mode_diag.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 mode_diag
!   #####################
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modd_surf_par, ONLY : xundef
USE modd_diag_n, ONLY : diag_t, diag_options_t
USE modd_diag_evap_isba_n, ONLY : diag_evap_isba_t
!
IMPLICIT NONE
!
CONTAINS 
!
SUBROUTINE alloc_surf_bud(DA,KLUA,KLUAC,KSWA)
!
TYPE(diag_t), INTENT(INOUT) :: DA
INTEGER, INTENT(IN) :: KLUA
INTEGER, INTENT(IN) :: KLUAC
INTEGER, INTENT(IN) :: KSWA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_SURF_BUD',0,zhook_handle)
!
ALLOCATE(da%XRN       (kluac))
ALLOCATE(da%XH        (kluac))
ALLOCATE(da%XLE       (kluac))
ALLOCATE(da%XLEI      (kluac))
ALLOCATE(da%XGFLUX    (kluac))
!
ALLOCATE(da%XSWD      (kluac))
ALLOCATE(da%XSWU      (kluac))
ALLOCATE(da%XLWD      (kluac))
ALLOCATE(da%XLWU      (kluac))
ALLOCATE(da%XFMU      (kluac))
ALLOCATE(da%XFMV      (kluac))
!
IF (kluac>0) THEN
  da%XRN      = xundef
  da%XH       = xundef
  da%XLE      = xundef
  da%XLEI     = xundef
  da%XGFLUX   = xundef

  da%XSWD     = xundef
  da%XSWU     = xundef
  da%XLWD     = xundef
  da%XLWU     = xundef
  da%XFMU     = xundef
  da%XFMV     = xundef
ENDIF
!
ALLOCATE(da%XSWBD     (klua,kswa))
ALLOCATE(da%XSWBU     (klua,kswa))
!
IF (klua>0) THEN
  da%XSWBD    = xundef
  da%XSWBU    = xundef
ENDIF
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_SURF_BUD',1,zhook_handle)
!
END SUBROUTINE alloc_surf_bud
!
SUBROUTINE alloc_n2m_bud(DA,KLUA)
!
TYPE(diag_t), INTENT(INOUT) :: DA
INTEGER, INTENT(IN) :: KLUA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_N2M_BUD',0,zhook_handle)
!
ALLOCATE(da%XRI           (klua))
ALLOCATE(da%XT2M          (klua))
ALLOCATE(da%XT2M_MIN      (klua))
ALLOCATE(da%XT2M_MAX      (klua))
ALLOCATE(da%XQ2M          (klua))
ALLOCATE(da%XHU2M         (klua))
ALLOCATE(da%XHU2M_MIN     (klua))
ALLOCATE(da%XHU2M_MAX     (klua))
ALLOCATE(da%XZON10M       (klua))
ALLOCATE(da%XMER10M       (klua))
ALLOCATE(da%XWIND10M      (klua))
ALLOCATE(da%XWIND10M_MAX  (klua))
!
IF (klua>0) THEN
  da%XRI      = xundef
  da%XT2M     = xundef
  da%XT2M_MIN = xundef
  da%XT2M_MAX = -xundef
  da%XQ2M     = xundef
  da%XHU2M    = xundef
  da%XHU2M_MIN= xundef
  da%XHU2M_MAX= -xundef
  da%XZON10M  = xundef
  da%XMER10M  = xundef
  da%XWIND10M = xundef
  da%XWIND10M_MAX = 0.0
ENDIF
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_N2M_BUD',1,zhook_handle)
!
END SUBROUTINE alloc_n2m_bud
!
SUBROUTINE alloc_coef_bud(DA,KLUA)
!
TYPE(diag_t), INTENT(INOUT) :: DA
INTEGER, INTENT(IN) :: KLUA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_COEF_BUD',0,zhook_handle)
!
ALLOCATE(da%XCD   (klua))
ALLOCATE(da%XCH   (klua))
ALLOCATE(da%XCE   (klua))
ALLOCATE(da%XZ0   (klua))
ALLOCATE(da%XZ0H  (klua))
!
IF (klua>0) THEN
  da%XCD      = xundef
  da%XCH      = xundef
  da%XCE      = xundef
  da%XZ0      = xundef
  da%XZ0H     = xundef
ENDIF

IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_COEF_BUD',1,zhook_handle)
!
END SUBROUTINE alloc_coef_bud
!
SUBROUTINE alloc_bud(DGO,DA,KLU,KSW)
!
TYPE(diag_options_t), INTENT(IN) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DA
INTEGER, INTENT(IN) :: KLU
INTEGER, INTENT(IN) :: KSW
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_BUD',0,zhook_handle)
!
IF (dgo%LSURF_BUDGET .OR. dgo%LSURF_BUDGETC) THEN
  CALL alloc_surf_bud(da,klu,klu,ksw)
ELSE
  CALL alloc_surf_bud(da,0,0,0)
END IF
!
!* parameters at 2m
!
IF (dgo%N2M>=1) THEN
  CALL alloc_n2m_bud(da,klu)
ELSE
  CALL alloc_n2m_bud(da,0)
END IF
!
!* transfer coefficients
!
IF (dgo%LCOEF) THEN
  CALL alloc_coef_bud(da,klu)
ELSE
  CALL alloc_coef_bud(da,0)
END IF
!
!* surface humidity
!
IF (dgo%LSURF_VARS) THEN
  ALLOCATE(da%XQS     (klu))
  da%XQS      = xundef
ELSE
  ALLOCATE(da%XQS     (0))
END IF
!
ALLOCATE(da%XTS(klu))
da%XTS = xundef
!
IF (lhook) CALL dr_hook('MODE_DIAG:ALLOC_BUD',1,zhook_handle)
!
END SUBROUTINE alloc_bud
!
SUBROUTINE init_bud(DGO,DA,DAC,PVAL)
!
TYPE(diag_options_t), INTENT(IN) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DA
TYPE(diag_t), INTENT(INOUT) :: DAC
REAL, INTENT(IN) :: PVAL
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_BUD',0,zhook_handle)
!
IF (dgo%LSURF_BUDGET) THEN
  CALL init_surf_bud(da,pval)
  da%XEVAP = pval
  da%XSUBL = pval
ENDIF
IF (dgo%LSURF_BUDGETC) THEN
  CALL init_surf_bud(dac,pval)
  dac%XEVAP = pval
  dac%XSUBL = pval  
ENDIF
IF (dgo%N2M>=1) CALL init_n2m_bud(da,pval)
IF (dgo%LCOEF) CALL init_coef_bud(da,pval)
da%XTS = pval
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_BUD',1,zhook_handle)
!
END SUBROUTINE init_bud
!
SUBROUTINE init_surf_bud(DA,PVAL)
!
TYPE(diag_t), INTENT(INOUT) :: DA
REAL, INTENT(IN) :: PVAL
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_SURF_BUD',0,zhook_handle)
!
da%XRN      = pval
da%XH       = pval
da%XLE      = pval
da%XLEI     = pval
da%XGFLUX   = pval
!
da%XSWD     = pval
da%XSWU     = pval
da%XLWD     = pval
da%XLWU     = pval
da%XFMU     = pval
da%XFMV     = pval
!
IF (SIZE(da%XSWBD)>0) THEN
  da%XSWBD    = pval
  da%XSWBU    = pval
ENDIF
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_SURF_BUD',1,zhook_handle)
!
END SUBROUTINE init_surf_bud
!
SUBROUTINE init_n2m_bud(DA,PVAL)
!
TYPE(diag_t), INTENT(INOUT) :: DA
REAL, INTENT(IN) :: PVAL
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_N2M_BUD',0,zhook_handle)
!
da%XRI      = pval
da%XT2M     = pval
da%XT2M_MIN = pval
da%XT2M_MAX = pval
da%XQ2M     = pval
da%XHU2M    = pval
da%XHU2M_MIN= pval
da%XHU2M_MAX= pval
da%XZON10M  = pval
da%XMER10M  = pval
da%XWIND10M = pval
da%XWIND10M_MAX = pval
da%XSFCO2       = pval
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_N2M_BUD',1,zhook_handle)
!
END SUBROUTINE init_n2m_bud
!
SUBROUTINE init_coef_bud(DA,PVAL)
!
TYPE(diag_t), INTENT(INOUT) :: DA
REAL, INTENT(IN) :: PVAL
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_COEF_BUD',0,zhook_handle)
!
da%XCD      = pval
da%XCH      = pval
da%XCE      = pval
da%XZ0      = pval
da%XZ0H     = pval
da%XZ0EFF   = pval
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_COEF_BUD',1,zhook_handle)
!
END SUBROUTINE init_coef_bud
!
SUBROUTINE diag_evap(DGO, DA, DAC, HPROGRAM, DAUP, DAUPC, KMASK )
!
IMPLICIT NONE
!
TYPE(diag_options_t), INTENT(IN) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DA
TYPE(diag_t), INTENT(INOUT) :: DAC
 CHARACTER(LEN=6),   INTENT(IN)  :: HPROGRAM ! program calling surf. schemes
TYPE(diag_t), INTENT(INOUT) :: DAUP
TYPE(diag_t), INTENT(INOUT) :: DAUPC
INTEGER, DIMENSION(:), INTENT(IN) :: KMASK
!
INTEGER :: JJ, ISIZE
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:DIAG_EVAP',0,zhook_handle)
!
isize = SIZE(kmask)
!
 CALL diag_cumul(dgo, da, dac, hprogram, daup, daupc, kmask)
!
IF (dgo%LSURF_BUDGET) THEN
  DO jj=1,isize
    daup%XEVAP    (kmask(jj))  = da%XEVAP     (jj)
    daup%XSUBL    (kmask(jj))  = da%XSUBL     (jj)
  ENDDO
END IF
!
IF (lhook) CALL dr_hook('MODE_DIAG:DIAG_EVAP',1,zhook_handle)
!
END SUBROUTINE diag_evap
!
SUBROUTINE diag_cumul(DGO, DA, DAC, HPROGRAM, DAUP, DAUPC, KMASK )
!
IMPLICIT NONE
!
TYPE(diag_options_t), INTENT(IN) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DA
TYPE(diag_t), INTENT(INOUT) :: DAC
 CHARACTER(LEN=6),   INTENT(IN)  :: HPROGRAM ! program calling surf. schemes
TYPE(diag_t), INTENT(INOUT) :: DAUP
TYPE(diag_t), INTENT(INOUT) :: DAUPC
INTEGER, DIMENSION(:), INTENT(IN) :: KMASK
!
INTEGER :: JJ, ISIZE
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:DIAG_CUMUL',0,zhook_handle)
!
isize = SIZE(kmask)
!
 CALL diag(dgo, da, hprogram, daup, kmask)
!
IF (dgo%LSURF_BUDGETC) THEN
  DO jj=1,isize
    daupc%XRN      (kmask(jj))  = dac%XRN       (jj)
    daupc%XH       (kmask(jj))  = dac%XH        (jj)
    daupc%XLE      (kmask(jj))  = dac%XLE       (jj)
    daupc%XLEI     (kmask(jj))  = dac%XLEI      (jj)
    daupc%XGFLUX   (kmask(jj))  = dac%XGFLUX    (jj)
    daupc%XEVAP    (kmask(jj))  = dac%XEVAP     (jj)
    daupc%XSUBL    (kmask(jj))  = dac%XSUBL     (jj)
    daupc%XSWD     (kmask(jj))  = dac%XSWD      (jj)
    daupc%XSWU     (kmask(jj))  = dac%XSWU      (jj)
    daupc%XLWD     (kmask(jj))  = dac%XLWD      (jj)
    daupc%XLWU     (kmask(jj))  = dac%XLWU      (jj)
    daupc%XFMU     (kmask(jj))  = dac%XFMU      (jj)
    daupc%XFMV     (kmask(jj))  = dac%XFMV      (jj)
  ENDDO        
END IF
!
IF (lhook) CALL dr_hook('MODE_DIAG:DIAG_CUMUL',1,zhook_handle)
!
END SUBROUTINE diag_cumul
!
SUBROUTINE diag(DGO, DA, HPROGRAM, DAUP, KMASK )
!
IMPLICIT NONE
!
TYPE(diag_options_t), INTENT(IN) :: DGO
TYPE(diag_t), INTENT(INOUT) :: DA
 CHARACTER(LEN=6),   INTENT(IN)  :: HPROGRAM ! program calling surf. schemes
TYPE(diag_t), INTENT(INOUT) :: DAUP
INTEGER, DIMENSION(:), INTENT(IN) :: KMASK
!
INTEGER :: JJ, ISIZE
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:DIAG',0,zhook_handle)
!
isize = SIZE(kmask)
!
IF (dgo%LSURF_BUDGET) THEN
  DO jj=1,isize
    daup%XRN      (kmask(jj))  = da%XRN       (jj)
    daup%XH       (kmask(jj))  = da%XH        (jj)
    daup%XLE      (kmask(jj))  = da%XLE       (jj)
    daup%XLEI     (kmask(jj))  = da%XLEI      (jj)
    daup%XGFLUX   (kmask(jj))  = da%XGFLUX    (jj)
    daup%XSWD     (kmask(jj))  = da%XSWD      (jj)
    daup%XSWU     (kmask(jj))  = da%XSWU      (jj)
    daup%XLWD     (kmask(jj))  = da%XLWD      (jj)
    daup%XLWU     (kmask(jj))  = da%XLWU      (jj)
    daup%XFMU     (kmask(jj))  = da%XFMU      (jj)
    daup%XFMV     (kmask(jj))  = da%XFMV      (jj)
    daup%XSWBD    (kmask(jj),:) = da%XSWBD     (jj,:)
    daup%XSWBU    (kmask(jj),:) = da%XSWBU     (jj,:)
  ENDDO
END IF
!
IF (dgo%N2M>=1 .OR. dgo%LSURF_BUDGET) THEN
  DO jj=1,isize
    daup%XTS       (kmask(jj))  = da%XTS      (jj)
  ENDDO        
ENDIF
!
IF (dgo%N2M>=1) THEN
  DO jj=1,isize
    daup%XRI      (kmask(jj))  = da%XRI       (jj)
    daup%XT2M     (kmask(jj))  = da%XT2M      (jj)
    daup%XT2M_MIN (kmask(jj))  = da%XT2M_MIN  (jj)
    daup%XT2M_MAX (kmask(jj))  = da%XT2M_MAX  (jj)
    daup%XQ2M     (kmask(jj))  = da%XQ2M      (jj)
    daup%XHU2M    (kmask(jj))  = da%XHU2M     (jj)
    daup%XHU2M_MIN(kmask(jj))  = da%XHU2M_MIN (jj)
    daup%XHU2M_MAX(kmask(jj))  = da%XHU2M_MAX (jj)
    daup%XZON10M  (kmask(jj))  = da%XZON10M   (jj)
    daup%XMER10M  (kmask(jj))  = da%XMER10M   (jj)
    daup%XWIND10M (kmask(jj))  = da%XWIND10M   (jj)
    daup%XWIND10M_MAX (kmask(jj))  = da%XWIND10M_MAX   (jj)
  ENDDO        
END IF 
!
IF (dgo%LCOEF) THEN
  DO jj=1,isize
    daup%XCD      (kmask(jj))  = da%XCD       (jj)
    daup%XCH      (kmask(jj))  = da%XCH       (jj)
    daup%XCE      (kmask(jj))  = da%XCE       (jj)
    daup%XZ0      (kmask(jj))  = da%XZ0       (jj)
    daup%XZ0H     (kmask(jj))  = da%XZ0H      (jj)
  ENDDO        
END IF
!
IF (dgo%LSURF_VARS) THEN
  DO jj=1,isize
    daup%XQS      (kmask(jj))  = da%XQS       (jj)
  ENDDO          
ENDIF
!
IF (lhook) CALL dr_hook('MODE_DIAG:DIAG',1,zhook_handle)
!
END SUBROUTINE diag
!
SUBROUTINE init_evap_bud(DEA)
!
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_EVAP_BUD',0,zhook_handle)
!
dea%XLEG        = 0.0
dea%XLEGI       = 0.0
dea%XLEV        = 0.0
dea%XLES        = 0.0
!
dea%XLESL       = 0.0
dea%XSNDRIFT    = 0.0
!
dea%XLER        = 0.0
dea%XLETR       = 0.0
!
dea%XDRAIN      = 0.0
dea%XRUNOFF     = 0.0
dea%XDRIP       = 0.0
dea%XRRVEG      = 0.0
dea%XMELT       = 0.0
!
dea%XIRRIG_FLUX = 0.0
!
dea%XGPP        = 0.0
dea%XRESP_AUTO  = 0.0
dea%XRESP_ECO   = 0.0  
!
dea%XQSB        = 0.0
dea%XHORT       = 0.0
!
dea%XIFLOOD    = 0.0
dea%XPFLOOD    = 0.0
dea%XLE_FLOOD  = 0.0
dea%XLEI_FLOOD = 0.0
!
dea%XRN_SN_FR    = 0.0
dea%XH_SN_FR     = 0.0
dea%XLEI_SN_FR   = 0.0
dea%XLE_SN_FR    = 0.0
dea%XGFLUX_SN_FR = 0.0
dea%XLEG_SN_FR   = 0.0
dea%XLEGI_SN_FR  = 0.0
dea%XLEV_SN_FR   = 0.0
dea%XLETR_SN_FR  = 0.0
dea%XUSTAR_SN_FR = 0.0
dea%XLER_SN_FR   = 0.0
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_EVAP_BUD',1,zhook_handle)
!
END SUBROUTINE init_evap_bud
!
SUBROUTINE init_meb_bud(DEA)
!
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_MEB_BUD',0,zhook_handle)
!
dea%XLELITTER      = 0.0
dea%XLELITTERI     = 0.0
dea%XDRIPLIT       = 0.0
dea%XRRLIT         = 0.0
!
dea%XLEV_CV        = 0.0
dea%XLES_CV        = 0.0
dea%XLETR_CV       = 0.0
dea%XLE_CV         = 0.0
dea%XH_CV          = 0.0
dea%XMELT_CV       = 0.0
dea%XFRZ_CV        = 0.0
!
dea%XLETR_GV       = 0.0
dea%XLER_GV        = 0.0
dea%XLE_GV         = 0.0
dea%XH_GV          = 0.0
!
dea%XLE_GN         = 0.0
dea%XEVAP_GN       = 0.0
dea%XH_GN          = 0.0
dea%XSR_GN         = 0.0
dea%XSWDOWN_GN     = 0.0
dea%XLWDOWN_GN     = 0.0
!
dea%XEVAP_G        = 0.0
dea%XLE_CA         = 0.0
dea%XH_CA          = 0.0
!
dea%XSWUP          = 0.0
dea%XLWUP          = 0.0
!
dea%XSWNET_V       = 0.0
dea%XSWNET_G       = 0.0
dea%XSWNET_N       = 0.0
dea%XSWNET_NS      = 0.0
dea%XLWNET_V       = 0.0
dea%XLWNET_G       = 0.0
dea%XLWNET_N       = 0.0
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_MEB_BUD',1,zhook_handle)
!
END SUBROUTINE init_meb_bud
!
SUBROUTINE init_water_bud(DEA)
!
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_WATER_BUD',0,zhook_handle)
!
dea%XDWG    = 0.0
dea%XDWGI   = 0.0
dea%XDWR    = 0.0
dea%XDSWE   = 0.0
dea%XWATBUD = 0.0
!
IF (lhook) CALL dr_hook('MODE_DIAG:INIT_WATER_BUD',1,zhook_handle)
!
END SUBROUTINE init_water_bud
!
SUBROUTINE avg_diag_tstep_evap(PTSTEP, DEAC, DEA)
!
REAL, INTENT(IN) :: PTSTEP
TYPE(diag_evap_isba_t), INTENT(IN) :: DEAC
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:AVG_DIAG_TSTEP_EVAP',0,zhook_handle)
!
dea%XLEG    (:) = deac%XLEG    (:)/ptstep
dea%XLEGI   (:) = deac%XLEGI   (:)/ptstep
dea%XLEV    (:) = deac%XLEV    (:)/ptstep
dea%XLES    (:) = deac%XLES    (:)/ptstep
!
dea%XLESL   (:) = deac%XLESL   (:)/ptstep
dea%XSNDRIFT(:) = deac%XSNDRIFT(:)/ptstep
!
dea%XLER    (:) = deac%XLER    (:)/ptstep
dea%XLETR   (:) = deac%XLETR   (:)/ptstep
dea%XDRAIN  (:) = deac%XDRAIN  (:)/ptstep
dea%XQSB    (:) = deac%XQSB    (:)/ptstep
dea%XRUNOFF (:) = deac%XRUNOFF (:)/ptstep
dea%XHORT   (:) = deac%XHORT   (:)/ptstep
dea%XDRIP   (:) = deac%XDRIP   (:)/ptstep
dea%XRRVEG  (:) = deac%XRRVEG  (:)/ptstep
dea%XMELT   (:) = deac%XMELT   (:)/ptstep
!
IF (lhook) CALL dr_hook('MODE_DIAG:AVG_DIAG_TSTEP_EVAP',1,zhook_handle)
!
END SUBROUTINE avg_diag_tstep_evap
!
SUBROUTINE avg_diag_tstep_water(PTSTEP, DEAC, DEA)
!
REAL, INTENT(IN) :: PTSTEP
TYPE(diag_evap_isba_t), INTENT(IN) :: DEAC
TYPE(diag_evap_isba_t), INTENT(INOUT) :: DEA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:AVG_DIAG_TSTEP_WATER',0,zhook_handle)
!
dea%XDWG   (:) = deac%XDWG   (:)/ptstep
dea%XDWGI  (:) = deac%XDWGI  (:)/ptstep
dea%XDWR   (:) = deac%XDWR   (:)/ptstep
dea%XDSWE  (:) = deac%XDSWE  (:)/ptstep
dea%XWATBUD(:) = deac%XWATBUD(:)/ptstep
!
IF (lhook) CALL dr_hook('MODE_DIAG:AVG_DIAG_TSTEP_WATER',1,zhook_handle)
!
END SUBROUTINE avg_diag_tstep_water
!
SUBROUTINE avg_diag_tstep_surf(PTSTEP, DAC, DA)
!
REAL, INTENT(IN) :: PTSTEP
TYPE(diag_t), INTENT(IN) :: DAC
TYPE(diag_t), INTENT(INOUT) :: DA
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('MODE_DIAG:AVG_DIAG_TSTEP_SURF',0,zhook_handle)
!
da%XRN   (:) = dac%XRN   (:)/ptstep
da%XH    (:) = dac%XH    (:)/ptstep
da%XLE   (:) = dac%XLE   (:)/ptstep
da%XLEI  (:) = dac%XLEI  (:)/ptstep
da%XGFLUX(:) = dac%XGFLUX(:)/ptstep
!
da%XEVAP (:) = dac%XEVAP (:)/ptstep
da%XSUBL (:) = dac%XSUBL (:)/ptstep
!
da%XSWD(:) = dac%XSWD(:)/ptstep
da%XSWU(:) = dac%XSWU(:)/ptstep
da%XLWD(:) = dac%XLWD(:)/ptstep
da%XLWU(:) = dac%XLWU(:)/ptstep
!
IF (lhook) CALL dr_hook('MODE_DIAG:AVG_DIAG_TSTEP_SURF',1,zhook_handle)
!
END SUBROUTINE avg_diag_tstep_surf
!
!-------------------------------------------------------------------------------
!
END MODULE mode_diag