equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/sofam.h

@@ -1,155 +0,0 @@
-#ifndef SOFAMHDEF
-#define SOFAMHDEF
-
-/*
-**  - - - - - - - -
-**   s o f a m . h
-**  - - - - - - - -
-**
-**  Macros used by SOFA library.
-**
-**  Please note that the constants defined below are to be used only in
-**  the context of the SOFA software, and have no other official IAU
-**  status.
-**
-**  This revision:   2012 February 23
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-
-#include "sofa.h"
-
-
-/* Pi */
-#define DPI (3.141592653589793238462643)
-
-/* 2Pi */
-#define D2PI (6.283185307179586476925287)
-
-/* Degrees to radians */
-#define DD2R (1.745329251994329576923691e-2)
-
-/* Radians to arcseconds */
-#define DR2AS (206264.8062470963551564734)
-
-/* Arcseconds to radians */
-#define DAS2R (4.848136811095359935899141e-6)
-
-/* Seconds of time to radians */
-#define DS2R (7.272205216643039903848712e-5)
-
-/* Arcseconds in a full circle */
-#define TURNAS (1296000.0)
-
-/* Milliarcseconds to radians */
-#define DMAS2R (DAS2R / 1e3)
-
-/* Length of tropical year B1900 (days) */
-#define DTY (365.242198781)
-
-/* Seconds per day. */
-#define DAYSEC (86400.0)
-
-/* Days per Julian year */
-#define DJY (365.25)
-
-/* Days per Julian century */
-#define DJC (36525.0)
-
-/* Days per Julian millennium */
-#define DJM (365250.0)
-
-/* Reference epoch (J2000.0), Julian Date */
-#define DJ00 (2451545.0)
-
-/* Julian Date of Modified Julian Date zero */
-#define DJM0 (2400000.5)
-
-/* Reference epoch (J2000.0), Modified Julian Date */
-#define DJM00 (51544.5)
-
-/* 1977 Jan 1.0 as MJD */
-#define DJM77 (43144.0)
-
-/* TT minus TAI (s) */
-#define TTMTAI (32.184)
-
-/* AU (m) */
-#define DAU (149597870e3)
-
-/* Speed of light (AU per day) */
-#define DC (DAYSEC / 499.004782)
-
-/* L_G = 1 - d(TT)/d(TCG) */
-#define ELG (6.969290134e-10)
-
-/* L_B = 1 - d(TDB)/d(TCB), and TDB (s) at TAI 1977/1/1.0 */
-#define ELB (1.550519768e-8)
-#define TDB0 (-6.55e-5)
-
-/* dint(A) - truncate to nearest whole number towards zero (double) */
-#define dint(A) ((A)<0.0?ceil(A):floor(A))
-
-/* dnint(A) - round to nearest whole number (double) */
-#define dnint(A) ((A)<0.0?ceil((A)-0.5):floor((A)+0.5))
-
-/* dsign(A,B) - magnitude of A with sign of B (double) */
-#define dsign(A,B) ((B)<0.0?-fabs(A):fabs(A))
-
-/* Reference ellipsoids */
-#define WGS84 1
-#define GRS80 2
-#define WGS72 3
-
-#endif
-
-/*----------------------------------------------------------------------
-**
-**  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/sp00.c

@@ -1,118 +0,0 @@
-#include "sofam.h"
-
-double iauSp00(double date1, double date2)
-/*
-**  - - - - - - - -
-**   i a u S p 0 0
-**  - - - - - - - -
-**
-**  The TIO locator s', positioning the Terrestrial Intermediate Origin
-**  on the equator of the Celestial Intermediate Pole.
-**
-**  Status:  canonical model.
-**
-**  Given:
-**     date1,date2  double    TT as a 2-part Julian Date (Note 1)
-**
-**  Returned (function value):
-**                  double    the TIO locator s' in 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 TIO locator s' is obtained from polar motion observations by
-**     numerical integration, and so is in essence unpredictable.
-**     However, it is dominated by a secular drift of about
-**     47 microarcseconds per century, which is the approximation
-**     evaluated by the present function.
-**
-**  Reference:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**  This revision:  2008 May 24
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double t, sp;
-
-
-/* Interval between fundamental epoch J2000.0 and current date (JC). */
-   t = ((date1 - DJ00) + date2) / DJC;
-
-/* Approximate s'. */
-   sp = -47e-6 * t * DAS2R;
-
-   return sp;
-
-/*----------------------------------------------------------------------
-**
-**  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/starpm.c

@@ -1,205 +0,0 @@
-#include "sofam.h"
-
-int iauStarpm(double ra1, double dec1,
-              double pmr1, double pmd1, double px1, double rv1,
-              double ep1a, double ep1b, double ep2a, double ep2b,
-              double *ra2, double *dec2,
-              double *pmr2, double *pmd2, double *px2, double *rv2)
-/*
-**  - - - - - - - - - -
-**   i a u S t a r p m
-**  - - - - - - - - - -
-**
-**  Star proper motion:  update star catalog data for space motion.
-**
-**  Status:  support function.
-**
-**  Given:
-**     ra1    double     right ascension (radians), before
-**     dec1   double     declination (radians), before
-**     pmr1   double     RA proper motion (radians/year), before
-**     pmd1   double     Dec proper motion (radians/year), before
-**     px1    double     parallax (arcseconds), before
-**     rv1    double     radial velocity (km/s, +ve = receding), before
-**     ep1a   double     "before" epoch, part A (Note 1)
-**     ep1b   double     "before" epoch, part B (Note 1)
-**     ep2a   double     "after" epoch, part A (Note 1)
-**     ep2b   double     "after" epoch, part B (Note 1)
-**
-**  Returned:
-**     ra2    double     right ascension (radians), after
-**     dec2   double     declination (radians), after
-**     pmr2   double     RA proper motion (radians/year), after
-**     pmd2   double     Dec proper motion (radians/year), after
-**     px2    double     parallax (arcseconds), after
-**     rv2    double     radial velocity (km/s, +ve = receding), after
-**
-**  Returned (function value):
-**            int        status:
-**                          -1 = system error (should not occur)
-**                           0 = no warnings or errors
-**                           1 = distance overridden (Note 6)
-**                           2 = excessive velocity (Note 7)
-**                           4 = solution didn't converge (Note 8)
-**                        else = binary logical OR of the above warnings
-**
-**  Notes:
-**
-**  1) The starting and ending TDB dates ep1a+ep1b and ep2a+ep2b are
-**     Julian Dates, apportioned in any convenient way between the two
-**     parts (A and B).  For example, JD(TDB)=2450123.7 could be
-**     expressed in any of these ways, among others:
-**
-**             epna          epnb
-**
-**         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) In accordance with normal star-catalog conventions, the object's
-**     right ascension and declination are freed from the effects of
-**     secular aberration.  The frame, which is aligned to the catalog
-**     equator and equinox, is Lorentzian and centered on the SSB.
-**
-**     The proper motions are the rate of change of the right ascension
-**     and declination at the catalog epoch and are in radians per TDB
-**     Julian year.
-**
-**     The parallax and radial velocity are in the same frame.
-**
-**  3) Care is needed with units.  The star coordinates are in radians
-**     and the proper motions in radians per Julian year, but the
-**     parallax is in arcseconds.
-**
-**  4) The RA proper motion is in terms of coordinate angle, not true
-**     angle.  If the catalog uses arcseconds for both RA and Dec proper
-**     motions, the RA proper motion will need to be divided by cos(Dec)
-**     before use.
-**
-**  5) Straight-line motion at constant speed, in the inertial frame,
-**     is assumed.
-**
-**  6) An extremely small (or zero or negative) parallax is interpreted
-**     to mean that the object is on the "celestial sphere", the radius
-**     of which is an arbitrary (large) value (see the iauStarpv
-**     function for the value used).  When the distance is overridden in
-**     this way, the status, initially zero, has 1 added to it.
-**
-**  7) If the space velocity is a significant fraction of c (see the
-**     constant VMAX in the function iauStarpv),  it is arbitrarily set
-**     to zero.  When this action occurs, 2 is added to the status.
-**
-**  8) The relativistic adjustment carried out in the iauStarpv function
-**     involves an iterative calculation.  If the process fails to
-**     converge within a set number of iterations, 4 is added to the
-**     status.
-**
-**  Called:
-**     iauStarpv    star catalog data to space motion pv-vector
-**     iauPvu       update a pv-vector
-**     iauPdp       scalar product of two p-vectors
-**     iauPvstar    space motion pv-vector to star catalog data
-**
-**  This revision:  2008 May 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double pv1[2][3], tl1, dt, pv[2][3], r2, rdv, v2, c2mv2, tl2,
-          pv2[2][3];
-   int j1, j2, j;
-
-
-/* RA,Dec etc. at the "before" epoch to space motion pv-vector. */
-   j1 = iauStarpv(ra1, dec1, pmr1, pmd1, px1, rv1, pv1);
-
-/* Light time when observed (days). */
-   tl1 = iauPm(pv1[0]) / DC;
-
-/* Time interval, "before" to "after" (days). */
-   dt = (ep2a - ep1a) + (ep2b - ep1b);
-
-/* Move star along track from the "before" observed position to the */
-/* "after" geometric position. */
-   iauPvu(dt + tl1, pv1, pv);
-
-/* From this geometric position, deduce the observed light time (days) */
-/* at the "after" epoch (with theoretically unneccessary error check). */
-   r2 = iauPdp(pv[0], pv[0]);
-   rdv = iauPdp(pv[0], pv[1]);
-   v2 = iauPdp(pv[1], pv[1]);
-   c2mv2 = DC*DC - v2;
-   if (c2mv2 <=  0) return -1;
-   tl2 = (-rdv + sqrt(rdv*rdv + c2mv2*r2)) / c2mv2;
-
-/* Move the position along track from the observed place at the */
-/* "before" epoch to the observed place at the "after" epoch. */
-   iauPvu(dt + (tl1 - tl2), pv1, pv2);
-
-/* Space motion pv-vector to RA,Dec etc. at the "after" epoch. */
-   j2 = iauPvstar(pv2, ra2, dec2, pmr2, pmd2, px2, rv2);
-
-/* Final status. */
-   j = (j2 == 0) ? j1 : -1;
-
-   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.
-**
-**--------------------------------------------------------------------*/
-}