Skip to content

Add an option to use liquid potential temperature in local mixing of GFSPBL #977

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

Closed
WeiguoWang-NOAA wants to merge 3 commits into from
Closed

Add an option to use liquid potential temperature in local mixing of GFSPBL #977

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
b6bdc1d
Add an option to use liquid potential temperature in temperature equa…
BinLiu-NOAA Jun 2, 2025
8102fe9
Merge remote-tracking branch 'remotes/origin/develop' into feature/gf…
BinLiu-NOAA Jun 2, 2025
b0c32b2
Merge remote-tracking branch 'remotes/origin/develop' into feature/gf…
BinLiu-NOAA Jun 21, 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
4 changes: 2 additions & 2 deletions .gitmodules
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -8,8 +8,8 @@
branch = main
[submodule "ccpp/physics"]
path = ccpp/physics
url = https://github.com/ufs-community/ccpp-physics
branch = ufs/dev
url = https://github.com/hafs-community/ccpp-physics
branch = feature/gfspbl_lpt
[submodule "upp"]
path = upp
url = https://github.com/NOAA-EMC/UPP
Expand Down
6 changes: 5 additions & 1 deletion ccpp/data/GFS_typedefs.F90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -1377,6 +1377,7 @@ module GFS_typedefs
real(kind=kind_phys) :: elmx !< maximum allowed dissipation mixing length in boundary layer mass flux scheme
integer :: sfc_rlm !< choice of near surface mixing length in boundary layer mass flux scheme
integer :: tc_pbl !< control for TC applications in the PBL scheme
integer :: use_lpt !< control for using Liquid Potential Temp for TC applications in the GFSPBL scheme

!--- parameters for canopy heat storage (CHS) parameterization
real(kind=kind_phys) :: h0facu !< CHS factor for sensible heat flux in unstable surface layer
Expand Down Expand Up @@ -3918,6 +3919,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
real(kind=kind_phys) :: elmx = 300. !< maximum allowed dissipation mixing length in boundary layer mass flux scheme
integer :: sfc_rlm = 0 !< choice of near surface mixing length in boundary layer mass flux scheme
integer :: tc_pbl = 0 !< control for TC applications in the PBL scheme
integer :: use_lpt = 0 !< control for using Liquid Potential Temp for TC applications in the GFSPBL scheme
Copy link
Collaborator

@grantfirl grantfirl Jun 13, 2025

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@WeiguoWang-NOAA If this can only ever be 0 or 1, can this be a logical instead?

Copy link
Author

@WeiguoWang-NOAA WeiguoWang-NOAA Jun 13, 2025

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@WeiguoWang-NOAA If this can only ever be 0 or 1, can this be a logical instead?

@grantfirl I would prefer to use an integer. Other different approximations could be developed in the future.


!--- parameters for canopy heat storage (CHS) parameterization
real(kind=kind_phys) :: h0facu = 0.25
Expand Down Expand Up @@ -4184,7 +4186,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
diag_flux, diag_log, &
! vertical diffusion
xkzm_m, xkzm_h, xkzm_s, xkzminv, moninq_fac, dspfac, &
bl_upfr, bl_dnfr, rlmx, elmx, sfc_rlm, tc_pbl, &
bl_upfr, bl_dnfr, rlmx, elmx, sfc_rlm, tc_pbl, use_lpt, &
!--- canopy heat storage parameterization
h0facu, h0facs, &
!--- cellular automata
Expand Down Expand Up @@ -5155,6 +5157,7 @@ subroutine control_initialize (Model, nlunit, fn_nml, me, master, &
Model%elmx = elmx
Model%sfc_rlm = sfc_rlm
Model%tc_pbl = tc_pbl
Model%use_lpt = use_lpt

!--- canopy heat storage parametrization
Model%h0facu = h0facu
Expand Down Expand Up @@ -7000,6 +7003,7 @@ subroutine control_print(Model)
print *, ' elmx : ', Model%elmx
print *, ' sfc_rlm : ', Model%sfc_rlm
print *, ' tc_pbl : ', Model%tc_pbl
print *, ' use_lpt : ', Model%use_lpt
print *, ' '
print *, 'parameters for canopy heat storage parametrization'
print *, ' h0facu : ', Model%h0facu
Expand Down
6 changes: 6 additions & 0 deletions ccpp/data/GFS_typedefs.meta
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -6112,6 +6112,12 @@
units = none
dimensions = ()
type = integer
[use_lpt]
standard_name = control_for_using_LPT_for_TC_applications_in_the_PBL_scheme
long_name = control for using LPT in TC applications in the PBL scheme
units = none
dimensions = ()
type = integer
[h0facu]
standard_name = multiplicative_tuning_parameter_for_reduced_surface_heat_fluxes_due_to_canopy_heat_storage
long_name = canopy heat storage factor for sensible heat flux in unstable surface layer
Expand Down
2 changes: 1 addition & 1 deletion ccpp/physics
Open in desktop
Submodule physics updated 2 files
+13 −1 physics/PBL/SATMEDMF/satmedmfvdifq.F
+7 −0 physics/PBL/SATMEDMF/satmedmfvdifq.meta