equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/d2tf.c

@@ -1,160 +0,0 @@
-#include "sofam.h"
-
-void iauD2tf(int ndp, double days, char *sign, int ihmsf[4])
-/*
-**  - - - - - - - -
-**   i a u D 2 t f
-**  - - - - - - - -
-**
-**  Decompose days to hours, minutes, seconds, fraction.
-**
-**  Status:  vector/matrix support function.
-**
-**  Given:
-**     ndp     int     resolution (Note 1)
-**     days    double  interval in days
-**
-**  Returned:
-**     sign    char    '+' or '-'
-**     ihmsf   int[4]  hours, minutes, seconds, fraction
-**
-**  Notes:
-**
-**  1) The argument ndp is interpreted as follows:
-**
-**     ndp         resolution
-**      :      ...0000 00 00
-**     -7         1000 00 00
-**     -6          100 00 00
-**     -5           10 00 00
-**     -4            1 00 00
-**     -3            0 10 00
-**     -2            0 01 00
-**     -1            0 00 10
-**      0            0 00 01
-**      1            0 00 00.1
-**      2            0 00 00.01
-**      3            0 00 00.001
-**      :            0 00 00.000...
-**
-**  2) The largest positive useful value for ndp is determined by the
-**     size of days, the format of double on the target platform, and
-**     the risk of overflowing ihmsf[3].  On a typical platform, for
-**     days up to 1.0, the available floating-point precision might
-**     correspond to ndp=12.  However, the practical limit is typically
-**     ndp=9, set by the capacity of a 32-bit int, or ndp=4 if int is
-**     only 16 bits.
-**
-**  3) The absolute value of days may exceed 1.0.  In cases where it
-**     does not, it is up to the caller to test for and handle the
-**     case where days is very nearly 1.0 and rounds up to 24 hours,
-**     by testing for ihmsf[0]=24 and setting ihmsf[0-3] to zero.
-**
-**  This revision:  2008 May 17
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   int nrs, n;
-   double rs, rm, rh, a, w, ah, am, as, af;
-
-
-/* Handle sign. */
-   *sign = (char) ( ( days >= 0.0 ) ? '+' : '-' );
-
-/* Interval in seconds. */
-   a = DAYSEC * fabs(days);
-
-/* Pre-round if resolution coarser than 1s (then pretend ndp=1). */
-   if (ndp < 0) {
-      nrs = 1;
-      for (n = 1; n <= -ndp; n++) {
-          nrs *= (n == 2 || n == 4) ? 6 : 10;
-      }
-      rs = (double) nrs;
-      w = a / rs;
-      a = rs * dnint(w);
-   }
-
-/* Express the unit of each field in resolution units. */
-   nrs = 1;
-   for (n = 1; n <= ndp; n++) {
-      nrs *= 10;
-   }
-   rs = (double) nrs;
-   rm = rs * 60.0;
-   rh = rm * 60.0;
-
-/* Round the interval and express in resolution units. */
-   a = dnint(rs * a);
-
-/* Break into fields. */
-   ah = a / rh;
-   ah = dint(ah);
-   a -= ah * rh;
-   am = a / rm;
-   am = dint(am);
-   a -= am * rm;
-   as = a / rs;
-   as = dint(as);
-   af = a - as * rs;
-
-/* Return results. */
-   ihmsf[0] = (int) ah;
-   ihmsf[1] = (int) am;
-   ihmsf[2] = (int) as;
-   ihmsf[3] = (int) af;
-
-   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/dat.c

@@ -1,289 +0,0 @@
-#include "sofam.h"
-
-int iauDat(int iy, int im, int id, double fd, double *deltat )
-/*
-**  - - - - - - -
-**   i a u D a t
-**  - - - - - - -
-**
-**  For a given UTC date, calculate delta(AT) = TAI-UTC.
-**
-**     :------------------------------------------:
-**     :                                          :
-**     :                 IMPORTANT                :
-**     :                                          :
-**     :  A new version of this function must be  :
-**     :  produced whenever a new leap second is  :
-**     :  announced.  There are four items to     :
-**     :  change on each such occasion:           :
-**     :                                          :
-**     :  1) A new line must be added to the set  :
-**     :     of statements that initialize the    :
-**     :     array "changes".                     :
-**     :                                          :
-**     :  2) The parameter IYV must be set to     :
-**     :     the current year.                    :
-**     :                                          :
-**     :  3) The "Latest leap second" comment     :
-**     :     below must be set to the new leap    :
-**     :     second date.                         :
-**     :                                          :
-**     :  4) The "This revision" comment, later,  :
-**     :     must be set to the current date.     :
-**     :                                          :
-**     :  Change (2) must also be carried out     :
-**     :  whenever the function is re-issued,     :
-**     :  even if no leap seconds have been       :
-**     :  added.                                  :
-**     :                                          :
-**     :  Latest leap second:  2012 June 30       :
-**     :                                          :
-**     :__________________________________________:
-**
-**  Status:  support function.
-**
-**  Given:
-**     iy     int      UTC:  year (Notes 1 and 2)
-**     im     int            month (Note 2)
-**     id     int            day (Notes 2 and 3)
-**     fd     double         fraction of day (Note 4)
-**
-**  Returned:
-**     deltat double   TAI minus UTC, seconds
-**
-**  Returned (function value):
-**            int      status (Note 5):
-**                       1 = dubious year (Note 1)
-**                       0 = OK
-**                      -1 = bad year
-**                      -2 = bad month
-**                      -3 = bad day (Note 3)
-**                      -4 = bad fraction (Note 4)
-**
-**  Notes:
-**
-**  1) UTC began at 1960 January 1.0 (JD 2436934.5) and it is improper
-**     to call the function with an earlier date.  If this is attempted,
-**     zero is returned together with a warning status.
-**
-**     Because leap seconds cannot, in principle, be predicted in
-**     advance, a reliable check for dates beyond the valid range is
-**     impossible.  To guard against gross errors, a year five or more
-**     after the release year of the present function (see parameter
-**     IYV) is considered dubious.  In this case a warning status is
-**     returned but the result is computed in the normal way.
-**
-**     For both too-early and too-late years, the warning status is
-**     j=+1.  This is distinct from the error status j=-1, which
-**     signifies a year so early that JD could not be computed.
-**
-**  2) If the specified date is for a day which ends with a leap second,
-**     the UTC-TAI value returned is for the period leading up to the
-**     leap second.  If the date is for a day which begins as a leap
-**     second ends, the UTC-TAI returned is for the period following the
-**     leap second.
-**
-**  3) The day number must be in the normal calendar range, for example
-**     1 through 30 for April.  The "almanac" convention of allowing
-**     such dates as January 0 and December 32 is not supported in this
-**     function, in order to avoid confusion near leap seconds.
-**
-**  4) The fraction of day is used only for dates before the
-**     introduction of leap seconds, the first of which occurred at the
-**     end of 1971.  It is tested for validity (0 to 1 is the valid
-**     range) even if not used;  if invalid, zero is used and status
-**     j=-4 is returned.  For many applications, setting fd to zero is
-**     acceptable;  the resulting error is always less than 3 ms (and
-**     occurs only pre-1972).
-**
-**  5) The status value returned in the case where there are multiple
-**     errors refers to the first error detected.  For example, if the
-**     month and day are 13 and 32 respectively, j=-2 (bad month)
-**     will be returned.
-**
-**  6) In cases where a valid result is not available, zero is returned.
-**
-**  References:
-**
-**  1) For dates from 1961 January 1 onwards, the expressions from the
-**     file ftp://maia.usno.navy.mil/ser7/tai-utc.dat are used.
-**
-**  2) The 5ms timestep at 1961 January 1 is taken from 2.58.1 (p87) of
-**     the 1992 Explanatory Supplement.
-**
-**  Called:
-**     iauCal2jd    Gregorian calendar to Julian Day number
-**
-**  This revision:  2012 June 14
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-/* Release year for this version of iauDat */
-#define IYV (2012)
-
-/* Reference dates (MJD) and drift rates (s/day), pre leap seconds */
-   static const double drift[][2] = {
-      { 37300.0, 0.0012960 },
-      { 37300.0, 0.0012960 },
-      { 37300.0, 0.0012960 },
-      { 37665.0, 0.0011232 },
-      { 37665.0, 0.0011232 },
-      { 38761.0, 0.0012960 },
-      { 38761.0, 0.0012960 },
-      { 38761.0, 0.0012960 },
-      { 38761.0, 0.0012960 },
-      { 38761.0, 0.0012960 },
-      { 38761.0, 0.0012960 },
-      { 38761.0, 0.0012960 },
-      { 39126.0, 0.0025920 },
-      { 39126.0, 0.0025920 }
-   };
-
-/* Number of Delta(AT) expressions before leap seconds were introduced */
-#define NERA1 ((int) (sizeof drift / sizeof (double) / 2))
-
-/* Dates and Delta(AT)s */
-   static const struct {
-      int iyear, month;
-      double delat;
-   } changes[] = {
-      { 1960,  1,  1.4178180 },
-      { 1961,  1,  1.4228180 },
-      { 1961,  8,  1.3728180 },
-      { 1962,  1,  1.8458580 },
-      { 1963, 11,  1.9458580 },
-      { 1964,  1,  3.2401300 },
-      { 1964,  4,  3.3401300 },
-      { 1964,  9,  3.4401300 },
-      { 1965,  1,  3.5401300 },
-      { 1965,  3,  3.6401300 },
-      { 1965,  7,  3.7401300 },
-      { 1965,  9,  3.8401300 },
-      { 1966,  1,  4.3131700 },
-      { 1968,  2,  4.2131700 },
-      { 1972,  1, 10.0       },
-      { 1972,  7, 11.0       },
-      { 1973,  1, 12.0       },
-      { 1974,  1, 13.0       },
-      { 1975,  1, 14.0       },
-      { 1976,  1, 15.0       },
-      { 1977,  1, 16.0       },
-      { 1978,  1, 17.0       },
-      { 1979,  1, 18.0       },
-      { 1980,  1, 19.0       },
-      { 1981,  7, 20.0       },
-      { 1982,  7, 21.0       },
-      { 1983,  7, 22.0       },
-      { 1985,  7, 23.0       },
-      { 1988,  1, 24.0       },
-      { 1990,  1, 25.0       },
-      { 1991,  1, 26.0       },
-      { 1992,  7, 27.0       },
-      { 1993,  7, 28.0       },
-      { 1994,  7, 29.0       },
-      { 1996,  1, 30.0       },
-      { 1997,  7, 31.0       },
-      { 1999,  1, 32.0       },
-      { 2006,  1, 33.0       },
-      { 2009,  1, 34.0       },
-      { 2012,  7, 35.0       }
-   };
-
-/* Number of Delta(AT) changes */
-   const int NDAT = sizeof changes / sizeof changes[0];
-
-/* Miscellaneous local variables */
-   int j, i, m;
-   double da, djm0, djm;
-
-
-/* Initialize the result to zero. */
-   *deltat = da = 0.0;
-
-/* If invalid fraction of a day, set error status and give up. */
-   if (fd < 0.0 || fd > 1.0) return -4;
-
-/* Convert the date into an MJD. */
-   j = iauCal2jd(iy, im, id, &djm0, &djm);
-
-/* If invalid year, month, or day, give up. */
-   if (j < 0) return j;
-
-/* If pre-UTC year, set warning status and give up. */
-   if (iy < changes[0].iyear) return 1;
-
-/* If suspiciously late year, set warning status but proceed. */
-   if (iy > IYV + 5) j = 1;
-
-/* Combine year and month to form a date-ordered integer... */
-   m = 12*iy + im;
-
-/* ...and use it to find the preceding table entry. */
-   for (i = NDAT-1; i >=0; i--) {
-      if (m >= (12 * changes[i].iyear + changes[i].month)) break;
-   }
-
-/* Get the Delta(AT). */
-   da = changes[i].delat;
-
-/* If pre-1972, adjust for drift. */
-   if (i < NERA1) da += (djm + fd - drift[i][0]) * drift[i][1];
-
-/* Return the Delta(AT) value. */
-   *deltat = da;
-
-/* Return the status. */
-   return j;
-
-/*----------------------------------------------------------------------
-**
-**  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/depend

@@ -1,25 +0,0 @@
-celes_core.o: celes_core.c sofa.h sofam.h a2af.c a2tf.c af2a.c \
-		anp.c anpm.c bi00.c bp00.c bp06.c bpn2xy.c c2i00a.c c2i00b.c \
-		c2i06a.c c2ibpn.c c2ixy.c c2ixys.c c2s.c c2t00a.c c2t00b.c \
-		c2t06a.c c2tcio.c c2teqx.c c2tpe.c c2txy.c cal2jd.c cp.c cpv.c \
-		cr.c d2dtf.c d2tf.c dat.c dtdb.c dtf2d.c ee00.c ee00a.c ee00b.c \
-		ee06a.c eect00.c eform.c eo06a.c eors.c epb.c epb2jd.c epj.c \
-		epj2jd.c epv00.c eqeq94.c era00.c fad03.c fae03.c faf03.c \
-		faju03.c fal03.c falp03.c fama03.c fame03.c fane03.c faom03.c \
-		fapa03.c fasa03.c faur03.c fave03.c fk52h.c fk5hip.c fk5hz.c \
-		fw2m.c fw2xy.c gc2gd.c gc2gde.c gd2gc.c gd2gce.c gmst00.c \
-		gmst06.c gmst82.c gst00a.c gst00b.c gst06.c gst06a.c gst94.c \
-		h2fk5.c hfk5z.c ir.c jd2cal.c jdcalf.c num00a.c num00b.c \
-		num06a.c numat.c nut00a.c nut00b.c nut06a.c nut80.c nutm80.c \
-		obl06.c obl80.c p06e.c p2pv.c p2s.c pap.c pas.c pb06.c pdp.c \
-		pfw06.c plan94.c pm.c pmat00.c pmat06.c pmat76.c pmp.c pn.c \
-		pn00.c pn00a.c pn00b.c pn06.c pn06a.c pnm00a.c pnm00b.c \
-		pnm06a.c pnm80.c pom00.c ppp.c ppsp.c pr00.c prec76.c pv2p.c \
-		pv2s.c pvdpv.c pvm.c pvmpv.c pvppv.c pvstar.c pvu.c pvup.c \
-		pvxpv.c pxp.c rm2v.c rv2m.c rx.c rxp.c rxpv.c rxr.c ry.c \
-		rz.c s00.c s00a.c s00b.c s06.c s06a.c s2c.c s2p.c s2pv.c \
-		s2xpv.c sepp.c seps.c sp00.c starpm.c starpv.c sxp.c sxpv.c \
-		taitt.c taiut1.c taiutc.c tcbtdb.c tcgtt.c tdbtcb.c tdbtt.c \
-		tf2a.c tf2d.c tr.c trxp.c trxpv.c tttai.c tttcg.c tttdb.c \
-		ttut1.c ut1tai.c ut1tt.c ut1utc.c utctai.c utcut1.c xy06.c \
-		xys00a.c xys00b.c xys06a.c zp.c zpv.c zr.c