diag_surf_budgetc.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_surf_budgetc (D, DC, PTSTEP, ONOTICE)  
!     #########################################################################
!
!!****  *DIAG_SURF_BUDGETC * - Computes cumulated diagnostics 
!!
!!    PURPOSE
!!    -------
!
!!**  METHOD
!!    ------
!!
!!    REFERENCE
!!    ---------
!!      
!!
!!    AUTHOR
!!    ------
!!     B. Decharme 
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    08/2009
!!------------------------------------------------------------------
! 
USE modd_diag_n, ONLY : diag_t
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
TYPE(diag_t), INTENT(INOUT) :: D
TYPE(diag_t), INTENT(INOUT) :: DC
!
REAL,               INTENT(IN) :: PTSTEP  
!
LOGICAL, INTENT(IN) :: ONOTICE
!
!*      0.2    declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------------
!
!* total incoming and outgoing SW
!
IF (lhook) CALL dr_hook('DIAG_SURF_BUDGETC_FLAKE',0,zhook_handle)
!
IF (onotice) dc%XSWD(:) = dc%XSWD(:) + d%XSWD(:) * ptstep
dc%XSWU(:) = dc%XSWU(:) + d%XSWU(:) * ptstep
!
!*incoming outgoing LW
!
IF (onotice) dc%XLWD(:) = dc%XLWD(:) + d%XLWD(:) * ptstep
dc%XLWU(:) = dc%XLWU(:) + d%XLWU(:) * ptstep
!
!* net radiation
!
dc%XRN(:) = dc%XRN(:) + d%XRN(:) * ptstep
!
!* sensible heat flux
!
dc%XH(:) = dc%XH(:) + d%XH(:) * ptstep 
!
IF (onotice) THEN
  !
  !* latent heat flux
  !
  dc%XLE (:) = dc%XLE (:) + d%XLE (:) * ptstep 
  dc%XLEI(:) = dc%XLEI(:) + d%XLEI(:) * ptstep 
  !
  !* evaporation and sublimation (kg/m2)
  !
  dc%XEVAP(:) = dc%XEVAP(:) + d%XEVAP(:) * ptstep
  dc%XSUBL(:) = dc%XSUBL(:) + d%XSUBL(:) * ptstep
  !
ENDIF
!
!* storage flux
!
dc%XGFLUX(:) = dc%XGFLUX(:) + d%XGFLUX(:) * ptstep 
!
!* wind stress
!
dc%XFMU(:) = dc%XFMU(:) + d%XFMU(:) * ptstep 
dc%XFMV(:) = dc%XFMV(:) + d%XFMV(:) * ptstep
!
IF (lhook) CALL dr_hook('DIAG_SURF_BUDGETC',1,zhook_handle)
!
!-------------------------------------------------------------------------------------
!
END SUBROUTINE diag_surf_budgetc