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