Work Description

Each of these directories contain pickled files for each model runs.
These are prepared using Python2.7 and can be read as shown in
read_pickle.py:

The runs in each directory are full year simulations of GITM
with different lower boundary (LB) conditions. For each run, diurnally
and zonally averaged parameters are provided in a dictionary data structure.
These can be read using the function load_gitm_pickle()

1. gitm_msis/ : G/MSIS simulation (NRLMSIS00 and HWM14)
2. gitm_noSAO/ : G/NOSAO simulation (AO, SAO Flags = 0)
3. gitm_waccm_lb_cesm2.0/ : G/WX (WACCM-X 2.0 as LB)
4. gitm_waccm_nudge/ : G/NUDGE (WACCMX 2.0 as LB with constrained dynamics)
5. gitm_waccm_eddy/ : G/EDDY (WACCM-X 2.0 as the LB with seasonally
varying Kzz)
6. gitm_waccm_lb_cesm2.1/ : G/WX 2.1 as the LB
7. pure_msis : MSIS empirical model output on the GITM grid

Point to note : These GITM simulations are 3D simulations with
lat x lon x alt as the coordinate system.

Each simulation/directory has quantities that are stored as pickled files
Some quantities are stored as dictionaries where the keys represent the
day of the year, 2010. And the values are daily and zonal averages with
shape, lats x alts. These quantities are as follows :

1. alti_daily.p : Altitude in m
2. N2_daily.p : Nitrogen number density(N2) in m-3
3. T_daily.p : Temperature(T) in K
4. O_daily.p : Atomic Oxygen number density(O) in m-3
5. rho_daily.p : Neutral Mass Density(rho) in kgm-3
6. U_daily.p : Zonal Wind Velocity(U) in m/s
7. V_daily.p : Meridional Wind Velocity(V) in m/s
8. W_daily.p : Total Vertical Wind Velocity(W) in m/s
9. W_O_daily.p : Total Vertical Wind Velocity(W) for O in m/s
10.W_O2_daily.p : Total Vertical Wind Velocity(W) for O2 in m/s
11.W_N2_daily.p : Total Vertical Wind Velocity(W) for N2 in m/s

The following quantities are stored as one-dimensional arrays :
12. lats : Latitude in radians
13. day_dtimes : Dates in datetime format

WACCM-X 2.0 daily and zonally averaged 2010 data is in the folder
waccmx_cesm2.0/ and can be read using the load_waccmx_pickle().
The output of the model is interpolated onto an altitude grid of GITM,
And the latitude grid can be defined as np.linspace(-90,90,96).

1. comb_waccmx_N2.p : Nitrogen number density(N2) in m-3
2. comb_waccmx_O.p : Atomic Oxygen number density(O) in m-3
3. comb_waccmx_O2.p : Molecular Oxygen number density(O) in m-3
4. comb_waccmx_U.p : Zonal Wind Velocity(U) in m/s
5. comb_waccmx_V.p : Meridional Wind Velocity(U) in m/s
6. comb_waccmx_W.p : Total Vertical Wind Velocity(U) in m/s
7. comb_waccmx_T.p : Temperature(T) in K

These quantities are used as the lower boundary conditions for different
GITM simulations.