run_comment = "Simulations of OSU Experiments"
log_file    = "lasslab.log"
basenm      = "lasslab_"

# This particular parfile is used as an example parfile that is
# described in detail in the users guide.

#########################
#   OUTPUT PARAMETERS   #
#########################

### Checkpoint Options  ###
checkpointFileIntervalTime =  1.0
checkpointFileIntervalStep = 10

### Plot Options ###
plotFileNumber       = 0
plotFileIntervalStep = 10
plot_var_1           = "dens"
plot_var_2           = "trag"
plot_var_3           = "tele"
plot_var_4           = "tion"
plot_var_5           = "cham "
plot_var_6           = "foil"

### 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 = "reflecting"
rt_mgdXrBoundaryType = "vacuum"
rt_mgdYlBoundaryType = "vacuum"
rt_mgdYrBoundaryType = "reflecting"
rt_mgdZlBoundaryType = "reflecting"
rt_mgdZrBoundaryType = "reflecting"

useOpacity     = .false.

### SET Target (CH) OPACITY OPTIONS ###
op_targAbsorb   = "op_tabpa"
op_targEmiss    = "op_tabpe"
op_targTrans    = "op_tabro"
op_targFileType = "ionmix4"
op_targFileName = "polystyrene-imx-001.cn4"

### SET Foil (Cu) OPACITY OPTIONS ###
op_FoilAbsorb   = "op_tabpa"
op_FoilEmiss    = "op_tabpe"
op_FoilTrans    = "op_tabro"
op_FoilFileType = "ionmix4"
op_FoilFileName = "al-imx-003.cn4"

### SET Chamber (He) OPACITY OPTIONS ###
op_chamAbsorb   = "op_tabpa"
op_chamEmiss    = "op_tabpe"
op_chamTrans    = "op_tabro"
op_chamFileType = "ionmix4"
op_chamFileName = "he-imx-005.cn4"


### 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 = "vacuum"
diff_eleYrBoundaryType = "outflow"
diff_eleZlBoundaryType = "outflow"
diff_eleZrBoundaryType = "outflow"

### Heat-Exchange parameters ###
useHeatexchange = .true.

### Material and EOS parameters ###
eosModeInit      = "dens_temp_gather"
eosMode          = "dens_ie_recal_gather"
smallt           = 1.0
smallx           = 1.0e-99
eos_useLogTables = .false.

############################
#                          #
#     HYDRO PARAMETERS     #
#                          #
############################

# Use second order hybrid solver with minmod slope limiter. This
# essentially eliminates any Carbuncle instability.
 
order            = 2        # Interpolation order (first/second/third/fifth order)
transOrder       = 1
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    = "hybrid" # Roe, HLL, HLLC, LLF, Marquina, hybrid
entropy          = .false.  # Entropy fix for the Roe solver
shockDetect      = .false.  # Shock Detect for numerical stability


# Hydro boundary conditions:
xl_boundary_type = "periodic"
xr_boundary_type = "periodic"
yl_boundary_type = "reflect"
yr_boundary_type = "reflect"
zl_boundary_type = "reflect"
zr_boundary_type = "reflect"


##########################
#   INITIAL CONDITIONS   #
##########################

sim_targetRadius = 200.0e-04
sim_targetHeight = 27.0e-04
sim_interfaceght = 2.0e-04
sim_vacuumHeight = 8.0e-04


sim_amplitudeght = 3.0e-05
sim_lambdaX = 30.e-4

# Foil material defaults set for Aluminum at 25 eV:
sim_rhoFoil  = 8.93
sim_teleFoil = 290000.0
sim_tionFoil = 290000.0
sim_tradFoil = 290000.0
ms_FoilA = 63.5
ms_FoilZ = 29.0 
sim_zminFoil =0.02
eos_FoilEosType = "eos_tab"
eos_FoilSubType = "ionmix4"
eos_FoilTableFile = "al-imx-003.cn4"

# Target material defaults set for Helium at 400 eV:
sim_rhoTarg  = 1.1
sim_teleTarg = 4640000.0
sim_tionTarg = 4640000.0
sim_tradTarg = 4640000.0
ms_TargA = 12.01
ms_TargZ = 6.0
sim_zminTarg =0.08
eos_TargEosType = "eos_tab"
eos_TargSubType = "ionmix4"
eos_TargTableFile = "polystyrene-imx-001.cn4"

# Chamber material defaults set for Helium at room temperature:
sim_rhoCham  = 1.0e-08
sim_teleCham = 290.11375
sim_tionCham = 290.11375
sim_tradCham = 290.11375
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.1
cfl            = 0.4
dt_diff_factor = 1.0e+100 # Disable diffusion dt
rt_dtFactor    = 0.02
hx_dtFactor    = 1.0e+100
tmax           = 1.0e-8
dtmin          = 0.1e-14
dtinit         = 0.1e-14
dtmax          = 1.0e-08
nend           = 8000


#######################
#   MESH PARAMETERS   #
#######################

geometry = "cartesian"

# Domain size:
xmin = 0.0e-04
xmax = 40.0e-04
ymin = -40.0e-04
ymax = 40.0e-04

# Total number of blocks:
nblockx = 1
nblocky = 1

lrefine_max = 9
lrefine_min = 1
refine_var_1 = "dens"
refine_var_3= "pres"