Add a new E-epsilon PBL scheme: (#1303)

TYPE: new feature KEYWORDS: E-epsilon, EEPS, TKE, PBL SOURCE: Chunxi Zhang: 1: (past) International Pacific Research Center, University of Hawai‘i at Mānoa 2: (past) Center for Analysis and Prediction of Storms, University of Oklahoma 3: (current) I. M. System Group, Inc. (IMSG); NOAA/EMC Yuqing Wang International Pacific Research Center, and Department of Atmospheric Sciences, University of Hawai‘i at Mānoa DESCRIPTION OF CHANGES: A turbulence kinetic energy (TKE) and TKE dissipation rate (ε) based 1.5-order closure PBL parameterization is added: 1) the file module_bl_eepsilon.F includes the module for this scheme 2) standard procedures for a new scheme in WRF: changed init, addtendc and driver 3) three more variables are declared in registry.EM_COMMON. e.g., pep_pbl (TKE dissipation rate, state) pek_adv (TKE advection, scalar) pep_adv (PEP advection, scalar) the files module_first_rd_step_part1.F and start_em.F are revised accordingly this scheme is denfied as eepsscheme bl_pbl_physics == 16 4) add module_bl_eepsilon in Makefile LIST OF MODIFIED FILES: phys/module_bl_eepsilon.F phys/module_pbl_driver.F phys/module_physics_init.F phys/module_physics_addtendc.F phys/Makefile Registry/Registry.EM_COMM dyn_em/module_first_rk_step_part1.F dyn_em/start_em.F run/README.namelist TESTS CONDUCTED: 1. Automated jenkins testing is all PASS. 2. The performances of the newly implemented EEPS scheme and the existing Yonsei University (YSU) scheme, the University of Washington (UW) scheme, and Mellor–Yamada–Nakanishi–Niino (MYNN) scheme are evaluated over the stratocumulus dominated southeast Pacific (SEP) and over the Southern Great Plains (SGP) where strong PBL diurnal variation is common. The simulations by these PBL parameterizations are compared with various observations from two field campaigns: the Variability of American Monsoon Systems Project (VAMOS) Ocean–Cloud–Atmosphere–Land Study (VOCALS) in 2008 over the SEP and the Land–Atmosphere Feedback Experiment (LAFE) in 2017 over the SGP. Results show that the EEPS and YSU schemes perform comparably over both regions, while the MYNN scheme performs differently in many aspects, especially over the SEP. The EEPS (MYNN) scheme slightly (significantly) underestimates liquid water path over the SEP. Compared with observations, the UW scheme produces the best PBL height over the SEP. The MYNN produces too high PBL height over the western part of the SEP while both the YSU and EEPS schemes produce too low PBL and cloud-top heights. The differences among the PBL schemes in simulating the PBL features over the SGP are relatively small. 3. See published paper: Zhang, C., Y. Wang, and M. Xue, 2020: Evaluation of an E–ε and Three Other Boundary Layer Parameterization Schemes in the WRF Model over the Southeast Pacific and the Southern Great Plains. Mon. Wea. Rev., 148, 1121–1145, https://doi.org/10.1175/MWR-D-19-0084.1. RELEASE NOTE: A turbulence kinetic energy (TKE) and TKE dissipation rate (ε) based 1.5-order closure PBL parameterization (E–ε, EEPS) is added (Zhang et al. 2020, MWR). Works with surface layer options, 1, 91, 2 and 5. (Thanks to Chunxi Zhang and Yuqing Wang of University of Hawaii.)
wrf-model · Feb 26, 2021 · 1430495 · 1430495
1 parent 0fff69a
commit 1430495
Show file tree

Hide file tree

Showing 9 changed files with 1,115 additions and 5 deletions.
diff --git a/Registry/Registry.EM_COMMON b/Registry/Registry.EM_COMMON
@@ -1027,6 +1027,11 @@ state   real   ZOL            ij     misc        1         -     -         "ZOL"
 state   real   WSTAR_YSU      ij     misc        1         -     -        "WSTAR_YSU"           "mixed-layer velocity scale from ysupbl" "m/s"
 state   real   DELTA_YSU      ij     misc        1         -     -        "DELTA_YSU"           "entrainment layer depth from ysupbl"    "m"
 
+# EEPS PBL  scheme
+state   real   pep_pbl        ikj    misc        1         -     hr        "PEP_PBL"            "TKE dissipation rate from EEPS"    "m2 s-3"
+state   real   pek_adv        ikjftb  scalar      1         -      i0rusdf=(bdy_interp:dt) "pek_adv"       "TKE from EEPS"      "m2 s-2"
+state   real   pep_adv        ikjftb  scalar      1         -      i0rusdf=(bdy_interp:dt) "pep_adv"       "TKE dissipating rate from EEPS"      "m2 s-3"
+
 # MYJ PBL variables; GBM PBL: EXCH_H, EXCH_M
 state    real   EXCH_H          ikj     misc        1         Z     r         "EXCH_H"               "SCALAR EXCHANGE COEFFICIENTS "                "m2 s-1"
 state    real   EXCH_M          ikj     misc        1         Z     r         "EXCH_M"               "EXCHANGE COEFFICIENTS "                       "m2 s-1"
@@ -3016,6 +3021,7 @@ package   camuwpblscheme bl_pbl_physics==9           -             state:tauresx
 package   temfpblscheme  bl_pbl_physics==10          -             state:te_temf,kh_temf,km_temf,shf_temf,qf_temf,uw_temf,vw_temf,wupd_temf,mf_temf,thup_temf,qlup_temf,qtup_temf,cf3d_temf,hd_temf,lcl_temf,hct_temf,cfm_temf
 package   shinhongscheme bl_pbl_physics==11          -             state:el_pbl,tke_pbl
 package   gbmpblscheme   bl_pbl_physics==12          -             state:exch_tke,el_pbl,tke_pbl
+package   eepsscheme     bl_pbl_physics==16          -             scalar:pek_adv,pep_adv;state:tke_pbl,pep_pbl
 package   mrfscheme      bl_pbl_physics==99          -             -
 
 package   mynn_tkebudget bl_mynn_tkebudget==1        -             state:qSHEAR,qBUOY,qDISS,qWT,dqke

diff --git a/dyn_em/module_first_rk_step_part1.F b/dyn_em/module_first_rk_step_part1.F
@@ -1055,6 +1055,10 @@ SUBROUTINE first_rk_step_part1 (   grid , config_flags              &
      &        ,spp_pbl=config_flags%spp_pbl                               &
      &        ,pattern_spp_pbl=grid%pattern_spp_pbl                       &
      &        ,restart=config_flags%restart,cycling=config_flags%cycling  &
+!EEPS for ARW
+     &        ,pep=grid%pep_pbl                                           &
+     &        ,PEK_ADV=scalar(ims,kms,jms,P_pek_adv)                      &!TKEadvection
+     &        ,PEP_ADV=scalar(ims,kms,jms,P_pep_adv)                      &!TKEadvection
 !GWD for ARW
      &        ,GWD_OPT=config_flags%gwd_opt &
      &        ,DTAUX3D=grid%dtaux3d,DTAUY3D=grid%dtauy3d &

diff --git a/dyn_em/start_em.F b/dyn_em/start_em.F
@@ -1243,6 +1243,7 @@ SUBROUTINE start_domain_em ( grid, allowed_to_read &
                       ,grid%xtime1, grid%PBLHAVG, grid%TKEAVG                            &
                       ,grid%ccn_conc                                          & ! RAS
                       ,grid%QKE                                               &!JOE-for mynn
+                      ,grid%pep_pbl                                           & ! EEPS
                       ,grid%landusef,grid%landusef2,grid%mosaic_cat_index                            & ! danli mosaic
                       ,grid%TSK_mosaic,grid%TSLB_mosaic,grid%SMOIS_mosaic,grid%SH2O_mosaic           & ! danli mosaic
                       ,grid%CANWAT_mosaic,grid%SNOW_mosaic,grid%SNOWH_mosaic,grid%SNOWC_mosaic       & ! danli mosaic

diff --git a/phys/Makefile b/phys/Makefile
@@ -48,6 +48,7 @@ MODULES = \
 	module_checkerror.o \
 	module_bl_camuwpbl_driver.o \
 	module_bl_mfshconvpbl.o \
+	module_bl_eepsilon.o \
 	module_shcu_camuwshcu_driver.o \
 	module_shcu_camuwshcu.o \
 	module_shcu_deng.o          \