equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/nut80.c

@@ -1,325 +0,0 @@
-#include "sofam.h"
-
-void iauNut80(double date1, double date2, double *dpsi, double *deps)
-/*
-**  - - - - - - - - -
-**   i a u N u t 8 0
-**  - - - - - - - - -
-**
-**  Nutation, IAU 1980 model.
-**
-**  Status:  canonical model.
-**
-**  Given:
-**     date1,date2   double    TT as a 2-part Julian Date (Note 1)
-**
-**  Returned:
-**     dpsi          double    nutation in longitude (radians)
-**     deps          double    nutation in obliquity (radians)
-**
-**  Notes:
-**
-**  1) The TT date date1+date2 is a Julian Date, apportioned in any
-**     convenient way between the two arguments.  For example,
-**     JD(TT)=2450123.7 could be expressed in any of these ways,
-**     among others:
-**
-**            date1          date2
-**
-**         2450123.7           0.0       (JD method)
-**         2451545.0       -1421.3       (J2000 method)
-**         2400000.5       50123.2       (MJD method)
-**         2450123.5           0.2       (date & time method)
-**
-**     The JD method is the most natural and convenient to use in
-**     cases where the loss of several decimal digits of resolution
-**     is acceptable.  The J2000 method is best matched to the way
-**     the argument is handled internally and will deliver the
-**     optimum resolution.  The MJD method and the date & time methods
-**     are both good compromises between resolution and convenience.
-**
-**  2) The nutation components are with respect to the ecliptic of
-**     date.
-**
-**  Called:
-**     iauAnpm      normalize angle into range +/- pi
-**
-**  Reference:
-**
-**     Explanatory Supplement to the Astronomical Almanac,
-**     P. Kenneth Seidelmann (ed), University Science Books (1992),
-**     Section 3.222 (p111).
-**
-**  This revision:  2008 September 30
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double t, el, elp, f, d, om, dp, de, arg, s, c;
-   int j;
-
-/* Units of 0.1 milliarcsecond to radians */
-   const double U2R = DAS2R / 1e4;
-
-/* ------------------------------------------------ */
-/* Table of multiples of arguments and coefficients */
-/* ------------------------------------------------ */
-
-/* The units for the sine and cosine coefficients are 0.1 mas and */
-/* the same per Julian century */
-
-   static const struct {
-      int nl,nlp,nf,nd,nom; /* coefficients of l,l',F,D,Om */
-      double sp,spt;        /* longitude sine, 1 and t coefficients */
-      double ce,cet;        /* obliquity cosine, 1 and t coefficients */
-   } x[] = {
-
-   /* 1-10 */
-      {  0,  0,  0,  0,  1, -171996.0, -174.2,  92025.0,    8.9 },
-      {  0,  0,  0,  0,  2,    2062.0,    0.2,   -895.0,    0.5 },
-      { -2,  0,  2,  0,  1,      46.0,    0.0,    -24.0,    0.0 },
-      {  2,  0, -2,  0,  0,      11.0,    0.0,      0.0,    0.0 },
-      { -2,  0,  2,  0,  2,      -3.0,    0.0,      1.0,    0.0 },
-      {  1, -1,  0, -1,  0,      -3.0,    0.0,      0.0,    0.0 },
-      {  0, -2,  2, -2,  1,      -2.0,    0.0,      1.0,    0.0 },
-      {  2,  0, -2,  0,  1,       1.0,    0.0,      0.0,    0.0 },
-      {  0,  0,  2, -2,  2,  -13187.0,   -1.6,   5736.0,   -3.1 },
-      {  0,  1,  0,  0,  0,    1426.0,   -3.4,     54.0,   -0.1 },
-
-   /* 11-20 */
-      {  0,  1,  2, -2,  2,    -517.0,    1.2,    224.0,   -0.6 },
-      {  0, -1,  2, -2,  2,     217.0,   -0.5,    -95.0,    0.3 },
-      {  0,  0,  2, -2,  1,     129.0,    0.1,    -70.0,    0.0 },
-      {  2,  0,  0, -2,  0,      48.0,    0.0,      1.0,    0.0 },
-      {  0,  0,  2, -2,  0,     -22.0,    0.0,      0.0,    0.0 },
-      {  0,  2,  0,  0,  0,      17.0,   -0.1,      0.0,    0.0 },
-      {  0,  1,  0,  0,  1,     -15.0,    0.0,      9.0,    0.0 },
-      {  0,  2,  2, -2,  2,     -16.0,    0.1,      7.0,    0.0 },
-      {  0, -1,  0,  0,  1,     -12.0,    0.0,      6.0,    0.0 },
-      { -2,  0,  0,  2,  1,      -6.0,    0.0,      3.0,    0.0 },
-
-   /* 21-30 */
-      {  0, -1,  2, -2,  1,      -5.0,    0.0,      3.0,    0.0 },
-      {  2,  0,  0, -2,  1,       4.0,    0.0,     -2.0,    0.0 },
-      {  0,  1,  2, -2,  1,       4.0,    0.0,     -2.0,    0.0 },
-      {  1,  0,  0, -1,  0,      -4.0,    0.0,      0.0,    0.0 },
-      {  2,  1,  0, -2,  0,       1.0,    0.0,      0.0,    0.0 },
-      {  0,  0, -2,  2,  1,       1.0,    0.0,      0.0,    0.0 },
-      {  0,  1, -2,  2,  0,      -1.0,    0.0,      0.0,    0.0 },
-      {  0,  1,  0,  0,  2,       1.0,    0.0,      0.0,    0.0 },
-      { -1,  0,  0,  1,  1,       1.0,    0.0,      0.0,    0.0 },
-      {  0,  1,  2, -2,  0,      -1.0,    0.0,      0.0,    0.0 },
-
-   /* 31-40 */
-      {  0,  0,  2,  0,  2,   -2274.0,   -0.2,    977.0,   -0.5 },
-      {  1,  0,  0,  0,  0,     712.0,    0.1,     -7.0,    0.0 },
-      {  0,  0,  2,  0,  1,    -386.0,   -0.4,    200.0,    0.0 },
-      {  1,  0,  2,  0,  2,    -301.0,    0.0,    129.0,   -0.1 },
-      {  1,  0,  0, -2,  0,    -158.0,    0.0,     -1.0,    0.0 },
-      { -1,  0,  2,  0,  2,     123.0,    0.0,    -53.0,    0.0 },
-      {  0,  0,  0,  2,  0,      63.0,    0.0,     -2.0,    0.0 },
-      {  1,  0,  0,  0,  1,      63.0,    0.1,    -33.0,    0.0 },
-      { -1,  0,  0,  0,  1,     -58.0,   -0.1,     32.0,    0.0 },
-      { -1,  0,  2,  2,  2,     -59.0,    0.0,     26.0,    0.0 },
-
-   /* 41-50 */
-      {  1,  0,  2,  0,  1,     -51.0,    0.0,     27.0,    0.0 },
-      {  0,  0,  2,  2,  2,     -38.0,    0.0,     16.0,    0.0 },
-      {  2,  0,  0,  0,  0,      29.0,    0.0,     -1.0,    0.0 },
-      {  1,  0,  2, -2,  2,      29.0,    0.0,    -12.0,    0.0 },
-      {  2,  0,  2,  0,  2,     -31.0,    0.0,     13.0,    0.0 },
-      {  0,  0,  2,  0,  0,      26.0,    0.0,     -1.0,    0.0 },
-      { -1,  0,  2,  0,  1,      21.0,    0.0,    -10.0,    0.0 },
-      { -1,  0,  0,  2,  1,      16.0,    0.0,     -8.0,    0.0 },
-      {  1,  0,  0, -2,  1,     -13.0,    0.0,      7.0,    0.0 },
-      { -1,  0,  2,  2,  1,     -10.0,    0.0,      5.0,    0.0 },
-
-   /* 51-60 */
-      {  1,  1,  0, -2,  0,      -7.0,    0.0,      0.0,    0.0 },
-      {  0,  1,  2,  0,  2,       7.0,    0.0,     -3.0,    0.0 },
-      {  0, -1,  2,  0,  2,      -7.0,    0.0,      3.0,    0.0 },
-      {  1,  0,  2,  2,  2,      -8.0,    0.0,      3.0,    0.0 },
-      {  1,  0,  0,  2,  0,       6.0,    0.0,      0.0,    0.0 },
-      {  2,  0,  2, -2,  2,       6.0,    0.0,     -3.0,    0.0 },
-      {  0,  0,  0,  2,  1,      -6.0,    0.0,      3.0,    0.0 },
-      {  0,  0,  2,  2,  1,      -7.0,    0.0,      3.0,    0.0 },
-      {  1,  0,  2, -2,  1,       6.0,    0.0,     -3.0,    0.0 },
-      {  0,  0,  0, -2,  1,      -5.0,    0.0,      3.0,    0.0 },
-
-   /* 61-70 */
-      {  1, -1,  0,  0,  0,       5.0,    0.0,      0.0,    0.0 },
-      {  2,  0,  2,  0,  1,      -5.0,    0.0,      3.0,    0.0 },
-      {  0,  1,  0, -2,  0,      -4.0,    0.0,      0.0,    0.0 },
-      {  1,  0, -2,  0,  0,       4.0,    0.0,      0.0,    0.0 },
-      {  0,  0,  0,  1,  0,      -4.0,    0.0,      0.0,    0.0 },
-      {  1,  1,  0,  0,  0,      -3.0,    0.0,      0.0,    0.0 },
-      {  1,  0,  2,  0,  0,       3.0,    0.0,      0.0,    0.0 },
-      {  1, -1,  2,  0,  2,      -3.0,    0.0,      1.0,    0.0 },
-      { -1, -1,  2,  2,  2,      -3.0,    0.0,      1.0,    0.0 },
-      { -2,  0,  0,  0,  1,      -2.0,    0.0,      1.0,    0.0 },
-
-   /* 71-80 */
-      {  3,  0,  2,  0,  2,      -3.0,    0.0,      1.0,    0.0 },
-      {  0, -1,  2,  2,  2,      -3.0,    0.0,      1.0,    0.0 },
-      {  1,  1,  2,  0,  2,       2.0,    0.0,     -1.0,    0.0 },
-      { -1,  0,  2, -2,  1,      -2.0,    0.0,      1.0,    0.0 },
-      {  2,  0,  0,  0,  1,       2.0,    0.0,     -1.0,    0.0 },
-      {  1,  0,  0,  0,  2,      -2.0,    0.0,      1.0,    0.0 },
-      {  3,  0,  0,  0,  0,       2.0,    0.0,      0.0,    0.0 },
-      {  0,  0,  2,  1,  2,       2.0,    0.0,     -1.0,    0.0 },
-      { -1,  0,  0,  0,  2,       1.0,    0.0,     -1.0,    0.0 },
-      {  1,  0,  0, -4,  0,      -1.0,    0.0,      0.0,    0.0 },
-
-   /* 81-90 */
-      { -2,  0,  2,  2,  2,       1.0,    0.0,     -1.0,    0.0 },
-      { -1,  0,  2,  4,  2,      -2.0,    0.0,      1.0,    0.0 },
-      {  2,  0,  0, -4,  0,      -1.0,    0.0,      0.0,    0.0 },
-      {  1,  1,  2, -2,  2,       1.0,    0.0,     -1.0,    0.0 },
-      {  1,  0,  2,  2,  1,      -1.0,    0.0,      1.0,    0.0 },
-      { -2,  0,  2,  4,  2,      -1.0,    0.0,      1.0,    0.0 },
-      { -1,  0,  4,  0,  2,       1.0,    0.0,      0.0,    0.0 },
-      {  1, -1,  0, -2,  0,       1.0,    0.0,      0.0,    0.0 },
-      {  2,  0,  2, -2,  1,       1.0,    0.0,     -1.0,    0.0 },
-      {  2,  0,  2,  2,  2,      -1.0,    0.0,      0.0,    0.0 },
-
-   /* 91-100 */
-      {  1,  0,  0,  2,  1,      -1.0,    0.0,      0.0,    0.0 },
-      {  0,  0,  4, -2,  2,       1.0,    0.0,      0.0,    0.0 },
-      {  3,  0,  2, -2,  2,       1.0,    0.0,      0.0,    0.0 },
-      {  1,  0,  2, -2,  0,      -1.0,    0.0,      0.0,    0.0 },
-      {  0,  1,  2,  0,  1,       1.0,    0.0,      0.0,    0.0 },
-      { -1, -1,  0,  2,  1,       1.0,    0.0,      0.0,    0.0 },
-      {  0,  0, -2,  0,  1,      -1.0,    0.0,      0.0,    0.0 },
-      {  0,  0,  2, -1,  2,      -1.0,    0.0,      0.0,    0.0 },
-      {  0,  1,  0,  2,  0,      -1.0,    0.0,      0.0,    0.0 },
-      {  1,  0, -2, -2,  0,      -1.0,    0.0,      0.0,    0.0 },
-
-   /* 101-106 */
-      {  0, -1,  2,  0,  1,      -1.0,    0.0,      0.0,    0.0 },
-      {  1,  1,  0, -2,  1,      -1.0,    0.0,      0.0,    0.0 },
-      {  1,  0, -2,  2,  0,      -1.0,    0.0,      0.0,    0.0 },
-      {  2,  0,  0,  2,  0,       1.0,    0.0,      0.0,    0.0 },
-      {  0,  0,  2,  4,  2,      -1.0,    0.0,      0.0,    0.0 },
-      {  0,  1,  0,  1,  0,       1.0,    0.0,      0.0,    0.0 }
-   };
-
-/* Number of terms in the series */
-   const int NT = (int) (sizeof x / sizeof x[0]);
-
-/*--------------------------------------------------------------------*/
-
-/* Interval between fundamental epoch J2000.0 and given date (JC). */
-   t = ((date1 - DJ00) + date2) / DJC;
-
-/* --------------------- */
-/* Fundamental arguments */
-/* --------------------- */
-
-/* Mean longitude of Moon minus mean longitude of Moon's perigee. */
-   el = iauAnpm(
-        (485866.733 + (715922.633 + (31.310 + 0.064 * t) * t) * t)
-        * DAS2R + fmod(1325.0 * t, 1.0) * D2PI);
-
-/* Mean longitude of Sun minus mean longitude of Sun's perigee. */
-   elp = iauAnpm(
-         (1287099.804 + (1292581.224 + (-0.577 - 0.012 * t) * t) * t)
-         * DAS2R + fmod(99.0 * t, 1.0) * D2PI);
-
-/* Mean longitude of Moon minus mean longitude of Moon's node. */
-   f = iauAnpm(
-       (335778.877 + (295263.137 + (-13.257 + 0.011 * t) * t) * t)
-       * DAS2R + fmod(1342.0 * t, 1.0) * D2PI);
-
-/* Mean elongation of Moon from Sun. */
-   d = iauAnpm(
-       (1072261.307 + (1105601.328 + (-6.891 + 0.019 * t) * t) * t)
-       * DAS2R + fmod(1236.0 * t, 1.0) * D2PI);
-
-/* Longitude of the mean ascending node of the lunar orbit on the */
-/* ecliptic, measured from the mean equinox of date. */
-   om = iauAnpm(
-        (450160.280 + (-482890.539 + (7.455 + 0.008 * t) * t) * t)
-        * DAS2R + fmod(-5.0 * t, 1.0) * D2PI);
-
-/* --------------- */
-/* Nutation series */
-/* --------------- */
-
-/* Initialize nutation components. */
-   dp = 0.0;
-   de = 0.0;
-
-/* Sum the nutation terms, ending with the biggest. */
-   for (j = NT-1; j >= 0; j--) {
-
-   /* Form argument for current term. */
-      arg = (double)x[j].nl  * el
-          + (double)x[j].nlp * elp
-          + (double)x[j].nf  * f
-          + (double)x[j].nd  * d
-          + (double)x[j].nom * om;
-
-   /* Accumulate current nutation term. */
-      s = x[j].sp + x[j].spt * t;
-      c = x[j].ce + x[j].cet * t;
-      if (s != 0.0) dp += s * sin(arg);
-      if (c != 0.0) de += c * cos(arg);
-   }
-
-/* Convert results from 0.1 mas units to radians. */
-   *dpsi = dp * U2R;
-   *deps = de * U2R;
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/nutm80.c

@@ -1,117 +0,0 @@
-#include "sofam.h"
-
-void iauNutm80(double date1, double date2, double rmatn[3][3])
-/*
-**  - - - - - - - - - -
-**   i a u N u t m 8 0
-**  - - - - - - - - - -
-**
-**  Form the matrix of nutation for a given date, IAU 1980 model.
-**
-**  Status:  support function.
-**
-**  Given:
-**     date1,date2    double          TDB date (Note 1)
-**
-**  Returned:
-**     rmatn          double[3][3]    nutation matrix
-**
-**  Notes:
-**
-**  1) The TT date date1+date2 is a Julian Date, apportioned in any
-**     convenient way between the two arguments.  For example,
-**     JD(TT)=2450123.7 could be expressed in any of these ways,
-**     among others:
-**
-**            date1          date2
-**
-**         2450123.7           0.0       (JD method)
-**         2451545.0       -1421.3       (J2000 method)
-**         2400000.5       50123.2       (MJD method)
-**         2450123.5           0.2       (date & time method)
-**
-**     The JD method is the most natural and convenient to use in
-**     cases where the loss of several decimal digits of resolution
-**     is acceptable.  The J2000 method is best matched to the way
-**     the argument is handled internally and will deliver the
-**     optimum resolution.  The MJD method and the date & time methods
-**     are both good compromises between resolution and convenience.
-**
-**  2) The matrix operates in the sense V(true) = rmatn * V(mean),
-**     where the p-vector V(true) is with respect to the true
-**     equatorial triad of date and the p-vector V(mean) is with
-**     respect to the mean equatorial triad of date.
-**
-**  Called:
-**     iauNut80     nutation, IAU 1980
-**     iauObl80     mean obliquity, IAU 1980
-**     iauNumat     form nutation matrix
-**
-**  This revision:  2008 May 12
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double dpsi, deps, epsa;
-
-
-/* Nutation components and mean obliquity. */
-   iauNut80(date1, date2, &dpsi, &deps);
-   epsa = iauObl80(date1, date2);
-
-/* Build the rotation matrix. */
-   iauNumat(epsa, dpsi, deps, rmatn);
-
-   return;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}

data/examples/celes_parts/obl06.c

@@ -1,118 +0,0 @@
-#include "sofam.h"
-
-double iauObl06(double date1, double date2)
-/*
-**  - - - - - - - - -
-**   i a u O b l 0 6
-**  - - - - - - - - -
-**
-**  Mean obliquity of the ecliptic, IAU 2006 precession model.
-**
-**  Status:  canonical model.
-**
-**  Given:
-**     date1,date2  double   TT as a 2-part Julian Date (Note 1)
-**
-**  Returned (function value):
-**                  double   obliquity of the ecliptic (radians, Note 2)
-**
-**  Notes:
-**
-**  1) The TT date date1+date2 is a Julian Date, apportioned in any
-**     convenient way between the two arguments.  For example,
-**     JD(TT)=2450123.7 could be expressed in any of these ways,
-**     among others:
-**
-**            date1          date2
-**
-**         2450123.7           0.0       (JD method)
-**         2451545.0       -1421.3       (J2000 method)
-**         2400000.5       50123.2       (MJD method)
-**         2450123.5           0.2       (date & time method)
-**
-**     The JD method is the most natural and convenient to use in
-**     cases where the loss of several decimal digits of resolution
-**     is acceptable.  The J2000 method is best matched to the way
-**     the argument is handled internally and will deliver the
-**     optimum resolution.  The MJD method and the date & time methods
-**     are both good compromises between resolution and convenience.
-**
-**  2) The result is the angle between the ecliptic and mean equator of
-**     date date1+date2.
-**
-**  Reference:
-**
-**     Hilton, J. et al., 2006, Celest.Mech.Dyn.Astron. 94, 351
-**
-**  This revision:  2009 March 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double t, eps0;
-
-
-/* Interval between fundamental date J2000.0 and given date (JC). */
-   t = ((date1 - DJ00) + date2) / DJC;
-
-/* Mean obliquity. */
-   eps0 = (84381.406     +
-          (-46.836769    +
-          ( -0.0001831   +
-          (  0.00200340  +
-          ( -0.000000576 +
-          ( -0.0000000434) * t) * t) * t) * t) * t) * DAS2R;
-
-   return eps0;
-
-/*----------------------------------------------------------------------
-**
-**  Celes is a wrapper of the SOFA Library for Ruby.
-**
-**  This file is redistributed and relicensed in accordance with 
-**  the SOFA Software License (http://www.iausofa.org/tandc.html).
-**
-**  The original library is available from IAU Standards of
-**  Fundamental Astronomy (http://www.iausofa.org/).
-**
-**
-**
-**
-**
-**  Copyright (C) 2013, Naoki Arita
-**  All rights reserved.
-**
-**  Redistribution and use in source and binary forms, with or without
-**  modification, are permitted provided that the following conditions
-**  are met:
-**
-**  1 Redistributions of source code must retain the above copyright
-**    notice, this list of conditions and the following disclaimer.
-**
-**  2 Redistributions in binary form must reproduce the above copyright
-**    notice, this list of conditions and the following disclaimer in
-**    the documentation and/or other materials provided with the
-**    distribution.
-**
-**  3 Neither the name of the Standards Of Fundamental Astronomy Board,
-**    the International Astronomical Union nor the names of its
-**    contributors may be used to endorse or promote products derived
-**    from this software without specific prior written permission.
-**
-**  THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS
-**  "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT
-**  LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
-**  FOR A PARTICULAR PURPOSE ARE DISCLAIMED.  IN NO EVENT SHALL THE
-**  COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
-**  INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING,
-**  BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-**  LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
-**  CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT
-**  LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN
-**  ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
-**  POSSIBILITY OF SUCH DAMAGE.
-**
-**--------------------------------------------------------------------*/
-}