cumul_diag_tebn.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 cumul_diag_teb_n (DMTC, DMT, GDDEC, GDDE, GRDEC, GRDE, TOP, PTSTEP)
!##################################
!
!
!!****  *CUMUL_DIAG_TEB_n*  
!!
!!    PURPOSE
!!    -------
!      Cumulates some diagnostics for TEB
!     
!!**  METHOD
!!    ------
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------ 
!!
!!      
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      C. de Munck       * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    02/2013
!!                  08/2013 (V. Masson) adds solar panels
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_diag_misc_teb_n, ONLY : diag_misc_teb_t
USE modd_teb_option_n, ONLY : teb_options_t
!
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
!
USE modd_surf_par,        ONLY :  xundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
TYPE(diag_misc_teb_t), INTENT(INOUT) :: DMTC
TYPE(diag_misc_teb_t), INTENT(INOUT) :: DMT
TYPE(teb_options_t), INTENT(INOUT) :: TOP
!
TYPE(diag_evap_isba_t), INTENT(INOUT) :: GDDE
TYPE(diag_evap_isba_t), INTENT(INOUT) :: GDDEC
TYPE(diag_evap_isba_t), INTENT(INOUT) :: GRDE
TYPE(diag_evap_isba_t), INTENT(INOUT) :: GRDEC
!
REAL,               INTENT(IN) :: PTSTEP            ! time step
!
!*      0.2    declarations of local variables
!
INTEGER :: JI 

REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
!       0.     Initialization
!              --------------
IF (lhook) CALL dr_hook('CUMUL_DIAG_TEB_N',0,zhook_handle)
!
!       1.     Time-cumulated diagnostics for TEB
!              ----------------------------------
!
DO ji=1,SIZE(dmt%XRUNOFF_ROOF,1)
!
 IF (top%LSOLAR_PANEL) THEN
    IF (dmt%XTHER_PROD_BLD(ji) .NE. xundef) THEN
     dmtc%XTHER_PROD_BLD(ji)     =  dmtc%XTHER_PROD_BLD(ji)     + dmt%XTHER_PROD_BLD(ji)  * ptstep
    ENDIF
    !
    IF (dmt%XPHOT_PROD_BLD(ji) .NE. xundef) THEN
     dmtc%XPHOT_PROD_BLD(ji)     =  dmtc%XPHOT_PROD_BLD(ji)     + dmt%XPHOT_PROD_BLD(ji)  * ptstep
    ENDIF
 END IF

 IF (top%CBEM == 'BEM') THEN
    IF (dmt%XHVAC_COOL(ji) .NE. xundef) THEN
     dmtc%XHVAC_COOL(ji)     =  dmtc%XHVAC_COOL(ji)       + dmt%XHVAC_COOL(ji)        * ptstep
    ENDIF
    !
    IF (dmt%XHVAC_HEAT(ji) .NE. xundef) THEN
     dmtc%XHVAC_HEAT(ji)     =  dmtc%XHVAC_HEAT(ji)       + dmt%XHVAC_HEAT(ji)        * ptstep
    ENDIF
 ENDIF
 !
 IF (dmt%XRUNOFF_TOWN(ji) .NE. xundef) THEN
  dmtc%XRUNOFF_TOWN(ji)      =  dmtc%XRUNOFF_TOWN(ji)     + dmt%XRUNOFF_TOWN(ji)      * ptstep
 ENDIF
 !
 IF (dmt%XRUNOFF_ROAD(ji) .NE. xundef) THEN
  dmtc%XRUNOFF_ROAD(ji)      =  dmtc%XRUNOFF_ROAD(ji)     + dmt%XRUNOFF_ROAD(ji)      * ptstep
 ENDIF
 !
 IF (dmt%XRUNOFF_ROOF(ji) .NE. xundef) THEN 
  dmtc%XRUNOFF_ROOF(ji)      =  dmtc%XRUNOFF_ROOF(ji)     + dmt%XRUNOFF_ROOF(ji)      * ptstep
 ENDIF
 !
 IF (dmt%XRUNOFF_STRLROOF(ji) .NE. xundef) THEN
  dmtc%XRUNOFF_STRLROOF(ji)  =  dmtc%XRUNOFF_STRLROOF(ji) + dmt%XRUNOFF_STRLROOF(ji)  * ptstep
 ENDIF
 !
 IF (dmt%XIRRIG_ROAD(ji) .NE. xundef) THEN
   dmtc%XIRRIG_ROAD(ji)      =  dmtc%XIRRIG_ROAD(ji)      + dmt%XIRRIG_ROAD(ji)       * ptstep
 ENDIF
 !
 IF (top%LGARDEN) THEN
   !
   IF (gdde%XRUNOFF(ji) .NE. xundef) THEN
     gddec%XRUNOFF(ji)    =  gddec%XRUNOFF(ji)   + gdde%XRUNOFF(ji)    * ptstep
   ENDIF 
   !
   IF (gdde%XDRAIN(ji) .NE. xundef) THEN
     gddec%XDRAIN(ji)    =  gddec%XDRAIN(ji)    + gdde%XDRAIN(ji)     * ptstep
   ENDIF
   !
   IF (dmt%XIRRIG_GARDEN(ji) .NE. xundef) THEN
     dmtc%XIRRIG_GARDEN(ji)    =  dmtc%XIRRIG_GARDEN(ji)    + dmt%XIRRIG_GARDEN(ji)     * ptstep
   ENDIF
   !
 ENDIF
 ! 
 IF (top%LGREENROOF) THEN 

    IF (grde%XRUNOFF(ji) .NE. xundef) THEN
     grdec%XRUNOFF(ji) =  grdec%XRUNOFF(ji)+ grde%XRUNOFF(ji) * ptstep
    ENDIF
    !
    IF (grde%XDRAIN(ji) .NE. xundef) THEN
     grdec%XDRAIN(ji)  =  grdec%XDRAIN(ji) + grde%XDRAIN(ji)  * ptstep
    ENDIF
    !
    IF (dmt%XIRRIG_GREENROOF(ji) .NE. xundef) THEN
     dmtc%XIRRIG_GREENROOF(ji)  =  dmtc%XIRRIG_GREENROOF(ji) + dmt%XIRRIG_GREENROOF(ji)  * ptstep
    ENDIF
 ENDIF
 !
ENDDO
!
IF (lhook) CALL dr_hook('CUMUL_DIAG_TEB_N',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE cumul_diag_teb_n