Program: MKJD

  mkjd..................calculates the julian day from year, month, day
 
  call: jd=mkjd(y,m,d);
 
        y: the year
        m: the month
        d: the day (with decimal fractions, where .0 is 00:00h)
 
  result: jd: the julian day, as astronomers know it
              .0 is 12:00h (UT)
          THE RESULT IS VALID FOR NON-NEGATIVE JD ONLY !!!!
 
  Algorithm: from
 
  Jean Meeus
  "Astronomische Berechnungen"
  2., durchges. Auflage
  Joh. Ambrosius Barth, Leipzig, Berlin, Heidleberg, 1994
 
  The algorithm takes into account that the day after the
  04.Oct.1582 was the 15.Oct1582 due to the invention of
  the gregorian calendar!
 
  Note that the beginning of a Julian Day is 12:00 UT (because astronomers
  wanted to avoid a date change during night - their observation time)
  Also note that in astronomical considerations, a year 0 exists which
  is identical to the year 1 before christ in the gregorianic calendar.
 
  Martin Knapmeyer, 21.08.1995

Read M-File Source Code


eof.