read_nam_pgd_orography.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 read_nam_pgd_orography(HPROGRAM, HZS, HFILETYPE, PUNIF_ZS, &
                                          HOROGTYPE, PENV, OIMP_ZS,&
                                  HSLOPE, HSLOPEFILETYPE, OEXPLICIT_SLOPE )  
!     ##############################################################
!
!!**** *READ_NAM_PGD_OROGRAPHY* reads namelist for Orography
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    B. Decharme        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    02/2010
!!    M Lafaysse 07/2013 : explicit slope
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_surf_par,       ONLY : xundef
!
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
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!                                   
 CHARACTER(LEN=6),    INTENT(IN)    :: HPROGRAM    ! Type of program
 CHARACTER(LEN=28),   INTENT(OUT)   :: HZS         ! file name for orography
 CHARACTER(LEN=6),    INTENT(OUT)   :: HFILETYPE   ! data file type
REAL,                INTENT(OUT)   :: PUNIF_ZS    ! uniform orography
 CHARACTER(LEN=3),    INTENT(OUT)   :: HOROGTYPE   ! orogpraphy type 
REAL,                INTENT(OUT)   :: PENV        ! parameter for enveloppe orography:
LOGICAL,             INTENT(OUT)   :: OIMP_ZS     ! Imposed orography from another PGD file
 CHARACTER(LEN=28),   INTENT(OUT),OPTIONAL   :: HSLOPE         ! file name for slope
 CHARACTER(LEN=6),    INTENT(OUT),OPTIONAL   :: HSLOPEFILETYPE   ! data file type
LOGICAL,             INTENT(OUT),OPTIONAL   :: OEXPLICIT_SLOPE ! Slope is computed from explicit ZS field and not subgrid orography
!                                  
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
INTEGER                           :: ILUOUT    ! output listing logical unit
INTEGER                           :: ILUNAM    ! namelist file logical unit
LOGICAL                           :: GFOUND    ! flag when namelist is present
!
!*    0.3    Declaration of namelists
!            ------------------------
!
 CHARACTER(LEN=28)        :: YZS         ! file name for orography
 CHARACTER(LEN=6)         :: YZSFILETYPE   ! data file type
 CHARACTER(LEN=28)        :: YSLOPE         ! file name for slope
 CHARACTER(LEN=6)         :: YSLOPEFILETYPE   ! data file type
REAL                     :: XUNIF_ZS    ! uniform orography
 CHARACTER(LEN=3)         :: COROGTYPE   ! orogpraphy type 
!                                       ! 'AVG' : average orography
!                                       ! 'SIL' : silhouette orography
!                                       ! 'ENV' : enveloppe orography
REAL                     :: XENV        ! parameter for enveloppe orography:
!                                       ! zs = avg_zs + XENV * SSO_STEDV
LOGICAL                  :: LIMP_ZS     ! Imposed orography from another PGD file
LOGICAL                  :: LEXPLICIT_SLOPE ! Slope is computed from explicit ZS field and not subgrid orography
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
NAMELIST/nam_zs/yzs, yzsfiletype, xunif_zs, corogtype, xenv, limp_zs , & 
                yslope, yslopefiletype, lexplicit_slope
!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations of defaults
!             ---------------------------
!
IF (lhook) CALL dr_hook('READ_NAM_PGD_OROGRAPHY',0,zhook_handle)
xunif_zs       = xundef
yzs            = '                          '
yzsfiletype    = '      '
yslope            = '                          '
yslopefiletype      = '      '
!
COROGTYPE      = 'ENV'
xenv           = 0.
limp_zs        = .false.
lexplicit_slope=.false.
!
 CALL get_luout(hprogram,iluout)
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of namelist
!             -------------------
!
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_ZS',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_zs)
!
 CALL close_namelist(hprogram,ilunam)
!
!-------------------------------------------------------------------------------
!
hzs       = yzs       ! file name for orography
hfiletype = yzsfiletype ! data file type
IF (PRESENT(hslope)) THEN
  hslope       = yslope       ! file name for slope
  hslopefiletype = yslopefiletype ! data file type
END IF
punif_zs  = xunif_zs  ! uniform orography
horogtype = corogtype ! orogpraphy type 
penv      = xenv      ! parameter for enveloppe orography:
oimp_zs   = limp_zs   ! Imposed orography from another PGD file
IF (PRESENT(oexplicit_slope)) THEN
    oexplicit_slope=lexplicit_slope ! Slope is computed from explicit ZS field and not subgrid orography
END IF
IF (lhook) CALL dr_hook('READ_NAM_PGD_OROGRAPHY',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_nam_pgd_orography