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 (DGCT, DGMT, 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_cumul_teb_n, ONLY : diag_cumul_teb_t
USE modd_diag_misc_teb_n, ONLY : diag_misc_teb_t
USE modd_teb_option_n, ONLY : teb_options_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_cumul_teb_t), INTENT(INOUT) :: dgct
TYPE(diag_misc_teb_t), INTENT(INOUT) :: dgmt
TYPE(teb_options_t), INTENT(INOUT) :: top
!
      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(dgmt%CUR%XRUNOFF_ROOF,1)
!
 IF (top%LSOLAR_PANEL) THEN
    IF (dgmt%CUR%XTHER_PROD_BLD(ji) .NE. xundef) THEN
     dgct%CUR%XTHER_PROD_BLDC(ji)     =  dgct%CUR%XTHER_PROD_BLDC(ji)     + dgmt%CUR%XTHER_PROD_BLD(ji)  * ptstep
    ENDIF
    !
    IF (dgmt%CUR%XPHOT_PROD_BLD(ji) .NE. xundef) THEN
     dgct%CUR%XPHOT_PROD_BLDC(ji)     =  dgct%CUR%XPHOT_PROD_BLDC(ji)     + dgmt%CUR%XPHOT_PROD_BLD(ji)  * ptstep
    ENDIF
 END IF

 IF (top%CBEM == 'BEM') THEN
    IF (dgmt%CUR%XHVAC_COOL(ji) .NE. xundef) THEN
     dgct%CUR%XHVACC_COOL(ji)     =  dgct%CUR%XHVACC_COOL(ji)       + dgmt%CUR%XHVAC_COOL(ji)        * ptstep
    ENDIF
    !
    IF (dgmt%CUR%XHVAC_HEAT(ji) .NE. xundef) THEN
     dgct%CUR%XHVACC_HEAT(ji)     =  dgct%CUR%XHVACC_HEAT(ji)       + dgmt%CUR%XHVAC_HEAT(ji)        * ptstep
    ENDIF
 ENDIF
 !
 IF (dgmt%CUR%XRUNOFF_TOWN(ji) .NE. xundef) THEN
  dgct%CUR%XRUNOFFC_TOWN(ji)      =  dgct%CUR%XRUNOFFC_TOWN(ji)     + dgmt%CUR%XRUNOFF_TOWN(ji)      * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XRUNOFF_GARDEN(ji) .NE. xundef) THEN
  dgct%CUR%XRUNOFFC_GARDEN(ji)    =  dgct%CUR%XRUNOFFC_GARDEN(ji)   + dgmt%CUR%XRUNOFF_GARDEN(ji)    * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XRUNOFF_ROAD(ji) .NE. xundef) THEN
  dgct%CUR%XRUNOFFC_ROAD(ji)      =  dgct%CUR%XRUNOFFC_ROAD(ji)     + dgmt%CUR%XRUNOFF_ROAD(ji)      * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XRUNOFF_ROOF(ji) .NE. xundef) THEN 
  dgct%CUR%XRUNOFFC_ROOF(ji)      =  dgct%CUR%XRUNOFFC_ROOF(ji)     + dgmt%CUR%XRUNOFF_ROOF(ji)      * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XRUNOFF_STRLROOF(ji) .NE. xundef) THEN
  dgct%CUR%XRUNOFFC_STRLROOF(ji)  =  dgct%CUR%XRUNOFFC_STRLROOF(ji) + dgmt%CUR%XRUNOFF_STRLROOF(ji)  * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XDRAIN_GARDEN(ji) .NE. xundef) THEN
   dgct%CUR%XDRAINC_GARDEN(ji)    =  dgct%CUR%XDRAINC_GARDEN(ji)    + dgmt%CUR%XDRAIN_GARDEN(ji)     * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XIRRIG_GARDEN(ji) .NE. xundef) THEN
   dgct%CUR%XIRRIGC_GARDEN(ji)    =  dgct%CUR%XIRRIGC_GARDEN(ji)    + dgmt%CUR%XIRRIG_GARDEN(ji)     * ptstep
 ENDIF
 !
 IF (dgmt%CUR%XIRRIG_ROAD(ji) .NE. xundef) THEN
   dgct%CUR%XIRRIGC_ROAD(ji)      =  dgct%CUR%XIRRIGC_ROAD(ji)      + dgmt%CUR%XIRRIG_ROAD(ji)       * ptstep
 ENDIF
 !
 IF (top%LGREENROOF) THEN 
    IF (dgmt%CUR%XRUNOFF_GREENROOF(ji) .NE. xundef) THEN
     dgct%CUR%XRUNOFFC_GREENROOF(ji) =  dgct%CUR%XRUNOFFC_GREENROOF(ji)+ dgmt%CUR%XRUNOFF_GREENROOF(ji) * ptstep
    ENDIF
    !
    IF (dgmt%CUR%XDRAIN_GREENROOF(ji) .NE. xundef) THEN
     dgct%CUR%XDRAINC_GREENROOF(ji)  =  dgct%CUR%XDRAINC_GREENROOF(ji) + dgmt%CUR%XDRAIN_GREENROOF(ji)  * ptstep
    ENDIF
    !
    IF (dgmt%CUR%XIRRIG_GREENROOF(ji) .NE. xundef) THEN
     dgct%CUR%XIRRIGC_GREENROOF(ji)  =  dgct%CUR%XIRRIGC_GREENROOF(ji) + dgmt%CUR%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