modd_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.
!     ######################
      MODULE modd_diag_teb_n
!     ######################
!
!!****  *MODD_DIAG_TEB - declaration of diagnostics for TEB scheme
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      V. Masson   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       01/2004
!!      Modified    01/2006 : sea flux parameterization.
!
!*       0.   DECLARATIONS
!             ------------
!
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE

TYPE diag_teb_t
!------------------------------------------------------------------------------
!
  REAL    :: XDIAG_TSTEP  ! time step for diagnostics writing
!
  INTEGER :: N2M          ! flag for 2 meters (and 10 meters) quantities
  LOGICAL :: L2M_MIN_ZS   ! flag for 2 meters quantities evaluated on
!                         ! the minimum orographyy of the grid      
  LOGICAL :: LSURF_BUDGET ! flag for surface energy budget
  LOGICAL :: LRAD_BUDGET  ! flag for radiative energy budget
  LOGICAL :: LCOEF        ! flag for transfer coefficients
  LOGICAL :: LSURF_VARS   ! flag for surface variables
!
  LOGICAL :: LPGD         ! flag for writing of PGD files
  LOGICAL :: LPGD_FIX     ! flag for writing of PGD files for time
!                           invariant field
!
!* averaged variables
!
  REAL, POINTER, DIMENSION(:)   :: XRI      ! Bulk-Richardson number           (-)
  REAL, POINTER, DIMENSION(:)   :: XCD      ! drag coefficient for wind        (W/s2)
  REAL, POINTER, DIMENSION(:)   :: XCH      ! drag coefficient for heat        (W/s)
  REAL, POINTER, DIMENSION(:)   :: XCE      ! drag coefficient for vapor       (W/s/K)
  REAL, POINTER, DIMENSION(:)   :: XZ0      ! roughness length for momentum    (m)
  REAL, POINTER, DIMENSION(:)   :: XZ0H     ! roughness length for heat        (m)
  REAL, POINTER, DIMENSION(:)   :: XRN      ! net radiation at surface         (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XH       ! sensible heat flux               (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XLE      ! latent heat flux                 (W/m2) 
  REAL, POINTER, DIMENSION(:)   :: XGFLUX   ! net soil-vegetation flux         (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XT2M     ! air temperature at 2 meters      (K)
  REAL, POINTER, DIMENSION(:)   :: XT2M_MIN ! Minimum air temperature at 2 meters      (K)
  REAL, POINTER, DIMENSION(:)   :: XT2M_MAX ! Maximum air temperature at 2 meters      (K)
  REAL, POINTER, DIMENSION(:)   :: XQ2M     ! air humidity at 2 meters         (kg/kg)
  REAL, POINTER, DIMENSION(:)   :: XHU2M    ! air relative humidity at 2 meters(-)
  REAL, POINTER, DIMENSION(:)   :: XHU2M_MIN! Minimum relative humidity at 2 meters    (-)
  REAL, POINTER, DIMENSION(:)   :: XHU2M_MAX! Maximum relative humidity at 2 meters    (-)
  REAL, POINTER, DIMENSION(:)   :: XQS      ! air humidity at surface          (kg/kg)
  REAL, POINTER, DIMENSION(:)   :: XZON10M  ! zonal wind at 10 meters          (m/s)
  REAL, POINTER, DIMENSION(:)   :: XMER10M  ! meridian wind at 10 meters       (m/s)
  REAL, POINTER, DIMENSION(:)   :: XWIND10M ! wind at 10 meters                (m/s)
  REAL, POINTER, DIMENSION(:)   :: XWIND10M_MAX! Maximum wind at 10 meters     (m/s)
  REAL, POINTER, DIMENSION(:)   :: XSFCO2   ! CO2 flux                         (m/s*kg_CO2/kg_air)
  REAL, POINTER, DIMENSION(:)   :: XLWD     ! downward long wave radiation     (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XLWU     ! upward long wave radiation       (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XSWD     ! downward short wave radiation    (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XSWU     ! upward short wave radiation      (W/m2)
  REAL, POINTER, DIMENSION(:,:) :: XSWBD    ! downward short wave radiation by spectral band   (W/m2)
  REAL, POINTER, DIMENSION(:,:) :: XSWBU    ! upward short wave radiation by spectral band (W/m2)
  REAL, POINTER, DIMENSION(:)   :: XFMU     ! horizontal momentum flux zonal   (m2/s2)
  REAL, POINTER, DIMENSION(:)   :: XFMV     ! horizontal momentum flux meridian (m2/s2)             
  REAL, POINTER, DIMENSION(:)   :: XDIAG_TS ! arithmetic mean of surface temperature (K)
!------------------------------------------------------------------------------
!

END TYPE diag_teb_t



 CONTAINS

!




SUBROUTINE diag_teb_init(YDIAG_TEB)
TYPE(diag_teb_t), INTENT(INOUT) :: ydiag_teb
REAL(KIND=JPRB) :: zhook_handle
IF (lhook) CALL dr_hook("MODD_DIAG_TEB_N:DIAG_TEB_INIT",0,zhook_handle)
  nullify(ydiag_teb%XRI)
  nullify(ydiag_teb%XCD)
  nullify(ydiag_teb%XCH)
  nullify(ydiag_teb%XCE)
  nullify(ydiag_teb%XZ0)
  nullify(ydiag_teb%XZ0H)
  nullify(ydiag_teb%XRN)
  nullify(ydiag_teb%XH)
  nullify(ydiag_teb%XLE)
  nullify(ydiag_teb%XGFLUX)
  nullify(ydiag_teb%XT2M)
  nullify(ydiag_teb%XT2M_MIN)
  nullify(ydiag_teb%XT2M_MAX)
  nullify(ydiag_teb%XQ2M)
  nullify(ydiag_teb%XHU2M)
  nullify(ydiag_teb%XHU2M_MIN)
  nullify(ydiag_teb%XHU2M_MAX)
  nullify(ydiag_teb%XQS)
  nullify(ydiag_teb%XZON10M)
  nullify(ydiag_teb%XMER10M)
  nullify(ydiag_teb%XWIND10M)
  nullify(ydiag_teb%XWIND10M_MAX)
  nullify(ydiag_teb%XSFCO2)
  nullify(ydiag_teb%XLWD)
  nullify(ydiag_teb%XLWU)
  nullify(ydiag_teb%XSWD)
  nullify(ydiag_teb%XSWU)
  nullify(ydiag_teb%XSWBD)
  nullify(ydiag_teb%XSWBU)
  nullify(ydiag_teb%XFMU)
  nullify(ydiag_teb%XFMV)
  nullify(ydiag_teb%XDIAG_TS)
ydiag_teb%XDIAG_TSTEP=0.
ydiag_teb%N2M=0
ydiag_teb%L2M_MIN_ZS=.false.
ydiag_teb%LSURF_BUDGET=.false.
ydiag_teb%LRAD_BUDGET=.false.
ydiag_teb%LCOEF=.false.
ydiag_teb%LSURF_VARS=.false.
ydiag_teb%LPGD=.false.
ydiag_teb%LPGD_FIX=.false.
IF (lhook) CALL dr_hook("MODD_DIAG_TEB_N:DIAG_TEB_INIT",1,zhook_handle)
END SUBROUTINE diag_teb_init


END MODULE modd_diag_teb_n