<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><br><div><div>On 23 Jan 2012, at 14:57, Jose Fiestas Iquira wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><div>Hello,<br><br>Regarding the Binary Black Hole coalescence example from the ETK.<br><br>In the BH diagnostics files are columns with quadrupole terms, which<br>are not constant. Do you know of some tool in ETK to calculate the GW<br>emission from these terms?<br></div></blockquote><div><br></div><div>No, not that I am aware of. These are the quadrupole moments of the individual apparent horizons, rather than the quadrupole moment of the system as a whole. The gravitational wave emission can be computed by using the WeylScal4 and Multipole thorns. These give you the spherical harmonic multipoles of Psi4 (e.g. l = 2, m = 2 is the dominant one), which can then be converted to strain by "two integrations in time" [this is actually quite subtle and hard to get right - see Reisswig and Pollney, <span class="Apple-style-span" style="font-family: arial, verdana, sans-serif; ">arXiv:1006.1632, for an overview and solution).</span></div><div><br></div><br><blockquote type="cite"><div><br>Thanks,<br>Jose<br><br>On Tue, Jan 17, 2012 at 5:18 AM, Ian Hinder <<a href="mailto:ian.hinder@aei.mpg.de">ian.hinder@aei.mpg.de</a>> wrote:<br><blockquote type="cite"><br></blockquote><blockquote type="cite">On 17 Jan 2012, at 13:37, Jose Fiestas Iquira wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The Einstein Analysis guide is very useful. I found here<br></blockquote><blockquote type="cite"><a href="http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech142.html#x159-2520000P8.4">http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech142.html#x159-2520000P8.4</a><br></blockquote><blockquote type="cite">the units used in 'Hydrobase'.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">If I assume these are the units for all thorns, I obtain from data in<br></blockquote><blockquote type="cite">BH_diagnostics.ah1.gp a unit of distance of about 1 Km. To estimate<br></blockquote><blockquote type="cite">the time I am using a mass of ~0.5 M_sun for both BHs (m_irreducible<br></blockquote><blockquote type="cite">?) , and I read a final time (cctk_time) ~ 1microsec. Are these<br></blockquote><blockquote type="cite">quantities in the usual range?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The length unit is G M / c^2 and the time unit is G M / c^3. For M = M_sun,<br></blockquote><blockquote type="cite">I get the length unit as 1.5 km and the time unit as 4.9 us. The final time<br></blockquote><blockquote type="cite">in the simulation is 100 M (i.e. 100 G M / c^3), which would be 0.5 ms.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Tue, Jan 17, 2012 at 4:04 AM, Jose Fiestas Iquira <<a href="mailto:jafiestas@lbl.gov">jafiestas@lbl.gov</a>><br></blockquote><blockquote type="cite">wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I sent the wrong file, here is the right one.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">$ hg tip<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> gives me:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-----------------------------------------------------<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">changeset: 3421:dff4df2b3687<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">tag: tip<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">user: Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu">schnetter@cct.lsu.edu</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">date: Thu Jan 12 14:11:27 2012 -0500<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">summary: CarpetLib: Allocate space for dstbox before<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">restricting/prolongating<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">[40514 refs]<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-----------------------------------------------------<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I copied this version from:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">./GetComponents -a<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="http://svn.einsteintoolkit.org/manifest/trunk/einsteintoolkit.th">http://svn.einsteintoolkit.org/manifest/trunk/einsteintoolkit.th</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Tue, Jan 17, 2012 at 1:13 AM, Ian Hinder <<a href="mailto:ian.hinder@aei.mpg.de">ian.hinder@aei.mpg.de</a>> wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On 16 Jan 2012, at 23:21, Jose Fiestas Iquira wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks, I am attaching the output file.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The output file looks fine - I can't see any problem with the simulation.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> According to the output file, it reached t = 100 M, and in the parameter<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">file that is listed as cctk_final_time, so all looks well!<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I did copy the development version. As I can see, it was already<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">fixed, so I should not get this warning, since I copied the<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">development version.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Yes, if you are using the current development version, then I don't know why<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">you get the warning. Can you check which version of Carpet you have by<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">doing<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">cd repos/carpet<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">hg tip<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">to be sure you are using the development version?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Is the puncture tracker not a better way to plot BH trajectories.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">It's cheaper to use the puncture tracker, because computing the horizon<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">centroid is more expensive, but it has to be done every iteration, whereas<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">the horizon can be found less frequently. For comparison with<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">point-particle results, I don't think it makes much difference; the results<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">are gauge dependent anyway so differences on the order of the horizon size<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">cannot be interpreted physically. Since the parameter file already has<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">horizon-finding built in, it is easier for you to use that.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">If you want to figure out how to add the puncture tracker to the simulation,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">you can start with these parameters:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">ActiveThorns = "PunctureTracker"<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PunctureTracker::track [0] = yes<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PunctureTracker::initial_x [0] = $xp<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PunctureTracker::which_surface_to_store_info[0] = 0<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PunctureTracker::track [1] = yes<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PunctureTracker::initial_x [1] = $xm<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">PunctureTracker::which_surface_to_store_info[1] = 1<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">IOASCII::out0D_vars = "PunctureTracker::pt_loc"<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Replace $xp and $xm with the initial x coordinates of the +x and -x<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">punctures (this will be \pm TwoPunctures::par_b). You then need to choose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">two spherical surfaces<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">(<a href="https://docs.einsteintoolkit.org/et-docs/Einstein_Toolkit_standards">https://docs.einsteintoolkit.org/et-docs/Einstein_Toolkit_standards</a>) which<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">are not being used currently to store the puncture locations in. If there<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">is already 0D output, just add the puncturetracker variable separated with a<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">space.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">did some plots, what are really the combinations xx, xy, xz, yy, zz ?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Are they related to the trajectories of one BH (projections)? What<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">about the second BH?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">You should have BH_diagnostics.ah1.gp and BH_diagnostics.ah2.gp. These<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">correspond to the first and second BHs. If you look at the top of these<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">files, you will see an explanation of the columns. More information about<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">each column is available in the AHFinderDirect documentation<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">(<a href="http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech128.html#x144-2264000O4">http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech128.html#x144-2264000O4</a>).<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> Specifically, information about the BH_diagnostics output files is<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">at <a href="http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech128.html#x144-2274000O4.5.5">http://einsteintoolkit.org/documentation/ThornGuide/ThornGuidech128.html#x144-2274000O4.5.5</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">To plot the trajectory of one of the BHs, you want to plot column 3 and<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">column 4 in a parametric plot.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Cheers,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Mon, Jan 16, 2012 at 4:39 AM, Ian Hinder <<a href="mailto:ian.hinder@aei.mpg.de">ian.hinder@aei.mpg.de</a>> wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On 14 Jan 2012, at 00:41, Jose Fiestas Iquira wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Hi Ian,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I did the changes but it started and finished early. I am attaching<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">the qc0-mclachlan.err file. At the end there is a warning in one of<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">the Carpet files (mask_test.c)<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The warning you see,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">WARNING level 1 in thorn CarpetReduce processor 0 host qb517<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> (line 120 of<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">/home/ettest24/Cactus/arrangements/Carpet/CarpetReduce/src/mask_test.c):<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> -> Simulation domain volume and reduction weight sum differ<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">is a known issue (<a href="https://trac.einsteintoolkit.org/ticket/434">https://trac.einsteintoolkit.org/ticket/434</a>) which was<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">fixed after the last release and the fix is only available in the<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">development version. I believe that the warning is not fatal, and it should<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">not have terminated the simulation. Can you also post the qc0-mclachlan.out<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">file?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Where can I find the puncture tracker?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The PunctureTracker thorn is part of the Einstein Toolkit thornlist, so you<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">should have it already. However, I think the best thing for you to do<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">initially would be to look at the BH_diagnostics files for the apparent<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">horizon centroid locations for tracking the locations of the BHs.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Fri, Jan 13, 2012 at 4:38 AM, Ian Hinder <<a href="mailto:ian.hinder@aei.mpg.de">ian.hinder@aei.mpg.de</a>> wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On 13 Jan 2012, at 12:32, Jose Fiestas Iquira wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Yes, trajectories and wave forms would be enough for now. Which are<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">the names of this files? Many of the .asc files have data like kxx,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">kxy, kxz, etc. Is there any description of them?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Not that I am aware of. Files with extension .x.asc are 1D ascii output of<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">grid functions from CarpetIOASCII. kxy is the variable name for the x,y<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">component of the extrinsic curvature tensor K. The waveform file should be<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">in mp_psi4_l2_m2_r*.asc where * corresponds to the radius of the extraction<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">sphere. However, looking at the qc0-mclachlan.par file, I see that there is<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">no wave extraction switched on. You should be able to enable wave<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">extraction by adding something like the following to the parameter file:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">ActiveThorns = "WeylScal4 CarpetInterp AEILocalInterp Multipole"<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Multipole::nradii = 3<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Multipole::radius[0] = 30<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Multipole::radius[1] = 40<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Multipole::radius[2] = 50<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Multipole::variables = "WeylScal4::Psi4r{sw=-2 cmplx=’WeylScal4::Psi4i’}"<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Multipole::l_max = 4<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">WeylScal4::fdOrder = 4<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">(I haven't tested this, so let me know if it doesn't work)<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I've just found that the usual PunctureTracker thorn is not enabled in this<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">example parameter file. So to look at trajectories you will have to use the<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">apparent horizon centroids. These are in the files BH_diagnostics.ah*.gp.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> The meanings of the columns are at the top of the files. You want the<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">centroid x and y columns to plot in GnuPlot.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Fri, Jan 13, 2012 at 12:41 AM, Ian Hinder <<a href="mailto:ian.hinder@aei.mpg.de">ian.hinder@aei.mpg.de</a>> wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On 12 Jan 2012, at 23:46, Jose Fiestas Iquira wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Dear all,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">I run the tutorial and found other interesting examples like the<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">McLachlan tool for Binary Black Holes.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Regarding the McLachlan tool,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Could somebody help me finding a description of the output? Is all<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">contained in the .asc files? I would like to visualize it.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Hi Jose,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">What sort of output do you want to visualise? Cactus thorns typically<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">output "lightweight" information such as BH trajectories and waveforms as<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">ASCII, and more heavyweight data -- like gridfunctions -- as HDF5. For the<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">trajectories (from thorn Puncture Tracker) and waveforms (from Multipole)<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">you can use gnuplot or anything that can read standard ASCII output to<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">visualise them. For gridfunctions in HDF5 format, I think the easiest way<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">to get started is to use VisIt, for which there is a Carpet plugin (though<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">installing VisIt is sometimes not so easy!).<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Does anyone know if the Carpet plugin is distributed with VisIt yet? I<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">believe that was the intention.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Cheers,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Wed, Jan 11, 2012 at 7:06 PM, Jose Fiestas Iquira <<a href="mailto:jafiestas@lbl.gov">jafiestas@lbl.gov</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks to Erik for advise. ETK example simulation run successfully.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Wed, Jan 11, 2012 at 11:49 AM, Jose Fiestas Iquira <<a href="mailto:jafiestas@lbl.gov">jafiestas@lbl.gov</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Right, it looks good now:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">./simfactory/bin/sim list-simulations<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"> static_tov [ACTIVE (RUNNING), restart 0000, job id 595490]<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Let see how it looks when it finishes. The example was set to 8 hours.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Thanks,<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Jose<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Wed, Jan 11, 2012 at 11:45 AM, Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu">schnetter@cct.lsu.edu</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">On Wed, Jan 11, 2012 at 2:42 PM, Jose Fiestas Iquira <<a href="mailto:jafiestas@lbl.gov">jafiestas@lbl.gov</a>><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">wrote:<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">The only thing is that now<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">$ ./simfactory/bin/sim list-simulations<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">returns<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Error: unknown command list-simulation<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Aborting Simfactory.<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Did you really type "list-simulations", and the error message says<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">"list-simulation" (without the final "s")?<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">-erik<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Erik Schnetter <<a href="mailto:schnetter@cct.lsu.edu">schnetter@cct.lsu.edu</a>> <a href="http://www.cct.lsu.edu/~eschnett/">http://www.cct.lsu.edu/~eschnett/</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">_______________________________________________<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Users mailing list<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="mailto:Users@einsteintoolkit.org">Users@einsteintoolkit.org</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="http://lists.einsteintoolkit.org/mailman/listinfo/users">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Ian Hinder<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="http://numrel.aei.mpg.de/people/hinder">http://numrel.aei.mpg.de/people/hinder</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Ian Hinder<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="http://numrel.aei.mpg.de/people/hinder">http://numrel.aei.mpg.de/people/hinder</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><qc0-mclachlan.err><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Ian Hinder<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="http://numrel.aei.mpg.de/people/hinder">http://numrel.aei.mpg.de/people/hinder</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><qc0-mclachlan.out><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">Ian Hinder<br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><a href="http://numrel.aei.mpg.de/people/hinder">http://numrel.aei.mpg.de/people/hinder</a><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite"><br></blockquote><blockquote type="cite">--<br></blockquote><blockquote type="cite">Ian Hinder<br></blockquote><blockquote type="cite"><a href="http://numrel.aei.mpg.de/people/hinder">http://numrel.aei.mpg.de/people/hinder</a><br></blockquote><blockquote type="cite"><br></blockquote></div></blockquote></div><br><div>
<span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: -webkit-auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; font-size: medium; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>-- </div><div>Ian Hinder</div><div><a href="http://numrel.aei.mpg.de/people/hinder">http://numrel.aei.mpg.de/people/hinder</a></div></div></span></div></span></div></span></span>
</div>
<br></body></html>