Sphere-In-Box issue with EOS selection

[Geoffrey Bryan]

Hello all,

 I have a problem with running the sphere-in-box simulation   "out-of-the-box". If I use phantomsetup to generate the setup file "sinb.in"  interactively the option chosen for the equation of state (eos) is  ieos=8 (8) Barotropic - which is fine as this is appropriate as it is used in the paper by Price and Bates.

However, if I then run phantom using  the command ./phantom sinb  I get the following error message:

FATAL ERROR! read_infile: only ieos=2 makes sense if storing thermal energy

I can change the value of ieos to 2 (adiabatic) and the code will run but this is the wrong equation of state!

Also the sinb.in input file generated by ./phantomsetup does not contain some of the vital values for the simulation (such as masstoflux) which means that their values default to the values stored in the code.  So can I add lines to the initialisation file (in this case "sinb.in") such as:

masstoflux = 10.0 ?

I'd be appreciative of any help here as this is the first time I've run phantom in earnest (I've been to the summer school in Shanghai but that was a while ago.....)

Cheers,

 Geoff

[James Wurster]

Hello Geoff

I assume you are running phantom using setup=sphereinbox.  To reproduce Price & Bate (2007), I suggest using setup=jet.  This will automatically turn on the relevant processes. 

For completeness, to use ieos=8, you need to compile with ISOTHERMAL=yes (ie. add this to the command line when you type make).

For the mass-to-flux ratio, you should be prompted for this value when you run phantomsetup; you can then modify it in sinb.setup and rerun phantomsetup if you choose.  You will only be prompted for this if you compile with MHD=yes, which is the default for setup=jet but not setup=sphereinbox.

Hope that this helps.

James

--
You received this message because you are subscribed to the Google Groups "Phantom users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to phantomsph-use...@googlegroups.com.
To post to this group, send email to phantoms...@googlegroups.com.
To view this discussion on the web, visit https://groups.google.com/d/msgid/phantomsph-users/817ce985-619b-4c5b-a9ee-19b9235aaf44%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Geoffrey Bryan]

Thanks James,

 That helped a lot!

 Cheers,