Skip to content
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

Soil-veg-atm coupling #885

Closed
wants to merge 24 commits into from
Closed

Soil-veg-atm coupling #885

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
24 commits
Select commit Hold shift + click to select a range
4cbea92
set up option for thermal roughness
Jan 30, 2022
22de66b
change czil
Jan 31, 2022
10fa17e
canopy height dependant czil
Feb 16, 2022
b90d4e2
add canopy heat storage and gvf impact on thermal conductivity
HelinWei-NOAA Mar 7, 2022
ae7ac42
add sfcdif3 as a separate subroutine
barlage Mar 7, 2022
c50f50a
change fveg to shdfac in sfcdif3 vege call
barlage Mar 7, 2022
53c0c7a
move trs options to sfcdif3
barlage Mar 7, 2022
f093f77
fix missing czil1 in vege_flux
barlage Mar 7, 2022
13a1d14
add some clean up to energy
barlage Mar 8, 2022
ebb4fa1
add some groundwater mods from ncar code
barlage Mar 8, 2022
b6e0573
Merge pull request #17 from barlage/sfcdif3_subroutine
HelinWei-NOAA Mar 8, 2022
41cf4ec
gvf impact on thermal conductivity limited to the first soil layer
HelinWei-NOAA Mar 8, 2022
b43c030
Merge pull request #18 from barlage/sfcdif3_subroutine
HelinWei-NOAA Mar 8, 2022
c1d813e
correct the reference height
HelinWei-NOAA Mar 8, 2022
2ae3f48
Merge branch 'NCAR:main' into lsm_upgrades_for_p8c
HelinWei-NOAA Mar 9, 2022
11b50ca
to read new hig-res ice climatology data
HelinWei-NOAA Mar 9, 2022
4ed3982
replace fveg by lai/laimax to be used for dependent
HelinWei-NOAA Mar 10, 2022
8e1b316
simplify the code with internal function maxval
HelinWei-NOAA Mar 10, 2022
70507a0
to avoid exception floating point
HelinWei-NOAA Mar 10, 2022
d33598b
revert the df1 change due to some negative impact on surface temperature
HelinWei-NOAA Mar 11, 2022
3095d71
correct the condition to avoid a divide by zero exception
HelinWei-NOAA Mar 14, 2022
27ea849
further refinement of the impact of vegetation on zvfun
HelinWei-NOAA Mar 15, 2022
c722905
replace shdfac by fveg for zvfun
HelinWei-NOAA Mar 16, 2022
4284846
modify the eddy diffusivity for heat at the top of the canopy
wzzheng90 Mar 18, 2022
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
24 changes: 19 additions & 5 deletions physics/module_sf_noahmp_glacier.f90
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -62,6 +62,7 @@ module noahmp_glacier_globals
INTEGER :: OPT_GLA != 1 !(suggested 1)

INTEGER :: OPT_SFC != 1 !(suggested 1)
INTEGER :: OPT_TRS != 1 !(suggested 2)

! adjustable parameters for snow processes

Expand Down Expand Up @@ -1129,8 +1130,10 @@ subroutine glacier_flux (nsoil ,nsnow ,emg ,isnow ,df ,dzsnso
real (kind=kind_phys) :: b !< temporary calculation
real (kind=kind_phys) :: t, tdc !< kelvin to degree celsius with limit -50 to +50
real (kind=kind_phys), dimension( 1:nsoil) :: sice !< soil ice
real (kind=kind_phys) :: czil !< calculate roughness length of heat

tdc(t) = min( 50., max(-50.,(t-tfrz)) )
czil=0.1

! -----------------------------------------------------------------
! initialization variables that do not depend on stability iteration
Expand All @@ -1155,10 +1158,18 @@ subroutine glacier_flux (nsoil ,nsnow ,emg ,isnow ,df ,dzsnso
fv = ur*vkc/log(zlvli/z0m)
reyni = fv*z0m/(1.5e-05) !introduction of fv dependent z0h for the iter

if (reyni .gt. 2.0) then
z0h = z0m/exp(2.46*(reyni)**0.25 - log(7.4)) !Brutsaert 1982
else
z0h = z0m/exp(-log(0.397)) !Brusaert 1982, table 4
if (opt_trs == 1) then
z0h = z0m
elseif (opt_trs == 2) then
z0h = z0m*exp(-czil*0.4*258.2*sqrt(fv*z0m))
elseif (opt_trs == 3) then
z0h = z0m*0.1
elseif (opt_trs == 4) then
if (reyni .gt. 2.0) then
z0h = z0m/exp(2.46*(reyni)**0.25 - log(7.4)) !Brutsaert 1982
else
z0h = z0m/exp(-log(0.397)) !Brusaert 1982, table 4
endif
endif

z0h_total = z0h
Expand Down Expand Up @@ -3328,7 +3339,8 @@ end subroutine error_glacier
! ==================================================================================================

!>\ingroup NoahMP_LSM
subroutine noahmp_options_glacier(iopt_alb ,iopt_snf ,iopt_tbot, iopt_stc, iopt_gla, iopt_sfc)
subroutine noahmp_options_glacier(iopt_alb ,iopt_snf ,iopt_tbot, iopt_stc, iopt_gla,&
iopt_sfc, iopt_trs)

implicit none

Expand All @@ -3339,6 +3351,7 @@ subroutine noahmp_options_glacier(iopt_alb ,iopt_snf ,iopt_tbot, iopt_stc, iop
!! 1 -> semi-implicit; 2 -> full implicit (original noah)
integer, intent(in) :: iopt_gla !< glacier option (1->phase change; 2->simple)
integer, intent(in) :: iopt_sfc !< sfc scheme option
integer, intent(in) :: iopt_trs !< thermal roughness option

! -------------------------------------------------------------------------------------------------

Expand All @@ -3348,6 +3361,7 @@ subroutine noahmp_options_glacier(iopt_alb ,iopt_snf ,iopt_tbot, iopt_stc, iop
opt_stc = iopt_stc
opt_gla = iopt_gla
opt_sfc = iopt_sfc
opt_trs = iopt_trs

end subroutine noahmp_options_glacier

Expand Down
Loading