/*
   julian - Numerator program to convert a month-day-year to a Julian date.
   
   Reference:  More Subroutine Sandwich, Grillo & Robertson, 1983, John
               Wiley and Sons, Inc., ppg 36

*/

"Enter year (e.g. 1989) "
get i
"Enter number of month  "
get j
"Enter day of month     "
get kl

if (j > 2)
   ma = j - 3
   ya = i
endif
if (j <= 2)
   ma = j + 9
   ya = i - 1
endif
ca = int (ya / 100)
da = ya - 100 * ca
n = int (146097 * ca / 4) + kl + int (1461 * da / 4) + 1721119 + int ((153 * ma + 2) / 5)
"Julian date is         "
n

