trad_body.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.
MODULE modi_trad_body
INTERFACE
FUNCTION trad_body(PSCA_SW, PREF_SW_FAC, PREF_SW_GRND,   &
                   pemit_lw_fac, pemit_lw_grnd, plw_rad, &
                   pbld, pbld_height, pwall_o_hor,       &
                   pdir_sw, pzenith) result(ptrad_body)
REAL, DIMENSION(:), INTENT(IN)  :: psca_sw       ! Diffuse solar radiation (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pref_sw_fac   ! Solar radiation reflected by facade [wall + glazing] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pref_sw_grnd  ! Solar radiation reflected by ground [road + garden] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pemit_lw_fac  ! Longwave radiation emitted by the facade [wall + glazing] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pemit_lw_grnd ! Longwave radiation emitted by the ground [road + garden] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: plw_rad       ! Atmospheric longwave radiation (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pbld          ! plan area density of building (m2(bld)/m2(urban))
REAL, DIMENSION(:), INTENT(IN)  :: pbld_height   ! building height (m)
REAL, DIMENSION(:), INTENT(IN)  :: pwall_o_hor   ! ratio between facade and urban horizontal surface (m2(facade)/m2(urban))
REAL, DIMENSION(:), INTENT(IN), OPTIONAL :: pdir_sw !Direct solar radiation (W/m2)
REAL, DIMENSION(:), INTENT(IN), OPTIONAL :: pzenith !solar zenithal angle (rad from vert.)
REAL, DIMENSION(SIZE(PSCA_SW)) :: ptrad_body
END FUNCTION trad_body
END INTERFACE
END MODULE modi_trad_body
!   ##########################################################################
FUNCTION trad_body(PSCA_SW, PREF_SW_FAC, PREF_SW_GRND, PEMIT_LW_FAC, PEMIT_LW_GRND, PLW_RAD,&
                   pbld, pbld_height, pwall_o_hor, pdir_sw, pzenith) result(ptrad_body)
!   ##########################################################################
!
!!****  *TRAD_BODY  
!!
!!    PURPOSE
!!    -------
!
!     Computes the radiant temperature equivalent to the total radiation
!     received by the human body
!     
!!**  METHOD
!     ------
!
!!    EXTERNAL
!!    --------
!!
!!    IMPLICIT ARGUMENTS
!!    ------------------
!! a supplement
!!    MODD_CST
!!
!!    REFERENCE
!!    ---------
!!   www.utci.org
!!      
!!    AUTHOR
!!    ------
!!
!!      G. Pigeon           * Meteo-France *
!!
!!    MODIFICATIONS
!!    -------------
!!      Original  03/2011
!!      V.MASSON   08/2014 : bug in road view factor in computation of Universal Thermal Climate Index (diagnostic only)
!-------------------------------------------------------------------------------
!
!*       0.     DECLARATIONS
!               ------------
!
USE modd_csts, ONLY : xstefan, xpi
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
IMPLICIT NONE
!
!*      0.1    declarations of arguments
REAL, DIMENSION(:), INTENT(IN)  :: psca_sw       ! Diffuse solar radiation (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pref_sw_fac   ! Solar radiation reflected by facade [wall + glazing] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pref_sw_grnd  ! Solar radiation reflected by ground [road + garden] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pemit_lw_fac  ! Longwave radiation emitted by the facade [wall + glazing] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pemit_lw_grnd ! Longwave radiation emitted by the ground [road + garden] (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: plw_rad       ! Atmospheric longwave radiation (W/m2)
REAL, DIMENSION(:), INTENT(IN)  :: pbld          ! plan area density of building (m2(bld)/m2(urban))
REAL, DIMENSION(:), INTENT(IN)  :: pbld_height   ! building height (m)
REAL, DIMENSION(:), INTENT(IN)  :: pwall_o_hor   ! ratio between facade and urban horizontal surface (m2(facade)/m2(urban))
REAL, DIMENSION(:), INTENT(IN), OPTIONAL :: pdir_sw !Direct solar radiation (W/m2)
REAL, DIMENSION(:), INTENT(IN), OPTIONAL :: pzenith !solar zenithal angle (rad from vert.)
!REAL, DIMENSION(:), INTENT(IN) :: PDIR_SW !Direct solar radiation (W/m2)
!REAL, DIMENSION(:), INTENT(IN) :: PZENITH !solar zenithal angle (rad from vert.)
REAL, DIMENSION(SIZE(PSCA_SW))    :: ptrad_body
!REAL, DIMENSION(:) :: PTRAD_BODY

!*      0.2    declarations of local variables
REAL :: zhb = 1.7 !average height of human person (m)
REAL :: zab = 0.7 !absorption coef of solar radiation by human body
REAL :: zeb = 0.97 !emissivity of human body
!REAL, DIMENSION(SIZE(PBLD)) :: TRAD_BODY
REAL, DIMENSION(SIZE(PBLD)) :: zwroad !width of the road (m)
REAL, DIMENSION(SIZE(PBLD)) :: zl1, zl2, zl4 !lengths for view factor calculation
REAL, DIMENSION(SIZE(PBLD)) :: zffac !facade view factor of human body
REAL, DIMENSION(SIZE(PBLD)) :: zfgrnd !ground view factor of human body
REAL, DIMENSION(SIZE(PBLD)) :: zfsky !sky view factor of human body
REAL, DIMENSION(SIZE(PBLD)) :: zdirswbody !solar radiation received by human body
REAL, DIMENSION(SIZE(PBLD)) :: zelev !solar elevation angle
REAL, DIMENSION(SIZE(PBLD)) :: zradbody !total radiation received by human body
INTEGER :: jj
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('TRAD_BODY',0,zhook_handle)
!
DO jj = 1, SIZE(pbld_height) 
  !
  !*  1 - calculation of view factors
  zwroad(jj) = pbld_height(jj) * 2. * (1. - pbld(jj)) / pwall_o_hor(jj)
  !
  zl1(jj) = sqrt(zhb**2                   + (zwroad(jj)/2.)**2)
  zl2(jj) = sqrt( pbld_height(jj)**2      + (zwroad(jj)/2.)**2)
  zl4(jj) = sqrt((pbld_height(jj)-zhb)**2 + (zwroad(jj)/2.)**2)
  !
  zffac(jj) = (zl1(jj) + zl2(jj) - zwroad(jj)/2. - zl4(jj)) / (2. * zhb)
  zfgrnd(jj) = 0.5*zwroad(jj)/zhb
  zfgrnd(jj) = 0.5 * (zfgrnd(jj) + 1. - sqrt(zfgrnd(jj)**2 + 1.))
  zfsky(jj) = 1. - zffac(jj) - zfgrnd(jj)
  !
  !*  2 - base calculation for both sun and shade
  zradbody(jj) = zab/zeb * &
              ( psca_sw(jj)*zfsky(jj) + pref_sw_fac(jj)*zffac(jj) + pref_sw_grnd(jj)*zfgrnd(jj) ) &
              + plw_rad(jj)*zfsky(jj) + pemit_lw_fac(jj)*zffac(jj) + pemit_lw_grnd(jj)*zfgrnd(jj)
  !
ENDDO
!

!*  3 - add direct contribution in case of sunny conditions 
IF (present(pdir_sw) .AND. present(pzenith)) THEN
  DO jj = 1, SIZE(pbld_height)
    zelev(jj) = xpi/2. - pzenith(jj)
    IF (zelev(jj) < 1e-6) zelev(jj) = 0.
    zdirswbody(jj) = pdir_sw(jj) * 0.308 * cos( zelev(jj)*(1-(zelev(jj)*180./xpi)**2)/48402. )
      ! the direct solar radiation is weighted by a projected area factor which can be expressed by this equation
      ! for a rotationally symmetric human being (Fanger, 1970)
    zradbody(jj) = zradbody(jj) + zab/zeb*zdirswbody(jj)
  ENDDO
ENDIF

!PTRAD_BODY(:) = (ZRADBODY(:)/XSTEFAN)**0.25
DO jj=1, SIZE(pbld_height)
   ptrad_body(jj) = (zradbody(jj)/xstefan)**0.25
ENDDO
!
IF (lhook) CALL dr_hook('TRAD_BODY',1,zhook_handle)
!
END FUNCTION trad_body