255 lines
11 KiB
C
255 lines
11 KiB
C
|
|
/************************************************************
|
|
$Header: /home/dieter/sweph/RCS/swephlib.h,v 1.74 2008/06/16 10:07:20 dieter Exp $
|
|
|
|
Authors: Dieter Koch and Alois Treindl, Astrodienst Zurich
|
|
|
|
************************************************************/
|
|
/* Copyright (C) 1997 - 2008 Astrodienst AG, Switzerland. All rights reserved.
|
|
|
|
License conditions
|
|
------------------
|
|
|
|
This file is part of Swiss Ephemeris.
|
|
|
|
Swiss Ephemeris is distributed with NO WARRANTY OF ANY KIND. No author
|
|
or distributor accepts any responsibility for the consequences of using it,
|
|
or for whether it serves any particular purpose or works at all, unless he
|
|
or she says so in writing.
|
|
|
|
Swiss Ephemeris is made available by its authors under a dual licensing
|
|
system. The software developer, who uses any part of Swiss Ephemeris
|
|
in his or her software, must choose between one of the two license models,
|
|
which are
|
|
a) GNU public license version 2 or later
|
|
b) Swiss Ephemeris Professional License
|
|
|
|
The choice must be made before the software developer distributes software
|
|
containing parts of Swiss Ephemeris to others, and before any public
|
|
service using the developed software is activated.
|
|
|
|
If the developer choses the GNU GPL software license, he or she must fulfill
|
|
the conditions of that license, which includes the obligation to place his
|
|
or her whole software project under the GNU GPL or a compatible license.
|
|
See http://www.gnu.org/licenses/old-licenses/gpl-2.0.html
|
|
|
|
If the developer choses the Swiss Ephemeris Professional license,
|
|
he must follow the instructions as found in http://www.astro.com/swisseph/
|
|
and purchase the Swiss Ephemeris Professional Edition from Astrodienst
|
|
and sign the corresponding license contract.
|
|
|
|
The License grants you the right to use, copy, modify and redistribute
|
|
Swiss Ephemeris, but only under certain conditions described in the License.
|
|
Among other things, the License requires that the copyright notices and
|
|
this notice be preserved on all copies.
|
|
|
|
Authors of the Swiss Ephemeris: Dieter Koch and Alois Treindl
|
|
|
|
The authors of Swiss Ephemeris have no control or influence over any of
|
|
the derived works, i.e. over software or services created by other
|
|
programmers which use Swiss Ephemeris functions.
|
|
|
|
The names of the authors or of the copyright holder (Astrodienst) must not
|
|
be used for promoting any software, product or service which uses or contains
|
|
the Swiss Ephemeris. This copyright notice is the ONLY place where the
|
|
names of the authors can legally appear, except in cases where they have
|
|
given special permission in writing.
|
|
|
|
The trademarks 'Swiss Ephemeris' and 'Swiss Ephemeris inside' may be used
|
|
for promoting such software, products or services.
|
|
*/
|
|
|
|
#define PREC_IAU_1976 1
|
|
#define PREC_IAU_2000 2
|
|
#define PREC_IAU_2006 3
|
|
#define PREC_BRETAGNON_2003 4
|
|
#define PREC_LASKAR_1986 5
|
|
#define PREC_SIMON_1994 6
|
|
#define PREC_WILLIAMS_1994 7
|
|
#define PREC_VONDRAK_2011 8
|
|
|
|
/* Precession coefficients for remote past and future.
|
|
* One of the following four defines must be true.
|
|
*/
|
|
#define USE_PREC_VONDRAK_2011 TRUE
|
|
#define USE_PREC_WILLIAMS_1994 FALSE
|
|
#define USE_PREC_SIMON_1994 FALSE
|
|
#define USE_PREC_LASKAR_1986 FALSE
|
|
#define USE_PREC_BRETAGNON_2003 FALSE
|
|
/* IAU precession 1976 or 2003 for recent centuries.
|
|
* only one of the following two defines may be TRUE */
|
|
#define USE_PREC_IAU_1976 FALSE
|
|
#define USE_PREC_IAU_2000 FALSE
|
|
#define USE_PREC_IAU_2006 FALSE /* precession model P03 */
|
|
#define PREC_IAU_1976_CTIES 2.0 /* J2000 +/- two centuries */
|
|
#define PREC_IAU_2000_CTIES 2.0 /* J2000 +/- two centuries */
|
|
/* we use P03 for whole ephemeris */
|
|
#define PREC_IAU_2006_CTIES 75.0 /* J2000 +/- 75 centuries */
|
|
|
|
/* choose between the following nutation models */
|
|
#define NUT_IAU_1980 FALSE
|
|
#define NUT_IAU_2000A FALSE /* very time consuming ! */
|
|
#define NUT_IAU_2000B TRUE /* fast, but precision of milli-arcsec */
|
|
/* Set TRUE, to include Herring's (1987) corrections to IAU 1980
|
|
* nutation series. AA (1996) neglects them. */
|
|
#define NUT_CORR_1987 FALSE
|
|
|
|
/* frame bias */
|
|
#define FRAME_BIAS_IAU2006 TRUE /* if false, frame bias iau2000 will be used
|
|
* difference is minimal. */
|
|
|
|
/* For reproducing JPL Horizons to 2 mas (SEFLG_JPLHOR):
|
|
* The user has to keep the following files up to date which contain
|
|
* the earth orientation parameters related to the IAU 1980 nutation
|
|
* theory.
|
|
* Download the file
|
|
* datacenter.iers.org/eop/-/somos/5Rgv/document/tx13iers.u24/eopc04_08.62-now
|
|
* and rename it as eop_1962_today.txt. For current data and estimations for
|
|
* the near future, also download maia.usno.navy.mil/ser7/finals.all and
|
|
* rename it as eop_finals.txt */
|
|
#define DPSI_DEPS_IAU1980_FILE_EOPC04 "eop_1962_today.txt"
|
|
#define DPSI_DEPS_IAU1980_FILE_FINALS "eop_finals.txt"
|
|
#define DPSI_DEPS_IAU1980_TJD0_HORIZONS 2437684.5
|
|
#define HORIZONS_TJD0_DPSI_DEPS_IAU1980 2437684.5
|
|
#define INCLUDE_CODE_FOR_DPSI_DEPS_IAU1980 TRUE
|
|
/* You can set the latter false if you do not want to compile the
|
|
* code required to reproduce JPL Horizons.
|
|
* Keep it TRUE in order to reproduce JPL Horizons following
|
|
* IERS Conventions 1996 (1992), p. 22. Call swe_calc_ut() with
|
|
* iflag|SEFLG_JPLHOR. This options runs only, if the files
|
|
* DPSI_DEPS_IAU1980_FILE_EOPC04 and DPSI_DEPS_IAU1980_FILE_FINALS
|
|
* are in the ephemeris path.
|
|
*/
|
|
|
|
/* If the above define INCLUDE_CODE_FOR_DPSI_DEPS_IAU1980 is FALSE or
|
|
* the software does not find the earth orientation files (see above)
|
|
* in the ephemeris path, then SEFLG_JPLHOR will run as
|
|
* SEFLG_JPLHOR_APPROX.
|
|
* The following define APPROXIMATE_HORIZONS_ASTRODIENST defines
|
|
* the handling of SEFLG_JPLHOR_APPROX.
|
|
* With this flag, planetary positions are always calculated
|
|
* using a recent precession/nutation model.
|
|
* If APPROXIMATE_HORIZONS_ASTRODIENST is FALSE, then the
|
|
* frame bias as recommended by IERS Conventions 2003 and 2010
|
|
* is *not* applied. Instead, dpsi_bias and deps_bias are added to
|
|
* nutation. This procedure is found in some older astronomical software.
|
|
* Equatorial apparent positions will be close to JPL Horizons
|
|
* (within a few mas) beetween 1962 and current years. Ecl. longitude
|
|
* will be good, latitude bad.
|
|
* If APPROXIMATE_HORIZONS_ASTRODIENST is TRUE, the approximation of
|
|
* JPL Horizons is even better. Frame bias matrix is applied with
|
|
* some correction to RA and another correction is added to epsilon.
|
|
*/
|
|
#define APPROXIMATE_HORIZONS_ASTRODIENST TRUE
|
|
|
|
#define USE_HORIZONS_METHOD_BEFORE_1980 TRUE /* Horizons method before 20-jan-1962 */
|
|
/* The latter, if combined with SEFLG_JPLHOR provides good agreement
|
|
* with JPL Horizons for 1800 - today. However, Horizons uses correct
|
|
* dpsi and deps only after 20-jan-1962. For all dates before that
|
|
* it uses dpsi and deps of 20-jan-1962, which provides a continuous
|
|
* ephemeris, but does not make sense otherwise.
|
|
* Before 1800, even this option does not provide agreement with Horizons,
|
|
* because Horizons uses a different precession model (Owen 1986)
|
|
* before 1800, which is not included in the Swiss Ephemeris.
|
|
* If this macro is FALSE then the program defaults to SEFLG_JPLHOR_APPROX
|
|
* outside the time range of correction data dpsi and deps.
|
|
* Note that this will result in a non-continuous ephemeris near
|
|
* 20-jan-1962 and current years.
|
|
*/
|
|
|
|
/* coordinate transformation */
|
|
extern void swi_coortrf(double *xpo, double *xpn, double eps);
|
|
|
|
/* coordinate transformation */
|
|
extern void swi_coortrf2(double *xpo, double *xpn, double sineps, double coseps);
|
|
|
|
/* cartesian to polar coordinates */
|
|
extern void swi_cartpol(double *x, double *l);
|
|
|
|
/* cartesian to polar coordinates with velocity */
|
|
extern void swi_cartpol_sp(double *x, double *l);
|
|
extern void swi_polcart_sp(double *l, double *x);
|
|
|
|
/* polar to cartesian coordinates */
|
|
extern void swi_polcart(double *l, double *x);
|
|
|
|
/* GCRS to J2000 */
|
|
extern void swi_bias(double *x, double tjd, int32 iflag, AS_BOOL backward);
|
|
extern void swi_get_eop_time_range(void);
|
|
/* GCRS to FK5 */
|
|
extern void swi_icrs2fk5(double *x, int32 iflag, AS_BOOL backward);
|
|
|
|
/* precession */
|
|
extern int swi_precess(double *R, double J, int32 iflag, int direction );
|
|
extern void swi_precess_speed(double *xx, double t, int32 iflag, int direction);
|
|
|
|
/* from sweph.c, light deflection, aberration, etc. */
|
|
extern void swi_deflect_light(double *xx, double dt, int32 iflag);
|
|
extern void swi_aberr_light(double *xx, double *xe, int32 iflag);
|
|
extern int swi_plan_for_osc_elem(int32 iflag, double tjd, double *xx);
|
|
extern int swi_trop_ra2sid_lon(double *xin, double *xout, double *xoutr, int32 iflag, char *serr);
|
|
extern int swi_trop_ra2sid_lon_sosy(double *xin, double *xout, double *xoutr, int32 iflag, char *serr);
|
|
extern int swi_get_observer(double tjd, int32 iflag,
|
|
AS_BOOL do_save, double *xobs, char *serr);
|
|
extern void swi_force_app_pos_etc();
|
|
|
|
/* obliquity of ecliptic */
|
|
extern void swi_check_ecliptic(double tjd, int32 iflag);
|
|
extern double swi_epsiln(double J, int32 iflag);
|
|
extern void swi_ldp_peps(double J, double *dpre, double *deps);
|
|
|
|
/* nutation */
|
|
extern void swi_check_nutation(double tjd, int32 iflag);
|
|
extern int swi_nutation(double J, int32 iflag, double *nutlo);
|
|
extern void swi_nutate(double *xx, int32 iflag, AS_BOOL backward);
|
|
|
|
extern void swi_mean_lunar_elements(double tjd,
|
|
double *node, double *dnode,
|
|
double *peri, double *dperi);
|
|
/* */
|
|
extern double swi_mod2PI(double x);
|
|
|
|
/* evaluation of chebyshew series and derivative */
|
|
extern double swi_echeb(double x, double *coef, int ncf);
|
|
extern double swi_edcheb(double x, double *coef, int ncf);
|
|
|
|
/* cross product of vectors */
|
|
extern void swi_cross_prod(double *a, double *b, double *x);
|
|
/* dot product of vecotrs */
|
|
extern double swi_dot_prod_unit(double *x, double *y);
|
|
|
|
extern double swi_angnorm(double x);
|
|
|
|
/* generation of SWISSEPH file names */
|
|
extern void swi_gen_filename(double tjd, int ipli, char *fname);
|
|
|
|
/* cyclic redundancy checksum (CRC), 32 bit */
|
|
extern uint32 swi_crc32(unsigned char *buf, int len);
|
|
|
|
extern int swi_cutstr(char *s, char *cutlist, char *cpos[], int nmax);
|
|
extern char *swi_right_trim(char *s);
|
|
|
|
extern double swi_kepler(double E, double M, double ecce);
|
|
|
|
extern char *swi_get_fict_name(int32 ipl, char *s);
|
|
|
|
extern void swi_FK4_FK5(double *xp, double tjd);
|
|
|
|
extern char *swi_strcpy(char *to, char *from);
|
|
extern char *swi_strncpy(char *to, char *from, size_t n);
|
|
|
|
extern double swi_deltat_ephe(double tjd_ut, int32 epheflag);
|
|
|
|
#ifdef TRACE
|
|
# define TRACE_COUNT_MAX 10000
|
|
extern FILE *swi_fp_trace_c;
|
|
extern FILE *swi_fp_trace_out;
|
|
extern void swi_open_trace(char *serr);
|
|
extern int32 swi_trace_count;
|
|
static char *fname_trace_c = "swetrace.c";
|
|
static char *fname_trace_out = "swetrace.txt";
|
|
#ifdef FORCE_IFLAG
|
|
static char *fname_force_flg = "force.flg";
|
|
#endif
|
|
#endif /* TRACE */
|