equationoftime 4.0.0 → 4.1.1

data/examples/celes_parts/faur03.c

@@ -1,99 +0,0 @@
-#include "sofam.h"
-
-double iauFaur03(double t)
-/*
-**  - - - - - - - - - -
-**   i a u F a u r 0 3
-**  - - - - - - - - - -
-**
-**  Fundamental argument, IERS Conventions (2003):
-**  mean longitude of Uranus.
-**
-**  Status:  canonical model.
-**
-**  Given:
-**     t     double    TDB, Julian centuries since J2000.0 (Note 1)
-**
-**  Returned  (function value):
-**           double    mean longitude of Uranus, radians (Note 2)
-**
-**  Notes:
-**
-**  1) Though t is strictly TDB, it is usually more convenient to use
-**     TT, which makes no significant difference.
-**
-**  2) The expression used is as adopted in IERS Conventions (2003) and
-**     is adapted from Simon et al. (1994).
-**
-**  References:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**     Simon, J.-L., Bretagnon, P., Chapront, J., Chapront-Touze, M.,
-**     Francou, G., Laskar, J. 1994, Astron.Astrophys. 282, 663-683
-**
-**  This revision:  2009 December 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double a;
-
-
-/* Mean longitude of Uranus (IERS Conventions 2003). */
-   a = fmod(5.481293872 + 7.4781598567 * t, D2PI);
-
-   return a;
-
-/*----------------------------------------------------------------------
-**
-**  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/fave03.c

@@ -1,102 +0,0 @@
-#include "sofam.h"
-
-double iauFave03(double t)
-/*
-**  - - - - - - - - - -
-**   i a u F a v e 0 3
-**  - - - - - - - - - -
-**
-**  Fundamental argument, IERS Conventions (2003):
-**  mean longitude of Venus.
-**
-**  Status:  canonical model.
-**
-**  Given:
-**     t     double    TDB, Julian centuries since J2000.0 (Note 1)
-**
-**  Returned (function value):
-**           double    mean longitude of Venus, radians (Note 2)
-**
-**  Notes:
-**
-**  1) Though t is strictly TDB, it is usually more convenient to use
-**     TT, which makes no significant difference.
-**
-**  2) The expression used is as adopted in IERS Conventions (2003) and
-**     comes from Souchay et al. (1999) after Simon et al. (1994).
-**
-**  References:
-**
-**     McCarthy, D. D., Petit, G. (eds.), IERS Conventions (2003),
-**     IERS Technical Note No. 32, BKG (2004)
-**
-**     Simon, J.-L., Bretagnon, P., Chapront, J., Chapront-Touze, M.,
-**     Francou, G., Laskar, J. 1994, Astron.Astrophys. 282, 663-683
-**
-**     Souchay, J., Loysel, B., Kinoshita, H., Folgueira, M. 1999,
-**     Astron.Astrophys.Supp.Ser. 135, 111
-**
-**  This revision:  2009 December 16
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double a;
-
-
-/* Mean longitude of Venus (IERS Conventions 2003). */
-   a = fmod(3.176146697 + 1021.3285546211 * t, D2PI);
-
-   return a;
-
-/*----------------------------------------------------------------------
-**
-**  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/fk52h.c

@@ -1,143 +0,0 @@
-#include "sofam.h"
-
-void iauFk52h(double r5, double d5,
-              double dr5, double dd5, double px5, double rv5,
-              double *rh, double *dh,
-              double *drh, double *ddh, double *pxh, double *rvh)
-/*
-**  - - - - - - - - -
-**   i a u F k 5 2 h
-**  - - - - - - - - -
-**
-**  Transform FK5 (J2000.0) star data into the Hipparcos system.
-**
-**  Status:  support function.
-**
-**  Given (all FK5, equinox J2000.0, epoch J2000.0):
-**     r5      double    RA (radians)
-**     d5      double    Dec (radians)
-**     dr5     double    proper motion in RA (dRA/dt, rad/Jyear)
-**     dd5     double    proper motion in Dec (dDec/dt, rad/Jyear)
-**     px5     double    parallax (arcsec)
-**     rv5     double    radial velocity (km/s, positive = receding)
-**
-**  Returned (all Hipparcos, epoch J2000.0):
-**     rh      double    RA (radians)
-**     dh      double    Dec (radians)
-**     drh     double    proper motion in RA (dRA/dt, rad/Jyear)
-**     ddh     double    proper motion in Dec (dDec/dt, rad/Jyear)
-**     pxh     double    parallax (arcsec)
-**     rvh     double    radial velocity (km/s, positive = receding)
-**
-**  Notes:
-**
-**  1) This function transforms FK5 star positions and proper motions
-**     into the system of the Hipparcos catalog.
-**
-**  2) The proper motions in RA are dRA/dt rather than
-**     cos(Dec)*dRA/dt, and are per year rather than per century.
-**
-**  3) The FK5 to Hipparcos transformation is modeled as a pure
-**     rotation and spin;  zonal errors in the FK5 catalog are not
-**     taken into account.
-**
-**  4) See also iauH2fk5, iauFk5hz, iauHfk5z.
-**
-**  Called:
-**     iauStarpv    star catalog data to space motion pv-vector
-**     iauFk5hip    FK5 to Hipparcos rotation and spin
-**     iauRxp       product of r-matrix and p-vector
-**     iauPxp       vector product of two p-vectors
-**     iauPpp       p-vector plus p-vector
-**     iauPvstar    space motion pv-vector to star catalog data
-**
-**  Reference:
-**
-**     F.Mignard & M.Froeschle, Astron. Astrophys. 354, 732-739 (2000).
-**
-**  This revision:  2009 December 17
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   int i;
-   double pv5[2][3], r5h[3][3], s5h[3], wxp[3], vv[3], pvh[2][3];
-
-
-/* FK5 barycentric position/velocity pv-vector (normalized). */
-   iauStarpv(r5, d5, dr5, dd5, px5, rv5, pv5);
-
-/* FK5 to Hipparcos orientation matrix and spin vector. */
-   iauFk5hip(r5h, s5h);
-
-/* Make spin units per day instead of per year. */
-   for ( i = 0; i < 3; s5h[i++] /= 365.25 );
-
-/* Orient the FK5 position into the Hipparcos system. */
-   iauRxp(r5h, pv5[0], pvh[0]);
-
-/* Apply spin to the position giving an extra space motion component. */
-   iauPxp(pv5[0], s5h, wxp);
-
-/* Add this component to the FK5 space motion. */
-   iauPpp(wxp, pv5[1], vv);
-
-/* Orient the FK5 space motion into the Hipparcos system. */
-   iauRxp(r5h, vv, pvh[1]);
-
-/* Hipparcos pv-vector to spherical. */
-   iauPvstar(pvh, rh, dh, drh, ddh, pxh, rvh);
-
-   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/fk5hip.c

@@ -1,126 +0,0 @@
-#include "sofam.h"
-
-void iauFk5hip(double r5h[3][3], double s5h[3])
-/*
-**  - - - - - - - - - -
-**   i a u F k 5 h i p
-**  - - - - - - - - - -
-**
-**  FK5 to Hipparcos rotation and spin.
-**
-**  Status:  support function.
-**
-**  Returned:
-**     r5h   double[3][3]  r-matrix: FK5 rotation wrt Hipparcos (Note 2)
-**     s5h   double[3]     r-vector: FK5 spin wrt Hipparcos (Note 3)
-**
-**  Notes:
-**
-**  1) This function models the FK5 to Hipparcos transformation as a
-**     pure rotation and spin;  zonal errors in the FK5 catalogue are
-**     not taken into account.
-**
-**  2) The r-matrix r5h operates in the sense:
-**
-**           P_Hipparcos = r5h x P_FK5
-**
-**     where P_FK5 is a p-vector in the FK5 frame, and P_Hipparcos is
-**     the equivalent Hipparcos p-vector.
-**
-**  3) The r-vector s5h represents the time derivative of the FK5 to
-**     Hipparcos rotation.  The units are radians per year (Julian,
-**     TDB).
-**
-**  Called:
-**     iauRv2m      r-vector to r-matrix
-**
-**  Reference:
-**
-**     F.Mignard & M.Froeschle, Astron. Astrophys. 354, 732-739 (2000).
-**
-**  This revision:  2009 March 14
-**
-**  Original version 2012-03-01
-**
-**  Copyright (C) 2013 Naoki Arita.  See notes at end.
-*/
-{
-   double v[3];
-
-/* FK5 wrt Hipparcos orientation and spin (radians, radians/year) */
-   double epx, epy, epz;
-   double omx, omy, omz;
-
-
-   epx = -19.9e-3 * DAS2R;
-   epy =  -9.1e-3 * DAS2R;
-   epz =  22.9e-3 * DAS2R;
-
-   omx = -0.30e-3 * DAS2R;
-   omy =  0.60e-3 * DAS2R;
-   omz =  0.70e-3 * DAS2R;
-
-/* FK5 to Hipparcos orientation expressed as an r-vector. */
-   v[0] = epx;
-   v[1] = epy;
-   v[2] = epz;
-
-/* Re-express as an r-matrix. */
-   iauRv2m(v, r5h);
-
-/* Hipparcos wrt FK5 spin expressed as an r-vector. */
-   s5h[0] = omx;
-   s5h[1] = omy;
-   s5h[2] = omz;
-
-   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.
-**
-**--------------------------------------------------------------------*/
-}