init_bemn.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 init_bem_n ( DGU, DTCO, UG, U, TM, &
                             kluout)
!     #############################################################
!
!!****  *INIT_BEM_n* - routine to initialize Building Energy Model
!!
!!    PURPOSE
!!    -------
!!
!!**  METHOD
!!    ------
!!
!!    EXTERNAL
!!    --------
!!
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!
!!    AUTHOR
!!    ------
!!      G. Pigeon   *Meteo France*
!!
!!    MODIFICATIONS
!!    -------------
!!      Original    04/2012
!-------------------------------------------------------------------------------
!
!*       0.    DECLARATIONS
!              ------------
!
!
USE modd_diag_surf_atm_n, ONLY : diag_surf_atm_t
USE modd_data_cover_n, ONLY : data_cover_t
USE modd_surf_atm_grid_n, ONLY : surf_atm_grid_t
USE modd_surf_atm_n, ONLY : surf_atm_t
USE modd_surfex_n, ONLY : teb_model_t
!
USE modi_convert_patch_teb
USE modi_window_data
USE modi_hvac_autosize
USE modi_bem_morpho
USE modi_stores_hvac_autosize
!
USE modi_abor1_sfx
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*       0.1   Declarations of arguments
!              -------------------------
!
!
TYPE(diag_surf_atm_t), INTENT(INOUT) :: dgu
TYPE(data_cover_t), INTENT(INOUT) :: dtco
TYPE(surf_atm_grid_t), INTENT(INOUT) :: ug
TYPE(surf_atm_t), INTENT(INOUT) :: u
TYPE(teb_model_t), INTENT(INOUT) :: tm
!
INTEGER, INTENT(IN) :: kluout ! logical unit of output listing
!
!*       0.2   Declarations of local variables
!              -------------------------------
!
INTEGER                         :: jj               ! counter
INTEGER                         :: ilu              ! sizes of TEB arrays
LOGICAL                         :: gprint           ! flag for warning prints in output file
!
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!         Initialisation for IO
!

IF (lhook) CALL dr_hook('INIT_BEM_N',0,zhook_handle)
!
!
!*       3.     Physiographic data fields from land cover:
!               -----------------------------------------
!
ilu = SIZE(tm%TOP%XCOVER,1)
IF (tm%TOP%CBEM=='DEF') ilu=0
!
ALLOCATE(tm%B%CUR%XHC_FLOOR    (ilu,tm%BOP%NFLOOR_LAYER))
ALLOCATE(tm%B%CUR%XTC_FLOOR    (ilu,tm%BOP%NFLOOR_LAYER))
ALLOCATE(tm%B%CUR%XD_FLOOR     (ilu,tm%BOP%NFLOOR_LAYER))
!
ALLOCATE(tm%B%CUR%XTCOOL_TARGET(ilu))
ALLOCATE(tm%B%CUR%XTHEAT_TARGET(ilu))
ALLOCATE(tm%B%CUR%XEFF_HEAT    (ilu))
ALLOCATE(tm%B%CUR%XSHGC        (ilu))
ALLOCATE(tm%B%CUR%XQIN         (ilu))
ALLOCATE(tm%B%CUR%XQIN_FRAD    (ilu))
ALLOCATE(tm%B%CUR%XSHGC_SH     (ilu))
ALLOCATE(tm%B%CUR%XU_WIN       (ilu))
ALLOCATE(tm%B%CUR%XTRAN_WIN    (ilu))
ALLOCATE(tm%B%CUR%XFLOOR_HEIGHT(ilu))
ALLOCATE(tm%B%CUR%XINF         (ilu))
!
ALLOCATE(tm%B%CUR%XQIN_FLAT    (ilu))
ALLOCATE(tm%B%CUR%XHR_TARGET   (ilu))
ALLOCATE(tm%B%CUR%XV_VENT      (ilu))
ALLOCATE(tm%B%CUR%XCAP_SYS_HEAT(ilu))
ALLOCATE(tm%B%CUR%XCAP_SYS_RAT (ilu))
ALLOCATE(tm%B%CUR%XT_ADP       (ilu))
ALLOCATE(tm%B%CUR%XM_SYS_RAT   (ilu))
ALLOCATE(tm%B%CUR%XCOP_RAT     (ilu))
ALLOCATE(tm%B%CUR%XT_SIZE_MAX  (ilu))
ALLOCATE(tm%B%CUR%XT_SIZE_MIN  (ilu))
ALLOCATE(tm%B%CUR%XF_WATER_COND(ilu))
ALLOCATE(tm%B%CUR%CNATVENT     (ilu))
ALLOCATE(tm%B%CUR%XNATVENT     (ilu))
!
ALLOCATE(tm%B%CUR%XABS_WIN (ilu))
ALLOCATE(tm%B%CUR%XUGG_WIN (ilu))
ALLOCATE(tm%B%CUR%LSHADE   (ilu))
ALLOCATE(tm%B%CUR%XSHADE   (ilu))
ALLOCATE(tm%B%CUR%LSHAD_DAY(ilu))
ALLOCATE(tm%B%CUR%LNATVENT_NIGHT(ilu))
ALLOCATE(tm%B%CUR%XAUX_MAX    (ilu))
ALLOCATE(tm%B%CUR%XN_FLOOR(ilu))
ALLOCATE(tm%B%CUR%XGLAZ_O_BLD(ilu))
ALLOCATE(tm%B%CUR%XMASS_O_BLD(ilu))
ALLOCATE(tm%B%CUR%XFLOOR_HW_RATIO(ilu))
ALLOCATE(tm%B%CUR%XF_FLOOR_MASS(ilu))
ALLOCATE(tm%B%CUR%XF_FLOOR_WALL(ilu))
ALLOCATE(tm%B%CUR%XF_FLOOR_WIN(ilu))
ALLOCATE(tm%B%CUR%XF_FLOOR_ROOF(ilu))
ALLOCATE(tm%B%CUR%XF_WALL_FLOOR(ilu))
ALLOCATE(tm%B%CUR%XF_WALL_MASS(ilu))
ALLOCATE(tm%B%CUR%XF_WALL_WIN(ilu))
ALLOCATE(tm%B%CUR%XF_WIN_FLOOR(ilu))
ALLOCATE(tm%B%CUR%XF_WIN_MASS(ilu))
ALLOCATE(tm%B%CUR%XF_WIN_WALL(ilu))
ALLOCATE(tm%B%CUR%XF_WIN_WIN(ilu))
ALLOCATE(tm%B%CUR%XF_MASS_FLOOR(ilu))
ALLOCATE(tm%B%CUR%XF_MASS_WALL(ilu))
ALLOCATE(tm%B%CUR%XF_MASS_WIN(ilu))

SELECT CASE(tm%TOP%CBEM)
!----------
 CASE("DEF")
!-----------
   !parameters that needs to be 0 for calculation
   tm%B%CUR%XGR  (:)         = 0.
   tm%B%CUR%XF_WASTE_CAN(:)  = 0.
!----------
 CASE("BEM")
!----------

  tm%B%CUR%XAUX_MAX(:) = 5.
  CALL convert_patch_teb(tm%BDD, tm%DTB, dtco, tm%DTT, tm%TOP, &
                         tm%TOP%XCOVER,tm%TOP%LCOVER,0.,                 &
                      phc_floor=tm%B%CUR%XHC_FLOOR, ptc_floor=tm%B%CUR%XTC_FLOOR, &
                      pd_floor=tm%B%CUR%XD_FLOOR, ptcool_target=tm%B%CUR%XTCOOL_TARGET, &
                      ptheat_target=tm%B%CUR%XTHEAT_TARGET, pf_waste_can=tm%B%CUR%XF_WASTE_CAN, &
                      peff_heat=tm%B%CUR%XEFF_HEAT, pqin=tm%B%CUR%XQIN, pqin_frad=tm%B%CUR%XQIN_FRAD, &
                      pshgc=tm%B%CUR%XSHGC, pu_win=tm%B%CUR%XU_WIN, pgr=tm%B%CUR%XGR,  &
                      pshgc_sh=tm%B%CUR%XSHGC_SH, pfloor_height=tm%B%CUR%XFLOOR_HEIGHT, &
                      pinf=tm%B%CUR%XINF, pf_water_cond=tm%B%CUR%XF_WATER_COND, &
                      pqin_flat=tm%B%CUR%XQIN_FLAT, phr_target=tm%B%CUR%XHR_TARGET, &
                      pv_vent=tm%B%CUR%XV_VENT, pcap_sys_heat=tm%B%CUR%XCAP_SYS_HEAT, &
                      pcap_sys_rat=tm%B%CUR%XCAP_SYS_RAT, pt_adp=tm%B%CUR%XT_ADP, &
                      pm_sys_rat=tm%B%CUR%XM_SYS_RAT, pcop_rat=tm%B%CUR%XCOP_RAT,     &
                      pt_size_max=tm%B%CUR%XT_SIZE_MAX, pt_size_min=tm%B%CUR%XT_SIZE_MIN,    &
                      pshade=tm%B%CUR%XSHADE, pnatvent=tm%B%CUR%XNATVENT)
   !
   !
   ! *.     indoor relative surf. and view factors
   !        --------------------------------------
   !
   CALL bem_morpho(tm%T%CUR%XBLD, tm%T%CUR%XWALL_O_HOR, tm%T%CUR%XBLD_HEIGHT, &
                   tm%B%CUR%XFLOOR_HEIGHT, tm%B%CUR%XGR, tm%B%CUR%XN_FLOOR, tm%T%CUR%XWALL_O_BLD, &
                   tm%B%CUR%XGLAZ_O_BLD, tm%B%CUR%XMASS_O_BLD, tm%B%CUR%XFLOOR_HW_RATIO, &
                   tm%B%CUR%XF_FLOOR_MASS, tm%B%CUR%XF_FLOOR_WALL, tm%B%CUR%XF_FLOOR_WIN, &
                   tm%B%CUR%XF_FLOOR_ROOF, tm%B%CUR%XF_WALL_FLOOR, tm%B%CUR%XF_WALL_MASS, &
                   tm%B%CUR%XF_WALL_WIN, tm%B%CUR%XF_WIN_FLOOR, tm%B%CUR%XF_WIN_MASS, &
                   tm%B%CUR%XF_WIN_WALL,  tm%B%CUR%XF_MASS_FLOOR, tm%B%CUR%XF_MASS_WALL, &
                   tm%B%CUR%XF_MASS_WIN, tm%B%CUR%XF_WASTE_CAN, tm%B%CUR%XF_WIN_WIN      )
   !
   ! *.     Window optical and thermal data
   !        -------------------------------
   !
   CALL window_data(ilu, tm%B%CUR%XSHGC, tm%B%CUR%XU_WIN, tm%B%CUR%XALB_WIN, &
                    tm%B%CUR%XABS_WIN, tm%B%CUR%XUGG_WIN, tm%B%CUR%XTRAN_WIN)
   gprint = .false.
   DO jj=1,SIZE(tm%B%CUR%XSHADE)
      IF (tm%B%CUR%XSHADE(jj) >= 0.0 .AND. tm%B%CUR%XSHADE(jj) < 0.5) THEN
         tm%B%CUR%LSHADE(jj) = .false.
      ELSEIF (tm%B%CUR%XSHADE(jj) >= 0.5 .AND. tm%B%CUR%XSHADE(jj) <= 1.0) THEN
         tm%B%CUR%LSHADE(jj) = .true.
      ELSE
       gprint = .true.
       tm%B%CUR%LSHADE(jj) = .false.
      ENDIF
   ENDDO
   IF (gprint) WRITE(kluout,*) &
   'TEB-BEM : Error in specifying shading devices for at least one point, no shading device for these points'
   tm%B%CUR%LSHAD_DAY(:) = .false.
   !
   ! *.     Nocturnal surventilation
   !        ------------------------
   gprint = .false.
   DO jj=1,SIZE(tm%B%CUR%XNATVENT)
      IF (tm%B%CUR%XNATVENT(jj) >= 0.0 .AND. tm%B%CUR%XNATVENT(jj) < 0.5) THEN
        tm%B%CUR%CNATVENT(jj) = 'NONE'
      ELSEIF (tm%B%CUR%XNATVENT(jj) >= 0.5 .AND. tm%B%CUR%XNATVENT(jj) < 1.5) THEN
        tm%B%CUR%CNATVENT(jj) = 'MANU'
      ELSEIF (tm%B%CUR%XNATVENT(jj) >= 1.5 .AND. tm%B%CUR%XNATVENT(jj) <= 2.5) THEN
        tm%B%CUR%CNATVENT(jj) = 'AUTO'        
      ELSEIF (tm%B%CUR%XNATVENT(jj) >= 2.5 .AND. tm%B%CUR%XNATVENT(jj) <= 3.5) THEN
        tm%B%CUR%CNATVENT(jj) = 'MECH'        
      ELSE
        gprint = .true.
        tm%B%CUR%CNATVENT(jj) = 'NONE'        
      ENDIF
    ENDDO
    IF (gprint) WRITE(kluout,*) 'TEB-BEM : Chosen option for surventilation is not yet implemented; None venting is kept instead'

   tm%B%CUR%LNATVENT_NIGHT(:) = .false.
   !
END SELECT
!
!-------------------------------------------------------------------------------
!
!*       8.     Building HVAC automatic sizing:
!               -------------------------------  
IF (tm%TOP%CBEM=='BEM' .AND. tm%BOP%LAUTOSIZE) THEN
  CALL hvac_autosize(tm%BDD, dtco, dgu, tm%B, tm%BOP, ug, u, tm%TG, tm%T, tm%TOP, &
                     ilu,kluout)
  !* stores the real systems characteristics in physiographic data 
  !  for further use
  CALL stores_hvac_autosize(tm%B, tm%BOP, tm%DTB)
ENDIF
!
!-------------------------------------------------------------------------------
!
!
IF (lhook) CALL dr_hook('INIT_BEM_N',1,zhook_handle)
!
!
END SUBROUTINE init_bem_n