GFS_v16_RRTMGP Suite

Overview

This is an experimental version of the Global Forecast System (GFS) in which the RRTMG radiation parameterization used operationally is replaced with the RTE-RRTMGP parameterization for research and development.

Suite Definition File

<?xml version="1.0" encoding="UTF-8"?>
 
<suite name="SCM_GFS_v16_RRTMGP" version="1">
  <group name="time_vary">
    <subcycle loop="1">
      <scheme>GFS_time_vary_pre</scheme>
      <scheme>GFS_rrtmgp_setup</scheme>
      <scheme>GFS_rad_time_vary</scheme>
      <scheme>GFS_phys_time_vary</scheme>
    </subcycle>
  </group>
  <group name="radiation">
    <subcycle loop="1">
      <scheme>GFS_suite_interstitial_rad_reset</scheme>
      <scheme>GFS_rrtmgp_pre</scheme>
      <scheme>GFS_radiation_surface</scheme>
      <scheme>GFS_rrtmgp_cloud_mp</scheme>
      <scheme>GFS_rrtmgp_cloud_overlap</scheme>
      <scheme>GFS_cloud_diagnostics</scheme>
      <scheme>rrtmgp_aerosol_optics</scheme>
      <scheme>rrtmgp_sw_main</scheme>
      <scheme>rrtmgp_lw_main</scheme>
      <scheme>GFS_rrtmgp_post</scheme>
    </subcycle>
  </group>
  <group name="physics">
    <subcycle loop="1">
      <scheme>GFS_suite_interstitial_phys_reset</scheme>
      <scheme>GFS_suite_stateout_reset</scheme>
      <scheme>get_prs_fv3</scheme>
      <scheme>GFS_suite_interstitial_1</scheme>
      <scheme>GFS_surface_generic_pre</scheme>
      <scheme>GFS_surface_composites_pre</scheme>
      <scheme>dcyc2t3</scheme>
      <scheme>GFS_surface_composites_inter</scheme>
      <scheme>GFS_suite_interstitial_2</scheme>
    </subcycle>
    <!-- Surface iteration loop -->
    <subcycle loop="2">
      <scheme>sfc_diff</scheme>
      <scheme>GFS_surface_loop_control_part1</scheme>
      <scheme>sfc_nst_pre</scheme>
      <scheme>sfc_nst</scheme>
      <scheme>sfc_nst_post</scheme>
      <scheme>lsm_noah</scheme>
      <scheme>sfc_sice</scheme>
      <scheme>GFS_surface_loop_control_part2</scheme>
    </subcycle>
    <!-- End of surface iteration loop -->
    <subcycle loop="1">
      <scheme>GFS_surface_composites_post</scheme>
      <scheme>sfc_diag</scheme>
      <scheme>sfc_diag_post</scheme>
      <scheme>GFS_surface_generic_post</scheme>
      <scheme>GFS_PBL_generic_pre</scheme>
      <scheme>satmedmfvdifq</scheme>
      <scheme>GFS_PBL_generic_post</scheme>
      <scheme>GFS_GWD_generic_pre</scheme>
      <scheme>cires_ugwp</scheme>
      <scheme>cires_ugwp_post</scheme>
      <scheme>GFS_GWD_generic_post</scheme>
      <scheme>GFS_suite_stateout_update</scheme>
      <scheme>h2ophys</scheme>
      <scheme>get_phi_fv3</scheme>
      <scheme>GFS_suite_interstitial_3</scheme>
      <scheme>GFS_DCNV_generic_pre</scheme>
      <scheme>samfdeepcnv</scheme>
      <scheme>GFS_DCNV_generic_post</scheme>
      <scheme>GFS_SCNV_generic_pre</scheme>
      <scheme>samfshalcnv</scheme>
      <scheme>GFS_SCNV_generic_post</scheme>
      <scheme>GFS_suite_interstitial_4</scheme>
      <scheme>cnvc90</scheme>
      <scheme>GFS_MP_generic_pre</scheme>
      <scheme>gfdl_cloud_microphys</scheme>
      <scheme>GFS_MP_generic_post</scheme>
      <scheme>maximum_hourly_diagnostics</scheme>
      <scheme>GFS_physics_post</scheme>
    </subcycle>
  </group>
</suite>

Namelist

• General physics options

  &gfs_physics_nml
      cal_pre       = .false.
      cdmbgwd       = 4.0, 0.15, 1.0, 1.0
      cnvcld        = .true.
      cnvgwd        = .true.
      debug         = .false.
      do_shum       = .false.
      do_skeb       = .false.
      do_spp        = .false.
      do_sppt       = .false.
      do_tofd       = .true.
      do_ugwp       = .false.
      dspheat       = .true.
      effr_in       = .true.
      fhcyc         = 0
      fhlwr         = 1200.0
      fhswr         = 1200.0
      fhzero        = 1.0
      h2o_phys      = .true.
      hybedmf       = .false.
      iaer          = 5111
      ialb          = 1
      iau_inc_files = ''
      icliq_sw      = 2
      ico2          = 2
      iems          = 1
      imfdeepcnv    = 2
      imfshalcnv    = 2
      imp_physics   = 11
      iopt_alb      = 2
      iopt_btr      = 1
      iopt_crs      = 1
      iopt_dveg     = 1
      iopt_frz      = 1
      iopt_inf      = 1
      iopt_rad      = 1
      iopt_run      = 1
      iopt_sfc      = 1
      iopt_snf      = 4
      iopt_stc      = 1
      iopt_tbot     = 2
      iovr          = 3
      isatmedmf     = 1
      isol          = 2
      isot          = 1
      isubc_lw      = 2
      isubc_sw      = 2
      ivegsrc       = 1
      ldiag3d       = .true.
      ldiag_ugwp    = .false.
      lgfdlmprad    = .true.
      lheatstrg     = .true.
      lndp_type     = 0
      lsm           = 1
      lsoil         = 4
      lwhtr         = .true.
      n_var_lndp    = 0
      n_var_spp     = 0
      nsfullradar_diag = 3600
      nst_anl       = .true.
      nstf_name     = 2, 1, 0, 0, 0
      oz_phys       = .false.
      oz_phys_2015  = .true.
      pdfcld        = .false.
      prautco       = 0.00015, 0.00015
      pre_rad       = .false.
      print_diff_pgr = .false.
      prslrd0       = 0.0
      psautco       = 0.0008, 0.0005
      qdiag3d       = .true.
      random_clds   = .false.
      redrag        = .true.
      satmedmf      = .true.
      shal_cnv      = .true.
      swhtr         = .true.
      trans_trac    = .true.
      use_ufo       = .true.
      do_RRTMGP          = .true.
      doGP_cldoptics_LUT = .true.
      doGP_lwscat        = .true.
      active_gases       = 'h2o_co2_o3_n2o_ch4_o2'
      ngases             = 6
      rrtmgp_root        = '../../ccpp/physics/physics/rte-rrtmgp/'
      lw_file_gas        = 'rrtmgp/data/rrtmgp-data-lw-g128-210809.nc'
      lw_file_clouds     = 'extensions/cloud_optics/rrtmgp-cloud-optics-coeffs-lw.nc'
      sw_file_gas        = 'rrtmgp/data/rrtmgp-data-sw-g112-210809.nc'
      sw_file_clouds     = 'extensions/cloud_optics/rrtmgp-cloud-optics-coeffs-sw.nc'
      rrtmgp_nGptsSW     = 112
      rrtmgp_nGptsLW     = 128
      rrtmgp_nBandsLW    = 16
      rrtmgp_nBandsSW    = 14
  /

• GFS_UGWP_v0_UGWP related options

  &cires_ugwp_nml
      knob_ugwp_azdir    = 2, 4, 4, 4
      knob_ugwp_doaxyz   = 1
      knob_ugwp_doheat   = 1
      knob_ugwp_dokdis   = 1
      knob_ugwp_effac    = 1, 1, 1, 1
      knob_ugwp_ndx4lh   = 1
      knob_ugwp_solver   = 2
      knob_ugwp_source   = 1, 1, 0, 0
      knob_ugwp_stoch    = 0, 0, 0, 0
      knob_ugwp_version  = 0
      knob_ugwp_wvspec   = 1, 25, 25, 25
      launch_level       = 27
  /

• nstf_name = \([2,0,0,0,0]^1 [2,1,0,0,0]^2\)
  • \(^1\) NSST is on and coupled with spin up off
  • \(^2\) NSST is on and coupled with spin up on
• Grid-spacing dependency of cdmbgwd
  • C48: cdmbgwd="0.071,2.1,1.0,1.0"
  • C96: cdmbgwd="0.14,1.8,1.0,1.0"
  • C192: cdmbgwd="0.23,1.5,1.0,1.0"
  • C384: cdmbgwd="1.1,0.72,1.0,1.0"
  • C768: cdmbgwd="4.0,0.15,1.0,1.0"
  • C1152: cdmbgwd="4.0,0.10,1.0,1.0"
  • C3072: cdmbgwd="4.0,0.05,1.0,1.0"