equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/eform.c

@@ -1,147 +0,0 @@
-#include "sofam.h"
-
-int iauEform ( int n, double *a, double *f )
-/*
-**  - - - - - - - - -
-**   i a u E f o r m
-**  - - - - - - - - -
-**
-**  Earth reference ellipsoids.
-**
-**  Status:  canonical.
-**
-**  Given:
-**     n    int         ellipsoid identifier (Note 1)
-**
-**  Returned:
-**     a    double      equatorial radius (meters, Note 2)
-**     f    double      flattening (Note 2)
-**
-**  Returned (function value):
-**          int         status:  0 = OK
-**                              -1 = illegal identifier (Note 3)
-**
-**  Notes:
-**
-**  1) The identifier n is a number that specifies the choice of
-**     reference ellipsoid.  The following are supported:
-**
-**        n    ellipsoid
-**
-**        1     WGS84
-**        2     GRS80
-**        3     WGS72
-**
-**     The n value has no significance outside the SOFA software.  For
-**     convenience, symbols WGS84 etc. are defined in sofam.h.
-**
-**  2) The ellipsoid parameters are returned in the form of equatorial
-**     radius in meters (a) and flattening (f).  The latter is a number
-**     around 0.00335, i.e. around 1/298.
-**
-**  3) For the case where an unsupported n value is supplied, zero a and
-**     f are returned, as well as error status.
-**
-**  References:
-**
-**     Department of Defense World Geodetic System 1984, National
-**     Imagery and Mapping Agency Technical Report 8350.2, Third
-**     Edition, p3-2.
-**
-**     Moritz, H., Bull. Geodesique 66-2, 187 (1992).
-**
-**     The Department of Defense World Geodetic System 1972, World
-**     Geodetic System Committee, May 1974.
-**
-**     Explanatory Supplement to the Astronomical Almanac,
-**     P. Kenneth Seidelmann (ed), University Science Books (1992),
-**     p220.
-**
-**  This revision:  2012 February 23
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-
-/* Look up a and f for the specified reference ellipsoid. */
-   switch ( n ) {
-
-   case WGS84:
-      *a = 6378137.0;
-      *f = 1.0 / 298.257223563;
-      break;
-
-   case GRS80:
-      *a = 6378137.0;
-      *f = 1.0 / 298.257222101;
-      break;
-
-   case WGS72:
-      *a = 6378135.0;
-      *f = 1.0 / 298.26;
-      break;
-
-   default:
-
-   /* Invalid identifier. */
-      *a = 0.0;
-      *f = 0.0;
-      return -1;
-
-   }
-
-/* OK status. */
-   return 0;
-
-/*----------------------------------------------------------------------
-**
-**  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/eo06a.c

@@ -1,131 +0,0 @@
-#include "sofam.h"
-
-double iauEo06a(double date1, double date2)
-/*
-**  - - - - - - - - -
-**   i a u E o 0 6 a
-**  - - - - - - - - -
-**
-**  Equation of the origins, IAU 2006 precession and IAU 2000A nutation.
-**
-**  Status:  support function.
-**
-**  Given:
-**     date1,date2  double    TT as a 2-part Julian Date (Note 1)
-**
-**  Returned (function value):
-**                  double    equation of the origins in 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 equation of the origins is the distance between the true
-**     equinox and the celestial intermediate origin and, equivalently,
-**     the difference between Earth rotation angle and Greenwich
-**     apparent sidereal time (ERA-GST).  It comprises the precession
-**     (since J2000.0) in right ascension plus the equation of the
-**     equinoxes (including the small correction terms).
-**
-**  Called:
-**     iauPnm06a    classical NPB matrix, IAU 2006/2000A
-**     iauBpn2xy    extract CIP X,Y coordinates from NPB matrix
-**     iauS06       the CIO locator s, given X,Y, IAU 2006
-**     iauEors      equation of the origins, Given NPB matrix and s
-**
-**  References:
-**
-**     Capitaine, N. & Wallace, P.T., 2006, Astron.Astrophys. 450, 855
-**
-**     Wallace, P.T. & Capitaine, N., 2006, Astron.Astrophys. 459, 981
-**
-**  This revision:  2008 May 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double r[3][3], x, y, s, eo;
-
-
-/* Classical nutation x precession x bias matrix. */
-   iauPnm06a(date1, date2, r);
-
-/* Extract CIP coordinates. */
-   iauBpn2xy(r, &x, &y);
-
-/* The CIO locator, s. */
-   s = iauS06(date1, date2, x, y);
-
-/* Solve for the EO. */
-   eo = iauEors(r, s);
-
-   return eo;
-
-/*----------------------------------------------------------------------
-**
-**  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/eors.c

@@ -1,108 +0,0 @@
-#include "sofam.h"
-
-double iauEors(double rnpb[3][3], double s)
-/*
-**  - - - - - - - -
-**   i a u E o r s
-**  - - - - - - - -
-**
-**  Equation of the origins, given the classical NPB matrix and the
-**  quantity s.
-**
-**  Status:  support function.
-**
-**  Given:
-**     rnpb  double[3][3]  classical nutation x precession x bias matrix
-**     s     double        the quantity s (the CIO locator)
-**
-**  Returned (function value):
-**           double        the equation of the origins in radians.
-**
-**  Notes:
-**
-**  1)  The equation of the origins is the distance between the true
-**      equinox and the celestial intermediate origin and, equivalently,
-**      the difference between Earth rotation angle and Greenwich
-**      apparent sidereal time (ERA-GST).  It comprises the precession
-**      (since J2000.0) in right ascension plus the equation of the
-**      equinoxes (including the small correction terms).
-**
-**  2)  The algorithm is from Wallace & Capitaine (2006).
-**
-** References:
-**
-**     Capitaine, N. & Wallace, P.T., 2006, Astron.Astrophys. 450, 855
-**
-**     Wallace, P. & Capitaine, N., 2006, Astron.Astrophys. 459, 981
-**
-**  This revision:  2008 May 26
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double x, ax, xs, ys, zs, p, q, eo;
-
-
-/* Evaluate Wallace & Capitaine (2006) expression (16). */
-   x = rnpb[2][0];
-   ax = x / (1.0 + rnpb[2][2]);
-   xs = 1.0 - ax * x;
-   ys = -ax * rnpb[2][1];
-   zs = -x;
-   p = rnpb[0][0] * xs + rnpb[0][1] * ys + rnpb[0][2] * zs;
-   q = rnpb[1][0] * xs + rnpb[1][1] * ys + rnpb[1][2] * zs;
-   eo = ((p != 0) || (q != 0)) ? s - atan2(q, p) : s;
-
-   return eo;
-
-/*----------------------------------------------------------------------
-**
-**  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/epb.c

@@ -1,96 +0,0 @@
-#include "sofam.h"
-
-double iauEpb(double dj1, double dj2)
-/*
-**  - - - - - - -
-**   i a u E p b
-**  - - - - - - -
-**
-**  Julian Date to Besselian Epoch.
-**
-**  Status:  support function.
-**
-**  Given:
-**     dj1,dj2    double     Julian Date (see note)
-**
-**  Returned (function value):
-**                double     Besselian Epoch.
-**
-**  Note:
-**
-**     The Julian Date is supplied in two pieces, in the usual SOFA
-**     manner, which is designed to preserve time resolution.  The
-**     Julian Date is available as a single number by adding dj1 and
-**     dj2.  The maximum resolution is achieved if dj1 is 2451545D0
-**     (J2000.0).
-**
-**  Reference:
-**
-**     Lieske,J.H., 1979. Astron.Astrophys.,73,282.
-**
-**  This revision:  2009 December 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-/* J2000.0 minus B1900.0 (2415019.81352) in Julian days */
-   static const double D1900 = 36524.68648;
-
-   double epb;
-
-
-   epb = 1900.0 + ((dj1 - DJ00) + (dj2 + D1900)) / DTY;
-
-   return epb;
-
-/*----------------------------------------------------------------------
-**
-**  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.
-**
-**--------------------------------------------------------------------*/
-}