equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/ut1utc.c

@@ -1,193 +0,0 @@
-#include "sofam.h"
-
-int iauUt1utc(double ut11, double ut12, double dut1,
-              double *utc1, double *utc2)
-/*
-**  - - - - - - - - - -
-**   i a u U t 1 u t c
-**  - - - - - - - - - -
-**
-**  Time scale transformation:  Universal Time, UT1, to Coordinated
-**  Universal Time, UTC.
-**
-**  Status:  canonical.
-**
-**  Given:
-**     ut11,ut12  double   UT1 as a 2-part Julian Date (Note 1)
-**     dut1       double   Delta UT1: UT1-UTC in seconds (Note 2)
-**
-**  Returned:
-**     utc1,utc2  double   UTC as a 2-part quasi Julian Date (Notes 3,4)
-**
-**  Returned (function value):
-**                int      status: +1 = dubious year (Note 5)
-**                                  0 = OK
-**                                 -1 = unacceptable date
-**
-**  Notes:
-**
-**  1) ut11+ut12 is Julian Date, apportioned in any convenient way
-**     between the two arguments, for example where ut11 is the Julian
-**     Day Number and ut12 is the fraction of a day.  The returned utc1
-**     and utc2 form an analogous pair, except that a special convention
-**     is used, to deal with the problem of leap seconds - see Note 3.
-**
-**  2) Delta UT1 can be obtained from tabulations provided by the
-**     International Earth Rotation and Reference Systems Service.  The
-**     value changes abruptly by 1s at a leap second;  however, close to
-**     a leap second the algorithm used here is tolerant of the "wrong"
-**     choice of value being made.
-**
-**  3) JD cannot unambiguously represent UTC during a leap second unless
-**     special measures are taken.  The convention in the present
-**     function is that the returned quasi JD day UTC1+UTC2 represents
-**     UTC days whether the length is 86399, 86400 or 86401 SI seconds.
-**
-**  4) The function iauD2dtf can be used to transform the UTC quasi-JD
-**     into calendar date and clock time, including UTC leap second
-**     handling.
-**
-**  5) The warning status "dubious year" flags UTCs that predate the
-**     introduction of the time scale and that are too far in the future
-**     to be trusted.  See iauDat for further details.
-**
-**  Called:
-**     iauJd2cal    JD to Gregorian calendar
-**     iauDat       delta(AT) = TAI-UTC
-**     iauCal2jd    Gregorian calendar to JD
-**
-**  References:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**     Explanatory Supplement to the Astronomical Almanac,
-**     P. Kenneth Seidelmann (ed), University Science Books (1992)
-**
-**  This revision:  2011 May 14
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   int big1;
-   int i, iy, im, id, js;
-   double duts, u1, u2, d1, dats1, d2, fd, dats2, ddats, us1, us2, du;
-
-
-/* UT1-UTC in seconds. */
-   duts = dut1;
-
-/* Put the two parts of the UT1 into big-first order. */
-   big1 = ( ut11 >= ut12 );
-   if ( big1 ) {
-      u1 = ut11;
-      u2 = ut12;
-   } else {
-      u1 = ut12;
-      u2 = ut11;
-   }
-
-/* See if the UT1 can possibly be in a leap-second day. */
-   d1 = u1;
-   dats1 = 0;
-   for ( i = -1; i <= 3; i++ ) {
-      d2 = u2 + (double) i;
-      if ( iauJd2cal(d1, d2, &iy, &im, &id, &fd) ) return -1;
-      js = iauDat(iy, im, id, 0.0, &dats2);
-      if ( js < 0 ) return -1;
-      if ( i == - 1 ) dats1 = dats2;
-      ddats = dats2 - dats1;
-      if ( fabs(ddats) >= 0.5 ) {
-
-      /* Yes, leap second nearby: ensure UT1-UTC is "before" value. */
-         if ( ddats * duts >= 0 ) duts -= ddats;
-
-      /* UT1 for the start of the UTC day that ends in a leap. */
-         if ( iauCal2jd(iy, im, id, &d1, &d2) ) return -1;
-         us1 = d1;
-         us2 = d2 - 1.0 + duts/DAYSEC;
-
-      /* Is the UT1 after this point? */
-         du = u1 - us1;
-         du += u2 - us2;
-         if ( du > 0 ) {
-
-         /* Yes:  fraction of the current UTC day that has elapsed. */
-            fd = du * DAYSEC / ( DAYSEC + ddats );
-
-         /* Ramp UT1-UTC to bring about SOFA's JD(UTC) convention. */
-            duts += ddats * ( fd <= 1.0 ? fd : 1.0 );
-         }
-
-      /* Done. */
-         break;
-      }
-      dats1 = dats2;
-   }
-
-/* Subtract the (possibly adjusted) UT1-UTC from UT1 to give UTC. */
-   u2 -= duts / DAYSEC;
-
-/* Result, safeguarding precision. */
-   if ( big1 ) {
-      *utc1 = u1;
-      *utc2 = u2;
-   } else {
-      *utc1 = u2;
-      *utc2 = u1;
-   }
-
-/* Status. */
-   return js;
-
-/*----------------------------------------------------------------------
-**
-**  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/utctai.c

@@ -1,163 +0,0 @@
-#include "sofam.h"
-
-int iauUtctai(double utc1, double utc2, double *tai1, double *tai2)
-/*
-**  - - - - - - - - - -
-**   i a u U t c t a i
-**  - - - - - - - - - -
-**
-**  Time scale transformation:  Coordinated Universal Time, UTC, to
-**  International Atomic Time, TAI.
-**
-**  Status:  canonical.
-**
-**  Given:
-**     utc1,utc2  double   UTC as a 2-part quasi Julian Date (Notes 1-4)
-**
-**  Returned:
-**     tai1,tai2  double   TAI as a 2-part Julian Date (Note 5)
-**
-**  Returned (function value):
-**                int      status: +1 = dubious year (Note 3)
-**                                  0 = OK
-**                                 -1 = unacceptable date
-**
-**  Notes:
-**
-**  1) utc1+utc2 is quasi Julian Date (see Note 2), apportioned in any
-**     convenient way between the two arguments, for example where utc1
-**     is the Julian Day Number and utc2 is the fraction of a day.
-**
-**  2) JD cannot unambiguously represent UTC during a leap second unless
-**     special measures are taken.  The convention in the present
-**     function is that the JD day represents UTC days whether the
-**     length is 86399, 86400 or 86401 SI seconds.
-**
-**  3) The warning status "dubious year" flags UTCs that predate the
-**     introduction of the time scale and that are too far in the future
-**     to be trusted.  See iauDat  for further details.
-**
-**  4) The function iauDtf2d converts from calendar date and time of day
-**     into 2-part Julian Date, and in the case of UTC implements the
-**     leap-second-ambiguity convention described above.
-**
-**  5) The returned TAI1,TAI2 are such that their sum is the TAI Julian
-**     Date.
-**
-**  Called:
-**     iauJd2cal    JD to Gregorian calendar
-**     iauDat       delta(AT) = TAI-UTC
-**     iauCal2jd    Gregorian calendar to JD
-**
-**  References:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**     Explanatory Supplement to the Astronomical Almanac,
-**     P. Kenneth Seidelmann (ed), University Science Books (1992)
-**
-**  This revision:  2010 September 10
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-**
-*/
-{
-   int big1;
-   int iy, im, id, js, iyt, imt, idt;
-   double u1, u2, fd, dats, fdt, datst, ddat, z1, z2, a2;
-
-
-/* Put the two parts of the UTC into big-first order. */
-   big1 = ( utc1 >= utc2 );
-   if ( big1 ) {
-      u1 = utc1;
-      u2 = utc2;
-   } else {
-      u1 = utc2;
-      u2 = utc1;
-   }
-
-/* Get TAI-UTC now. */
-   if ( iauJd2cal(u1, u2, &iy, &im, &id, &fd) ) return -1;
-   js = iauDat(iy, im, id, fd, &dats);
-   if ( js < 0 ) return -1;
-
-/* Get TAI-UTC tomorrow. */
-   if ( iauJd2cal(u1+1.5, u2-fd, &iyt, &imt, &idt, &fdt) ) return -1;
-   js = iauDat(iyt, imt, idt, fdt, &datst);
-   if ( js < 0 ) return -1;
-
-/* If today ends in a leap second, scale the fraction into SI days. */
-   ddat = datst - dats;
-   if ( fabs(ddat) > 0.5 ) fd += fd * ddat / DAYSEC;
-
-/* Today's calendar date to 2-part JD. */
-   if ( iauCal2jd(iy, im, id, &z1, &z2) ) return -1;
-
-/* Assemble the TAI result, preserving the UTC split and order. */
-   a2 = z1 - u1;
-   a2 += z2;
-   a2 += fd + dats / DAYSEC;
-   if ( big1 ) {
-      *tai1 = u1;
-      *tai2 = a2;
-   } else {
-      *tai1 = a2;
-      *tai2 = u1;
-   }
-
-/* Status. */
-   return js;
-
-/*----------------------------------------------------------------------
-**
-**  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/utcut1.c

@@ -1,151 +0,0 @@
-#include "sofam.h"
-
-int iauUtcut1(double utc1, double utc2, double dut1,
-              double *ut11, double *ut12)
-/*
-**  - - - - - - - - - -
-**   i a u U t c u t 1
-**  - - - - - - - - - -
-**
-**  Time scale transformation:  Coordinated Universal Time, UTC, to
-**  Universal Time, UT1.
-**
-**  Status:  canonical.
-**
-**  Given:
-**     utc1,utc2  double   UTC as a 2-part quasi Julian Date (Notes 1-4)
-**     dut1       double   Delta UT1 = UT1-UTC in seconds (Note 5)
-**
-**  Returned:
-**     ut11,ut12  double   UT1 as a 2-part Julian Date (Note 6)
-**
-**  Returned (function value):
-**                int      status: +1 = dubious year (Note 7)
-**                                  0 = OK
-**                                 -1 = unacceptable date
-**
-**  Notes:
-**
-**  1) utc1+utc2 is quasi Julian Date (see Note 2), apportioned in any
-**     convenient way between the two arguments, for example where utc1
-**     is the Julian Day Number and utc2 is the fraction of a day.
-**
-**  2) JD cannot unambiguously represent UTC during a leap second unless
-**     special measures are taken.  The convention in the present
-**     function is that the JD day represents UTC days whether the
-**     length is 86399, 86400 or 86401 SI seconds.
-**
-**  3) The warning status "dubious year" flags UTCs that predate the
-**     introduction of the time scale and that are too far in the future
-**     to be trusted.  See iauDat  for further details.
-**
-**  4) The function iauDtf2d  converts from calendar date and time of
-**     day into 2-part Julian Date, and in the case of UTC implements
-**     the leap-second-ambiguity convention described above.
-**
-**  5) Delta UT1 can be obtained from tabulations provided by the
-**     International Earth Rotation and Reference Systems Service.  It
-**     It is the caller's responsibility to supply a DUT argument
-**     containing the UT1-UTC value that matches the given UTC.
-**
-**  6) The returned ut11,ut12 are such that their sum is the UT1 Julian
-**     Date.
-**
-**  7) The warning status "dubious year" flags UTCs that predate the
-**     introduction of the time scale and that are too far in the future
-**     to be trusted.  See iauDat for further details.
-**
-**  References:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**     Explanatory Supplement to the Astronomical Almanac,
-**     P. Kenneth Seidelmann (ed), University Science Books (1992)
-**
-**  Called:
-**     iauJd2cal    JD to Gregorian calendar
-**     iauDat       delta(AT) = TAI-UTC
-**     iauUtctai    UTC to TAI
-**     iauTaiut1    TAI to UT1
-**
-**  This revision:  2010 May 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   int iy, im, id, js, jw;
-   double w, dat, dta, tai1, tai2;
-
-
-/* Look up TAI-UTC. */
-   if ( iauJd2cal(utc1, utc2, &iy, &im, &id, &w) ) return -1;
-   js = iauDat ( iy, im, id, 0.0, &dat);
-   if ( js < 0 ) return -1;
-
-/* Form UT1-TAI. */
-   dta = dut1 - dat;
-
-/* UTC to TAI to UT1. */
-   jw = iauUtctai(utc1, utc2, &tai1, &tai2);
-   if ( jw < 0 ) {
-      return -1;
-   } else if ( jw > 0 ) {
-      js = jw;
-   }
-   if ( iauTaiut1(tai1, tai2, dta, ut11, ut12) ) return -1;
-
-/* Status. */
-   return js;
-
-/*----------------------------------------------------------------------
-**
-**  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.
-**
-**--------------------------------------------------------------------*/
-}