<div dir="ltr">Klaus,<div><br></div><div>It's failing in the actual Eos routine.  I can trace the abort message that you see in the log below back to ../multiTemp/Multitype/Eos.F90.  It's being given input of mode MODE_DENS_ TEMP_ELE when it's called in rt_computeDt.F90.  The select case statement which starts on line 495 doesn't have a case specification for that mode, so it goes to the default case which aborts.<br><div><br></div><div><u>Setup command line</u>: </div><div>./setup Burst_Rad -auto -1d -nxb=10 +spherical +mtmmmt +mgd +uhd3t mgd_meshgroups=6 species=air<br></div><div><br></div><div><u>Relevant units to my question</u>:</div><div>Multispecies/MultispeciesMain<br>physics/Diffuse/DiffuseMain/Unsplit<br>physics/Diffuse/localAPI<br>physics/Eos/EosMain/Tabulated/Hdf5TableRead<br>physics/Eos/EosMain/multiTemp/Gamma/Ye<br>physics/Eos/EosMain/multiTemp/Multitype<br>physics/Eos/localAPI<br>physics/Hydro/HydroMain/unsplit/Hydro_Unsplit<br>physics/Hydro/HydroMain/unsplit/multiTemp<br>physics/Hydro/localAPI<br>physics/RadTrans/RadTransMain/MGD<br>physics/RadTrans/localAPI<br>physics/materialProperties/Opacity/OpacityMain/Multispecies/method/Constcm2g<br>physics/materialProperties/Opacity/localAPI<br></div><div><br></div><div><u>Relevant flash.par flags</u>:</div><div># EoS parameters (nitrogen)<br>ms_airA = 14.0067<br>ms_airZ = 7.0<br>eos_airEosType = "eos_gam"<br><br># Diffusion parameters<br>useDiffuse = .true.<br><br># Radiation (multi-group diffusion) parameters<br>useRadTrans = .true.<br>rt_useMGD = .true.<br>rt_mgdNumGroups = 6<br>rt_mgdBounds_1 = 0.000000e+00<br>rt_mgdBounds_2 = 9.992797e-01<br>rt_mgdBounds_3 = 9.992797e+00<br>rt_mgdBounds_4 = 9.992797e+01<br>rt_mgdBounds_5 = 9.992797e+02<br>rt_mgdBounds_6 = 9.992797e+03<br>rt_mgdBounds_7 = 9.992797e+04<br>rt_mgdFlMode = "fl_levermorepomraning1981"<br>rt_computeDt = .true.<br><br># Opacity parameters (general)<br>useOpacity      = .true.<br><br># Opacity parameters (air)<br>op_airAbsorb = "op_constcm2g"<br>op_airEmiss = "op_constcm2g"<br>op_airTrans = "op_constcm2g"<br>op_airAbsorbConstant = 1.e-1<br>op_airEmissConstant = 1.e-1<br>op_airTransConstant = 1.e-1<br><br># Heat exchange parameters<br>useHeatexchange = .false.<br></div><div><br></div><div><u>Error message from log</u>:</div><div>[ 12-11-2019  07:47:52.744 ] [Driver_evolveFlash]: Entering evolution loop<br> [ 12-11-2019  07:47:52.744 ] step: n=1 t=0.000000E+00 dt=1.000000E-15<br> [ 12-11-2019  07:47:52.824 ] [DRIVER_ABORT]: Driver_abort() called by PE           0<br> [ 12-11-2019  07:47:52.824 ] abort_message: [Eos] Unrecognized input mode given to Eos</div><div><br></div><div>Keith</div><div></div></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 10, 2019 at 5:52 PM Klaus Weide <<a href="mailto:klaus@flash.uchicago.edu">klaus@flash.uchicago.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">On Tue, 10 Dec 2019, Keith Galvin wrote:<br>
<br>
> Hi,<br>
> <br>
> I'm performing a 3T simulation using the Multitemp/Multitype EoS and I seem<br>
> to be running into issues when trying to compute a radiation time step.<br>
> The code seems to crash at the first time step inside of the EoS call when<br>
> computing the rad time step.  Specifically, it's complaining about an<br>
> invalid EoS mode.  After looking through the source code, it looks as if<br>
> the Multitemp/Multitype EoS routine isn't designed to handle the mode input<br>
> of MODE_DENS_TEMP_ELE, which is the mode specified in the routine that<br>
> calculates the rad time step.  Is the rt_computeDt code only designed to be<br>
> used with the RadFLAH versions of the EoS routines, or is there some way to<br>
> use it for full 3T calculations with multiple types?  Thanks!<br>
<br>
<br>
Keith,<br>
<br>
You should give more information on your setup.<br>
<br>
What is your setup command line?<br>
<br>
What is the list of units (as written by the 'setup' command)?<br>
<br>
If you can tell, which is the routine that makes the failing Eos or <br>
Eos_wrapped call?<br>
<br>
Klaus<br>
</blockquote></div>