pgd_teb_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_teb_par (DTCO, DGU, UG, U, USS, BDD, DTT, DTI, TG, &
                              hprogram,ogarden,ogreenroof,hbld_atype)
!     ##############################################################
!
!!**** *PGD_TEB_PAR* monitor for averaging and interpolations of cover fractions
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    10/12/97
!!
!!       Modified 08/12/05, P. Le Moigne: user defined fields
!!    G. Pigeon      09/2012: add ROUGH_WALL/ROUGH_ROOF for outdoor convection
!!    V. Masson      08/2013: adds solar panels
!!    V. Masson      10/2013: adds residential fraction
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
!
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_surf_atm_sso_n, ONLY : surf_atm_sso_t
USE modd_bld_description_n, ONLY : bld_desc_t
USE modd_data_teb_n, ONLY : data_teb_t
USE modd_data_isba_n, ONLY : data_isba_t
USE modd_teb_grid_n, ONLY : teb_grid_t
!
USE modd_surf_par,   ONLY : xundef, nundef
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
USE modi_ini_var_from_data_0d
USE modi_ini_var_from_data
USE modi_test_nam_var_surf
USE modi_read_csvdata_teb
USE modi_bldcode
!
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_cover_t), INTENT(INOUT) :: dtco
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(surf_atm_grid_t), INTENT(INOUT) :: ug
TYPE(surf_atm_t), INTENT(INOUT) :: u
TYPE(surf_atm_sso_t), INTENT(INOUT) :: uss
TYPE(bld_desc_t), INTENT(INOUT) :: bdd
TYPE(data_teb_t), INTENT(INOUT) :: dtt
TYPE(data_isba_t), INTENT(INOUT) :: dti
TYPE(teb_grid_t), INTENT(INOUT) :: tg
!
 CHARACTER(LEN=6),    INTENT(IN)    :: hprogram     ! Type of program
LOGICAL,             INTENT(IN)    :: ogarden      ! T if urban green areas
LOGICAL,             INTENT(IN)    :: ogreenroof   ! T if greenroofs option is activated
 CHARACTER(LEN=3),    INTENT(OUT)   :: hbld_atype    ! Type of building averaging
!
!
!*    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
!
REAL, DIMENSION(TG%NDIM) :: zwork
REAL                  :: zunif     ! temporary variable
!
!*    0.3    Declaration of namelists
!            ------------------------
!
INTEGER, PARAMETER :: nroof_max  = 9
INTEGER, PARAMETER :: nroad_max  = 9
INTEGER, PARAMETER :: nwall_max  = 9
INTEGER            :: npar_road_layer ! number of road layers
INTEGER            :: npar_roof_layer ! number of roof layers
INTEGER            :: npar_wall_layer ! number of wall layers
!
! Geometric Parameters:
!
INTEGER                                 :: nunif_bldtype
 CHARACTER(LEN=28)                       :: cfnam_bldtype
 CHARACTER(LEN=6)                        :: cftyp_bldtype
INTEGER                                 :: nunif_bld_age
 CHARACTER(LEN=28)                       :: cfnam_bld_age
 CHARACTER(LEN=6)                        :: cftyp_bld_age
 CHARACTER(LEN=28)                       :: ccsvdatafile
INTEGER                                 :: nunif_usetype
 CHARACTER(LEN=28)                       :: cfnam_usetype
 CHARACTER(LEN=6)                        :: cftyp_usetype
 CHARACTER(LEN=3)                        :: cbld_atype         ! type of averaging for buildings

!
REAL                                    :: xunif_bld          ! fraction of buildings            (-)
REAL                                    :: xunif_bld_height   ! buildings height 'h'             (m)
REAL                                    :: xunif_wall_o_hor   ! wall surf. / hor. surf.          (-)
REAL                                    :: xunif_z0_town      ! roughness length for momentum    (m)
REAL                                    :: xunif_garden       ! fraction of veg in the streets   (-)
REAL                                    :: xunif_greenroof    ! fraction of greenroofs on roofs  (-)
REAL                                    :: xunif_road_dir     ! road direction (deg from North, clockwise)
 CHARACTER(LEN=28)                       :: cfnam_bld          ! file name for BLD 
 CHARACTER(LEN=28)                       :: cfnam_bld_height   ! file name for BLD_HEIGHT
 CHARACTER(LEN=28)                       :: cfnam_wall_o_hor   ! file name for WALL_O_HOR
 CHARACTER(LEN=28)                       :: cfnam_z0_town      ! file name for Z0_TOWN
 CHARACTER(LEN=28)                       :: cfnam_garden       ! file name for GARDEN  
 CHARACTER(LEN=28)                       :: cfnam_greenroof    ! file name for GREENROOF
 CHARACTER(LEN=28)                       :: cfnam_road_dir     ! file name for ROAD_DIR  
 CHARACTER(LEN=6)                        :: cftyp_bld          ! file type for BLD 
 CHARACTER(LEN=6)                        :: cftyp_bld_height   ! file type for BLD_HEIGHT
 CHARACTER(LEN=6)                        :: cftyp_wall_o_hor   ! file type for WALL_O_HOR
 CHARACTER(LEN=6)                        :: cftyp_z0_town      ! file type for Z0_TOWN
 CHARACTER(LEN=6)                        :: cftyp_garden       ! file type for GARDEN  
 CHARACTER(LEN=6)                        :: cftyp_greenroof    ! file type for GREENROOF
 CHARACTER(LEN=6)                        :: cftyp_road_dir     ! file type for ROAD_DIR
!
! Roof parameters
!
REAL                                    :: xunif_alb_roof     ! roof albedo                      (-)
REAL                                    :: xunif_emis_roof    ! roof emissivity                  (-)
 CHARACTER(LEN=28)                       :: cfnam_alb_roof     ! file name for ALB_ROOF
 CHARACTER(LEN=28)                       :: cfnam_emis_roof    ! file name for EMIS_ROOF
 CHARACTER(LEN=6)                        :: cftyp_alb_roof     ! file name for ALB_ROOF   
 CHARACTER(LEN=6)                        :: cftyp_emis_roof    ! file name for EMIS_ROOF  
REAL, DIMENSION(NROOF_MAX)              :: xunif_hc_roof      ! roof layers heat capacity        (J/K/m3)
REAL, DIMENSION(NROOF_MAX)              :: xunif_tc_roof      ! roof layers thermal conductivity (W/K/m)
REAL, DIMENSION(NROOF_MAX)              :: xunif_d_roof       ! depth of roof layers             (m)
 CHARACTER(LEN=28), DIMENSION(NROOF_MAX) :: cfnam_hc_roof      ! file name for HC_ROOF   
 CHARACTER(LEN=28), DIMENSION(NROOF_MAX) :: cfnam_tc_roof      ! file name for TC_ROOF
 CHARACTER(LEN=28), DIMENSION(NROOF_MAX) :: cfnam_d_roof       ! file name for D_ROOF
 CHARACTER(LEN=6),  DIMENSION(NROOF_MAX) :: cftyp_hc_roof      ! file type for HC_ROOF   
 CHARACTER(LEN=6),  DIMENSION(NROOF_MAX) :: cftyp_tc_roof      ! file type for TC_ROOF
 CHARACTER(LEN=6),  DIMENSION(NROOF_MAX) :: cftyp_d_roof       ! file type for D_ROOF
REAL                                    :: xunif_rough_roof  ! roof roughness coef
 CHARACTER(LEN=28)                       :: cfnam_rough_roof  ! file name for ROUGH_ROOF
 CHARACTER(LEN=6)                        :: cftyp_rough_roof  ! file type for ROUGH_ROOF
!
!
! Road parameters
!
REAL                                    :: xunif_alb_road     ! road albedo                      (-)
REAL                                    :: xunif_emis_road    ! road emissivity                  (-)
 CHARACTER(LEN=28)                       :: cfnam_alb_road     ! file name for ALB_ROAD
 CHARACTER(LEN=28)                       :: cfnam_emis_road    ! file name for EMIS_ROAD
 CHARACTER(LEN=6)                        :: cftyp_alb_road     ! file type for ALB_ROAD
 CHARACTER(LEN=6)                        :: cftyp_emis_road    ! file type for EMIS_ROAD
REAL, DIMENSION(NROAD_MAX)              :: xunif_hc_road      ! road layers heat capacity        (J/K/m3)
REAL, DIMENSION(NROAD_MAX)              :: xunif_tc_road      ! road layers thermal conductivity (W/K/m)
REAL, DIMENSION(NROAD_MAX)              :: xunif_d_road       ! depth of road layers             (m)
 CHARACTER(LEN=28), DIMENSION(NROAD_MAX) :: cfnam_hc_road      ! file name for HC_ROAD   
 CHARACTER(LEN=28), DIMENSION(NROAD_MAX) :: cfnam_tc_road      ! file name for TC_ROAD
 CHARACTER(LEN=28), DIMENSION(NROAD_MAX) :: cfnam_d_road       ! file name for D_ROAD
 CHARACTER(LEN=6),  DIMENSION(NROAD_MAX) :: cftyp_hc_road      ! file type for HC_ROAD   
 CHARACTER(LEN=6),  DIMENSION(NROAD_MAX) :: cftyp_tc_road      ! file type for TC_ROAD
 CHARACTER(LEN=6),  DIMENSION(NROAD_MAX) :: cftyp_d_road       ! file type for D_ROAD
!
! Wall parameters
!
REAL                                    :: xunif_alb_wall     ! wall albedo                      (-)
REAL                                    :: xunif_emis_wall    ! wall emissivity                  (-)
 CHARACTER(LEN=28)                       :: cfnam_alb_wall     ! file name for ALB_WALL
 CHARACTER(LEN=28)                       :: cfnam_emis_wall    ! file name for EMIS_WALL
 CHARACTER(LEN=6)                        :: cftyp_alb_wall     ! file type for ALB_WALL
 CHARACTER(LEN=6)                        :: cftyp_emis_wall    ! file type for EMIS_WALL
REAL, DIMENSION(NWALL_MAX)              :: xunif_hc_wall      ! wall layers heat capacity        (J/K/m3)
REAL, DIMENSION(NWALL_MAX)              :: xunif_tc_wall      ! wall layers thermal conductivity (W/K/m)
REAL, DIMENSION(NWALL_MAX)              :: xunif_d_wall       ! depth of wall layers             (m)
 CHARACTER(LEN=28), DIMENSION(NWALL_MAX) :: cfnam_hc_wall      ! file name for HC_WALL   
 CHARACTER(LEN=28), DIMENSION(NWALL_MAX) :: cfnam_tc_wall      ! file name for TC_WALL
 CHARACTER(LEN=28), DIMENSION(NWALL_MAX) :: cfnam_d_wall       ! file name for D_WALL
 CHARACTER(LEN=6),  DIMENSION(NWALL_MAX) :: cftyp_hc_wall      ! file type for HC_WALL   
 CHARACTER(LEN=6),  DIMENSION(NWALL_MAX) :: cftyp_tc_wall      ! file type for TC_WALL
 CHARACTER(LEN=6),  DIMENSION(NWALL_MAX) :: cftyp_d_wall       ! file type for D_WALL
REAL                                    :: xunif_rough_wall  ! wall roughness coef
 CHARACTER(LEN=28)                       :: cfnam_rough_wall  ! file name for ROUGH_WALL
 CHARACTER(LEN=6)                        :: cftyp_rough_wall  ! file type for ROUGH_WALL
REAL                                    :: xunif_residential ! residential fraction
 CHARACTER(LEN=28)                       :: cfnam_residential ! file name for RESIDENTIAL
 CHARACTER(LEN=6)                        :: cftyp_residential ! file type for RESIDENTIAL
!
! anthropogenic fluxes
!
REAL                                    :: xunif_h_traffic    ! anthropogenic sensible
!                                                             ! heat fluxes due to traffic       (W/m2)
REAL                                    :: xunif_le_traffic   ! anthropogenic latent
!                                                             ! heat fluxes due to traffic       (W/m2)
REAL                                    :: xunif_h_industry   ! anthropogenic sensible                   
!                                                             ! heat fluxes due to factories     (W/m2)
REAL                                    :: xunif_le_industry  ! anthropogenic latent
!                                                             ! heat fluxes due to factories     (W/m2)
 CHARACTER(LEN=28)                       :: cfnam_h_traffic    ! file name for H_TRAFFIC
 CHARACTER(LEN=28)                       :: cfnam_le_traffic   ! file name for LE_TRAFFIC
 CHARACTER(LEN=28)                       :: cfnam_h_industry   ! file name for H_INDUSTRY
 CHARACTER(LEN=28)                       :: cfnam_le_industry  ! file name for LE_INDUSTRY
 CHARACTER(LEN=6)                        :: cftyp_h_traffic    ! file type for H_TRAFFIC
 CHARACTER(LEN=6)                        :: cftyp_le_traffic   ! file type for LE_TRAFFIC
 CHARACTER(LEN=6)                        :: cftyp_h_industry   ! file type for H_INDUSTRY
 CHARACTER(LEN=6)                        :: cftyp_le_industry  ! file type for LE_INDUSTRY
!
! Solar panels parameters
!
REAL                                    :: xunif_emis_panel    ! emissivity of solar panel       (-)
REAL                                    :: xunif_alb_panel     ! albedo     of solar panel       (-)
REAL                                    :: xunif_eff_panel     ! efficiency of solar panel       (-)
REAL                                    :: xunif_frac_panel    ! fraction   of solar panel       (-)
 CHARACTER(LEN=28)                       :: cfnam_emis_panel   ! file name for EMIS_PANEL
 CHARACTER(LEN=28)                       :: cfnam_alb_panel    ! file name for ALB_PANEL
 CHARACTER(LEN=28)                       :: cfnam_eff_panel    ! file name for EFF_PANEL
 CHARACTER(LEN=28)                       :: cfnam_frac_panel   ! file name for FRAC_PANEL
 CHARACTER(LEN=6)                        :: cftyp_emis_panel   ! file type for EMIS_PANEL
 CHARACTER(LEN=6)                        :: cftyp_alb_panel    ! file type for ALB_PANEL
 CHARACTER(LEN=6)                        :: cftyp_eff_panel    ! file type for EFF_PANEL
 CHARACTER(LEN=6)                        :: cftyp_frac_panel   ! file type for FRAC_PANEL

REAL(KIND=JPRB) :: zhook_handle
!

namelist/nam_data_teb/      npar_roof_layer, npar_road_layer, npar_wall_layer,&
                              cbld_atype,                                     &
                              nunif_bldtype, cfnam_bldtype, cftyp_bldtype,    &
                              nunif_bld_age, cfnam_bld_age, cftyp_bld_age,    &
                              ccsvdatafile,                                   &
                              nunif_usetype, cfnam_usetype, cftyp_usetype,    &
                              xunif_alb_roof,                                 &
                              xunif_emis_roof, xunif_hc_roof, xunif_tc_roof,  &
                              xunif_d_roof, xunif_alb_road, xunif_emis_road,  &
                              xunif_hc_road, xunif_tc_road, xunif_d_road,     &
                              xunif_alb_wall, xunif_emis_wall, xunif_hc_wall, &
                              xunif_tc_wall, xunif_d_wall,                    &
                              xunif_z0_town, xunif_bld, xunif_bld_height,     &
                              xunif_wall_o_hor,                               &
                              xunif_h_traffic, xunif_le_traffic,              &
                              xunif_h_industry, xunif_le_industry,            &
                              xunif_garden, xunif_greenroof,                  &
                              xunif_road_dir,                                 &
                              cfnam_alb_roof,                                 &
                              cfnam_emis_roof, cfnam_hc_roof, cfnam_tc_roof,  &
                              cfnam_d_roof, cfnam_alb_road, cfnam_emis_road,  &
                              cfnam_hc_road, cfnam_tc_road, cfnam_d_road,     &
                              cfnam_alb_wall, cfnam_emis_wall, cfnam_hc_wall, &
                              cfnam_tc_wall, cfnam_d_wall,                    &
                              cfnam_z0_town, cfnam_bld, cfnam_bld_height,     &
                              cfnam_wall_o_hor,                               &
                              cfnam_h_traffic, cfnam_le_traffic,              &
                              cfnam_h_industry, cfnam_le_industry,            &
                              cfnam_garden, cfnam_road_dir, cfnam_greenroof,  &
                              cftyp_alb_roof,                                 &
                              cftyp_emis_roof, cftyp_hc_roof, cftyp_tc_roof,  &
                              cftyp_d_roof, cftyp_alb_road, cftyp_emis_road,  &
                              cftyp_hc_road, cftyp_tc_road, cftyp_d_road,     &
                              cftyp_alb_wall, cftyp_emis_wall, cftyp_hc_wall, &
                              cftyp_tc_wall, cftyp_d_wall,                    &
                              cftyp_z0_town, cftyp_bld, cftyp_bld_height,     &
                              cftyp_wall_o_hor,                               &
                              cftyp_h_traffic, cftyp_le_traffic,              &
                              cftyp_h_industry, cftyp_le_industry,            &
                              cftyp_garden, cftyp_road_dir, cftyp_greenroof,  &
                              xunif_rough_roof, cfnam_rough_roof, cftyp_rough_roof, &
                              xunif_rough_wall, cfnam_rough_wall, cftyp_rough_wall, &
                              xunif_residential,cfnam_residential,cftyp_residential,&
                              xunif_emis_panel, cfnam_emis_panel, cftyp_emis_panel, &
                              xunif_alb_panel,  cfnam_alb_panel,  cftyp_alb_panel,  &
                              xunif_eff_panel,  cfnam_eff_panel,  cftyp_eff_panel,  &
                              xunif_frac_panel, cfnam_frac_panel, cftyp_frac_panel

!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations
!             ---------------
!
IF (lhook)   CALL dr_hook('PGD_TEB_PAR',0,zhook_handle)
npar_roof_layer=0
npar_road_layer=0
npar_wall_layer=0
 cbld_atype ='MAJ'
nunif_bldtype      = nundef
nunif_bld_age      = nundef
nunif_usetype      = nundef
xunif_bld          = xundef
xunif_bld_height   = xundef
xunif_wall_o_hor   = xundef
xunif_z0_town      = xundef
xunif_alb_roof     = xundef
xunif_emis_roof    = xundef
xunif_hc_roof      = xundef
xunif_tc_roof      = xundef
xunif_d_roof       = xundef
xunif_alb_road     = xundef
xunif_emis_road    = xundef
xunif_hc_road      = xundef
xunif_tc_road      = xundef
xunif_d_road       = xundef
xunif_alb_wall     = xundef
xunif_emis_wall    = xundef
xunif_hc_wall      = xundef
xunif_tc_wall      = xundef
xunif_d_wall       = xundef
xunif_h_traffic    = xundef
xunif_le_traffic   = xundef
xunif_h_industry   = xundef
xunif_le_industry  = xundef
xunif_garden       = xundef
xunif_greenroof    = xundef
xunif_road_dir     = xundef
xunif_rough_roof   = xundef
xunif_rough_wall   = xundef
xunif_residential  = xundef
xunif_emis_panel   = xundef
xunif_alb_panel    = xundef
xunif_eff_panel    = xundef
xunif_frac_panel   = xundef

 cfnam_bldtype      = '                            '
 cfnam_bld_age      = '                            '
 cfnam_usetype      = '                            '
 ccsvdatafile       ='                            '
 cfnam_bld          = '                            '
 cfnam_bld_height   = '                            '
 cfnam_wall_o_hor   = '                            '
 cfnam_z0_town      = '                            '

 cfnam_alb_roof(:) = '                            '
 cfnam_emis_roof(:) = '                            '
 cfnam_hc_roof(:) = '                            '
 cfnam_tc_roof(:) = '                            '
 cfnam_d_roof(:) = '                            '
 cfnam_rough_roof(:) = '                            '
 cfnam_rough_wall(:) = '                            '
 cfnam_residential(:)= '                            '
 cfnam_alb_road(:) = '                            '
 cfnam_emis_road(:) = '                            '
 cfnam_hc_road(:) = '                            '
 cfnam_tc_road(:) = '                            '
 cfnam_d_road(:) = '                            '
 cfnam_alb_wall(:) = '                            '
 cfnam_emis_wall(:) = '                            '
 cfnam_hc_wall(:) = '                            '
 cfnam_tc_wall(:) = '                            '
 cfnam_d_wall(:) = '                            '

 cfnam_h_traffic    = '                            '
 cfnam_le_traffic   = '                            '
 cfnam_h_industry   = '                            '
 cfnam_le_industry  = '                            '

 cfnam_garden       = '                            '
 cfnam_greenroof    = '                            '
 cfnam_road_dir     = '                            '

 cfnam_emis_panel   = '                            '
 cfnam_alb_panel    = '                            '
 cfnam_eff_panel    = '                            '
 cfnam_frac_panel   = '                            '

 cftyp_bldtype      = '      '
 cftyp_bld_age      = '      '
 cftyp_usetype      = '      '
 cftyp_bld          = '      '
 cftyp_bld_height   = '      '
 cftyp_wall_o_hor   = '      '
 cftyp_z0_town      = '      '
 cftyp_alb_roof(:)  = '      '
 cftyp_emis_roof(:) = '      '
 cftyp_hc_roof(:)   = '      '
 cftyp_tc_roof(:)   = '      '
 cftyp_d_roof(:)    = '      '
 cftyp_rough_roof(:)    = '      '
 cftyp_rough_wall(:)    = '      '
 cftyp_residential(:)   = '      '
 cftyp_alb_road(:)  = '      '
 cftyp_emis_road(:) = '      '
 cftyp_hc_road(:)   = '      '
 cftyp_tc_road(:)   = '      '
 cftyp_d_road(:)    = '      '
 cftyp_alb_wall(:)  = '      '
 cftyp_emis_wall(:) = '      '
 cftyp_hc_wall(:)   = '      '
 cftyp_tc_wall(:)   = '      '
 cftyp_d_wall(:)    = '      '
 cftyp_h_traffic    = '      '
 cftyp_le_traffic   = '      '
 cftyp_h_industry   = '      '
 cftyp_le_industry  = '      '
 cftyp_garden       = '      '
 cftyp_greenroof    = '      '
 cftyp_road_dir     = '      '
!
 cftyp_emis_panel   = '      '
 cftyp_alb_panel    = '      '
 cftyp_eff_panel    = '      '
 cftyp_frac_panel   = '      '
!
!-------------------------------------------------------------------------------
!
!*    2.      Input file for cover types
!             --------------------------
!
 CALL get_luout(hprogram,iluout)
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_DATA_TEB',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_data_teb)
!
 CALL close_namelist(hprogram,ilunam)
!
 CALL test_nam_var_surf(iluout,'CBLD_ATYPE',cbld_atype,'ARI','MAJ')
!
dtt%NPAR_ROOF_LAYER = npar_roof_layer
dtt%NPAR_ROAD_LAYER = npar_road_layer
dtt%NPAR_WALL_LAYER = npar_wall_layer
!
hbld_atype = cbld_atype
!-------------------------------------------------------------------------------
!
!* coherence check
!
IF ((     any(xunif_hc_road/=xundef) .OR. any(len_trim(cfnam_hc_road)>0) &
     .OR. any(xunif_tc_road/=xundef) .OR. any(len_trim(cfnam_tc_road)>0) &
     .OR. any(xunif_d_road /=xundef) .OR. any(len_trim(cfnam_d_road )>0) &
    ) .AND. npar_road_layer<1                                  ) THEN
  CALL abor1_sfx('In order to initialize road thermal quantities, please specify NPAR_ROAD_LAYER in namelist NAM_DATA_TEB')
END IF
!
IF ((     any(xunif_hc_roof/=xundef) .OR. any(len_trim(cfnam_hc_roof)>0) &
     .OR. any(xunif_tc_roof/=xundef) .OR. any(len_trim(cfnam_tc_roof)>0) &
     .OR. any(xunif_d_roof /=xundef) .OR. any(len_trim(cfnam_d_roof )>0) &
    ) .AND. npar_roof_layer<1                                  ) THEN
  CALL abor1_sfx('In order to initialize ROOF thermal quantities, please specify NPAR_ROOF_LAYER in namelist NAM_DATA_TEB')
END IF
!
IF ((     any(xunif_hc_wall/=xundef) .OR. any(len_trim(cfnam_hc_wall)>0) &
     .OR. any(xunif_tc_wall/=xundef) .OR. any(len_trim(cfnam_tc_wall)>0) &
     .OR. any(xunif_d_wall /=xundef) .OR. any(len_trim(cfnam_d_wall )>0) &
    ) .AND. npar_wall_layer<1                                  ) THEN
  CALL abor1_sfx('In order to initialize WALL thermal quantities, please specify NPAR_WALL_LAYER in namelist NAM_DATA_TEB')
END IF
!-------------------------------------------------------------------------------
ALLOCATE(dtt%NPAR_BLDTYPE     (tg%NDIM))
ALLOCATE(dtt%NPAR_BLD_AGE     (tg%NDIM))
ALLOCATE(dtt%NPAR_BLDCODE     (tg%NDIM))
ALLOCATE(dtt%NPAR_USETYPE     (tg%NDIM))
ALLOCATE(dtt%XPAR_Z0_TOWN     (tg%NDIM))
ALLOCATE(dtt%XPAR_ALB_ROOF    (tg%NDIM))
ALLOCATE(dtt%XPAR_EMIS_ROOF   (tg%NDIM))
ALLOCATE(dtt%XPAR_ALB_ROAD    (tg%NDIM))
ALLOCATE(dtt%XPAR_EMIS_ROAD   (tg%NDIM))
ALLOCATE(dtt%XPAR_ALB_WALL    (tg%NDIM))
ALLOCATE(dtt%XPAR_EMIS_WALL   (tg%NDIM))
ALLOCATE(dtt%XPAR_BLD         (tg%NDIM))
ALLOCATE(dtt%XPAR_BLD_HEIGHT  (tg%NDIM))
ALLOCATE(dtt%XPAR_WALL_O_HOR  (tg%NDIM))
ALLOCATE(dtt%XPAR_H_TRAFFIC   (tg%NDIM))
ALLOCATE(dtt%XPAR_LE_TRAFFIC  (tg%NDIM))
ALLOCATE(dtt%XPAR_H_INDUSTRY  (tg%NDIM))
ALLOCATE(dtt%XPAR_LE_INDUSTRY (tg%NDIM))
ALLOCATE(dtt%XPAR_GARDEN      (tg%NDIM))
ALLOCATE(dtt%XPAR_GREENROOF   (tg%NDIM))
ALLOCATE(dtt%XPAR_ROAD_DIR    (tg%NDIM))
!
ALLOCATE(dtt%XPAR_HC_ROOF     (tg%NDIM,npar_roof_layer))
ALLOCATE(dtt%XPAR_TC_ROOF     (tg%NDIM,npar_roof_layer))
ALLOCATE(dtt%XPAR_D_ROOF      (tg%NDIM,npar_roof_layer))
ALLOCATE(dtt%XPAR_HC_ROAD     (tg%NDIM,npar_road_layer))
ALLOCATE(dtt%XPAR_TC_ROAD     (tg%NDIM,npar_road_layer))
ALLOCATE(dtt%XPAR_D_ROAD      (tg%NDIM,npar_road_layer))
ALLOCATE(dtt%XPAR_HC_WALL     (tg%NDIM,npar_wall_layer))
ALLOCATE(dtt%XPAR_TC_WALL     (tg%NDIM,npar_wall_layer))
ALLOCATE(dtt%XPAR_D_WALL      (tg%NDIM,npar_wall_layer))
ALLOCATE(dtt%XPAR_ROUGH_ROOF    (tg%NDIM))
ALLOCATE(dtt%XPAR_ROUGH_WALL    (tg%NDIM))
ALLOCATE(dtt%XPAR_RESIDENTIAL   (tg%NDIM))
!
ALLOCATE(dtt%XPAR_EMIS_PANEL  (tg%NDIM))
ALLOCATE(dtt%XPAR_ALB_PANEL   (tg%NDIM))
ALLOCATE(dtt%XPAR_EFF_PANEL   (tg%NDIM))
ALLOCATE(dtt%XPAR_FRAC_PANEL  (tg%NDIM))
!
!-------------------------------------------------------------------------------
IF (nroof_max < npar_roof_layer) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) 'Please update pgd_teb_par.f90 routine :      '
  WRITE(iluout,*) 'The maximum number of ROOF LAYER             '
  WRITE(iluout,*) 'in the declaration of the namelist variables '
  WRITE(iluout,*) 'must be increased to : ', npar_roof_layer
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx('PGD_TEB_PAR: MAXIMUM NUMBER OF NROOF_LAYER MUST BE INCREASED')
ENDIF
!-------------------------------------------------------------------------------
IF (nroad_max < npar_road_layer) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) 'Please update pgd_teb_par.f90 routine :      '
  WRITE(iluout,*) 'The maximum number of ROAD LAYER             '
  WRITE(iluout,*) 'in the declaration of the namelist variables '
  WRITE(iluout,*) 'must be increased to : ', npar_road_layer
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx('PGD_TEB_PAR: MAXIMUM NUMBER OF NROAD_LAYER MUST BE INCREASED')
ENDIF
!-------------------------------------------------------------------------------
IF (nwall_max < npar_wall_layer) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) 'Please update pgd_teb_par.f90 routine :      '
  WRITE(iluout,*) 'The maximum number of WALL LAYER             '
  WRITE(iluout,*) 'in the declaration of the namelist variables '
  WRITE(iluout,*) 'must be increased to : ', npar_wall_layer
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx('PGD_TEB_PAR: MAXIMUM NUMBER OF NWALL_LAYER MUST BE INCREASED')
ENDIF

!-------------------------------------------------------------------------------
!
!*    3.      user defined fields are prescribed
!             ----------------------------------
!
!* building's type
zunif = xundef
IF (nunif_bldtype/=nundef) zunif=float(nunif_bldtype)
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'MAJ','BLDTYPE    ','TWN', cfnam_bldtype,cftyp_bldtype,zunif,&
        zwork(:),dtt%LDATA_BLDTYPE )
IF (.NOT. dtt%LDATA_BLDTYPE) THEN
  DEALLOCATE(dtt%NPAR_BLDTYPE)
ELSE
  dtt%NPAR_BLDTYPE = nint(zwork)
END IF
!
!* building's age
zunif = xundef
IF (nunif_bld_age/=nundef) zunif=float(nunif_bld_age)
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','BLD_AGE    ','TWN', cfnam_bld_age,cftyp_bld_age,zunif,&
        zwork(:),dtt%LDATA_BLD_AGE )
IF (.NOT. dtt%LDATA_BLD_AGE) THEN
  DEALLOCATE(dtt%NPAR_BLD_AGE)
ELSE
  dtt%NPAR_BLD_AGE = nint(zwork)
END IF
!
IF (dtt%LDATA_BLDTYPE .AND. .NOT. dtt%LDATA_BLD_AGE) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) ' You chose to define building types :        '
  IF (nunif_bldtype/=nundef) THEN
    WRITE(iluout,*) ' NUNIF_BLDTYPE=', nunif_bldtype
  ELSE
    WRITE(iluout,*) ' CFNAM_BLDTYPE =',cfnam_bldtype 
    WRITE(iluout,*) ' CFTYP_BLDTYPE =',cftyp_bldtype 
  END IF
  WRITE(iluout,*) ' But            '
  WRITE(iluout,*) " You did not chose to define building's age"
  WRITE(iluout,*) '- - - - - - - - - - - - - - - - - - - - - - -'
  WRITE(iluout,*) ' Please define the construction date of the buildings. '
  WRITE(iluout,*) ' To do so, use either :'
  WRITE(iluout,*) ' NUNIF_BLD_AGE   (to have a uniform construction date for all buildings'
  WRITE(iluout,*) ' or CFNAM_BLD_AGE and CFTYP_BLD_AGE (to incorporate spatial data ) '
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx("PGD_TEB_PAR: Building's age data is missing")
END IF
!
!* building's use
zunif = xundef
IF (nunif_usetype/=nundef) zunif=float(nunif_usetype)
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'MAJ','USETYPE    ','TWN', cfnam_usetype,cftyp_usetype,zunif,&
        zwork(:),dtt%LDATA_USETYPE )
IF (.NOT. dtt%LDATA_USETYPE) THEN
  DEALLOCATE(dtt%NPAR_USETYPE)
ELSE
  dtt%NPAR_USETYPE = nint(zwork)
END IF
!
IF (dtt%LDATA_BLDTYPE .OR. dtt%LDATA_BLD_AGE .OR. dtt%LDATA_USETYPE)  CALL read_csvdata_teb(bdd, &
                                                                                            hprogram,ccsvdatafile)
!
!* building's code
IF (ASSOCIATED(dtt%NPAR_BLDTYPE)) dtt%NPAR_BLDCODE(:) = bldcode(bdd, dtt%NPAR_BLDTYPE,dtt%NPAR_BLD_AGE)
!
!
!* other building parameters
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','BLD        ','TWN', cfnam_bld,cftyp_bld,xunif_bld,dtt%XPAR_BLD,dtt%LDATA_BLD )
IF (.NOT.dtt%LDATA_BLD) DEALLOCATE(dtt%XPAR_BLD)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','BLD_HEIGHT ','TWN',cfnam_bld_height,cftyp_bld_height,xunif_bld_height,&
        dtt%XPAR_BLD_HEIGHT,dtt%LDATA_BLD_HEIGHT)
IF (.NOT.dtt%LDATA_BLD_HEIGHT) DEALLOCATE(dtt%XPAR_BLD_HEIGHT)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','WALL_O_HOR ','TWN',cfnam_wall_o_hor,cftyp_wall_o_hor,xunif_wall_o_hor,&
        dtt%XPAR_WALL_O_HOR,dtt%LDATA_WALL_O_HOR)
IF (.NOT.dtt%LDATA_WALL_O_HOR) DEALLOCATE(dtt%XPAR_WALL_O_HOR)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'CDN','Z0_TOWN    ','TWN',cfnam_z0_town,cftyp_z0_town,xunif_z0_town,&
        dtt%XPAR_Z0_TOWN,dtt%LDATA_Z0_TOWN)
IF (.NOT.dtt%LDATA_Z0_TOWN) DEALLOCATE(dtt%XPAR_Z0_TOWN)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'ALB_ROOF   ','TWN',cfnam_alb_roof,cftyp_alb_roof,xunif_alb_roof  ,&
        dtt%XPAR_ALB_ROOF,dtt%LDATA_ALB_ROOF)
IF (.NOT.dtt%LDATA_ALB_ROOF) DEALLOCATE(dtt%XPAR_ALB_ROOF)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'EMIS_ROOF  ','TWN',cfnam_emis_roof,cftyp_emis_roof,xunif_emis_roof ,&
        dtt%XPAR_EMIS_ROOF,dtt%LDATA_EMIS_ROOF)
IF (.NOT.dtt%LDATA_EMIS_ROOF) DEALLOCATE(dtt%XPAR_EMIS_ROOF)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'HC_ROOF  ','TWN',cfnam_hc_roof,cftyp_hc_roof, &
        xunif_hc_roof,dtt%XPAR_HC_ROOF,dtt%LDATA_HC_ROOF ) 
IF (.NOT.dtt%LDATA_HC_ROOF) DEALLOCATE(dtt%XPAR_HC_ROOF)
! 
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'TC_ROOF  ','TWN',cfnam_tc_roof,cftyp_tc_roof, &
                 xunif_tc_roof ,dtt%XPAR_TC_ROOF, dtt%LDATA_TC_ROOF ) 
IF (.NOT.dtt%LDATA_TC_ROOF) DEALLOCATE(dtt%XPAR_TC_ROOF)
! 
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'D_ROOF   ','TWN',cfnam_d_roof,cftyp_d_roof, &
                 xunif_d_roof  ,dtt%XPAR_D_ROOF , dtt%LDATA_D_ROOF ) 
IF (.NOT.dtt%LDATA_D_ROOF) DEALLOCATE(dtt%XPAR_D_ROOF)
! 
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','ALB_ROAD   ','TWN',cfnam_alb_road  ,cftyp_alb_road  ,xunif_alb_road  ,&
        dtt%XPAR_ALB_ROAD, dtt%LDATA_ALB_ROAD  )
IF (.NOT.dtt%LDATA_ALB_ROAD) DEALLOCATE(dtt%XPAR_ALB_ROAD)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','EMIS_ROAD  ','TWN',cfnam_emis_road ,cftyp_emis_road ,xunif_emis_road ,&
        dtt%XPAR_EMIS_ROAD, dtt%LDATA_EMIS_ROAD )
IF (.NOT.dtt%LDATA_EMIS_ROAD) DEALLOCATE(dtt%XPAR_EMIS_ROAD)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'HC_ROAD  ','TWN',cfnam_hc_road ,cftyp_hc_road , &
                   xunif_hc_road ,dtt%XPAR_HC_ROAD, dtt%LDATA_HC_ROAD  )  
IF (.NOT.dtt%LDATA_HC_ROAD) DEALLOCATE(dtt%XPAR_HC_ROAD)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'TC_ROAD  ','TWN',cfnam_tc_road ,cftyp_tc_road , &
                   xunif_tc_road ,dtt%XPAR_TC_ROAD, dtt%LDATA_TC_ROAD  )  
IF (.NOT.dtt%LDATA_TC_ROAD) DEALLOCATE(dtt%XPAR_TC_ROAD)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,'ARI','D_ROAD   ','TWN',cfnam_d_road  ,cftyp_d_road  , &
                   xunif_d_road  ,dtt%XPAR_D_ROAD , dtt%LDATA_D_ROAD  )
IF (.NOT.dtt%LDATA_D_ROAD) DEALLOCATE(dtt%XPAR_D_ROAD)
!  
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'ALB_WALL   ','TWN',cfnam_alb_wall  ,cftyp_alb_wall  ,xunif_alb_wall  ,&
        dtt%XPAR_ALB_WALL, dtt%LDATA_ALB_WALL   )
IF (.NOT.dtt%LDATA_ALB_WALL) DEALLOCATE(dtt%XPAR_ALB_WALL)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'EMIS_WALL  ','TWN',cfnam_emis_wall ,cftyp_emis_wall ,xunif_emis_wall ,&
        dtt%XPAR_EMIS_WALL, dtt%LDATA_EMIS_WALL  )
IF (.NOT.dtt%LDATA_EMIS_WALL) DEALLOCATE(dtt%XPAR_EMIS_WALL)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'HC_WALL  ','TWN',cfnam_hc_wall ,cftyp_hc_wall , &
                   xunif_hc_wall ,dtt%XPAR_HC_WALL, dtt%LDATA_HC_WALL  ) 
IF (.NOT.dtt%LDATA_HC_WALL) DEALLOCATE(dtt%XPAR_HC_WALL)
! 
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'TC_WALL  ','TWN',cfnam_tc_wall ,cftyp_tc_wall , &
                   xunif_tc_wall ,dtt%XPAR_TC_WALL, dtt%LDATA_TC_WALL  ) 
IF (.NOT.dtt%LDATA_TC_WALL) DEALLOCATE(dtt%XPAR_TC_WALL)
! 
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,cbld_atype,'D_WALL   ','TWN',cfnam_d_wall  ,cftyp_d_wall  , &
                   xunif_d_wall  ,dtt%XPAR_D_WALL , dtt%LDATA_D_WALL  ) 
IF (.NOT.dtt%LDATA_D_WALL) DEALLOCATE(dtt%XPAR_D_WALL)
! 
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','H_TRAFFIC  ','TWN',cfnam_h_traffic  ,cftyp_h_traffic  ,xunif_h_traffic  ,&
        dtt%XPAR_H_TRAFFIC, dtt%LDATA_H_TRAFFIC   )
IF (.NOT.dtt%LDATA_H_TRAFFIC) DEALLOCATE(dtt%XPAR_H_TRAFFIC)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','LE_TRAFFIC ','TWN',cfnam_le_traffic ,cftyp_le_traffic ,xunif_le_traffic ,&
        dtt%XPAR_LE_TRAFFIC, dtt%LDATA_LE_TRAFFIC  )
IF (.NOT.dtt%LDATA_LE_TRAFFIC) DEALLOCATE(dtt%XPAR_LE_TRAFFIC)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','H_INDUSTRY ','TWN',cfnam_h_industry ,cftyp_h_industry ,xunif_h_industry ,&
        dtt%XPAR_H_INDUSTRY, dtt%LDATA_H_INDUSTRY  )
IF (.NOT.dtt%LDATA_H_INDUSTRY) DEALLOCATE(dtt%XPAR_H_INDUSTRY)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','LE_INDUSTRY','TWN',cfnam_le_industry,cftyp_le_industry,xunif_le_industry,&
        dtt%XPAR_LE_INDUSTRY, dtt%LDATA_LE_INDUSTRY )
IF (.NOT.dtt%LDATA_LE_INDUSTRY) DEALLOCATE(dtt%XPAR_LE_INDUSTRY)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'ROUGH_ROOF','TWN',cfnam_rough_roof,cftyp_rough_roof,xunif_rough_roof ,&
        dtt%XPAR_ROUGH_ROOF,dtt%LDATA_ROUGH_ROOF)
IF (.NOT.dtt%LDATA_ROUGH_ROOF) DEALLOCATE(dtt%XPAR_ROUGH_ROOF)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'ROUGH_WALL','TWN',cfnam_rough_wall,cftyp_rough_wall,xunif_rough_wall ,&
        dtt%XPAR_ROUGH_WALL,dtt%LDATA_ROUGH_WALL)
IF (.NOT.dtt%LDATA_ROUGH_WALL) DEALLOCATE(dtt%XPAR_ROUGH_WALL)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'RESIDENTIAL','TWN',cfnam_residential,cftyp_residential,xunif_residential ,&
        dtt%XPAR_RESIDENTIAL,dtt%LDATA_RESIDENTIAL)
IF (.NOT.dtt%LDATA_RESIDENTIAL) DEALLOCATE(dtt%XPAR_RESIDENTIAL)
!-------------------------------------------------------------------------------
!
!* coherence checks
!
 CALL coherence_thermal_data('ROAD',dtt%LDATA_HC_ROAD,dtt%LDATA_TC_ROAD,dtt%LDATA_D_ROAD)
 CALL coherence_thermal_data('ROOF',dtt%LDATA_HC_ROOF,dtt%LDATA_TC_ROOF,dtt%LDATA_D_ROOF)
 CALL coherence_thermal_data('WALL',dtt%LDATA_HC_WALL,dtt%LDATA_TC_WALL,dtt%LDATA_D_WALL)

!-------------------------------------------------------------------------------
!
!* road directions
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','ROAD_DIR   ','TWN',cfnam_road_dir  ,cftyp_road_dir    ,xunif_road_dir   ,&
        dtt%XPAR_ROAD_DIR, dtt%LDATA_ROAD_DIR    )
IF (.NOT.dtt%LDATA_ROAD_DIR) DEALLOCATE(dtt%XPAR_ROAD_DIR)
!
!-------------------------------------------------------------------------------
!
!* solar panels
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','EMIS_PANEL ','BLD',cfnam_emis_panel,cftyp_emis_panel,xunif_emis_panel,&
       dtt%XPAR_EMIS_PANEL, dtt%LDATA_EMIS_PANEL    )
IF (.NOT.dtt%LDATA_EMIS_PANEL) DEALLOCATE(dtt%XPAR_EMIS_PANEL)
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','ALB_PANEL  ','BLD',cfnam_alb_panel ,cftyp_alb_panel ,xunif_alb_panel ,&
       dtt%XPAR_ALB_PANEL , dtt%LDATA_ALB_PANEL     )
IF (.NOT.dtt%LDATA_ALB_PANEL ) DEALLOCATE(dtt%XPAR_ALB_PANEL )
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','EFF_PANEL  ','BLD',cfnam_eff_panel ,cftyp_eff_panel ,xunif_eff_panel ,&
       dtt%XPAR_EFF_PANEL , dtt%LDATA_EFF_PANEL     )
IF (.NOT.dtt%LDATA_EFF_PANEL ) DEALLOCATE(dtt%XPAR_EFF_PANEL )
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','FRAC_PANEL ','BLD',cfnam_frac_panel,cftyp_frac_panel,xunif_frac_panel,&
       dtt%XPAR_FRAC_PANEL, dtt%LDATA_FRAC_PANEL    )
IF (.NOT.dtt%LDATA_FRAC_PANEL) DEALLOCATE(dtt%XPAR_FRAC_PANEL)
!
!-------------------------------------------------------------------------------
!
!* greenroof fraction
!
IF (ogreenroof) THEN
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,cbld_atype,'GREENROOF','BLD',cfnam_greenroof,cftyp_greenroof,xunif_greenroof ,&
        dtt%XPAR_GREENROOF,dtt%LDATA_GREENROOF)
  IF (.NOT.dtt%LDATA_GREENROOF) DEALLOCATE(dtt%XPAR_GREENROOF)
ELSE IF ( (xunif_greenroof/=0. .AND. xunif_greenroof/=xundef) .OR. len_trim(cfnam_greenroof)/=0) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) ' You chose not to include greenroofs in urban areas : LGREENROOF=.FALSE.     '
  WRITE(iluout,*) ' But            '
  IF (xunif_greenroof/=0. .AND. xunif_greenroof/=xundef) THEN
    WRITE(iluout,*) ' You also chose a greenroof fraction that is not zero : XUNIF_GREENROOF=',xunif_greenroof
  ELSE
    WRITE(iluout,*) ' You also chose a greenroof fraction that is not zero : CFNAM_GREENROOF=',cfnam_greenroof
  END IF
  WRITE(iluout,*) '- - - - - - - - - - - - - - - - - - - - - - -'
  WRITE(iluout,*) ' Please choose either:'
  WRITE(iluout,*) ' LGREENROOF=.TRUE. or set GREENROOF fraction to zero (XUNIF_GREENROOF=0.) in namelist PGD_TEB_PAR'
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx('PGD_TEB_PAR: GREENROOF flag and GREENROOF fraction not coherent')
END IF
!
!-------------------------------------------------------------------------------
!
!* gardens
!
IF (ogarden) THEN
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','GARDEN     ','TWN',cfnam_garden    ,cftyp_garden    ,xunif_garden    ,&
        dtt%XPAR_GARDEN, dtt%LDATA_GARDEN    )
  IF (.NOT.dtt%LDATA_GARDEN) DEALLOCATE(dtt%XPAR_GARDEN)
ELSE IF ( (xunif_garden/=0. .AND. xunif_garden/=xundef) .OR. len_trim(cfnam_garden)/=0) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) ' You chose not to include gardens in urban areas : LGARDEN=.FALSE.     '
  WRITE(iluout,*) ' But            '
  IF (xunif_garden/=0. .AND. xunif_garden/=xundef) THEN
    WRITE(iluout,*) ' You also chose a garden fraction that is not zero : XUNIF_GARDEN=',xunif_garden
  ELSE
    WRITE(iluout,*) ' You also chose a garden fraction that is not zero : CFNAM_GARDEN=',cfnam_garden
  END IF
  WRITE(iluout,*) '- - - - - - - - - - - - - - - - - - - - - - -'
  WRITE(iluout,*) ' Please choose either:'
  WRITE(iluout,*) ' LGARDEN=.TRUE. or set GARDEN fraction to zero (XUNIF_GARDEN=0.) in namelist PGD_TEB_PAR'
  WRITE(iluout,*) '- - - - - - - - - - - - - - - - - - - - - - -'
  WRITE(iluout,*) ' Beware that in this case, it may be necessary to change the'
  WRITE(iluout,*) ' road fraction if you want to keep the same canyon aspect ratio'
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx('PGD_TEB_PAR: GARDEN flag and GARDEN fraction not coherent')
END IF
!
!
!-------------------------------------------------------------------------------
IF (lhook)   CALL dr_hook('PGD_TEB_PAR',1,zhook_handle)
!-------------------------------------------------------------------------------
 CONTAINS
SUBROUTINE coherence_thermal_data(HTYPE,ODATA_HC,ODATA_TC,ODATA_D)
 CHARACTER(LEN=4), INTENT(IN) :: htype
LOGICAL,          INTENT(IN) :: odata_hc
LOGICAL,          INTENT(IN) :: odata_tc
LOGICAL,          INTENT(IN) :: odata_d
!
IF (odata_hc .OR. odata_tc .OR. odata_d) THEN
  IF (.NOT. (odata_hc .AND. odata_tc .AND. odata_d)) THEN
    WRITE(iluout,*) '*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*'
    WRITE(iluout,*) 'When specifying data for thermal ',trim(htype),' characteristics,'
    WRITE(iluout,*) 'All three parameters MUST be defined:'
    WRITE(iluout,*) 'Heat capacity, Thermal conductivity and depths of layers'
    WRITE(iluout,*) ' '
    WRITE(iluout,*) 'In your case :'
    IF (odata_hc) THEN
      WRITE(iluout,*) 'Heat capacity is defined'
    ELSE
      WRITE(iluout,*) 'Heat capacity is NOT defined'
    END IF
    IF (odata_tc) THEN
      WRITE(iluout,*) 'Thermal conductivity is defined'
    ELSE
      WRITE(iluout,*) 'Thermal conductivity is NOT defined'
    END IF
    IF (odata_d) THEN
      WRITE(iluout,*) 'Depths of layers are defined'
    ELSE
      WRITE(iluout,*) 'Depths of layers are NOT defined'
    END IF
    WRITE(iluout,*) '*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*'
    CALL abor1_sfx('Heat capacity, Thermal conductivity and depths of layers MUST all be defined for '//htype)
  END IF
END IF
END SUBROUTINE coherence_thermal_data
!-------------------------------------------------------------------------------
!
END SUBROUTINE pgd_teb_par