SURFEX v7.3
General documentation of Surfex
 All Classes Files Functions Variables Typedefs
/home/dasprezs/EXPORT_v7_3/src/SURFEX/facade_e_budget.F90
Go to the documentation of this file.
00001 SUBROUTINE FACADE_E_BUDGET(HWALL_OPT, HBEM,                                    &
00002                            PT_WALL_A, PT_WALL_B, PTSTEP,                       &
00003                            PHC_WALL, PTC_WALL, PD_WALL,                        &
00004                            PDN_ROAD, PRHOA, PAC_WALL,                          &
00005                            PAC_BLD,                                            &
00006                            PTI_BLD, PLW_RAD, PPS, PEXNS,                       &
00007                            PABS_SW_WALL_A, PABS_SW_WALL_B, PT_CANYON,          &
00008                            PTS_ROAD, PTSNOW_ROAD, PTS_GARDEN, PU_WIN, PGR,     &
00009                            PABS_SW_WIN, PF_WIN_MASS, PF_WIN_WALL,              &
00010                            PF_WIN_FLOOR, PTS_MASS, PTS_FLOOR, PEMIS_WALL,      &
00011                            PLW_WA_TO_WB,                                       &
00012                            PLW_R_TO_WA, PLW_R_TO_WB, PLW_G_TO_WA, PLW_G_TO_WB, &
00013                            PLW_S_TO_WA, PLW_S_TO_WB, PLW_NR_TO_WA,PLW_NR_TO_WB,&
00014                            PLW_WIN_TO_WA, PLW_WIN_TO_WB,                       &
00015                            PLW_S_TO_WIN, PLW_WA_TO_WIN, PLW_WB_TO_WIN,         &
00016                            PLW_R_TO_WIN, PLW_G_TO_WIN, PLW_NR_TO_WIN, PT_WIN1, &
00017                            PFLX_BLD_WALL_A, PDQS_WALL_A,                       &
00018                            PFLX_BLD_WALL_B, PDQS_WALL_B,                       &
00019                            PT_WIN2, PEMIT_LW_FAC,                              &
00020                            PABS_LW_WALL_A, PABS_LW_WALL_B, PABS_LW_WIN,        &
00021                            PH_WALL_A, PH_WALL_B,                               &
00022                            PIMB_WALL, PF_WALL_MASS, PF_WALL_FLOOR,             &
00023                            PF_WALL_WIN, PRADHT_IN, PRAD_ROOF_WALL,             &
00024                            PRAD_ROOF_WIN, PRAD_WALL_FLOOR, PRAD_WALL_MASS,     &
00025                            PRAD_WIN_FLOOR, PRAD_WIN_MASS, PCONV_WALL_BLD,      &
00026                            PCONV_WIN_BLD, PAC_WIN,  PLOAD_IN_WALL, PLOAD_IN_WIN )
00027 !
00028 !   ##########################################################################
00029 !
00030 !!****  *FACADE_E_BUDGET*  
00031 !!
00032 !!    PURPOSE
00033 !!    -------
00034 !
00035 !     Computes the evolution of wall and window temperature from
00036 !     wall_layer_e_budget et window_e_budget
00037 !         
00038 !     
00039 !!**  METHOD
00040 !     ------
00041 !
00042 !!    EXTERNAL
00043 !!    --------
00044 !!
00045 !!
00046 !!    IMPLICIT ARGUMENTS
00047 !!    ------------------
00048 !!
00049 !!    MODD_CST
00050 !!
00051 !!      
00052 !!    REFERENCE
00053 !!    ---------
00054 !!
00055 !!      
00056 !!    AUTHOR
00057 !!    ------
00058 !!
00059 !!      G. Pigeon           * Meteo-France *
00060 !!
00061 !!    MODIFICATIONS
00062 !!    -------------
00063 !!      Original    10/11/2011 
00064 !!      G. Pigeon     /09/2012: new conv. coef for indoor/outdoor 
00065 !!      G. Pigeon     /10/2012: new arg. : solar heat load for indoor
00066 !-------------------------------------------------------------------------------
00067 !
00068 !*       0.     DECLARATIONS
00069 !               ------------
00070 !
00071 USE MODD_SURF_PAR, ONLY : XUNDEF
00072 USE MODD_CSTS, ONLY : XCPD
00073 !
00074 USE MODI_WALL_LAYER_E_BUDGET
00075 USE MODI_WINDOW_E_BUDGET
00076 USE YOMHOOK   ,ONLY : LHOOK,   DR_HOOK
00077 USE PARKIND1  ,ONLY : JPRB
00078 !
00079 IMPLICIT NONE
00080 !
00081 !*      0.1    declarations of arguments
00082 !
00083  CHARACTER(LEN=4), INTENT(IN)      :: HWALL_OPT     ! TEB option for walls
00084                                                    ! 'UNIF' : uniform walls
00085                                                    ! 'TWO ' : two different opposite walls
00086  CHARACTER(LEN=3), INTENT(IN)      :: HBEM         ! Building Energy model 'DEF' or 'BEM'
00087 REAL, DIMENSION(:,:), INTENT(INOUT) :: PT_WALL_A  ! wall A layers temperatures
00088 REAL, DIMENSION(:,:), INTENT(INOUT) :: PT_WALL_B  ! wall B layers temperatures
00089 REAL,               INTENT(IN)    :: PTSTEP       ! time step
00090 REAL, DIMENSION(:,:), INTENT(IN)  :: PHC_WALL     ! heat capacity for wall layers
00091 REAL, DIMENSION(:,:), INTENT(IN)  :: PTC_WALL     ! thermal conductivity for wall layers
00092 REAL, DIMENSION(:,:), INTENT(IN)  :: PD_WALL      ! depth of wall layers
00093 REAL, DIMENSION(:), INTENT(IN)    :: PDN_ROAD     ! snow-covered fraction on roads
00094 REAL, DIMENSION(:), INTENT(IN)    :: PRHOA        ! rho
00095 REAL, DIMENSION(:), INTENT(IN)    :: PAC_WALL     ! aerodynamical conductance
00096 !                                                 ! between wall and canyon
00097 REAL, DIMENSION(:), INTENT(IN)    :: PAC_BLD      ! aerodynamical conductance
00098                                                   ! inside the building itself
00099 REAL, DIMENSION(:), INTENT(IN)    :: PTI_BLD      ! inside building temperature
00100 REAL, DIMENSION(:), INTENT(IN)    :: PLW_RAD      ! atmospheric infrared radiation
00101 REAL, DIMENSION(:), INTENT(IN)    :: PPS          ! pressure at the surface
00102 REAL, DIMENSION(:), INTENT(IN)    :: PEXNS        ! surface Exner function
00103 REAL, DIMENSION(:), INTENT(IN)    :: PABS_SW_WALL_A! absorbed solar radiation
00104 REAL, DIMENSION(:), INTENT(IN)    :: PABS_SW_WALL_B! absorbed solar radiation
00105 REAL, DIMENSION(:), INTENT(IN)    :: PT_CANYON    ! air canyon temperature
00106 REAL, DIMENSION(:), INTENT(IN)    :: PTS_ROAD     ! road surface temperature
00107 REAL, DIMENSION(:), INTENT(IN)    :: PTSNOW_ROAD  ! road snow temperature
00108 REAL, DIMENSION(:), INTENT(IN)    :: PTS_GARDEN   ! green area surface temperature
00109 REAL, DIMENSION(:), INTENT(IN)    :: PU_WIN       ! window U-factor [W m-2 K-1]
00110 REAL, DIMENSION(:), INTENT(IN)    :: PGR          ! glazing ratio
00111 REAL, DIMENSION(:), INTENT(IN)    :: PABS_SW_WIN  ! window absorbed shortwave radiation [W m-2] 
00112 REAL, DIMENSION(:),   INTENT(IN)  :: PF_WIN_MASS  ! window mass view factor
00113 REAL, DIMENSION(:),   INTENT(IN)  :: PF_WIN_WALL  ! window wall view factor
00114 REAL, DIMENSION(:),   INTENT(IN)  :: PF_WIN_FLOOR ! window floor view factor
00115 REAL, DIMENSION(:), INTENT(IN)    :: PTS_MASS     ! surface mass temperature  [K]
00116 REAL, DIMENSION(:), INTENT(IN)    :: PTS_FLOOR    ! floor layers temperatures [K]
00117 REAL, DIMENSION(:), INTENT(IN)    :: PEMIS_WALL   ! wall emissivity
00118 !
00119 REAL, DIMENSION(:), INTENT(IN)    :: PLW_WA_TO_WB ! LW interactions wall  -> wall (both versions)
00120 REAL, DIMENSION(:), INTENT(IN)    :: PLW_R_TO_WA  ! LW interactions road -> wall ; DEF formulation 
00121 REAL, DIMENSION(:), INTENT(IN)    :: PLW_R_TO_WB  ! LW interactions road -> wall ; DEF formulation 
00122 REAL, DIMENSION(:), INTENT(IN)    :: PLW_G_TO_WA  ! LW interactions garden -> wall ; DEF formulation
00123 REAL, DIMENSION(:), INTENT(IN)    :: PLW_G_TO_WB  ! LW interactions garden -> wall ; DEF formulation
00124 REAL, DIMENSION(:), INTENT(IN)    :: PLW_S_TO_WA  ! LW interactions sky   -> wall 
00125 REAL, DIMENSION(:), INTENT(IN)    :: PLW_S_TO_WB  ! LW interactions sky   -> wall 
00126 REAL, DIMENSION(:), INTENT(IN)    :: PLW_NR_TO_WA ! LW interactions road(snow) -> wall 
00127 REAL, DIMENSION(:), INTENT(IN)    :: PLW_NR_TO_WB ! LW interactions road(snow) -> wall 
00128 REAL, DIMENSION(:), INTENT(IN)    :: PLW_WA_TO_WIN! Radiative heat trasfer coeff wall-window 
00129                                                   ! [W K-1 m-2]
00130 REAL, DIMENSION(:), INTENT(IN)    :: PLW_WB_TO_WIN! Radiative heat trasfer coeff wall-window 
00131                                                   ! [W K-1 m-2]
00132 REAL, DIMENSION(:), INTENT(IN)    :: PLW_S_TO_WIN ! Radiative heat trasfer coeff window-sky 
00133                                                   ! [W K-1 m-2]
00134 REAL, DIMENSION(:), INTENT(IN)    :: PLW_WIN_TO_WA! Radiative heat trasfer coeff window-wall
00135                                                   ! [W K-1 m-2] 
00136 REAL, DIMENSION(:), INTENT(IN)    :: PLW_WIN_TO_WB! Radiative heat trasfer coeff window-wall
00137                                                   ! [W K-1 m-2] 
00138 REAL, DIMENSION(:), INTENT(IN)    :: PLW_R_TO_WIN ! Radiative heat trasfer coeff window-road 
00139                                                   ! [W K-1 m-2]
00140 REAL, DIMENSION(:), INTENT(IN)    :: PLW_G_TO_WIN ! Radiative heat trasfer coeff window-garden
00141 REAL, DIMENSION(:), INTENT(IN)    :: PLW_NR_TO_WIN! Radiative heat trasfer coeff window-snow
00142 !
00143 REAL, DIMENSION(:), INTENT(INOUT) :: PT_WIN1      ! outdoor window temperature [K]
00144 REAL, DIMENSION(:), INTENT(OUT)   :: PFLX_BLD_WALL_A! flux from bld to wall
00145 REAL, DIMENSION(:), INTENT(INOUT) :: PDQS_WALL_A    ! heat storage inside the wall 
00146 REAL, DIMENSION(:), INTENT(OUT)   :: PFLX_BLD_WALL_B! flux from bld to wall
00147 REAL, DIMENSION(:), INTENT(INOUT) :: PDQS_WALL_B    ! heat storage inside the wall 
00148 REAL, DIMENSION(:), INTENT(INOUT) :: PT_WIN2      ! indoor window temperature [K]
00149 REAL, DIMENSION(:), INTENT(OUT)   :: PEMIT_LW_FAC ! LW flux emitted by the facade (W/m² of facade)
00150 REAL, DIMENSION(:), INTENT(OUT)   :: PABS_LW_WALL_A ! absorbed infrared rad.
00151 REAL, DIMENSION(:), INTENT(OUT)   :: PABS_LW_WALL_B ! absorbed infrared rad.
00152 REAL, DIMENSION(:), INTENT(OUT)   :: PABS_LW_WIN  ! window absorbed shortwave radiation [W m-2] 
00153 REAL, DIMENSION(:), INTENT(OUT)   :: PH_WALL_A    ! Sensible heat flux from wall to air [W/m²(wall)]
00154                                                   ! wall = facade - glazing
00155 REAL, DIMENSION(:), INTENT(OUT)   :: PH_WALL_B    ! Sensible heat flux from wall to air [W/m²(wall)]
00156                                                   ! wall = facade - glazing
00157 REAL, DIMENSION(:), INTENT(OUT)   :: PIMB_WALL    ! wall residual energy imbalance 
00158                                                   ! for verification [W m-2]
00159 REAL, DIMENSION(:),   INTENT(IN)  :: PF_WALL_FLOOR ! View factor wall-floor
00160 REAL, DIMENSION(:),   INTENT(IN)  :: PF_WALL_MASS  ! View factor wall-mass
00161 REAL, DIMENSION(:),   INTENT(IN)  :: PF_WALL_WIN   ! View factor wall-win
00162 REAL, DIMENSION(:),   INTENT(IN)  :: PRADHT_IN     ! Indoor radiant heat transfer coefficient
00163                                                     ! [W K-1 m-2]
00164 REAL, DIMENSION(:), INTENT(IN)    :: PRAD_ROOF_WALL ! rad. fluxes from roof to wall [W m-2(roof)]
00165 REAL, DIMENSION(:), INTENT(IN)    :: PRAD_ROOF_WIN  ! rad. fluxes from roof to win [W m-2(roof)]
00166 REAL, DIMENSION(:), INTENT(OUT)   :: PRAD_WALL_FLOOR! rad. fluxes from wall to floor [W m-2(wall)]
00167 REAL, DIMENSION(:), INTENT(OUT)   :: PRAD_WALL_MASS ! rad. fluxes from wall to mass [W m-2(wall)]
00168 REAL, DIMENSION(:), INTENT(OUT)   :: PRAD_WIN_FLOOR ! rad. fluxes from window to floor [W m-2(win)]
00169 REAL, DIMENSION(:), INTENT(OUT)   :: PRAD_WIN_MASS  ! rad. fluxes from window to mass [W m-2(win)]
00170 REAL, DIMENSION(:), INTENT(OUT)   :: PCONV_WALL_BLD ! conv. fluxes from wall to bld [W m-2(wall)]
00171 REAL, DIMENSION(:), INTENT(OUT)   :: PCONV_WIN_BLD  ! conv. fluxes from window to bld [W m-2(wind.)]
00172 REAL, DIMENSION(:), INTENT(IN)    :: PAC_WIN        ! window aerodynamic conductance
00173 !
00174 REAL, DIMENSION(:), INTENT(IN)    :: PLOAD_IN_WALL  ! solar + inter. heat gains W/m² [wall]
00175 REAL, DIMENSION(:), INTENT(IN)    :: PLOAD_IN_WIN   ! solar + inter. heat gains W/m² [win]
00176 !*      0.2    declarations of local variables
00177 !
00178 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZRAD_WALL_WIN   ! rad. fluxes between averaged wall and win [W m-2(wall)]
00179 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZRAD_WALL_A_WIN ! rad. fluxes between one wall and win [W m-2(wall)]
00180 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZEMIS_WIN     ! window emissivity
00181 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZEMIT_LW_WIN  ! Longwave radiation emitted by the window [W m-2(window)]
00182 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZEMIT_LW_WALL_A ! Longwave radiation emitted by the wall [W m-2(wall)]
00183 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZEMIT_LW_WALL_B ! Longwave radiation emitted by the wall [W m-2(wall)]
00184 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZTS_WALL_A    ! surface temperature of wall A at previous time-step
00185 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZTS_WALL_B    ! surface temperature of wall B at previous time-step
00186 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZTI_WALL_A    ! internal temperature of wall A at previous time-step
00187 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZTI_WALL_B    ! internal temperature of wall B at previous time-step
00188 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZTS_WALL      ! surface temperature of averaged wall at new time-step
00189 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZIMB_WALL     ! wall residual energy imbalance for verification [W m-2]
00190 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZRAD_WALL_FLOOR! rad. fluxes from wall to floor [W m-2(wall)]
00191 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZRAD_WALL_MASS ! rad. fluxes from wall to mass [W m-2(wall)]
00192 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZCONV_WALL_BLD ! conv. fluxes from wall to bld [W m-2(wall)]
00193 REAL, DIMENSION(SIZE(PTI_BLD)) :: ZLW_W_TO_WIN   ! Radiative heat trasfer coeff wall-window [W K-1 m-2]
00194 INTEGER                        :: IWALL
00195 REAL(KIND=JPRB) :: ZHOOK_HANDLE
00196 !-------------------------------------------------------------------------------
00197 IF (LHOOK) CALL DR_HOOK('FACADE_E_BUDGET',0,ZHOOK_HANDLE)
00198 !-------------------------------------------------------------------------------
00199 !
00200 !* surface temperature of the opposite wall
00201 !  ----------------------------------------
00202 !
00203 ZTS_WALL_A(:) = PT_WALL_A(:,1)
00204 ZTS_WALL_B(:) = PT_WALL_B(:,1)
00205 IWALL = SIZE(PT_WALL_A,2)
00206 ZTI_WALL_A(:) = PT_WALL_A(:,IWALL)
00207 ZTI_WALL_B(:) = PT_WALL_B(:,IWALL)
00208 !
00209 ! *Convection heat transfer coefficients [W m-2 K-1] from EP Engineering Reference
00210 ! --------------------------------------------------------------------------------
00211 !
00212 !
00213 ! *opaque balance of the facade for wall A
00214 !  ---------------------------------------
00215 !
00216  CALL WALL_LAYER_E_BUDGET(HBEM, PT_WALL_A, ZTS_WALL_B, ZTI_WALL_B, PTSTEP,                     &
00217                          PHC_WALL, PTC_WALL, PD_WALL, PDN_ROAD, PRHOA, PAC_WALL,  &
00218                          PAC_BLD, PTI_BLD, PLW_RAD, PPS, PEXNS, PABS_SW_WALL_A, PT_CANYON,    &
00219                          PTS_ROAD, PTSNOW_ROAD, PTS_GARDEN, PTS_MASS, PTS_FLOOR, PEMIS_WALL,  &
00220                          PLW_WA_TO_WB, PLW_R_TO_WA, PLW_G_TO_WA, PLW_NR_TO_WA,                &
00221                          PLW_WIN_TO_WA, PLW_S_TO_WA,                                          &
00222                          PT_WIN1, PFLX_BLD_WALL_A, PDQS_WALL_A, PT_WIN2, PABS_LW_WALL_A,      &
00223                          ZEMIT_LW_WALL_A, PH_WALL_A, ZIMB_WALL, PF_WALL_MASS,                 &
00224                          PF_WALL_FLOOR, PF_WALL_WIN, PRADHT_IN, PRAD_ROOF_WALL,               &
00225                          ZRAD_WALL_A_WIN, ZRAD_WALL_FLOOR, ZRAD_WALL_MASS, ZCONV_WALL_BLD,    &
00226                          PLOAD_IN_WALL)
00227 !
00228 PIMB_WALL       = 0.5 * ZIMB_WALL
00229 PRAD_WALL_FLOOR = 0.5 * ZRAD_WALL_FLOOR
00230 PRAD_WALL_MASS  = 0.5 * ZRAD_WALL_MASS 
00231 PCONV_WALL_BLD  = 0.5 * ZCONV_WALL_BLD 
00232 ZRAD_WALL_WIN   = 0.5 * ZRAD_WALL_A_WIN
00233 !
00234 ! *opaque balance of the facade for wall B
00235 !  ---------------------------------------
00236 !
00237 IF (HWALL_OPT/='UNIF') THEN
00238  CALL WALL_LAYER_E_BUDGET(HBEM, PT_WALL_B, ZTS_WALL_A, ZTI_WALL_A, PTSTEP,                     &
00239                          PHC_WALL, PTC_WALL, PD_WALL, PDN_ROAD, PRHOA, PAC_WALL,              &
00240                          PAC_BLD, PTI_BLD, PLW_RAD, PPS, PEXNS, PABS_SW_WALL_B, PT_CANYON,    &
00241                          PTS_ROAD, PTSNOW_ROAD, PTS_GARDEN, PTS_MASS, PTS_FLOOR, PEMIS_WALL,  &
00242                          PLW_WA_TO_WB, PLW_R_TO_WB, PLW_G_TO_WB, PLW_NR_TO_WB,                &
00243                          PLW_WIN_TO_WB, PLW_S_TO_WB,                                          &
00244                          PT_WIN1, PFLX_BLD_WALL_B, PDQS_WALL_B, PT_WIN2, PABS_LW_WALL_B,      &
00245                          ZEMIT_LW_WALL_B, PH_WALL_B, ZIMB_WALL, PF_WALL_MASS,                 &
00246                          PF_WALL_FLOOR, PF_WALL_WIN, PRADHT_IN, PRAD_ROOF_WALL,               &
00247                          ZRAD_WALL_A_WIN, ZRAD_WALL_FLOOR, ZRAD_WALL_MASS, ZCONV_WALL_BLD,    & 
00248                          PLOAD_IN_WALL                                                        )
00249 ELSE
00250   PT_WALL_B        = PT_WALL_A
00251   PH_WALL_B        = PH_WALL_A
00252   PABS_LW_WALL_B   = PABS_LW_WALL_A
00253   PDQS_WALL_B      = PDQS_WALL_A
00254   PFLX_BLD_WALL_B  = PFLX_BLD_WALL_A
00255   ZEMIT_LW_WALL_B  = ZEMIT_LW_WALL_A
00256 END IF
00257 !
00258 PIMB_WALL       = PIMB_WALL       + 0.5 * ZIMB_WALL
00259 PRAD_WALL_FLOOR = PRAD_WALL_FLOOR + 0.5 * ZRAD_WALL_FLOOR
00260 PRAD_WALL_MASS  = PRAD_WALL_MASS  + 0.5 * ZRAD_WALL_MASS 
00261 PCONV_WALL_BLD  = PCONV_WALL_BLD  + 0.5 * ZCONV_WALL_BLD 
00262 ZRAD_WALL_WIN   = ZRAD_WALL_WIN   + 0.5 * ZRAD_WALL_A_WIN
00263 !
00264 !-------------------------------------------------------------------------------
00265 !
00266 ! *Energy Balance for windows (averaged on both walls)
00267 !  ---------------------------------------------------
00268 !
00269 IF (HBEM == 'DEF') ZEMIT_LW_WIN(:) = 0.
00270 IF (HBEM == 'BEM') THEN
00271         ZEMIS_WIN(:) = 0.84 ! From Engineering Reference documentation p219 
00272                             ! (see http://apps1.eere.energy.gov/buildings/energyplus/
00273    ZLW_W_TO_WIN(:) = PLW_WA_TO_WIN(:) + PLW_WB_TO_WIN(:)
00274    ZTS_WALL(:) = XUNDEF
00275    WHERE (ZLW_W_TO_WIN(:)>0.) &
00276    ZTS_WALL(:) = ( PLW_WA_TO_WIN(:)*PT_WALL_A(:,1)+PLW_WB_TO_WIN(:)*PT_WALL_B(:,1) ) / ZLW_W_TO_WIN(:)
00277    CALL WINDOW_E_BUDGET(PU_WIN, ZEMIS_WIN, PF_WIN_FLOOR, PF_WIN_WALL, PF_WIN_MASS,                &
00278                         ZLW_W_TO_WIN, PLW_R_TO_WIN, PLW_G_TO_WIN, PLW_NR_TO_WIN, PLW_S_TO_WIN,    &
00279                         PRAD_ROOF_WIN, ZRAD_WALL_WIN, PABS_SW_WIN, PLW_RAD, PAC_WIN,              &
00280                         PRADHT_IN, PTI_BLD, PTS_MASS, PTS_FLOOR, PRHOA, PDN_ROAD,                 &
00281                         PT_CANYON, ZTS_WALL, PTS_ROAD, PTSNOW_ROAD, PTS_GARDEN,                   &
00282                         PT_WIN1, PT_WIN2, PRAD_WIN_FLOOR, PRAD_WIN_MASS, PCONV_WIN_BLD,           &
00283                         ZEMIT_LW_WIN, PABS_LW_WIN, PLOAD_IN_WIN                                   )
00284 ENDIF
00285 
00286 !*        wall, and win emitted LW radiation on snow-free surfaces
00287 !         ----------------------------------------------------------------
00288 !
00289 !
00290 PEMIT_LW_FAC (:) = 0.5*(ZEMIT_LW_WALL_A(:)+ZEMIT_LW_WALL_B(:)) * (1-PGR(:)) +  ZEMIT_LW_WIN(:) * PGR(:)
00291 !
00292 !-------------------------------------------------------------------------------
00293 IF (LHOOK) CALL DR_HOOK('FACADE_E_BUDGET',1,ZHOOK_HANDLE)
00294 END SUBROUTINE FACADE_E_BUDGET