This is the repository for RTE+RRTMGP, a set of codes for computing radiative fluxes in planetary atmospheres. RTE+RRTMGP is described in a paper in Journal of Advances in Modeling Earth Systems.
RRTMGP uses a k-distribution to provide an optical description (absorption and possibly Rayleigh optical depth) of the gaseous atmosphere, along with the relevant source functions, on a pre-determined spectral grid given temperatures, pressures, and gas concentration. The k-distribution currently distributed with this package is applicable to the Earth's atmosphere under present-day, pre-industrial, and 4xCO2 conditions.
RTE computes fluxes given spectrally-resolved optical descriptions and source functions. The fluxes are normally summarized or reduced via a user extensible class.
Example programs and documentation are evolving - please see examples/ in the repo and Wiki on the project's Github page. Suggestions are welcome. Meanwhile for questions please contact Robert Pincus and Eli Mlawer at rrtmgp@aer.com.
- The default method for solution for longwave problems that include scattering has been changed from 2-stream methods to a re-scaled and refined no-scattering calculation following Tang et al. 2018.
- In RRTMGP gas optics, the spectrally-resolved solar source function in can be adjusted by specifying the total solar irradiance (
gas_optics%set_tsi(tsi)
) and/or the facular and sunspot indicies (gas_optics%set_solar_variability(mg_index, sb_index, tsi)
)from the NRLSSI2 model of solar variability. rte_lw()
now includes optional arguments for computing the Jacobian (derivative) of broadband flux with respect to changes in surface temperature. In calculations neglecting scattering only the Jacobian of upwelling flux is computed. When using re-scaling to account for scattering the Jacobians of both up- and downwelling flux are computed.- A new module,
mo_rte_config
, contains two logical variables that indicate whether arguments to routines are to be checked for correct extents and/or valid values. These variables can be changed via calls torte_config_checks()
. Setting the values to.false.
removes the checks. Invalid values may cause incorrect results, crashes, or other mayhem
Relative to commit 69d36c9
to master
on Apr 20, 2020, the required arguments to both the longwave and shortwave versions of ty_gas_optics_rrtmgp%load()
have changed.
- Set environment variables
FC
(the Fortran 2003 compiler) andFCFLAGS
(compiler flags). Examples are provided in theCompiler-flags.md
file. - Set environment variables
RRTMGP_ROOT
to the top-level RTE+RRTMGP directory andRTE_KERNELS
toopenacc
if you want the OpenACC/OpenMP kernels rather than the default. make libs
in the top-level directory will make the RTE and RRTMGP libraries.- The examples and testing codes use netCDF. Set the variables
NCHOME
andNFHOME
to the roots of the C and Fortran netCDF installations, thenmake tests
to build and run these. (A few files need to be downloaded forexamples/rfmaip-clear-sky
. The default is to download these withwget
but a Python script is also available.) - Evaluating the results of the tests requires
Python
with thexarray
package and its depdencies installed. Comparisons can be made withmake check
in the top level directory. make
invoked without a target in the top level attempts all three steps.
Two examples are provided in examples/
, one for clear skies and one including clouds. Directory tests/
contains regression testing (e.g. to ensure that answers are independent of orientation) and unit testing (to be sure all the code paths are tested). See the README file and codes in each directory for further information.