Skip to content

Commit

Permalink
Bug fix in physics/module_sf_noahmp_glacier.F90: only compute glacier…
Browse files Browse the repository at this point in the history
… albedo in snowalb_bats_glacier if cosz > 0
  • Loading branch information
climbfuji committed Oct 8, 2024
1 parent da75531 commit 78889db
Showing 1 changed file with 14 additions and 12 deletions.
26 changes: 14 additions & 12 deletions physics/SFC_Models/Land/Noahmp/module_sf_noahmp_glacier.F90
Original file line number Diff line number Diff line change
Expand Up @@ -971,18 +971,20 @@ subroutine snowalb_bats_glacier (nband,cosz,fage,albsnd,albsni)
albsni(1: nband) = 0.

! when cosz > 0

sl=2.0
sl1=1./sl
sl2=2.*sl
cf1=((1.+sl1)/(1.+sl2*cosz)-sl1)
fzen=amax1(cf1,0.)

albsni(1)=0.95*(1.-c1*fage)
albsni(2)=0.65*(1.-c2*fage)

albsnd(1)=albsni(1)+0.4*fzen*(1.-albsni(1)) ! vis direct
albsnd(2)=albsni(2)+0.4*fzen*(1.-albsni(2)) ! nir direct
! https://github.com/ufs-community/ccpp-physics/issues/227
if (cosz>0) then
sl=2.0
sl1=1./sl
sl2=2.*sl
cf1=((1.+sl1)/(1.+sl2*cosz)-sl1)
fzen=amax1(cf1,0.)

albsni(1)=0.95*(1.-c1*fage)
albsni(2)=0.65*(1.-c2*fage)

albsnd(1)=albsni(1)+0.4*fzen*(1.-albsni(1)) ! vis direct
albsnd(2)=albsni(2)+0.4*fzen*(1.-albsni(2)) ! nir direct
endif

end subroutine snowalb_bats_glacier
! ==================================================================================================
Expand Down

0 comments on commit 78889db

Please sign in to comment.