# FLASH3 Configuration file for 1D Coronal Loop with proton-electron
# plasma, conduction, and radiative cooling. No magnetic field
# is included so the ordinary hydro PPM solver is used.
# Uses plane parallel gravity with changes to implement semi-circular
# loop geometry
# Use Ionize without any additional SPECIES to get just a proton+electron
# plasma.  Config in Simulation/SimulationComposition/Ionize contains
# just e and h 

REQUIRES Driver
REQUIRES physics/Hydro
REQUIRES physics/Eos/EosMain/Gamma
REQUIRES physics/sourceTerms/Diffuse/DiffuseMain
REQUIRES physics/sourceTerms/Heat/HeatMain/StatPlusGauss
REQUIRES physics/sourceTerms/Cool/CoolMain/Radloss
REQUIRES physics/materialProperties/Conductivity/ConductivityMain/Spitzer
REQUIRES physics/materialProperties/Viscosity/ViscosityMain/Spitzer
REQUIRES Simulation/SimulationComposition/Ionize
REQUIRES physics/Gravity/GravityMain/PlanePar
REQUIRES IO
REQUIRES Grid/GridBoundaryConditions/OneRow

# Values of the physical parameters are taken from Table 1. of
# Chae & Poland (1997). Assume a plasma upflow and use an average 
# value of .4 km/sec from Table 1 (Chae&Poland). Units are cgs. 
# The coronal input energy is given as a percentage value 
# between 0 and 1.

D sim_plasmaVel   The plasma inflow velocity for solar model, parallel to x-axis
D sim_neCor       Particle density in corona
D sim_TChr        Temperature of chromosphere
D sim_TCor        Temperature of corona
D sim_p0Chr       Pressure at base of chromosphere
D sim_Gsun        Solar gravity at the surface
D sim_Gamma       Ideal gas constant gamma
 
PARAMETER  sim_plasmaVel    REAL    0.0e+00
PARAMETER  sim_neCor        REAL    1.0e+09
PARAMETER  sim_TChr         REAL    8000.0
PARAMETER  sim_TCor         REAL    1.0e+06  
PARAMETER  sim_p0Chr        REAL    0.4 
PARAMETER  sim_Gsun         REAL    2.74e+04
PARAMETER  sim_Gamma        REAL    1.66666667