trip_diag_run.F90

Go to the documentation of this file.

!     #########
      SUBROUTINE trip_diag_run (TPDG, TPG, &
                                KLISTING,KLON,KLAT,PRUNTIME)  
!     #####################################################
!
!!****  *TRIP_DIAG_RUN*  
!!
!!    PURPOSE
!!    -------
!
!     TRIP river routing run mean outputs.
!     
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    28/05/05 
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
USE modd_trip_diag, ONLY : trip_diag_t
USE modd_trip_grid, ONLY : trip_grid_t
!
USE modn_trip,       ONLY : cgroundw, cvit, lflood
USE modn_trip_run,   ONLY : ldiag_misc
USE modd_trip_oasis, ONLY : lcpl_land
!                           
!
USE mode_rw_trip
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
!
!
TYPE(trip_diag_t), INTENT(INOUT) :: TPDG
TYPE(trip_grid_t), INTENT(INOUT) :: TPG
!
INTEGER, INTENT(IN)             :: KLISTING
INTEGER, INTENT(IN)             :: KLON
INTEGER, INTENT(IN)             :: KLAT
!
REAL, INTENT(IN)                :: PRUNTIME 
!
!*      0.2    declarations of local variables
!
 CHARACTER(LEN=18), PARAMETER         :: YRUN         ='TRIP_DIAG_RUN.nc'
 CHARACTER(LEN=50)                    :: YFILE
 CHARACTER(LEN=10)                    :: YVNAME
!
REAL,   DIMENSION(KLON,KLAT)         :: ZWRITE
LOGICAL,DIMENSION(KLON,KLAT)         :: LMASK
LOGICAL,DIMENSION(KLON,KLAT)         :: LMASK_GW
!
INTEGER :: ITNUM, ITVAL
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
IF (lhook) CALL dr_hook('TRIP_DIAG_RUN',0,zhook_handle)
!
! * Trip mask
!
lmask(:,:) = tpg%GMASK(:,:)
!
! * Groundwater specific mask
!
IF(cgroundw/='DEF')THEN
  lmask_gw(:,:) = tpg%GMASK_GW(:,:)        
ENDIF
!
! * Recup diag file
!
itnum = 1
itval = 0
yfile =yrun
!
! * Store output in diag file
!
yvname = 'SURF_STO'
zwrite = tpdg%TDIAG_RUN%XSURF_STO / pruntime
 CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!                
yvname = 'QDIS'
zwrite = tpdg%TDIAG_RUN%XQDIS / pruntime 
 CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
IF(ldiag_misc)THEN                
  yvname = 'QSIN'
  zwrite = tpdg%TDIAG_RUN%XQIN / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval) 
ENDIF
!
IF(lcpl_land.AND.ldiag_misc)THEN                
  yvname = 'RUNOFF'
  zwrite = tpdg%TDIAG_RUN%XRUNOFF !mm of water
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)
  yvname = 'DRAIN'
  zwrite = tpdg%TDIAG_RUN%XDRAIN !mm of water
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)
ENDIF
!
IF(cgroundw/='DEF')THEN
!
  yvname = 'QGF'
  zwrite = tpdg%TDIAG_RUN%XQGF / pruntime 
  CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)
!
  yvname = 'GROUND_STO'
  zwrite = tpdg%TDIAG_RUN%XGROUND_STO / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)        
!
ENDIF
!
IF(cgroundw=='DIF')THEN
!
  yvname = 'HGROUND'
  zwrite = tpdg%TDIAG_RUN%XHGROUND / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)
!
  yvname = 'FWTD'
  zwrite = tpdg%TDIAG_RUN%XFWTD / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)
!
  yvname = 'WTD'
  zwrite = tpdg%TDIAG_RUN%XWTD / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)    
!
  IF(ldiag_misc)THEN                
!
    yvname = 'QGCELL'
    zwrite = tpdg%TDIAG_RUN%XQGCELL / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)
!
    yvname = 'HGHRIV'
    zwrite = tpdg%TDIAG_RUN%XHGHS / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask_gw,zwrite,itnum,itval)
!
  ENDIF
!
ENDIF
!
IF(cvit=='VAR')THEN
!
  yvname = 'VEL'
  zwrite = tpdg%TDIAG_RUN%XVEL / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
  yvname = 'HSTREAM'
  zwrite = tpdg%TDIAG_RUN%XHS / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
ENDIF
!
IF(lflood)THEN
!
  yvname = 'FFLOOD'
  zwrite = tpdg%TDIAG_RUN%XFF / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
  yvname = 'FLOOD_STO'
  zwrite = tpdg%TDIAG_RUN%XFLOOD_STO / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
  yvname = 'HFLOOD'
  zwrite = tpdg%TDIAG_RUN%XHF / pruntime
  CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
  IF(ldiag_misc)THEN
!
    yvname = 'FSOURCE'
    zwrite = tpdg%TDIAG_RUN%XSOURCE !mm of water
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)
!
    yvname = 'QFR'
    zwrite = tpdg%TDIAG_RUN%XQFR / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!                
    yvname = 'QRF'
    zwrite = tpdg%TDIAG_RUN%XQRF / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
    yvname = 'VFIN'
    zwrite = tpdg%TDIAG_RUN%XVFIN / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
    yvname = 'VFOUT'
    zwrite = tpdg%TDIAG_RUN%XVFOUT / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
    yvname = 'HSF'
    zwrite = tpdg%TDIAG_RUN%XHSF / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
    yvname = 'WF'
    zwrite = tpdg%TDIAG_RUN%XWF / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)        
!
    yvname = 'LF'
    zwrite = tpdg%TDIAG_RUN%XLF / pruntime
    CALL write_trip(klisting,yfile,yvname,lmask,zwrite,itnum,itval)       
!
  ENDIF       
!
ENDIF
!
IF (lhook) CALL dr_hook('TRIP_DIAG_RUN',1,zhook_handle)
!
!-------------------------------------------------------------------------------
END SUBROUTINE trip_diag_run