modd_trip.F90

Go to the documentation of this file.

!###############
MODULE modd_trip
!###############
!
!!****  *MODD_TRIP - declaration of surface variable for TRIP RRM
!!
!!    PURPOSE
!!    -------
!
!!
!!**  IMPLICIT ARGUMENTS
!!    ------------------
!!      None 
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!      B. Decharme   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original       21/05/08
!!      B. Decharme 10/2016  bug surface/groundwater coupling
!
!*       0.   DECLARATIONS
!             ------------
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE

TYPE trip_t
!
!-------------------------------------------------------------------------------
!
! Input river geometry Parameters :
!
REAL, POINTER, DIMENSION(:,:) :: xtaug          ! ground water transfer time   [s]
REAL, POINTER, DIMENSION(:,:) :: xslopebed      ! river bed slopes             [m/m]
REAL, POINTER, DIMENSION(:,:) :: xwidth         ! river widths                 [m]
REAL, POINTER, DIMENSION(:,:) :: xn             ! Manning roughness coeficient [-] (0.03 to 0.065)
REAL, POINTER, DIMENSION(:,:) :: xn_flood       ! Manning coeficient over floodplains  [-] (currently 0.1)
REAL, POINTER, DIMENSION(:,:) :: xhc_bed        ! River bed depth              [m]
REAL, POINTER, DIMENSION(:,:) :: xweff          ! Porosité efficace
REAL, POINTER, DIMENSION(:,:) :: xtrans         ! Transmissivité
REAL, POINTER, DIMENSION(:,:) :: xnum_aqui      ! Numéro aquifère
REAL, POINTER, DIMENSION(:,:) :: xnear_aqui     ! Nearest aquifer
REAL, POINTER, DIMENSION(:,:) :: xtopo_riv      ! River elevation              [m]
!
!-------------------------------------------------------------------------------
!
! Time varing variables :
!
REAL, POINTER, DIMENSION(:,:) :: xsurf_sto        ! river channel storage        [kg]
REAL, POINTER, DIMENSION(:,:) :: xground_sto      ! groundwater storage          [kg]
REAL, POINTER, DIMENSION(:,:) :: xflood_sto       ! Floodplain water storage     [kg]
REAL, POINTER, DIMENSION(:,:) :: xhground         ! Groudwater height            [m]
REAL, POINTER, DIMENSION(:,:) :: xhflood          ! Floodplain water depth       [m]
REAL, POINTER, DIMENSION(:,:) :: xfflood          ! Floodplain grid-cell fraction [-]
REAL, POINTER, DIMENSION(:,:) :: xwflood          ! Floodplain width             [m]
REAL, POINTER, DIMENSION(:,:) :: xflood_len       ! Floodplain lenght            [m]
!
!-------------------------------------------------------------------------------
!
! Floodplain fonctions :
!
REAL, POINTER, DIMENSION(:,:,:) :: xtab_f         ! Flood fraction array
REAL, POINTER, DIMENSION(:,:,:) :: xtab_h         ! Topo height array
REAL, POINTER, DIMENSION(:,:,:) :: xtab_vf        ! Flood volume array
!
!-------------------------------------------------------------------------------
!
! Groundwater fonctions :
!
REAL, POINTER, DIMENSION(:,:,:) :: xtabgw_f       ! Groundwater fraction array
REAL, POINTER, DIMENSION(:,:,:) :: xtabgw_h       ! Topo height array
!
!-------------------------------------------------------------------------------
!
! Coupling variable with other models :
!
REAL, POINTER, DIMENSION(:,:) :: xcpl_fwtd         ! grid-cell fraction of water table to rise
REAL, POINTER, DIMENSION(:,:) :: xcpl_wtd          ! Water table depth            [m]
!
REAL, POINTER, DIMENSION(:,:) :: xcpl_fflood       ! Flood fraction    [-]
REAL, POINTER, DIMENSION(:,:) :: xcpl_piflood      ! Floodplains potential infiltration  [kg/m2]
!
REAL, POINTER, DIMENSION(:,:) :: xcpl_rivdis       ! River discharges             [kg/m2]
REAL, POINTER, DIMENSION(:,:) :: xcpl_calvgre      ! Calving flux over greenland  [kg/m2]
REAL, POINTER, DIMENSION(:,:) :: xcpl_calvant      ! Calving flux over antarctica [kg/m2]
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
END TYPE trip_t
!
 CONTAINS
!
SUBROUTINE trip_init(YTRIP)
TYPE(trip_t), INTENT(INOUT) :: YTRIP
REAL(KIND=JPRB) :: ZHOOK_HANDLE
IF (lhook) CALL dr_hook("MODD_TRIP:TRIP_INIT",0,zhook_handle)
!
  NULLIFY(ytrip%XTAUG)
  NULLIFY(ytrip%XSLOPEBED)
  NULLIFY(ytrip%XWIDTH)
  NULLIFY(ytrip%XN)
  NULLIFY(ytrip%XN_FLOOD)
  NULLIFY(ytrip%XHC_BED)
  NULLIFY(ytrip%XSURF_STO)
  NULLIFY(ytrip%XGROUND_STO)
  NULLIFY(ytrip%XFLOOD_STO)
  NULLIFY(ytrip%XHFLOOD)
  NULLIFY(ytrip%XFFLOOD)
  NULLIFY(ytrip%XWFLOOD)
  NULLIFY(ytrip%XFLOOD_LEN)
  NULLIFY(ytrip%XWEFF)
  NULLIFY(ytrip%XTRANS)
  NULLIFY(ytrip%XNUM_AQUI)
  NULLIFY(ytrip%XNEAR_AQUI)
  NULLIFY(ytrip%XTOPO_RIV)
  NULLIFY(ytrip%XTABGW_H)
  NULLIFY(ytrip%XHGROUND)
  NULLIFY(ytrip%XTAB_F)
  NULLIFY(ytrip%XTAB_H)
  NULLIFY(ytrip%XTAB_VF)
  NULLIFY(ytrip%XTABGW_F)
  NULLIFY(ytrip%XTABGW_H)
!
  NULLIFY(ytrip%XCPL_WTD)
  NULLIFY(ytrip%XCPL_FWTD)
  NULLIFY(ytrip%XCPL_RIVDIS)
  NULLIFY(ytrip%XCPL_FFLOOD)
  NULLIFY(ytrip%XCPL_PIFLOOD)
  NULLIFY(ytrip%XCPL_CALVGRE)
  NULLIFY(ytrip%XCPL_CALVANT)
!
IF (lhook) CALL dr_hook("MODD_TRIP:TRIP_INIT",1,zhook_handle)
END SUBROUTINE trip_init

!
END MODULE modd_trip