[FLASH-USERS] ERROR: bad size of IONMIX4 energy group grid
Arthur Hirsch-Passicos
a.hirsch at hzdr.de
Mon Aug 12 07:52:56 EDT 2024
Dear Flash users,
I'm a new user of FLASH and am currently trying to launch my first few
non-example simulations. I am basing this case on the LaserSlab example
and am confronted to an error with the EOS tables of my targets. I tried
with several .cn4 files, but they all give me this error:
[eos_tabBrowseIonmix4Tables] IONMIX4 file found: he-imx-005.cn4
[eos_tabBrowseIonmix4Tables] IONMIX4 file found: h-imx-1grp.cn4
in eos_inittabulated, tableName = he-imx-005.cn4
in eos_inittabulated, groupName = -none-
in eos_inittabulated, tableName = h-imx-1grp.cn4
in eos_inittabulated, groupName = -none-
RadTrans initialized
[EnergyDeposition_init] INFO: Using ed_irradVar= 16
ed_irradVarName=lase
Source terms initialized
iteration, no. not moved = 0 0
refined: total leaf blocks = 2
refined: total blocks = 2
[amr_morton_process]: Initializing surr_blks using standard orrery
implementation
iteration, no. not moved = 0 0
refined: total leaf blocks = 8
refined: total blocks = 10
iteration, no. not moved = 0 1
iteration, no. not moved = 1 0
refined: total leaf blocks = 32
refined: total blocks = 42
iteration, no. not moved = 0 11
iteration, no. not moved = 1 0
refined: total leaf blocks = 92
refined: total blocks = 122
Finished with Grid_initDomain, no restart
1 6
Driver_abort called. See log file for details.
Error message is [op_readIonmix4Tables] ERROR: bad size of IONMIX4
energy group grid
I have no more idea of how to fix it. If anyone has encountered this,
how did you fix it? Thank you in advance for your answer.
Best regards,
Arthur Hirsch-Passicos
Post-doctoral Fellow
Helmhotz-Zentrum Dresden Rossendorf
-------------- next part --------------
A non-text attachment was scrubbed...
Name: lasjet.log
Type: text/x-log
Size: 144693 bytes
Desc: not available
URL: <http://flash.rochester.edu/pipermail/flash-users/attachments/20240812/4ce1fec4/attachment-0001.bin>
-------------- next part --------------
run_comment = "Laser Slab Example Simulation"
log_file = "lasjet.log"
basenm = "lasjet_"
# This particular parfile is used as an example that is described in
# detail in the users guide.
##########################
# #
# I/O PARAMETERS #
# #
##########################
### Checkpoint Options ###
checkpointFileIntervalTime = 1.0
checkpointFileIntervalStep = 1000
### Plot Options ###
plotFileNumber = 0
plotFileIntervalStep = 100
plotFileIntervalTime = 0.5e-12
plot_var_1 = "dens"
plot_var_2 = "depo"
plot_var_3 = "tele"
plot_var_4 = "tion"
plot_var_5 = "trad"
plot_var_6 = "ye "
plot_var_7 = "sumy"
plot_var_8 = "cham"
plot_var_9 = "targ"
plot_var_10 = "pres"
plot_var_11 = "bdry"
### Restart Options ###
restart = .false.
checkpointFileNumber = 0
########################################
# #
# RADIATION/OPACITY PARAMETERS #
# #
########################################
rt_useMGD = .true.
rt_mgdNumGroups = 6
rt_mgdBounds_1 = 1.0e-01
rt_mgdBounds_2 = 1.0e+00
rt_mgdBounds_3 = 1.0e+01
rt_mgdBounds_4 = 1.0e+02
rt_mgdBounds_5 = 1.0e+03
rt_mgdBounds_6 = 1.0e+04
rt_mgdBounds_7 = 1.0e+05
rt_mgdFlMode = "fl_harmonic"
rt_mgdFlCoef = 1.0
rt_mgdXlBoundaryType = "vacuum"
rt_mgdXrBoundaryType = "vacuum"
rt_mgdYlBoundaryType = "vacuum"
rt_mgdYrBoundaryType = "vacuum"
rt_mgdZlBoundaryType = "vacuum"
rt_mgdZrBoundaryType = "vacuum"
useOpacity = .true.
### SET CHAMBER (HELIUM) OPACITY OPTIONS ###
op_chamAbsorb = "op_tabpa"
op_chamEmiss = "op_tabpe"
op_chamTrans = "op_tabro"
op_chamFileType = "ionmix4"
op_chamFileName = "he-imx-005.cn4"
### SET TARGET (ALUMINUM) OPACITY OPTIONS ###
op_targAbsorb = "op_tabpa"
op_targEmiss = "op_tabpe"
op_targTrans = "op_tabro"
op_targFileType = "ionmix4"
op_targFileName = "h-imx-1grp.cn4"
############################
# #
# LASER PARAMETERS #
# #
############################
useEnergyDeposition = .true.
ed_maxRayCount = 10000
ed_gradOrder = 2
### LASER IO OPTIONS ###
ed_useLaserIO = .false.
ed_laserIOMaxNumberOfPositions = 10000
ed_laserIOMaxNumberOfRays = 256
### SETUP LASER PULSES ###
ed_numberOfPulses = 1
# Define Pulse 1:
ed_numberOfSections_1 = 128
ed_time_1_1 = 0.000000e+00
ed_power_1_1 = 3.481276e+05
ed_time_1_2 = 1.181102e-15
ed_power_1_2 = 5.982787e+05
ed_time_1_3 = 2.362205e-15
ed_power_1_3 = 1.019379e+06
ed_time_1_4 = 3.543307e-15
ed_power_1_4 = 1.722009e+06
ed_time_1_5 = 4.724409e-15
ed_power_1_5 = 2.884046e+06
ed_time_1_6 = 5.905512e-15
ed_power_1_6 = 4.788904e+06
ed_time_1_7 = 7.086614e-15
ed_power_1_7 = 7.883832e+06
ed_time_1_8 = 8.267717e-15
ed_power_1_8 = 1.286784e+07
ed_time_1_9 = 9.448819e-15
ed_power_1_9 = 2.082291e+07
ed_time_1_10 = 1.062992e-14
ed_power_1_10 = 3.340752e+07
ed_time_1_11 = 1.181102e-14
ed_power_1_11 = 5.313912e+07
ed_time_1_12 = 1.299213e-14
ed_power_1_12 = 8.380148e+07
ed_time_1_13 = 1.417323e-14
ed_power_1_13 = 1.310256e+08
ed_time_1_14 = 1.535433e-14
ed_power_1_14 = 2.031084e+08
ed_time_1_15 = 1.653543e-14
ed_power_1_15 = 3.121526e+08
ed_time_1_16 = 1.771654e-14
ed_power_1_16 = 4.756343e+08
ed_time_1_17 = 1.889764e-14
ed_power_1_17 = 7.185327e+08
ed_time_1_18 = 2.007874e-14
ed_power_1_18 = 1.076186e+09
ed_time_1_19 = 2.125984e-14
ed_power_1_19 = 1.598067e+09
ed_time_1_20 = 2.244094e-14
ed_power_1_20 = 2.352719e+09
ed_time_1_21 = 2.362205e-14
ed_power_1_21 = 3.434096e+09
ed_time_1_22 = 2.480315e-14
ed_power_1_22 = 4.969606e+09
ed_time_1_23 = 2.598425e-14
ed_power_1_23 = 7.130150e+09
ed_time_1_24 = 2.716535e-14
ed_power_1_24 = 1.014244e+10
ed_time_1_25 = 2.834646e-14
ed_power_1_25 = 1.430388e+10
ed_time_1_26 = 2.952756e-14
ed_power_1_26 = 2.000010e+10
ed_time_1_27 = 3.070866e-14
ed_power_1_27 = 2.772540e+10
ed_time_1_28 = 3.188976e-14
ed_power_1_28 = 3.810576e+10
ed_time_1_29 = 3.307087e-14
ed_power_1_29 = 5.192431e+10
ed_time_1_30 = 3.425197e-14
ed_power_1_30 = 7.014843e+10
ed_time_1_31 = 3.543307e-14
ed_power_1_31 = 9.395772e+10
ed_time_1_32 = 3.661417e-14
ed_power_1_32 = 1.247712e+11
ed_time_1_33 = 3.779528e-14
ed_power_1_33 = 1.642718e+11
ed_time_1_34 = 3.897638e-14
ed_power_1_34 = 2.144269e+11
ed_time_1_35 = 4.015748e-14
ed_power_1_35 = 2.774998e+11
ed_time_1_36 = 4.133858e-14
ed_power_1_36 = 3.560519e+11
ed_time_1_37 = 4.251969e-14
ed_power_1_37 = 4.529300e+11
ed_time_1_38 = 4.370079e-14
ed_power_1_38 = 5.712368e+11
ed_time_1_39 = 4.488189e-14
ed_power_1_39 = 7.142799e+11
ed_time_1_40 = 4.606299e-14
ed_power_1_40 = 8.854986e+11
ed_time_1_41 = 4.724409e-14
ed_power_1_41 = 1.088365e+12
ed_time_1_42 = 4.842520e-14
ed_power_1_42 = 1.326260e+12
ed_time_1_43 = 4.960630e-14
ed_power_1_43 = 1.602322e+12
ed_time_1_44 = 5.078740e-14
ed_power_1_44 = 1.919279e+12
ed_time_1_45 = 5.196850e-14
ed_power_1_45 = 2.279259e+12
ed_time_1_46 = 5.314961e-14
ed_power_1_46 = 2.683592e+12
ed_time_1_47 = 5.433071e-14
ed_power_1_47 = 3.132611e+12
ed_time_1_48 = 5.551181e-14
ed_power_1_48 = 3.625465e+12
ed_time_1_49 = 5.669291e-14
ed_power_1_49 = 4.159951e+12
ed_time_1_50 = 5.787402e-14
ed_power_1_50 = 4.732383e+12
ed_time_1_51 = 5.905512e-14
ed_power_1_51 = 5.337511e+12
ed_time_1_52 = 6.023622e-14
ed_power_1_52 = 5.968496e+12
ed_time_1_53 = 6.141732e-14
ed_power_1_53 = 6.616957e+12
ed_time_1_54 = 6.259843e-14
ed_power_1_54 = 7.273088e+12
ed_time_1_55 = 6.377953e-14
ed_power_1_55 = 7.925865e+12
ed_time_1_56 = 6.496063e-14
ed_power_1_56 = 8.563310e+12
ed_time_1_57 = 6.614173e-14
ed_power_1_57 = 9.172842e+12
ed_time_1_58 = 6.732283e-14
ed_power_1_58 = 9.741669e+12
ed_time_1_59 = 6.850394e-14
ed_power_1_59 = 1.025723e+13
ed_time_1_60 = 6.968504e-14
ed_power_1_60 = 1.070764e+13
ed_time_1_61 = 7.086614e-14
ed_power_1_61 = 1.108218e+13
ed_time_1_62 = 7.204724e-14
ed_power_1_62 = 1.137165e+13
ed_time_1_63 = 7.322835e-14
ed_power_1_63 = 1.156882e+13
ed_time_1_64 = 7.440945e-14
ed_power_1_64 = 1.166868e+13
ed_time_1_65 = 7.559055e-14
ed_power_1_65 = 1.166868e+13
ed_time_1_66 = 7.677165e-14
ed_power_1_66 = 1.156882e+13
ed_time_1_67 = 7.795276e-14
ed_power_1_67 = 1.137165e+13
ed_time_1_68 = 7.913386e-14
ed_power_1_68 = 1.108218e+13
ed_time_1_69 = 8.031496e-14
ed_power_1_69 = 1.070764e+13
ed_time_1_70 = 8.149606e-14
ed_power_1_70 = 1.025723e+13
ed_time_1_71 = 8.267717e-14
ed_power_1_71 = 9.741669e+12
ed_time_1_72 = 8.385827e-14
ed_power_1_72 = 9.172842e+12
ed_time_1_73 = 8.503937e-14
ed_power_1_73 = 8.563310e+12
ed_time_1_74 = 8.622047e-14
ed_power_1_74 = 7.925865e+12
ed_time_1_75 = 8.740157e-14
ed_power_1_75 = 7.273088e+12
ed_time_1_76 = 8.858268e-14
ed_power_1_76 = 6.616957e+12
ed_time_1_77 = 8.976378e-14
ed_power_1_77 = 5.968496e+12
ed_time_1_78 = 9.094488e-14
ed_power_1_78 = 5.337511e+12
ed_time_1_79 = 9.212598e-14
ed_power_1_79 = 4.732383e+12
ed_time_1_80 = 9.330709e-14
ed_power_1_80 = 4.159951e+12
ed_time_1_81 = 9.448819e-14
ed_power_1_81 = 3.625465e+12
ed_time_1_82 = 9.566929e-14
ed_power_1_82 = 3.132611e+12
ed_time_1_83 = 9.685039e-14
ed_power_1_83 = 2.683592e+12
ed_time_1_84 = 9.803150e-14
ed_power_1_84 = 2.279259e+12
ed_time_1_85 = 9.921260e-14
ed_power_1_85 = 1.919279e+12
ed_time_1_86 = 1.003937e-13
ed_power_1_86 = 1.602322e+12
ed_time_1_87 = 1.015748e-13
ed_power_1_87 = 1.326260e+12
ed_time_1_88 = 1.027559e-13
ed_power_1_88 = 1.088365e+12
ed_time_1_89 = 1.039370e-13
ed_power_1_89 = 8.854986e+11
ed_time_1_90 = 1.051181e-13
ed_power_1_90 = 7.142799e+11
ed_time_1_91 = 1.062992e-13
ed_power_1_91 = 5.712368e+11
ed_time_1_92 = 1.074803e-13
ed_power_1_92 = 4.529300e+11
ed_time_1_93 = 1.086614e-13
ed_power_1_93 = 3.560519e+11
ed_time_1_94 = 1.098425e-13
ed_power_1_94 = 2.774998e+11
ed_time_1_95 = 1.110236e-13
ed_power_1_95 = 2.144269e+11
ed_time_1_96 = 1.122047e-13
ed_power_1_96 = 1.642718e+11
ed_time_1_97 = 1.133858e-13
ed_power_1_97 = 1.247712e+11
ed_time_1_98 = 1.145669e-13
ed_power_1_98 = 9.395772e+10
ed_time_1_99 = 1.157480e-13
ed_power_1_99 = 7.014843e+10
ed_time_1_100 = 1.169291e-13
ed_power_1_100 = 5.192431e+10
ed_time_1_101 = 1.181102e-13
ed_power_1_101 = 3.810576e+10
ed_time_1_102 = 1.192913e-13
ed_power_1_102 = 2.772540e+10
ed_time_1_103 = 1.204724e-13
ed_power_1_103 = 2.000010e+10
ed_time_1_104 = 1.216535e-13
ed_power_1_104 = 1.430388e+10
ed_time_1_105 = 1.228346e-13
ed_power_1_105 = 1.014244e+10
ed_time_1_106 = 1.240157e-13
ed_power_1_106 = 7.130150e+09
ed_time_1_107 = 1.251969e-13
ed_power_1_107 = 4.969606e+09
ed_time_1_108 = 1.263780e-13
ed_power_1_108 = 3.434096e+09
ed_time_1_109 = 1.275591e-13
ed_power_1_109 = 2.352719e+09
ed_time_1_110 = 1.287402e-13
ed_power_1_110 = 1.598067e+09
ed_time_1_111 = 1.299213e-13
ed_power_1_111 = 1.076186e+09
ed_time_1_112 = 1.311024e-13
ed_power_1_112 = 7.185327e+08
ed_time_1_113 = 1.322835e-13
ed_power_1_113 = 4.756343e+08
ed_time_1_114 = 1.334646e-13
ed_power_1_114 = 3.121526e+08
ed_time_1_115 = 1.346457e-13
ed_power_1_115 = 2.031084e+08
ed_time_1_116 = 1.358268e-13
ed_power_1_116 = 1.310256e+08
ed_time_1_117 = 1.370079e-13
ed_power_1_117 = 8.380148e+07
ed_time_1_118 = 1.381890e-13
ed_power_1_118 = 5.313912e+07
ed_time_1_119 = 1.393701e-13
ed_power_1_119 = 3.340752e+07
ed_time_1_120 = 1.405512e-13
ed_power_1_120 = 2.082291e+07
ed_time_1_121 = 1.417323e-13
ed_power_1_121 = 1.286784e+07
ed_time_1_122 = 1.429134e-13
ed_power_1_122 = 7.883832e+06
ed_time_1_123 = 1.440945e-13
ed_power_1_123 = 4.788904e+06
ed_time_1_124 = 1.452756e-13
ed_power_1_124 = 2.884046e+06
ed_time_1_125 = 1.464567e-13
ed_power_1_125 = 1.722009e+06
ed_time_1_126 = 1.476378e-13
ed_power_1_126 = 1.019379e+06
ed_time_1_127 = 1.488189e-13
ed_power_1_127 = 5.982787e+05
ed_time_1_128 = 1.500000e-13
ed_power_1_128 = 0.000000e+00
### SETUP LASER BEAM ###
ed_numberOfBeams = 1
# Setup Gaussian Beam:
ed_lensX_1 = -5.1e-03
ed_lensY_1 = 0.0
ed_lensSemiAxisMajor_1 = 3.0e-04
ed_targetX_1 = 0.0
ed_targetY_1 = 0.0
ed_targetSemiAxisMajor_1 = 3.0e-04
ed_pulseNumber_1 = 1
ed_wavelength_1 = 0.8
ed_crossSectionFunctionType_1 = "gaussian1D"
ed_gaussianExponent_1 = 1.0
ed_gaussianRadiusMajor_1 = 1.501403e-04
ed_numberOfRays_1 = 1024
ed_gridType_1 = "regular1D"
#################################
# #
# CONDUCTION PARAMETERS #
# #
#################################
useDiffuse = .true.
useConductivity = .true.
diff_useEleCond = .true.
diff_eleFlMode = "fl_larsen"
diff_eleFlCoef = 0.06
diff_thetaImplct = 1.0
diff_eleXlBoundaryType = "outflow"
diff_eleXrBoundaryType = "outflow"
diff_eleYlBoundaryType = "outflow"
diff_eleYrBoundaryType = "outflow"
diff_eleZlBoundaryType = "outflow"
diff_eleZrBoundaryType = "outflow"
####################################
# #
# HEAT EXCHANGE PARAMETERS #
# #
####################################
useHeatexchange = .true.
##########################
# #
# EOS PARAMETERS #
# #
##########################
eosModeInit = "dens_temp_gather"
smallt = 1.0
smallx = 1.0e-99
eos_useLogTables = .false.
############################
# #
# HYDRO PARAMETERS #
# #
############################
useHydro = .true.
order = 3 # Interpolation order (first/second/third/fifth order)
slopeLimiter = "minmod" # Slope limiters (minmod, mc, vanLeer, hybrid, limited)
LimitedSlopeBeta = 1. # Slope parameter for the "limited" slope by Toro
charLimiting = .true. # Characteristic limiting vs. Primitive limiting
use_avisc = .true. # use artificial viscosity (originally for PPM)
cvisc = 0.1 # coefficient for artificial viscosity
use_flattening = .false. # use flattening (dissipative) (originally for PPM)
use_steepening = .false. # use contact steepening (originally for PPM)
use_upwindTVD = .false. # use upwind biased TVD slope for PPM (need nguard=6)
RiemannSolver = "hllc" # Roe, HLL, HLLC, LLF, Marquina, hybrid
entropy = .false. # Entropy fix for the Roe solver
shockDetect = .false. # Shock Detect for numerical stability
use_hybridOrder = .true. # Enforce Riemann density jump
# Hydro boundary conditions:
xl_boundary_type = "outflow"
xr_boundary_type = "outflow"
yl_boundary_type = "outflow"
yr_boundary_type = "outflow"
zl_boundary_type = "outflow"
zr_boundary_type = "outflow"
##############################
# #
# INITIAL CONDITIONS #
# #
##############################
sim_targetRadius = 2.2e-0.3
sim_targetHeight = 1
sim_vacuumHeight = 0
# Target material defaults set for H at Cryo temperature:
sim_rhoTarg = 2.7
sim_teleTarg = 5.0
sim_tionTarg = 5.0
sim_tradTarg = 5.0
ms_targA = 1.00784
ms_targZ = 1.0
ms_targZMin = 0.1
eos_targEosType = "eos_tab"
eos_targSubType = "ionmix4"
eos_targTableFile = "h-imx-1grp.cn4"
# Chamber material defaults set for Helium at pressure 1.6 mbar:
sim_rhoCham = 1.0e-06
sim_teleCham = 5.0
sim_tionCham = 5.0
sim_tradCham = 5.0
ms_chamA = 4.002602
ms_chamZ = 2.0
eos_chamEosType = "eos_tab"
eos_chamSubType = "ionmix4"
eos_chamTableFile = "he-imx-005.cn4"
###########################
# #
# TIME PARAMETERS #
# #
###########################
tstep_change_factor = 1.10
cfl = 0.4
dt_diff_factor = 1.0e+100 # Disable diffusion dt
rt_dtFactor = 0.2
hx_dtFactor = 1.0e+100
tmax = 10.0e-12
dtmin = 1.0e-20
dtinit = 1.0e-16
dtmax = 3.0e-13
nend = 10000000
###########################
# #
# MESH PARAMETERS #
# #
###########################
geometry = "cartesian"
# Domain size:
xmin = -5.0e-03
xmax = 5.0e-03
ymin = -5.0e-03
ymax = 5.0e-03
# Total number of blocks:
nblockx = 16
nblocky = 16
lrefine_max = 3
lrefine_min = 1
refine_var_1 = "dens"
refine_var_2 = "tele"
# Advanced Parameters
hy_eosModeAfter = "dens_ie_gather"
hy_3TMode = "ragelike"
dr_posdefVar_1 = "eion"
dr_usePosdefComputeDt = TRUE
dr_posdefVar_2 = "eele"
dr_posdefVar_3 = "tele"
dr_posdefVar_4 = "tion"
nblockx = 1
nblocky = 2
lrefine_max = 4
lrefine_min = 1
refine_var_1 = "dens"
refine_var_2 = "tele"
More information about the flash-users
mailing list