alloc_trip_diag.F90

Go to the documentation of this file.

!     #########
      SUBROUTINE alloc_trip_diag (TPDG,KLON,KLAT)
!     #######################################################################
!
!!****  *ALLOC_TRIP_DIAG*  
!!
!!    PURPOSE
!!    -------
!
!     Define the name and unit of each trip output variable.
!     
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    01/2017 
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
!
USE modd_trip_diag, ONLY : trip_diag_t
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!
!*      0.1    declarations of arguments
!
TYPE(trip_diag_t), INTENT(INOUT) :: TPDG
!
INTEGER, INTENT(IN)          :: KLON, KLAT
!
!*      0.2    declarations of output variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('ALLOC_TRIP_DIAG',0,zhook_handle)
!
! * Allocate and initialyse diagnostic variables
!
 CALL alloc_diag(tpdg%TDIAG)
!
! * Allocate and initialyse run mean diagnostic variables
!
 CALL alloc_diag(tpdg%TDIAG_RUN)
!
IF (lhook) CALL dr_hook('ALLOC_TRIP_DIAG',1,zhook_handle)
!
!-------------------------------------------------------------------------------
 CONTAINS
!-------------------------------------------------------------------------------
!
SUBROUTINE alloc_diag(TPDIAG)
!
USE modn_trip,       ONLY : cgroundw, cvit, lflood
USE modn_trip_run,   ONLY : ldiag_misc
USE modd_trip_oasis, ONLY : lcpl_land, lcpl_flood
!
USE modd_type_diag
!
IMPLICIT NONE
!
TYPE(diag) :: TPDIAG
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
IF (lhook) CALL dr_hook('ALLOC_TRIP_DIAG:ALLOC_DIAG',0,zhook_handle)
!
ALLOCATE(tpdiag%XSURF_STO(klon,klat))
ALLOCATE(tpdiag%XQDIS    (klon,klat))
tpdiag%XSURF_STO(:,:) = 0.0
tpdiag%XQDIS    (:,:) = 0.0

IF(ldiag_misc)THEN
  ALLOCATE(tpdiag%XQIN     (klon,klat))
  tpdiag%XQIN      (:,:) = 0.0
ELSE
  ALLOCATE(tpdiag%XQIN(0,0))
ENDIF
!
IF(lcpl_land.AND.ldiag_misc)THEN
  ALLOCATE(tpdiag%XRUNOFF(klon,klat))
  ALLOCATE(tpdiag%XDRAIN (klon,klat))
  tpdiag%XRUNOFF(:,:) = 0.0
  tpdiag%XDRAIN (:,:) = 0.0
ELSE
  ALLOCATE(tpdiag%XRUNOFF(0,0))
  ALLOCATE(tpdiag%XDRAIN (0,0))
ENDIF

IF(cgroundw/='DEF')THEN
  ALLOCATE(tpdiag%XGROUND_STO (klon,klat))
  ALLOCATE(tpdiag%XQGF        (klon,klat))
  tpdiag%XGROUND_STO (:,:) = 0.0
  tpdiag%XQGF        (:,:) = 0.0
ELSE
  ALLOCATE(tpdiag%XGROUND_STO(0,0))
  ALLOCATE(tpdiag%XQGF       (0,0))
ENDIF

IF(cgroundw=='DIF')THEN
!
  ALLOCATE(tpdiag%XHGROUND(klon,klat))
  ALLOCATE(tpdiag%XFWTD   (klon,klat))
  ALLOCATE(tpdiag%XWTD    (klon,klat))
  tpdiag%XHGROUND(:,:) = 0.0
  tpdiag%XFWTD   (:,:) = 0.0
  tpdiag%XWTD    (:,:) = 0.0
!  
  IF(ldiag_misc)THEN
    ALLOCATE(tpdiag%XQGCELL (klon,klat))
    ALLOCATE(tpdiag%XHGHS   (klon,klat))
    tpdiag%XQGCELL (:,:) = 0.0
    tpdiag%XHGHS   (:,:) = 0.0
  ELSE
    ALLOCATE(tpdiag%XQGCELL (0,0))
    ALLOCATE(tpdiag%XHGHS   (0,0))          
  ENDIF
!
ELSE
  ALLOCATE(tpdiag%XFWTD   (0,0))
  ALLOCATE(tpdiag%XWTD    (0,0))
  ALLOCATE(tpdiag%XHGROUND(0,0))
  ALLOCATE(tpdiag%XQGCELL (0,0))
  ALLOCATE(tpdiag%XHGHS   (0,0))
ENDIF


IF(cvit=='VAR')THEN
  ALLOCATE(tpdiag%XHS (klon,klat))
  ALLOCATE(tpdiag%XVEL(klon,klat))
  tpdiag%XHS (:,:) = 0.0
  tpdiag%XVEL(:,:) = 0.0
ELSE
  ALLOCATE(tpdiag%XHS (0,0))
  ALLOCATE(tpdiag%XVEL(0,0))
ENDIF
!
!
!
IF (lcpl_flood.AND.ldiag_misc) THEN
   ALLOCATE(tpdiag%XSOURCE(klon,klat))
   tpdiag%XSOURCE(:,:) = 0.0
ELSE
   ALLOCATE(tpdiag%XSOURCE(0,0))
ENDIF
!
!
!
IF(lflood)THEN
  ALLOCATE(tpdiag%XFLOOD_STO (klon,klat))
  ALLOCATE(tpdiag%XFF        (klon,klat))
  ALLOCATE(tpdiag%XHF        (klon,klat))
  tpdiag%XFLOOD_STO (:,:) = 0.0
  tpdiag%XFF        (:,:) = 0.0
  tpdiag%XHF        (:,:) = 0.0
  IF(ldiag_misc)THEN
    ALLOCATE(tpdiag%XQFR   (klon,klat))
    ALLOCATE(tpdiag%XQRF   (klon,klat))
    ALLOCATE(tpdiag%XVFIN  (klon,klat))
    ALLOCATE(tpdiag%XVFOUT (klon,klat))
    ALLOCATE(tpdiag%XWF    (klon,klat))
    ALLOCATE(tpdiag%XLF    (klon,klat))
    ALLOCATE(tpdiag%XHSF   (klon,klat))
    tpdiag%XQFR   (:,:) = 0.0
    tpdiag%XQRF   (:,:) = 0.0
    tpdiag%XVFIN  (:,:) = 0.0
    tpdiag%XVFOUT (:,:) = 0.0
    tpdiag%XWF    (:,:) = 0.0
    tpdiag%XLF    (:,:) = 0.0
    tpdiag%XHSF   (:,:) = 0.0
  ELSE
    ALLOCATE(tpdiag%XQFR   (0,0))
    ALLOCATE(tpdiag%XQRF   (0,0))
    ALLOCATE(tpdiag%XVFIN  (0,0))
    ALLOCATE(tpdiag%XVFOUT (0,0))
    ALLOCATE(tpdiag%XWF    (0,0))
    ALLOCATE(tpdiag%XLF    (0,0))
    ALLOCATE(tpdiag%XHSF   (0,0))
  ENDIF
ELSE
  ALLOCATE(tpdiag%XFLOOD_STO (0,0))
  ALLOCATE(tpdiag%XHF        (0,0))
  ALLOCATE(tpdiag%XFF        (0,0))
  ALLOCATE(tpdiag%XQFR       (0,0))
  ALLOCATE(tpdiag%XQRF       (0,0))
  ALLOCATE(tpdiag%XVFIN      (0,0))
  ALLOCATE(tpdiag%XVFOUT     (0,0))
  ALLOCATE(tpdiag%XWF        (0,0))
  ALLOCATE(tpdiag%XLF        (0,0))
  ALLOCATE(tpdiag%XHSF       (0,0))
ENDIF
!
IF (lhook) CALL dr_hook('ALLOC_TRIP_DIAG:ALLOC_DIAG',1,zhook_handle)
!
END SUBROUTINE alloc_diag
!
!-------------------------------------------------------------------------------
END SUBROUTINE alloc_trip_diag