trip_diag_init.F90

Go to the documentation of this file.

SUBROUTINE trip_diag_init(PSOUT,PSIN,PVEL,PHS,PGOUT,PGNEG,PHG_OLD,   &
                          PWTD,PFWTD,PQGCELL,PHGHS,                  &
                          PQFR,PQRF,PVFIN,PVFOUT,PHSF,PDISCHARGE,    &
                          PGSTO_ALL,PGSTO2_ALL,PGIN_ALL,PGOUT_ALL    )  
!     #####################################################
!
!!****  *TRIP_DIAG_INIT*  
!!
!!    PURPOSE
!!    -------
!
!     TRIP diag compuation
!     
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    12/12/13 
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modn_trip_run,  ONLY : ldiag_misc
USE modn_trip,      ONLY : cgroundw, cvit, lflood
USE modd_trip_par,  ONLY : xundef
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
REAL, DIMENSION(:,:), INTENT(OUT) :: PSOUT      !streamflow                    [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PSIN       !grid-cell input streamflow    [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PVEL       !river velocity                [m/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PHS        !River heigh                   [m]
REAL, DIMENSION(:,:), INTENT(OUT) :: PGOUT      !Groundwater outflow           [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PGNEG      !Groundwater intflow (neg)     [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PWTD       !Water table depth for coupling[m]
REAL, DIMENSION(:,:), INTENT(OUT) :: PFWTD      !Fraction of water table to rise
REAL, DIMENSION(:,:), INTENT(OUT) :: PQGCELL    !lateral groundwater exchanges [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PHGHS      !groundwater minus river heigh [m]
REAL, DIMENSION(:,:), INTENT(OUT) :: PQFR       !floodplains to river exchange [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PQRF       !river to floodplains exchange [kg/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PVFIN      !QRF velocity                  [m/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PVFOUT     !QFR velocity                  [m/s]
REAL, DIMENSION(:,:), INTENT(OUT) :: PHSF       !river minus flodd heigh       [m]
REAL, DIMENSION(:,:), INTENT(OUT) :: PDISCHARGE !Cumulated river discharges    [kg]
REAL, DIMENSION(:,:), INTENT(OUT) :: PHG_OLD    ! Water table depth before     
!                                                 water mass conservation      [m]
!
REAL                , INTENT(OUT) :: PGSTO_ALL  !Global groundwater storage at t    [kg]
REAL                , INTENT(OUT) :: PGSTO2_ALL !Global groundwater storage at t-1  [kg]
REAL                , INTENT(OUT) :: PGIN_ALL   !Global gw recharge + lateral input [kg/m2/s]
REAL                , INTENT(OUT) :: PGOUT_ALL  !Global gw outflow                  [kg/m2/s]
!
!*      0.2    declarations of local variables
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('TRIP_DIAG_INIT',0,zhook_handle)
!
!
pdischarge(:,:) = 0.0
psout(:,:) = 0.0
psin(:,:) = 0.0
pgout(:,:) = 0.0
pgneg(:,:) = 0.0
!
pvel(:,:) = xundef        
phs(:,:) = xundef    
pqfr(:,:) = xundef  
pqrf(:,:) = xundef
pvfin(:,:) = xundef
pvfout(:,:) = xundef
phsf(:,:) = xundef
pqgcell(:,:) = xundef
pwtd(:,:) = xundef
pfwtd(:,:) = xundef
phghs(:,:) = xundef
phg_old(:,:) = xundef
!
IF(cgroundw=='DIF')THEN
   IF(ldiag_misc)THEN
     pqgcell(:,:) = 0.0
     phghs(:,:) = 0.0
   ENDIF
ENDIF
!
IF(cvit=='VAR')THEN
  pvel(:,:) = 0.0
  phs(:,:) = 0.0
ENDIF
!
IF(lflood)THEN 
  pqfr(:,:) = 0.0
  pqrf(:,:) = 0.0
  pvfin(:,:) = 0.0
  pvfout(:,:) = 0.0
  phsf(:,:) = 0.0
ENDIF
!
pgsto_all  = xundef
pgsto2_all = xundef
pgin_all   = xundef
pgout_all  = xundef
!
!
IF (lhook) CALL dr_hook('TRIP_DIAG_INIT',1,zhook_handle)
!-------------------------------------------------------------------------------
!
END SUBROUTINE trip_diag_init