Skip to content

Commit

Permalink
Add scale-aware to the new Tiedtke cu scheme as an option (#1806)
Browse files Browse the repository at this point in the history 
Add scale-aware to the new Tiedtke scheme as an option

TYPE: new feature

KEYWORDS: new Tiedtke cu scheme, scale aware option

SOURCE: internal

DESCRIPTION OF CHANGES:
A scale-aware or grid-distance-dependent option is added to the new Tiedtke scheme. The code is based on the paper by Wang (2022), and can be turned on by using namelist option ntiedtke_dx_opt = 1. The option becomes active when grid sizes is less than 15 km.
Wang, W., 2022: Forecasting Convection with a “Scale-Aware” Tiedtke Cumulus Parameterization Scheme at Kilometer Scales. _Wea. Forecasting_, 37/8, 1491 - 1507, 10.1175/WAF-D-21-0179.1.

LIST OF MODIFIED FILES: 
M       Registry/Registry.EM_COMMON
M       dyn_em/module_first_rk_step_part1.F
M       phys/module_cu_ntiedtke.F
M       phys/module_cumulus_driver.F
M       run/README.namelist

TESTS CONDUCTED: 
1. Do mods fix problem? How can that be demonstrated, and was that test conducted?
2. The Jenkins tests have passed.

RELEASE NOTE: A scale-aware or grid-distance-dependent option is added to the new Tiedtke scheme. The code is based on the paper by Wang (2022), and can be turned on by using namelist option ntiedtke_dx_opt = 1. The option becomes active when grid sizes is less than 15 km.
  • Loading branch information
@weiwangncar
weiwangncar authored Jan 24, 2023
1 parent 2a29582 commit 7639166
Show file tree
Hide file tree
Showing 5 changed files with 85 additions and 10 deletions.
1 change: 1 addition & 0 deletions Registry/Registry.EM_COMMON
View file
Original file line number Diff line number Diff line change
Expand Up @@ -2449,6 +2449,7 @@ rconfig integer cu_diag namelist,physics max_domains 0
rconfig integer kf_edrates namelist,physics max_domains 0 rh "kf_edrates" "output entrainment/detrainment rates and convective timescale for KF schemes" ""
rconfig integer kfeta_trigger namelist,physics 1 1 rh "KFETA Trigger function" "" ""
rconfig integer nsas_dx_factor namelist,physics 1 0 rh "NSAS DX-dependent option" "" ""
rconfig integer ntiedtke_dx_opt namelist,physics 1 0 rh "nTiedtke DX-dependent option" "" ""
rconfig real CUDT namelist,physics max_domains 0 h "CUDT" "" ""
rconfig real GSMDT namelist,physics max_domains 0 h "GSMDT" "" ""
rconfig integer ISFFLX namelist,physics 1 1 irh "ISFFLX" "" ""
Expand Down
1 change: 1 addition & 0 deletions dyn_em/module_first_rk_step_part1.F
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1474,6 +1474,7 @@ SUBROUTINE first_rk_step_part1 ( grid , config_flags &
!BSINGH -ENDS
& ,KFETA_TRIGGER=config_flags%kfeta_trigger &
& ,NSAS_DX_FACTOR=config_flags%nsas_dx_factor &
& ,NTIEDTKE_DX_OPT=config_flags%ntiedtke_dx_opt &
! Dimension arguments
& ,IDS=ids,IDE=ide, JDS=jds,JDE=jde, KDS=kds,KDE=kde &
& ,IMS=ims,IME=ime, JMS=jms,JME=jme, KMS=kms,KME=kme &
Expand Down
86 changes: 76 additions & 10 deletions phys/module_cu_ntiedtke.F
View file
Original file line number Diff line number Diff line change
Expand Up @@ -150,6 +150,7 @@ subroutine cu_ntiedtke( &
,u3d,v3d,w,t3d,qv3d,qc3d,qi3d,pi3d,rho3d &
,qvften,thften &
,dz8w,pcps,p8w,xland,cu_act_flag,dx &
,ntiedtke_dx_opt &
,ids,ide, jds,jde, kds,kde &
,ims,ime, jms,jme, kms,kme &
,its,ite, jts,jte, kts,kte &
Expand Down Expand Up @@ -191,6 +192,7 @@ subroutine cu_ntiedtke( &
!-- dz8w dz between full levels (m)
!-- qfx upward moisture flux at the surface (kg/m^2/s)
!-- hfx upward heat flux at the surface (w/m^2)
!-- ntiedtke_dx_opt whether the schemem is scale-aware
!-- dt time step (s)
!-- ids start index for i in domain
!-- ide end index for i in domain
Expand All @@ -217,6 +219,8 @@ subroutine cu_ntiedtke( &
itimestep, &
stepcu

integer, intent(in) :: ntiedtke_dx_opt

real, intent(in) :: &
dt
real, dimension(ims:ime, jms:jme), intent(in) :: &
Expand Down Expand Up @@ -426,7 +430,7 @@ subroutine cu_ntiedtke( &
!
!########################################################################
call tiecnvn(u1,v1,t1,q1,q2,q3,q1b,t1b,ghtl,ghti,omg,prsl,prsi,evap,heatflux, &
rn,slimsk,im,kx,kx1,delt,dx2d)
rn,slimsk,im,kx,kx1,delt,dx2d,ntiedtke_dx_opt)

do i=its,ite
raincv(i,j)=rn(i)/stepcu
Expand Down Expand Up @@ -555,7 +559,7 @@ end subroutine ntiedtkeinit
! level 1 subroutine 'tiecnvn'
!-----------------------------------------------------------------
subroutine tiecnvn(pu,pv,pt,pqv,pqc,pqi,pqvf,ptf,poz,pzz,pomg, &
& pap,paph,evap,hfx,zprecc,lndj,lq,km,km1,dt,dx)
& pap,paph,evap,hfx,zprecc,lndj,lq,km,km1,dt,dx,ntiedtke_dx_opt)
!-----------------------------------------------------------------
! this is the interface between the model and the mass
! flux convection module
Expand All @@ -577,6 +581,7 @@ subroutine tiecnvn(pu,pv,pt,pqv,pqc,pqi,pqvf,ptf,poz,pzz,pomg, &
& zqsat(lq,km), pqc(lq,km), pqi(lq,km), zrain(lq)
real pqvf(lq,km), ptf(lq,km)
real dx(lq)
integer ntiedtke_dx_opt

integer icbot(lq), ictop(lq), ktype(lq), lndj(lq)
logical locum(lq)
Expand All @@ -585,6 +590,18 @@ subroutine tiecnvn(pu,pv,pt,pqv,pqc,pqi,pqvf,ptf,poz,pzz,pomg, &
integer i,j,k,lq,km,km1
real dt,ztpp1
real zew,zqs,zcor
real scale_fac(lq), dxref
!
! set scale-dependency factor when dx is < 15 km
!
dxref = 15000.
do j=1,lq
if (dx(j).lt.dxref .and. ntiedtke_dx_opt.eq.1) then
scale_fac(j) = (1.+log(dxref/dx(j)))**3
else
scale_fac(j) = 1.+1.33e-5*dx(j)
end if
end do
!
ztmst=dt
!
Expand Down Expand Up @@ -638,7 +655,8 @@ subroutine tiecnvn(pu,pv,pt,pqv,pqc,pqi,pqvf,ptf,poz,pzz,pomg, &
& pssfc, locum, &
& ktype, icbot, ictop, ztu, zqu, &
& zlu, zlude, zmfu, zmfd, zrain,&
& pcte, phhfl, lndj, pgeoh, dx)
& pcte, phhfl, lndj, pgeoh, dx, &
& scale_fac, ntiedtke_dx_opt )
!
! to include the cloud water and cloud ice detrained from convection
!
Expand Down Expand Up @@ -693,7 +711,8 @@ subroutine cumastrn &
& pssfc, ldcum, &
& ktype, kcbot, kctop, ptu, pqu,&
& plu, plude, pmfu, pmfd, prain,&
& pcte, phhfl, lndj, zgeoh, dx)
& pcte, phhfl, lndj, zgeoh, dx, &
& scale_fac, ntiedtke_dx_opt )
implicit none
!
!***cumastrn* master routine for cumulus massflux-scheme
Expand Down Expand Up @@ -801,6 +820,8 @@ subroutine cumastrn &
& ktype(klon), lndj(klon)
logical ldcum(klon)
logical loddraf(klon), llo1, llo2(klon)
real scale_fac(klon)
integer ntiedtke_dx_opt

! local varaiables
real zcons,zcons2,zqumqe,zdqmin,zdh,zmfmax
Expand Down Expand Up @@ -901,7 +922,8 @@ subroutine cumastrn &
& zuu, zvu, pmfu, zmfub,&
& zmfus, zmfuq, zmful, plude, zdmfup,&
& kcbot, kctop, ictop0, icum, ztmst,&
& zqsenh, zlglac, lndj, wup, wbase, kdpl, pmfude_rate )
& zqsenh, zlglac, lndj, wup, wbase, kdpl, pmfude_rate,&
& scale_fac, ntiedtke_dx_opt )

!* (b) check cloud depth and change entrainment rate accordingly
! calculate precipitation rate (for downdraft calculation)
Expand Down Expand Up @@ -1017,10 +1039,14 @@ subroutine cumastrn &
if(ldcum(jl).and.ktype(jl).eq.1) then
ikb = kcbot(jl)
ikt = kctop(jl)
ztau = ztauc(jl) * (1.+1.33e-5*dx(jl))
ztau = ztauc(jl) * scale_fac(jl)
ztau = max(ztmst,ztau)
ztau = max(360.,ztau)
ztau = min(10800.,ztau)
if(ntiedtke_dx_opt.eq.1) then
ztau = min(43200.,ztau)
else
ztau = min(10800.,ztau)
end if
if(nonequil) then
zcape2(jl)= max(0.,zcape2(jl))
zcape(jl) = max(0.,min(zcape1(jl)-zcape2(jl),5000.))
Expand Down Expand Up @@ -1060,6 +1086,9 @@ subroutine cumastrn &
zmfub1(jl) = zmfub(jl)
end if
zmfub1(jl) = min(zmfub1(jl),zmfmax)
if(ntiedtke_dx_opt.eq.1) then
zmfub1(jl) = zmfub1(jl)/(scale_fac(jl))**(.5)
end if
end if

!* 6.3 mid-level convection - nothing special
Expand Down Expand Up @@ -2102,7 +2131,8 @@ subroutine cuascn &
& puu, pvu, pmfu, pmfub, &
& pmfus, pmfuq, pmful, plude, pdmfup,&
& kcbot, kctop, kctop0, kcum, ztmst,&
& pqsenh, plglac, lndj, wup, wbase, kdpl, pmfude_rate)
& pqsenh, plglac, lndj, wup, wbase, kdpl, pmfude_rate,&
& scale_fac, ntiedtke_dx_opt )
implicit none
! this routine does the calculations for cloud ascents
! for cumulus parameterization
Expand Down Expand Up @@ -2223,12 +2253,18 @@ subroutine cuascn &
real zrnew,zz,zdmfeu,zdmfdu,dp
real zfac,zbuoc,zdkbuo,zdken,zvv,zarg,zchange,zxe,zxs,zdshrd
real atop1,atop2,abot
real scale_fac(klon)
integer ntiedtke_dx_opt
!--------------------------------
!* 1. specify parameters
!--------------------------------
zcons2=3./(g*ztmst)
zfacbuo = 0.5/(1.+0.5)
if(ntiedtke_dx_opt.eq.1) then
zprcdgw = cprcon*zrg/(scale_fac(jl)**(1./3.))
else
zprcdgw = cprcon*zrg
end if
z_cldmax = 5.e-3
z_cwifrac = 0.5
z_cprc2 = 0.5
Expand Down Expand Up @@ -2307,7 +2343,7 @@ subroutine cuascn &
& pgeo, pgeoh, ldcum, ktype, klab, zlrain,&
& pmfu, pmfub, kcbot, ptu,&
& pqu, plu, puu, pvu, pmfus,&
& pmfuq, pmful, pdmfup)
& pmfuq, pmful, pdmfup, ntiedtke_dx_opt )
is = 0
jlm = 0
do jl = 1,klon
Expand Down Expand Up @@ -3709,7 +3745,7 @@ subroutine cubasmcn &
& pgeo, pgeoh, ldcum, ktype, klab, plrain,&
& pmfu, pmfub, kcbot, ptu,&
& pqu, plu, puu, pvu, pmfus,&
& pmfuq, pmful, pdmfup)
& pmfuq, pmful, pdmfup, ntiedtke_dx_opt )
implicit none
! m.tiedtke e.c.m.w.f. 12/89
! c.zhang iprc 05/2012
Expand Down Expand Up @@ -3746,9 +3782,38 @@ subroutine cubasmcn &
! local variabels
integer jl,kk,klev,klon,klevp1,klevm1
real zzzmb
integer ntiedtke_dx_opt
!--------------------------------------------------------
!* 1. calculate entrainment and detrainment rates
! -------------------------------------------------------
if(ntiedtke_dx_opt.eq.1) then
do jl=1,klon
if(.not.ldcum(jl) .and. klab(jl,kk+1).eq.0) then
if(lmfmid .and. pqen(jl,kk) .gt. 0.80*pqsen(jl,kk).and. &
! ww (no mid level conv, if the air is saturated)
pqen(jl,kk) .le. pqsen(jl,kk).and. &
pgeo(jl,kk)*zrg .gt. 5.0e2 .and. &
& pgeo(jl,kk)*zrg .lt. 1.0e4 ) then
ptu(jl,kk+1)=(cpd*pten(jl,kk)+pgeo(jl,kk)-pgeoh(jl,kk+1))&
& *rcpd
pqu(jl,kk+1)=pqen(jl,kk)
plu(jl,kk+1)=0.
zzzmb=max(cmfcmin,-pverv(jl,kk)*zrg)
zzzmb=min(zzzmb,cmfcmax)
pmfub(jl)=zzzmb
pmfu(jl,kk+1)=pmfub(jl)
pmfus(jl,kk+1)=pmfub(jl)*(cpd*ptu(jl,kk+1)+pgeoh(jl,kk+1))
pmfuq(jl,kk+1)=pmfub(jl)*pqu(jl,kk+1)
pmful(jl,kk+1)=0.
pdmfup(jl,kk+1)=0.
kcbot(jl)=kk
klab(jl,kk+1)=1
plrain(jl,kk+1)=0.0
ktype(jl)=3
end if
end if
end do
else
do jl=1,klon
if(.not.ldcum(jl) .and. klab(jl,kk+1).eq.0) then
if(lmfmid .and. pqen(jl,kk) .gt. 0.80*pqsen(jl,kk).and. &
Expand All @@ -3773,6 +3838,7 @@ subroutine cubasmcn &
end if
end if
end do
end if
return
end subroutine cubasmcn
!---------------------------------------------------------
Expand Down
4 changes: 4 additions & 0 deletions phys/module_cumulus_driver.F
View file
Original file line number Diff line number Diff line change
Expand Up @@ -125,6 +125,7 @@ SUBROUTINE cumulus_driver(grid &
! Optional trigger function activation variable
,kfeta_trigger &
,nsas_dx_factor &
,ntiedtke_dx_opt &
#if ( WRF_DFI_RADAR == 1 )
! Optional CAP suppress option --- 3.2 CLEANUP TODO -- THESE SHOULD BE OPTIONAL, NOT #IF/#ENDIF
,do_capsuppress &
Expand Down Expand Up @@ -324,6 +325,7 @@ SUBROUTINE cumulus_driver(grid &
!-- W0AVG average vertical velocity, (for KF scheme) (m/s)
!-- kfeta_trigger namelist for KF trigger (=1, default; =2, moisture-advection-dependent trigger)
!-- nsas_dx_factor namelist for NSAS deep scheme to have some dependency on grid sizes
!-- ntiedtke_dx_opt namelist for nTiedtke cu scheme to have some dependency on grid sizes
!-- rho density (kg/m^3)
!-- CONVCLD Convective cloud (for BMJ scheme) (kg/m^2)
!-- CCLDFRA convective cloud fraction (for BMJ scheme)
Expand Down Expand Up @@ -525,6 +527,7 @@ SUBROUTINE cumulus_driver(grid &
INTEGER, INTENT(IN ), OPTIONAL :: kfeta_trigger
INTEGER, INTENT(IN ), OPTIONAL :: nsas_dx_factor
INTEGER, INTENT(IN ) :: ntiedtke_dx_opt
REAL, INTENT(IN ) :: DT, DX
REAL, DIMENSION( ims:ime , jms:jme ), OPTIONAL, &
Expand Down Expand Up @@ -1423,6 +1426,7 @@ SUBROUTINE cumulus_driver(grid &
,QV3D=QV_CURR,QC3D=QC_CURR,QI3D=QI_CURR &
,DZ8W=dz8w,PCPS=p,P8W=p8w,XLAND=XLAND,DX=dx2d &
,CU_ACT_FLAG=CU_ACT_FLAG &
,NTIEDTKE_DX_OPT=NTIEDTKE_DX_OPT &
,IDS=ids,IDE=ide,JDS=jds,JDE=jde,KDS=kds,KDE=kde &
,IMS=ims,IME=ime,JMS=jms,JME=jme,KMS=kms,KME=kme &
,ITS=its,ITE=ite,JTS=jts,JTE=jte,KTS=kts,KTE=kte &
Expand Down
3 changes: 3 additions & 0 deletions run/README.namelist
View file
Original file line number Diff line number Diff line change
Expand Up @@ -852,6 +852,9 @@ Namelist variables for controlling the adaptive time step option:
= 3, for small grid distances (DX < 5 km)
nsas_dx_factor = 0, ! default option
= 1, NSAS grid-distance dependent option
ntiedtke_dx_opt = 0, default option
= 1, new Tiedtke grid-distance dependent option (scale-aware)
when grid size falls below 15 km
For KF-CuP scheme: recommended to use with cu_rad_feedback
shallowcu_forced_ra(max_dom) radiative impact of shallow Cu by a prescribed maximum cloud fraction
= .false., option off, default
Expand Down

0 comments on commit 7639166

Please sign in to comment.