pgd_bem_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_bem_par (DTCO, DGU, UG, U, USS, DTB, DTI, TG, &
                              hprogram,oautosize)
!     ##########################################
!
!!**** *PGD_BEM_PAR* monitor for averaging and interpolations of BEM input data
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    G. Pigeon        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    08/2011
!!    G. Pigeon   09/2012, NPAR_FLOOR_LAYER default to 1
!!
!!
!----------------------------------------------------------------------------
!
!*    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_data_bem_n, ONLY : data_bem_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
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_get_luout
USE modi_open_namelist
USE modi_close_namelist
USE modi_test_nam_var_surf
USE modi_ini_var_from_data_0d
USE modi_ini_var_from_data
USE modi_abor1_sfx
!
USE mode_pos_surf
!
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(data_bem_t), INTENT(INOUT) :: dtb
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) :: oautosize    ! T for automatic determination
!                                            ! of HVAC systems charcateristics
!
!
!*    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               :: jlayer    ! loop counter on layers
!
!
!*    0.3    Declaration of namelists
!            ------------------------
!
! bem options
!
INTEGER, PARAMETER                      :: nfloor_max  = 9
INTEGER                                 :: npar_floor_layer ! number of floor layers
!
REAL                                    :: xunif_shade    ! 
 CHARACTER(LEN=28)                       :: cfnam_shade    ! 
 CHARACTER(LEN=6)                        :: cftyp_shade    ! 
REAL                                    :: xunif_natvent    ! 
 CHARACTER(LEN=28)                       :: cfnam_natvent    ! 
 CHARACTER(LEN=6)                        :: cftyp_natvent    ! 
!
! Floor parameters
!
REAL, DIMENSION(NFLOOR_MAX)             :: xunif_hc_floor     ! floor layers heat capacity (J/K/m3)
REAL, DIMENSION(NFLOOR_MAX)             :: xunif_tc_floor     ! floor layers thermal conduc (W/K/m)
REAL, DIMENSION(NFLOOR_MAX)             :: xunif_d_floor      ! depth of floor layers      (m)
REAL                                    :: xunif_floor_height ! building floor height [m]
 CHARACTER(LEN=28), DIMENSION(NFLOOR_MAX):: cfnam_hc_floor     ! file name for HC_FLOOR   
 CHARACTER(LEN=28), DIMENSION(NFLOOR_MAX):: cfnam_tc_floor     ! file name for TC_FLOOR
 CHARACTER(LEN=28), DIMENSION(NFLOOR_MAX):: cfnam_d_floor      ! file name for D_FLOOR
 CHARACTER(LEN=28)                       :: cfnam_floor_height ! file name for FLOOR_HEIGHT
 CHARACTER(LEN=6),  DIMENSION(NFLOOR_MAX):: cftyp_hc_floor     ! file type for HC_FLOOR   
 CHARACTER(LEN=6),  DIMENSION(NFLOOR_MAX):: cftyp_tc_floor     ! file type for TC_FLOOR
 CHARACTER(LEN=6),  DIMENSION(NFLOOR_MAX):: cftyp_d_floor      ! file type for D_FLOOR
 CHARACTER(LEN=6)                        :: cftyp_floor_height ! file type for FLOOR_HEIGHT
!
! AC systems parameters
!
REAL                                    :: xunif_tcool_target  !cooling setpoint
REAL                                    :: xunif_theat_target  !heating setpoint
REAL                                    :: xunif_f_waste_can  ! fraction of waste heat into the canyon
REAL                                    :: xunif_eff_heat     ! efficiency of the heating system
REAL                                    :: xunif_hr_target    ! Relative humidity setpoint
REAL                                    :: xunif_cap_sys_heat ! Capacity of the heating system 
REAL                                    :: xunif_cap_sys_rat  ! Rated capacity of the cooling system
REAL                                    :: xunif_t_adp        ! Apparatus dewpoint temperature of the
REAL                                    :: xunif_m_sys_rat    ! Rated HVAC mass flow rate 
REAL                                    :: xunif_cop_rat      ! Rated COP of the cooling system
REAL                                    :: xunif_f_water_cond ! fraction of evaporation of condensers
 CHARACTER(LEN=28)                       :: cfnam_tcool_target ! file name for TCOOL_TARGET
 CHARACTER(LEN=28)                       :: cfnam_theat_target ! file name for THEAT_TARGET
 CHARACTER(LEN=28)                       :: cfnam_f_waste_can  ! file name for F_WASTE_CAN
 CHARACTER(LEN=28)                       :: cfnam_eff_heat     ! file name for EFF_HEAT
 CHARACTER(LEN=28)                       :: cfnam_hr_target    ! Relative humidity setpoint
 CHARACTER(LEN=28)                       :: cfnam_cap_sys_heat ! Capacity of the heating system 
 CHARACTER(LEN=28)                       :: cfnam_cap_sys_rat  ! Rated capacity of the cooling system
 CHARACTER(LEN=28)                       :: cfnam_t_adp        ! Apparatus dewpoint temperature of the
 CHARACTER(LEN=28)                       :: cfnam_m_sys_rat    ! Rated HVAC mass flow rate 
 CHARACTER(LEN=28)                       :: cfnam_cop_rat      ! Rated COP of the cooling system
 CHARACTER(LEN=28)                       :: cfnam_f_water_cond ! fraction of evaporation of condensers
 CHARACTER(LEN=6)                        :: cftyp_tcool_target ! file type for TCOOL_TARGET
 CHARACTER(LEN=6)                        :: cftyp_theat_target ! file type for THEAT_TARGET
 CHARACTER(LEN=6)                        :: cftyp_f_waste_can  ! file type for F_WASTE_CAN
 CHARACTER(LEN=6)                        :: cftyp_eff_heat     ! file type for EFF_HEAT
 CHARACTER(LEN=6)                        :: cftyp_hr_target    ! Relative humidity setpoint
 CHARACTER(LEN=6)                        :: cftyp_cap_sys_heat ! Capacity of the heating system 
 CHARACTER(LEN=6)                        :: cftyp_cap_sys_rat  ! Rated capacity of the cooling system
 CHARACTER(LEN=6)                        :: cftyp_t_adp        ! Apparatus dewpoint temperature of the
 CHARACTER(LEN=6)                        :: cftyp_m_sys_rat    ! Rated HVAC mass flow rate 
 CHARACTER(LEN=6)                        :: cftyp_cop_rat      ! Rated COP of the cooling system
 CHARACTER(LEN=6)                        :: cftyp_f_water_cond ! fraction of evaporation of condensers
!
! Internal heat gains
REAL                                    :: xunif_qin          ! internal heat gains [W m-2(floor)]
REAL                                    :: xunif_qin_frad     ! radiant fraction of int heat gains
REAL                                    :: xunif_qin_flat     ! Latent franction of internal heat gains
 CHARACTER(LEN=28)                       :: cfnam_qin          ! file name for QIN
 CHARACTER(LEN=28)                       :: cfnam_qin_frad     ! file name for QIN_FRAD
 CHARACTER(LEN=28)                       :: cfnam_qin_flat     ! Latent franction of internal heat gains
 CHARACTER(LEN=6)                        :: cftyp_qin          ! file type for QIN
 CHARACTER(LEN=6)                        :: cftyp_qin_frad     ! file type for QIN_FRAD
 CHARACTER(LEN=6)                        :: cftyp_qin_flat     ! Latent franction of internal heat gains
!
! window parameters
REAL                                    :: xunif_gr           ! glazing ratio
REAL                                    :: xunif_shgc         ! solar transmitance of windows
REAL                                    :: xunif_shgc_sh      ! solar transmitance of windows + shading
REAL                                    :: xunif_u_win        ! glazing thermal resistance[K m W-2]
 CHARACTER(LEN=28)                       :: cfnam_gr           ! file name for GR
 CHARACTER(LEN=28)                       :: cfnam_shgc         ! file name for SHGC
 CHARACTER(LEN=28)                       :: cfnam_shgc_sh      ! file name for SHGC_SH
 CHARACTER(LEN=28)                       :: cfnam_u_win        ! file name for U_WIN
 CHARACTER(LEN=6)                        :: cftyp_gr           ! file type for GR
 CHARACTER(LEN=6)                        :: cftyp_shgc         ! file type for SHGC
 CHARACTER(LEN=6)                        :: cftyp_shgc_sh      ! file type for SHGC
 CHARACTER(LEN=6)                        :: cftyp_u_win        ! file type for U_WIN
!
! air renewal
REAL                                    :: xunif_inf          ! infiltration/ventilation flow rate [AC/H]
REAL                                    :: xunif_v_vent       ! Ventilation flow rate [AC/H]
 CHARACTER(LEN=28)                       :: cfnam_inf          ! file name for INF
 CHARACTER(LEN=28)                       :: cfnam_v_vent       ! Ventilation flow rate [AC/H]
 CHARACTER(LEN=6)                        :: cftyp_inf          ! file type for INF
 CHARACTER(LEN=6)                        :: cftyp_v_vent       ! Ventilation flow rate [AC/H]
!
! parameters for autosize calculation of the AC systems
REAL                                    :: xunif_t_size_max   !
REAL                                    :: xunif_t_size_min    !
 CHARACTER(LEN=28)                       :: cfnam_t_size_max    ! 
 CHARACTER(LEN=28)                       :: cfnam_t_size_min    ! 
 CHARACTER(LEN=6)                        :: cftyp_t_size_max    !
 CHARACTER(LEN=6)                        :: cftyp_t_size_min    !
!
REAL, DIMENSION(TG%NDIM) :: zwork
REAL(KIND=JPRB) :: zhook_handle
!
namelist/nam_data_bem/ npar_floor_layer,                          &
                  xunif_hc_floor, xunif_tc_floor, xunif_d_floor,  &
                  xunif_floor_height,                             &
                  xunif_tcool_target, xunif_theat_target,         &
                  xunif_f_waste_can, xunif_eff_heat,              &
                  xunif_f_water_cond, xunif_hr_target,            &
                  xunif_qin, xunif_qin_frad, xunif_qin_flat,      &   
                  xunif_shgc, xunif_u_win, xunif_gr,xunif_shgc_sh,&
                  xunif_inf, xunif_v_vent,                        &
                  xunif_cap_sys_heat,                             &
                  xunif_cap_sys_rat, xunif_t_adp, xunif_m_sys_rat,&
                  xunif_cop_rat, xunif_t_size_max,                &
                  xunif_t_size_min,                               &
                  xunif_shade, cfnam_shade, cftyp_shade,          &
                  xunif_natvent, cfnam_natvent, cftyp_natvent,    &
                  cfnam_hc_floor, cfnam_tc_floor, cfnam_d_floor,  &
                  cfnam_floor_height,                             &
                  cfnam_tcool_target, cfnam_theat_target,         &
                  cfnam_f_waste_can, cfnam_eff_heat,              &
                  cfnam_f_water_cond, cfnam_hr_target,            &
                  cfnam_qin, cfnam_qin_frad, cfnam_qin_flat,      &   
                  cfnam_shgc, cfnam_u_win, cfnam_gr,              &
                  cfnam_shgc_sh, cfnam_inf, cfnam_v_vent,         &
                  cfnam_cap_sys_heat,                             &
                  cfnam_cap_sys_rat, cfnam_t_adp, cfnam_m_sys_rat,&
                  cfnam_cop_rat, cfnam_t_size_max,                &
                  cfnam_t_size_min,                               &
                  cftyp_hc_floor, cftyp_tc_floor, cftyp_d_floor,  &
                  cftyp_floor_height,                             &
                  cftyp_tcool_target, cftyp_theat_target,         &
                  cftyp_f_waste_can, cftyp_eff_heat,              &
                  cftyp_f_water_cond, cftyp_hr_target,            &
                  cftyp_qin, cftyp_qin_frad, cftyp_qin_flat,      &   
                  cftyp_shgc, cftyp_u_win, cftyp_gr,              &
                  cftyp_shgc_sh, cftyp_inf, cftyp_v_vent,         &
                  cftyp_cap_sys_heat,                             &
                  cftyp_cap_sys_rat, cftyp_t_adp, cftyp_m_sys_rat,&
                  cftyp_cop_rat, cftyp_t_size_max,                &
                  cftyp_t_size_min
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('PGD_BEM_PAR',0,zhook_handle)
!
!*    1.      Initializations
!             ---------------
!
npar_floor_layer = 1
xunif_shade   = xundef
xunif_natvent = xundef
xunif_hc_floor     = xundef
xunif_tc_floor     = xundef
xunif_d_floor      = xundef
xunif_tcool_target = xundef
xunif_theat_target = xundef
xunif_f_waste_can  = xundef
xunif_eff_heat     = xundef
xunif_qin          = xundef
xunif_qin_frad     = xundef
xunif_shgc         = xundef
xunif_u_win        = xundef
xunif_gr           = xundef
xunif_shgc_sh      = xundef
xunif_floor_height = xundef
xunif_inf          = xundef
xunif_f_water_cond = xundef
xunif_qin_flat     = xundef
xunif_hr_target    = xundef
xunif_v_vent       = xundef
xunif_cap_sys_heat = xundef
xunif_cap_sys_rat  = xundef
xunif_t_adp        = xundef
xunif_m_sys_rat    = xundef
xunif_cop_rat      = xundef
xunif_t_size_max   = xundef
xunif_t_size_min   = xundef
!
 cfnam_shade        = '                            '
 cfnam_natvent      = '                            '
 cfnam_hc_floor(:) = '                            '
 cfnam_tc_floor(:) = '                            '
 cfnam_d_floor(:) = '                            '
 cfnam_tcool_target = '                            '
 cfnam_theat_target = '                            '
 cfnam_f_waste_can  = '                            '
 cfnam_eff_heat     = '                            '
 cfnam_qin          = '                            '
 cfnam_qin_frad     = '                            '
 cfnam_shgc         = '                            '
 cfnam_u_win        = '                            '
 cfnam_gr           = '                            '
 cfnam_shgc_sh      = '                            '
 cfnam_floor_height = '                            '
 cfnam_inf          = '                            '
 cfnam_f_water_cond = '                            '
 cfnam_qin_flat     = '                            '
 cfnam_hr_target    = '                            '
 cfnam_v_vent       = '                            '
 cfnam_cap_sys_heat = '                            '
 cfnam_cap_sys_rat  = '                            '
 cfnam_t_adp        = '                            '
 cfnam_m_sys_rat    = '                            '
 cfnam_cop_rat      = '                            '
 cfnam_t_size_max   = '                            '
 cfnam_t_size_min   = '                            '
!
 cftyp_shade        = '      '
 cftyp_natvent      = '      '
 cftyp_hc_floor(:)  = '      '
 cftyp_tc_floor(:)  = '      '
 cftyp_d_floor(:)   = '      '
 cftyp_tcool_target = '      '
 cftyp_theat_target = '      '
 cftyp_f_waste_can  = '      '
 cftyp_eff_heat     = '      '
 cftyp_qin          = '      '
 cftyp_qin_frad     = '      '
 cftyp_shgc         = '      '
 cftyp_u_win        = '      '
 cftyp_gr           = '      '
 cftyp_shgc_sh      = '      '
 cftyp_floor_height = '      '
 cftyp_inf          = '      '
 cftyp_f_water_cond = '      '
 cftyp_qin_flat     = '      '
 cftyp_hr_target    = '      '
 cftyp_v_vent       = '      '
 cftyp_cap_sys_heat = '      '
 cftyp_cap_sys_rat  = '      '
 cftyp_t_adp        = '      '
 cftyp_m_sys_rat    = '      '
 cftyp_cop_rat      = '      '
 cftyp_t_size_max   = '      '
 cftyp_t_size_min   = '      '
!
!-------------------------------------------------------------------------------
!
!*    2.      Input file for cover types
!             --------------------------
!
 CALL get_luout(hprogram,iluout)
 CALL open_namelist(hprogram,ilunam)
!
 CALL posnam(ilunam,'NAM_DATA_BEM',gfound,iluout)
IF (gfound) READ(unit=ilunam,nml=nam_data_bem)
!
 CALL close_namelist(hprogram,ilunam)
!
!-------------------------------------------------------------------------------
!
dtb%NPAR_FLOOR_LAYER = npar_floor_layer
!
!-------------------------------------------------------------------------------
!
!* coherence check
!
IF ((     any(xunif_hc_floor/=xundef) .OR. any(len_trim(cfnam_hc_floor)>0) &
     .OR. any(xunif_tc_floor/=xundef) .OR. any(len_trim(cfnam_tc_floor)>0) &
     .OR. any(xunif_d_floor /=xundef) .OR. any(len_trim(cfnam_d_floor )>0) &
    ) .AND. npar_floor_layer<1                                  ) THEN
  CALL abor1_sfx('In order to initialize FLOOR thermal quantities, please specify NPAR_FLOOR_LAYER in namelist NAM_DATA_BEM')
END IF
!
!-------------------------------------------------------------------------------
!
ALLOCATE(dtb%XPAR_HC_FLOOR    (tg%NDIM,npar_floor_layer))
ALLOCATE(dtb%XPAR_TC_FLOOR    (tg%NDIM,npar_floor_layer))
ALLOCATE(dtb%XPAR_D_FLOOR     (tg%NDIM,npar_floor_layer))
!
ALLOCATE(dtb%XPAR_TCOOL_TARGET  (tg%NDIM))
ALLOCATE(dtb%XPAR_THEAT_TARGET  (tg%NDIM))
ALLOCATE(dtb%XPAR_F_WASTE_CAN   (tg%NDIM))
ALLOCATE(dtb%XPAR_EFF_HEAT      (tg%NDIM))
ALLOCATE(dtb%XPAR_QIN           (tg%NDIM))
ALLOCATE(dtb%XPAR_QIN_FRAD      (tg%NDIM))
ALLOCATE(dtb%XPAR_SHGC          (tg%NDIM))
ALLOCATE(dtb%XPAR_U_WIN         (tg%NDIM))
ALLOCATE(dtb%XPAR_GR            (tg%NDIM))
ALLOCATE(dtb%XPAR_SHGC_SH       (tg%NDIM))
ALLOCATE(dtb%XPAR_FLOOR_HEIGHT  (tg%NDIM))
ALLOCATE(dtb%XPAR_INF           (tg%NDIM))
!
ALLOCATE(dtb%XPAR_F_WATER_COND (tg%NDIM))
ALLOCATE(dtb%XPAR_QIN_FLAT     (tg%NDIM))
ALLOCATE(dtb%XPAR_HR_TARGET    (tg%NDIM))
ALLOCATE(dtb%XPAR_V_VENT       (tg%NDIM))
ALLOCATE(dtb%XPAR_CAP_SYS_HEAT (tg%NDIM))
ALLOCATE(dtb%XPAR_CAP_SYS_RAT  (tg%NDIM))
ALLOCATE(dtb%XPAR_T_ADP        (tg%NDIM))
ALLOCATE(dtb%XPAR_M_SYS_RAT    (tg%NDIM))
ALLOCATE(dtb%XPAR_COP_RAT      (tg%NDIM))
ALLOCATE(dtb%XPAR_T_SIZE_MAX   (tg%NDIM))
ALLOCATE(dtb%XPAR_T_SIZE_MIN   (tg%NDIM))
!
ALLOCATE(dtb%XPAR_SHADE        (tg%NDIM))
ALLOCATE(dtb%XPAR_NATVENT      (tg%NDIM))
!
!-------------------------------------------------------------------------------
IF (nfloor_max < npar_floor_layer) THEN
  WRITE(iluout,*) '---------------------------------------------'
  WRITE(iluout,*) 'Please update pgd_bem_par.F90 routine :      '
  WRITE(iluout,*) 'The maximum number of FLOOR LAYER             '
  WRITE(iluout,*) 'in the declaration of the namelist variables '
  WRITE(iluout,*) 'must be increased to : ', npar_floor_layer
  WRITE(iluout,*) '---------------------------------------------'
  CALL abor1_sfx('PGD_BEM_PAR: MAXIMUM NUMBER OF NPAR_FLOOR_LAYER MUST BE INCREASED')
END IF
!-------------------------------------------------------------------------------
!
!*    3.      user defined fields are prescribed
!             ----------------------------------
!
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'MAJ','SHADE      ','TWN', cfnam_shade, cftyp_shade, xunif_shade, &
        dtb%XPAR_SHADE, dtb%LDATA_SHADE )
IF (.NOT.dtb%LDATA_SHADE) DEALLOCATE(dtb%XPAR_SHADE)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'MAJ','NATVENT   ','TWN', cfnam_natvent, cftyp_natvent, xunif_natvent, &
        dtb%XPAR_NATVENT, dtb%LDATA_NATVENT )
IF (.NOT.dtb%LDATA_NATVENT) DEALLOCATE(dtb%XPAR_NATVENT)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,'INV','HC_FLOOR  ','TWN',cfnam_hc_floor,cftyp_hc_floor, &
        xunif_hc_floor,dtb%XPAR_HC_FLOOR,dtb%LDATA_HC_FLOOR )
IF (.NOT.dtb%LDATA_HC_FLOOR) DEALLOCATE(dtb%XPAR_HC_FLOOR)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,'ARI','TC_FLOOR  ','TWN',cfnam_tc_floor,cftyp_tc_floor, &
                 xunif_tc_floor ,dtb%XPAR_TC_FLOOR, dtb%LDATA_TC_FLOOR )
IF (.NOT.dtb%LDATA_TC_FLOOR) DEALLOCATE(dtb%XPAR_TC_FLOOR)
!
 CALL ini_var_from_data(dtco, dgu, ug, u, uss, dti, &
                        hprogram,'ARI','D_FLOOR   ','TWN',cfnam_d_floor,cftyp_d_floor, &
                 xunif_d_floor  ,dtb%XPAR_D_FLOOR , dtb%LDATA_D_FLOOR )
IF (.NOT.dtb%LDATA_D_FLOOR) DEALLOCATE(dtb%XPAR_D_FLOOR)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','TCOOL_TARGET','TWN',cfnam_tcool_target, cftyp_tcool_target, xunif_tcool_target, &
        dtb%XPAR_TCOOL_TARGET, dtb%LDATA_TCOOL_TARGET)
IF (.NOT.dtb%LDATA_TCOOL_TARGET) DEALLOCATE(dtb%XPAR_TCOOL_TARGET)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','THEAT_TARGET','TWN',cfnam_theat_target, cftyp_theat_target, xunif_theat_target, &
        dtb%XPAR_THEAT_TARGET, dtb%LDATA_THEAT_TARGET)
IF (.NOT.dtb%LDATA_THEAT_TARGET) DEALLOCATE(dtb%XPAR_THEAT_TARGET)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','F_WASTE_CAN','TWN',cfnam_f_waste_can, cftyp_f_waste_can, xunif_f_waste_can, &
        dtb%XPAR_F_WASTE_CAN, dtb%LDATA_F_WASTE_CAN)
IF (.NOT.dtb%LDATA_F_WASTE_CAN) DEALLOCATE(dtb%XPAR_F_WASTE_CAN)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','EFF_HEAT','TWN',cfnam_eff_heat, cftyp_eff_heat, xunif_eff_heat, &
        dtb%XPAR_EFF_HEAT, dtb%LDATA_EFF_HEAT)
IF (.NOT.dtb%LDATA_EFF_HEAT) DEALLOCATE(dtb%XPAR_EFF_HEAT)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','QIN','TWN',cfnam_qin, cftyp_qin, xunif_qin, dtb%XPAR_QIN, dtb%LDATA_QIN)
IF (.NOT.dtb%LDATA_QIN) DEALLOCATE(dtb%XPAR_QIN)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','QIN_FRAD','TWN',cfnam_qin_frad, cftyp_qin_frad, xunif_qin_frad, &
        dtb%XPAR_QIN_FRAD, dtb%LDATA_QIN_FRAD)
IF (.NOT.dtb%LDATA_QIN_FRAD) DEALLOCATE(dtb%XPAR_QIN_FRAD)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','SHGC','TWN',cfnam_shgc, cftyp_shgc, xunif_shgc, dtb%XPAR_SHGC, dtb%LDATA_SHGC)
IF (.NOT.dtb%LDATA_SHGC) DEALLOCATE(dtb%XPAR_SHGC)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','U_WIN','TWN',cfnam_u_win, cftyp_u_win, xunif_u_win, dtb%XPAR_U_WIN, dtb%LDATA_U_WIN)
IF (.NOT.dtb%LDATA_U_WIN) DEALLOCATE(dtb%XPAR_U_WIN)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','GR','TWN',cfnam_gr, cftyp_gr, xunif_gr, dtb%XPAR_GR, dtb%LDATA_GR)
IF (.NOT.dtb%LDATA_GR) DEALLOCATE(dtb%XPAR_GR)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','SHGC_SH','TWN',cfnam_shgc_sh, cftyp_shgc_sh, xunif_shgc_sh, &
        dtb%XPAR_SHGC_SH, dtb%LDATA_SHGC_SH)
IF (.NOT.dtb%LDATA_SHGC_SH) DEALLOCATE(dtb%XPAR_SHGC_SH)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','FLOOR_HEIGHT','TWN',cfnam_floor_height, cftyp_floor_height, xunif_floor_height, &
        dtb%XPAR_FLOOR_HEIGHT, dtb%LDATA_FLOOR_HEIGHT)
IF (.NOT.dtb%LDATA_FLOOR_HEIGHT) DEALLOCATE(dtb%XPAR_FLOOR_HEIGHT)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','INF','TWN',cfnam_inf, cftyp_inf, xunif_inf, dtb%XPAR_INF, dtb%LDATA_INF)
IF (.NOT.dtb%LDATA_INF) DEALLOCATE(dtb%XPAR_INF)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','F_WATER_COND','TWN',cfnam_f_water_cond, cftyp_f_water_cond, xunif_f_water_cond, &
        dtb%XPAR_F_WATER_COND, dtb%LDATA_F_WATER_COND)
IF (.NOT.dtb%LDATA_F_WATER_COND) DEALLOCATE(dtb%XPAR_F_WATER_COND)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','QIN_FLAT','TWN',cfnam_qin_flat, cftyp_qin_flat, xunif_qin_flat, &
        dtb%XPAR_QIN_FLAT, dtb%LDATA_QIN_FLAT)
IF (.NOT.dtb%LDATA_QIN_FLAT) DEALLOCATE(dtb%XPAR_QIN_FLAT)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','HR_TARGET','TWN',cfnam_hr_target, cftyp_hr_target, xunif_hr_target, &
        dtb%XPAR_HR_TARGET, dtb%LDATA_HR_TARGET)
IF (.NOT.dtb%LDATA_HR_TARGET) DEALLOCATE(dtb%XPAR_HR_TARGET)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','V_VENT','TWN',cfnam_v_vent, cftyp_v_vent, &
                           xunif_v_vent, dtb%XPAR_V_VENT, dtb%LDATA_V_VENT)
IF (.NOT.dtb%LDATA_V_VENT) DEALLOCATE(dtb%XPAR_V_VENT)
!
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','T_SIZE_MAX','TWN',cfnam_t_size_max, cftyp_t_size_max, xunif_t_size_max, &
        dtb%XPAR_T_SIZE_MAX, dtb%LDATA_T_SIZE_MAX)
IF (.NOT.dtb%LDATA_T_SIZE_MAX) DEALLOCATE(dtb%XPAR_T_SIZE_MAX)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','T_SIZE_MIN','TWN',cfnam_t_size_min, cftyp_t_size_min, xunif_t_size_min, &
        dtb%XPAR_T_SIZE_MIN, dtb%LDATA_T_SIZE_MIN)
IF (.NOT.dtb%LDATA_T_SIZE_MIN) DEALLOCATE(dtb%XPAR_T_SIZE_MIN)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','CAP_SYS_HEAT','TWN',cfnam_cap_sys_heat, cftyp_cap_sys_heat, xunif_cap_sys_heat, &
        dtb%XPAR_CAP_SYS_HEAT, dtb%LDATA_CAP_SYS_HEAT)
IF (.NOT.dtb%LDATA_CAP_SYS_HEAT) DEALLOCATE(dtb%XPAR_CAP_SYS_HEAT)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','CAP_SYS_RAT','TWN',cfnam_cap_sys_rat, cftyp_cap_sys_rat, xunif_cap_sys_rat, &
        dtb%XPAR_CAP_SYS_RAT, dtb%LDATA_CAP_SYS_RAT)
IF (.NOT.dtb%LDATA_CAP_SYS_RAT) DEALLOCATE(dtb%XPAR_CAP_SYS_RAT)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','M_SYS_RAT','TWN',cfnam_m_sys_rat, cftyp_m_sys_rat, xunif_m_sys_rat, &
        dtb%XPAR_M_SYS_RAT, dtb%LDATA_M_SYS_RAT)
IF (.NOT.dtb%LDATA_M_SYS_RAT) DEALLOCATE(dtb%XPAR_M_SYS_RAT)
  !
IF (oautosize) THEN
  IF (dtb%LDATA_CAP_SYS_HEAT .OR. dtb%LDATA_CAP_SYS_RAT .OR. dtb%LDATA_M_SYS_RAT) THEN
    WRITE(iluout,*) '==> You choose LAUTOSIZE=T <=='
    WRITE(iluout,*) 'Therefore HVAC systems characteristics will be computed automatically'
    IF (dtb%LDATA_CAP_SYS_HEAT) THEN
      WRITE(iluout,*) 'Data you provided for CAP_SYS_HEAT are then discarded.'
      DEALLOCATE(dtb%XPAR_CAP_SYS_HEAT)
    END IF
    IF (dtb%LDATA_CAP_SYS_RAT ) THEN
      WRITE(iluout,*) 'Data you provided for CAP_SYS_RAT  are then discarded.'
      DEALLOCATE(dtb%XPAR_CAP_SYS_RAT)
    END IF
    IF (dtb%LDATA_M_SYS_RAT   ) THEN
      WRITE(iluout,*) 'Data you provided for M_SYS_RAT    are then discarded.'
      DEALLOCATE(dtb%XPAR_M_SYS_RAT)
    END IF
  END IF
  dtb%LDATA_CAP_SYS_HEAT = .false.
  dtb%LDATA_CAP_SYS_RAT  = .false.
  dtb%LDATA_M_SYS_RAT    = .false.
ELSE
  IF (dtb%LDATA_T_SIZE_MAX .OR. dtb%LDATA_T_SIZE_MAX) THEN
    WRITE(iluout,*) '==> You choose LAUTOSIZE=F <=='
    WRITE(iluout,*) 'Therefore HVAC systems characteristics are specified'
    WRITE(iluout,*) 'and you do not need the minimal and maximum temperatures'
    WRITE(iluout,*) 'that would be used if you have chosen an automatic calibration.'
    IF (dtb%LDATA_T_SIZE_MAX) THEN
      WRITE(iluout,*) 'Data you provided for T_SIZE_MAX are then discarded.'
      DEALLOCATE(dtb%XPAR_T_SIZE_MAX)
    END IF
    IF (dtb%LDATA_T_SIZE_MIN) THEN
      WRITE(iluout,*) 'Data you provided for T_SIZE_MIN are then discarded.'
      DEALLOCATE(dtb%XPAR_T_SIZE_MIN)
    END IF
  END IF
  dtb%LDATA_T_SIZE_MAX = .false.
  dtb%LDATA_T_SIZE_MIN = .false.
END IF
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','T_ADP','TWN',cfnam_t_adp, cftyp_t_adp, xunif_t_adp, dtb%XPAR_T_ADP, dtb%LDATA_T_ADP)
IF (.NOT.dtb%LDATA_T_ADP) DEALLOCATE(dtb%XPAR_T_ADP)
!
 CALL ini_var_from_data_0d(dtco, dgu, ug, u, uss, &
                           hprogram,'ARI','COP_RAT','TWN',cfnam_cop_rat, cftyp_cop_rat, xunif_cop_rat, &
        dtb%XPAR_COP_RAT, dtb%LDATA_COP_RAT)
IF (.NOT.dtb%LDATA_COP_RAT) DEALLOCATE(dtb%XPAR_COP_RAT)

!-------------------------------------------------------------------------------
!
!* coherence checks
!
 CALL coherence_thermal_data_fl('FLOOR',dtb%LDATA_HC_FLOOR,dtb%LDATA_TC_FLOOR,dtb%LDATA_D_FLOOR)
!
!-------------------------------------------------------------------------------
IF (lhook) CALL dr_hook('PGD_BEM_PAR',1,zhook_handle)
!-------------------------------------------------------------------------------
 CONTAINS
SUBROUTINE coherence_thermal_data_fl(HTYPE,ODATA_HC,ODATA_TC,ODATA_D)
 CHARACTER(LEN=5), 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_fl
!-------------------------------------------------------------------------------
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE pgd_bem_par