<html>
<head>
<meta http-equiv="content-type" content="text/html; charset=UTF-8">
</head>
<body>
<p><font color="#0c0068">Hello all</font></p>
<p><font color="#0c0068">We have compiled and ran a customised
version Laserslab with the example3d.par and everything seems to
work. On the contrary when we want to use hydrogen instead of
helium and/or aluminum we encounter a number of errors,
different ones each time we use a different EOS tab file.<br>
</font></p>
<ul>
<li><font color="#0c0068">IONMIX h-100grp-lte.imx</font></li>
</ul>
<p><font color="#0c0068">We have set it up with the ./setup
custom_Laserslab -auto -3d +hdf5typeio species=cham,targ +mtmmmt
+laser +uhd3t +mgd mgd_meshgroups=100 -parfile=custom_3d.par
+pm4dev</font></p>
<p><font color="#0c0068">and I take a succesfull make.</font></p>
<p><font color="#0c0068">In the custom_3d.par i have declared 100</font></p>
<p><font color="#0c0068">rt_mgdNumGroups = 100<br>
#rt_mgdBoundEntry = "grbd_manual"<br>
rt_mgdBounds_1 = 0.0<br>
rt_mgdBounds_2 = 0.00139292180089<br>
rt_mgdBounds_3 = 0.001564702535</font></p>
<p><font color="#0c0068">...</font></p>
<p><font color="#0c0068">rt_mgdBounds_100 = 124.0/<br>
</font></p>
<p><font color="#0c0068">rt_mgdBounds_101 = 10000.</font></p>
<p><font color="#0c0068">acquired from the
FLASH4.6.1/source/Simulation/SimulationMain/radflaHD/SupernovaRad1D/flash_working.par<br>
</font></p>
<p><font color="#0c0068">Then when I ran it I receive the following
error, relevant to IONMIX / Opacity energy group mismatch, bold
in the output flow! Lasslab.log attached for further
information!</font></p>
<p><font color="#5a5a5a">[itazes@login01 object]$ mpirun -np 2
/users/pr007/itazes/FLASH4.6.1/object/flash4<br>
...</font></p>
<p><font color="#5a5a5a">Grid_init: resolution based on runtime
params:<br>
lrefine dx dy dz<br>
3 0.003 0.003 0.003<br>
4 0.001 0.001 0.001<br>
5 6.250E-04 6.250E-04 6.250E-04<br>
MaterialProperties initialized<br>
Cosmology initialized<br>
[eos_tabBrowseIonmixTables] IONMIX file found: h-100grp-lte.imx<br>
Grid_init: resolution based on runtime params:<br>
lrefine dx dy dz<br>
3 0.003 0.003 0.003<br>
4 0.001 0.001 0.001<br>
5 6.250E-04 6.250E-04 6.250E-04<br>
MaterialProperties initialized<br>
Cosmology initialized<br>
[eos_tabBrowseIonmixTables] IONMIX file found: h-100grp-lte.imx<br>
in eos_inittabulated, tableName = h-100grp-lte.imx<br>
in eos_inittabulated, groupName = -none-<br>
in eos_inittabulated, tableName = h-100grp-lte.imx<br>
in eos_inittabulated, groupName = -none-<br>
[eos_tabBrowseIonmixTables] IONMIX file found: h-100grp-lte.imx<br>
[eos_tabBrowseIonmixTables] IONMIX file found: h-100grp-lte.imx<br>
in eos_inittabulated, tableName = h-100grp-lte.imx<br>
in eos_inittabulated, groupName = -none-<br>
in eos_inittabulated, tableName = h-100grp-lte.imx<br>
in eos_inittabulated, groupName = -none-<br>
RadTrans initialized<br>
[EnergyDeposition_init] INFO: Using ed_irradVar= 16
ed_irradVarName=lase<br>
RadTrans initialized<br>
[EnergyDeposition_init] INFO: Using ed_irradVar= 16
ed_irradVarName=lase<br>
Source terms initialized<br>
Source terms initialized<br>
iteration, no. not moved = 0 0<br>
refined: total leaf blocks = 8<br>
refined: total blocks = 9<br>
[amr_morton_process]: Initializing surr_blks using standard
orrery implementation<br>
iteration, no. not moved = 0 0<br>
refined: total leaf blocks = 8<br>
refined: total blocks = 9<br>
[amr_morton_process]: Initializing surr_blks using standard
orrery implementation<br>
INFO: Grid_fillGuardCells is ignoring masking.<br>
INFO: Grid_fillGuardCells is ignoring masking.<br>
iteration, no. not moved = 0 0<br>
refined: total leaf blocks = 64<br>
refined: total blocks = 73<br>
iteration, no. not moved = 0 0<br>
refined: total leaf blocks = 64<br>
refined: total blocks = 73<br>
Finished with Grid_initDomain, no restart<br>
<br>
Driver_abort called. See log file for details.<br>
<b> Error message is [op_readIonmixTables] ERROR: IONMIX /
Opacity energy group mismatch<br>
Calling MPI_Abort() for shutdown in 2 seconds!</b><br>
<br>
Finished with Grid_initDomain, no restart<br>
<br>
Driver_abort called. See log file for details.<br>
<b> Error message is [op_readIonmixTables] ERROR: IONMIX /
Opacity energy group mismatch</b><b><br>
</b><b> Calling MPI_Abort() for shutdown in 2 seconds</b>!<br>
<br>
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0<br>
--------------------------------------------------------------------------<br>
Primary job terminated normally, but 1 process returned<br>
a non-zero exit code. Per user-direction, the job has been
aborted.<br>
--------------------------------------------------------------------------<br>
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0<br>
--------------------------------------------------------------------------<br>
mpirun detected that one or more processes exited with non-zero
status, thus causing<br>
the job to be terminated. The first process to do so was:<br>
<br>
Process name: [[54963,1],0]<br>
Exit code: 1<br>
--------------------------------------------------------------------------<br>
[itazes@login01 object]$</font></p>
<p><font color="#5a5a5a"><br>
</font></p>
<ul>
<li><font color="#0c0068"><font color="#0c0068">IONMIX4
h-imx-1grp.cn4</font></font></li>
</ul>
<p><font color="#0c0068"><font color="#0c0068"></font>Furthermore,
we decided to use the IONMIX4 hydrogen tabulated EOS file,
h-imx-1grp.cn4 setting it up with <br>
</font></p>
<p><font color="#0c0068">We have set it up with the ./setup
custom_Laserslab -auto -3d +hdf5typeio species=cham,targ +mtmmmt
+laser +uhd3t +mgd mgd_meshgroups=1 -parfile=custom_3d.par
+pm4dev</font></p>
<p><font color="#0c0068">In this case we managed to surpass the
IONMIX/ Opacity energy group mismatch error in this case only 2
boundaries in our .par where used as following:<br>
</font></p>
<p><font color="#0c0068">rt_mgdNumGroups = 1<br>
#rt_mgdBoundEntry = "grbd_manual"<br>
rt_mgdBounds_1 = 1.0e-2<br>
rt_mgdBounds_2 = 1<br>
</font></p>
<p><font color="#0c0068">and we encounter another error:</font></p>
<p><font color="#0c0068">Error message is [Diffuse]: computed dt is
not positive! Aborting!<br>
Calling MPI_Abort() for shutdown in 2 seconds,</font></p>
<p><font color="#0c0068"> bold in the output flow! Lasslab.log
attached for further information!</font></p>
<p><font color="#0c0068">We have already tried to play with the CFL
value and lower it but changed in out output</font></p>
<p><font color="#0c0068">also we changed the value of <br>
</font></p>
<p><font color="#0c0068">tdt_diff_factor = 1.0e+100 # Disable
diffusion dt</font></p>
<p><font color="#0c0068">but we do not really understand its
purpose. FYI we attach the flash.par we used<br>
</font></p>
<p><font color="#5a5a5a">[itazes@login01 object]$ mpirun -np 1
/users/pr007/itazes/FLASH4.6.1/object/flash4<br>
RuntimeParameters_read: ignoring unknown parameter
"iGridSize"...<br>
RuntimeParameters_read: ignoring unknown parameter
"jGridSize"...<br>
RuntimeParameters_read: ignoring unknown parameter
"kGridSize"...<br>
Grid_init: resolution based on runtime params:<br>
lrefine dx dy dz<br>
3 0.003 0.003 0.003<br>
4 0.001 0.001 0.001<br>
5 6.250E-04 6.250E-04 6.250E-04<br>
MaterialProperties initialized<br>
Cosmology initialized<br>
[eos_tabBrowseIonmix4Tables] IONMIX4 file found: h-imx-1grp.cn4<br>
[eos_tabBrowseIonmix4Tables] IONMIX4 file found: h-imx-1grp.cn4<br>
in eos_inittabulated, tableName = h-imx-1grp.cn4<br>
in eos_inittabulated, groupName = -none-<br>
in eos_inittabulated, tableName = h-imx-1grp.cn4<br>
in eos_inittabulated, groupName = -none-<br>
RadTrans initialized<br>
[EnergyDeposition_init] INFO: Using ed_irradVar= 16
ed_irradVarName=lase<br>
Source terms initialized<br>
iteration, no. not moved = 0 0<br>
refined: total leaf blocks = 8<br>
refined: total blocks = 9<br>
[amr_morton_process]: Initializing surr_blks using standard
orrery implementation<br>
INFO: Grid_fillGuardCells is ignoring masking.<br>
iteration, no. not moved = 0 0<br>
refined: total leaf blocks = 64<br>
refined: total blocks = 73<br>
Finished with Grid_initDomain, no restart<br>
Ready to call Hydro_init<br>
[Hydro_init] INFO: Turning shockDetect on because RiemannSolver
is hybrid.<br>
Hydro initialized<br>
Gravity initialized<br>
<br>
Driver_abort called. See log file for details.<br>
<b> Error message is [Diffuse]: computed dt is not positive!
Aborting!</b><b><br>
</b><b> Calling MPI_Abort() for shutdown in 2 seconds!</b><br>
<br>
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0<br>
--------------------------------------------------------------------------<br>
Primary job terminated normally, but 1 process returned<br>
a non-zero exit code. Per user-direction, the job has been
aborted.<br>
--------------------------------------------------------------------------<br>
--------------------------------------------------------------------------<br>
mpirun detected that one or more processes exited with non-zero
status, thus causing<br>
the job to be terminated. The first process to do so was:<br>
<br>
Process name: [[52177,1],0]<br>
Exit code: 1<br>
--------------------------------------------------------------------------<br>
[itazes@login01 object]$</font><br>
</p>
<p>Stelios</p>
<pre class="moz-signature" cols="72">--
Stylianos Passalidis
Phd candidate
Laboratory of Physical Chemistry - Matter and Radiation
Sorbonne Université
4 place Jussieu, 75005 Paris
Bureau 111, Corridor 43-44</pre>
<p><br>
</p>
</body>
</html>