modd_teb_greenroofn.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.
!##################
MODULE modd_teb_greenroof_n
!##################
!
!!****  *MODD_TEB_GREENROOF - declaration of ISBA scheme packed surface parameters for urban green roofs
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      A. Lemonsu *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       09/2009
!!      C. de Munck     06/2011 
!!      V. Masson       06/2013  splits module in 4
!!
!-------------------------------------------------------------------------------
!
!*       0.   DECLARATIONS
!             ------------
!
USE modd_type_snow
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE

TYPE teb_greenroof_1p_t
!-------------------------------------------------------------------------------
!
! Prognostic variables:
!
! - Snow Cover:
!
  TYPE(surf_snow)                :: TSNOW            ! snow state: 
                                                     ! scheme type/option                      (-)
                                                     ! number of layers                        (-)
                                                     ! snow (& liq. water) content             (kg/m2)
                                                     ! heat content                            (J/m2)
                                                     ! temperature                             (K)
                                                     ! density                                 (kg m-3)
!
! - Soil and vegetation heat and water:
!
  REAL, POINTER, DIMENSION(:)    :: XWR              ! liquid water retained on the
                                                     ! foliage of the vegetation
                                                     ! canopy                                  (kg/m2)
  REAL, POINTER, DIMENSION(:,:)  :: XTG              ! surface and sub-surface soil 
                                                     ! temperature profile                     (K)
  REAL, POINTER, DIMENSION(:,:)  :: XWG              ! soil volumetric water content profile   (m3/m3)
  REAL, POINTER, DIMENSION(:,:)  :: XWGI             ! soil liquid water equivalent volumetric 
                                                     ! ice content profile                     (m3/m3)
  REAL, POINTER, DIMENSION(:)    :: XRESA            ! aerodynamic resistance                  (s/m)
!

! - Vegetation: Ags Prognostic (YPHOTO = 'AGS', 'LAI', 'AST', 'LST', 'NIT', 'NCB')
!
  REAL, POINTER, DIMENSION(:)    :: XAN              ! net CO2 assimilation                    (mg/m2/s)
  REAL, POINTER, DIMENSION(:)    :: XANDAY           ! daily net CO2 assimilation              (mg/m2)
  REAL, POINTER, DIMENSION(:)    :: XANFM            ! maximum leaf assimilation               (mg/m2/s)
  REAL, POINTER, DIMENSION(:)    :: XLE              ! evapotranspiration                      (W/m2)
  REAL, POINTER, DIMENSION(:)     :: XFAPARC       ! Fapar of vegetation (cumul)
  REAL, POINTER, DIMENSION(:)     :: XFAPIRC       ! Fapir of vegetation (cumul)
  REAL, POINTER, DIMENSION(:)     :: XLAI_EFFC     ! Effective LAI (cumul)
  REAL, POINTER, DIMENSION(:)     :: XMUS          ! cos zenithal angle (cumul)    
!
! - Vegetation: Ags Prognostic (YPHOTO = 'NIT', 'NCB')
!
  REAL, POINTER, DIMENSION(:,:)  :: XRESP_BIOMASS    ! daily cumulated respiration of 
                                                     ! biomass                              (kg/m2/s)
  REAL, POINTER, DIMENSION(:,:)  :: XBIOMASS         ! biomass of previous day              (kg/m2) 
!
! - SGH scheme
!                                                     
  REAL, POINTER, DIMENSION(:)    :: XKSAT_ICE       ! hydraulic conductivity at saturation
                                                    ! over frozen area (m s-1)                                     
!-------------------------------------------------------------------------------
!
! - Snow and flood fractions and total albedo at time t:
!
  REAL, POINTER, DIMENSION(:)    :: XPSNG              ! Snow fraction over ground
  REAL, POINTER, DIMENSION(:)    :: XPSNV              ! Snow fraction over vegetation
  REAL, POINTER, DIMENSION(:)    :: XPSNV_A            ! Snow fraction over vegetation
  REAL, POINTER, DIMENSION(:)    :: XPSN               ! Total Snow fraction
! 
  REAL, POINTER, DIMENSION(:)    :: XSNOWFREE_ALB      ! snow free albedo                        (-)
  REAL, POINTER, DIMENSION(:)    :: XSNOWFREE_ALB_VEG  ! snow free albedo for vegetation         (-)
  REAL, POINTER, DIMENSION(:)    :: XSNOWFREE_ALB_SOIL ! snow free albedo for soil               (-)
!
!-------------------------------------------------------------------------------
!
END TYPE teb_greenroof_1p_t
!
TYPE teb_greenroof_t
  !
  TYPE(teb_greenroof_1p_t), POINTER :: ALP(:) => NULL()
  TYPE(teb_greenroof_1p_t), POINTER :: CUR => NULL()
  !
END TYPE teb_greenroof_t
!
 CONTAINS
!
!
SUBROUTINE teb_greenroof_goto_patch(YTEB_GREENROOF,KTO_PATCH)
        TYPE(teb_greenroof_t), INTENT(INOUT) :: yteb_greenroof
INTEGER, INTENT(IN) :: kto_patch
REAL(KIND=JPRB) :: zhook_handle
!
! Current patch is set to patch KTO_PATCH
IF (lhook) CALL dr_hook('MODD_TEB_GREENROOF_N:TEB_GREENROOF_GOTO_PATCH',0,zhook_handle)

yteb_greenroof%CUR => yteb_greenroof%ALP(kto_patch)

IF (lhook) CALL dr_hook('MODD_TEB_GREENROOF_N:TEB_GREENROOF_GOTO_PATCH',1,zhook_handle)
!
END SUBROUTINE teb_greenroof_goto_patch
!
SUBROUTINE teb_greenroof_init(YTEB_GREENROOF,KPATCH)
TYPE(teb_greenroof_t), INTENT(INOUT) :: yteb_greenroof
INTEGER, INTENT(IN) :: kpatch
INTEGER :: jp
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_TEB_GREENROOF_N:TEB_GREENROOF_INIT",0,zhook_handle)
 ALLOCATE(yteb_greenroof%ALP(kpatch))
 yteb_greenroof%CUR => yteb_greenroof%ALP(1)
DO jp=1,kpatch
  nullify(yteb_greenroof%ALP(jp)%XWR)
  nullify(yteb_greenroof%ALP(jp)%XTG)
  nullify(yteb_greenroof%ALP(jp)%XWG)
  nullify(yteb_greenroof%ALP(jp)%XWGI)
  nullify(yteb_greenroof%ALP(jp)%XRESA)
  nullify(yteb_greenroof%ALP(jp)%XAN)
  nullify(yteb_greenroof%ALP(jp)%XANDAY)
  nullify(yteb_greenroof%ALP(jp)%XANFM)
  nullify(yteb_greenroof%ALP(jp)%XLE)
  nullify(yteb_greenroof%ALP(jp)%XFAPARC)
  nullify(yteb_greenroof%ALP(jp)%XFAPIRC)
  nullify(yteb_greenroof%ALP(jp)%XLAI_EFFC)  
  nullify(yteb_greenroof%ALP(jp)%XMUS)     
  nullify(yteb_greenroof%ALP(jp)%XRESP_BIOMASS)
  nullify(yteb_greenroof%ALP(jp)%XBIOMASS)
  nullify(yteb_greenroof%ALP(jp)%XKSAT_ICE)
  nullify(yteb_greenroof%ALP(jp)%XPSNG)
  nullify(yteb_greenroof%ALP(jp)%XPSNV)
  nullify(yteb_greenroof%ALP(jp)%XPSNV_A)
  nullify(yteb_greenroof%ALP(jp)%XPSN)
  nullify(yteb_greenroof%ALP(jp)%XSNOWFREE_ALB)
  nullify(yteb_greenroof%ALP(jp)%XSNOWFREE_ALB_VEG)
  nullify(yteb_greenroof%ALP(jp)%XSNOWFREE_ALB_SOIL)
ENDDO 
IF (lhook) CALL dr_hook("MODD_TEB_GREENROOF_N:TEB_GREENROOF_INIT",1,zhook_handle)
END SUBROUTINE teb_greenroof_init


END MODULE modd_teb_greenroof_n