<div dir="ltr">let me add the AHFinderDirect correctly finds an AH after merger (the one that should be stored in surface 2).<div><br></div><div>Moreover, Federico also run a sim in which AHFinderDirect was writing on surfaces 0, 1, and 2 while PunctureTracker was writing on surfaces 3 and 4. In this case Outflow was giving the same error on all three surfaces 0, 1, and 2.</div><div><br></div><div>My naive impression is that we are missing something and AHFinderDirect is not saving the surfaces properly and therefore Outflow is not able to read them.</div><div><br></div><div>I had a look at previous par files from my BNS simulations (where Outflow worked) and the only difference I saw was in the number of points in theta and phi for the surfaces.</div><div><br></div><div>Thanks,</div><div>Bruno</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 20 feb 2020 alle ore 18:37 Federico Cattorini <<a href="mailto:f.cattorini@campus.unimib.it">f.cattorini@campus.unimib.it</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div><span style="font-family:arial,sans-serif"><font size="2">Dear all,</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">I am having issues with the thorn Outflow running a BHB simulation with WhiskyMHD (parfile attached). Outflow correctly computes the flows on the single orbiting BHs during the inspiral, but after the merger it fails to compute it on the newly formed BH. <br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">I get the following warning in the .out:</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">WARNING level 1 from host <a href="http://r149c18s02.marconi.cineca.it" target="_blank">r149c18s02.marconi.cineca.it</a> process 0<br> while executing schedule bin CCTK_ANALYSIS, routine Outflow::outflow<br> in thorn Outflow, file /marconi/home/userexternal/fcattori/EinsteinToolkit/ET_2018_09/Cactus/arrangements/EinsteinAnalysis/Outflow/src/outflow.c:958:<br> -> didn't find valid detector surface for sn=2, det=2</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">In my parfile I implement Outflow as</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">outflow::compute_every = 64 <br>outflow::num_detectors = 3 <br>outflow::surface_index[0] = 0 <br>outflow::surface_index[1] = 1 <br>outflow::surface_index[2] = 2 <br>outflow::interpolator_name = "Lagrange polynomial interpolation" <br>outflow::interpolator_pars = "order=4" <br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2">AHFinderDirect uses three spherical surfaces (n. 0,1,2). PunctureTracker tracks and write its values on n. 0 and 1.</font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><pre><span style="font-family:arial,sans-serif"><font size="2">Does anyone know how to address this issue?<br><br></font></span></pre><pre><span style="font-family:arial,sans-serif"><font size="2">Thanks a lot,<br><br></font></span></pre><pre><span style="font-family:arial,sans-serif"><font size="2">Federico<br></font></span></pre></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div><div><span style="font-family:arial,sans-serif"><font size="2"><br></font></span></div></div>
_______________________________________________<br>
Users mailing list<br>
<a href="mailto:Users@einsteintoolkit.org" target="_blank">Users@einsteintoolkit.org</a><br>
<a href="http://lists.einsteintoolkit.org/mailman/listinfo/users" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/mailman/listinfo/users</a><br>
</blockquote></div><br clear="all"><div><br></div>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><p><font color="#000000">Pr<font face="arial, sans-serif">of. Bruno Giacomazzo<br>Department of Physics<br>University of Milano-Bicocca<br></font></font><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">Piazza della Scienza 3<br></font></span><span style="color:rgb(0,0,0)"><font face="arial, sans-serif">20126 Milano<br></font></span><span style="font-size:12.8px;color:rgb(0,0,0)">Italy</span></p><p><font color="#000000"><span style="font-size:12.8px">email: </span><span style="font-size:12.8px"><a href="mailto:bruno.giacomazzo@unimib.it" target="_blank">bruno.giacomazzo@unimib.it</a><br></span><span style="font-size:12.8px">phone: (+39) 02 6448 2321</span></font><br><font color="#000000"><span style="font-size:12.8px">web: </span></font><span style="color:rgb(0,0,0);font-size:12.8px"><a href="http://www.brunogiacomazzo.org/" style="font-size:12.8px" target="_blank">http://www.brunogiacomazzo.org</a></span></p><p><font color="#000000">----------------------------------------------------------------------<br><span style="font-size:12.8px">There are only 10 types of people in the world:<br></span><span style="font-size:12.8px">Those who understand binary, and those who don't<br></span><span style="font-size:12.8px">----------------------------------------------------------------------<br></span></font></p></div></div></div></div></div></div></div>