sw_daycycle.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 sw_daycycle(KI, PZENITH, PTOT_SW)
!     #############################################################
!
USE yomhook   ,ONLY : lhook,   dr_hook
USE parkind1  ,ONLY : jprb
!
!http://www.apesimulator.it/help/models/solarradiation/Calculating_extra-terrestrial_solar_radiation.html
!
INTEGER,             INTENT(IN)  :: ki        ! number of points
REAL, DIMENSION(KI), INTENT(IN)  :: pzenith     ! Zenith angle (rad)
REAL, DIMENSION(KI), INTENT(OUT) :: ptot_sw     ! Extraterrestrial solar radiation (W m-2)
!
!RJ: missing declaration
INTEGER :: jj
!
REAL      :: zc_sol  ! Solar constant [W/m2]
REAL      :: zd_corr ! Correction to acutal solar distance at any specific day of the year
INTEGER   :: zjulian ! Julian day of the year
REAL(KIND=JPRB) :: zhook_handle
!
IF (lhook) CALL dr_hook('SW_DAYCYCLE',0,zhook_handle)
!
zjulian = 197
zc_sol  = 4.921*1e6/3600                ! W m-2
zd_corr = 1 + 0.0334 * cos(0.01721 * zjulian - 0.0552)
!
DO jj=1,ki
  ptot_sw(jj) = zc_sol * zd_corr * cos(pzenith(jj))
END DO
!
IF (lhook) CALL dr_hook('SW_DAYCYCLE',1,zhook_handle)
!
END SUBROUTINE sw_daycycle