fractdew.f90

来自「CLM集合卡曼滤波数据同化算法」· F90 代码 · 共 52 行

  SUBROUTINE fractdew ( sigf,lai,sai,dewmx,ldew,fwet,fdry )
       
!=======================================================================
!      Source file: fractdew.f90
! Original version: Yongjiu Dai, September 15, 1999
!
! determine fraction of foliage covered by water and
! fraction of foliage that is dry and transpiring
!
!=======================================================================

  IMPLICIT NONE

! Dummy argument
  real, INTENT(in) :: &
        sigf, &! fraction of veg cover, excluding snow-covered veg [-]
        lai,  &! leaf area index  [-]
        sai,  &! stem area index  [-]
        dewmx,&! maximum allowed dew [0.1 mm]
        ldew   ! depth of water on foliage [kg/m2/s]

  real, INTENT(out) :: &
        fwet, &! fraction of foliage covered by water [-]
        fdry   ! fraction of foliage that is green and dry [-]
! Local
  real lsai,  &! lai + sai
       dewmxi,&! inverse of maximum allowed dew [1/mm]
       vegt    ! sigf*lsai
!
!-----------------------------------------------------------------------
! Fwet is the fraction of all vegetation surfaces which are wet 
! including stem area which contribute to evaporation
      lsai = lai + sai
      dewmxi  = 1.0/dewmx
      vegt   =  sigf*lsai

      fwet = 0
      if(ldew > 0.) then
         fwet = ((dewmxi/vegt)*ldew)**.666666666666

! Check for maximum limit of fwet
         fwet = min(fwet,1.0)

      end if

! fdry is the fraction of lai which is dry because only leaves can 
! transpire. Adjusted for stem area which does not transpire
      fdry = (1.-fwet)*lai/lsai
!
!-----------------------------------------------------------------------
  END SUBROUTINE fractdew