pgd_tsz0_par.F90

Go to the documentation of this file.

!SFX_LIC Copyright 1994-2014 CNRS, Meteo-France and Universite Paul Sabatier
!SFX_LIC This is part of the SURFEX software governed by the CeCILL-C licence
!SFX_LIC version 1. See LICENSE, CeCILL-C_V1-en.txt and CeCILL-C_V1-fr.txt  
!SFX_LIC for details. version 1.
!     #########
      SUBROUTINE pgd_tsz0_par (DTZ, &
                               hprogram)
!     ##############################################################
!
!!**** *PGD_TSZ0_PAR* monitor for averaging and interpolations of sst
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    P. Le Moigne        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original     09/2007
!!    P. Le Moigne 03/2015 tsz0 time management
!!
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE modd_data_tsz0_n, ONLY : data_tsz0_t
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
!
USE mode_pos_surf
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(data_tsz0_t), INTENT(INOUT) :: dtz
!
 CHARACTER(LEN=6),    INTENT(IN)    :: hprogram     ! Type of program
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER               :: iluout    ! output listing logical unit
INTEGER               :: ilunam    ! namelist file  logical unit
LOGICAL               :: gfound    ! true if namelist is found
!
INTEGER               :: jtime     ! loop counter on time
!
!*    0.3    Declaration of namelists
!            ------------------------
!
INTEGER            :: ntime
INTEGER, PARAMETER :: ntime_max    = 37
!
REAL, DIMENSION(NTIME_MAX)     :: xunif_dts
REAL, DIMENSION(NTIME_MAX)     :: xunif_dhugrd
 CHARACTER(LEN=28), DIMENSION(NTIME_MAX)  :: cfnam_dts
 CHARACTER(LEN=28), DIMENSION(NTIME_MAX)  :: cfnam_dhugrd
 CHARACTER(LEN=6), DIMENSION(NTIME_MAX)   :: cftyp_dts
 CHARACTER(LEN=6), DIMENSION(NTIME_MAX)   :: cftyp_dhugrd
!
! name of files containing data
!
REAL(KIND=JPRB) :: zhook_handle
!
namelist/nam_data_tsz0/ntime, xunif_dts, xunif_dhugrd  
!-------------------------------------------------------------------------------
!
!*    1.      Initializations
!             ---------------
!
IF (lhook) CALL dr_hook('PGD_TSZ0_PAR',0,zhook_handle)
ntime             = 25
xunif_dts(:)     = -0.250
xunif_dhugrd(:)   = 0.0
 cfnam_dts(:) = '                            '
 cfnam_dhugrd(:) = '                            '
 cftyp_dts(:) = '      '
 cftyp_dhugrd(:) = '      '
!
!-------------------------------------------------------------------------------
!
!*    2.      Input file for cover types
!             --------------------------
!
 CALL get_luout(hprogram,iluout)
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_DATA_TSZ0',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_data_tsz0)
!
 CALL close_namelist(hprogram,ilunam)
!
IF (ntime > ntime_max) THEN
   WRITE(iluout,*)'NTIME SHOULD NOT EXCEED',ntime_max
   CALL abor1_sfx('PGD_TSZ0_PAR: NTIME TOO BIG')
ENDIF
!
ALLOCATE(dtz%XDATA_DTS    (ntime))
ALLOCATE(dtz%XDATA_DHUGRD (ntime))
!
!-------------------------------------------------------------------------------
!
!*    3.      Uniform fields are prescribed
!             -----------------------------
!
IF (ntime==1) THEN
  dtz%XDATA_DTS   (:) = xunif_dts(1)
  dtz%XDATA_DHUGRD(:) = xunif_dhugrd(1)
ELSE
  DO jtime=1,ntime
    dtz%XDATA_DTS   (jtime) = xunif_dts(jtime)
    dtz%XDATA_DHUGRD(jtime) = xunif_dhugrd(jtime)
  END DO
ENDIF
!
IF (lhook) CALL dr_hook('PGD_TSZ0_PAR',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE pgd_tsz0_par