#	Configuration file for the Sod shock-tube problem
#	(Sod, G. A., 1978, J. Comp. Phys., 27, 1)

REQUIRES Driver
REQUIRES physics/Hydro

USESETUPVARS ThreeT

IF ThreeT
   # If this is a 3T simulation (i.e. includes the +3t setup shortcut)
   # then use the multitemp gamma law EOS
   REQUIRES physics/Eos/EosMain/multiTemp/Gamma
   REQUESTS physics/Diffuse/DiffuseMain/Split
   REQUESTS physics/Hydro/HydroMain/split/PPM/multiTemp
	
   REQUIRES physics/sourceTerms/Heatexchange/HeatexchangeMain/Spitzer
   REQUIRES physics/materialProperties/Conductivity/ConductivityMain/SpitzerHighZ
   REQUIRES physics/materialProperties/Viscosity/ViscosityMain/Spitzer/
ELSE
   REQUIRES physics/Eos
ENDIF

USESETUPVARS withParticles

IF withParticles
   PARTICLETYPE passive INITMETHOD lattice MAPMETHOD quadratic ADVMETHOD rungekutta
   REQUIRES Particles/ParticlesMain
   REQUESTS Particles/ParticlesMapping/Quadratic
   REQUESTS Particles/ParticlesInitialization/Lattice
   REQUESTS IO/IOMain/
   REQUESTS IO/IOParticles
ENDIF

D sim_confGeometry Geometry of the physical problem configuration (initial condition).
D & If empty, "none", or "default", use the Grid (coordinate) geometry for the physical geometry.
PARAMETER sim_confGeometry STRING "none" [ cartesian, polar ,cylindrical,spherical,"","none","NONE","default","DEFAULT" ]

D sim_rhoLeft Density in the left part of the grid 
PARAMETER sim_rhoLeft	REAL    1.       [0 to ]

D sim_rhoRight Density in the right part of the grid 
PARAMETER sim_rhoRight	REAL    0.125    [0 to ]

D sim_pLeft Pressure  in the left part of the grid 
PARAMETER sim_pLeft	REAL    1.       [0 to ]

D sim_pRight Pressure  in the righ part of the grid 
PARAMETER sim_pRight	REAL    0.1      [0 to ]

D sim_uLeft fluid velocity in the left part of the grid 
PARAMETER sim_uLeft	REAL    0.       

D sim_uRight fluid velocity in the right part of the grid 
PARAMETER sim_uRight	REAL    0.       

D sim_xangle  Angle made by diaphragm normal w/x-axis (deg)
PARAMETER sim_xangle	REAL	0.       [0 to 360]

D sim_ yangle Angle made by diaphragm normal w/y-axis (deg)
PARAMETER sim_yangle	REAL	90.      [0 to 360]

D sim_posnR  Point of intersection between the shock plane and the x-axis
PARAMETER sim_posn	    REAL    0.5

PARAMETER sim_sigma	    REAL    0.5

# For a variant of the Sod simulation with two materials, add TwoMaterials=True
# to the setup command line.
USESETUPVARS TwoMaterials
IF TwoMaterials
   PPDEFINE SIMULATION_TWO_MATERIALS
   SPECIES LEFT
   SPECIES RGHT

   D sim_abarLeft  ion mean molecular weight of material on left
   PARAMETER sim_abarLeft   REAL   1.        [0 to ]
   D sim_zbarLeft  ion average charge for material on left
   PARAMETER sim_zbarLeft   REAL   1.        [0 to ]
   D sim_abarRight ion mean molecular weight of material on right
   PARAMETER sim_abarRight  REAL   1.        [0 to ]
   D sim_zbarRight ion average charge for material on right
   PARAMETER sim_zbarRight  REAL   1.        [0 to ]
ENDIF

#SPECIES H1

# Add some options for 3T simulations:
IF ThreeT
   D gamma       Ratio of specific heats for gas - for initialization
   PARAMETER gamma         REAL    1.6666666666666667  [0.0 to ] 

   D gammaIon       Ratio of specific heats for ion component, should be 5./3.
   PARAMETER gammaIon         REAL    1.6666666666666667 [0.0 to ]

   # This line is needed for entropy advection to work:
   MASS_SCALAR sele EOSMAP: SELE

   PARAMETER sim_pionLeft  REAL -1.0
   PARAMETER sim_pionRight REAL -1.0
   PARAMETER sim_peleLeft  REAL -1.0
   PARAMETER sim_peleRight REAL -1.0
   PARAMETER sim_pradLeft  REAL -1.0
   PARAMETER sim_pradRight REAL -1.0
ENDIF

USESETUPVARS AddtlVelVars
IF AddtlVelVars
   LINKIF Grid_computeUserVars.F90  Grid
   D vlcx_variable velocity componenent in the X-direction of sim_confGeometry
   D vlcy_variable velocity componenent in the Y-direction of sim_confGeometry
   D vlcz_variable velocity componenent in the Z-direction of sim_confGeometry
   VARIABLE vlcx	TYPE: PER_MASS
   VARIABLE vlcy	TYPE: PER_MASS
   VARIABLE vlcz	TYPE: PER_MASS
ENDIF
