Skip to content

Commit

Permalink
Merge pull request NOAA-GFDL#48 from spencerkclark/noahmp/update-snow…
Browse files Browse the repository at this point in the history 
…-fraction-parameters

Update parameters used for snow cover fraction scheme
  • Loading branch information
@lharris4
lharris4 authored Jun 18, 2024
2 parents 573746a + e60641b commit fb99cae
Show file tree
Hide file tree
Showing 3 changed files with 16 additions and 5 deletions.
3 changes: 2 additions & 1 deletion gsmphys/module_sf_noahmplsm.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -172,6 +172,7 @@ module module_sf_noahmplsm
real (kind=kind_phys) :: den !tree density (no. of trunks per m2)
real (kind=kind_phys) :: rc !tree crown radius (m)
real (kind=kind_phys) :: mfsno !snowmelt m parameter ()
real (kind=kind_phys) :: scffac !snow cover factor (m)
real (kind=kind_phys) :: saim(12) !monthly stem area index, one-sided
real (kind=kind_phys) :: laim(12) !monthly leaf area index, one-sided
real (kind=kind_phys) :: sla !single-side leaf area per kg [m2/kg]
Expand Down Expand Up @@ -1657,7 +1658,7 @@ subroutine energy (parameters,ice ,vegtyp ,ist ,nsnow ,nsoil , & !in
if(snowh.gt.0.) then
bdsno = sneqv / snowh
fmelt = (bdsno/100.)**parameters%mfsno
fsno = tanh( snowh /(2.5* z0 * fmelt))
fsno = tanh( snowh /(parameters%scffac * fmelt))
! print*,'bdsno=',bdsno,sneqv,snowh,parameters%mfsno,fmelt,fsno
endif

Expand Down
17 changes: 13 additions & 4 deletions gsmphys/noahmp_tables.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -88,11 +88,20 @@ module noahmp_tables
& 0.30, 0.30, 0.00, 0.00, 0.00, 0.00, &
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 /

! C. He 12/17/2020: optimized MFSNO values dependent on land type based on evaluation with SNOTEL SWE and MODIS SCF, surface albedo
real (kind=kind_phys) :: mfsno_table(mvt) !snowmelt curve parameter ()
data ( mfsno_table(i),i=1,mvt) / 2.50, 2.50, 2.50, 2.50, 2.50, 2.50, &
& 2.50, 2.50, 2.50, 2.50, 2.50, 2.50, &
& 2.50, 2.50, 2.50, 2.50, 2.50, 2.50, &
& 2.50, 2.50, 0.00, 0.00, 0.00, 0.00, &
data ( mfsno_table(i),i=1,mvt) / 1.00, 1.00, 1.00, 1.00, 1.00, 2.00, &
& 2.00, 2.00, 2.00, 2.00, 3.00, 3.00, &
& 4.00, 4.00, 2.50, 3.00, 3.00, 3.50, &
& 3.50, 3.50, 0.00, 0.00, 0.00, 0.00, &
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 /

! C. He 12/17/2020: optimized snow cover factor (m) in SCF formulation to replace original constant 2.5*z0,z0=0.002m, based on evaluation with SNOTEL SWE and MODIS SCF, surface albedo
real (kind=kind_phys) :: scffac_table(mvt) !snow cover factor (m)
data ( scffac_table(i),i=1,mvt) / 0.005, 0.005, 0.005, 0.005, 0.005, 0.008, &
& 0.008, 0.010, 0.010, 0.010, 0.010, 0.007, &
& 0.021, 0.013, 0.015, 0.008, 0.015, 0.015, &
& 0.015, 0.015, 0.00, 0.00, 0.00, 0.00, &
& 0.00, 0.00, 0.00, 0.00, 0.00, 0.00 /

!
Expand Down
1 change: 1 addition & 0 deletions gsmphys/sfc_noahmp_drv.f
View file
Original file line number Diff line number Diff line change
Expand Up @@ -963,6 +963,7 @@ subroutine transfer_mp_parameters (vegtype,soiltype,slopetype, &
parameters%den = den_table(vegtype) !tree density (no. of trunks per m2)
parameters%rc = rc_table(vegtype) !tree crown radius (m)
parameters%mfsno = mfsno_table(vegtype) !snowmelt m parameter ()
parameters%scffac = scffac_table(vegtype) !snow cover factor (m)
parameters%saim = saim_table(vegtype,:) !monthly stem area index, one-sided
parameters%laim = laim_table(vegtype,:) !monthly leaf area index, one-sided
parameters%sla = sla_table(vegtype) !single-side leaf area per kg [m2/kg]
Expand Down

0 comments on commit fb99cae

Please sign in to comment.