BROWSER/out_doc80/diag__teb__initn_8F90_source.html

 !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_teb_init_n (DGT, DGUT, &

                                   hprogram,klu,ksw)

 !     #####################

 !

 !!****  *DIAG_TEB_INIT_n* - routine to initialize TEB diagnostic variables

 !!

 !!    PURPOSE

 !!    -------

 !!

 !!**  METHOD

 !!    ------

 !!

 !!    EXTERNAL

 !!    --------

 !!

 !!

 !!    IMPLICIT ARGUMENTS

 !!    ------------------

 !!

 !!    REFERENCE

 !!    ---------

 !!

 !!

 !!    AUTHOR

 !!    ------

 !!      V. Masson   *Meteo France*

 !!

 !!    MODIFICATIONS

 !!    -------------

 !!      Original    01/2004

 !!       V. Masson  10/2013 Adds integrated UTCI diagnostics

 !       B. decharme 04/2013 : Add DIAG_TS

 !-------------------------------------------------------------------------------

 !

 !*       0.    DECLARATIONS

 !              ------------

 !

 !

 USE modd_diag_teb_n, ONLY : diag_teb_t

 USE modd_diag_utci_teb_n, ONLY : diag_utci_teb_t

 !

 USE modd_surf_par,   ONLY : xundef

 USE modd_type_date_surf

 !

 USE modd_utci,              ONLY : nutci_stress


 !

 USE modi_read_surf

 !

 !

 USE yomhook   ,ONLY : lhook,   dr_hook

 USE parkind1  ,ONLY : jprb

 !

 IMPLICIT NONE

 !

 !*       0.1   Declarations of arguments

 !              -------------------------

 !

 !

 TYPE(diag_teb_t), INTENT(INOUT) :: dgt

 TYPE(diag_utci_teb_t), INTENT(INOUT) :: dgut

 !

 INTEGER, INTENT(IN) :: klu   ! size of arrays

 INTEGER, INTENT(IN) :: ksw   ! spectral bands

  CHARACTER(LEN=6), INTENT(IN):: hprogram  ! program calling

 !

 !*       0.2   Declarations of local variables

 !              -------------------------------

 !

 INTEGER           :: iresp          ! IRESP  : return-code if a problem appears

  CHARACTER(LEN=12) :: yrec           ! Name of the article to be read

 REAL(KIND=JPRB) :: zhook_handle

 !

 !-------------------------------------------------------------------------------

 !

 !* surface energy budget

 !

 IF (lhook) CALL dr_hook('DIAG_TEB_INIT_N',0,zhook_handle)

 !

 ALLOCATE(dgt%XDIAG_TS(klu))

 dgt%XDIAG_TS = xundef

 !

 IF (dgt%LSURF_BUDGET) THEN

   ALLOCATE(dgt%XRN     (klu))

   ALLOCATE(dgt%XH      (klu))

   ALLOCATE(dgt%XLE     (klu))

   ALLOCATE(dgt%XGFLUX  (klu))

   ALLOCATE(dgt%XSWD    (klu))

   ALLOCATE(dgt%XSWU    (klu))

   ALLOCATE(dgt%XSWBD   (klu,ksw))

   ALLOCATE(dgt%XSWBU   (klu,ksw))

   ALLOCATE(dgt%XLWD    (klu))

   ALLOCATE(dgt%XLWU    (klu))

   ALLOCATE(dgt%XFMU    (klu))

   ALLOCATE(dgt%XFMV    (klu))

   ALLOCATE(dgt%XSFCO2  (klu))

   !

   dgt%XRN      = xundef

   dgt%XH       = xundef

   dgt%XLE      = xundef

   dgt%XGFLUX   = xundef

   dgt%XSWD     = xundef

   dgt%XSWU     = xundef

   dgt%XSWBD    = xundef

   dgt%XSWBU    = xundef

   dgt%XLWD     = xundef

   dgt%XLWU     = xundef

   dgt%XFMU     = xundef

   dgt%XFMV     = xundef

   dgt%XSFCO2   = xundef

 ELSE

   ALLOCATE(dgt%XRN     (0))

   ALLOCATE(dgt%XH      (0))

   ALLOCATE(dgt%XLE     (0))

   ALLOCATE(dgt%XGFLUX  (0))

   ALLOCATE(dgt%XSWD    (0))

   ALLOCATE(dgt%XSWU    (0))

   ALLOCATE(dgt%XSWBD   (0,0))

   ALLOCATE(dgt%XSWBU   (0,0))

   ALLOCATE(dgt%XLWD    (0))

   ALLOCATE(dgt%XLWU    (0))

   ALLOCATE(dgt%XFMU    (0))

   ALLOCATE(dgt%XFMV    (0))

   ALLOCATE(dgt%XSFCO2  (0))

 END IF

 !

 !* parameters at 2m

 !

 IF (dgt%N2M>=1) THEN

   ALLOCATE(dgt%XRI     (klu))

   ALLOCATE(dgt%XT2M    (klu))

   ALLOCATE(dgt%XT2M_MIN (klu))

   ALLOCATE(dgt%XT2M_MAX (klu))

   ALLOCATE(dgt%XQ2M    (klu))

   ALLOCATE(dgt%XHU2M   (klu))

   ALLOCATE(dgt%XHU2M_MIN(klu))

   ALLOCATE(dgt%XHU2M_MAX(klu))

   ALLOCATE(dgt%XZON10M (klu))

   ALLOCATE(dgt%XMER10M (klu))

   ALLOCATE(dgt%XWIND10M (klu))

   ALLOCATE(dgt%XWIND10M_MAX(klu))

   !

   dgt%XRI      = xundef

   dgt%XT2M     = xundef

   dgt%XT2M_MIN = xundef

   dgt%XT2M_MAX = -xundef

   dgt%XQ2M     = xundef

   dgt%XHU2M    = xundef

   dgt%XHU2M_MIN= xundef

   dgt%XHU2M_MAX=-xundef

   dgt%XZON10M  = xundef

   dgt%XMER10M  = xundef

   dgt%XWIND10M = xundef

   dgt%XWIND10M_MAX = -xundef

 ELSE

   ALLOCATE(dgt%XRI      (0))

   ALLOCATE(dgt%XT2M     (0))

   ALLOCATE(dgt%XT2M_MIN (0))

   ALLOCATE(dgt%XT2M_MAX (0))

   ALLOCATE(dgt%XQ2M     (0))

   ALLOCATE(dgt%XHU2M    (0))

   ALLOCATE(dgt%XHU2M_MIN(0))

   ALLOCATE(dgt%XHU2M_MAX(0))

   ALLOCATE(dgt%XZON10M  (0))

   ALLOCATE(dgt%XMER10M  (0))

   ALLOCATE(dgt%XWIND10M (0))

   ALLOCATE(dgt%XWIND10M_MAX(0))

 END IF

 !!

 !* miscellaneous fields

 !

 IF (dgt%N2M>0 .AND. dgut%LUTCI) THEN

   !

   ALLOCATE(dgut%XUTCI_IN       (klu))

   ALLOCATE(dgut%XUTCI_OUTSUN   (klu))

   ALLOCATE(dgut%XUTCI_OUTSHADE (klu))

   ALLOCATE(dgut%XTRAD_SUN      (klu))

   ALLOCATE(dgut%XTRAD_SHADE    (klu))

   ALLOCATE(dgut%XUTCIC_IN      (klu,nutci_stress))

   ALLOCATE(dgut%XUTCIC_OUTSUN  (klu,nutci_stress))

   ALLOCATE(dgut%XUTCIC_OUTSHADE(klu,nutci_stress))

   !

   dgut%XUTCI_IN        = xundef

   dgut%XUTCI_OUTSUN    = xundef

   dgut%XUTCI_OUTSHADE  = xundef

   dgut%XTRAD_SUN       = xundef

   dgut%XTRAD_SHADE     = xundef

   dgut%XUTCIC_IN       = 0.

   dgut%XUTCIC_OUTSUN   = 0.

   dgut%XUTCIC_OUTSHADE = 0.

   !

 ELSE

   ALLOCATE(dgut%XUTCI_IN       (0))

   ALLOCATE(dgut%XUTCI_OUTSUN   (0))

   ALLOCATE(dgut%XUTCI_OUTSHADE (0))

   ALLOCATE(dgut%XTRAD_SUN      (0))

   ALLOCATE(dgut%XTRAD_SHADE    (0))

   ALLOCATE(dgut%XUTCIC_IN      (0,0))

   ALLOCATE(dgut%XUTCIC_OUTSUN  (0,0))

   ALLOCATE(dgut%XUTCIC_OUTSHADE(0,0))

 ENDIF

 !

 !* transfer coefficients

 !

 IF (dgt%LCOEF) THEN

   ALLOCATE(dgt%XCD     (klu))

   ALLOCATE(dgt%XCH     (klu))

   ALLOCATE(dgt%XCE     (klu))

   ALLOCATE(dgt%XZ0     (klu))

   ALLOCATE(dgt%XZ0H    (klu))

   !

   dgt%XCD      = xundef

   dgt%XCH      = xundef

   dgt%XCE      = xundef

   dgt%XZ0      = xundef

   dgt%XZ0H     = xundef

 ELSE

   ALLOCATE(dgt%XCD     (0))

   ALLOCATE(dgt%XCH     (0))

   ALLOCATE(dgt%XCE     (0))

   ALLOCATE(dgt%XZ0     (0))

   ALLOCATE(dgt%XZ0H    (0))

 END IF

 !

 !

 !* surface humidity

 !

 IF (dgt%LSURF_VARS) THEN

   ALLOCATE(dgt%XQS     (klu))

   !

   dgt%XQS      = xundef

 ELSE

   ALLOCATE(dgt%XQS     (0))

 END IF

 IF (lhook) CALL dr_hook('DIAG_TEB_INIT_N',1,zhook_handle)

 !

 !-------------------------------------------------------------------------------

 !

 END SUBROUTINE diag_teb_init_n

modi_read_surf
Definition: read_surf.F90:6

modd_diag_utci_teb_n::diag_utci_teb_t
Definition: modd_diag_utci_tebn.F90:41

modd_utci
Definition: modd_utci.F90:6

modd_diag_teb_n
Definition: modd_diag_tebn.F90:6

modd_diag_utci_teb_n
Definition: modd_diag_utci_tebn.F90:6

modd_surf_par
Definition: modd_surf_par.F90:6

modd_type_date_surf
Definition: modd_type_date_surf.F90:6

diag_teb_init_n
subroutine diag_teb_init_n(DGT, DGUT, HPROGRAM, KLU, KSW)
Definition: diag_teb_initn.F90:6

modd_diag_teb_n::diag_teb_t
Definition: modd_diag_tebn.F90:41