CCPP SciDoc  v6.0.0
Common Community Physics Package Developed at DTC
GFS RRTMG Shortwave/Longwave Radiation Scheme


Radiative processes are among the most complex and computationally intensive parts of all model physics. As an essential component of modeling the atmosphere, radiation directly and indirectly connects all physics processes with model dynamics, and it regulates the overall earth-atmosphere energy exchanges and transformations.

The schematic radiation module structure is shown in Figure 1.

Figure 1: Schematic Radiation Module Structure

GFS radiation package is intended to provide a fast and accurate method of determining the total radiative flux at any given location. These calculations provide both the total radiative flux at the ground surface, which is needed to establish the surface energy budget, and the vertical radiative flux divergence, which is used to calculate the radiative heating and cooling rates of a given atmospheric layer. The magnitude of the terms in the surface energy budget can set the stage for moist deep convection and are crucial to the formation of low-level clouds. In addition, the vertical radiative flux divergence can produce substantial cooling, particularly at the tops of clouds, which can have strong dynamical effects on cloud evolution.

Since 2007, the longwave (LW) and the shortwave (SW) radiation parameterizations in NCEP's operational GFS are both modified and optimized versions of the Rapid Radiative Transfer Model for GCMs (RRTMG_LW v2.3 and RRTMG_SW v2.3, respectively) developed at AER inc. (Atmospheric and Environmental Research) (Iacono et al.(2008) [97]; Mlawer et al.(1997) [136]; Iacono et al.(2000) [96]; Clough et al. 2005 [41]). The LW algorithm contains 140 unevenly distributed g-points (quadrature points) in 16 broad spectral bands, while the SW algorithm includes 112 g-points in 14 bands. In addition to the major atmospheric absorbing gases of ozone, water vapor, and carbon dioxide, the algorithm also includes various minor absorbing species such as methane, nitrous oxide, oxygen, and in the longwave up to four types of halocarbons (CFCs). To represent statistically the unresolved subgrid cloud variability when dealing multi layered clouds, a Monte-Carlo Independent Column Approximation (McICA) method is used in the RRTMG radiative transfer. A maximum-random cloud overlap method is used in both LW and SW radiation calculations.

Radiative fields from model outputs ( \(W m^{-2}\)) include:

  • At surface total sky
    • DLWRFsfc: Downward LW
    • DSWRFsfc: Downward SW
    • ULWRFsfc: Upward LW
    • USWRFsfc: Upward SW
    • NBDSFsfc: Near IR beam downward
    • NDDSFsfc: Near IR diffuse downward
    • VBDSFsfc: UV+Visible beam downward
    • VDDSFsfc: UV+Visible diffuse downward
    • DUVBsfc: UV-B downward flux
  • At surface clear sky
    • CSDLFsfc: Downward LW
    • CSDSFsfc: Downward SW
    • CSULFsfc: Upward LW
    • CSULFsfc: Upward LW
    • CSUSFsfc: Upward sw
    • CDUVBsfc: UV-B downward flux
  • At TOA total sky
    • DSWRFtoa: Downward SW
    • ULWRFtoa: Upward LW
    • USWRFtoa: Upward SW
  • At TOA clear sky:
    • CSULFtoa: Upward LW
    • CSUSFtoa: Upward SW

CCPP Physics Updates

CCPP v6.0.0

Requests have been made by many physics developers and users to rewrite the cloud routines (routines progcld) for radiation computation in the program radiation_clouds.f. Those cloud subroutines are very similar, and have many lines of common code. We modified the radiation_clouds.f module, and includes all the calculations of the cloud properties to a new subroutine radiation_clouds_prop. We also moved the common code from subroutines progcld to this new subroutine. Subroutine radiation_clouds_prop can be called by RRTMG and RRTMGP. A single call to the subroutine radiation_clouds_prop can connect to the calculations of the cloud radiation properties for all the microphysics schemes.

Summary of the major changes:

  • radiation_clouds.f

    A new subroutine “radiation_clouds_prop” was added to radiation_clouds.f. This new subroutine calculates all cloud radiation properties for all the microphysics schemes. Subroutines "progcld*" were renamed based on the input variables imp_physics, and the inactive subroutines were removed from file radiation_clouds.f

    • 'progcld1' — > progcld_zhao_carr
    • 'progcld3' — > progcld_zhao_carr_pdf
    • 'progcld4' — > progcld_gfdl_lin
    • 'progcld5' — > progcld_fer_hires
    • 'progcld6' — > progcld_thompson_wsm6
    • 'progclduni' — > progclduni
    • 'progcld_thompson'— > progcld_thompson
  • GFS_rrtmg_pre.F90

    Removed the “progcld” subroutine calls, and replaced them with a single subroutine call to “radiation_clouds_prop”.

  • radiation_cloud_overlap.F90

    Replaced subroutine “get_alpha_exp” with “get_alpha_exper”. The new subroutine revises alpha for exponential random cloud overlap option. This new subroutine is used in programs GFS_rrtmgp_cloud_overlap_pre.F90 and GFS_rrtmgp_gfdlmp_pre.F90.

  • Subroutine getml() has been modified. The subroutine computes low, mid, high, total and boundary clouds, and is used in GFS_cloud_diagnostics.F90.

Intraphysics Communication

General Algorithm