============================================================================================================================== Included are 8 python readable IDL save files. (import scipy as sp, import scipy.io as spio, Python command: spio.readsav('$PATH/figureX.sav',verbose=True), where $PATH is the IDL save file local path and X is the figure number 1-8.) (IDL command: restore,'$PATH/figureX.sav',description=description, where $PATH is the IDL save file local path and Xis the figure number 1-8.) The variables for each are file are described below. ============================================================================================================================== 'figure1.sav' Description = Global average mass density and their IAVs relative to their global annual averages, and the change in the IAVs from the modified versus the standard simulation. Model Simulation(s) = TIME-GCM Standard, Reduced O+O+M, and Double H+O3 doy = day of year (days) ndoy = number of days masden_stnd = Total neutral mass density at 400 km from the Standard simulation (x10e-15 g cm-3) masden_red = Total neutral mass density at 400 km from the Reduced O+O+M simulation (x10e-15 g cm-3) masden_dbl = Total neutral mass density at 400 km from the Double H+O3 simulation (x10e-15 g cm-3) per_masden_stnd = IAV in total neutral mass density at 400 km from the Standard simulation (%) per_masden_red = IAV in total neutral mass density at 400 km from the Reduced O+O+M simulation (%) per_masden_dbl = IAV in total neutral mass density at 400 km from the Double H+O3 simulation (%) deliav_red = Difference in total neutral mass density IAV, Reduced O+O+M - Standard simulation (%) deliav_dbl = Difference in total neutral mass density IAV, Double H+O3 - Standard simulation (%) ============================================================================================================================== 'figure2.sav' Description = Global annual average vertical profile of O/Ox (O+O3) ratio from the Standard simulation. Model Simulation(s) = TIME-GCM Standard logp_mid = model log-pressure level (logp_mid to hpa = p0*exp(-logp_mid(k))*1.e-3, p0=5.e-4 microbars) np_mid = number of model log-pressure levels oox_ratio = atomic oxygen to total odd oxygen (O+O3) ratio ============================================================================================================================== 'figure3.sav' Description = Terms of the Ox and O2 continuity equation (4) as a function of day of year. Model Simulation = TIME-GCM Standard doy = day of year (days) ndoy = number of days ox_chem_r2 = total odd oxygen chemical term in Region 2 (% day-1) ox_dpsidt_r2 = time rate of change of odd oxygen in Region 2 (% day-1) ox_tvertrans_r2ub = total vertical transport of odd oxygen at Region 2 upper boundary (% day-1) ox_tvertrans_r2lb = total vertical transport of odd oxygen at Region 2 lower boundary (% day-1) ox_chem_r1 = total odd oxygen chemical term in Region 1 (% day-1) ox_dpsidt_r1 = time rate of change of odd oxygen in Region 1 (% day-1) ox_tvertrans_r1ub = total vertical transport of odd oxygen at Region 1 upper boundary (% day-1) ox_tvertrans_r1lb = total vertical transport of odd oxygen at Region 1 lower boundary (% day-1) o2_chem_r2 = total molecular oxygen chemical term in Region 2 (% day-1) o2_dpsidt_r2 = time rate of change of molecular oxygen in Region 2 (% day-1) o2_tvertrans_r2ub = total vertical transport of molecular oxygen at Region 2 upper boundary (% day-1) o2_tvertrans_r2lb = total vertical transport of molecular oxygen at Region 2 lower boundary (% day-1) o2_chem_r1 = total molecular oxygen chemical term in Region 1 (% day-1) o2_dpsidt_r1 = time rate of change of molecular oxygen in Region 1 (% day-1) o2_tvertrans_r1ub = total vertical transport of molecular oxygen at Region 1 upper boundary (% day-1) o2_tvertrans_r1lb = total vertical transport of molecular oxygen at Region 1 lower boundary (% day-1) ============================================================================================================================== 'figure4.sav' Description = Chemical production and loss terms of Ox and O2 from equation (4) in Region 2. Model Simulation(s) = TIME-GCM Standard doy = day of year (days) ndoy = number of days o2_chem_r2 = total molecular oxygen chemical term in Region 2 (% day-1) ox_chem_r2 = total odd oxygen chemical term in Region 2 (% day-1) pdo2_losso2_r2 = photochemical loss of molecular oxygen from photodissociation in Region 2 (% day-1) pdo2_prodox_r2 = photochemical production of odd oxygen (o) from photodissociation of molecular oxygen in Region 2 (% day-1) ============================================================================================================================== 'figure5.sav' Description = Chemical production and loss terms of Ox from equation (4) in Region 1 as a function of day of year. Model Simulation(s) = TIME-GCM Standard doy = day of year (days) ndoy = number of days p_lox_r1 = net production-loss of odd oxygen in Region 1 (% day-1) prod_ox_r1 = odd oxygen production in Region 1 (% day-1) loss_ox_r1 = odd oxygen loss in Region 1 (% day-1) rkm20_oom_r1 = O + O + M loss rate (% day-1) rkm25_ooh_r1 = O + OH loss rate (% day-1) rkm26_oho2_r1 = O + HO2 loss rate (% day-1) rkm37_ho3_r1 = H + O3 loss rate (% day-1) ============================================================================================================================== 'figure6.sav' Description = Vertical profiles of Ox chemical and vertical transport terms from equation 5, including the relative changes between the Double H+O3 and Standard simulations at March equinox and June solstice. Model Simulation(s) = TIME-GCM Standard, Reduced O+O+M, and Double H+O3 logp_mid = model log-pressure level (logp_mid to hpa = p0*exp(-logp_mid(k))*1.e-3, p0=5.e-4 microbars) np_mid = number of model log-pressure levels ox_chemdz_stnd_eqnx = total odd oxygen chemical term from the Standard simulation on day of year 80 (x10e2 % day-1) ox_tvertransdz_stnd_eqnx = total odd oxygen vertical transport term from the Standard simulation on day of year 80 (x10e2 % day-1) ox_chemdz_stnd_sol = total odd oxygen chemical term from the Standard simulation on day of year 172 (x10e2 % day-1) ox_tvertransdz_stnd_sol = total odd oxygen vertical transport term from the Standard simulation on day of year 172 (x10e2 % day-1) ox_chemdz_red_eqnx = total odd oxygen chemical term from the Reduced O+O+M simulation on day of year 80 (x10e2 % day-1) ox_tvertransdz_red_eqnx = total odd oxygen vertical transport term from the Reduced O+O+M simulation on day of year 80 (x10e2 % day-1) ox_chemdz_red_sol = total odd oxygen chemical term from the Reduced O+O+M simulation on day of year 172 (x10e2 % day-1) ox_tvertransdz_red_sol = total odd oxygen vertical transport term from the Reduced O+O+M simulation on day of year 172 (x10e2 % day-1) ox_chemdz_dbl_eqnx = total odd oxygen chemical term from the Double H+O3 simulation simulation on day of year 80 (x10e2 % day-1) ox_tvertransdz_dbl_eqnx = total odd oxygen vertical transport term from the Double H+O3 simulation on day of year 80 (x10e2 % day-1) ox_chemdz_dbl_sol = total odd oxygen chemical term from the Double H+O3 simulation on day of year 172 (x10e2 % day-1) ox_tvertransdz_dbl_sol = total odd oxygen vertical transport term from the Double H+O3 simulation on day of year 172 (x10e2 % day-1) ox_chemdz_perchg_dbl_eqnx = relative change in the total odd oxygen chemical term on day of year 80 ((Double H+O3 - Standard)/MOx_Stnd, % day-1) ox_wbulkdz_perchg_dbl_eqnx = relative change in the odd oxygen vertical transport by the bulk wind on day of year 80 (% day-1) ox_weddydz_perchg_dbl_eqnx = relative change in the odd oxygen vertical transport by eddy diffusion on day of year 80 (% day-1) ox_chemdz_perchg_dbl_sol = relative change in the total odd oxygen chemical term on day of year 172 (% day-1) ox_wbulkdz_perchg_dbl_sol = relative change in the odd oxygen vertical transport by the bulk wind day of year 172 (% day-1) ox_weddydz_perchg_dbl_sol = relative change in the odd oxygen vertical transport by eddy diffusion on day of year 172 (% day-1) ============================================================================================================================== 'figure7.sav' Description = 'Geometric height of the [O] peak in the TIME-GCM standard simulation (in km as a function of latitude lat and doy) and the o mass mixing ratios at log-pressure level -7 from the different TIME-GCM simulations (ommr_plev7_***, units mmr x 10-3) and their sao amplitudes (sao_ommr_***, units x 10-3) as a function of day of year (doy). doy = day of year (days) ndoy = number of days lat = latitude (deg) ny = number of latitudes zg_omax = geometric height of the [O] peak (km) ommr_plev7_stnd = Zonally and diurnally averaged O mass mixing ratio at log-pressure level -7 from the Standard simulation (mmr x 10-3) ommr_plev7_dbl = Zonally and diurnally averaged O mass mixing ratio at log-pressure level -7 from the Double H+O3 simulation (mmr x 10-3) sao_ommr_stnd = Global average O mass mixing SAO from the Standard simulation (mmr x 10-3) sao_ommr_dbl = Global average O mass mixing SAO from the Double H+O3 simulation (mmr x 10-3) ============================================================================================================================== 'figure8.sav' Description = 'IAVs relative to their global annual averages in mass density at 400 km from TIE-GCM with and without MSIS composition at the model lower boundary. AO Residual mass densities at 400 km from TIE-GCM with and without MSIS composition at the model lower boundary. Day-of-year vs. altitude representation of [O] IAVs relative to their global annual averages in 5 km intervals from a TIE-GCM simulation with MSIS composition at the model lower boundary. Model Simulation(s) = TIE-GCM MSIS+HWM U,V,T,Z and TIE-GCM MSIS+HWM U,V,T,Z + Comp. doy = day of year (days) ndoy = number of days hgt = altitude (km) nhgt = number of altitudes per_masden_msishwm = IAV in total neutral mass density at 400 km from the TIE-GCM MSIS+HWM U,V,T,Z simulation (%) per_masden_msishwm_comp = IAV in total neutral mass density at 400 km from the TIE-GCM MSIS+HWM U,V,T,Z + Comp. simulation (%) per_aoresid_msishwm = AO residuals in total neutral mass density at 400 km from the TIE-GCM MSIS+HWM U,V,T,Z simulation (%) per_aoresid_msishwm_comp = AO residuals in total neutral mass density at 400 km from the TIE-GCM MSIS+HWM U,V,T,Z + Comp. simulation (%) per_o1_msishwm_comp = IAV in [O] as a function of altitude from the TIE-GCM MSIS+HWM U,V,T,Z + Comp. simulation (%) ==============================================================================================================================