averaged_albedo_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 averaged_albedo_teb(HBEM, HROAD_DIR, HWALL_OPT,     &
                       pzenith, pazim,                               &
                       pbld, pgarden, proad_dir, proad, pfrac_gr,    &
                       pfrac_panel, palb_panel,                      &
                       pwall_o_hor, pcan_hw_ratio,                   &
                       palb_roof,                                    &
                       palb_road, psvf_road,                         &
                       palb_wall, psvf_wall,                         &
                       palb_garden, psvf_garden,                     &
                       palb_greenroof,                               &
                       tsnow_roof, tsnow_road,                       &
                       pgr, pshgc, pshgc_sh, pabs_win, palb_win,     &
                       oshad_day,                                    &
                       pdir_alb_town, psca_alb_town, ptran_win       )  
!     ###################################################
!
!!**** *AVERAGED_ALBEDO_TEB* computes averaged albedo for TEB scheme
!!
!!    PURPOSE
!!    -------
!!
!!    METHOD
!!    ------
!!   
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!!
!!    REFERENCE
!!    ---------
!!
!!    AUTHOR
!!    ------
!!
!!    V. Masson        Meteo-France
!!
!!    MODIFICATION
!!    ------------
!!
!!    Original    01/2004
!     C. de Munck & A. Lemonsu   09/2011 Greenroofs
!!    G. Pigeon                  09/2012 PTRAN_WIN as arguments
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
USE modd_surf_par, ONLY : xundef
USE modd_type_snow
!
USE modi_urban_solar_abs
USE mode_surf_snow_frac
!
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*    0.1    Declaration of arguments
!            ------------------------
!
 CHARACTER(LEN=3),   INTENT(IN) :: hbem          ! Building Energy model 'DEF' or 'BEM'
 CHARACTER(LEN=4),   INTENT(IN) :: hroad_dir     ! road direction options
              ! 'UNIF' : classical TEB version, all roads uniformely present
              ! 'ORIE' : specified road ORIEntation 
              !          (in fact many roads  direction
              !          possible because there can be many TEB patches)
 CHARACTER(LEN=4),   INTENT(IN) :: hwall_opt     ! wall options
              ! 'UNIF' : classical TEB version, all walls are uniform
              ! 'TWO ' : the two opposite walls receive different solar energy
REAL, DIMENSION(:), INTENT(IN) :: pzenith       ! zenithal solar angle
REAL, DIMENSION(:), INTENT(IN) :: pazim         ! solar azimuthal angle
!                                               ! (radian from N, clockwise)
!
REAL, DIMENSION(:), INTENT(IN) :: pbld          ! building fraction
REAL, DIMENSION(:), INTENT(IN) :: pgarden       ! green area fraction
REAL, DIMENSION(:), INTENT(IN) :: proad_dir     ! Road direction
                                                ! (deg from N, clockwise)
REAL, DIMENSION(:), INTENT(IN) :: proad         ! road fraction
REAL, DIMENSION(:), INTENT(IN) :: pfrac_gr      ! fraction of greenroof
REAL, DIMENSION(:), INTENT(IN) :: pfrac_panel   ! fraction of solar panels
REAL, DIMENSION(:), INTENT(IN) :: palb_panel    ! albedo   of solar panels
REAL, DIMENSION(:), INTENT(IN) :: pwall_o_hor   ! vertical surf. / horizontal surf.
REAL, DIMENSION(:), INTENT(IN) :: psvf_road     ! sky-view-factor from roads
REAL, DIMENSION(:), INTENT(IN) :: psvf_wall     ! sky-view-factor from walls
REAL, DIMENSION(:), INTENT(IN) :: psvf_garden   ! sky-view-factor from green areas
REAL, DIMENSION(:), INTENT(IN) :: pcan_hw_ratio ! canyon height/width ratio
!
REAL, DIMENSION(:), INTENT(IN) :: palb_roof     ! roof albedo
REAL, DIMENSION(:), INTENT(IN) :: palb_road     ! road albedo
REAL, DIMENSION(:), INTENT(IN) :: palb_wall     ! wall albedo
REAL, DIMENSION(:), INTENT(IN) :: palb_garden   ! green areas albedo
REAL, DIMENSION(:), INTENT(IN) :: palb_greenroof! green roof albedo
TYPE(surf_snow),    INTENT(IN) :: tsnow_roof    ! snow on roofs
TYPE(surf_snow),    INTENT(IN) :: tsnow_road    ! snow on roads
!
REAL, DIMENSION(:), INTENT(IN) :: pgr           ! Glazing ratio
REAL, DIMENSION(:), INTENT(IN) :: pshgc         ! Window solar transmittance
REAL, DIMENSION(:), INTENT(IN) :: pshgc_sh      ! Window + shading solar heat gain coef.
!
REAL, DIMENSION(:), INTENT(IN) :: pabs_win      ! Window solar absortance
REAL, DIMENSION(:), INTENT(OUT):: palb_win      ! Window solar reflectance
LOGICAL, DIMENSION(:),INTENT(IN) :: oshad_day   ! has shading been necessary this day ?
!
REAL, DIMENSION(:), INTENT(OUT):: pdir_alb_town ! direct albedo
REAL, DIMENSION(:), INTENT(OUT):: psca_alb_town ! diffuse albedo
REAL, DIMENSION(:), INTENT(IN) :: ptran_win     ! window transmittance
!
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL, DIMENSION(SIZE(PBLD)) :: zdn_roof       ! snow fraction 
REAL, DIMENSION(SIZE(PBLD)) :: zdn_road       ! on the surface
REAL, DIMENSION(SIZE(PBLD)) :: zdf_roof       ! free-snow fraction 
REAL, DIMENSION(SIZE(PBLD)) :: zdf_road       ! on the surface
LOGICAL, DIMENSION(SIZE(PBLD)) :: gmask       ! .false. (= no snow precip.)
!
!
REAL, DIMENSION(SIZE(PBLD)) :: zdir_sw        ! direct and diffuse shortwave radiation
REAL, DIMENSION(SIZE(PBLD)) :: zsca_sw        ! to mimic radiation behaviour of town
!
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_wall_a     ! shortwave absorbed by walls
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_wall_b     ! shortwave absorbed by walls
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_road       ! shortwave absorbed by roads
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_roof       ! shortwave absorbed by roofs
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_garden     ! shortwave absorbed by green areas
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_greenroof  ! shortwave absorbed by green roofs
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_snow_road  ! shortwave absorbed by snow
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_snow_roof  ! on roads, roofs,
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_panel      ! shortwave absorbed by solar panels
!
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_road       ! shortwave received by roads
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_wall_a     ! shortwave received by walls
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_wall_b     ! shortwave received by walls
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_garden     ! shortwave received by green areas
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_snow_road  ! shortwave received by snow on roads
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_roof       ! shortwave received by roofs
!
REAL, DIMENSION(SIZE(PBLD)) :: zsw_rad_garden ! total solar radiation reaching green areas
REAL, DIMENSION(SIZE(PBLD)) :: zabs_sw_win       ! shortwave absorbed by walls
REAL, DIMENSION(SIZE(PBLD)) :: zrec_sw_win       ! shortwave received by walls
REAL, DIMENSION(SIZE(PBLD)) :: zref_sw_grnd      !
REAL, DIMENSION(SIZE(PBLD)) :: zref_sw_fac       !
REAL, DIMENSION(SIZE(PBLD)) :: ztr_sw_win        !
REAL, DIMENSION(SIZE(PBLD)) :: ze_shading        !
LOGICAL, DIMENSION(SIZE(PBLD)) :: gshad_day
REAL(KIND=JPRB) :: zhook_handle
!-------------------------------------------------------------------------------
!
!* snow fractions
!  --------------
!
IF (lhook) CALL dr_hook('AVERAGED_ALBEDO_TEB',0,zhook_handle)
gmask(:) = .false.
 CALL snow_frac_road(tsnow_road%WSNOW(:,1,1),gmask,zdn_road,zdf_road)
 CALL snow_frac_roof(tsnow_roof%WSNOW(:,1,1),gmask,zdn_roof,zdf_roof)
!
!
!* town  direct and diffuse albedo
!  -------------------------------
!
zdir_sw=1.
zsca_sw=1.
gshad_day=.false.
IF (SIZE(oshad_day)>0) gshad_day=oshad_day
!
 CALL urban_solar_abs(hbem, hroad_dir, hwall_opt,                   &
                     zdir_sw, zsca_sw, pzenith, pazim,             &
                     pbld, pgarden, proad_dir, proad, pfrac_gr,    &
                     pwall_o_hor, pcan_hw_ratio,                   &
                     palb_roof,                                    &
                     palb_road, psvf_road, palb_wall, psvf_wall,   &
                     pfrac_panel, palb_panel,                      &
                     palb_garden, psvf_garden,                     &
                     palb_greenroof,                               &
                     tsnow_roof%ALB(:,1), tsnow_road%ALB(:,1),     &
                     zdn_roof, zdf_roof, zdn_road, zdf_road,       &
                     pgr, pabs_win, pshgc, pshgc_sh, palb_win,     &
                     zabs_sw_roof, zabs_sw_road,                   &
                     zabs_sw_wall_a, zabs_sw_wall_b,               &
                     zabs_sw_garden, zabs_sw_greenroof,            &
                     zabs_sw_snow_roof, zabs_sw_snow_road,         &
                     zabs_sw_panel,                                &
                     zrec_sw_road,  zrec_sw_snow_road,             &
                     zrec_sw_wall_a, zrec_sw_wall_b,               &
                     zrec_sw_garden, zrec_sw_roof,                 &
                     pdir_alb_town, psca_alb_town,                 &
                     zsw_rad_garden, zabs_sw_win, zrec_sw_win,     &
                     ptran_win,                                    &
                     zref_sw_grnd, zref_sw_fac, ztr_sw_win,        &
                     ze_shading, gshad_day, gmask                  )
IF (lhook) CALL dr_hook('AVERAGED_ALBEDO_TEB',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE averaged_albedo_teb