!$Author: fredfox $
!$Date: 2002-09-09 16:09:16 $
!$Revision: 1.3 $
!$Source: /weru/cvs/weps/weps.src/util/solar/declination.for,v $

      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