CMR Header

CCMC Metadata Registry (CMR)

Menu Page

Go to:

CMR: View Simulation Run Info

SAMI3/3.22/ROR/SAMI3-HWM-01_2019-08-TP-01_082823_IT_1

Run Description

CCMC ROR run for SAMI3 3.22 with RORID SAMI3-HWM-01_2019-08-TP-01_082823_IT_1

Simulation Time:

Simulation Start Time: 2019-08-04T00:00:00Z
Simulation End Time: 2019-08-08T00:00:00Z
Time Step in second: 900
Simulation Time Description:
Simulation Time Caveats:

Temporal Dependence (whether the outputs of the simulation code are time-dependent?): true

Regions:

Earth.NearSurface.Ionosphere

Simulation Domain :

Contacts :

Min-Yang.Chou, ModelUser

List of Run Input Parameter Grouping and their associated properties :

Name Description Caveats Properties
All settings for SAMI3 model
All settings used to run SAMI3 model
Name Description Caveats Units Default Value or Actual Value Valid Min Valid Max
vw
constant vertical wind (cm/s)
-best to leave vw=0
cm/s 0
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
imonth_1
month of starting day-of-year
iday_1
day (in month) of starting day-of-year
imonth_2
month of last day-of-year
iday 2
day (in month) of last day-of-year
rmin
used to specify grid along the magnetic field
- best to leave rmin = 90
90.
altmin
used to specify grid along the magnetic field
- best to leave altmin = 85 km
85.
fbar
F10.7A
- 81 day average of F10.7
f10p7
F10.7
- solar radio flux at 10.7 cm
- used as a measure of solar activity
Ap
Ap index
- The Ap index is a daily measure of magnetic activity
iyear
Year of the event
nday
Day-of-the-year (DOY) of the event
- the first value is the starting day-of-year
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
not used (leave 0)
0
veh
not used (leave 0)
0
gams
adjusts grid spacing along the magnetic field line
- best to leave gams = .1
.1
gamp
adjusts grid spacing in latitude
- best to leave gams = 2.5
2.5
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
alt_crit_avg
used to specify values at high latitude boundary
- best to leave alt_crit_avg = 6370.
6370
blat_min
latitude of the base of the first field line
- best to leave blat_min = 6.85
6.85
blat_max
latitude of the base of the higest field line
- nominal values 60-89 (degree)
- maximum value is 89 (degree)
degree 89.
psmooth
smoothing factor (not used)
1.
hall
Hall conductivity
- hall = true (include the Hall and Pedersen conductances in the potential solver)
- hall = false (include only the Pedersen conductance in the potential solver)
- by toggling this switch the influence of the Hall conductance on the electric field is elucidated
- generally, hall = true should be used
true
restart
restart a run
- restart = false (do not restart)
- restart = true (restart a run; needs rst files)
false
storm_ti
for Volland-Stern potential (not used)
storm_tf
for Volland-Stern potential (not used)
vexb_max
maximum E X B velocity (cm/s)
- nominally use 40 km/s; doesn't affect ionosphere dynnamics
cm/s 4000000
lmadala
potential solver.
- lmadala = true (solve for the neutral wind dynamo field)
- lmadala = false (does not solve for the dynamo field: E X B drift is zero)
true
lcr
corotation potential.
- lcr = true (adds the corotation potential)
- lcr = false (does not add the corotation potenial_
- best to leave lcr = false
false
lvs
Volland-Stern high latitude potential.
- best to leave lvs = false
false
lweimer
Weimer high latitude potential.
- lweimer = true (use Weimer high latitude potential); needs OMNI.txt file linked to year of Weimer data file.
- lweimer = false  (do not use Weimer high latitude potential)
true
decay_time
enhanced loss of H+ and He+.
- invoked at distance greater than pcrit . RE where RE is the earth radius
1
pcrit
used for emhanced loss (see decay_time).
- nominally, pcrit = 8
8
lhwm93
HWM93 neutral wind model.
- lhwm93 = true (use HWM93)
- lhwm93 = false (do not use HWM93)
- nominally use lhwm93 = false
false
lhwm14
HWM14 neutral wind model.
- lhwm14 = true (use HWM14)
- lhwm14 = false (do not use HWM14)
- nominally use lhwm14 = true
true
anu_drag0
anomalous ion drag.
- nominally, use anu_drag0 = 4.e3
- used to prevent high velocity nighttime flows into the ionosphere from the plasmasphere
4.e3
kp
Kp index.
- used in the Hardy model 
- valid values are 1 to 6
xhardy
modifies ionization associated with high latitude precipitation.
- nominal value xhardy = 1
1
tphi
how often the potential equation is solved (sec).
- nominal value tphi = 60
second 60
euv_fac
modifies EUVAC spectrum.
- nominal value euv_fac = 1
- euv_fac < 1 decreases EUV intensity
- euv_fac > 1 increases EUV intensity
1
tmax
maximum electron temperature (K).
- nominal range 10^4 - 10^5
K 1e4
lfism
FISM EUV spectrum.
- lfism = true (use FISM daily spectrum); needs fism linked to FISM year data (up to 2018)
- lfism = false (use EUVAC spectrum)
false
neutral_atmosphere_model
Neutral atmosphere model used in the SAMI3 simulation (NRLMSIS2.1/HWM14 or TIEGCM/WEIMER or ICON-TIEGCM)
Note the allowed interval of ICON-TIEGCM is based on the availability of ICON level 4 data (2019/12/22-2021/12/30) NRLMSIS2.1/HWM14

Run Access Information :

Keyword(s) :

validation

Linked to Parent Spase Resource(s) (example: SimulationModel) :

spase://CCMC/SimulationModel/SAMI3/3.22

List of Publications Referencing this Run :

CMR Footer

Curator: Chiu Wiegand | NASA Official: Dr. Masha Kuznetsova | Privacy and Security Notices | Accessibility | CCMC Data Collection Consent Agreement