convert_patch_teb.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 convert_patch_teb (BDD, DTB, DTCO, DTT, TOP, &
                                   PDEF_ROAD_DIR, T, B, TPN )
!     ##############################################################
!
!!**** *CONVERT_PATCH_TEB* compilation of conver_cover_isba, pgd_isba_par and
!!      init_from_data_isba_parn
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    S. Faroux        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    16/11/10
!!    G. Pigeon     /09/12: add ROUGH_ROOF and ROUGH_WALL coef for outdoor conv. coef
!!    V. Masson      08/13  adds solar panels
!!
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
!
USE modd_bld_description_n, ONLY : bld_desc_t
USE modd_data_bem_n, ONLY : data_bem_t
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_data_teb_n, ONLY : data_teb_t
USE modd_teb_option_n, ONLY : teb_options_t
!
USE modd_teb_n, ONLY : teb_t
USE modd_bem_n, ONLY : bem_t
USE modd_teb_panel_n, ONLY : teb_panel_t
!
USE modd_surf_par,       ONLY : xundef
!
USE modd_data_cover_par, ONLY : ndata_road_layer, ndata_wall_layer,               &
                                ndata_roof_layer, ndata_floor_layer
USE modd_data_cover,     ONLY : xdata_z0_town, xdata_alb_roof,                    &
                                xdata_emis_roof, xdata_hc_roof, xdata_tc_roof,    &
                                xdata_d_roof, xdata_alb_road, xdata_emis_road,    &
                                xdata_hc_road, xdata_tc_road, xdata_d_road,       &
                                xdata_alb_wall, xdata_emis_wall, xdata_hc_wall,   &
                                xdata_tc_wall, xdata_d_wall, xdata_bld_height,    &
                                xdata_h_traffic, xdata_le_traffic,                &
                                xdata_h_industry, xdata_le_industry,              &
                                xdata_hc_floor, xdata_tc_floor, xdata_d_floor,    &
                                xdata_tcool_target, xdata_theat_target,           &
                                xdata_f_waste_can, xdata_eff_heat, xdata_qin,     &
                                xdata_qin_frad, xdata_shgc, xdata_u_win, xdata_gr,&
                                xdata_shgc_sh, xdata_floor_height, xdata_inf,     &
                                xdata_f_water_cond, xdata_qin_flat,               &   
                                xdata_hr_target, xdata_v_vent, xdata_cap_sys_heat,&
                                xdata_cap_sys_rat, xdata_t_adp, xdata_m_sys_rat,  &
                                xdata_cop_rat, xdata_t_size_max, xdata_t_size_min,&
                                xdata_shade, xdata_natvent, xdata_rough_roof,     &
                                xdata_rough_wall, xdata_residential,xdata_frac_gr,&
                                xdata_emis_panel,xdata_alb_panel, xdata_eff_panel,&
                                xdata_frac_panel
!
!
USE modi_thermal_layers_conf
USE modi_ini_data_param_teb
USE modi_av_pgd
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
USE modi_abor1_sfx
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
!
TYPE(bld_desc_t), INTENT(INOUT) :: BDD
TYPE(data_bem_t), INTENT(INOUT) :: DTB
TYPE(data_cover_t), INTENT(INOUT) :: DTCO
TYPE(data_teb_t), INTENT(INOUT) :: DTT
TYPE(teb_options_t), INTENT(INOUT) :: TOP
!
REAL,                   INTENT(IN)    :: PDEF_ROAD_DIR ! default road direction
!
TYPE(teb_t), INTENT(INOUT), OPTIONAL :: T
TYPE(bem_t), INTENT(INOUT), OPTIONAL :: B
TYPE(teb_panel_t), INTENT(INOUT), OPTIONAL :: TPN
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL, DIMENSION(:,:), ALLOCATABLE :: ZTC, ZHC, ZD
INTEGER               :: JLAYER    ! loop counter on layers
INTEGER               :: IL        ! number of points
!
 CHARACTER(LEN=3)      :: YAVG      ! Type of averaging for buildings
 CHARACTER(LEN=4)      :: YAREA     ! Area where parameter is averaged
!
!*    0.3    Declaration of namelists
!            ------------------------
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!
!-------------------------------------------------------------------------------
!
!*    1.      Initializations
!             ---------------
!
IF (lhook) CALL dr_hook('CONVERT_PATCH_TEB',0,zhook_handle)
!
IF (top%CBLD_ATYPE=='ARI') THEN
  yavg=top%CBLD_ATYPE
  yarea='BLD'
ELSE
  yavg=top%CBLD_ATYPE
  yarea='BLV'
END IF
!
il = SIZE(top%XCOVER,1)
!-------------------------------------------------------------------------------
!
IF (PRESENT(t)) THEN
!       ROAD DIRECTION 
!       --------------
  IF (dtt%LDATA_ROAD_DIR) THEN
    t%XROAD_DIR=dtt%XPAR_ROAD_DIR
  ELSE
    t%XROAD_DIR=pdef_road_dir
  ENDIF
!
!       GARDEN fraction
!       ---------------
  IF (dtt%LDATA_GARDEN) THEN
    t%XGARDEN=dtt%XPAR_GARDEN
  ELSE
    CALL av_pgd(dtco, t%XGARDEN, top%XCOVER, dtco%XDATA_GARDEN(:),'TWN','ARI',top%LCOVER)
  ENDIF
!
!     GREENROOF FRACTION
!     ------------------
  IF (dtt%LDATA_GREENROOF) THEN
    t%XGREENROOF=dtt%XPAR_GREENROOF
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pgreenroof=t%XGREENROOF)
  ELSE
    CALL av_pgd(dtco, t%XGREENROOF,top%XCOVER ,xdata_frac_gr(:),'BLD','ARI',top%LCOVER)
  ENDIF
!
!       BLD fraction
!       -------------
  IF (dtt%LDATA_BLD) THEN
    t%XBLD=dtt%XPAR_BLD  
  ELSE
    CALL av_pgd(dtco, t%XBLD, top%XCOVER, dtco%XDATA_BLD(:),'TWN','ARI',top%LCOVER)
  ENDIF
!
!
!      Z0_TOWN
!      ----------
  IF (dtt%LDATA_Z0_TOWN) THEN
    t%XZ0_TOWN=dtt%XPAR_Z0_TOWN  
  ELSE
    CALL av_pgd(dtco, t%XZ0_TOWN ,top%XCOVER ,xdata_z0_town(:),'TWN','CDN',top%LCOVER)  
  ENDIF  
!
!      BLD Height
!      ----------
  IF (dtt%LDATA_BLD_HEIGHT) THEN
    t%XBLD_HEIGHT=dtt%XPAR_BLD_HEIGHT
  ELSE
    CALL av_pgd(dtco, t%XBLD_HEIGHT ,top%XCOVER ,xdata_bld_height(:),'BLD','ARI',top%LCOVER)  
  ENDIF  
!
!      WALL O HOR
!      ----------
  IF (dtt%LDATA_WALL_O_HOR) THEN
    t%XWALL_O_HOR=dtt%XPAR_WALL_O_HOR 
  ELSE
    CALL av_pgd(dtco, t%XWALL_O_HOR ,top%XCOVER ,dtco%XDATA_WALL_O_HOR (:),'BLD','ARI',top%LCOVER)  
  ENDIF  
!
!      ALB_ROOF
!      ----------
  IF (dtt%LDATA_ALB_ROOF) THEN
    t%XALB_ROOF=dtt%XPAR_ALB_ROOF
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,palb_roof=t%XALB_ROOF)
  ELSE
    CALL av_pgd(dtco, t%XALB_ROOF ,top%XCOVER ,xdata_alb_roof(:),yarea,yavg,top%LCOVER)  
  ENDIF  
!
!      EMIS_ROOF
!      ----------
  IF (dtt%LDATA_EMIS_ROOF) THEN
    t%XEMIS_ROOF=dtt%XPAR_EMIS_ROOF
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pemis_roof=t%XEMIS_ROOF)
  ELSE
    CALL av_pgd(dtco, t%XEMIS_ROOF ,top%XCOVER ,xdata_emis_roof(:),yarea,yavg,top%LCOVER)  
  ENDIF  
!
!      HC_ROOF, TC_ROOF, D_ROOF
!      ------------------------
!
  IF (dtt%LDATA_HC_ROOF) THEN
    CALL allocate_thermal_work(dtt%NPAR_ROOF_LAYER)
    zhc=dtt%XPAR_HC_ROOF
    ztc=dtt%XPAR_TC_ROOF
    zd =dtt%XPAR_D_ROOF
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL allocate_thermal_work(bdd%NDESC_ROOF_LAYER)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,phc_roof=zhc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,ptc_roof=ztc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pd_roof =zd )
  ELSE
    CALL allocate_thermal_work(ndata_roof_layer)
    DO jlayer=1,ndata_roof_layer
      CALL av_pgd(dtco, zhc(:,jlayer), top%XCOVER, xdata_hc_roof(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, ztc(:,jlayer), top%XCOVER ,xdata_tc_roof(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, zd(:,jlayer), top%XCOVER ,xdata_d_roof(:,jlayer),yarea,yavg,top%LCOVER)
    ENDDO
  ENDIF 
  CALL thermal_layers_conf('ROOF ',zd,t%XD_ROOF,zhc,t%XHC_ROOF,ztc,t%XTC_ROOF) 
  CALL deallocate_thermal_work
!
!      ALB_ROAD
!      ----------
  IF (dtt%LDATA_ALB_ROAD) THEN
    t%XALB_ROAD=dtt%XPAR_ALB_ROAD
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd,  dtt%NPAR_BLDCODE,palb_road=t%XALB_ROAD)
  ELSE
    CALL av_pgd(dtco, &
                 t%XALB_ROAD ,top%XCOVER ,xdata_alb_road(:),'STR','ARI',top%LCOVER)  
  ENDIF  
!
!      EMIS_ROAD
!      ----------
  IF (dtt%LDATA_EMIS_ROAD) THEN
    t%XEMIS_ROAD=dtt%XPAR_EMIS_ROAD
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pemis_road=t%XEMIS_ROAD)
  ELSE
    CALL av_pgd(dtco, &
                 t%XEMIS_ROAD ,top%XCOVER ,xdata_emis_road(:),'STR','ARI',top%LCOVER)  
  ENDIF  
!
!      HC_ROAD, TC_ROAD, D_ROAD
!      ------------------------
!
  IF (dtt%LDATA_HC_ROAD) THEN
    CALL allocate_thermal_work(dtt%NPAR_ROAD_LAYER)
    zhc=dtt%XPAR_HC_ROAD
    ztc=dtt%XPAR_TC_ROAD
    zd =dtt%XPAR_D_ROAD
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL allocate_thermal_work(bdd%NDESC_ROAD_LAYER)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,phc_road=zhc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,ptc_road=ztc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pd_road =zd )
  ELSE
    CALL allocate_thermal_work(ndata_road_layer)
    DO jlayer=1,ndata_road_layer
      CALL av_pgd(dtco, zhc(:,jlayer), top%XCOVER, xdata_hc_road(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, ztc(:,jlayer), top%XCOVER ,xdata_tc_road(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, zd(:,jlayer), top%XCOVER ,xdata_d_road(:,jlayer),yarea,yavg,top%LCOVER)
      !CALL AV_PGD (T%XHC_ROAD(:,JLAYER), TOP%XCOVER, XDATA_HC_ROAD (:,JLAYER),YAREA,YAVG)
      !CALL AV_PGD (T%XTC_ROAD(:,JLAYER), TOP%XCOVER ,XDATA_TC_ROAD (:,JLAYER),YAREA,YAVG)
      !CALL AV_PGD (T%XD_ROAD(:,JLAYER), TOP%XCOVER ,XDATA_D_ROAD  (:,JLAYER),YAREA,YAVG)
    ENDDO
  ENDIF 
  CALL thermal_layers_conf('ROAD ',zd,t%XD_ROAD,zhc,t%XHC_ROAD,ztc,t%XTC_ROAD) 
  CALL deallocate_thermal_work
!
!      ALB_WALL
!      ----------
  IF (dtt%LDATA_ALB_WALL) THEN
    t%XALB_WALL=dtt%XPAR_ALB_WALL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,palb_wall=t%XALB_WALL)
  ELSE
    CALL av_pgd(dtco, t%XALB_WALL ,top%XCOVER ,xdata_alb_wall(:),yarea,yavg,top%LCOVER)  
  ENDIF  
!
!      EMIS_WALL
!      ----------
  IF (dtt%LDATA_EMIS_WALL) THEN
    t%XEMIS_WALL=dtt%XPAR_EMIS_WALL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pemis_wall=t%XEMIS_WALL)
  ELSE
    CALL av_pgd(dtco, t%XEMIS_WALL ,top%XCOVER ,xdata_emis_wall(:),yarea,yavg,top%LCOVER)  
  ENDIF  
!
!      HC_WALL, TC_WALL, D_WALL
!      ------------------------
!
  IF (dtt%LDATA_HC_WALL) THEN
    CALL allocate_thermal_work(dtt%NPAR_WALL_LAYER)
    zhc=dtt%XPAR_HC_WALL
    ztc=dtt%XPAR_TC_WALL
    zd =dtt%XPAR_D_WALL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL allocate_thermal_work(bdd%NDESC_WALL_LAYER)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,phc_wall=zhc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,ptc_wall=ztc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pd_wall =zd )
  ELSE
    CALL allocate_thermal_work(ndata_wall_layer)
    DO jlayer=1,ndata_wall_layer
      CALL av_pgd(dtco, zhc(:,jlayer), top%XCOVER, xdata_hc_wall(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, ztc(:,jlayer), top%XCOVER ,xdata_tc_wall(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, zd(:,jlayer), top%XCOVER ,xdata_d_wall(:,jlayer),yarea,yavg,top%LCOVER)
    ENDDO
  ENDIF 
  CALL thermal_layers_conf('WALL ',zd,t%XD_WALL,zhc,t%XHC_WALL,ztc,t%XTC_WALL) 
  CALL deallocate_thermal_work
!
!     H_TRAFFIC
!     ---------
  IF (dtt%LDATA_H_TRAFFIC) THEN
    t%XH_TRAFFIC=dtt%XPAR_H_TRAFFIC
  ELSE
    CALL av_pgd(dtco, t%XH_TRAFFIC ,top%XCOVER ,xdata_h_traffic(:),'TWN','ARI',top%LCOVER)
  ENDIF 
!
!     LE_TRAFFIC
!     ----------
  IF (dtt%LDATA_LE_TRAFFIC) THEN
    t%XLE_TRAFFIC=dtt%XPAR_LE_TRAFFIC
  ELSE
    CALL av_pgd(dtco, t%XLE_TRAFFIC ,top%XCOVER ,xdata_le_traffic(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     H_INDUSTRY
!     ----------
  IF (dtt%LDATA_H_INDUSTRY) THEN
    t%XH_INDUSTRY=dtt%XPAR_H_INDUSTRY
  ELSE
    CALL av_pgd(dtco, t%XH_INDUSTRY ,top%XCOVER ,xdata_h_industry(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     LE_INDUSTRY
!     -----------
  IF (dtt%LDATA_LE_INDUSTRY) THEN
    t%XLE_INDUSTRY=dtt%XPAR_LE_INDUSTRY
  ELSE
    CALL av_pgd(dtco, t%XLE_INDUSTRY ,top%XCOVER ,xdata_le_industry(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!
!      ROUGH_ROOF
!      ----------
  IF (dtt%LDATA_ROUGH_ROOF) THEN
    t%XROUGH_ROOF=dtt%XPAR_ROUGH_ROOF
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,prough_roof=t%XROUGH_ROOF)
  ELSE
    CALL av_pgd(dtco, t%XROUGH_ROOF ,top%XCOVER ,xdata_rough_roof(:),yarea,yavg,top%LCOVER)  
  ENDIF  
!
!      ROUGH_WALL
!      ----------
  IF (dtt%LDATA_ROUGH_WALL) THEN
    t%XROUGH_WALL=dtt%XPAR_ROUGH_WALL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,prough_wall=t%XROUGH_WALL)
  ELSE
    CALL av_pgd(dtco, t%XROUGH_WALL ,top%XCOVER ,xdata_rough_wall(:),yarea,yavg,top%LCOVER)  
  ENDIF  
!
!      RESIDENTIAL use fraction
!      ------------------------
  IF (dtt%LDATA_RESIDENTIAL) THEN
    t%XRESIDENTIAL=dtt%XPAR_RESIDENTIAL
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,presidential=t%XRESIDENTIAL)
  ELSE
    CALL av_pgd(dtco, t%XRESIDENTIAL ,top%XCOVER ,xdata_residential(:),yarea,yavg,top%LCOVER)  
  ENDIF
!
ENDIF
!      HC_FLOOR, TC_FLOOR, D_FLOOR
!      ------------------------
!
IF (PRESENT(b)) THEN

  IF (dtb%LDATA_HC_FLOOR) THEN
    CALL allocate_thermal_work(dtb%NPAR_FLOOR_LAYER)
    zhc=dtb%XPAR_HC_FLOOR
    ztc=dtb%XPAR_TC_FLOOR
    zd =dtb%XPAR_D_FLOOR
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL allocate_thermal_work(bdd%NDESC_FLOOR_LAYER)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,phc_floor=zhc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,ptc_floor=ztc)
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pd_floor =zd )
  ELSE
    CALL allocate_thermal_work(ndata_floor_layer)
    DO jlayer=1,ndata_floor_layer
      CALL av_pgd(dtco, zhc(:,jlayer), top%XCOVER, xdata_hc_floor(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, ztc(:,jlayer), top%XCOVER ,xdata_tc_floor(:,jlayer),yarea,yavg,top%LCOVER)
      CALL av_pgd(dtco, zd(:,jlayer), top%XCOVER ,xdata_d_floor(:,jlayer),yarea,yavg,top%LCOVER)
    ENDDO
  ENDIF 
  CALL thermal_layers_conf('FLOOR',zd,b%XD_FLOOR,zhc,b%XHC_FLOOR,ztc,b%XTC_FLOOR) 
  CALL deallocate_thermal_work
!
!
!     TCOOL_TARGET
!     -----------
  IF (dtb%LDATA_TCOOL_TARGET) THEN
    b%XTCOOL_TARGET=dtb%XPAR_TCOOL_TARGET
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,ptcool_target=b%XTCOOL_TARGET)
  ELSE
    CALL av_pgd(dtco, b%XTCOOL_TARGET ,top%XCOVER ,xdata_tcool_target(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     THEAT_TARGET
!     -----------
  IF (dtb%LDATA_THEAT_TARGET) THEN
    b%XTHEAT_TARGET=dtb%XPAR_THEAT_TARGET
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,ptheat_target=b%XTHEAT_TARGET)
  ELSE
    CALL av_pgd(dtco, b%XTHEAT_TARGET ,top%XCOVER ,xdata_theat_target(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     -----------
  IF (dtb%LDATA_F_WASTE_CAN) THEN
    b%XF_WASTE_CAN=dtb%XPAR_F_WASTE_CAN
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pf_waste_can=b%XF_WASTE_CAN)
  ELSE
    CALL av_pgd(dtco, b%XF_WASTE_CAN ,top%XCOVER ,xdata_f_waste_can(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     EFF_HEAT
!     -----------
  IF (dtb%LDATA_EFF_HEAT) THEN
    b%XEFF_HEAT=dtb%XPAR_EFF_HEAT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,peff_heat=b%XEFF_HEAT)
  ELSE
    CALL av_pgd(dtco, b%XEFF_HEAT ,top%XCOVER ,xdata_eff_heat(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     QIN
!     -----------
  IF (dtb%LDATA_QIN) THEN
    b%XQIN=dtb%XPAR_QIN
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,pqin=b%XQIN)
  ELSE
    CALL av_pgd(dtco, b%XQIN ,top%XCOVER ,xdata_qin(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     QIN_FRAD
!     -----------
  IF (dtb%LDATA_QIN_FRAD) THEN
    b%XQIN_FRAD=dtb%XPAR_QIN_FRAD
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,pqin_frad=b%XQIN_FRAD)
  ELSE
    CALL av_pgd(dtco, b%XQIN_FRAD ,top%XCOVER ,xdata_qin_frad(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     SHGC
!     -----------
  IF (dtb%LDATA_SHGC) THEN
    b%XSHGC=dtb%XPAR_SHGC
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pshgc=b%XSHGC)
  ELSE
    CALL av_pgd(dtco, b%XSHGC ,top%XCOVER ,xdata_shgc(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     U_WIN
!     -----------
  IF (dtb%LDATA_U_WIN) THEN
    b%XU_WIN=dtb%XPAR_U_WIN
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pu_win=b%XU_WIN)
  ELSE
    CALL av_pgd(dtco, b%XU_WIN ,top%XCOVER ,xdata_u_win(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     GR
!     -----------
  IF (dtb%LDATA_GR) THEN
    b%XGR=dtb%XPAR_GR
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pgr=b%XGR)
  ELSE
    CALL av_pgd(dtco, b%XGR ,top%XCOVER ,xdata_gr(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     SHGC_SH
!     -----------
  IF (dtb%LDATA_SHGC_SH) THEN
    b%XSHGC_SH=dtb%XPAR_SHGC_SH
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,pshgc_sh=b%XSHGC_SH)
  ELSE
    CALL av_pgd(dtco, b%XSHGC_SH ,top%XCOVER ,xdata_shgc_sh(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     FLOOR_HEIGHT
!     -----------
  IF (dtb%LDATA_FLOOR_HEIGHT) THEN
    b%XFLOOR_HEIGHT=dtb%XPAR_FLOOR_HEIGHT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pfloor_height=b%XFLOOR_HEIGHT)
  ELSE
    CALL av_pgd(dtco, b%XFLOOR_HEIGHT ,top%XCOVER ,xdata_floor_height(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     INF
!     -----------
  IF (dtb%LDATA_INF) THEN
    b%XINF=dtb%XPAR_INF
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pinf=b%XINF)
  ELSE
    CALL av_pgd(dtco, b%XINF ,top%XCOVER ,xdata_inf(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     F_WATER_COND
!     -----------
  IF (dtb%LDATA_F_WATER_COND) THEN
    b%XF_WATER_COND=dtb%XPAR_F_WATER_COND
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pf_water_cond=b%XF_WATER_COND)
  ELSE
    CALL av_pgd(dtco, b%XF_WATER_COND ,top%XCOVER ,xdata_f_water_cond(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     QIN_FLAT
!     -----------
  IF (dtb%LDATA_QIN_FLAT) THEN
    b%XQIN_FLAT=dtb%XPAR_QIN_FLAT
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,pqin_flat=b%XQIN_FLAT)
  ELSE
    CALL av_pgd(dtco, b%XQIN_FLAT ,top%XCOVER ,xdata_qin_flat(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     HR_TARGET
!     -----------
  IF (dtb%LDATA_HR_TARGET) THEN
    b%XHR_TARGET=dtb%XPAR_HR_TARGET
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,phr_target=b%XHR_TARGET)
  ELSE
    CALL av_pgd(dtco, b%XHR_TARGET ,top%XCOVER ,xdata_hr_target(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     V_VENT
!     -----------
  IF (dtb%LDATA_V_VENT) THEN
    b%XV_VENT=dtb%XPAR_V_VENT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pv_vent=b%XV_VENT)
  ELSE
    CALL av_pgd(dtco, b%XV_VENT ,top%XCOVER ,xdata_v_vent(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     CAP_SYS_HEAT
!     -----------
  IF (dtb%LDATA_CAP_SYS_HEAT) THEN
    b%XCAP_SYS_HEAT=dtb%XPAR_CAP_SYS_HEAT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pcap_sys_heat=b%XCAP_SYS_HEAT)
  ELSE
    CALL av_pgd(dtco, b%XCAP_SYS_HEAT ,top%XCOVER ,xdata_cap_sys_heat(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     CAP_SYS_RAT
!     -----------
  IF (dtb%LDATA_CAP_SYS_RAT) THEN
    b%XCAP_SYS_RAT=dtb%XPAR_CAP_SYS_RAT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pcap_sys_rat=b%XCAP_SYS_RAT)
  ELSE
    CALL av_pgd(dtco, b%XCAP_SYS_RAT ,top%XCOVER ,xdata_cap_sys_rat(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     T_ADP
!     -----------
  IF (dtb%LDATA_T_ADP) THEN
    b%XT_ADP=dtb%XPAR_T_ADP
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pt_adp=b%XT_ADP)
  ELSE
    CALL av_pgd(dtco, b%XT_ADP ,top%XCOVER ,xdata_t_adp(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     M_SYS_RAT
!     -----------
  IF (dtb%LDATA_M_SYS_RAT) THEN
    b%XM_SYS_RAT=dtb%XPAR_M_SYS_RAT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pm_sys_rat=b%XM_SYS_RAT)
  ELSE
    CALL av_pgd(dtco, b%XM_SYS_RAT ,top%XCOVER ,xdata_m_sys_rat(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     COP_RAT
!     -----------
  IF (dtb%LDATA_COP_RAT) THEN
    b%XCOP_RAT=dtb%XPAR_COP_RAT
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pcop_rat=b%XCOP_RAT)
  ELSE
    CALL av_pgd(dtco, b%XCOP_RAT ,top%XCOVER ,xdata_cop_rat(:),'TWN','ARI',top%LCOVER)
  ENDIF   
!
!     T_SIZE_MAX
!     -----------
  IF (dtb%LDATA_T_SIZE_MAX) THEN
    b%XT_SIZE_MAX=dtb%XPAR_T_SIZE_MAX
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pt_size_max=b%XT_SIZE_MAX)
  ELSE
    CALL av_pgd(dtco, b%XT_SIZE_MAX ,top%XCOVER ,xdata_t_size_max(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     T_SIZE_MIN
!     -----------
  IF (dtb%LDATA_T_SIZE_MIN) THEN
    b%XT_SIZE_MIN=dtb%XPAR_T_SIZE_MIN
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pt_size_min=b%XT_SIZE_MIN)
  ELSE
    CALL av_pgd(dtco, b%XT_SIZE_MIN ,top%XCOVER ,xdata_t_size_min(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     SHADE
!     -----
  IF (dtb%LDATA_SHADE) THEN
    b%XSHADE(:) = dtb%XPAR_SHADE(:)
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,pshade=b%XSHADE)
  ELSE
    CALL av_pgd(dtco, b%XSHADE,top%XCOVER ,xdata_shade(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
!     NATVENT
!     -------
  IF (dtb%LDATA_NATVENT) THEN
    b%XNATVENT(:) = dtb%XPAR_NATVENT(:)
  ELSEIF (dtt%LDATA_USETYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_USETYPE,pnatvent=b%XNATVENT)
  ELSE
    CALL av_pgd(dtco, b%XNATVENT,top%XCOVER,xdata_natvent(:),'TWN','ARI',top%LCOVER)
  ENDIF  
!
ENDIF
!
IF (PRESENT(tpn)) THEN
!
!     SOLAR PANELS
!     ------------
  IF (dtt%LDATA_EMIS_PANEL) THEN
    tpn%XEMIS_PANEL=dtt%XPAR_EMIS_PANEL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pemis_panel=tpn%XEMIS_PANEL)
  ELSE
    CALL av_pgd(dtco, tpn%XEMIS_PANEL ,top%XCOVER ,xdata_emis_panel(:),'BLD','ARI',top%LCOVER)
  ENDIF  
!
  IF (dtt%LDATA_ALB_PANEL) THEN
    tpn%XALB_PANEL=dtt%XPAR_ALB_PANEL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,palb_panel=tpn%XALB_PANEL)
  ELSE
    CALL av_pgd(dtco, tpn%XALB_PANEL ,top%XCOVER ,xdata_alb_panel(:),'BLD','ARI',top%LCOVER)
  ENDIF  
!
  IF (dtt%LDATA_EFF_PANEL) THEN
    tpn%XEFF_PANEL=dtt%XPAR_EFF_PANEL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,peff_panel=tpn%XEFF_PANEL)
  ELSE
    CALL av_pgd(dtco, tpn%XEFF_PANEL ,top%XCOVER ,xdata_eff_panel(:),'BLD','ARI',top%LCOVER)
  ENDIF  
!
  IF (dtt%LDATA_FRAC_PANEL) THEN
    tpn%XFRAC_PANEL=dtt%XPAR_FRAC_PANEL
  ELSEIF (dtt%LDATA_BLDTYPE) THEN
    CALL ini_data_param_teb(bdd, dtt%NPAR_BLDCODE,pfrac_panel=tpn%XFRAC_PANEL)
  ELSE
    CALL av_pgd(dtco, tpn%XFRAC_PANEL ,top%XCOVER ,xdata_frac_panel(:),'BLD','ARI',top%LCOVER)
  ENDIF  
!
ENDIF
!
IF (lhook) CALL dr_hook('CONVERT_PATCH_TEB',1,zhook_handle)
!
!-------------------------------------------------------------------------------
CONTAINS
!
SUBROUTINE allocate_thermal_work(KPAR)
INTEGER, INTENT(IN) :: KPAR
ALLOCATE(zd(il,kpar))
ALLOCATE(zhc(il,kpar))
ALLOCATE(ztc(il,kpar))
END SUBROUTINE allocate_thermal_work
SUBROUTINE deallocate_thermal_work
DEALLOCATE(zd)
DEALLOCATE(zhc)
DEALLOCATE(ztc)
END SUBROUTINE deallocate_thermal_work
!-------------------------------------------------------------------------------
!
END SUBROUTINE convert_patch_teb