isba_fluxes_meb.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     ##########################################################################
      SUBROUTINE isba_fluxes_meb(KK, PK, PEK, DK, DEK, DMK, PRHOA, PLTT, PSIGMA_F,PSIGMA_FN, &
                                 PRN_V, PRN_G, PLWNET_V_DTV, PLWNET_V_DTG, PLWNET_V_DTN, &
                                 PLWNET_G_DTV, PLWNET_G_DTG, PLWNET_G_DTN, PLWNET_N_DTV, &
                                 PLWNET_N_DTG, PLWNET_N_DTN, PTHRMA_TA, PTHRMB_TA,       &
                                 PTHRMA_TC, PTHRMB_TC, PTHRMA_TG, PTHRMB_TG, PTHRMA_TV,  &
                                 PTHRMB_TV, PTHRMA_TN, PTHRMB_TN, PQSAT_G, PQSAT_V,      &
                                 PQSATI_N, PPSNA, PPSNCV, PFROZEN1,PLEG_DELTA,           &
                                 PLEGI_DELTA, PHUGI, PHVG, PHVN, PFLXC_CA, PFLXC_GV,     &
                                 PFLXC_VG_C, PFLXC_VN_C, PFLXC_GN, PFLXC_N_A, PFLXC_MOM, &
                                 PFLXC_CV, PHVGS, PHVNS, PTG, PDQSAT_G, PDQSAT_V,        &
                                 PDQSATI_N, PTA_IC, PQA_IC, PDELTA_V, PDELTAT_G,         &
                                 PDELTAT_V, PDELTAT_N, PSW_UP, PSW_RAD, PLW_RAD, PLW_UP, &
                                 PH_N_A, PEVAP_C_A, PEVAP_N_A, PLEG, PLEGI, PLES, PLEL,  &
                                 PEVAPN, PEMIS                             )
!     ##########################################################################
!
!!****  *ISBA_FLXUES_MEB*  
!!
!!    PURPOSE
!!    -------
!
!     Calculates the implicit fluxes for implicit or explicit atmospheric
!     coupling and fluxes needed by hydrology, soil and snow routines.
!     finally, compute soil phase changes.
!         
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    none
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!
!!      
!!    REFERENCE
!!    ---------
!!
!!    Noilhan and Planton (1989)
!!    Belair (1995)
!!    * to be done * (2011)
!!      
!!    AUTHOR
!!    ------
!!
!!      A. Boone           * Meteo-France *
!!      P. Samuelsson      * SMHI *
!!      S. Gollvik         * SMHI * 
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    22/01/11 
!!
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_isba_n, ONLY : isba_k_t, isba_pe_t, isba_p_t
USE modd_diag_n, ONLY : diag_t
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
USE modd_diag_misc_isba_n, ONLY : diag_misc_isba_t
!
USE modd_isba_par,       ONLY : xemissoil, xemisveg
USE modd_csts,           ONLY : xlvtt, xlstt, xstefan
!
USE modi_isba_emis_meb
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
TYPE(isba_k_t), INTENT(INOUT) :: KK
TYPE(isba_p_t), INTENT(INOUT) :: PK
TYPE(isba_pe_t), INTENT(INOUT) :: PEK
TYPE(diag_t), INTENT(INOUT) :: DK
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEK
TYPE(diag_misc_isba_t), INTENT(INOUT) :: DMK
!
!*      0.1    declarations of arguments
!
REAL, DIMENSION(:),   INTENT(IN)   :: PRHOA, PLTT
!                                     PRHOA = reference level air density (kg m-3)
!                                     PLTT  = latent heat normalization factor (J kg-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PSIGMA_F, PSIGMA_FN
!
REAL, DIMENSION(:),   INTENT(IN)   :: PLWNET_V_DTV, PLWNET_V_DTG, PLWNET_V_DTN
!                                     PLWNET_V_DTV, PLWNET_V_DTG, PLWNET_V_DTN = Vegetation canopy net LW radiation 
!                                     derivatives w/r/t surface temperature(s) (W m-2 K-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PLWNET_G_DTV, PLWNET_G_DTG, PLWNET_G_DTN
!                                     PLWNET_G_DTV, PLWNET_G_DTG, PLWNET_G_DTN = Understory-ground net LW radiation 
!                                          derivatives w/r/t surface temperature(s) (W m-2 K-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PLWNET_N_DTV, PLWNET_N_DTG, PLWNET_N_DTN
!                                     PLWNET_N_DTV, PLWNET_N_DTG, PLWNET_N_DTN = Ground-based snow net LW radiation 
!                                          derivatives w/r/t surface temperature(s) (W m-2 K-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PTHRMA_TA, PTHRMB_TA, PTHRMA_TC, PTHRMB_TC,                     &
                                      PTHRMA_TG, PTHRMB_TG, PTHRMA_TV, PTHRMB_TV, PTHRMA_TN, PTHRMB_TN
!                                     PTHRMA_TA                                                    (J kg-1 K-1)
!                                     PTHRMB_TA = linear transform coefficinets for atmospheric
!                                                 thermal variable for lowest atmospheric level.   (J kg-1)
!                                                 Transform T to dry static energy or enthalpy.
!                                     PTHRMA_TC                                                    (J kg-1 K-1)
!                                     PTHRMB_TC = linear transform coefficinets for atmospheric
!                                                 thermal variable for canopy air                  (J kg-1)
!                                                 Transform T to dry static energy or enthalpy.
!                                     PTHRMA_TG,V,N                                                (J kg-1 K-1)
!                                     PTHRMB_TG,V,N = linear transform coefficinets for atmospheric
!                                                 thermal variable for surfaces (G, V, and N)      (J kg-1)
!                                                 Transform T to dry static energy or enthalpy.
!
REAL, DIMENSION(:),   INTENT(IN)   :: PQSAT_G, PQSAT_V, PQSATI_N
!                                     PQSAT_G  = saturation specific humidity for understory surface    (kg kg-1)
!                                     PQSAT_V  = saturation specific humidity for the vegetation canopy (kg kg-1)
!                                     PQSATI_N = saturation specific humidity over ice for the snowpack (kg kg-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PPSNA, PPSNCV, PFROZEN1
!                                     PPSNA    = fraction of vegetation canopy buried by ground-based snowpack (-)
!                                     PPSNCV   = fraction of vegetation canopy covered by intercepted snow     (-)
!                                     PFROZEN1 = frozen fraction of surface ground layer                       (-)
!
!
REAL, DIMENSION(:),   INTENT(IN)   :: PLEG_DELTA, PLEGI_DELTA, PHUGI, PHVG, PHVN
!
!                                     PHVG = Halstead coefficient of non-buried (snow) canopy vegetation       (-)                         
!                                     PHVN = Halstead coefficient of paritally-buried (snow) canopy vegetation (-)                         
!
REAL, DIMENSION(:),   INTENT(IN)   :: PFLXC_CA, PFLXC_GV, PFLXC_VG_C, PFLXC_VN_C, PFLXC_GN, PFLXC_N_A,   &
                                      PFLXC_CV, PFLXC_MOM
!                                     PFLXC_CA  = Flux form heat transfer coefficient: canopy air to atmosphere (kg m-2 s-1)
!                                     PFLXC_GV  = As above, but for : ground-understory to canopy air           (kg m-2 s-1)
!                                     PFLXC_VG_C = As above, but for : non-snow buried canopy to canopy air      (kg m-2 s-1)
!                                     PFLXC_VN_C = As above, but for : partially snow-buried canopy air to canopy 
!                                                  air                                                           (kg m-2 s-1)
!                                     PFLXC_CV  = As above, but for : bulk vegetation canopy to canopy air      (kg m-2 s-1)
!                                     PFLXC_GN  = As above, but for : ground-based snow to atmosphere           (kg m-2 s-1)
!                                     PFLXC_N_A  = As above, but for : ground-based snow to canopy air           (kg m-2 s-1)
!                                     PFLXC_MOM  = flux form drag transfer coefficient: canopy air to atmosphere (kg m-2 s-1)
!
REAL, DIMENSION(:,:), INTENT(IN)   :: PTG
!                                     PTG    = Soil temperature profile (K)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PDQSAT_G, PDQSAT_V, PDQSATI_N
!                                     PQSAT_G  = saturation specific humidity derivative for understory 
!                                                surface               (kg kg-1 K-1)
!                                     PQSAT_V  = saturation specific humidity derivative for the vegetation 
!                                                canopy                (kg kg-1 K-1)
!                                     PQSATI_N = saturation specific humidity derivative over ice for the 
!                                                ground-based snowpack (kg kg-1 K-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PHVGS, PHVNS
!                                     PHVGS     = Dimensionless pseudo humidity factor for computing vapor
!                                                 fluxes from the non-buried part of the canopy to the canopy air    (-)
!                                     PHVNS     = Dimensionless pseudo humidity factor for computing vapor
!                                                 fluxes from the partly-buried part of the canopy to the canopy air (-)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PTA_IC, PQA_IC
!                                     PTA_IC    = Near-ground air temperature        (K)
!                                     PQA_IC    = Near-ground air specific humidity  (kg kg-1)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PSW_UP, PSW_RAD, PLW_RAD
!                                     PSW_UP  = total upwelling shortwave radiation from the surface at the atmosphere (W m-2)
!                                     PSW_RAD = downwelling shortwave radiation from the atmosphere above the canopy  (W m-2)
!                                     PLW_RAD = downwelling longwave radiation from the atmosphere above the canopy  (W m-2)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PDELTA_V
!                                     PDELTA_V = Explicit canopy interception fraction                 (-)
!
REAL, DIMENSION(:),   INTENT(IN)   :: PDELTAT_V, PDELTAT_N, PDELTAT_G
!                                     PDELTAT_V = Time change in vegetation canopy temperature (K)
!                                     PDELTAT_N = Time change in snowpack surface temperature  (K)
!                                     PDELTAT_G = Time change in soil surface temperature      (K)
!
REAL, DIMENSION(:),   INTENT(OUT)  :: PRN_V, PRN_G
!                                     PRN_G = Understory-ground net radiation (W m-2)
!                                     PRN_V = Vegetation canopy net radiation (W m-2)
!
REAL, DIMENSION(:),   INTENT(OUT)  :: PLW_UP  
!                                     PLW_UP  = total net longwave upwelling radiation to the atmosphere  (W m-2)
!
REAL, DIMENSION(:),   INTENT(OUT)  :: PH_N_A
!                                     PH_N_A  = Sensible heat flux: ground based snowpack to overlying atmosphere (W m-2)
!
REAL, DIMENSION(:),   INTENT(OUT)  :: PEVAP_C_A,  PEVAP_N_A, PEVAPN
!                                     PEVAP_C_A = Water flux: canopy air space to overlying atmosphere (kg m-2 s-1)
!                                     PEVAP_N_A = Water flux: ground based snowpack to overlying atmosphere (kg m-2 s-1)
!                                     PEVAPN    = Water flux: ground based snowpack to both canopy air space and overlying atmosphere (kg m-2 s-1)
!
REAL, DIMENSION(:),   INTENT(OUT)  :: PLEG, PLEGI, PLES, PLEL
!                                     PLEG        = Latent heat flux: baresoil evaporation (W m-2)
!                                     PLEGI       = Latent heat flux: baresoil sublimation (W m-2)
!                                     PLES        = Latent heat flux: net sublimation from ground-based snowpack to canopy air and overlying atmosphere (W m-2)
!                                     PLEL        = Latent heat flux: net evaporation from ground-based snowpack to canopy air and overlying atmosphere (W m-2)
!
REAL, DIMENSION(:),   INTENT(OUT)  :: PEMIS
!                                     PEMIS       = effective (aggregated) net surface emissivity (-)
!
!
!*      0.2    declarations of local variables
!
!
REAL, DIMENSION(SIZE(PEK%XTV(:)))         :: ZWORK
!                                     ZWORK = working array
!
REAL, DIMENSION(SIZE(PEK%XTV(:)))         :: ZSAIR, ZSAIRC
!                                     ZSAIR   = atmospheric value of the therodynamic variable
!                                     ZSAIRC  = canopy air value of the therodynamic variable
!
REAL, DIMENSION(SIZE(PEK%XTV(:)))         :: ZEVAP_CV
!                                     ZEVAP_CV = Water flux: Evapotranspiration vapor flux from the vegetation canopy (kg m-2 s-1)
!
REAL, DIMENSION(SIZE(PEK%XTV(:)))         :: ZQSATN_V, ZQSATIN_N, ZQSATN_G
!                                     ZQSATN_V  = saturation specific humidity (over water) for the vegetation canopy (kg kg-1)
!                                     ZQSATIN_N = saturation specific humidity (over ice) for the snow (kg kg-1)
!                                                 NOTE that liquid water can only exist when the snowpack T=XTT in the model, 
!                                                 and at the freezing point, the value is the same over ice and water, therefore
!                                                 over snow, we do not need to explicitly consider a "ZQSATN_N"
!                                     ZQSATN_G  = saturation specific humidity (over water) for the understory (kg kg-1)
!                
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
!*       0.     Initialization:
!               ---------------
!
IF (lhook) CALL dr_hook('ISBA_FLUXES_MEB',0,zhook_handle)
!-------------------------------------------------------------------------------
!
!*       1.     Radiative Fluxes
!               ----------------
!
! LWnet: transform from explicit to implicit (i.e. at time t+dt)
!
dek%XLWNET_V(:)  = dek%XLWNET_V(:) + plwnet_v_dtv(:)*pdeltat_v(:)   &
                                   + plwnet_v_dtg(:)*pdeltat_g(:)   &
                                   + plwnet_v_dtn(:)*pdeltat_n(:)

dek%XLWNET_G(:)  = dek%XLWNET_G(:) + plwnet_g_dtv(:)*pdeltat_v(:)   &
                                   + plwnet_g_dtg(:)*pdeltat_g(:)   &
                                   + plwnet_g_dtn(:)*pdeltat_n(:)

dek%XLWNET_N(:)  = dek%XLWNET_N(:) + plwnet_n_dtv(:)*pdeltat_v(:)   &
                                   + plwnet_n_dtg(:)*pdeltat_g(:)   &
                                   + plwnet_n_dtn(:)*pdeltat_n(:)
!
! LWup at t+dt
!
plw_up(:)   = plw_rad(:) - (dek%XLWNET_V(:) + dek%XLWNET_G(:) + dek%XLWNET_N(:))
!
!
! Effective emissivity:
!
 CALL isba_emis_meb(pek%XPSN, ppsna, psigma_f, psigma_fn, pek%TSNOW%EMIS, pemis   )
!
! Now compute the effective radiative temperature while
! imposing the constraint: 
!
!    LW_RAD * (1 - EMIS ) + EMIS * XSTEFAN * TS_RAD**4 = LWUP
!
! Using the effective emissivity ensures that the upwelling radiation from the surface (RHS)
! model will be equal to the upwelling radiation computed in the atmospheric model (LHS)
! (i.e. LWUP is consistent with EMIS & TS_RAD), thereby insuring energy conservation from
! the surface to the atmosphere. Solving the above equation for
! the radiative T gives:
!
dk%XTSRAD(:)    = ((plw_up(:) - plw_rad(:)*(1.0-pemis(:)))/(xstefan*pemis(:)))**0.25
!
!
! Rnet (t+dt)
!
prn_v(:)  = dek%XSWNET_V(:) + dek%XLWNET_V(:)
!
prn_g(:)  = dek%XSWNET_G(:) + dek%XLWNET_G(:)
!
dmk%XRNSNOW(:) = dek%XSWNET_N(:) + dek%XLWNET_N(:)
!
!
! total Rnet (t+dt):
!
dk%XRN(:)    = prn_g(:) + prn_v(:) + dmk%XRNSNOW(:)
!
!
!*       2.a    Implicit (Turbulent) Sensible Heat Fluxes
!               -----------------------------------------

! First get input thermo variable (could be enthalpy (air heat capacity x potential temperature or dry static energy)

zsair(:)  = pthrmb_ta(:) + pthrma_ta(:)*pta_ic(:)
zsairc(:) = pthrmb_tc(:) + pthrma_tc(:)*pek%XTC(:)

! Sensible heat fluxes (W m-2):
! - Canopy air to atmosphere, vegetation canopy to canopy air (implicitly includes from canopy intercepted snow), 
!   understory-ground to canopy air,
!   ground-based snow to canopy air, ground-based snow to atmosphere:

dek%XH_CA(:)  = pflxc_ca(:) *( zsairc(:)                                      - zsair(:))*(1.0 - pek%XPSN(:)*ppsna(:))
dek%XH_CV(:)  = pflxc_cv(:) *( pthrmb_tv(:) + pthrma_tv(:)*pek%XTV(:)         - zsairc(:)) 
dek%XH_GV(:)  = pflxc_gv(:) *( pthrmb_tg(:) + pthrma_tg(:)*ptg(:,1)           - zsairc(:))*(1.0 - pek%XPSN(:))
dek%XH_GN(:)  = pflxc_gn(:) *( pthrmb_tn(:) + pthrma_tn(:)*dmk%XSNOWTEMP(:,1) - zsairc(:))*       pek%XPSN(:) *(1.0-ppsna(:)) 
ph_n_a(:)  = pflxc_n_a(:) *( pthrmb_tn(:) + pthrma_tn(:)*dmk%XSNOWTEMP(:,1) - zsair(:))*       pek%XPSN(:) *     ppsna(:) 

! - Net sensible heat flux from ground-based snow (to the canopy and the atmosphere (from
!   the buried-vegetation canopy fraction)) (W m-2) 

dmk%XHSNOW(:) = dek%XH_GN(:) + ph_n_a(:) 

! FINAL sensible heat flux to the atmosphere (W m-2):

dk%XH(:)      = dek%XH_CA(:) + ph_n_a(:)

!
!*       2.b    Implicit (Turbulent) Vapor and Latent Heat Fluxes
!               -------------------------------------------------
! Note, to convert any of the latent heat fluxes back to vapor fluxes, 
! simply divide by XLVTT, even sublimation fluxes as XLSTT already accounted for.

! - first get 'new' surface specific humidities, qsatn, at time t+dt:

zqsatn_g(:)  =        pqsat_g(:)  + pdqsat_g(:)  * pdeltat_g(:)
zqsatn_v(:)  =        pqsat_v(:)  + pdqsat_v(:)  * pdeltat_v(:)
zqsatin_n(:) =        pqsati_n(:) + pdqsati_n(:) * pdeltat_n(:)

! - Evaporation and Sublimation latent heat fluxes from the soil, respectively:
! (W m-2)

zwork(:)      = (1.-pek%XPSN(:)-kk%XFF(:)) * pflxc_gv(:)
pleg(:)       = zwork(:)*pleg_delta(:) *( dk%XHUG(:) *zqsatn_g(:) - pek%XQC(:) )*(1.-pfrozen1(:))*xlvtt
plegi(:)      = zwork(:)*plegi_delta(:)*(   phugi(:) *zqsatn_g(:) - pek%XQC(:) )*    pfrozen1(:) *xlstt

! - Latent heat flux from frozen and unfrozen flooded zones (W m-2)

zwork(:)      = kk%XFF(:) * pflxc_gv(:)*( zqsatn_g(:) - pek%XQC(:) )
dek%XLE_FLOOD(:)  = zwork(:) * (1.-kk%XFFROZEN(:))* xlvtt 
dek%XLEI_FLOOD(:) = zwork(:) *     kk%XFFROZEN(:) * xlstt 

! - Evapotranspiration vapor flux from the vegetation canopy (kg m-2 s-1)

zevap_cv(:) = (1.-ppsncv(:)) * phvgs(:) * pflxc_cv(:)*( zqsatn_v(:) - pek%XQC(:) ) * (xlvtt/pltt(:))

! - Latent heat flux from the canopy (liquid) water interception reservoir (W m-2)

dek%XLER_CV(:) = ( (1.-ppsna(:))*pek%XPSN(:) * pflxc_vn_c(:)    +              &
                             (1.-pek%XPSN(:))* pflxc_vg_c(:)  ) *              &
                 xlvtt * (1.-ppsncv(:))* pdelta_v(:) * ( zqsatn_v(:) - pek%XQC(:) )
!
! - Total latent heat flux (evapotranspiration) from the vegetation to the canopy air space (W m-2)
!   *without* sublimation (for TOTAL evapotranspiration and sublimation, add PLESC here)

dek%XLEV_CV(:)  = pltt(:) * zevap_cv(:)
!
! - latent heat flux from transpiration from the canopy (W m-2)

dek%XLETR_CV(:) = dek%XLEV_CV(:) - dek%XLER_CV(:) 

! Snow sublimation and evaporation latent heat flux from canopy-intercepted snow (W m-2)

dek%XLES_CV(:)  =  ppsncv(:) * xlstt * phvns(:) * pflxc_cv(:)*( zqsatn_v(:) - pek%XQC(:) )

! - Total latent heat flux from vegetation canopy overstory to canopy air space
!   (including transpiration, liquid water store, canopy snow sublimation):

dek%XLE_CV(:)   = dek%XLEV_CV(:) + dek%XLES_CV(:)

! - Vapor flux from the ground-based snowpack to the canopy air (kg m-2 s-1):

zwork(:)        = pflxc_gn(:)*(zqsatin_n(:) - pek%XQC(:))*pek%XPSN(:)*(1.0-ppsna(:))
dek%XEVAP_GN(:) = zwork(:)*(xlstt/pltt(:))
dek%XLE_GN(:)   = zwork(:)* xlstt                              ! W m-2

! - latent heat flux from transpiration from canopy veg (evapotranspiration)

dek%XLETR(:)     = dek%XLETR_CV(:)

! Total latent heat flux from transpiration from understory veg and canopy veg (evapotranspiration and sublimation)
!   and intercepted water on both reservoirs (W m-2) 

dek%XLEV(:)      = dek%XLETR(:) + dek%XLER_CV(:) 

! Total latent heat flux from intercepted water (canopy and understory vegetation):
! (does not include intercepted snow sublimation): W m-2

dek%XLER(:)      = dek%XLER_CV(:)

! - Vapor flux from the ground-based snowpack (part burying the canopy vegetation) to the atmosphere (kg m-2 s-1):

zwork(:)     = pflxc_n_a(:) *( zqsatin_n(:) - pqa_ic(:))*       pek%XPSN(:)*     ppsna(:)
pevap_n_a(:) = zwork(:) *(xlstt/pltt)

! - Net Snow (groud-based) sublimation latent heat flux (W m-2) to the canopy air space and the overlying atmosphere:
ples(:)      = ( pflxc_gn(:) *( zqsatin_n(:) - pek%XQC(:))*       pek%XPSN(:)*(1.0-ppsna(:)) + zwork(:) ) * xlstt

! - Net Snow evaporation (liquid water) latent heat flux (W m-2)

plel(:)      = pltt(:)*(dek%XEVAP_GN(:) + pevap_n_a(:)) - ples(:)

! - Total mass flux from ground-based snowpack (kg m-2 s-1)

pevapn(:)    = (plel(:) + ples(:))/pltt(:)

! - Total snow-free vapor flux from the understory (flooded areas, baresoil and understory vegetation)
!   to the canopy air space (W m-2 and kg m-2 s-1, respectively):

dek%XLE_GV(:)   = dek%XLE_FLOOD(:) + dek%XLEI_FLOOD(:) + plegi(:) + pleg(:) 

dek%XEVAP_G(:)  = dek%XLE_GV(:)/pltt(:) 

! - Net vapor flux from canopy air to the atmosphere (kg m-2 s-1)

pevap_c_a(:) = pflxc_ca(:) *( pek%XQC(:) - pqa_ic(:))*(1.0 - pek%XPSN(:)*ppsna(:))

dek%XLE_CA(:)   = pltt(:) * pevap_c_a(:) ! W m-2

! FINAL net vapor flux from the surface to the Atmosphere:
! - Net vapor flux from canopy air and exposed ground based snow (from part of snow 
!   burying the vegetation canopy) to the atmosphere (kg m-2 s-1)
!
dk%XEVAP(:)     = pevap_c_a(:) + pevap_n_a(:)   
! 
! Total latent heat flux of surface/snow/vegetation: W m-2
!
pek%XLE(:)       = dk%XEVAP(:)*pltt(:)                  
!
! Total sublimation from the surface/snow/vegetation: W m-2
!
dk%XLEI(:)      = ples(:) + plegi(:) + dek%XLEI_FLOOD(:)
!
! Total sublimation from the surface/snow/vegetation: kg m-2 s-1
!
dk%XSUBL(:)     = dk%XLEI(:)/ pltt(:)
!
IF (lhook) CALL dr_hook('ISBA_FLUXES_MEB',1,zhook_handle)
!
END SUBROUTINE isba_fluxes_meb