Overview

The HAFS_thompson suite is one of the physics suites in the operational HAFS, currently applied in the National Hurricane Center's Northern Atlantic basin. Wang et al. (2024) [135] documented the details of physics schemes used in HAFS, including a few TC-specific modifications such as roughness length type and PBL scheme which can be configured through the namelist options.

The HAFS_thompson suite uses the parameterizations in the following order:

Suite Definition File

<?xml version="1.0" encoding="UTF-8"?>
 
<suite name="FV3_HAFS_v1_thompson_nonsst" version="1">
  <!-- <init></init> -->
  <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_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="phys_ps">
    <subcycle loop="1">
      <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_ocean</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>unified_ugwp</scheme>
      <scheme>unified_ugwp_post</scheme>
      <scheme>GFS_GWD_generic_post</scheme>
      <scheme>GFS_suite_stateout_update</scheme>
    </subcycle>
  </group>
  <group name="phys_ts">
    <subcycle loop="1">
      <scheme>GFS_photochemistry</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>mp_thompson_pre</scheme>
      <scheme>mp_thompson</scheme>
      <scheme>mp_thompson_post</scheme>
      <scheme>GFS_MP_generic_post</scheme>
      <scheme>maximum_hourly_diagnostics</scheme>
      <scheme>GFS_physics_post</scheme>
    </subcycle>
  </group>
  <group name="stochastics">
    <subcycle loop="1">
      <scheme>GFS_stochastics</scheme>
    </subcycle>
  </group>
  <!-- <finalize></finalize> -->
</suite>

Namelist

General physics options
&gfs_physics_nml

fhzero = 3.

ldiag3d = .false.

lradar = .true.

avg_max_length = 10800.

h2o_phys = .true.

fhcyc = 0.

use_ufo = .true.

pre_rad = .false.

imp_physics = 8

iovr = 3

ltaerosol = .false.

ttendlim = -999

dt_inner = 45

sedi_semi = .true.

decfl = 10

pdfcld = .false.

min_lakeice = 0.15

min_seaice = 0.15

use_cice_alb = .false

fhswr = 900.

fhlwr = 900.

ialb = 1

iems = 1

iaer = 1011

ico2 = 2

isubc_sw = 2

isubc_lw = 2

isol = 2

lwhtr = .true.

swhtr = .true.

cnvgwd = .false.

cal_pre = .false.

redrag = .true.

dspheat = .true.

hybedmf = .false.

hurr_pbl = .true. ! HWRF moninedmf

moninq_fac = -1.0 ! HWRF_moninedmf

satmedmf = .true.

isatmedmf = 1

lseaspray = .false.

rlmx = 300.

elmx = 300.

sfc_rlm = 1

tc_pbl = 0

random_clds = .false.

trans_trac = .true.

cnvcld = .true.

shal_cnv = .true. !Shallow convection

do_deep = .true.

imfshalcnv = 2

imfdeepcnv = 2

clam_deep = 0.15

betal_deep = 0.05

betas_deep = 0.05

c1_deep = 0.002

c1_shal = 0.0005

progsigma = .false.

betascu = 4.0

betamcu = 1.0

betadcu = 1.0

cdmbgwd = 1.0,1.0,1.0,1.0

sfc_z0_type = 6

prslrd0 = 0.

ivegsrc = 1

isot = 1

lsoil = 4

lsm = 1

iopt_dveg = 4

iopt_crs = 2

iopt_btr = 1

iopt_run = 1

iopt_sfc = 3

iopt_frz = 1

iopt_inf = 1

iopt_rad = 3

iopt_alb = 1

iopt_snf = 4

iopt_tbot = 2

iopt_stc = 3

iopt_trs = 2

lsoil_lsm = 4

debug = .false.

oz_phys = .false.

oz_phys_2015 = .true.

nst_anl = .true.

nstf_name = 2,1,0,0,0

cplflx = .false.

cplice = .false.

cplocn2atm = .true.

icplocn2atm = 1

cplwav = .false.

cplwav2atm = .false.

cpl_imp_mrg = .true.

psautco = 0.0008, 0.0005

prautco = 0.00015, 0.00015

iau_delthrs = 6

iaufhrs = 0

iau_inc_files = ,

iau_drymassfixer = .false.

lgfdlmprad = .false.

effr_in = .true.

do_sppt = .false.

do_shum = .false.

do_skeb = .false.

lndp_type = 0

n_var_lndp = 0

do_ca = .false.

ca_sgs = .false.

nca = 1

ncells = 5

nlives = 12

nseed = 1

nfracseed = 0.5

ca_trigger = .false.

ca_entr = .false.

ca_closure = .false.

ca_global = .false.

nca_g = 1

ncells_g = 1

nlives_g = 100

nseed_g = 100

ca_smooth = .false.

nspinup = 1

iseed_ca = 1

nsmooth = 100

ca_amplitude = 0.35

print_diff_pgr = .true.

do_ugwp = .false.

do_tofd = .false.

gwd_opt = 2

do_ugwp_v0 = .false.

do_ugwp_v1 = .false.

do_ugwp_v0_orog_only = .false.

do_ugwp_v0_nst_only = .true.

do_ugwp_v1_w_gsldrag = .false.

do_ugwp_v1_orog_only = .false.

do_gsl_drag_ls_bl = .true.

do_gsl_drag_ss = .true.

do_gsl_drag_tofd = .true.

/