GFS_v16beta Suite

Overview

Version 16 of the Global Forecast System (GFS) will be implemented operationally by the NOAA National Centers for Environmental Prediction (NCEP) in 2021. GFS_v16beta is a prototype of the GFS_v16 suite. The main difference between the GFS_v15p2 and GFS_v16beta suites is the replacement of the K-based EDMF PBL scheme with a moist TKE based one.

The GFS_v16beta physics 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 Rayleigh Damping Calculation
• 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_v16beta" lib="ccppphys" ver="5">
  <!-- <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>rrtmg_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>rayleigh_damp</scheme>
      <scheme>GFS_suite_stateout_update</scheme>
      <scheme>ozphys_2015</scheme>
      <scheme>h2ophys</scheme>
      <scheme>GFS_DCNV_generic_pre</scheme>
      <scheme>get_phi_fv3</scheme>
      <scheme>GFS_suite_interstitial_3</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>
    </subcycle>
  </group>
  <group name="stochastics">
    <subcycle loop="1">
      <scheme>GFS_stochastics</scheme>
    </subcycle>
  </group>
  <!-- <finalize></finalize> -->
</suite>

Namelist

&gfs_physics_nml
  fhzero            = 6
  h2o_phys          = .true.
  ldiag3d           = .false.
  fhcyc             = 24
  use_ufo           = .true.
  pre_rad           = .false.
  ncld              = 5
  imp_physics       = 11
  pdfcld            = .false.
  fhswr             = 3600.
  fhlwr             = 3600.
  ialb              = 1
  iems              = 1
  iaer              = 5111
  icliq_sw          = 2
  iovr_lw           = 3
  iovr_sw           = 3
  ico2              = 2
  isubc_sw          = 2
  isubc_lw          = 2
  isol              = 2
  lwhtr             = .true.
  swhtr             = .true.
  cnvgwd            = .true.
  shal_cnv          = .true.
  cal_pre           = .false.
  redrag            = .true.
  dspheat           = .true.
  hybedmf           = .false.
  satmedmf          = .true.
  isatmedmf         = 1
  lheatstrg         = .true.
  random_clds       = .false.
  trans_trac        = .true.
  cnvcld            = .true.
  imfshalcnv        = 2
  imfdeepcnv        = 2
  cdmbgwd           = 4.0,0.15,1.0,1.0
  prslrd0           = 0.
  ivegsrc           = 1
  isot              = 1
  lsoil             = 4
  lsm               = 1
  iopt_dveg         = 1
  iopt_crs          = 1
  iopt_btr          = 1
  iopt_run          = 1
  iopt_sfc          = 1
  iopt_frz          = 1
  iopt_inf          = 1
  iopt_rad          = 1
  iopt_alb          = 2
  iopt_snf          = 4
  iopt_tbot         = 2
  iopt_stc          = 1
  debug             = .false.
  oz_phys           = .F.
  oz_phys_2015      = .T.
  nstf_name         = @[NSTF_NAME]
  nst_anl           = .true.
  psautco           = 0.0008,0.0005
  prautco           = 0.00015,0.00015
  lgfdlmprad        = .true.
  effr_in           = .true.
  ldiag_ugwp        = .false.
  do_ugwp           = .false.
  do_tofd           = .true.
  do_sppt           = .false.
  do_shum           = .false.
  do_skeb           = .false.
  do_sfcperts       = .false.
/

&gfdl_cloud_microphysics_nml
  sedi_transport    = .true.
  do_sedi_heat      = .false.
  rad_snow          = .true.
  rad_graupel       = .true.
  rad_rain          = .true.
  const_vi          = .F.
  const_vs          = .F.
  const_vg          = .F.
  const_vr          = .F.
  vi_max            = 1.
  vs_max            = 2.
  vg_max            = 12.
  vr_max            = 12.
  qi_lim            = 1.
  prog_ccn          = .false.
  do_qa             = .true.
  fast_sat_adj      = .true.
  tau_l2v           = 225.
  tau_v2l           = 150.
  tau_g2v           = 900.
  rthresh           = 10.e-6  
  dw_land           = 0.16
  dw_ocean          = 0.10
  ql_gen            = 1.0e-3
  ql_mlt            = 1.0e-3
  qi0_crt           = 8.0E-5
  qs0_crt           = 1.0e-3
  tau_i2s           = 1000.
  c_psaci           = 0.05
  c_pgacs           = 0.01
  rh_inc            = 0.30
  rh_inr            = 0.30
  rh_ins            = 0.30
  ccn_l             = 300.
  ccn_o             = 100.
  c_paut            = 0.5
  c_cracw           = 0.8
  use_ppm           = .false.
  use_ccn           = .true.
  mono_prof         = .true.
  z_slope_liq       = .true.
  z_slope_ice       = .true.
  de_ice            = .false.
  fix_negative      = .true.
  icloud_f          = 1
  mp_time           = 150.
  reiflag           = 2
/

&cires_ugwp_nml
  knob_ugwp_solver  = 2
  knob_ugwp_source  = 1,1,0,0
  knob_ugwp_wvspec  = 1,25,25,25
  knob_ugwp_azdir   = 2,4,4,4
  knob_ugwp_stoch   = 0,0,0,0
  knob_ugwp_effac   = 1,1,1,1
  knob_ugwp_doaxyz  = 1
  knob_ugwp_doheat  = 1
  knob_ugwp_dokdis  = 1
  knob_ugwp_ndx4lh  = 1
  knob_ugwp_version = 0
  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