!$Author: joelevin $
!$Date: 2011-03-24 11:33:26 -0500 (Thu, 24 Mar 2011) $
!$Revision: 11724 $
!$HeadURL: https://svn.weru.ksu.edu/weru/weps1/trunk/weps.src/src/util/lib_solar/declination.for $

      real function   declination(idoy)

!     + + + PURPOSE + + +
!     This function calculates the declination of the earth with respect
!     the sun based on the day of the year

!     + + + KEYWORDS + + +
!     solar declination

!     + + + ARGUMENT DECLARATIONS + + +
      integer idoy

!     + + + ARGUMENT DEFINITIONS + + +
!     idoy   - Day of year

!     + + + LOCAL VARAIBLES + + +
      real b

!     + + + LOCAL DEFINITIONS + + +
!     b      - sub calculation (time of year, radians)

!     + + + COMMON BLOCKS + + +
      include 'p1unconv.inc'

!     + + + END SPECIFICATIONS + + +

!     Calculate declination angle (dec)
      b = (360.0/365.0)*(idoy-81.25) * degtorad				!h-55
      declination = 23.45*sin(b)							!h-58

      return
      end