diag_teb_veg_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_veg_init_n(DK, DEK, DECK, DMK, KLU, KSNOW_LAYER)
              !     #####################
!
USE modd_surf_par, ONLY : xundef
!
USE mode_diag
!
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 yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
TYPE(diag_t), INTENT(INOUT) :: DK
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEK
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DECK
TYPE(diag_misc_isba_t), INTENT(INOUT) :: DMK
INTEGER, INTENT(IN) :: KLU
INTEGER, INTENT(IN) :: KSNOW_LAYER
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('DIAG_TEB_VEG_INIT',0,zhook_handle)
!
! Main diag
!
!general
ALLOCATE(dk%XRN   (klu))
ALLOCATE(dk%XH    (klu))
ALLOCATE(dk%XLE   (klu))
ALLOCATE(dk%XLEI  (klu))
ALLOCATE(dk%XGFLUX(klu))
ALLOCATE(dk%XEVAP (klu))
ALLOCATE(dk%XSUBL (klu))
!
dk%XRN   (:) = xundef
dk%XH    (:) = xundef
dk%XLE   (:) = xundef
dk%XLEI  (:) = xundef
dk%XGFLUX(:) = xundef
dk%XEVAP (:) = xundef
dk%XSUBL (:) = xundef
!
!
ALLOCATE(dk%XRI(klu))
!
dk%XRI(:) = xundef
!
!
ALLOCATE(dk%XCD (klu))
ALLOCATE(dk%XCH (klu))
ALLOCATE(dk%XZ0 (klu))
ALLOCATE(dk%XZ0H(klu))
!
dk%XCD (:) = xundef
dk%XCH (:) = xundef
dk%XZ0 (:) = xundef
dk%XZ0H(:) = xundef
!
!
ALLOCATE(dk%XZ0EFF(klu))
ALLOCATE(dk%XCDN  (klu))
ALLOCATE(dk%XHUG  (klu))
ALLOCATE(dk%XHU   (klu))
ALLOCATE(dk%XQS   (klu))
!
dk%XZ0EFF(:) = xundef
dk%XCDN  (:) = xundef
dk%XHUG  (:) = xundef
dk%XHU   (:) = xundef
dk%XQS   (:) = xundef
!
!
ALLOCATE(dk%XTS   (klu))
ALLOCATE(dk%XTSRAD(klu))
ALLOCATE(dk%XALBT (klu))
!
dk%XTS   (:) = xundef
dk%XTSRAD(:) = xundef
dk%XALBT (:) = xundef
!
!---------------------------
!
! Evaporation diag
!
! general
ALLOCATE(dek%XLEG (klu))
ALLOCATE(dek%XLEGI(klu))
ALLOCATE(dek%XLEV (klu))
ALLOCATE(dek%XLES (klu))
ALLOCATE(dek%XLESL(klu))
ALLOCATE(dek%XLER (klu))
ALLOCATE(dek%XLETR(klu))
!
dek%XLEG (:) = xundef
dek%XLEGI(:) = xundef
dek%XLEV (:) = xundef
dek%XLES (:) = xundef
dek%XLESL(:) = xundef
dek%XLER (:) = xundef
dek%XLETR(:) = xundef
!
!
ALLOCATE(dek%XDRAIN (klu))
ALLOCATE(dek%XRUNOFF(klu))
ALLOCATE(dek%XHORT  (klu))
ALLOCATE(dek%XQSB   (klu))
ALLOCATE(dek%XIRRIG_FLUX(klu))
!
dek%XDRAIN (:) = xundef
dek%XRUNOFF(:) = xundef
dek%XHORT  (:) = xundef
dek%XQSB   (:) = xundef
dek%XIRRIG_FLUX(:) = xundef
!
!
ALLOCATE(dek%XDRIP (klu))
ALLOCATE(dek%XRRVEG(klu))
!
dek%XDRIP (:) = xundef
dek%XRRVEG(:) = xundef
!
!
ALLOCATE(dek%XMELT    (klu))
ALLOCATE(dek%XMELTADV (klu))
ALLOCATE(dek%XRESTORE (klu))
ALLOCATE(dek%XSNDRIFT (klu))
ALLOCATE(dek%XSWNET_N (klu))
ALLOCATE(dek%XSWNET_NS(klu))
ALLOCATE(dek%XLWNET_N (klu))
!
dek%XMELT    (:) = xundef
dek%XMELTADV (:) = xundef
dek%XRESTORE (:) = xundef
dek%XSNDRIFT (:) = xundef
dek%XSWNET_N (:) = xundef
dek%XSWNET_NS(:) = xundef
dek%XLWNET_N (:) = xundef
!
!
ALLOCATE(dek%XLE_FLOOD (klu))
ALLOCATE(dek%XLEI_FLOOD(klu))
ALLOCATE(dek%XPFLOOD   (klu))
ALLOCATE(dek%XIFLOOD   (klu))
!
dek%XLE_FLOOD (:) = xundef
dek%XLEI_FLOOD(:) = xundef
dek%XPFLOOD   (:) = xundef
dek%XIFLOOD   (:) = xundef
!
!
ALLOCATE(dek%XGPP      (klu))
ALLOCATE(dek%XRESP_ECO (klu))
ALLOCATE(dek%XRESP_AUTO(klu))
!
dek%XGPP      (:) = xundef
dek%XRESP_ECO (:) = xundef
dek%XRESP_AUTO(:) = xundef
!
!
ALLOCATE(dek%XRN_SN_FR   (klu))
ALLOCATE(dek%XH_SN_FR    (klu))
ALLOCATE(dek%XLEI_SN_FR  (klu))
ALLOCATE(dek%XLE_SN_FR   (klu))
ALLOCATE(dek%XGFLUX_SN_FR(klu))
ALLOCATE(dek%XLEG_SN_FR  (klu))
ALLOCATE(dek%XLEGI_SN_FR (klu))
ALLOCATE(dek%XLEV_SN_FR  (klu))
ALLOCATE(dek%XLETR_SN_FR (klu))
ALLOCATE(dek%XUSTAR_SN_FR(klu))
ALLOCATE(dek%XLER_SN_FR  (klu))
!
dek%XRN_SN_FR   (:) = xundef
dek%XH_SN_FR    (:) = xundef
dek%XLEI_SN_FR  (:) = xundef
dek%XLE_SN_FR   (:) = xundef
dek%XGFLUX_SN_FR(:) = xundef
dek%XLEG_SN_FR  (:) = xundef
dek%XLEGI_SN_FR (:) = xundef
dek%XLEV_SN_FR  (:) = xundef
dek%XLETR_SN_FR (:) = xundef
dek%XUSTAR_SN_FR(:) = xundef
dek%XLER_SN_FR  (:) = xundef
!
!
ALLOCATE(deck%XDRAIN (klu))
ALLOCATE(deck%XRUNOFF(klu))
!
deck%XDRAIN (:) = 0.
deck%XRUNOFF(:) = 0.
!
!--------------------------------------
!
! Misc diag
!
ALLOCATE(dmk%XC1  (klu))
ALLOCATE(dmk%XC2  (klu))
ALLOCATE(dmk%XWGEQ(klu))
ALLOCATE(dmk%XCG  (klu))
ALLOCATE(dmk%XCT  (klu))
ALLOCATE(dmk%XRS  (klu))
ALLOCATE(dmk%XHV  (klu))
ALLOCATE(dmk%XGRNDFLUX (klu))
ALLOCATE(dmk%XSNOWTEMP (klu,ksnow_layer))
ALLOCATE(dmk%XSNOWHMASS(klu))
ALLOCATE(dmk%XSNOWLIQ  (klu,ksnow_layer))
ALLOCATE(dmk%XSNOWDZ   (klu,ksnow_layer))
ALLOCATE(dmk%XSRSFC    (klu))
ALLOCATE(dmk%XRRSFC    (klu))
ALLOCATE(dmk%XRNSNOW   (klu))
ALLOCATE(dmk%XHSNOW    (klu))
ALLOCATE(dmk%XGFLUXSNOW(klu))
ALLOCATE(dmk%XHPSNOW   (klu))
ALLOCATE(dmk%XUSTARSNOW(klu))
ALLOCATE(dmk%XCDSNOW   (klu))
ALLOCATE(dmk%XCHSNOW   (klu))
!
dmk%XC1        = xundef
dmk%XC2        = xundef
dmk%XWGEQ      = xundef
dmk%XCG        = xundef
dmk%XCT        = xundef
dmk%XRS        = xundef
dmk%XHV        = xundef
dmk%XGRNDFLUX  = xundef
dmk%XSNOWTEMP  = xundef
dmk%XSNOWHMASS = xundef
dmk%XSNOWLIQ   = xundef
dmk%XSNOWDZ    = xundef
dmk%XSRSFC     = xundef
dmk%XRRSFC     = xundef
dmk%XRNSNOW    = xundef
dmk%XHSNOW     = xundef
dmk%XGFLUXSNOW = xundef
dmk%XHPSNOW    = xundef
dmk%XUSTARSNOW = xundef
dmk%XCDSNOW    = xundef
dmk%XCHSNOW    = xundef
!
! LPROSNOW = .FALSE. for GARDEN GREENROOF
ALLOCATE(dmk%XSNOWDEND(0,0))
!
IF (lhook) CALL dr_hook('DIAG_TEB_VEG_INIT',1,zhook_handle)
!
END SUBROUTINE diag_teb_veg_init_n