summaryrefslogtreecommitdiff
path: root/time/solar_declination.c
diff options
context:
space:
mode:
Diffstat (limited to 'time/solar_declination.c')
-rw-r--r--time/solar_declination.c77
1 files changed, 0 insertions, 77 deletions
diff --git a/time/solar_declination.c b/time/solar_declination.c
deleted file mode 100644
index dd89028..0000000
--- a/time/solar_declination.c
+++ /dev/null
@@ -1,77 +0,0 @@
-/*
- * (C)2012 Michael Duane Rice All rights reserved.
- *
- * Redistribution and use in source and binary forms, with or without
- * modification, are permitted provided that the following conditions are
- * met:
- *
- * Redistributions of source code must retain the above copyright notice, this
- * list of conditions and the following disclaimer. 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. Neither the name of the copyright holders
- * nor the names of 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 OWNER 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.
- */
-
-/* $Id: solar_declination.c 2369 2013-04-28 14:19:35Z swfltek $ */
-
-/*
- Were it not for the eccentricity of Earths orbit, this would be a trivial function.
-
- We compute the Earths orbital position with respect to perihelion, from which we derive a
- 'velocity correction factor'. We then compute the orbital angle with respect to the
- December solstice, as 'modulated' by that correction factor.
-
- Due to the accumulation of rounding errors, the computed December solstice of 2135 will lag
- the actual solstice by many hours. A fudge factor, 'LAG', distributes the error across
- the 136 year range of this library.
-*/
-
-#include <time.h>
-#include <math.h>
-#include "ephemera_common.h"
-
-#define LAG 38520
-
-double
-solar_declination(const time_t * timer)
-{
-
- uint32_t fT, oV;
- double dV, dT;
-
- /* Determine orbital angle relative to perihelion of January 1999 */
- oV = *timer % ANOM_YEAR;
- oV += PERIHELION;
- dV = oV;
- dV /= ANOM_CYCLE;
-
- /* Derive velocity correction factor from the perihelion angle */
- dV = sin(dV);
- dV *= DELTA_V;
-
- /* Determine orbital angle relative to solstice of December 1999 */
- fT = *timer % TROP_YEAR;
- fT += SOLSTICE + LAG;
- dT = fT;
- dT /= TROP_CYCLE;
- dT += dV;
-
- /* Finally having the solstice angle, we can compute the declination */
- dT = cos(dT) * INCLINATION;
-
- return -dT;
-}