julian.f

Go to the documentation of this file.

      integer function julian(DAY, MONTH, YEAR)
C****************************************************************
C
C     Converts Gregorian DAY/MONTH/YEAR into Julian day.
C
C     Input parameters:
C     ----------------
C     DAY:    I   Day of MONTH in YEAR; 1 to 31
C     MONTH:  I   Month in YEAR; 1 to 12
C     YEAR:   I   Year in which to convert DAY/MONTH
C
C     Function:
C     --------
C     LEAP:   L   SP GEN_UTILS
C
C     Created by Michael Darzi, GSC, 9/89.
C               To UNIX: BD Schieber, SAIC/GSC, 1993
C****************************************************************
      implicit none
      logical leap
      integer day, month, year, jdays(12,2), ii
      data jdays/0,31,59,90,120,151,181,212,243,273,304,334,
     .           0,31,60,91,121,152,182,213,244,274,305,335/
 
      ii = 1
      if (leap(year)) ii = 2
      julian = min(max(1,day),31) + jdays(min(max(1,month),12),ii)
 
      return
      end