From 71dbe4604f4569029b56242d02f656317513bdc8 Mon Sep 17 00:00:00 2001 From: Momme Date: Sat, 14 Mar 2015 09:37:22 +0000 Subject: [PATCH] Fixed lon,lat units in zenith angle computation. --- zenith_angle.F90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/zenith_angle.F90 b/zenith_angle.F90 index de9e0b0..d1a23af 100644 --- a/zenith_angle.F90 +++ b/zenith_angle.F90 @@ -81,10 +81,10 @@ subroutine do_surface(self,_ARGUMENTS_DO_SURFACE_) _GET_HORIZONTAL_(self%id_lat,lat) ! Sun hour angle : - zenithA = (h-12._rk)*15._rk*deg2rad + lon + zenithA = ((h-12._rk)*15._rk + lon)*deg2rad ! Cosine of the solar zenith angle : - zenithA =sin(lat)*sin(sundec)+cos(lat)*cos(sundec)*cos(zenithA) + zenithA =sin(lat*deg2rad)*sin(sundec)+cos(lat*deg2rad)*cos(sundec)*cos(zenithA) zenithA =acos(max(0._rk,zenithA))/deg2rad _SET_HORIZONTAL_DIAGNOSTIC_(self%id_zenith_angle,zenithA)