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. The GFS_v16 suite uses the parameterizations in the following order:

• GFS RRTMG Shortwave/Longwave Radiation Scheme
• GFS Surface Layer Scheme
• GFS Near-Surface Sea Temperature Scheme
• GFS Simple Ocean Scheme
• GFS Noah Land Surface Model
• GFS Sea Ice Scheme
• GFS Scale-aware TKE-based Moist Eddy-Diffusion Mass-Flux (EDMF) PBL and Free Atmospheric Turbulence Scheme
• Unified Gravity Wave Physics Scheme - Version 0
• GFS Ozone Photochemistry (2015) Scheme
• GFS 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="FV3_GFS_v16" version="1">
  <!-- <init></init> -->
  <group name="fast_physics">
    <subcycle loop="1">
      <scheme>fv_sat_adj</scheme>
    </subcycle>
  </group>
  <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_pre</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>ozphys_2015</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>phys_tend</scheme>
    </subcycle>
  </group>
  <group name="stochastics">
    <subcycle loop="1">
      <scheme>GFS_stochastics</scheme>
    </subcycle>
  </group>
  <!-- <finalize></finalize> -->
</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        = .true.
    do_sedi_heat = .false.
    dw_land      = 0.16
    dw_ocean     = 0.1
    fast_sat_adj = .true.
    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_tofd      = .true.
    do_ugwp      = .false.
    dspheat      = .true.
    effr_in      = .true.
    fhcyc        = 0.0
    fhlwr        = 3600.0
    fhswr        = 3600.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
    iopt_trs     = 2
    iovr         = 3
    isatmedmf    = 1
    isol         = 2
    isot         = 1
    isubc_lw     = 2
    isubc_sw     = 2
    ivegsrc      = 1
    ldiag3d      = .false.
    ldiag_ugwp   = .false.
    lgfdlmprad   = .true.
    lheatstrg    = .true.
    lsm          = 1
    lsoil        = 4
    lwhtr        = .true.
    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
    random_clds  = .false.
    redrag       = .true.
    satmedmf     = .true.
    sfclay_compute_flux = .false.
    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"