SUBROUTINE JULIAN(JDAY) C C Routine to return the current julian day C IMPLICIT INTEGER (A-Z) DIMENSION MONTHS(0:12) DATA MONTHS/0,31,59,90,120,151,181,212,243,273,304,334,365/ C CALL IDATE(MON,DAY,YEAR) ! Get current date JDAY = MONTHS(MON-1)+DAY IF (MOD(YEAR,4).EQ.0.AND.MON.GT.2) JDAY = JDAY+1 ! Leap year RETURN ! That's all END