Skip to content

Add tuning parameters for the similarity profiles of the surface layer #4

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
wlanghans wants to merge 9 commits into
base: wolfgang/rebase_sofar-071124_on_main
Choose a base branch
from
Open

Add tuning parameters for the similarity profiles of the surface layer #4

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
9 commits
Select commit Hold shift + click to select a range
355f5ff
Added alpha tuning parameters, one for stable one for unstable regime…
May 13, 2025
21a177d
Added parameter tune_ocean_surface_layer to switch on/off the tuning …
May 14, 2025
ccd9ba3
Add missing description of new parameters
May 14, 2025
2fb1409
Added what was missing because of auto-save deactivated
May 14, 2025
34c9d18
Added description of new input arguments to MO subroutine
May 14, 2025
4571c80
Make the PBL schemes (not all of them) consistent with the implemente…
May 14, 2025
87358fa
Fixed the atan functions from python lingo to fortran
May 15, 2025
c78cc6d
Fixed some up
May 16, 2025
8c5be0e
removed extra anpersand
May 16, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
17 changes: 12 additions & 5 deletions GFS_layer/GFS_physics_driver.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1202,7 +1202,9 @@ subroutine GFS_physics_driver &
Sfcprop%ffmm, Sfcprop%ffhh, Sfcprop%uustar, &
wind, Tbd%phy_f2d(1,Model%num_p2d), fm10, fh2, &
sigmaf, vegtype, Sfcprop%shdmax, Model%ivegsrc, &
tsurf, flag_iter) !, Model%redrag, Model%z0s_max, &
tsurf, flag_iter, &
Model%alpha_stable, Model%alpha_unstable, Model%tune_ocean_surface_layer)
! , Model%redrag, Model%z0s_max, &
!Model%do_z0_moon, Model%do_z0_hwrf15, &
!Model%do_z0_hwrf17, Model%do_z0_hwrf17_hwonly, &
!Model%wind_th_hwrf)
Expand All @@ -1225,7 +1227,8 @@ subroutine GFS_physics_driver &
tsurf, flag_iter, Model%redrag, Model%z0s_max, &
Model%do_z0_moon, Model%do_z0_hwrf15, &
Model%do_z0_hwrf17, Model%do_z0_hwrf17_hwonly, &
Model%wind_th_hwrf)
Model%wind_th_hwrf, &
Model%alpha_stable, Model%alpha_unstable, Model%tune_ocean_surface_layer)
else
! GFS original sfc_diff modified by kgao
call sfc_diff (im,Statein%pgr, Statein%ugrs, Statein%vgrs,&
Expand Down Expand Up @@ -1709,7 +1712,8 @@ subroutine GFS_physics_driver &
Model%dspheat, dusfc1, dvsfc1, dtsfc1, dqsfc1, Diag%hpbl,&
gamt, gamq, dkt, kinver, Model%xkzm_m, Model%xkzm_h, &
Model%xkzm_s, lprnt, ipr, &
Model%xkzminv, Model%moninq_fac)
Model%xkzminv, Model%moninq_fac, Model%alpha_stable, &
Model%alpha_unstable, Model%tune_ocean_surface_layer)
! if (lprnt) write(0,*)' dtdtm=',(dtdt(ipr,k),k=1,15)
! if (lprnt) write(0,*)' dqdtm=',(dqdt(ipr,k,1),k=1,15)

Expand All @@ -1731,7 +1735,9 @@ subroutine GFS_physics_driver &
Model%xkzm_ml, Model%xkzm_hl, Model%xkzm_mi, Model%xkzm_hi, &
Model%xkzm_s, Model%xkzminv, Model%do_dk_hb19, &
Model%xkzm_lim, Model%xkgdx, &
Model%rlmn, Model%rlmx, Model%cap_k0_land, dkt)
Model%rlmn, Model%rlmx, Model%cap_k0_land, dkt, &
Model%alpha_stable, Model%alpha_unstable, &
Model%tune_ocean_surface_layer)

elseif (Model%isatmedmf == 1) then
do i=1,im
Expand All @@ -1758,7 +1764,8 @@ subroutine GFS_physics_driver &
Model%xkzm_s, Model%xkzminv, Model%rlmx, Model%zolcru, &
Model%cs0, Model%do_dk_hb19, Model%xkgdx, &
Model%dspfac, Model%bl_upfr, Model%bl_dnfr, dkt, &
flux_cg, flux_en) !cg as up and en as down
flux_cg, flux_en, Model%alpha_stable, Model%alpha_unstable, &
Model%tune_ocean_surface_layer) !cg as up and en as down
endif

elseif (Model%ysupbl) then
Expand Down
15 changes: 14 additions & 1 deletion GFS_layer/GFS_typedefs.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -680,6 +680,9 @@ module GFS_typedefs
logical :: do_z0_hwrf17 !< flag for using z0 scheme in 2017 HWRF (kgao)
logical :: do_z0_hwrf17_hwonly !< flag for using z0 scheme in 2017 HWRF only under high wind (kgao)
real(kind=kind_phys) :: wind_th_hwrf !< wind speed threshold when z0 level off as in HWRF (kgao)
real(kind=kind_phys) :: alpha_stable !< tuning parameter for the dimensionless momentum and scalar gradient function in the surface layer
real(kind=kind_phys) :: alpha_unstable !< tuning parameter for the dimensionless momentum gradient function in the surface layer
logical :: tune_ocean_surface_layer !< if true, use alphas above to modify the surface layer dimensionless gradients over oceans
logical :: hybedmf !< flag for hybrid edmf pbl scheme
logical :: myj_pbl !< flag for NAM MYJ tke scheme
logical :: ysupbl !< flag for ysu pbl scheme (version in WRFV3.8)
Expand Down Expand Up @@ -2401,6 +2404,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
logical :: do_z0_hwrf17 = .false. !< flag for using z0 scheme in 2017 HWRF
logical :: do_z0_hwrf17_hwonly = .false. !< flag for using z0 scheme in 2017 HWRF only under high wind
real(kind=kind_phys) :: wind_th_hwrf = 33. !< wind speed threshold when z0 level off as in HWRF
real(kind=kind_phys) :: alpha_stable = 5. !< tuning parameter for the dimensionless momentum and scalar gradient function in the surface layer
real(kind=kind_phys) :: alpha_unstable = 16. !< tuning parameter for the dimensionless momentum gradient function in the surface layer
logical :: tune_ocean_surface_layer = .false. !< if true, use alphas above to modify the surface layer dimensionless gradients over oceans
logical :: hybedmf = .false. !< flag for hybrid edmf pbl scheme
logical :: myj_pbl = .false. !< flag for NAM MYJ tke-based scheme
logical :: ysupbl = .false. !< flag for hybrid edmf pbl scheme
Expand Down Expand Up @@ -2633,7 +2639,8 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
cscnv, cal_pre, do_aw, do_shoc, shocaftcnv, shoc_cld, &
h2o_phys, pdfcld, shcnvcw, redrag, sfc_gfdl, z0s_max, &
sfc_coupled, do_z0_moon, do_z0_hwrf15, do_z0_hwrf17, &
do_z0_hwrf17_hwonly, wind_th_hwrf, &
do_z0_hwrf17_hwonly, wind_th_hwrf, alpha_stable, &
alpha_unstable, tune_ocean_surface_layer, &
hybedmf, dspheat, lheatstrg, hour_canopy, afac_canopy, &
cnvcld, no_pbl, xkzm_lim, xkzm_fac, xkgdx, &
rlmn, rlmx, zolcru, cs0, &
Expand Down Expand Up @@ -2878,6 +2885,9 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
Model%do_z0_hwrf17 = do_z0_hwrf17
Model%do_z0_hwrf17_hwonly = do_z0_hwrf17_hwonly
Model%wind_th_hwrf = wind_th_hwrf
Model%alpha_stable = alpha_stable
Model%alpha_unstable = alpha_unstable
Model%tune_ocean_surface_layer = tune_ocean_surface_layer
Model%hybedmf = hybedmf
Model%myj_pbl = myj_pbl
Model%ysupbl = ysupbl
Expand Down Expand Up @@ -3585,6 +3595,9 @@ subroutine control_print(Model)
print *, ' do_z0_hwrf17 : ', Model%do_z0_hwrf17
print *, ' do_z0_hwrf17_hwonly : ', Model%do_z0_hwrf17_hwonly
print *, ' wind_th_hwrf : ', Model%wind_th_hwrf
print *, ' alpha_stable : ', Model%alpha_stable
print *, ' alpha_unstable : ', Model%alpha_unstable
print *, ' tune_ocean_surface_layer : ', Model%tune_ocean_surface_layer
print *, ' hybedmf : ', Model%hybedmf
print *, ' myj_pbl : ', Model%myj_pbl
print *, ' ysupbl : ', Model%ysupbl
Expand Down
30 changes: 10 additions & 20 deletions gsmphys/moninedmf.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -93,7 +93,8 @@ subroutine moninedmf(ix,im,km,ntrac,ntcw,dv,du,tau,rtg, &
& prsi,del,prsl,prslk,phii,phil,delt,dspheat, &
& dusfc,dvsfc,dtsfc,dqsfc,hpbl,hgamt,hgamq,dkt, &
& kinver,xkzm_m,xkzm_h,xkzm_s,lprnt,ipr, &
& xkzminv,moninq_fac)
& xkzminv,moninq_fac, alpha_stable, alpha_unstable, &
& tune_ocean_surface_layer)
!
use machine , only : kind_phys
use funcphys , only : fpvs
Expand All @@ -103,11 +104,12 @@ subroutine moninedmf(ix,im,km,ntrac,ntcw,dv,du,tau,rtg, &
!
! arguments
!
logical lprnt
logical lprnt, tune_ocean_surface_layer
integer ipr
integer ix, im, km, ntrac, ntcw, kpbl(im), kinver(im)
!
real(kind=kind_phys) delt, xkzm_m, xkzm_h, xkzm_s
real(kind=kind_phys) delt, xkzm_m, xkzm_h, xkzm_s, alpha_stable,
& alpha_unstable
real(kind=kind_phys) dv(im,km), du(im,km), &
& tau(im,km), rtg(im,km,ntrac), &
& u1(ix,km), v1(ix,km), &
Expand Down Expand Up @@ -512,24 +514,12 @@ subroutine moninedmf(ix,im,km,ntrac,ntcw,dv,du,tau,rtg, &
!! \f[
!! w_s = (u_*^3 + 7\epsilon k w_*^3)^{1/3}
!! \f]
call get_similarity_parameters(im, rbsoil, fm, fh, rimin,
& sfcflg, zfmin, sfcfrac, hpbl, zl(1:im,1),
& aphi16, aphi5, alpha_unstable, alpha_stable,
& tune_ocean_surface_layer, phim, phih, zol)

do i=1,im
zol(i) = max(rbsoil(i)*fm(i)*fm(i)/fh(i),rimin)
if(sfcflg(i)) then
zol(i) = min(zol(i),-zfmin)
else
zol(i) = max(zol(i),zfmin)
endif
zol1 = zol(i)*sfcfrac*hpbl(i)/zl(i,1)
if(sfcflg(i)) then
! phim(i) = (1.-aphi16*zol1)**(-1./4.)
! phih(i) = (1.-aphi16*zol1)**(-1./2.)
tem = 1.0 / (1. - aphi16*zol1)
phih(i) = sqrt(tem)
phim(i) = sqrt(phih(i))
else
phim(i) = 1. + aphi5*zol1
phih(i) = phim(i)
endif
wscale(i) = ustar(i)/phim(i)
ustmin(i) = ustar(i)/aphi5
wscale(i) = max(wscale(i),ustmin(i))
Expand Down
88 changes: 67 additions & 21 deletions gsmphys/satmedmfvdiff.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -36,7 +36,8 @@ subroutine satmedmfvdif(ix,im,km,ntrac,ntcw,ntiw,ntke,
& dspheat,dusfc,dvsfc,dtsfc,dqsfc,hpbl,
& kinver,xkzm_mo,xkzm_ho,xkzm_ml,xkzm_hl,xkzm_mi,xkzm_hi,
& xkzm_s,xkzinv,do_dk_hb19,xkzm_lim,xkgdx,
& rlmn, rlmx, cap_k0_land, dkt_out)
& rlmn, rlmx, cap_k0_land, dkt_out,
& alpha_stable, alpha_unstable, tune_ocean_surface_layer)
!
use machine , only : kind_phys
use funcphys , only : fpvs
Expand All @@ -53,7 +54,8 @@ subroutine satmedmfvdif(ix,im,km,ntrac,ntcw,ntiw,ntke,
!
real(kind=kind_phys) delt, xkzm_s, xkzm_lim,
& xkzm_mo, xkzm_ho, xkzm_ml, xkzm_hl,
& xkzm_mi, xkzm_hi
& xkzm_mi, xkzm_hi,
& alpha_stable, alpha_unstable
real(kind=kind_phys) dv(im,km), du(im,km),
& tdt(im,km), rtg(im,km,ntrac),
& u1(ix,km), v1(ix,km),
Expand All @@ -77,7 +79,7 @@ subroutine satmedmfvdif(ix,im,km,ntrac,ntcw,ntiw,ntke,
! kgao note - q1 and rtg are local var now

!
logical dspheat, cap_k0_land, do_dk_hb19
logical dspheat, cap_k0_land, do_dk_hb19, tune_ocean_surface_layer
! flag for tke dissipative heating
real(kind=kind_phys),dimension(1:im,1:km),intent(OUT)::dkt_out

Expand Down Expand Up @@ -592,24 +594,10 @@ subroutine satmedmfvdif(ix,im,km,ntrac,ntcw,ntiw,ntke,
!
! compute similarity parameters
!
do i=1,im
zol(i) = max(rbsoil(i)*fm(i)*fm(i)/fh(i),rimin)
if(sfcflg(i)) then
zol(i) = min(zol(i),-zfmin)
else
zol(i) = max(zol(i),zfmin)
endif
!
zol1 = zol(i)*sfcfrac*hpbl(i)/zl(i,1)
if(sfcflg(i)) then
tem = 1.0 / (1. - aphi16*zol1)
phih(i) = sqrt(tem)
phim(i) = sqrt(phih(i))
else
phim(i) = 1. + aphi5*zol1
phih(i) = phim(i)
endif
enddo
call get_similarity_parameters(im, rbsoil, fm, fh, rimin,
& sfcflg, zfmin, sfcfrac, hpbl, zl(1:im,1),
& aphi16, aphi5, alpha_unstable, alpha_stable,
& tune_ocean_surface_layer, phim, phih, zol)
!
do i=1,im
if(pblflg(i)) then
Expand Down Expand Up @@ -1564,3 +1552,61 @@ subroutine tridit(l,n,nt,cl,cm,cu,rt,au,at)
c-----------------------------------------------------------------------
return
end

subroutine get_similarity_parameters(
& im, rbsoil, fm, fh, rimin,
& sfcflg, zfmin, sfcfrac, hpbl, zl,
& aphi16, aphi5, alpha_unstable, alpha_stable,
& tune_ocean_surface_layer, phim, phih, zol
& )

use machine, only : kind_phys

implicit none
! input
integer, intent(in) :: im
real(kind=kind_phys), dimension(1:im), intent(in) ::
& rbsoil, fm, fh, hpbl, zl
real(kind=kind_phys), intent(in) :: rimin, zfmin, aphi16, aphi5,
& alpha_unstable, alpha_stable, sfcfrac
logical, dimension(1:im), intent(in) :: sfcflg
logical, intent(in) :: tune_ocean_surface_layer

! output
real(kind=kind_phys), dimension(1:im), intent(out) :: phim, phih,
& zol

! local
integer i
real(kind=kind_phys) :: zol1, tem

do i=1,im
zol(i) = max(rbsoil(i)*fm(i)*fm(i)/fh(i),rimin)
if(sfcflg(i)) then
zol(i) = min(zol(i),-zfmin)
else
zol(i) = max(zol(i),zfmin)
endif
!
zol1 = zol(i)*sfcfrac*hpbl(i)/zl(i)
if(sfcflg(i)) then
tem = 1.0 / (1. - aphi16*zol1)
phih(i) = sqrt(tem)
if (tune_ocean_surface_layer) then
tem = 1.0 / (1. - alpha_unstable*zol1)
phim(i) = sqrt(sqrt(tem))
else
phim(i) = sqrt(phih(i))
end if
else
if (tune_ocean_surface_layer) then
phim(i) = 1. + alpha_stable*zol1
else
phim(i) = 1. + aphi5*zol1
end if
phih(i) = phim(i)
endif
enddo
return
end subroutine

30 changes: 9 additions & 21 deletions gsmphys/satmedmfvdifq.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -56,7 +56,8 @@ subroutine satmedmfvdifq(ix,im,km,ntrac,ntcw,ntiw,ntke,
& kinver,xkzm_mo,xkzm_ho,xkzm_ml,xkzm_hl,xkzm_mi,xkzm_hi,
& xkzm_s,xkzinv,rlmx,zolcru,cs0,
& do_dk_hb19,xkgdx,dspfac,bl_upfr,bl_dnfr,dkt_out,
& flux_up, flux_dn)
& flux_up, flux_dn,
& alpha_stable, alpha_unstable, tune_ocean_surface_layer)
!
use machine , only : kind_phys
use funcphys , only : fpvs
Expand All @@ -73,7 +74,8 @@ subroutine satmedmfvdifq(ix,im,km,ntrac,ntcw,ntiw,ntke,
!
real(kind=kind_phys) delt, xkzm_mo, xkzm_ho, xkzm_s, dspfac,
& bl_upfr, bl_dnfr, xkzm_ml, xkzm_hl,
& xkzm_mi, xkzm_hi
& xkzm_mi, xkzm_hi, alpha_stable,
& alpha_unstable
real(kind=kind_phys) dv(im,km), du(im,km),
& tdt(im,km), rtg(im,km,ntrac),
& u1(ix,km), v1(ix,km),
Expand All @@ -97,7 +99,7 @@ subroutine satmedmfvdifq(ix,im,km,ntrac,ntcw,ntiw,ntke,
& rtg_in(im,km,ntrac)
! kgao note - q1 and rtg are local var now
!
logical dspheat, do_dk_hb19
logical dspheat, do_dk_hb19, tune_ocean_surface_layer
! flag for tke dissipative heating
real(kind=kind_phys)::dkt_out(im,km),flux_up(im,km),flux_dn(im,km)
!
Expand Down Expand Up @@ -619,24 +621,10 @@ subroutine satmedmfvdifq(ix,im,km,ntrac,ntcw,ntiw,ntke,
!
! compute similarity parameters
!
do i=1,im
zol(i) = max(rbsoil(i)*fm(i)*fm(i)/fh(i),rimin)
if(sfcflg(i)) then
zol(i) = min(zol(i),-zfmin)
else
zol(i) = max(zol(i),zfmin)
endif
!
zol1 = zol(i)*sfcfrac*hpbl(i)/zl(i,1)
if(sfcflg(i)) then
tem = 1.0 / (1. - aphi16*zol1)
phih(i) = sqrt(tem)
phim(i) = sqrt(phih(i))
else
phim(i) = 1. + aphi5*zol1
phih(i) = phim(i)
endif
enddo
call get_similarity_parameters(im, rbsoil, fm, fh, rimin,
& sfcflg, zfmin, sfcfrac, hpbl, zl(1:im,1),
& aphi16, aphi5, alpha_unstable, alpha_stable,
& tune_ocean_surface_layer, phim, phih, zol)
!
do i=1,im
if(pblflg(i)) then
Expand Down
18 changes: 13 additions & 5 deletions gsmphys/sfc_diff_coupled.f
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -4,7 +4,9 @@ subroutine sfc_diff_coupled(im,ps,u1,v1,t1,q1,z1,
& stress,fm,fh,
& ustar,wind,ddvel,fm10,fh2,
& sigmaf,vegtype,shdmax,ivegsrc,
& tsurf,flag_iter) !,redrag,
& tsurf,flag_iter, alpha_stable, alpha_unstable,
& tune_ocean_surface_layer)
! ,redrag,
! & z0s_max)
! & do_z0_moon, do_z0_hwrf15, do_z0_hwrf17,
! & do_z0_hwrf17_hwonly, wind_th_hwrf)
Expand All @@ -26,10 +28,12 @@ subroutine sfc_diff_coupled(im,ps,u1,v1,t1,q1,z1,
&, prsl1, prslki, stress
&, fm, fh, ustar, wind, ddvel
&, fm10, fh2, sigmaf, shdmax
&, tsurf, snwdph
&, tsurf, snwdph, alpha_stable
&, alpha_unstable

integer, dimension(im) ::vegtype, islimsk

logical flag_iter(im)
logical flag_iter(im), tune_ocean_surface_layer
! logical redrag
! logical do_z0_moon, do_z0_hwrf15, do_z0_hwrf17 ! kgao
! &, do_z0_hwrf17_hwonly ! kgao
Expand Down Expand Up @@ -143,7 +147,9 @@ subroutine sfc_diff_coupled(im,ps,u1,v1,t1,q1,z1,
! --- call similarity

call monin_obukhov_similarity
& (z1(i), snwdph(i), thv1, wind(i), z0max, ztmax, tvs,
& (islimsk(i), z1(i), snwdph(i), thv1, wind(i), z0max,
Comment thread
wlanghans marked this conversation as resolved.
& ztmax, tvs,
& alpha_stable, alpha_unstable, tune_ocean_surface_layer,
& rb(i), fm(i), fh(i), fm10(i), fh2(i),
& cm(i), ch(i), stress(i), ustar(i))

Expand All @@ -165,7 +171,9 @@ subroutine sfc_diff_coupled(im,ps,u1,v1,t1,q1,z1,
! --- call similarity
! kgao: use z0/zt to do sfc diagnosis
call monin_obukhov_similarity
& (z1(i), snwdph(i), thv1, wind(i), z0max, ztmax, tvs,
& (islimsk(i), z1(i), snwdph(i), thv1, wind(i), z0max,
& ztmax, tvs,
& alpha_stable, alpha_unstable, tune_ocean_surface_layer,
& rb(i), fm(i), fh(i), fm10(i), fh2(i),
& cm(i), ch(i), tem1, tem2) !stress(i), ustar(i))

Expand Down
Loading