modd_teb_gardenn.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_garden_n
!##################
!
!!****  *MODD_TEB_GARDEN - declaration of packed surface parameters for ISBA scheme
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      A. Lemonsu   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       01/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_garden_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) 
!
!
!-------------------------------------------------------------------------------
!
! - 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_garden_1p_t
!
TYPE teb_garden_t
  !
  TYPE(teb_garden_1p_t), POINTER :: ALP(:) => NULL()
  TYPE(teb_garden_1p_t), POINTER :: CUR => NULL()
  !
END TYPE teb_garden_t
!


 CONTAINS


!


!

SUBROUTINE teb_garden_goto_patch(YTEB_GARDEN,KTO_PATCH)
TYPE(teb_garden_t), INTENT(INOUT) :: yteb_garden
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_GARDEN_N:TEB_GARDEN_GOTO_PATCH',0,zhook_handle)

yteb_garden%CUR => yteb_garden%ALP(kto_patch)

IF (lhook) CALL dr_hook('MODD_TEB_GARDEN_N:TEB_GARDEN_GOTO_PATCH',1,zhook_handle)
!
END SUBROUTINE teb_garden_goto_patch

SUBROUTINE teb_garden_init(YTEB_GARDEN,KPATCH)
TYPE(teb_garden_t), INTENT(INOUT) :: yteb_garden
INTEGER, INTENT(IN) :: kpatch
INTEGER :: jp
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_TEB_GARDEN_N:TEB_GARDEN_INIT",0,zhook_handle)
 ALLOCATE(yteb_garden%ALP(kpatch))
 yteb_garden%CUR => yteb_garden%ALP(1)
DO jp=1,kpatch
  nullify(yteb_garden%ALP(jp)%XWR)
  nullify(yteb_garden%ALP(jp)%XTG)
  nullify(yteb_garden%ALP(jp)%XWG)
  nullify(yteb_garden%ALP(jp)%XWGI)
  nullify(yteb_garden%ALP(jp)%XRESA)
  nullify(yteb_garden%ALP(jp)%XAN)
  nullify(yteb_garden%ALP(jp)%XANDAY)
  nullify(yteb_garden%ALP(jp)%XANFM)
  nullify(yteb_garden%ALP(jp)%XLE)
  nullify(yteb_garden%ALP(jp)%XFAPARC)
  nullify(yteb_garden%ALP(jp)%XFAPIRC)
  nullify(yteb_garden%ALP(jp)%XLAI_EFFC)  
  nullify(yteb_garden%ALP(jp)%XMUS)    
  nullify(yteb_garden%ALP(jp)%XRESP_BIOMASS)
  nullify(yteb_garden%ALP(jp)%XBIOMASS)
  nullify(yteb_garden%ALP(jp)%XPSNG)
  nullify(yteb_garden%ALP(jp)%XPSNV)
  nullify(yteb_garden%ALP(jp)%XPSNV_A)
  nullify(yteb_garden%ALP(jp)%XPSN)
  nullify(yteb_garden%ALP(jp)%XSNOWFREE_ALB)
  nullify(yteb_garden%ALP(jp)%XSNOWFREE_ALB_VEG)
  nullify(yteb_garden%ALP(jp)%XSNOWFREE_ALB_SOIL)
ENDDO 
IF (lhook) CALL dr_hook("MODD_TEB_GARDEN_N:TEB_GARDEN_INIT",1,zhook_handle)
END SUBROUTINE teb_garden_init



END MODULE modd_teb_garden_n