Skip to content

Commit

Permalink
Add changes for running v17p8_mynn
Browse files Browse the repository at this point in the history 
 - new field table, field_table_thompson_aero_tke_progsigma
 - updated namelist parameters in config.fcst
  • Loading branch information
@jkhender
jkhender committed Jan 11, 2024
1 parent 982a1f0 commit c4a41de
Show file tree
Hide file tree
Showing 2 changed files with 121 additions and 4 deletions.
55 changes: 51 additions & 4 deletions parm/config/gfs/config.fcst
View file
Original file line number Diff line number Diff line change
Expand Up @@ -46,6 +46,7 @@ export FCSTEXEC="ufs_model.x"
# Model configuration
export TYPE="nh"
export MONO="non-mono"
#JKHexport range_warn=".false." ## JKH

# Use stratosphere h2o physics
export h2o_phys=".true."
Expand Down Expand Up @@ -118,8 +119,21 @@ fi

# PBL/turbulance schemes
export hybedmf=".false."
export satmedmf=".true."
export isatmedmf=1
if [[ "$CCPP_SUITE" == "FV3_GFS_v17_p8_ugwpv1_mynn" || "$CCPP_SUITE" == "FV3_GFS_v17_p8_mynn" || "$CCPP_SUITE" == "FV3_GFS_v17_p8_c3_mynn" ]] ; then
export satmedmf=".false."
export isatmedmf=0
export shal_cnv=".false."
export do_mynnedmf=".true."
export do_mynnsfclay=".false."
export icloud_bl=1
export bl_mynn_tkeadvect=".true."
export bl_mynn_edmf=1
export bl_mynn_edmf_mom=1
export lcnorm=".true." ## JKH
else
export satmedmf=".true."
export isatmedmf=1
fi
tbf=""
if [[ "${satmedmf}" == ".true." ]]; then tbf="_satmedmf" ; fi

Expand Down Expand Up @@ -147,6 +161,32 @@ export doGP_lwscat=.false.
export iopt_sfc="3"
export iopt_trs="2"

# Convection Options: 2-SASAS, 3-GF
export progsigma=".true."
if [[ "$CCPP_SUITE" == "FV3_GFS_v17_p8_c3_mynn" ]] ; then
export imfdeepcnv=5
export imfshalcnv=-1 ## JKH - no shallow GF
elif [[ "$CCPP_SUITE" == "FV3_GFS_v17_p8_c3" ]] ; then
export progsigma=.false.
export imfdeepcnv=5
export imfshalcnv=5
elif [[ "$CCPP_SUITE" == "FV3_GFS_v17_p8_gf" ]] ; then
export imfdeepcnv=3
export imfshalcnv=3
else
export imfdeepcnv=2
if [[ "$CCPP_SUITE" == "FV3_GFS_v17_p8_ugwpv1_mynn" || "$CCPP_SUITE" == "FV3_GFS_v17_p8_mynn" ]] ; then
export imfshalcnv=-1
else
export imfshalcnv=2
fi
fi

#Convection schemes ### JKH - affects field table name
tbp=""
if [ "$progsigma" = ".true." ]; then tbp="_progsigma" ; fi


# Microphysics configuration
export dnats=0
export cal_pre=".true."
Expand All @@ -166,13 +206,20 @@ case ${imp_physics} in
;;
8) # Thompson
export ncld=2
export FIELD_TABLE="${HOMEgfs}/parm/ufs/fv3/field_table_thompson_noaero_tke${tbp}"
export nwat=6

if [[ "$CCPP_SUITE" == "FV3_GFS_v17_p8_ugwpv1_mynn" || "$CCPP_SUITE" == "FV3_GFS_v17_p8_mynn" || "$CCPP_SUITE" == "FV3_GFS_v17_p8_c3_mynn" ||
"$CCPP_SUITE" == "FV3_GFS_v17_p8_thompson" ]] ; then
export ltaerosol=".true."
export FIELD_TABLE="$HOMEgfs/parm/ufs/fv3/field_table_thompson_aero_tke${tbp}"
else
export ltaerosol=".false."
export FIELD_TABLE="${HOMEgfs}/parm/ufs/fv3/field_table_thompson_noaero_tke${tbp}"
fi

export cal_pre=".false."
export random_clds=".false."
export effr_in=".true."
export ltaerosol=".false."
export lradar=".true."
export ttendlim="-999"
export dt_inner=$((DELTIM/2))
Expand Down
70 changes: 70 additions & 0 deletions parm/ufs/fv3/field_table_thompson_aero_tke_progsigma
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,70 @@
# added by FRE: sphum must be present in atmos
# specific humidity for moist runs
"TRACER", "atmos_mod", "sphum"
"longname", "specific humidity"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=3.e-6" /
# prognostic cloud water mixing ratio
"TRACER", "atmos_mod", "liq_wat"
"longname", "cloud water mixing ratio"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=1.e30" /
# prognostic ice water mixing ratio
"TRACER", "atmos_mod", "ice_wat"
"longname", "cloud ice mixing ratio"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=1.e30" /
# prognostic rain water mixing ratio
"TRACER", "atmos_mod", "rainwat"
"longname", "rain water mixing ratio"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=1.e30" /
# prognostic snow water mixing ratio
"TRACER", "atmos_mod", "snowwat"
"longname", "snow water mixing ratio"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=1.e30" /
# prognostic Grau water mixing ratio
"TRACER", "atmos_mod", "graupel"
"longname", "graupel mixing ratio"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=1.e30" /
# prognostic cloud water number concentration
"TRACER", "atmos_mod", "water_nc"
"longname", "cloud liquid water number concentration"
"units", "/kg"
"profile_type", "fixed", "surface_value=0.0" /
# prognostic cloud ice number concentration
"TRACER", "atmos_mod", "ice_nc"
"longname", "cloud ice water number concentration"
"units", "/kg"
"profile_type", "fixed", "surface_value=0.0" /
# prognostic rain number concentration
"TRACER", "atmos_mod", "rain_nc"
"longname", "rain number concentration"
"units", "/kg"
"profile_type", "fixed", "surface_value=0.0" /
# prognostic ozone mixing ratio tracer
"TRACER", "atmos_mod", "o3mr"
"longname", "ozone mixing ratio"
"units", "kg/kg"
"profile_type", "fixed", "surface_value=1.e30" /
# water- and ice-friendly aerosols (Thompson)
"TRACER", "atmos_mod", "liq_aero"
"longname", "water-friendly aerosol number concentration"
"units", "/kg"
"profile_type", "fixed", "surface_value=0.0" /
"TRACER", "atmos_mod", "ice_aero"
"longname", "ice-friendly aerosol number concentration"
"units", "/kg"
"profile_type", "fixed", "surface_value=0.0" /
# prognostic subgrid scale turbulent kinetic energy
"TRACER", "atmos_mod", "sgs_tke"
"longname", "subgrid scale turbulent kinetic energy"
"units", "m2/s2"
"profile_type", "fixed", "surface_value=0.0" /
# prognotsitc sigmab tracer
"TRACER", "atmos_mod", "sigmab"
"longname", "sigma fraction"
"units", "fraction"
"profile_type", "fixed", "surface_value=0.0" /

0 comments on commit c4a41de

Please sign in to comment.