averaged_albedo_emis_isba.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 averaged_albedo_emis_isba (IO, S, NK, NP, NPE, &
                                 PZENITH, PTG1, PSW_BANDS, PDIR_ALB, PSCA_ALB, &
                                 PEMIS, PTSRAD, PTSURF, PDIR_SW, PSCA_SW        )
!     ###################################################
!
!!**** ** computes radiative fields used in ISBA
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    01/2004
!!     A. Bogatchev 09/2005 EBA snow option
!!     B. Decharme  2008    The fraction of vegetation covered by snow must be
!                            <= to ZSNG
!!     B. Decharme  2013    new coupling variable and optimization    
!!     P. Samuelsson 10/2014 MEB
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE modd_isba_options_n, ONLY : isba_options_t
USE modd_isba_n, ONLY: isba_s_t, isba_nk_t, isba_np_t, isba_npe_t, isba_k_t, &
                       isba_p_t, isba_pe_t
!
USE modd_surf_par,  ONLY : xundef
!
USE modd_type_snow
!
USE modd_csts,      ONLY : xstefan
USE mode_meb,       ONLY : mebpalphan
!
USE modi_albedo
USE modi_average_rad
USE modi_update_rad_isba_n
USE modi_isba_lwnet_meb
USE modi_unpack_same_rank
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(isba_options_t), INTENT(INOUT) :: IO
TYPE(isba_s_t), INTENT(INOUT) :: S
TYPE(isba_nk_t), INTENT(INOUT) :: NK
TYPE(isba_np_t), INTENT(INOUT) :: NP
TYPE(isba_npe_t), INTENT(INOUT) :: NPE
!
REAL, DIMENSION(:,:),   INTENT(IN)   :: PTG1        ! soil surface temperature
REAL, DIMENSION(:),     INTENT(IN)   :: PZENITH     
REAL, DIMENSION(:),     INTENT(IN)   :: PSW_BANDS   ! middle wavelength of each band
!
REAL, DIMENSION(:,:),   INTENT(OUT)  :: PDIR_ALB    ! averaged direct albedo  (per wavelength)
REAL, DIMENSION(:,:),   INTENT(OUT)  :: PSCA_ALB    ! averaged diffuse albedo (per wavelength)
REAL, DIMENSION(:),     INTENT(OUT)  :: PEMIS       ! averaged emissivity
REAL, DIMENSION(:),     INTENT(OUT)  :: PTSRAD      ! averaged radiaitve temp.
REAL, DIMENSION(:),     INTENT(OUT)  :: PTSURF      ! surface effective temperature         (K)
!
REAL, DIMENSION(:,:),   INTENT(IN), OPTIONAL   :: PDIR_SW ! Downwelling direct SW radiation
REAL, DIMENSION(:,:),   INTENT(IN), OPTIONAL   :: PSCA_SW ! Downwelling diffuse SW radiation
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
!
TYPE(isba_k_t), POINTER :: KK
TYPE(isba_p_t), POINTER :: PK
TYPE(isba_pe_t), POINTER :: PEK
!
REAL, DIMENSION(SIZE(PZENITH),SIZE(PSW_BANDS),IO%NPATCH) :: ZDIR_ALB_PATCH 
!                                                     ! direct albedo
REAL, DIMENSION(SIZE(PZENITH),SIZE(PSW_BANDS),IO%NPATCH) :: ZSCA_ALB_PATCH 
!                                                     ! diffuse albedo
REAL, DIMENSION(SIZE(PZENITH),IO%NPATCH) :: ZEMIS_PATCH   ! emissivity with snow-flood
REAL, DIMENSION(SIZE(PZENITH),IO%NPATCH) :: ZTSRAD_PATCH  ! Tsrad
REAL, DIMENSION(SIZE(PZENITH),IO%NPATCH) :: ZTSURF_PATCH  ! Tsurf
REAL, DIMENSION(SIZE(PTG1,1)) :: ZEMIS         ! emissivity with flood
!
REAL, DIMENSION(SIZE(PTG1,1)) :: ZSNOWDEPTH    ! Total snow depth
REAL, DIMENSION(SIZE(PTG1,1)) :: ZPALPHAN      ! Snow/canopy ratio factor 
REAL, DIMENSION(SIZE(PTG1,1)) :: ZLW_RAD       ! Fake downwelling LW rad
REAL, DIMENSION(SIZE(PTG1,1)) :: ZLW_UP        ! Upwelling LW rad
REAL, DIMENSION(SIZE(PTG1,1)) :: ZLWNET_N      ! LW net for snow surface
REAL, DIMENSION(SIZE(PTG1,1)) :: ZLWNET_V      ! LW net for canopy veg
REAL, DIMENSION(SIZE(PTG1,1)) :: ZLWNET_G      ! LW net for ground
REAL, DIMENSION(SIZE(PTG1,1)) :: ZDUMMY
REAL, DIMENSION(SIZE(PTG1,1)) :: ZEMISF
REAL, DIMENSION(SIZE(PTG1,1)) :: ZFF
!
LOGICAL :: LEXPLICIT_SNOW ! snow scheme key
!
INTEGER :: JP, JI,ISIZE ! loop on patches
INTEGER :: IMASK ! loop on patches
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
!*    0.      Init
!             ----
!
IF (lhook) CALL dr_hook('AVERAGED_ALBEDO_EMIS_ISBA',0,zhook_handle)
!
pdir_alb(:,:)=0.
psca_alb(:,:)=0.
pemis(:)  =0.
ptsrad(:)  =0.
ptsurf(:)  =0.
!
zdir_alb_patch(:,:,:)=0.
zsca_alb_patch(:,:,:)=0.
zemis_patch(:,:  )=0.
!
lexplicit_snow = (npe%AL(1)%TSNOW%SCHEME=='3-L'.OR.npe%AL(1)%TSNOW%SCHEME=='CRO')
!
ztsrad_patch(:,:) = 0.
ztsurf_patch(:,:) = 0.
DO jp = 1,io%NPATCH
  DO ji = 1,np%AL(jp)%NSIZE_P
    imask = np%AL(jp)%NR_P(ji)
    ztsrad_patch(imask,jp) = ptg1(ji,jp)
    ztsurf_patch(imask,jp) = ptg1(ji,jp)
  ENDDO
ENDDO
!
!
!*    1.      averaged albedo on natural continental surfaces (except prognostic snow)
!             -----------------------------------------------
!
DO jp = 1,io%NPATCH
  CALL albedo(io%CALBEDO, npe%AL(jp) ) 
ENDDO
!
!*    2.      averaged albedo and emis. on natural continental surfaces (with prognostic snow)
!             ---------------------------------------------------------
!
! A dummy downwelling LW radiation can be used for calculation of radiative surface temp 
!
zlw_rad(:) = 300.0
!    
!* Initialization of albedo for each wavelength, emissivity and snow/flood fractions
!
DO jp = 1,io%NPATCH
  !
  IF(PRESENT(pdir_sw))THEN
    !
    ! For the case when MEB patch albedo is requested downweeling SW is needed
    !
    CALL update_rad_isba_n(io, s, nk%AL(jp), np%AL(jp), npe%AL(jp), jp, pzenith, psw_bands,   &
                           zdir_alb_patch(:,:,jp), zsca_alb_patch(:,:,jp), zemis_patch(:,jp), &
                           pdir_sw, psca_sw    )
  ELSE
    !
    ! For cases when MEB patch albedo is not requested no downweeling SW is needed
    !
    CALL update_rad_isba_n(io, s, nk%AL(jp), np%AL(jp), npe%AL(jp), jp, pzenith, psw_bands, &
                           zdir_alb_patch(:,:,jp), zsca_alb_patch(:,:,jp), zemis_patch(:,jp))
    !
  ENDIF
  !
ENDDO
!
!
!* radiative surface temperature
!
DO jp  =  1,io%NPATCH
  !
  pek => npe%AL(jp)
  pk => np%AL(jp)
  kk => nk%AL(jp)
  !
  isize = pk%NSIZE_P
  !
  IF(io%LMEB_PATCH(jp))THEN  ! MEB patches
    !
    !   ZPALPHAN is needed as input to ISBA_LWNET_MEB
    !
    zsnowdepth(1:isize) = sum(pek%TSNOW%WSNOW(:,:)/pek%TSNOW%RHO(:,:),2)
    zpalphan(1:isize) = mebpalphan(zsnowdepth(1:isize),pek%XH_VEG(:))
    !
    !   ZLWNET_N,ZLWNET_V,ZLWNET_G are needed for ZLW_UP and ZTSRAD_PATCH
    !
    IF(io%LFLOOD)THEN
      zemisf(1:isize) = kk%XEMISF(:)
      zff(1:isize) = kk%XFF   (:)
    ELSE
      zemisf(1:isize) = xundef
      zff(1:isize) = 0.0
    ENDIF
    !
    CALL isba_lwnet_meb(pek%XLAI, pek%XPSN, zpalphan(1:isize), pek%TSNOW%EMIS, &
                        zemisf(1:isize), zff(1:isize),           &
                        pek%XTV, ptg1(1:isize,jp), pek%TSNOW%TS,       &
                        zlw_rad(1:isize), zlwnet_n(1:isize), zlwnet_v(1:isize), zlwnet_g(1:isize),   &
                        zdummy(1:isize), zdummy(1:isize), zdummy(1:isize), zdummy(1:isize), &
                        zdummy(1:isize), zdummy(1:isize), zdummy(1:isize), zdummy(1:isize), &
                        zdummy(1:isize), zdummy(1:isize), zdummy(1:isize), zdummy(1:isize)   )
    !
    zlw_up(1:isize)  = zlw_rad(1:isize) - (zlwnet_v(1:isize) + zlwnet_g(1:isize) + zlwnet_n(1:isize))
    !
    !   MEB patch radiative temperature
    !
    DO ji = 1,pk%NSIZE_P
      imask = pk%NR_P(ji)
      IF (zemis_patch(imask,jp)/=0.) THEN
        ztsrad_patch(imask,jp) = ((zlw_up(ji) - zlw_rad(ji)*(1.0-zemis_patch(imask,jp)))/ &
                                 (xstefan*zemis_patch(imask,jp)))**0.25
      ENDIF
    END DO
    !
  ELSE   ! Non-MEB patches
    !
    zemis(1:isize) = pek%XEMIS(:)
    !
    IF(io%LFLOOD.AND.lexplicit_snow)THEN
      WHERE(pek%XPSN(:)<1.0.AND.pek%XEMIS(:)/=xundef)          
        zemis(1:isize) = ((1.-kk%XFF(:)-pek%XPSN(:))*pek%XEMIS(:) + kk%XFF(:)*kk%XEMISF(:)) /(1.-pek%XPSN(:))
      ENDWHERE   
    ENDIF
    !
    IF (.NOT.lexplicit_snow) THEN
      CALL unpack_same_rank(pk%NR_P,ptg1(1:pk%NSIZE_P,jp),ztsrad_patch(:,jp),0.)
    ELSE IF (lexplicit_snow) THEN
      DO ji = 1,pk%NSIZE_P
        imask = pk%NR_P(ji)
        IF (pek%XEMIS(ji)/=xundef .AND. zemis_patch(imask,jp)/=0.) THEN
          ztsrad_patch(imask,jp) =( ( (1.-pek%XPSN(ji)) * zemis(ji)*ptg1(ji,jp)**4            &
                                     +    pek%XPSN(ji) *pek%TSNOW%EMIS(ji)*pek%TSNOW%TS(ji)**4 ) )**0.25  &
                                 / zemis_patch(imask,jp)**0.25  
        END IF
      ENDDO
    END IF

  ENDIF
!
END DO
!
!* averaged radiative fields
!
 CALL average_rad(s%XPATCH,                                                  &
                  zdir_alb_patch, zsca_alb_patch, zemis_patch, ztsrad_patch, &
                  pdir_alb,       psca_alb,       pemis,       ptsrad        )  
!
!* averaged effective temperature
!
IF(lexplicit_snow)THEN
  DO jp = 1,io%NPATCH
    pek => npe%AL(jp)
    pk => np%AL(jp)
    DO ji = 1,pk%NSIZE_P
      imask = pk%NR_P(ji)
      ztsurf_patch(imask,jp) = ptg1(ji,jp)*(1.-pek%XPSN(ji)) + pek%TSNOW%TS(ji)*pek%XPSN(ji)
    ENDDO
  ENDDO
ENDIF
!
DO jp=1,io%NPATCH
  DO ji=1,np%AL(jp)%NSIZE_P
    imask = np%AL(jp)%NR_P(ji)
    ptsurf(imask) = ptsurf(imask) + np%AL(jp)%XPATCH(ji) * ztsurf_patch(imask,jp)
  ENDDO
ENDDO
!
IF (lhook) CALL dr_hook('AVERAGED_ALBEDO_EMIS_ISBA',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE averaged_albedo_emis_isba