diag_teb_greenroof_initn.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 diag_teb_greenroof_init_n (DGMTO, DGGR, TGRO, TVG, &
                                            hprogram,klu,ksw)
!     #####################
!
!!****  *DIAG_TEB_GREENROOF_INIT_n* - routine to initialize TEB-ISBA diagnostic variables
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!    Based on diag_teb_garden_initn
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original                  02/2003 
!!      modified                  11/2003  by P. LeMoigne: surface cumulated energy budget
!!      modified                  10/2004  by P. LeMoigne: surface miscellaneous fields
!!      B. Decharme               2008     New diag for water budget and allow to reset
!!                                         cumulatives variables at the beginning of a run
!       C. de Munck & A. Lemonsu  09/2011  Greenroofs
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE modd_diag_misc_teb_option_n, ONLY : diag_misc_teb_options_t
USE modd_diag_teb_greenroof_n, ONLY : diag_teb_greenroof_t
USE modd_teb_greenroof_option_n, ONLY : teb_greenroof_options_t
USE modd_teb_veg_n, ONLY : teb_veg_options_t
!
USE modd_surf_par,             ONLY : xundef
USE modd_type_date_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(diag_misc_teb_options_t), INTENT(INOUT) :: dgmto
TYPE(diag_teb_greenroof_t), INTENT(INOUT) :: dggr
TYPE(teb_greenroof_options_t), INTENT(INOUT) :: tgro
TYPE(teb_veg_options_t), INTENT(INOUT) :: tvg
!
INTEGER, INTENT(IN)         :: klu       ! size of arrays
INTEGER, INTENT(IN)         :: ksw       ! spectral bands
 CHARACTER(LEN=6), INTENT(IN):: hprogram  ! program calling
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER           :: iresp          ! IRESP  : return-code if a problem appears
 CHARACTER(LEN=12) :: yrec           ! Name of the article to be read
REAL(KIND=JPRB) :: zhook_handle
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('DIAG_TEB_GREENROOF_INIT_N',0,zhook_handle)
dggr%XCG         = xundef
dggr%XC1         = xundef
dggr%XC2         = xundef
dggr%XWGEQ       = xundef
dggr%XCT         = xundef
dggr%XRS         = xundef
dggr%XHU         = xundef
dggr%XHUG        = xundef
dggr%XHV         = xundef
dggr%XRESTORE    = xundef
dggr%XRI         = xundef
dggr%XUSTAR      = xundef
dggr%XRN         = xundef
dggr%XH          = xundef
dggr%XGFLUX      = xundef
dggr%XSNOWTEMP   = xundef
dggr%XSNOWLIQ    = xundef
dggr%XSNOWDZ     = xundef
dggr%XSNOWHMASS  = xundef
dggr%XMELTADV    = xundef
IF (tvg%CPHOTO/='NON') THEN
  dggr%XIACAN    = xundef
END IF
dggr%XCD         = xundef
dggr%XCDN        = xundef
dggr%XCH         = xundef
dggr%XQS         = xundef
dggr%XLEI        = xundef
dggr%XLEG        = xundef
dggr%XLEGI       = xundef
dggr%XLEV        = xundef
dggr%XLES        = xundef
dggr%XLER        = xundef
dggr%XLETR       = xundef
dggr%XEVAP       = xundef
dggr%XDRAIN      = xundef
dggr%XRUNOFF     = xundef
dggr%XHORT       = xundef
dggr%XDRIP       = xundef
dggr%XRRVEG      = xundef
dggr%XMELT       = xundef
dggr%XALBT       = xundef
dggr%XEMIST      = xundef
!
!* surface energy budget
!
!IF (LSURF_BUDGET) THEN
  !
  ALLOCATE(dggr%XSWD      (klu))
  ALLOCATE(dggr%XSWU      (klu))
  ALLOCATE(dggr%XSWBD     (klu,ksw))
  ALLOCATE(dggr%XSWBU     (klu,ksw))
  ALLOCATE(dggr%XLWD      (klu))
  ALLOCATE(dggr%XLWU      (klu))
  ALLOCATE(dggr%XFMU      (klu))
  ALLOCATE(dggr%XFMV      (klu))
  !
  dggr%XSWD     = xundef
  dggr%XSWU     = xundef
  dggr%XSWBD    = xundef
  dggr%XSWBU    = xundef
  dggr%XLWD     = xundef
  dggr%XLWU     = xundef
  dggr%XFMU     = xundef
  dggr%XFMV     = xundef
  !
!END IF
!
!* surface temperature and parameters at 2m
!
ALLOCATE(dggr%XTS    (klu))
dggr%XTS     = xundef
ALLOCATE(dggr%XTSRAD (klu))
dggr%XTSRAD  = xundef
!
!* miscellaneous surface fields
!
IF (dgmto%LSURF_MISC_BUDGET) THEN
  !
  ALLOCATE(dggr%XSWI    (klu,tgro%NLAYER_GR))
  ALLOCATE(dggr%XTSWI   (klu,tgro%NLAYER_GR))
  ALLOCATE(dggr%XTWSNOW (klu))
  ALLOCATE(dggr%XTDSNOW (klu))
  dggr%XSWI     = xundef
  dggr%XTSWI    = xundef
  dggr%XTWSNOW  = xundef
  dggr%XTDSNOW  = xundef
ENDIF


  ALLOCATE(dggr%XALBT   (klu))
  ALLOCATE(dggr%XGPP    (klu))
  ALLOCATE(dggr%XRESP_AUTO  (klu))
  ALLOCATE(dggr%XRESP_ECO   (klu))
  !
  dggr%XALBT    = xundef
  dggr%XGPP     = xundef
  dggr%XRESP_AUTO   = xundef
  dggr%XRESP_ECO    = xundef  
  !
!END IF
!
!* transfer coefficients
!
!IF (LCOEF) THEN
  !
  ALLOCATE(dggr%XCE            (klu))
  ALLOCATE(dggr%XZ0_WITH_SNOW  (klu))
  ALLOCATE(dggr%XZ0H_WITH_SNOW (klu))
  ALLOCATE(dggr%XZ0EFF         (klu))
  !
  dggr%XCE            = xundef
  dggr%XZ0_WITH_SNOW  = xundef
  dggr%XZ0H_WITH_SNOW = xundef
  dggr%XZ0EFF         = xundef
!END IF
!
!
!* surface humidity
!
!IF (LSURF_VARS) THEN
  ALLOCATE(dggr%XQS            (klu))
  !
  dggr%XQS            = xundef
!END IF
!
IF (lhook) CALL dr_hook('DIAG_TEB_GREENROOF_INIT_N',1,zhook_handle)
!
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE diag_teb_greenroof_init_n