CCPP SciDoc  v6.0.0
Common Community Physics Package Developed at DTC
Supported Namelist Options

The SCM and the UFS Atmosphere, the atmospheric component of the UFS Weather Model, access runtime configurations from file input.nml. This file contains various namelists records that control aspects of the I/O, dynamics, physics etc. Most physics-related options are in records &gfs_physics_nml. Some schemes have their own namelist records as described below.

  • Namelist &gfs_physics_nml pertains to all of the suites used, but some of the variables are only relevant for specific parameterizations. Its variables are defined in file GFS_typedefs.F90 in the host model.
  • Namelist &gfdl_cloud_microphysics_nml is only relevant when the GFDL microphysics is used, and its variables are defined in module_gfdl_cloud_microphys.F90.
  • Namelist &cires_ugwp_nml specifies options for the use of CIRES Unified Gravity Wave Physics Version 0.
  • Namelist &nam_sfcperts specifies whether and how stochastic perturbations are used in the Noah Land Surface Model.
  • Namelist &nam_sppperts specifies options for the use of Stochastically-Perturbed Parameterizations (SPP)

Both the SDF and the input.nml contain information about how to specify the physics suite. Some of this information is redundant, and the user must make sure they are compatible.The safest practice is to use the SDF and namelist provided for each suite, since those are supported configurations. Changes to the SDF must be accompanied by corresponding changes to the namelist. While there is not a one-to-one correspondence between the namelist and the SDF, the tables below show some variables in the namelist that must match the SDF.

NML Description
Option CCPP scheme or interstitial Description Default Value
General options
fhzero gfs_phys_time_vary hour between clearing of diagnostic buckets 0.0
h2o_phys h2ophys flag for stratosphere h2o scheme .false.
ldiag3d see GFS_typedefs.F90 flag for 3D diagnostic fields .false.
qdiag3d see GFS_typedefs.F90 flag for 3D tracer diagnostic fields .false.
lssav see GFS_typedefs.F90 flag for storing diagnostics .false.
cplflx see GFS_typedefs.F90 flag for using fluxes provided by an external model .false.
cplwav see GFS_typedefs.F90 flag for using information produced by an external ocean wave model .false.
cplchm see GFS_typedefs.F90 flag for coupled chemistry diagnostics .false.
cplwav2atm see GFS_typedefs.F90 flag for wave to atm coupling .false.
oz_phys_2015 ozphys_2015 flag for new (2015) ozone physics .false.
fhcyc gfs_phys_time_vary frequency for surface data cycling in hours 0.0
use_ufo gfs_phys_time_vary flag for using unfiltered orography surface option .false.
ncld see GFS_typedefs.F90 number of hydrometeors 1
do_mynnsfclay see GFS_typedefs.F90 flag to activate MYNN-SFCLAY scheme .false.
do_sppt gfs_stochastics flag for stochastic SPPT option .false.
do_shum gfs_stochastics flag for stochastic SHUM option .false.
do_skeb gfs_stochastics flag for stochastic SKEB option .false.
do_sfcperts gfs_rrtmg_pre flag for stochastic surface perturbations option .false.
imp_physics choice of microphysics scheme choice of microphysics scheme:

  • 11: GFDL microphysics scheme
  • 8: Thompson microphysics scheme
  • 10: Morrison-Gettelman microphysics scheme
  • 17: NSSL microphysics scheme with background CCN
  • 18: NSSL microphysics scheme with predicted CCN (compatibility)
99
Parameters related to radiation scheme options
pdfcld gfs_rrtmg_pre flag for PDF clouds .false.
fhswr rrtmg_sw frequency for shortwave radiation (secs) 3600.
fhlwr rrtmg_lw frequency for longwave radiation (secs) 3600.
levr gfs_rrtmg_setup number of vertical levels for radiation calculations -99
nfxr gfs_rrtmg_pre second dimension of radiation input/output array fluxr 39+6
iflip gfs_rrtmg_setup control flag for vertical index direction
  • 0: index from TOA to surface
  • 1: index from surface to TOA
1
icliq_sw rrtmg_sw sw optical property for liquid clouds
  • 0: input cloud optical depth, ignoring iswcice setting
  • 1: cloud optical property scheme based on Hu and Stamnes (1993) [95] method
  • 2: cloud optical property scheme based on Hu and Stamnes (1993) [95] - updated
1
gfs_typedefs::gfs_control_type::iovr rrtmg_sw control flag for cloud overlap in SW & LW radiation
  • 0: random overlapping clouds
  • 1: max/ran overlapping clouds
  • 2: maximum overlap clouds (mcica only)
  • 3: decorrelation-length overlap (mcica only)
  • 4: exponential overlapping method
  • 5: exponential-random overlapping method
1
ictm gfs_rrtmg_setup external data time/date control flag
  • -2: same as 0, but superimpose seasonal cycle from climatology data set
  • -1: use user provided external data for the forecast time, no extrapolation
  • 0: use data at initial condition time, if not available, use latest, no extrapolation
  • 1: use data at the forecast time, if not available, use latest and extrapolation
  • yyyy0: use yyyy data for the forecast time, no further data extrapolation
  • yyyy1: use yyyy data for the forecast. if needed, do extrapolation to match the fcst time
1
crick_proof gfs_rrtmg_setup control flag for eliminating CRICK
  • .true.: apply layer smoothing to eliminate CRICK
  • .false.: do not apply layer smoothing
.false.
ccnorm gfs_rrtmg_setup control flag for in-cloud condensate mixing ratio
  • .true.: normalize cloud condensate
  • .false.: not normalize cloud condensate
.false.
norad_precip gfs_rrtmg_setup control flag for not using precip in radiation (Ferrier scheme)

  • .true.: snow/rain has no impact on radiation
  • .false.: snow/rain has impact on radiation
.false.
ialb gfs_rrtmg_setup SW surface albedo control flag:
  • 0: using climatology surface albedo scheme for SW
  • 1: using MODIS based land surface albedo for SW
0
iems gfs_rrtmg_setup LW surface emissivity control flag (ab 2-digit integer) :
  • a: =0 set surface air/ground t same for LW radiation
  • =1 set surface air/ground t diff for LW radiation
  • b: =0 use fixed surface emissivity = 1.0 (black-body)
  • =1 use varying climatology surface emissivity (veg based)
  • =2 future development (not yet)
0
iaer gfs_rrtmg_setup 4-digit aerosol flag (dabc for aermdl, volcanic, LW, SW):
  • d:tropospheric aerosol model scheme flag
    =0 or none, opac-climatology aerosol scheme
    =1 use gocart climatology aerosol scheme
    =2 use gocart prognostic aerosol scheme
    =5 opac-clim new spectral mapping
  • a:=0 use background stratospheric aerosol
    =1 include stratospheric volcanic aerosol
  • b:=0 no tropospheric aerosol in LW radiation
    =1 include tropospheric aerosol in LW
  • c:=0 no tropospheric aerosol in SW radiation
    =1 include tropospheric aerosol in SW
1
ico2 gfs_rrtmg_setup \(CO_2\) data source control flag:
  • 0: prescribed value (380 ppmv)
  • 1: yearly global averaged annual mean from observations
  • 2: monthly 15 degree horizontal resolution from observations
0
isubc_sw rrtmg_sw subgrid cloud approximation control flag in SW radiation:
  • 0: no McICA approximation in SW radiation
  • 1: use McICA with prescribed permutation seeds (test mode)
  • 2: use McICA with randomly generated permutation seeds
0
isubc_lw rrtmg_lw subgrid cloud approximation control flag in LW radiation:
  • 0: no McICA approximation in LW radiation
  • 1: use McICA with prescribed permutatition seeds (test mode)
  • 2: use McICA with randomly generated permutation seeds
0
isol gfs_rrtmg_setup solar constant scheme control flag:
  • 0: fixed value = 1366.0 \(W m^{-2}\) (old standard)
  • 10: fixed value = 1360.8 \(W m^{-2}\) (new standard)
  • 1: NOAA ABS-scale TSI table (yearly) with 11-yr cycle approximation
  • 2: NOAA TIM-scale TSI table (yearly) with 11-yr cycle approximation
  • 3: CMIP5 TIM-scale TSI table (yearly) with 11-yr cycle approximation
  • 4: CMIP5 TIM-scale TSI table (monthly) with 11-yr cycle approximation
0
lwhtr rrtmg_lw flag for output of longwave heating rate .true.
swhtr rrtmg_sw flag for output of shortwave heating rate .true.
nhfrad gfs_time_vary_pre number of timesteps for which to call radiation on physics timestep (coldstarts) 0
Parameters related to cumulus schemes
imfshalcnv choice of shallow convective scheme flag for mass flux shallow convective scheme:
  • 1:July 2010 version of mass-flux shallow convective scheme (operational as of 2016)
  • 2: scale- & aerosol-aware mass-flux shallow convective scheme (2017)
  • 3: scale- & aerosol-aware Grell-Freitas scheme
  • 4: new Tiedtke scheme (CAPS)
  • 0: modified Tiedtke's eddy-diffusion shallow convective scheme
  • -1: no shallow convection used
1
imfdeepcnv choice of deep convective scheme flag for mass-flux deep convective scheme:
  • -1: Chikira-Sugiyama deep convection (with cscnv = .T.)
  • 1: July 2010 version of SAS convective scheme (operational version as of 2016)
  • 2: scale- & aerosol-aware mass-flux deep convective scheme (2017)
  • 3: scale- & aerosol-aware Grell-Freitas scheme (GSD)
  • 4: new Tiedtke scheme (CAPS)
1
do_deep see GFS_typedefs.F90 consistency check for deep convection .true.
shal_cnv gfs_suite_interstitial flag for calling shallow convection .false.
lmfshal gfs_rrtmg_pre flag for mass-flux shallow convection scheme in the cloud fraction calculation shal_cnv .and. (imfshalcnv > 0)
lmfdeep2 gfs_rrtmg_pre flag for mass-flux deep convection scheme in the cloud fraction calculation imfdeepcnv == 2 .or. 3 .or.4
random_clds gfs_phys_time_vary flag for whether clouds are random .false.
trans_trac gfs_suite_interstitial flag for convective transport of tracers .false.
cal_pre gfs_phys_time_vary or gfs_MP_generic flag for calling precipitation type algorithm .false.
shcnvcw samfshalcnv flag for shallow convective cloud .false.
Parameters related to PBL scheme options
do_mynnedmf mynnedmf_wrapper flag to activate MYNN-EDMF scheme .false.
dspheat satmedmfvdifq flag for using TKE dissipative heating to temperature tendency in hybrid EDMF and TKE-EDMF schemes .false.
satmedmf satmedmfvdifq flag for calling scale-ware TKE-based EDMF PBL scheme .false.
isatmedmf satmedmfvdifq flag for scale-aware TKE-based moist EDMF scheme
  • 0: initial version of satmedmf (Nov.2018)
  • 1: updated version of satmedmf (as of May 2019)
0
do_ysu see GFS_typedefs.F90 flag for YSU PBL scheme .false.
debug see GFS_typedefs.F90 flag for debug printout .false.
xkzm_h satmedmfvdifq background vertical diffusion for heat and q 1.0d0
xkzm_m satmedmfvdifq background vertical diffusion for momentum 1.0d0
xkzm_s satmedmfvdifq sigma threshold for background mom. diffusion 1.0d0
dspfac satmedmfvdifq TKE dissipative heating factor 1.0
bl_upfr satmedmfvdifq updraft fraction in boundary layer mass flux scheme 0.13
bl_dnfr satmedmfvdifq downdraft fraction in boundary layer mass flux scheme 0.1
grav_settling mynnedmf_wrapper flag to activate gravitational settling of cloud droplets as described in Nakanishi (2000) [142] 0
bl_mynn_mixlength mynnedmf_wrapper flag for different version of mixing length formulation
  • 0: Original form from Nakanishi and Niino (2009) [141] . NO scale-awareness is applied to the master mixing length, regardless of "scaleware" setting
  • 1: HRRR operational form 201609-201807. Designed to work without the mass-flux scheme. Uses BouLac mixing length in free atmosphere.
  • 2: HRRR operational form 201807-present. Designed to be compatible with mass-flux scheme activated (default)
2
bl_mynn_edmf mynnedmf_wrapper flag to activate the mass-flux scheme
  • 0: deactivate mass-flux scheme
  • 1: activate dynamic multiplume mass-flux scheme

0
bl_mynn_edmf_mom mynnedmf_wrapper flag to activate the transport of momentum
  • 0: deactivate momentum transport in mass-flux scheme
  • 1: activate momentum transport in dynamic multiplume mass-flux scheme. bl_mynn_edmf must be set to 1
1
bl_mynn_edmf_tke mynnedmf_wrapper flag to activate the transport of TKE
  • 0: deactivate TKE transport in mass-flux scheme
  • 1: activate TKE transport in dynamic multiplume mass-flux scheme. bl_mynn_edmf must be set to 1
0
bl_mynn_edmf_part mynnedmf_wrapper flag to partitioning the MF and ED areas 0
bl_mynn_tkeadvect mynnedmf_wrapper activate computation of TKE advection (not yet in use for FV3)
  • false: deactivate TKE advection
  • true: activate TKE advection

.false.
bl_mynn_tkebudget mynnedmf_wrapper flag to activate TKE budget 0
bl_mynn_cloudpdf mynnedmf_wrapper flag to determine which cloud PDF to use
  • 0: use Sommeria-Deardorff subgrid cloud PDF
  • 1: use Kuwano-Yoshida subgrid cloud PDF
  • 2: use modified Chaboureau-Bechtold subgrid cloud PDF
2
bl_mynn_edmf_cloudmix mynnedmf_wrapper flag to activate mixing of cloud species
  • 0: deactivate the mixing of any water species mixing ratios
  • 1: activate the mixing of all water species mixing ratios
1
bl_mynn_mixqt mynnedmf_wrapper flag to mix total water or individual species
  • 0: mix individual water species separately
  • 1: DO NOT USE
0
icloud_bl mynnedmf_wrapper flag to coupling SGS clouds to radiation
  • 0: deactivate coupling subgrid clouds to radiation
  • 1: activate subgrid cloud coupling to radiation (highly suggested)
1
num_dfi_radar cu_gf_driver number of timespans with radar-prescried temperature tendencies 0
fh_dfi_radar cu_gf_driver begin+end of timespans to receive radar-prescribed temperature tendencies -2e10
do_cap_suppress cu_gf_driver enable convection suppression in GF scheme if fh_dfi_radar is specified .true.
ix_dfi_radar cu_gf_driver index within dfi_radar_tten of each timespan (-1 means "none") -1
dfi_radar_max_intervals cu_gf_driver number of radar-derived temperature tendency and/or convection suppression intervals 4
Parameters related to surface perturbation options
nsfcpert gfs_surface_generic_pre number of weights for stochastic surface perturbation 0
pertz0 gfs_surface_generic_pre magnitude of perturbation of momentum roughness length -999.
pertzt gfs_surface_generic_pre magnitude of perturbation of heat to momentum roughness length ratio -999.
pertshc gfs_surface_generic_pre magnitude of perturbation of soil hydraulic conductivity -999.
pertlai gfs_surface_generic_pre magnitude of perturbation of leaf area index -999.
pertalb gfs_surface_generic_pre magnitude of surface albedo perturbation -999.
pertvegf gfs_surface_generic_pre magnitude of perturbation of vegetation fraction -999.
Parameters related to Cellular Automata options
do_ca samfdeepcnv, gfs_stochastics cellular automata main switch .false.
ca_closure samfdeepcnv logical switch for CA on closure .false
ca_entr samfdeepcnv logical swith for CA on entrainment .false
ca_trigger samfdeepcnv logical switch for CA on trigger .false.
Parameters related to microphysics scheme options
lradar gfdl_cloud_microphys flag for computing radar reflectivity in Thompson MP scheme .false.
sedi_transport gfdl_cloud_microphys flag for turning on horizontal momentum transport during sedimentation .true.
do_sedi_w gfdl_cloud_microphys .true. to turn on vertical motion transport during sedimentation. (not supported in GFS physics) .false.
do_sedi_heat gfdl_cloud_microphys flag for turning on horizontal heat transport during sedimentation .true.
rad_snow gfdl_cloud_microphys flag for considering snow in cloud fraction calculation .true.
rad_graupel gfdl_cloud_microphys flag for considering graupel in cloud fraction calculation .true.
rad_rain gfdl_cloud_microphys flag for considering rain in cloud fraction calculation .true.
cld_min gfdl_cloud_microphys minimum cloud fraction. If total cloud condensate exceeds 1.0e-6 kg/kg, cloud fraction cannot be less than cld_min 0.05
const_vi gfdl_cloud_microphys flag for using constant cloud ice fall speed .false.
const_vs gfdl_cloud_microphys flag for using constant snow fall speed .false.
const_vg gfdl_cloud_microphys flag for using constant graupel fall speed .false.
const_vr gfdl_cloud_microphys flag for using constant rain fall speed .false.
vi_fac gfdl_cloud_microphys tunable factor for cloud ice fall or the constant cloud ice fall speed when const_vi is .true. 1.
vr_fac gfdl_cloud_microphys tunable factor for rain fall or the constant rain fall speed when const_vr is .true. 1.
vs_fac gfdl_cloud_microphys tunable factor for snow fall or the constant snow fall speed when const_vs is .true. 1.
vg_fac gfdl_cloud_microphys tunable factor for graupel fall or the constant graupel fall speed when const_vg is .true. 1.
vi_max gfdl_cloud_microphys maximum fall speed for cloud ice 0.5
vs_max gfdl_cloud_microphys maximum fall speed for snow 5.0
vg_max gfdl_cloud_microphys maximum fall speed for graupel 8.0
vr_max gfdl_cloud_microphys maximum fall speed for rain 12.0
qi_lim gfdl_cloud_microphys cloud ice limiter to prevent large ice built up in cloud ice freezing and deposition 1.
prog_ccn gfdl_cloud_microphys flag for activating prognostic CCN (not supported in GFS Physics) .false.
do_qa gfdl_cloud_microphys .true. to activate inline cloud fraction diagnosis in fast saturation adjustment. .false. to activate inline cloud fraction diagnosis in major cloud microphysics .true.
fast_sat_adj gfdl_cloud_microphys flag for adjusting cloud water evaporation (cloud water -> water vapor), cloud water freezing (cloud water -> cloud ice), cloud ice deposition (water vapor -> cloud ice) when fast saturation adjustment is activated (do_sat_adj = .true. in fv_core_nml block) .true.
tau_l2v gfdl_cloud_microphys time scale for evaporation of cloud water to water vapor. Increasing(decreasing) tau_l2v can decrease(boost) deposition of cloud water to water vapor 300.
tau_v2l gfdl_cloud_microphys time scale for condensation of water vapor to cloud water. Increasing(decreasing) tau_v2l can decrease(boost) condensation of water vapor to cloud water 150.
tau_g2v gfdl_cloud_microphys time scale for sublimation of graupel to water vapor. Increasing(decreasing) tau_g2v can decrease(boost) sublimation of graupel to water vapor 900.
tau_g2r gfdl_cloud_microphys time scale for graupel melting. Increasing(decreasing) tau_g2r can decrease(boost) melting of graupel to rain (graupel-> rain) 600.
tau_v2g gfdl_cloud_microphys time scale for deposition of water vapor to graupel. Increasing(decreasing) tau_v2g can decrease(boost) deposition of water vapor to graupel (water vapor -> graupel) 21600.
tau_l2r gfdl_cloud_microphys time scale for autoconversion of cloud water to rain. Increasing(decreasing) tau_l2r can decrese(boost) autoconversion of cloud water to rain (cloud water -> rain) 900.
tau_r2g gfdl_cloud_microphys time scale for freezing of rain to graupel. Increasing(decreasing) tau_r2g can decrease(boost) freezing of rain to graupel (rain->graupel) 900.
tau_i2s gfdl_cloud_microphys time scale for autoconversion of cloud ice to snow. Increasing(decreasing) tau_i2s can decrease(boost) autoconversion of cloud ice to snow (cloud ice -> snow) 1000.
tau_imlt gfdl_cloud_microphys time scale for cloud ice melting. Increasing(decreasing) tau_imlt can decrease(boost) melting of cloud ice to cloud water or rain (cloud ice -> cloud water or rain) 600.
tau_smlt gfdl_cloud_microphys time scale for snow melting. Increasing(decreasing) tau_smlt can decrease(boost) melting of snow to cloud water or rain (snow-> cloud water or rain) 900.
rthresh gfdl_cloud_microphys critical cloud water radius for autoconversion (cloud water -> rain). Increasing(decreasing) of rthresh makes the autoconversion harder(easier) 10.0e-6
dw_land gfdl_cloud_microphys base value for subgrid deviation/variability over land 0.20
dw_ocean gfdl_cloud_microphys base value for subgrid deviation/variability over ocean 0.10
ql_gen gfdl_cloud_microphys maximum value for cloud water generated from condensation of water vapor (water vapor-> cloud water) 1.0e-3
qi_gen gfdl_cloud_microphys maximum value of cloud ice generated from deposition of water vapor (water vapor->cloud ice) or freezing(cloud water -> cloud ice). Increasing(decreasing) qi_gen can increas(decrease) cloud ice 1.82e-6
ql_mlt gfdl_cloud_microphys maximum value of cloud water allowed from melted cloud ice (cloud ice -> cloud water or rain). Exceedance of which will become rain. Increasing(decreasing) ql_mlt can increase(decrease) cloud water and decrease(increase) rain 2.0e-3
qs_mlt gfdl_cloud_microphys maximum value of cloud water allowed from melted snow (snow -> cloud water or rain). Exceedance of which will become rain. Increasing(decreasing) qs_mlt can increas(decrease) cloud water and decrease (increase) rain 1.0e-6
ql0_max gfdl_cloud_microphys threshold of cloud water to rain autoconversion (cloud water -> rain). Increasing(decreasing) ql0_max can increase(decrease) rain and decrease(increase) cloud water 2.0e-3
qi0_max gfdl_cloud_microphys maximum value of cloud ice generated from other sources like convection. Exceedance of which will become snow. Increasing(decreasing) qi0_max can increase(decrease) cloud ice and decrease(increase) snow 1.0e-4
qi0_crt gfdl_cloud_microphys threshold of cloud ice to snow autoconversion (cloud ice -> snow). Increasing(decreasing) qi0_crt can increase(decrease) cloud ice and decrease(increase) snow 1.0e-4
qs0_crt gfdl_cloud_microphys threshold of snow to graupel autoconversion (snow -> graupel). Increasing(decreasing) qs0_crt can increase(decrease) snow and decrease(increase) graupel 1.0e-3
qc_crt gfdl_cloud_microphys minimum value of cloud condensate to allow partial cloudiness. Partial cloud can only exist when total cloud condensate exceeds qc_crt 5.0e-8
c_psaci gfdl_cloud_microphys accretion efficiency of cloud ice to snow (cloud ice -> snow). Increasing(decreasing) of c_psaci can boost(decrease) the accretion of cloud ice to snow 0.02
c_pgacs gfdl_cloud_microphys accretion efficiency of snow to graupel (snow -> graupel). Increasing(decreasing) of c_pgacs can boost(decrease) the accretion of snow to graupel 2.0e-3
rh_inc gfdl_cloud_microphys relative humidity increment for complete evaporation of cloud water and cloud ice 0.25
rh_inr gfdl_cloud_microphys relative humidity increment for sublimation of snow 0.25
rh_ins gfdl_cloud_microphys relative humidity increment for minimum evaporation of rain 0.25
rthresh gfdl_cloud_microphys critical cloud water radius for autoconversion(cloud water->rain). Increasing(decreasing) of rthresh makes the autoconversion harder(easier) 1.0e-5
ccn_l gfdl_cloud_microphys base CCN over land. Increasing(decreasing) ccn_l can on the one hand boost(decrease) the autoconversion of cloud water to rain, on the other hand make the autoconversion harder(easier). The unit is \(cm^{-3}\) 270.
ccn_o gfdl_cloud_microphys base CCN over ocean. Increasing(decreasing) ccn_o can on the one hand boost(decrease) the autoconversion of cloud water to rain, on the other hand make the autoconversion harder(easier). The unit is \(cm^{-3}\) 90.
c_paut gfdl_cloud_microphys autoconversion efficiency of cloud water to rain (cloud water -> rain). Increasing(decreasing) of c_paut can boost(decrease) the autoconversion of cloud water to rain 0.55
c_cracw gfdl_cloud_microphys accretion efficiency of cloud water to rain (cloud water -> rain). Increasing(decreasing) of c_cracw can boost(decrease) the accretion of cloud water to rain 0.9
sat_adj0 gfdl_cloud_microphys adjust factor for condensation of water vapor to cloud water (water vapor->cloud water) and deposition of water vapor to cloud ice 0.9
use_ppm gfdl_cloud_microphys true to use PPM fall scheme; false to use time-implicit monotonic fall scheme .false.
use_ccn gfdl_cloud_microphys true to compute prescribed CCN. It should be .true. when prog_ccn = .false. .false.
mono_prof gfdl_cloud_microphys true to turn on terminal fall with monotonic PPM scheme. This is used together with use_ppm=.true. .true.
z_slope_liq gfdl_cloud_microphys true to turn on vertically subgrid linear monotonic slope for autoconversion of cloud water to rain .true.
z_slope_ice gfdl_cloud_microphys true to turn on vertically subgrid linear monotonic slope for autoconversion of cloud ice to snow .false.
de_ice gfdl_cloud_microphys true to convert excessive cloud ice to snow to prevent ice over-built from other sources like convection scheme (not supported in GFS physics) .false.
fix_negative gfdl_cloud_microphys true to fix negative water species using nearby points .false.
icloud_f gfdl_cloud_microphys flag (0,1,or 2) for cloud fraction diagnostic scheme 0
irain_f gfdl_cloud_microphys flag (0 or 1) for cloud water autoconversion to rain scheme. 0: with subgrid variability; 1: no subgrid variability 0
mp_time gfdl_cloud_microphys time step of GFDL cloud microphysics (MP). If mp_time isn't divisible by physics time step or is larger than physics time step, the actual MP time step becomes dt/NINT[dt/MIN(dt,mp_time)] 150.
alin gfdl_cloud_microphys parameter a in Lin et al.(1983). Constant in empirical formula for \(U_R\). Increasing(decreasing) alin can boost(decrease) accretion of cloud water by rain and rain evaporation 842.
clin gfdl_cloud_microphys parameter c in Lin et al.(1983). Constant in empirical formula for \(U_S\). Increasing(decreasing) clin can boost(decrease) accretion of cloud water by snow, accretion of cloud ice by snow, snow sublimation and deposition, and snow melting 4.8
t_min gfdl_cloud_microphys temperature threshold for instant deposition. Deposit all water vapor to cloud ice when temperature is lower than t_min 178.
t_sub gfdl_cloud_microphys temperature threshold for sublimation. Cloud ice, snow or graupel stops(starts) sublimation when temperature is lower(higher) then t_sub 184.
mp_print gfdl_cloud_microphys .true. to turn on GFDL cloud microphysics debugging print out. (not supported in GFS physics) .false.
ltaerosol mp_thompson flag for using aerosol climotology in Thompson MP scheme .false.
ttendlim mp_thompson temperature tendency limiter per time step in K/s, set to < 0 to deactivate -999.0
ext_diag_thompson mp_thompson flag for extended diagnostic output from Thompson MP .false.
thompson_ext_ndiag3d mp_thompson number of 3d arrays for extended diagnostic output from Thompson MP 37
dt_inner mp_thompson time step for the inner loop in second -999.0
sedi_semi mp_thompson flag for semi Lagrangian sedi of rain .false.
decfl mp_thompson deformed CFL factor 8
effr_in gfdl_cloud_microphys, mp_thompson flag for using input cloud effective radii calculation .false.
cnvcld see GFS_typedefs.F90 flag for convective cloud .false.
lgfdlmprad gfs_rrtmg_pre flag for GFDL mp scheme and radiation consistency .false.
nssl_cccn mp_nssl CCN concentration (m^-3) 0.6e9
nssl_alphah mp_nssl graupel shape parameter 0.0
nssl_alphahl mp_nssl hail shape parameter 1.0
nssl_hail_on mp_nssl NSSL flag to activate the hail category .false.
nssl_ccn_on mp_nssl NSSL flag to activate the CCN category .true.
nssl_invertccn mp_nssl NSSL flag to treat CCN as activated or unactivated .true.
Parameters related to gravity drag scheme options
knob_ugwp_version cires_ugwp parameter selects a version of the UGWP implementation in FV3GFS-127L

  • 0: default version delivered to EMC in Jan 2019 for implementation
  • 1: version of UGWP under development that plans to consider the physics-based sources of NGWs (knob_ugwp_wvspec [2:4]), options for stochastic and deterministic excitation of waves (knob_ugwp_stoch), and switches between different UGWP schemes (knob_ugwp_solver)
0
knob_ugwp_doaxyz cires_ugwp parameter controls application of the momentum deposition for NGW-schemes
  • 0: the momentum tendencies due to NGWs are calculated, but tendencies do not change the horizontal winds
  • 1: default value; it changes the horizontal momentum tendencies and horizontal winds
1
knob_ugwp_doheat cires_ugwp parameter controls application of the heat deposition for NGW-schemes
  • 0: the temperature tendencies due to NGWs are calculated but tendencies do not change the temperature state
  • 1: default value; it changes the temperature tendencies and kinetic temperature
1
knob_ugwp_dokdis cires_ugwp parameter controls application of the eddy diffusion due to instability of NGWs
  • 0: the eddy diffusion tendencies due to NGWs are calculated but tendencies do not change the model state vector
  • 1: it computes eddy diffusion coefficient due to instability of NGWs; in UGWP v0, eddy viscosity, heat conductivity and tracer diffusion are not activated
0
knob_ugwp_solver cires_ugwp parameter controls the selection of UGWP-solvers(wave propagation, dissipation and wave breaking) for NGWs
  • 1: represents the discrete multi-wave solver with background dissipation and linear wave saturation
  • 2: represents the spectral deterministic solver with background dissipation and spectral saturation
  • 3: represents the discrete multi-wave solver with the background dissipation, extension of Alexander sand Dunkerton (1999)
  • 4: represents the spectral solver with background dissipation, extension of Doppler Spread Theory of Hines (1997)
1
knob_ugwp_ndx4lh cires_ugwp parameter controls the selection of the horizontal wavenumber(wavelength) for NGW schemes
  • 1: selects the \(4xdx\) sub-grid wavelength, where dx is the horizontal resolution of the model configuration (C96-400km; C768-52km)
2
knob_ugwp_wvspec cires_ugwp four-dimensional array defines number of waves in each arimuthal propagation (as defined by knob_ugwp_azdir) for GWs excited due to the following four sources:
(1) sub-grid orography (knob_ugwp_wvspec[1]=1),
(2) convective (knob_ugwp_wvspec[2]=25),
(3) frontal (knob_ugwp_wvspec[3]=25) activity,
(4) knob_ugwp_wvspec[4] represents number of wave excited by dynamical imbalances that may mimic both convective and front-jet mechanisms of GW triggering.
In UGWP v0, first two elements of the array, knob_ugwp_wvspec(1:2), control number of waves for stationary (OGW) and nonstationary waves (NGWs). 		1,32,32,32
knob_ugwp_azdir cires_ugwp four-dimensional array that defines number of azimuths for propagation of GWs triggered by four types of physics-based sources (orography, convection, front-jets, and dynamical imbalance). In UGWP v0, first two elements of the array, knob_ugwp_azdir(1:2), control number of azimuths for OGW and NGWs respectively.
2,4,4,4
knob_ugwp_stoch cires_ugwp four-dimensional array that control stochastic selection of GWs triggered by four types of physics-based sources.
Default values:0,0,0,0 - reflect determinstic selection of GW parameters without stochastic selection
0,0,0,0
knob_ugwp_effac cires_ugwp four-dimensional array that control efficiency of GWs triggerd by four types of physics-based sources.
Default values: 1.,1.,1.,1. - reflect that calculated GW-tendencies will be applied for the model state. 		1.,1.,1.,1.
launch_level cires_ugwp parameter has been introduced by EMC during implementation. It defines the interface model level from the surface at which NGWs are launched.
Default value for FV3GFS-64L, launch_level=25 and for FV3GFS-128L, launch_level=52. 		55
ldiag_ugwp cires_ugwp flag for CIRES UGWP diagnostics .false.
do_ugwp cires_ugwp flag for CIRES UGWP revised OGW
  • .T.: revised gwdps_v0
  • .F.: GFS operational orographic gwdps
.false.; The CIRES Unified Gravity Wave Physics (cires_ugwp) scheme is used in GFSv15p2 and GFSv16beta SDFs with do_ugwp=F in the namelist. In this setting, the cires_ugwp calls the operational GFS v15.2 orographic gravity wave drag (gwdps) scheme. When do_ugwp=.T., the cires_ugwp scheme calls an experimental orographic gravity wave (gwdps_v0)
do_tofd cires_ugwp flag for turbulent orographic form drag .false.
cnvgwd cires_ugwp flag for convective gravity wave drag scheme dependent on maxval(cdmbgwd(3:4) == 0.0) .false.
cdmbgwd(4) cires_ugwp multiplication factors for mountain blocking(1), orographic gravity wave drag(2)
  • [1]: GWDPS mountain blocking
  • [2]: GWDPS orographic gravity wave drag
  • [3]: the modulation total momentum flux of NGWs by intensities of the total precipitation
  • [4]: TKE for future tests and applications
2.0,0.25,1.0,1.0
nmtvr cires_ugwp number of topographic variables such as variance etc used in the GWD parameterization-10 more added if GSL orographic drag scheme is used 14
cgwf cires_ugwp, unified_ugwp multiplication factor for convective GWD 0.5d0,0.05d0
do_gwd see GFS_typedefs.F90 flag for gravity wave drag maxval(cdmbgwd) > 0.0
gwd_opt drag_suite flag for GWD scheme
  • 1: original GFS GWD
  • 2: Unified UGWP GWD
  • 22: Unified UGWP GWD with extra output
  • 3: GSL drag suite
  • 33: GSL drag suite with extra output
1
do_ugwp_v0 unified_ugwp flag for version 0 UGWP GWD .true.
do_ugwp_v0_orog_only unified_ugwp flag for version 0 UGWP GWD (orographic drag only) .false.
do_ugwp_v0_nst_only unified_ugwp flag for version 0 UGWP GWD (non-stationary GWD only) .false.
do_gsl_drag_ls_bl unified_ugwp, drag_suite flag for GSL drag (large-scale GWD and blocking only) .false.
do_gsl_drag_ss unified_ugwp, drag_suite flag for GSL drag (small-scale GWD only) .false.
do_gsl_drag_tofd unified_ugwp, drag_suite flag for GSL drag (turbulent orog form drag only) .false.
Parameters related to LSM options
lsm see GFS_typedefs.F90 flag for land surface model to use
  • 1: Noah LSM
  • 2: NoahMP LSM
  • 3: RUC LSM
1
lsoil lsm_noah number of soil layers 4
rdlai lsm_ruc flag to read leaf area index from input files .false.
ivegsrc lsm_noah, lsm_ruc, noahmpdrv, sfc_diff flag for vegetation type dataset choice:
  • 0: USGS
  • 1: IGBP(20 category): IGBP must be selected if NoahMP is used
  • 2: UMD (13 category)
2
isot lsm_noah, lsm_ruc, noahmpdrv flag for soil type dataset choice:
  • 0: Zobler soil type (9 category)
  • 1: STATSGO soil type (19 category): STATSGO must be selected if NoahMP is used

0

iopt_dveg noahmpdrv options for dynamic vegetation
  • 1: off (use table LAI; use FVEG = SHDFAC from input)
  • 2: on (together with iopt_crs = 1)
  • 3: off (use table LAI; calculate FVEG)
  • 4: off (use table LAI; use maximum vegetation fraction)
  • 5: on (use maximum vegetation fraction)
  • 6: on (use FVEG = SHDFAC from input)
  • 7: off (use input LAI; use FVEG = SHDFAC from input)
  • 8: off (use input LAI; calculate FVEG)
  • 9: off (use input LAI; use maximum vegetation fraction)
4
iopt_crs noahmpdrv options for canopy stomatal resistance
  • 1: Ball-Berry
  • 2: Jarvis
1
iopt_btr noahmpdrv options for soil moisture factor for stomatal resistance
  • 1: Noah (soil moisture)
  • 2: CLM (matric potential)
  • 3: SSIB (matric potential)
1
iopt_run noahmpdrv options for runoff and groundwater
  • 1: TOPMODEL with groundwater (Niu et al. 2007 [146])
  • 2: TOPMODEL with an equilibrium water table (Niu et al. 2005 [145])
  • 3: original surface and subsurface runoff (free drainage)
  • 4: bats surface and subsurface runoff (free drainage)
  • 5: Miguez-macho&Fan groundwater scheme (Miguez-Macho et al. 2007 [133]; Fan et al. 2007 [58]; needs further testing for public use)
3
iopt_sfc noahmpdrv options for surface layer drag coeff (CH&CM)
  • 1: m-o
  • 2: original Noah (Chen et al. 1997 [36])
1
iopt_frz noahmpdrv options for supercooled liquid water (or ice fraction)
  • 1: no interation (Niu and Yang (2006) [144] )
  • 2: Koren's iteration
1
iopt_inf noahmpdrv options for frozen soil permeability
  • 1: linear effects, more permeable (Niu and Yang (2006) [144])
  • 2: nonlinear effects, less permeable (old)
1
iopt_rad noahmpdrv options for radiation transfer
  • 1: modified two-stream (gap = f(solar angle, 3d structure ...)<1-fveg)
  • 2: two-stream applied to grid-cell (gap = 0)
  • 3: two-stream applied to vegetated fraction (gap=1-FVEG)
3
iopt_alb noahmpdrv options for ground snow surface albedo
  • 1: BATS
  • 2: CLASS
2
iopt_snf noahmpdrv options for partitioning precipitation into rainfall and snowfall
  • 1: Jordan (1991)
  • 2: BATS: when sfctmp < tfrz+2.2
  • 3: sfctmp < tfrz
  • 4: use WRF microphysics output
1
iopt_tbot noahmpdrv options for lower boundary condition of soil temperature
  • 1: zero heat flux from bottom (zbot and tbot not used)
  • 2: tbot at zbot (8m) read from a file (original Noah)
2
iopt_stc noahmpdrv options for snow/soil temperature time scheme (only layer 1)
  • 1: semi-implicit; flux top boundary condition
  • 2: full implicit (original Noah); temperature top boundary condition
  • 3: same as 1, but fsno for ts calculation (generally improve snow; v3.7)
1
iopt_trs noahmpdrv options for thermal roughness scheme:
  • 1: z0h=z0m
  • 2: czil
  • 3: ec
  • 4: kb inversed
2
Parameters related to other surface scheme options
nstf_name(5) sfc_nst NSST related paramters:

  • nstf_name(1): 0=NSST off, 1= NSST on but uncoupled, 2= NSST on and coupled
  • nstf_name(2): 1=NSST spin up on, 0=NSST spin up off
  • nstf_name(3): 1=NSST analysis on, 0=NSST analysis off
  • nstf_name(4): zsea1 in mm
  • nstf_name(5): zesa2 in mm
/0,0,1,0,5/
nst_anl gfs_phys_time_vary flag for NSST analysis in gcycle/sfcsub .false.
frac_grid fractional grid flag for fractional grid .false.
min_lakeice fractional grid minimum lake ice value 0.15d0
min_seaice fractional grid minimum sea ice value 1.0d-11
min_lake_height fractional grid minimum lake height value 250.0
sfc_z0_type sfc_diff surface roughness options over ocean
  • 0: no change
  • 6: areodynamical roughness over water with input 10-m wind
  • 7: slightly decrease Cd for higher wind speed compared to 6
  • negative when cplwav2atm=.true. - i.e. two way wave coupling
0
redrag sfc_diff flag for applying reduced drag coefficient for high wind over sea in GFS surface layer scheme .false.
lheatstrg gfs_surface_generic_post flag for canopy heat storage parameterization .false.
z0fac gfs_surface_generic_post surface roughness fraction factor 0.3
e0fac gfs_surface_generic_post latent heat flux fraction factor relative to sensible heat flux,e.g., e0fac=0.5 indicates that canopy heat storage for latent heat flux is 50% of that for sensible heat flux 0.5