fmtout |
fmtout |
|
|
False |
|
|
maxstep |
maximum number of time steps
- for a normal run (e.g., 48 hrs), this is a large integer, e.g., maxstep = 10000000
- for testing/debugging it can be set to a smaller number, e.g., maxstep = 10 |
|
|
1000000 |
|
|
hrmax |
runtime (hrs)
- maximum runtime is hrmax = 168 (7 days)
- nominal minimum runtime is hrmax = 36
(it takes ∼ 12 hrs for the ionosphere to develop) |
|
|
48. |
|
|
dthr |
output cadence (hrs)
- nominal output cadence is dthr = 0.25 (15 min)
- for testing/debugging it can be reduced to output data more often |
|
hours |
0.25 |
|
|
hrpr |
time data starts to be output
- nominally hrpr = 0 and data is output from beginning of the run
- alternatively, if hrmax = 48 and hrpr = 24, data is output starting after the first day
(this reduces the size of the output files) |
|
hours |
24. |
|
|
dt0 |
maximum time step (sec)
- should be in the range 2 - 8 sec
- the larger dt0, the faster the code runs in general
- however, depending on geophysical conditions and the grid, the code can give un- physical
results if initially too large (e.g., dt0 = 8) |
|
seconds |
30. |
|
|
grad_in |
grad_in |
|
|
300. |
|
|
glat_in |
glat_in |
|
|
0. |
|
|
glon_in |
glon_in |
|
|
0. |
|
|
fejer |
False for Sinusoidal Model and True for Fejer/ Scherliess Empirical Model |
|
|
False |
|
|
rmin |
used to specify grid along the magnetic field
- best to leave rmin = 90 |
|
|
150 |
|
|
rmax |
rmax |
|
|
20000 |
|
|
altmin |
used to specify grid along the magnetic field
- best to leave altmin = 85 km |
|
|
85 |
|
|
f10.7_avg |
F10.7A
- 81 day average of F10.7 |
|
|
70 |
|
|
f10.7 |
F10.7
- solar radio flux at 10.7 cm
- used as a measure of solar activity |
|
|
70 |
|
|
Ap |
Ap index
- The Ap index is a daily measure of magnetic activity |
|
|
6 |
|
|
year |
Year of the event |
|
|
2019 |
|
|
day |
Day-of-the-year (DOY) of the event |
|
|
250 |
|
|
mmass |
mass number used in nrlmsise00
- mmass = 48 (do not change) |
|
|
48 |
|
|
nion1 |
first ion specie
- nion1 = 1 (H+) (do not change) |
|
|
1 |
|
|
nion2 |
last ion specie
- nion1 = 7 (N+) (do not change) |
|
|
7 |
|
|
hrinit |
starting time (UT hrs)
- nominal value hrinit = 0 (i.e., midnight GMT) - best not to change |
|
UT hours |
0. |
|
|
tvn0 |
modifies the magnitude of the neutral wind speed
- nominally tvn0 = 1
- use caution when changing, e.g., tvn0 = 2 will double the neutral wind speed |
|
|
1 |
|
|
tvexb0 |
modifies the magnitude of the E × B velocity
- nominally tvexb0 = 1
- use caution when changing, e.g., tvexb0 = 0 will set the E × B drift to 0 |
|
|
1. |
|
|
ver |
Max ExB |
|
|
15 |
|
|
veh |
not used (leave 0) |
|
|
|
|
|
gams |
adjusts grid spacing along the magnetic field line
- best to leave gams = .1 |
|
|
3 |
|
|
gamp |
adjusts grid spacing in latitude
- best to leave gams = 2.5 |
|
|
3 |
|
|
snn |
multiplier for neutral densities
- the 7 entries correspond to the 7 neutral species in parameter mod.f90
- for example, if snn = 1.,2.,1.,1.,1.,1.,1., then the oxygen neutral density is doubled |
|
|
1.,1.,1.,1.,1.,1.,1. |
|
|
stn |
multiplier for the neutral temperature
- this will change the ion temperatures because of ion-neutral coupling - it will not change the scale height of the neutrals |
|
|
1. |
|
|
denmin |
minimum ion density allowed (cm−3)
- best to leave denmin = 1.e-6 |
|
|
1.e-6 |
|
|
alt_crit |
alt_crit |
|
|
150. |
|
|
cqe |
cqe |
|
|
6e-14 |
|
|