c$Author: fredfox $
c$Date: 2002-03-13 21:15:49 $
c$Revision: 1.1 $
c$Source: /weru/cvs/weps/weps.src/util/solar/declination.for,v $

      real function   declination(idoy)

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

c     + + + KEYWORDS + + +
c     solar declination

c     + + + ARGUMENT DECLARATIONS + + +
      integer idoy

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

c     + + + LOCAL VARAIBLES + + +
      real b, dec

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

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

c     + + + END SPECIFICATIONS + + +

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

      return
      end