# Example 2.5 - A model for Mauna Kea based on median profiles computed from
# NASA MERRA-2 reanalysis data for the period 2006-2015.
#
# Modified to output Trj and Jacobian with respect to trop. h2o.
#
#
# MERRA references:
#   M. M. Rienecker et al. (2011), MERRA: NASA's Modern-Era Retrospective
#   Analysis for Research and Applications. J. Climate 24:3624.
#
#   A. Molod et al (2015), Development of the GEOS-5 atmospheric general
#   circulation model: evolution from MERRA to MERRA2.  Geoscience Model
#   Development 8:1339. 
#
# MERRA is produced by the NASA/Goddard Global Modeling and Assimilation
# Office (GMAO).  MERRA-2 data are archived and distributed by the Goddard
# Earth Sciences Data and Information Services Center (GES DISC).
#
# Data reference:
#   Global Modeling and Assimilation Office (GMAO) (2015),
#   MERRA-2 inst3_3d_asm_Np: 3d,3-Hourly,Instantaneous,Pressure-Level,
#   Assimilation,Assimilated Meteorological Fields,version 5.12.4,
#   Greenbelt, MD, USA:Goddard Space Flight Center Distributed Active
#   Archive Center (GSFC DAAC),
#   Accessed 2016 Feb 10 at doi:10.5067/QBZ6MG944HW0
#
?
? Usage: am example2.5.trop_h2o_Jacobian.amc fmin fmax df zenith_angle Tground trop_h2o_scale
?
? Example: am example2.5.trop_h2o_Jacobian.amc  170 GHz  390 GHz  2 MHz  0 deg  277 K  1.0
?
? Tground is the surface temperature.  At the summit of Mauna Kea, this is
? usually close to the median value of 277 K for the free atmosphere at the
? same altitude.
? 
? trop_h2o_scale is a scaling factor on the median tropospheric water vapor
? profile.  Spectra for other percentile profiles can be approximated by
? scaling the median profile to obtain the same total precipitable water
? vapor (pwv) column as follows:
?
?  H2O profile percentile |    5  |   25  |   50  |   75  |   95
?          trop_h2o_scale | 0.206 | 0.513 | 1.00  | 2.07  | 4.60
?                pwv [mm] | 0.384 | 0.952 | 1.86  | 3.84  | 8.54
?  225 GHz zenith opacity | 0.029 | 0.051 | 0.087 | 0.170 | 0.367
? 
? Setting df larger than about 100 kHz may result in warnings of unresolved
? lines on the lowest-pressure levels in this model.  These warnings may
? be ignored if detailed resolution of line cores is not needed.
?

f %1 %2  %3 %4  %5 %6
output f GHz  Trj K
jacobian Trj
za %7 %8
tol 1e-4

Nscale troposphere h2o %11 1.0

T0 2.7 K

layer mesosphere
Pbase 0.1 mbar
Tbase 222.3 K
lineshape Voigt-Kielkopf
column dry_air vmr
column h2o vmr 6.46e-06
column o3 vmr 1.77e-06

layer mesosphere
Pbase 0.3 mbar
Tbase 243.3 K
lineshape Voigt-Kielkopf
column dry_air vmr
column h2o vmr 6.60e-06
column o3 vmr 1.65e-06

layer mesosphere
Pbase 0.4 mbar
Tbase 249.2 K
lineshape Voigt-Kielkopf
column dry_air vmr
column h2o vmr 6.75e-06
column o3 vmr 1.61e-06

layer mesosphere
Pbase 0.5 mbar
Tbase 254.6 K
lineshape Voigt-Kielkopf
column dry_air vmr
column h2o vmr 6.75e-06
column o3 vmr 1.80e-06

layer mesosphere
Pbase 0.7 mbar
Tbase 262.1 K
lineshape Voigt-Kielkopf
column dry_air vmr
column h2o vmr 6.72e-06
column o3 vmr 2.12e-06

layer stratosphere
Pbase 1 mbar
Tbase 264.7 K
lineshape Voigt-Kielkopf
column dry_air vmr
column h2o vmr 6.63e-06
column o3 vmr 2.66e-06

layer stratosphere
Pbase 2 mbar
Tbase 255.9 K
column dry_air vmr
column h2o vmr 6.34e-06
column o3 vmr 4.01e-06

layer stratosphere
Pbase 3 mbar
Tbase 248.4 K
column dry_air vmr
column h2o vmr 5.97e-06
column o3 vmr 5.88e-06

layer stratosphere
Pbase 4 mbar
Tbase 243.9 K
column dry_air vmr
column h2o vmr 5.69e-06
column o3 vmr 7.28e-06

layer stratosphere
Pbase 5 mbar
Tbase 240.3 K
column dry_air vmr
column h2o vmr 5.47e-06
column o3 vmr 8.13e-06

layer stratosphere
Pbase 7 mbar
Tbase 234.9 K
column dry_air vmr
column h2o vmr 5.25e-06
column o3 vmr 8.75e-06

layer stratosphere
Pbase 10 mbar
Tbase 230.0 K
column dry_air vmr
column h2o vmr 5.01e-06
column o3 vmr 9.10e-06

layer stratosphere
Pbase 20 mbar
Tbase 222.3 K
column dry_air vmr
column h2o vmr 4.72e-06
column o3 vmr 8.01e-06

layer stratosphere
Pbase 30 mbar
Tbase 216.6 K
column dry_air vmr
column h2o vmr 4.49e-06
column o3 vmr 5.67e-06

layer stratosphere
Pbase 40 mbar
Tbase 212.1 K
column dry_air vmr
column h2o vmr 4.37e-06
column o3 vmr 3.60e-06

layer stratosphere
Pbase 50 mbar
Tbase 208.5 K
column dry_air vmr
column h2o vmr 4.30e-06
column o3 vmr 2.23e-06

layer stratosphere
Pbase 70 mbar
Tbase 202.1 K
column dry_air vmr
column h2o vmr 4.30e-06
column o3 vmr 1.23e-06

layer troposphere
Pbase 100 mbar
Tbase 198.5 K
column dry_air vmr
column h2o vmr 4.83e-06
column o3 vmr 4.72e-07

layer troposphere
Pbase 150 mbar
Tbase 207.2 K
column dry_air vmr
column h2o vmr 1.04e-05
column o3 vmr 1.39e-07

layer troposphere
Pbase 200 mbar
Tbase 219.0 K
column dry_air vmr
column h2o vmr 3.40e-05
column o3 vmr 6.74e-08

layer troposphere
Pbase 250 mbar
Tbase 229.9 K
column dry_air vmr
column h2o vmr 9.15e-05
column o3 vmr 5.25e-08

layer troposphere
Pbase 300 mbar
Tbase 239.5 K
column dry_air vmr
column h2o vmr 1.79e-04
column o3 vmr 4.85e-08

layer troposphere
Pbase 350 mbar
Tbase 247.7 K
column dry_air vmr
column h2o vmr 2.76e-04
column o3 vmr 4.69e-08

layer troposphere
Pbase 400 mbar
Tbase 254.7 K
column dry_air vmr
column h2o vmr 3.99e-04
column o3 vmr 4.60e-08

layer troposphere
Pbase 450 mbar
Tbase 260.9 K
column dry_air vmr
column h2o vmr 5.71e-04
column o3 vmr 4.52e-08

layer troposphere
Pbase 500 mbar
Tbase 266.2 K
column dry_air vmr
column h2o vmr 7.94e-04
column o3 vmr 4.45e-08

layer troposphere
Pbase 550 mbar
Tbase 270.9 K
column dry_air vmr
column h2o vmr 1.09e-03
column o3 vmr 4.38e-08

layer troposphere
Pbase 600 mbar
Tbase 275.0 K
column dry_air vmr
column h2o vmr 1.46e-03
column o3 vmr 4.30e-08

layer troposphere
Pbase 625 mbar
Tbase %9  %10
column dry_air vmr
column h2o vmr 1.87e-03
column o3 vmr 4.23e-08