GFS_v16 Suite

Overview

Version 16 of the Global Forecast System (GFS) was implemented operationally by the NOAA National Centers for Environmental Prediction (NCEP) in 2021. While GFS version 16 did not make use of CCPP physics suites, the GFS_v16 suite is meant to emulate the physics used in that operational implementation. The GFS_v16 suite uses the parameterizations in the following order:

• RRTMG Shortwave/Longwave Radiation Scheme
• GFS Surface Layer Scheme
• Near-Surface Sea Temperature Scheme
• GFS Simple Ocean Scheme
• Noah Land Surface Model
• Sea Ice Scheme
• GFS Scale-aware TKE-based Moist Eddy-Diffusion Mass-Flux (EDMF) PBL and Free Atmospheric Turbulence Scheme
• Unified Gravity Wave Physics (GWP) Scheme - Version 0
• Ozone Photochemistry (2015) Scheme
• Stratospheric H2O Scheme
• GFS Scale-Aware Simplified Arakawa-Schubert (sa-SAS) Deep Convection Scheme
• GFS SAS-based Mass-Flux Scheme for Shallow convection (sa-MF)
• GFDL Cloud Microphysics Scheme

Suite Definition File

<?xml version="1.0" encoding="UTF-8"?>
 
<suite name="SCM_GFS_v16" version="1">
  <group name="time_vary">
    <subcycle loop="1">
      <scheme>GFS_time_vary_pre</scheme>
      <scheme>GFS_rrtmg_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_rrtmg_pre</scheme>
      <scheme>GFS_radiation_surface</scheme>
      <scheme>rad_sw_pre</scheme>
      <scheme>rrtmg_sw</scheme>
      <scheme>rrtmg_sw_post</scheme>
      <scheme>rrtmg_lw</scheme>
      <scheme>rrtmg_lw_post</scheme>
      <scheme>GFS_rrtmg_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

GFDL Cloud Microphysics Scheme namelist options

&gfdl_cloud_microphysics_nml
    c_cracw       = 0.8
    c_paut        = 0.5
    c_pgacs       = 0.01
    c_psaci       = 0.05
    ccn_l         = 300.0
    ccn_o         = 100.0
    const_vg      = .false.
    const_vi      = .false.
    const_vr      = .false.
    const_vs      = .false.
    de_ice        = .false.
    do_qa         = .false.
    do_sedi_heat  = .false.
    dw_land       = 0.16
    dw_ocean      = 0.1
    fast_sat_adj  = .false.
    fix_negative  = .true.
    icloud_f      = 1
    mono_prof     = .true.
    mp_time       = 150.0
    prog_ccn      = .false.
    qi0_crt       = 8e-05
    qi_lim        = 1.0
    ql_gen        = 0.001
    ql_mlt        = 0.001
    qs0_crt       = 0.001
    rad_graupel   = .true.
    rad_rain      = .true.
    rad_snow      = .true.
    reiflag       = 2
    rh_inc        = 0.3
    rh_inr        = 0.3
    rh_ins        = 0.3
    rthresh       = 1e-05
    sedi_transport = .true.
    tau_g2v       = 900.0
    tau_i2s       = 1000.0
    tau_l2v       = 225.0
    tau_v2l       = 150.0
    use_ccn       = .true.
    use_ppm       = .false.
    vg_max        = 12.0
    vi_max        = 1.0
    vr_max        = 12.0
    vs_max        = 2.0
    z_slope_ice   = .true.
    z_slope_liq   = .true.
/

Other namelist 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.
/

• 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"