swe4r 0.0.0

data/ext/swe4r/sweph.h

@@ -0,0 +1,556 @@
+/************************************************************
+   $Header: /home/dieter/sweph/RCS/sweph.h,v 1.74 2008/06/16 10:07:20 dieter Exp $
+   definitions and constants SWISSEPH
+
+  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.
+*/
+
+/*
+ * move over from swephexp.h
+ */
+
+#define SE_VERSION      "1.77.00"
+
+#define J2000           2451545.0  	/* 2000 January 1.5 */
+#define B1950           2433282.42345905  	/* 1950 January 0.923 */
+#define J1900           2415020.0  	/* 1900 January 0.5 */
+
+#define MPC_CERES       1
+#define MPC_PALLAS      2
+#define MPC_JUNO        3
+#define MPC_VESTA       4
+#define MPC_CHIRON      2060
+#define MPC_PHOLUS      5145
+
+#define SE_NAME_SUN             "Sun"
+#define SE_NAME_MOON            "Moon"
+#define SE_NAME_MERCURY         "Mercury"
+#define SE_NAME_VENUS           "Venus"
+#define SE_NAME_MARS            "Mars"
+#define SE_NAME_JUPITER         "Jupiter"
+#define SE_NAME_SATURN          "Saturn"
+#define SE_NAME_URANUS          "Uranus"
+#define SE_NAME_NEPTUNE         "Neptune"
+#define SE_NAME_PLUTO           "Pluto"
+#define SE_NAME_MEAN_NODE       "mean Node"
+#define SE_NAME_TRUE_NODE       "true Node"
+#define SE_NAME_MEAN_APOG       "mean Apogee"
+#define SE_NAME_OSCU_APOG       "osc. Apogee"
+#define SE_NAME_INTP_APOG       "intp. Apogee"
+#define SE_NAME_INTP_PERG       "intp. Perigee"
+#define SE_NAME_EARTH           "Earth"
+#define SE_NAME_CERES           "Ceres"
+#define SE_NAME_PALLAS          "Pallas"
+#define SE_NAME_JUNO            "Juno"
+#define SE_NAME_VESTA           "Vesta"
+#define SE_NAME_CHIRON          "Chiron"
+#define SE_NAME_PHOLUS          "Pholus"
+
+ 
+#define SE_NAME_CUPIDO          "Cupido"
+#define SE_NAME_HADES           "Hades"
+#define SE_NAME_ZEUS            "Zeus"
+#define SE_NAME_KRONOS          "Kronos"
+#define SE_NAME_APOLLON         "Apollon"
+#define SE_NAME_ADMETOS         "Admetos"
+#define SE_NAME_VULKANUS        "Vulkanus"
+#define SE_NAME_POSEIDON        "Poseidon"
+#define SE_NAME_ISIS            "Isis"
+#define SE_NAME_NIBIRU          "Nibiru"
+#define SE_NAME_HARRINGTON      "Harrington"
+#define SE_NAME_NEPTUNE_LEVERRIER       "Leverrier"
+#define SE_NAME_NEPTUNE_ADAMS   "Adams"
+#define SE_NAME_PLUTO_LOWELL    "Lowell"
+#define SE_NAME_PLUTO_PICKERING "Pickering"
+#define SE_NAME_VULCAN          "Vulcan"
+#define SE_NAME_WHITE_MOON      "White Moon"
+
+/* for delta t: intrinsic tidal acceleration in the mean motion of the moon,
+ * not given in the parameters list of the ephemeris files but computed
+ * by Chapront/Chapront-Touzé/Francou A&A 387 (2002), p. 705.
+ */
+#define SE_TIDAL_DE200          (-23.8946)
+#define SE_TIDAL_DE403          (-25.580)  /* was (-25.8) until V. 1.76.2 */
+#define SE_TIDAL_DE404          (-25.580)  /* was (-25.8) until V. 1.76.2 */
+#define SE_TIDAL_DE405          (-25.826)  /* was (-25.7376) until V. 1.76.2 */
+#define SE_TIDAL_DE406          (-25.826)  /* was (-25.7376) until V. 1.76.2 */
+
+#define SE_TIDAL_26             (-26.0)
+
+#define SE_TIDAL_DEFAULT        SE_TIDAL_DE406
+
+/*
+ * earlier content
+ */
+
+#define PI              M_PI	/* 3.14159265358979323846, math.h */
+#define TWOPI           (2.0 * PI)
+
+#define ENDMARK         -99
+
+#define SEI_EPSILON     -2
+#define SEI_NUTATION    -1
+#define SEI_EMB		0	
+#define SEI_EARTH	0	
+#define SEI_SUN  	0	
+#define	SEI_MOON	1	
+#define	SEI_MERCURY	2	
+#define	SEI_VENUS	3	
+#define	SEI_MARS	4	
+#define	SEI_JUPITER	5	
+#define	SEI_SATURN	6	
+#define	SEI_URANUS	7	
+#define	SEI_NEPTUNE	8	
+#define	SEI_PLUTO	9	
+#define	SEI_SUNBARY	10	/* barycentric sun */
+#define	SEI_ANYBODY	11	/* any asteroid */
+#define	SEI_CHIRON	12	
+#define	SEI_PHOLUS	13	
+#define	SEI_CERES	14	
+#define	SEI_PALLAS	15	
+#define	SEI_JUNO	16	
+#define	SEI_VESTA	17	
+
+#define SEI_NPLANETS    18
+
+#define SEI_MEAN_NODE   0
+#define SEI_TRUE_NODE   1
+#define SEI_MEAN_APOG   2
+#define SEI_OSCU_APOG   3
+#define SEI_INTP_APOG   4
+#define SEI_INTP_PERG   5
+
+#define SEI_NNODE_ETC    6
+
+#define SEI_FLG_HELIO   1
+#define SEI_FLG_ROTATE  2
+#define SEI_FLG_ELLIPSE 4
+#define SEI_FLG_EMBHEL  8   	/* TRUE, if heliocentric earth is given
+				 * instead of barycentric sun 
+				 * i.e. bary sun is computed from 
+				 * barycentric and heliocentric earth */
+
+#define SEI_FILE_PLANET	  0
+#define SEI_FILE_MOON	  1 
+#define SEI_FILE_MAIN_AST 2
+#define SEI_FILE_ANY_AST  3
+#define SEI_FILE_FIXSTAR  4
+
+#if 0
+#define SEI_FILE_TEST_ENDIAN     (97L * 65536L + 98L * 256L + 99L) /*abc*/
+#endif
+#define SEI_FILE_TEST_ENDIAN     (0x616263L) 	/* abc*/ 
+#define SEI_FILE_BIGENDIAN	0
+#define SEI_FILE_NOREORD	0
+#define SEI_FILE_LITENDIAN	1
+#define SEI_FILE_REORD  	2
+
+#define SEI_FILE_NMAXPLAN	50
+#define SEI_FILE_EFPOSBEGIN      500
+
+#define SE_FILE_SUFFIX	"se1"
+
+#define SEI_NEPHFILES   7
+#define SEI_CURR_FPOS   -1
+
+/* Chiron's orbit becomes chaotic 
+ * before 720 AD and after 4606 AD, because of close encounters
+ * with Saturn. Accepting a maximum error of 5 degrees, 
+ * the ephemeris is good between the following dates:
+ */
+#define CHIRON_START    1958470.5  	/* 1.1.650 */
+#define CHIRON_END      3419437.5  	/* 1.1.4650 */
+
+/* Pholus's orbit is unstable as well, because he sometimes
+ * approaches Saturn.
+ * Accepting a maximum error of 5 degrees,
+ * the ephemeris is good after the following date:
+ */
+#define PHOLUS_START    314845.5  	/* 1.1.-3850 */
+
+#define MOSHPLEPH_START	 625000.5
+#define MOSHPLEPH_END  	2818000.5
+#define MOSHLUEPH_START	 625000.5
+#define MOSHLUEPH_END  	2818000.5
+#define MOSHNDEPH_START	-254900.5	/* 14 Feb -5410 00:00 ET jul.cal.*/
+#define MOSHNDEPH_END  	3697000.5	/* 11 Dec 5409 00:00 ET, greg. cal */
+/*
+#define MOSHPLEPH_START	 -225000.5
+#define MOSHPLEPH_END  	3600000.5
+#define MOSHLUEPH_START	 -225000.5
+#define MOSHLUEPH_END  	3600000.5
+*/
+#if FALSE	/*	Alois commented out, not used anywhere  */
+#define JPLEPH_START	 625307.5	/* about -3000 (DE406) */
+#define JPLEPH_END	2816848.5	/* about  3000 (DE406) */
+#define SWIEPH_START	 625614.927151
+#define SWIEPH_END	2813641.5
+#define ALLEPH_START	MOSHPLEPH_START
+#define ALLEPH_END	MOSHPLEPH_END
+#define BEG_YEAR       (-3000)
+#define END_YEAR       3000
+#endif
+
+#define MAXORD          40
+
+#define NCTIES         6.0     /* number of centuries per eph. file */
+
+#define OK (0)
+#define ERR (-1)
+#define NOT_AVAILABLE (-2)
+#define BEYOND_EPH_LIMITS (-3)
+
+#define J_TO_J2000   	1
+#define J2000_TO_J   	-1
+
+/* we always use Astronomical Almanac constants, if available */
+#define MOON_MEAN_DIST  384400000.0		/* in m, AA 1996, F2 */
+#define MOON_MEAN_INCL  5.1453964		/* AA 1996, D2 */
+#define MOON_MEAN_ECC   0.054900489		/* AA 1996, F2 */
+/* #define SUN_EARTH_MRAT  328900.561400           Su/(Ea+Mo) AA 2006 K7 */
+#define SUN_EARTH_MRAT  332946.050895           /* Su / (Ea only) AA 2006 K7 */   
+#define EARTH_MOON_MRAT (1 / 0.0123000383)	/* AA 2006, K7 */
+#if 0
+#define EARTH_MOON_MRAT 81.30056		/* de406 */
+#endif
+#define AUNIT       	1.49597870691e+11  	/* au in meters, AA 2006 K6 */
+#define CLIGHT       	2.99792458e+8   	/* m/s, AA 1996 K6 */
+#if 0
+#define HELGRAVCONST    1.32712438e+20		/* G * M(sun), m^3/sec^2, AA 1996 K6 */
+#endif
+#define HELGRAVCONST    1.32712440017987e+20	/* G * M(sun), m^3/sec^2, AA 2006 K6 */
+#define GEOGCONST       3.98600448e+14 		/* G * M(earth) m^3/sec^2, AA 1996 K6 */
+#define KGAUSS		0.01720209895		/* Gaussian gravitational constant K6 */
+#define SUN_RADIUS      (959.63 / 3600 * DEGTORAD)  /*  Meeus germ. p 391 */
+#define EARTH_RADIUS	6378136.6		/* AA 2006 K6 */
+/*#define EARTH_OBLATENESS (1.0/ 298.257223563)	 * AA 1998 K13 */
+#define EARTH_OBLATENESS (1.0/ 298.25642)	/* AA 2006 K6 */
+#define EARTH_ROT_SPEED (7.2921151467e-5 * 86400) /* in rad/day, expl. suppl., p 162 */
+
+#define LIGHTTIME_AUNIT  (499.0047838061/3600/24) 	/* 8.3167 minutes (days), AA 2006 K6 */
+
+/* node of ecliptic measured on ecliptic 2000 */
+#define SSY_PLANE_NODE_E2000    (107.582569 * DEGTORAD)
+/* node of ecliptic measured on solar system rotation plane */
+#define SSY_PLANE_NODE          (107.58883388 * DEGTORAD)
+/* inclination of ecliptic against solar system rotation plane */
+#define SSY_PLANE_INCL          (1.578701 * DEGTORAD)
+
+#define KM_S_TO_AU_CTY	 21.095			/* km/s to AU/century */
+#define MOON_SPEED_INTV  0.00005 		/* 4.32 seconds (in days) */
+#define PLAN_SPEED_INTV  0.0001 	        /* 8.64 seconds (in days) */
+#define MEAN_NODE_SPEED_INTV  0.001		
+#define NODE_CALC_INTV  0.0001		
+#define NODE_CALC_INTV_MOSH   0.1		
+#define NUT_SPEED_INTV   0.0001
+#define DEFL_SPEED_INTV  0.0000005
+
+#define SE_LAPSE_RATE        0.0065  /* deg K / m, for refraction */
+
+#define square_sum(x)   (x[0]*x[0]+x[1]*x[1]+x[2]*x[2])
+#define dot_prod(x,y)   (x[0]*y[0]+x[1]*y[1]+x[2]*y[2])
+
+#define PNOINT2JPL {J_EARTH, J_MOON, J_MERCURY, J_VENUS, J_MARS, J_JUPITER, J_SATURN, J_URANUS, J_NEPTUNE, J_PLUTO, J_SUN, }
+
+/* planetary radii in meters */
+#define NDIAM  (SE_VESTA + 1)
+static const double pla_diam[NDIAM] = {1392000000.0, /* Sun */
+                           3476300.0, /* Moon */
+                           2439000.0 * 2, /* Mercury */
+                           6052000.0 * 2, /* Venus */
+                           3397200.0 * 2, /* Mars */
+                          71398000.0 * 2, /* Jupiter */
+                          60000000.0 * 2, /* Saturn */
+                          25400000.0 * 2, /* Uranus */
+                          24300000.0 * 2, /* Neptune */
+                           2500000.0 * 2, /* Pluto */
+                           0, 0, 0, 0,    /* nodes and apogees */
+                           6378140.0 * 2, /* Earth */
+                                 0.0, /* Chiron */
+                                 0.0, /* Pholus */
+                            913000.0, /* Ceres */
+                            523000.0, /* Pallas */
+                            244000.0, /* Juno */
+                            501000.0, /* Vesta */
+                        };
+
+
+/* Ayanamsas */
+struct aya_init {double t0, ayan_t0;};
+static const struct aya_init ayanamsa[] = {
+    {2433282.5, 24.042044444},	/* 0: Fagan/Bradley (Default) */
+    {J1900, 360 - 337.53953},   /* 1: Lahiri (Robert Hand) */
+    {J1900, 360 - 333.58695},   /* 2: De Luce (Robert Hand) */
+    {J1900, 360 - 338.98556},   /* 3: Raman (Robert Hand) */
+    {J1900, 360 - 341.33904},   /* 4: Ushashashi (Robert Hand) */
+    {J1900, 360 - 337.636111},  /* 5: Krishnamurti (Robert Hand) */
+    {J1900, 360 - 333.0369024}, /* 6: Djwhal Khool; (Graham Dawson)  
+                                 *    Aquarius entered on 1 July 2117 */
+    {J1900, 360 - 338.917778},  /* 7: Yukteshwar; (David Cochrane) */
+    {J1900, 360 - 338.634444},  /* 8: JN Bhasin; (David Cochrane) */
+    {1684532.5, -3.36667},      /* 9: Babylonian, Kugler 1 */
+    {1684532.5, -4.76667},      /*10: Babylonian, Kugler 2 */
+    {1684532.5, -5.61667},      /*11: Babylonian, Kugler 3 */
+    {1684532.5, -4.56667},      /*12: Babylonian, Huber */
+    {1673941, -5.079167},       /*13: Babylonian, Mercier;
+                                 *    eta Piscium culminates with zero point */
+    {1684532.5, -4.44088389},   /*14: t0 is defined by Aldebaran at 15 Taurus */
+    {1674484, -9.33333},        /*15: Hipparchos */
+    {1927135.8747793, 0},       /*16: Sassanian */
+    /*{1746443.513, 0},          *17: Galactic Center at 0 Sagittarius */
+    {1746447.518, 0},           /*17: Galactic Center at 0 Sagittarius */
+    {J2000, 0},	                /*18: J2000 */
+    {J1900, 0},	                /*19: J1900 */
+    {B1950, 0},	                /*20: B1950 */
+	};
+
+#define PLAN_DATA struct plan_data
+
+/* obliquity of ecliptic */
+struct epsilon {
+  double teps, eps, seps, ceps; 	/* jd, eps, sin(eps), cos(eps) */
+};
+
+/*
+extern struct epsilon oec2000;
+extern struct epsilon oec;
+*/
+
+struct plan_data {
+  /* the following data are read from file only once, immediately after 
+   * file has been opened */
+  int ibdy;		/* internal body number */
+  int32 iflg;		/* contains several bit flags describing the data:	
+			 * SEI_FLG_HELIO: true if helio, false if bary
+			 * SEI_FLG_ROTATE: TRUE if coefficients are referred 
+			 *      to coordinate system of orbital plane 
+			 * SEI_FLG_ELLIPSE: TRUE if reference ellipse */
+  int ncoe;		/* # of coefficients of ephemeris polynomial,
+			   is polynomial order + 1  */
+  /* where is the segment index on the file */
+  int32 lndx0;   	/* file position of begin of planet's index */
+  int32 nndx;		/* number of index entries on file: computed */
+  double tfstart;	/* file contains ephemeris for tfstart thru tfend */
+  double tfend;         /*      for this particular planet !!!            */
+  double dseg;		/* segment size (days covered by a polynomial)  */
+  /* orbital elements: */
+  double telem;		/* epoch of elements */
+  double prot;		
+  double qrot; 
+  double dprot;
+  double dqrot; 
+  double rmax;		/* normalisation factor of cheby coefficients */
+  /* in addition, if reference ellipse is used: */
+  double peri;		
+  double dperi; 
+  double *refep;	/* pointer to cheby coeffs of reference ellipse,
+			 * size of data is 2 x ncoe */
+  /* unpacked segment information, only updated when a segment is read: */
+  double tseg0, tseg1;	/* start and end jd of current segment */
+  double *segp;         /* pointer to unpacked cheby coeffs of segment;
+			 * the size is 3 x ncoe */
+  int neval;		/* how many coefficients to evaluate. this may
+			 * be less than ncoe */
+  /* result of most recent data evaluation for this body: */
+  double teval;		/* time for which previous computation was made */
+  int32 iephe;            /* which ephemeris was used */
+  double x[6];		/* position and speed vectors equatorial J2000 */
+  int32 xflgs;		/* hel., light-time, aberr., prec. flags etc. */
+  double xreturn[24];   /* return positions:
+			 * xreturn+0	ecliptic polar coordinates
+			 * xreturn+6	ecliptic cartesian coordinates
+			 * xreturn+12	equatorial polar coordinates
+			 * xreturn+18	equatorial cartesian coordinates
+			 */
+};
+
+/*
+ * stuff exported from swemplan.c and swemmoon.c 
+ * and constants used inside these functions.
+************************************************************/
+
+#define STR             4.8481368110953599359e-6 /* radians per arc second */
+
+/* moon, s. moshmoon.c */
+extern int swi_mean_node(double jd, double *x, char *serr);
+extern int swi_mean_apog(double jd, double *x, char *serr);
+extern int swi_moshmoon(double tjd, AS_BOOL do_save, double *xpm, char *serr) ;
+extern int swi_moshmoon2(double jd, double *x);
+extern int swi_intp_apsides(double J, double *pol, int ipli);
+
+/* planets, s. moshplan.c */
+extern int swi_moshplan(double tjd, int ipli, AS_BOOL do_save, double *xpret, double *xeret, char *serr);
+extern int swi_moshplan2(double J, int iplm, double *pobj);
+extern int swi_osc_el_plan(double tjd, double *xp, int ipl, int ipli, double *xearth, double *xsun, char *serr);
+extern FILE *swi_fopen(int ifno, char *fname, char *ephepath, char *serr);
+
+/* nutation */
+struct nut {
+  double tnut;
+  double nutlo[2];	/* nutation in longitude and obliquity */
+  double snut, cnut;	/* sine and cosine of nutation in obliquity */
+  double matrix[3][3];	
+};
+
+struct plantbl {
+  char max_harmonic[9];
+  char max_power_of_t;
+  signed char *arg_tbl;
+  double *lon_tbl;
+  double *lat_tbl;
+  double *rad_tbl;
+  double distance;
+};
+
+struct file_data {
+  char fnam[AS_MAXCH];	/* ephemeris file name */
+  int fversion;		/* version number of file */
+  char astnam[50];	/* asteroid name, if asteroid file */ 
+  int32 sweph_denum;     /* DE number of JPL ephemeris, which this file
+			 * is derived from. */
+  FILE *fptr;		/* ephemeris file pointer */
+  double tfstart;       /* file may be used from this date */
+  double tfend;         /*      through this date          */
+  int32 iflg; 		/* byte reorder flag and little/bigendian flag */
+  short npl;		/* how many planets in file */
+  int ipl[SEI_FILE_NMAXPLAN];	/* planet numbers */
+};
+ 
+struct gen_const {
+ double clight, 
+	aunit, 
+	helgravconst, 
+	ratme, 
+	sunradius;
+};
+
+struct save_positions {
+  int ipl;
+  double tsave;
+  int32 iflgsave;
+  /* position at t = tsave,
+   * in ecliptic polar (offset 0),
+   *    ecliptic cartesian (offset 6), 
+   *    equatorial polar (offset 12),
+   *    and equatorial cartesian coordinates (offset 18).
+   * 6 doubles each for position and speed coordinates.
+   */
+  double xsaves[24];    
+};
+
+struct node_data {
+  /* result of most recent data evaluation for this body: */
+  double teval;		/* time for which last computation was made */
+  int32 iephe;            /* which ephemeris was used */
+  double x[6];		/* position and speed vectors equatorial J2000 */
+  int32 xflgs;		/* hel., light-time, aberr., prec. flags etc. */
+  double xreturn[24];   /* return positions: 
+			 * xreturn+0	ecliptic polar coordinates
+			 * xreturn+6	ecliptic cartesian coordinates
+			 * xreturn+12	equatorial polar coordinates
+			 * xreturn+18	equatorial cartesian coordinates
+			 */
+};
+
+struct topo_data {
+  double geolon, geolat, geoalt;
+  double teval;
+  double tjd_ut;
+  double xobs[6];
+};
+
+struct sid_data {
+  int32 sid_mode;
+  double ayan_t0;
+  double t0;
+};
+
+struct swe_data {
+  AS_BOOL ephe_path_is_set;
+  short jpl_file_is_open;
+  FILE *fixfp;		/* fixed stars file pointer */
+  char ephepath[AS_MAXCH];
+  char jplfnam[AS_MAXCH];
+  short jpldenum;
+  AS_BOOL geopos_is_set;
+  AS_BOOL ayana_is_set;
+  AS_BOOL is_old_starfile;
+  struct file_data fidat[SEI_NEPHFILES];
+  struct gen_const gcdat;
+  struct plan_data pldat[SEI_NPLANETS];
+#if 0
+  struct node_data nddat[SEI_NNODE_ETC];
+#else
+  struct plan_data nddat[SEI_NNODE_ETC];
+#endif
+  struct save_positions savedat[SE_NPLANETS+1];
+  struct epsilon oec;
+  struct epsilon oec2000;
+  struct nut nut;
+  struct nut nut2000;
+  struct nut nutv;
+  struct topo_data topd;
+  struct sid_data sidd;
+  char astelem[AS_MAXCH * 2];
+  double ast_G;
+  double ast_H;
+  double ast_diam;
+  int i_saved_planet_name;
+  char saved_planet_name[80];
+};
+
+extern struct swe_data FAR swed;