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(TOP, T, TPN, B,  PZENITH, PAZIM, &
                                     PALB_GARDEN, PALB_GREENROOF, PDIR_ALB_TOWN, PSCA_ALB_TOWN)
!     ###################################################
!
!!**** *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 B%XTRAN_WIN as arguments
!----------------------------------------------------------------------------
!
!*    0.     DECLARATION
!            -----------
!
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_diag_misc_teb_n, ONLY : diag_misc_teb_t
!
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
!            ------------------------
!
TYPE(teb_options_t), INTENT(INOUT) :: TOP
TYPE(teb_t), INTENT(INOUT) :: T
TYPE(bem_t), INTENT(INOUT) :: B
TYPE(teb_panel_t), INTENT(INOUT) :: TPN
!
REAL, DIMENSION(:), INTENT(IN) :: PZENITH       ! zenithal solar angle
REAL, DIMENSION(:), INTENT(IN) :: PAZIM         ! solar azimuthal angle
!                                               ! (radian from N, clockwise)
!
REAL, DIMENSION(:), INTENT(IN) :: PALB_GARDEN   ! green areas albedo
REAL, DIMENSION(:), INTENT(IN) :: PALB_GREENROOF! green roof albedo
REAL, DIMENSION(:), INTENT(OUT):: PDIR_ALB_TOWN ! direct albedo
REAL, DIMENSION(:), INTENT(OUT):: PSCA_ALB_TOWN ! diffuse albedo
!
!*    0.2    Declaration of local variables
!            ------------------------------
!
REAL, DIMENSION(SIZE(T%XBLD)) :: ZDN_ROOF       ! snow fraction 
REAL, DIMENSION(SIZE(T%XBLD)) :: ZDN_ROAD       ! on the surface
REAL, DIMENSION(SIZE(T%XBLD)) :: ZDF_ROOF       ! free-snow fraction 
REAL, DIMENSION(SIZE(T%XBLD)) :: ZDF_ROAD       ! on the surface
LOGICAL, DIMENSION(SIZE(T%XBLD)) :: GMASK       ! .false. (= no snow precip.)
!
!
REAL, DIMENSION(SIZE(T%XBLD)) :: ZDIR_SW        ! direct and diffuse shortwave radiation
REAL, DIMENSION(SIZE(T%XBLD)) :: ZSCA_SW        ! to mimic radiation behaviour of town
!
REAL, DIMENSION(SIZE(T%XBLD)) :: ZABS_SW_PANEL      ! shortwave absorbed by solar panels
!
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_ROAD       ! shortwave received by roads
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_WALL_A     ! shortwave received by walls
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_WALL_B     ! shortwave received by walls
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_GARDEN     ! shortwave received by green areas
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_SNOW_ROAD  ! shortwave received by snow on roads
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_ROOF       ! shortwave received by roofs
!
REAL, DIMENSION(SIZE(T%XBLD)) :: ZSW_RAD_GARDEN ! total solar radiation reaching green areas
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREC_SW_WIN       ! shortwave received by walls
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREF_SW_GRND      !
REAL, DIMENSION(SIZE(T%XBLD)) :: ZREF_SW_FAC       !
REAL, DIMENSION(SIZE(T%XBLD)) :: ZE_SHADING        !
LOGICAL, DIMENSION(SIZE(T%XBLD)) :: GSHAD_DAY
!
TYPE(diag_misc_teb_t) :: YDMT
!
REAL(KIND=JPRB) :: ZHOOK_HANDLE
!-------------------------------------------------------------------------------
!
!* snow fractions
!  --------------
!
IF (lhook) CALL dr_hook('AVERAGED_ALBEDO_TEB',0,zhook_handle)
gmask(:) = .false.
 CALL snow_frac_road(t%TSNOW_ROAD%WSNOW(:,1),gmask,zdn_road,zdf_road)
 CALL snow_frac_roof(t%TSNOW_ROOF%WSNOW(:,1),gmask,zdn_roof,zdf_roof)
!
!
!* town  direct and diffuse albedo
!  -------------------------------
!
zdir_sw=1.
zsca_sw=1.
gshad_day=.false.
IF (SIZE(b%LSHAD_DAY)>0) gshad_day=b%LSHAD_DAY
!
 CALL urban_solar_abs(top, t, b, ydmt, zdir_sw, zsca_sw, pzenith, pazim,  &
                      tpn%XFRAC_PANEL, tpn%XALB_PANEL, palb_garden,        &
                      t%XSVF_GARDEN, palb_greenroof, zdn_roof, zdf_roof,   &
                      zdn_road, zdf_road,  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, zrec_sw_win, zref_sw_grnd,           &
                      zref_sw_fac, ze_shading, gshad_day, gmask,           &
                      oalb_only=.true. )
!
IF (lhook) CALL dr_hook('AVERAGED_ALBEDO_TEB',1,zhook_handle)
!
!-------------------------------------------------------------------------------
!
END SUBROUTINE averaged_albedo_teb