| 1 | !================================================================================
|
|---|
| 2 | ! AstroBEAR physics.data file.
|
|---|
| 3 | ! Input data describing built-in physics-related parameters of the code.
|
|---|
| 4 | ! Created for problem: Template
|
|---|
| 5 | ! Default values for variables in [brackets]
|
|---|
| 6 | !================================================================================
|
|---|
| 7 |
|
|---|
| 8 | &EOSData
|
|---|
| 9 | !================================================================================
|
|---|
| 10 | ! Equation of State stuff
|
|---|
| 11 | !================================================================================
|
|---|
| 12 | iEOS = 1 ! Equation of State [0]-ideal gas, 1-MultiSpecies gas, 3-Van der Waals gas, 4-Isothermal
|
|---|
| 13 | gamma = 1.666667d0 ! Adiabatic index used for ideal gas and Van der Waals gas [1.66666666667]
|
|---|
| 14 | IsoTemp = 10d0 ! Isothermal temperature for Isothermal EOS
|
|---|
| 15 | CriticalDens = 30.693677103 ! maximum density in Van der Waals EOS
|
|---|
| 16 | AttractionParam = 1.38d6 ! attraction pram in Van der Waals EOS
|
|---|
| 17 | /
|
|---|
| 18 | &AbundanceData
|
|---|
| 19 | HeFrac = 0.085d0
|
|---|
| 20 | xH = 0.01d0
|
|---|
| 21 | xHe = 0d0
|
|---|
| 22 | xHeII = 0d0
|
|---|
| 23 | ZFrac = 0.001165d0
|
|---|
| 24 | Zmu = 16.48d0
|
|---|
| 25 | Xmu = 0d0
|
|---|
| 26 | xZ = 4.13127d-4
|
|---|
| 27 | /
|
|---|
| 28 |
|
|---|
| 29 | &PhysicsData
|
|---|
| 30 |
|
|---|
| 31 | !================================================================================
|
|---|
| 32 | ! Field based refinement control
|
|---|
| 33 | !================================================================================
|
|---|
| 34 | InterpOpts = 1,1,1,1,1,0,0,0,1,1,1,1,1 ! [0]-constant, 1-minmod, 2-superbee, 3-vanLeer, 4-mc, 5-Parabolic (not conservative), 6-Linear (not limited). It is recommended to use constant interpolation for any fields in q that correspond to aux fields (iBx, iBy, iBz, etc...)
|
|---|
| 35 | refineVariableFactor = 1d0,1d0,1d0,0d0,1d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0,0d0 ! weight factors for each field used in triggering refinement. Combined with qtolerance for that level. [0d0]
|
|---|
| 36 |
|
|---|
| 37 | !================================================================================
|
|---|
| 38 | ! Global source term switches
|
|---|
| 39 | !================================================================================
|
|---|
| 40 | lSelfGravity = .false. ! Turns on self-gravity if true [.false.]
|
|---|
| 41 | UniformGravity = 0d0 ! Gravitational acceleration in the -y direction [0d0]
|
|---|
| 42 | OmegaRot = 0d0 ! Shifts to a rotating reference frame (enables coriolis and centripetal force terms)
|
|---|
| 43 | iCylindrical = 2 ! [0]-No cylindrical geometry, 1-Cylindrical with no angular momentum, 2-Cylindrical with angular momentum
|
|---|
| 44 | iCooling = 4 ! Cooling module [NoCool=0, AnalyticCool=1, DMCool=2, IICool=3, ZCool=4]
|
|---|
| 45 | SrcPrecision = 1d-3 ! Relative precision desired for source terms. [1d-3]
|
|---|
| 46 | lTrackHydrogen = .true.
|
|---|
| 47 | lTrackHelium = .true.
|
|---|
| 48 | ne_min=1d-3
|
|---|
| 49 | FloorTemp=1d3
|
|---|
| 50 | !================================================================================
|
|---|
| 51 | ! Explicit solver switches
|
|---|
| 52 | !================================================================================
|
|---|
| 53 | lResistive = .false. ! Turns on resistivity [.false.]
|
|---|
| 54 | ResType = 2 ! Resistivity type. 1-constant, 2-user defined distribution, 3-Spitzer. [1]
|
|---|
| 55 | resistivity = 1.2d-1 ! resistivity value [0d0]
|
|---|
| 56 | lConductive = .false. ! Turns on mhd heat conduction [.false.]
|
|---|
| 57 | ConType = 1 ! Conduction type. 1-constant, 2-Spitzer. [1]
|
|---|
| 58 | conductivity = 8d-2 ! conductivity value [0d0]
|
|---|
| 59 | lThermalFluxLimiter= .false. ! Turns on saturation heat flux. (only works for ConType==2). [.false.]
|
|---|
| 60 | lIsoDiff = .false. ! Turns on isotropic heat conduction. Allows isotropicity to be in the range (0,1). [.false.]
|
|---|
| 61 | lViscous = .false. ! Turns on viscosity (2D only) [.false.]
|
|---|
| 62 | VisType = 1 ! viscosity type. 1-classic, 2-Braginskii [1]
|
|---|
| 63 | bviscosity = 0d0 ! viscosity value [0d0]
|
|---|
| 64 |
|
|---|
| 65 | !================================================================================
|
|---|
| 66 | ! Density Protection Options
|
|---|
| 67 | !================================================================================
|
|---|
| 68 | lRestartOnDensityProtections = .false. ! Do density protections trigger restarts? [.false.]
|
|---|
| 69 | iDensityProtect = 2 ! 0-Set to MinDensity, 1-Set to minimum nearby density, [2]-Average nearby densities
|
|---|
| 70 | iMomentumProtect = 2 ! 0-Conserve momentum, 1-Set to zero, [2]-Average nearby velocities
|
|---|
| 71 | MinDensity = 1d-4 ! Minimum computational density before protection is triggered [1d-10]
|
|---|
| 72 |
|
|---|
| 73 | !================================================================================
|
|---|
| 74 | ! Pressure Protection Options
|
|---|
| 75 | !================================================================================
|
|---|
| 76 | lRestartOnPressureProtections = .false. ! Do pressure protections trigger restarts? [.false.]
|
|---|
| 77 | iPressureProtect = 2 ! 0-Set to MinTemp, 1-Set to minimim nearby pressure, [2]-Set to average nearby pressure, 3-Set to minimum nearby temperature, 4-Set to average nearby temperature
|
|---|
| 78 | MinTemp = 1d-10 ! [1d-10] minimum allowed temperature for the system in Kelvin before protection is triggered
|
|---|
| 79 |
|
|---|
| 80 | !================================================================================
|
|---|
| 81 | ! Other Protection Options
|
|---|
| 82 | !================================================================================
|
|---|
| 83 | lRestartOnHypreMaxIters = T ! Restart steps if hypre does not converge
|
|---|
| 84 |
|
|---|
| 85 | !================================================================================
|
|---|
| 86 | ! Description of various scaling parameters
|
|---|
| 87 | ! Define one of each of the following: [nScale/rScale], [TempScale/pScale], and set the other to 0d0.
|
|---|
| 88 | !================================================================================
|
|---|
| 89 | nScale = 1d2, ! number density scale parameter [particles/cc]
|
|---|
| 90 | rScale = 0, ! density scale [g/cc]
|
|---|
| 91 | TempScale = 1d4, ! temperature scale parameter [Kelvin]
|
|---|
| 92 | pScale = 0d0, ! pressure scale [dynes/cm^2]
|
|---|
| 93 | lScale = 5e15, ! length scale parameter [cm] (AU=1.49598e13, pc=3.08568025e18, R_sun=6.955e10
|
|---|
| 94 |
|
|---|
| 95 | !================================================================================
|
|---|
| 96 | ! MHD related section
|
|---|
| 97 | !================================================================================
|
|---|
| 98 | lMHD = .false. ! Magnetic Fields present? [.false.]
|
|---|
| 99 | lCheckDivergence = .false. ! Turn on divergence checking [.false.]
|
|---|
| 100 | /
|
|---|
