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
• CIRES Unified Gravity Wave Physics Scheme - Version 0
• GFS Rayleigh Damping Calculation
• GFS Ozone Photochemistry (2015) Scheme
• GFS Stratospheric H2O Scheme
• GFS SAS-based Mass-Flux Deep Convection Scheme
• GFS SAS-based Mass-Flux Shallow Convection Scheme
• GFDL Cloud Microphysics Scheme
• GFS Precipitation Type Diagnosis Scheme

Suite Definition File

• For NEMSIO initialization data: suite_FV3_GFS_v16beta.xml

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 [1.1,0.72,1.0,1.0] [0.23,1.5,1.0,1.0] [0.14,1.8,1.0,1.0] ! [C768] [C384] [C192] [C96]L64
  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

Note

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