%N mobliqu %F $MIR/src/subs/ephem.for
%D Mean obliquity of the ecliptic
%P rjs
%: utilities
%B
        double precision function mobliq(jday)

        implicit none
        double precision jday

  Return the mean obliquity of the ecliptic.

  Input:
    jday       Julian day.
  Output:
    mobliq     Mean obliquity of the ecliptic, in radians.

  Reference:
    Explanatory Supplement ... page 114.
