init_restart_trip.F90

Go to the documentation of this file.

!     #########
      SUBROUTINE init_restart_trip (TPG, &
                                     KLISTING,HFILE,KLON,KLAT,HTITLE,HTIMEUNIT,OTIME)
!     #######################################################################
!
!!****  *INIT_RESTART_TRIP*  
!!
!!    PURPOSE
!!    -------
!
!     Define the name and unit of each trip restart variables.
!     
!!    REFERENCE
!!    ---------
!!      
!!    AUTHOR
!!    ------
!!      B. Decharme
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    27/05/08 
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
!
!
USE modd_trip_grid, ONLY : trip_grid_t
!
USE mode_trip_netcdf
!
USE modn_trip,     ONLY : cgroundw, lflood
USE modd_trip_par, ONLY : xundef, lncprint
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE netcdf
!
USE modi_get_lonlat_trip
!
IMPLICIT NONE
!
!
!*      0.1    declarations of arguments
!
!
!
TYPE(trip_grid_t), INTENT(INOUT) :: TPG
!
 CHARACTER(LEN=*), INTENT(IN) :: HFILE, HTITLE, HTIMEUNIT
!
INTEGER, INTENT(IN)          :: KLISTING, KLON, KLAT
!
LOGICAL, INTENT(IN)          :: OTIME
!
!*      0.2    declarations of restart variables
!
 CHARACTER(LEN=NF90_MAX_NAME), DIMENSION(:), ALLOCATABLE :: YVNAME  !Name of each restart variable
 CHARACTER(LEN=NF90_MAX_NAME), DIMENSION(:), ALLOCATABLE :: YVLNAME !Long name of each restart variables
 CHARACTER(LEN=NF90_MAX_NAME), DIMENSION(:), ALLOCATABLE :: YUNIT   !Unit of each restart variable
!
 CHARACTER(LEN=NF90_MAX_NAME) :: YFILE,YTITLE,YTIMEUNIT
!
REAL,    DIMENSION(:), ALLOCATABLE ::  ZLON
REAL,    DIMENSION(:), ALLOCATABLE ::  ZLAT
LOGICAL, DIMENSION(:), ALLOCATABLE ::  LDOUBLE
!
INTEGER :: IND, INCID, INUM
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
!
! * Number of restart variable
!
IF (lhook) CALL dr_hook('INIT_RESTART_TRIP',0,zhook_handle)
ind = 1
IF(cgroundw=='CST'.OR.cgroundw=='DIF') ind = ind + 1
IF(lflood)ind = ind + 3
!
! * Allocate netcdf file attributs
!
ALLOCATE(yvname(ind))
ALLOCATE(yvlname(ind))
ALLOCATE(yunit(ind))
ALLOCATE(ldouble(ind))
ldouble(:)=.true.
!
ALLOCATE(zlon(klon))
ALLOCATE(zlat(klat))
!
! * Initialyse netcdf file attributs
!
yvname(1) = 'SURF_STO                  '
yvlname(1) = 'River storage             '
yunit(1) = 'kg                        '
!
inum = 1
!
IF(cgroundw=='CST')THEN
!        
inum = inum + 1
yvname(inum) = 'GROUND_STO                '
yvlname(inum) = 'Groundwater storage       '
yunit(inum) = 'kg                        '
!
ELSEIF(cgroundw=='DIF')THEN
!        
inum = inum + 1
yvname(inum) = 'HGROUND                   '
yvlname(inum) = 'Groundwater height        '
yunit(inum) = 'm                         '
!
ENDIF
!
IF(lflood)THEN
!
inum = inum + 1
yvname(inum) = 'FLOOD_STO                 '
yvlname(inum) = 'Floodplain storage        '
yunit(inum) = 'kg                        '

inum = inum + 1
yvname(inum) = 'FFLOOD                    '
yvlname(inum) = 'TRIP flooded fraction     '
yunit(inum) = '-                         '
!
inum = inum + 1
yvname(inum) = 'HFLOOD                    '
yvlname(inum) = 'Flood depth               '
yunit(inum) = 'm                         '
!
ENDIF
!
! * Create netcdf file
!
yfile     = hfile(1:len_trim(hfile))
ytitle    = htitle(1:len_trim(htitle))
ytimeunit = htimeunit(1:len_trim(htimeunit))
!
 CALL get_lonlat_trip(tpg, &
                     klon,klat,zlon,zlat)
!
 CALL nccreate(klisting,yfile,ytitle,ytimeunit,yvname,yvlname,yunit,zlon,zlat, &
              xundef,lncprint,incid,otime,odouble=ldouble)
!
 CALL ncclose(klisting,lncprint,yfile,incid)
!
! * Deallocate netcdf file attributs
!
DEALLOCATE(yvname  )
DEALLOCATE(yvlname )
DEALLOCATE(yunit   )
DEALLOCATE(zlon    )
DEALLOCATE(zlat    )
IF (lhook) CALL dr_hook('INIT_RESTART_TRIP',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!-------------------------------------------------------------------------------
END SUBROUTINE init_restart_trip