fem.verboserank: 0

# PROBLEM SELECTION 
#------------------
# possible: "sod" , "withman", "withmansmooth", "smooth1d" , "ffs" , "diffractio
problem: sod
# problemflag: 0
# riemanndata: [1., -1., 1.,    1., 1., 1.]

# number of threads used in OMP program
fem.parallel.numberofthreads: 4

# coarsest level that should be present
fem.adaptation.coarsestLevel: 0
# finest level that should be present
fem.adaptation.finestLevel: 4
fem.eoc.steps: 4
fem.timeprovider.factor: 0.35
femdg.limiter.indicator: jump

#femdg.stepper.fixedtimestep: 0.005
#femdg.nonblockingcomm: false

femdg.stepper.printcount: 1000
fem.io.outputformat: none
fem.io.savestep: 10.0

fem.ode.odesolver: EX
fem.ode.order: 3
fem.timeprovider.factor: 0.35
fem.timeprovider.updatestep: 1

fem.io.macroGridFile_2d:
fem.io.macroGridFile_2d:  ../../grids/grid2d_str1d.dgf
#fem.io.macroGridFile_2d:  ../../grids/unitcube2_per.dgf

fem.eoc.outputpath: ./data
#femdg.fixedtimestepeocloopfactor: 1

# GENERAL
#--------
paramfile: ../../parameter/paramBase



# PROBLEM SELECTION 
#------------------
# possible: "sod" , "withman", "withmansmooth", "smooth1d" , "ffs" , "diffraction" , "shockbubble"
problem: sod 
# problemflag: 0
# riemanndata: [1., -1., 1.,    1., 1., 1.]

# LIMITER SETTINGS
#-----------------
# 0 = only dg solution | 1 = only reconstruction | 2 = both 
femdg.limiter.admissiblefunctions: 1 
# tolerance for shock indicator 
femdg.limiter.tolerance: 1
# threshold for avoiding over-excessive limitation 
femdg.limiter.epsilon: 1e-8
# add indicator to outputvariables 
femdg.limiter.indicatoroutput: true

femdg.limiter.indicator: none

# DATA WRITER
#------------
fem.io.datafileprefix: RC
fem.io.savestep: 0.01 # SaveStep (write data every `saveStep' time period, <=0 deactivates)
fem.io.savecount: -1 # SaveCount (write data every saveCount time steps, <=0 deactivates)


# GRID SOLUTION 
#--------------
gridsol.savestep: 0.5
gridsol.firstwrite: 0.5
gridsol.filename: straka-checkpoint

# coarsest level that should be present (also initial refinement)
fem.adaptation.coarsestLevel: 0

# STEPPER
#--------
femdg.stepper.starttime: 0.
femdg.stepper.endtime: 0.15
#femdg.stepper.maxtimestep: 0.1


# PROBLEM SETUP
#--------------
paramfile: ../../parameter/paramPhysicalConstants

# choises are: LLF, HLL, HLLC, LLF2
dgadvectionflux.method: EULER-HLLC

g: 1. # gravity force [m/ss]
hSpeed0: 1. # x-axis speed [m/s]
vSpeed0: 1.25 # z-axis speed [m/s]
p0: 1. # surface pressure [Pa]
N: 0.01 # Brunt-Vaisala frequency
Re: 1. # Reynold's number
Pr: 1. # Prandtl's number
rho_a: 1.
rho_b: 0.05

delta_T: -15 # temperature perturbation peak [K]
x_c: 0.25 # x-center of perturbation ball
y_c: 0.25 # y-center of perturbation ball
z_c: 0.25 # z-center of perturbation ball
r: 0.25 # radius of perturbation ball

# coarsest level that should be present
fem.adaptation.coarsestLevel: 0
# finest level that should be present
fem.adaptation.finestLevel: 6


# DOMAIN SETUP
#-------------
fem.io.macroGridFile_1d:  ../../grids/unitcube1.dgf
fem.io.macroGridFile_2d:  ../../grids/grid2d_str1d.dgf
#fem.io.macroGridFile_2d:  ../../grids/grid2d_nonaffine.dgf
fem.io.macroGridFile_3d:  ../../grids/unitcube3.dgf


# SOLVER CONFIGURATION
#---------------------
fem.ode.odesolver: EX
fem.ode.order: 3
#fem.ode.order: 2
fem.ode.verbose: none # ode output: none, cfl, full
fem.ode.cflincrease: 1.05
fem.ode.miniterations: 95
fem.ode.maxiterations: 105
fem.ode.cflStart: 1.
#fem.ode.cflMax: 5
fem.timeprovider.factor: 0.25
fem.timeprovider.updatestep: 1

# parameter for the implicit solvers 
fem.solver.verbose: false
fem.solver.gmres.restart: 15
fem.solver.newton.verbose: false
fem.solver.newton.linear.verbose: false
fem.solver.newton.maxlineariterations: 1000
fem.solver.newton.tolerance: 1e-10

# GENERAL
#--------
paramfile: ../../parameter/paramBase


paramfile: ../../parameter/paramSolver
