<div dir="ltr"><div dir="ltr">Dear Roland,<div><br></div><div>Firstly, thank you for your considerations and detailed response. Also, thank you for forwarding the message. I will first try to use TOVSolverHot and make it work as you have suggested.</div><div><br></div><div>Best regards,</div><div><br></div><div>Beyhan.</div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Jun 21, 2020 at 9:56 PM Roland Haas <<a href="mailto:rhaas@illinois.edu">rhaas@illinois.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Hello Beyhan,<br>
<br>
> I am using Turing release on Fedora 31. I have been trying to<br>
> perform a binary neutron star merger simulation via finite<br>
> temperature equation of state (might be misleading) with GRHydro. As<br>
> far as I can see, this question was replied by Roland previously in<br>
> <a href="http://lists.einsteintoolkit.org/pipermail/users/2019-February/006756.html" rel="noreferrer" target="_blank">http://lists.einsteintoolkit.org/pipermail/users/2019-February/006756.html</a>.<br>
> <br>
> Parameters were added as has been suggested but If Meudon_BNS_ID does<br>
> not refer to source code of Meudon_Bin_NS then I could not find where<br>
> Meudon_BNS_ID or Meudon_Bin_NS ID is located in Cactus, if it does I<br>
> will try to work on it.<br>
Meudon_BNS_ID refers to the binary neutron star ID reader for LORENE<br>
that is found in the toolkit. Indeed the correct name is Meudon_Bin_NS .<br>
<br>
> In this setting, is it enough to use IOUtil parameter<br>
> -filereader_ID_vars = "all" - to allow initial temperature, entropy and<br>
> electron fraction to be read from the file?<br>
Not quite. filereader reads in data from Cactus 3d output files.<br>
However that is not what is done with the Meudon thorns. Instead they<br>
are regular initial data thorns that read in data from LORENE's output<br>
files. <br>
<br>
Looking at the file Meudon_Bin_NS/src/Bin_NS.cc then it will set<br>
HydroBase::temperature and HydroBase::Y_e if those variables have<br>
storage eg if you add<br>
<br>
HydroBase::initial_temperature = "read from file"<br>
HydroBase::initial_Y_e = "read from file"<br>
<br>
Really though the Meudon thorn should extend the valid ranges of those<br>
keywords to include "Meudon_Bin_NS" the way tit does for the other<br>
initial_XXX parameters.<br>
<br>
I noticed that there is a TOVSolverHot thorn in<br>
<a href="http://svn.einsteintoolkit.org/incoming/TOVSolverHot/" rel="noreferrer" target="_blank">http://svn.einsteintoolkit.org/incoming/TOVSolverHot/</a> that lets you<br>
construct initial data for an isothermal TOV star and that may be a an<br>
easier test before going for a full binary simulation.<br>
<br>
> I have tried to use APR eos with SFH0 eos table from<br>
> <a href="http://stellarcollapse.org" rel="noreferrer" target="_blank">stellarcollapse.org</a> (<br>
> <a href="https://stellarcollapse.org/~evanoc/Hempel_SFHoEOS_rho222_temp180_ye60_version_1.1_20120817.h5.bz2" rel="noreferrer" target="_blank">https://stellarcollapse.org/~evanoc/Hempel_SFHoEOS_rho222_temp180_ye60_version_1.1_20120817.h5.bz2</a><br>
> ) and got "EOS errori n prim2con_hot" in this setting and do not know<br>
> how this might be overcomed.<br>
There is not usually a really good general solution for this. It boils<br>
down to gut feeling and experience with this particular type of<br>
simulations. I would give the isothermal single TOV a try first and<br>
once that works work my way up to a binary.<br>
<br>
> Besides that, I suppose WhiskyTHC allows nuclear eos table to be used<br>
> with Lorene initial data by providing this table in pizza format and<br>
> uses some conversion of units but do not know what will be needed for<br>
> Einstein Toolkit.<br>
WhiskyTHC has a full example on how to use it in its tarball<br>
<a href="http://personal.psu.edu/dur566/downloads/WhiskyTHC.tar" rel="noreferrer" target="_blank">http://personal.psu.edu/dur566/downloads/WhiskyTHC.tar</a>. Its output is<br>
in the regular ET HydroBase variables so not really different from when<br>
using GRHydro.<br>
<br>
Yours,<br>
Roland<br>
<br>
-- <br>
My email is as private as my paper mail. I therefore support encrypting<br>
and signing email messages. Get my PGP key from <a href="http://pgp.mit.edu" rel="noreferrer" target="_blank">http://pgp.mit.edu</a> .<br>
</blockquote></div></div>