%N todayjul %F $MIR/src/subs/julday.for
%D Format the current day into a Julian day.
%P jm
%: julian-day, date, utilities
%B
      subroutine todayjul(julian)

      implicit none
      double precision julian

  Convert the current date to Julian date.  This is not to high
  accuracy, but it achieves the accuracy required. See "Astronomical
  Formulae for Calculators", Jean Meeus (Wiillmann-Bell Inc).
  The day is assumed to begin at 0 hours UT.

  Input:
    none

  Output:
    julian      Julian date.
