read_nam_pgd_teb_greenroof.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_teb_greenroof(HPROGRAM, KTIME_GR,KLAYER_GR,HTYP_GR,                   &
                                            punif_om_gr, punif_clay_gr, punif_sand_gr, punif_lai_gr,&
                                            hfnam_om_gr, hfnam_clay_gr, hfnam_sand_gr, hfnam_lai_gr,&
                                            hftyp_om_gr, hftyp_clay_gr, hftyp_sand_gr, hftyp_lai_gr )
!     ##############################################################
!
!!**** *READ_NAM_PGD_TEB_GREENROOF* reading of greenroof namelist
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!    A. Lemonsu / C. de Munck        : TEB GreenRoof namelist
!!
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original                     07/2011
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_surf_par,             ONLY : xundef
USE modd_teb_veg,              ONLY : ntime_gr_max, nlayer_gr_max
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
!
USE modi_abor1_sfx
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
INTEGER,                         INTENT(OUT)   :: ktime_gr
INTEGER,                         INTENT(OUT)   :: klayer_gr
 CHARACTER(LEN=5),                INTENT(OUT)   :: htyp_gr
REAL,DIMENSION(:),               INTENT(OUT)   :: punif_om_gr
REAL,DIMENSION(:),               INTENT(OUT)   :: punif_clay_gr
REAL,DIMENSION(:),               INTENT(OUT)   :: punif_sand_gr
REAL,DIMENSION(:),               INTENT(OUT)   :: punif_lai_gr
 CHARACTER(LEN=28),DIMENSION(:),  INTENT(OUT)   :: hfnam_om_gr
 CHARACTER(LEN=28),DIMENSION(:),  INTENT(OUT)   :: hfnam_clay_gr
 CHARACTER(LEN=28),DIMENSION(:),  INTENT(OUT)   :: hfnam_sand_gr
 CHARACTER(LEN=28),DIMENSION(:),  INTENT(OUT)   :: hfnam_lai_gr
 CHARACTER(LEN=6),DIMENSION(:),   INTENT(OUT)   :: hftyp_om_gr
 CHARACTER(LEN=6),DIMENSION(:),   INTENT(OUT)   :: hftyp_clay_gr
 CHARACTER(LEN=6),DIMENSION(:),   INTENT(OUT)   :: hftyp_sand_gr
 CHARACTER(LEN=6),DIMENSION(:),   INTENT(OUT)   :: hftyp_lai_gr

!*    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
!
!*    0.3    Declaration of namelists
!            ------------------------
!
INTEGER                                    :: nlayer_gr      ! 
INTEGER                                    :: ntime_gr       ! 
 CHARACTER(LEN=5)                           :: ctyp_gr        ! type of green roof
!
! uniform value
!
REAL,DIMENSION(NLAYER_GR_MAX)              :: xunif_om_gr      ! fraction of organic matter (OM) in green roof layer
REAL,DIMENSION(NLAYER_GR_MAX)              :: xunif_clay_gr    ! fraction of clay for the non-OM part of the green roof layer
REAL,DIMENSION(NLAYER_GR_MAX)              :: xunif_sand_gr    ! fraction of sand for the non-OM part of the green roof layer
REAL,DIMENSION(NTIME_GR_MAX)               :: xunif_lai_gr     ! LAI of green roof vegetation
!
! name of files containing data
!
 CHARACTER(LEN=28),DIMENSION(NLAYER_GR_MAX) :: cfnam_om_gr      ! fraction of organic matter (OM) in green roof layer
 CHARACTER(LEN=28),DIMENSION(NLAYER_GR_MAX) :: cfnam_clay_gr    ! fraction of clay for the non-OM part of the green roof layer
 CHARACTER(LEN=28),DIMENSION(NLAYER_GR_MAX) :: cfnam_sand_gr    ! fraction of sand for the non-OM part of the green roof layer
 CHARACTER(LEN=28),DIMENSION(NTIME_GR_MAX)  :: cfnam_lai_gr     ! LAI  of green roof
!
! type of files containing data
!
 CHARACTER(LEN=6 ),DIMENSION(NLAYER_GR_MAX) :: cftyp_om_gr      ! fraction of organic matter (OM) in green roof layer
 CHARACTER(LEN=6 ),DIMENSION(NLAYER_GR_MAX) :: cftyp_clay_gr    ! fraction of clay for the non-OM part of the green roof layer
 CHARACTER(LEN=6 ),DIMENSION(NLAYER_GR_MAX) :: cftyp_sand_gr    ! fraction of sand for the non-OM part of the green roof layer
 CHARACTER(LEN=6 ),DIMENSION(NTIME_GR_MAX)  :: cftyp_lai_gr     ! LAI  of green roof
!
REAL(KIND=JPRB) :: zhook_handle
!
namelist/nam_data_teb_greenroof/ ntime_gr,nlayer_gr,                                      &
                                 ctyp_gr,                                                 & 
                                 xunif_om_gr, xunif_clay_gr, xunif_sand_gr, xunif_lai_gr, &
                                 cfnam_om_gr, cfnam_clay_gr, cfnam_sand_gr, cfnam_lai_gr, &
                                 cftyp_om_gr, cftyp_clay_gr, cftyp_sand_gr, cftyp_lai_gr
!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations
!             ---------------
!
IF (lhook) CALL dr_hook('PGD_TEB_GREENROOF_PAR',0,zhook_handle)
!
ntime_gr         = 1
nlayer_gr        = 6
 ctyp_gr          = 'GRASS'           ! Grasses - graminoïds
!
xunif_om_gr      = xundef
xunif_clay_gr    = xundef
xunif_sand_gr    = xundef
xunif_lai_gr     = xundef
!
 cfnam_om_gr      = '                            '
 cfnam_clay_gr    = '                            '
 cfnam_sand_gr    = '                            '
 cfnam_lai_gr     = '                            '
!
 cftyp_om_gr      = '      '
 cftyp_clay_gr    = '      '
 cftyp_sand_gr    = '      '
 cftyp_lai_gr     = '      '
!
!-------------------------------------------------------------------------------
!
!*    2.      Reading of green roof namelist
!             ------------------------------
!
 CALL get_luout(hprogram,iluout)
 CALL open_namelist(hprogram,ilunam)
 CALL posnam(ilunam,'NAM_DATA_TEB_GREENROOF',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_data_teb_greenroof)
!
 CALL close_namelist(hprogram,ilunam)
!
!--------------------------------------------------
IF (ntime_gr/=1 .AND. ntime_gr/=12) THEN
  CALL abor1_sfx('NTIME_GR must be either equal to 1 (uniform LAI) or 12 (monthly LAI)')
END IF
!--------------------------------------------------
!
ktime_gr           = ntime_gr
klayer_gr          = nlayer_gr
htyp_gr            = ctyp_gr
punif_om_gr        = xunif_om_gr  
punif_clay_gr      = xunif_clay_gr
punif_sand_gr      = xunif_sand_gr
punif_lai_gr       = xunif_lai_gr 
hfnam_om_gr        = cfnam_om_gr  
hfnam_clay_gr      = cfnam_clay_gr
hfnam_sand_gr      = cfnam_sand_gr
hfnam_lai_gr       = cfnam_lai_gr 
hftyp_om_gr        = cftyp_om_gr  
hftyp_clay_gr      = cftyp_clay_gr
hftyp_sand_gr      = cftyp_sand_gr
hftyp_lai_gr       = cftyp_lai_gr 
!
IF (lhook) CALL dr_hook('READ_NAM_PGD_TEB_GREENROOF',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE read_nam_pgd_teb_greenroof