run_comment = "Laser Slab Example Simulation"
log_file    = "lasslab.log"
basenm      = "lasslab_"

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

##########################
#                        #
#     I/O PARAMETERS     #
#                        #
##########################

### Checkpoint Options  ###
checkpointFileIntervalTime = 2.0
checkpointFileIntervalStep = 1000

### Plot Options ###
plotFileNumber       = 0
plotFileIntervalStep = 20
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"

### 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"

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 = "al-imx-003.cn4"


############################
#                          #
#     LASER PARAMETERS     #
#                          #
############################
useEnergyDeposition = .true.
ed_maxRayCount      = 10000
ed_gradOrder        = 2

### LASER IO OPTIONS ###
ed_useLaserIO                  = .false.
ed_laserIOMaxNumberOfPositions = 10000
ed_laserIOMaxNumberOfRays      = 128

### SETUP LASER PULSES ###
ed_numberOfPulses = 1

# Define Pulse 1:
ed_numberOfSections_1 = 25

ed_time_1_1 = 4e-12
ed_time_1_2 = 8e-12
ed_time_1_3 = 1.2e-11
ed_time_1_4 = 1.6e-11
ed_time_1_5 = 2e-11
ed_time_1_6 = 2.4e-11
ed_time_1_7 = 2.8e-11
ed_time_1_8 = 3.2e-11
ed_time_1_9 = 3.6e-11
ed_time_1_10 = 4e-11
ed_time_1_11 = 4.3999999999999997e-11
ed_time_1_12 = 4.8e-11
ed_time_1_13 = 5.2e-11
ed_time_1_14 = 5.6e-11
ed_time_1_15 = 6e-11
ed_time_1_16 = 6.4e-11
ed_time_1_17 = 6.8e-11
ed_time_1_18 = 7.2e-11
ed_time_1_19 = 7.6e-11
ed_time_1_20 = 8e-11
ed_time_1_21 = 8.4e-11
ed_time_1_22 = 8.799999999999999e-11
ed_time_1_23 = 9.199999999999999e-11
ed_time_1_24 = 9.6e-11
ed_time_1_25 = 1e-10
ed_power_1_1 = 15046970018.468273
ed_power_1_2 = 27692071980.63451
ed_power_1_3 = 48214697723.15642
ed_power_1_4 = 79418383890.12585
ed_power_1_5 = 123759983402.13225
ed_power_1_6 = 182455533862.22116
ed_power_1_7 = 254478700244.9149
ed_power_1_8 = 335786658096.36926
ed_power_1_9 = 419172743215.27844
ed_power_1_10 = 495039933608.52893
ed_power_1_11 = 553101750652.935
ed_power_1_12 = 584638528706.2346
ed_power_1_13 = 584638528706.2346
ed_power_1_14 = 553101750652.9352
ed_power_1_15 = 495039933608.529
ed_power_1_16 = 419172743215.2786
ed_power_1_17 = 335786658096.3693
ed_power_1_18 = 254478700244.91498
ed_power_1_19 = 182455533862.22128
ed_power_1_20 = 123759983402.13232
ed_power_1_21 = 79418383890.12592
ed_power_1_22 = 48214697723.15646
ed_power_1_23 = 27692071980.63457
ed_power_1_24 = 15046970018.468273
ed_power_1_25 = 7734998962.633265


### SETUP LASER BEAM ###
ed_numberOfBeams = 1

ed_lensX_1                    = -1000.0e-07
ed_targetX_1                  =  6000.0e-07
ed_pulseNumber_1              =  1
ed_wavelength_1               =  1.053
ed_crossSectionFunctionType_1 = "uniform"
ed_numberOfRays_1             =  1
ed_gridType_1                 = "regular1D"
ed_gridnRadialTics_1          =  512


#################################
#                               #
#     CONDUCTION PARAMETERS     #
#                               #
#################################
useDiffuse      = .true.
useConductivity = .true.
diff_useEleCond = .true.
diff_eleFlMode  = "fl_larsen"
diff_eleFlCoef  = 0.06
diff_thetaImplct = 1.0

diff_eleXlBoundaryType = "neumann"
diff_eleXrBoundaryType = "neumann"


####################################
#                                  #
#     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"


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

sim_targetRadius = 200.0e-07
sim_targetHeight = 100.0e-07
sim_vacuumHeight = 6000.0e-07

# Target material defaults set for Aluminum at room temperature:
sim_rhoTarg  = 2.7
sim_teleTarg = 0.00005
sim_tionTarg = 0.00005
sim_tradTarg = 0.00005
ms_targA = 26.9815386
ms_targZ = 13.0
ms_targZMin = 0.02
eos_targEosType = "eos_tab"
eos_targSubType = "ionmix4"
eos_targTableFile = "al-imx-003.cn4"

# Chamber material defaults set for Helium at pressure 1.6 mbar:
sim_rhoCham  = 1.0e-06
sim_teleCham = 0.00005
sim_tionCham = 0.00005
sim_tradCham = 0.00005
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           = 150.0e-12
dtmin          = 1.0e-16
dtinit         = 1.0e-15
dtmax          = 20.0e-9
nend           = 10000000


###########################
#                         #
#     MESH PARAMETERS     #
#                         #
###########################
geometry = "cartesian"

# Domain size:
xmin = 0.0
xmax = 6500.0e-7

# Total number of blocks:
nblockx = 3

lrefine_max = 4
lrefine_min = 1
refine_var_1 = "dens"
refine_var_2 = "tele"
refine_var_3 = "bdry"