gergely/swe-glib
1
0
Fork 0
Code Issues Pull Requests Releases Wiki Activity
243eac9de8
swe-glib/swe/src/swetest.c

2943 lines
99 KiB
C
Raw Blame History

/*
$Header: /users/dieter/sweph/RCS/swetest.c,v 1.78 2010/06/25 07:22:10 dieter Exp $
swetest.c A test program
Authors: Dieter Koch and Alois Treindl, Astrodienst Zuerich
**************************************************************/
/* 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.
*/
/* attention: Microsoft Compiler does not accept strings > 2048 char */
static char *infocmd0 = "\n\
Swetest computes a complete set of geocentric planetary positions,\n\
for a given date or a sequence of dates.\n\
Input can either be a date or an absolute julian day number.\n\
0:00 (midnight).\n\
With the proper options, swetest can be used to output a printed\n\
ephemeris and transfer the data into other programs like spreadsheets\n\
for graphical display.\n\
Version: $Header: /users/dieter/sweph/RCS/swetest.c,v 1.78 2010/06/25 07:22:10 dieter Exp $\n\
\n";
static char *infocmd1 = "\n\
Command line options:\n\
help commands:\n\
-?, -h display whole info\n\
-hcmd display commands\n\
-hplan display planet numbers\n\
-hform display format characters\n\
-hdate display input date format\n\
-hexamp display examples\n\
input time formats:\n\
-bDATE begin date; e.g. -b1.1.1992 if\n\
Note: the date format is day month year (European style).\n\
-bj... begin date as an absolute Julian day number; e.g. -bj2415020.5\n\
-j... same as -bj\n\
-tHH.MMSS input time (ephemeris time)\n\
-ut input date is universal time\n\
-utHH:MM:SS input time\n\
-utHH.MMSS input time\n\
output time for eclipses, occultations, risings/settings is UT by default\n\
-lmt output date/time is LMT (with -geopos)\n\
-lat output date/time is LAT (with -geopos)\n\
object, number of steps, step with\n\
-pSEQ planet sequence to be computed.\n\
See the letter coding below.\n\
-dX differential ephemeris: print differential ephemeris between\n\
body X and each body in list given by -p\n\
example: -p2 -d0 -fJl -n366 -b1.1.1992 prints the longitude\n\
distance between SUN (planet 0) and MERCURY (planet 2)\n\
for a full year starting at 1 Jan 1992.\n\
-DX midpoint ephemeris, works the same way as the differential\n\
mode -d described above, but outputs the midpoint position.\n\
-nN output data for N consecutive days; if no -n option\n\
is given, the default is 1. If the option -n without a\n\
number is given, the default is 20.\n\
-sN timestep N days, default 1. This option is only meaningful\n\
when combined with option -n.\n\
";
static char *infocmd2 = "\
output format:\n\
-fSEQ use SEQ as format sequence for the output columns;\n\
default is PLBRS.\n\
-head don\'t print the header before the planet data. This option\n\
is useful when you want to paste the output into a\n\
spreadsheet for displaying graphical ephemeris.\n\
+head header before every step (with -s..) \n\
-gPPP use PPP as gap between output columns; default is a single\n\
blank. -g followed by white space sets the\n\
gap to the TAB character; which is useful for data entry\n\
into spreadsheets.\n\
astrological house system:\n\
-house[long,lat,hsys] \n\
include house cusps. The longitude, latitude (degrees with\n\
DECIMAL fraction) and house system letter can be given, with\n\
commas separated, + for east and north. If none are given,\n\
Greenwich UK and Placidus is used: 0.00,51.50,p.\n\
The output lists 12 house cusps, Asc, MC, ARMC and Vertex.\n\
Houses can only be computed if option -ut is given.\n\
A equal\n\
E equal\n\
B Alcabitius\n\
C Campanus\n\
G 36 Gauquelin sectors\n\
H horizon / azimut\n\
K Koch\n\
M Morinus\n\
O Porphyry\n\
P Placidus\n\
R Regiomontanus\n\
T Polich/Page (\"topocentric\")\n\
U Krusinski-Pisa-Goelzer\n\
V equal Vehlow\n\
W equal, whole sign\n\
X axial rotation system/ Meridian houses\n\
Y APC houses\n\
-hsy[hsys] \n\
house system to be used (for house positions of planets)\n\
for long, lat, hsys, see -house\n\
";
static char *infocmd3 = "\
-geopos[long,lat,elev] \n\
Geographic position. Can be used for azimuth and altitude\n\
or topocentric or house cups calculations.\n\
The longitude, latitude (degrees with DECIMAL fraction)\n\
and elevation (meters) can be given, with\n\
commas separated, + for east and north. If none are given,\n\
Greenwich is used: 0,51.5,0\n\
sidereal astrology:\n\
-ay.. ayanamsa, with number of method, e.g. ay0 for Fagan/Bradley\n\
-sid.. sidereal, with number of method; 'sid0' for Fagan/Bradley\n\
'sid1' for Lahiri\n\
-sidt0.. sidereal, projection on ecliptic of t0 \n\
-sidsp.. sidereal, projection on solar system plane \n\
ephemeris specifications:\n\
-edirPATH change the directory of the ephemeris files \n\
-eswe swiss ephemeris\n\
-ejpl jpl ephemeris (DE431), or with ephemeris file name\n\
-ejplde200.eph \n\
-emos moshier ephemeris\n\
-true true positions\n\
-noaberr no aberration\n\
-nodefl no gravitational light deflection\n\
-noaberr -nodefl astrometric positions\n\
-j2000 no precession (i.e. J2000 positions)\n\
-icrs ICRS (use Internat. Celestial Reference System)\n\
-nonut no nutation \n\
-speed calculate high precision speed \n\
-speed3 'low' precision speed from 3 positions \n\
do not use this option. -speed parameter\n\
is faster and preciser \n\
-iXX force iflag to value XX\n\
-testaa96 test example in AA 96, B37,\n\
i.e. venus, j2450442.5, DE200.\n\
attention: use precession IAU1976\n\
and nutation 1980 (s. swephlib.h)\n\
-testaa95\n\
-testaa97\n\
-roundsec round to seconds\n\
-roundmin round to minutes\n\
";
static char *infocmd4 = "\
observer position:\n\
-hel compute heliocentric positions\n\
-bary compute barycentric positions (bar. earth instead of node) \n\
-topo[long,lat,elev] \n\
topocentric positions. The longitude, latitude (degrees with\n\
DECIMAL fraction) and elevation (meters) can be given, with\n\
commas separated, + for east and north. If none are given,\n\
Zuerich is used: 8.55,47.38,400\n\
\n\
special events:\n\
-solecl solar eclipse\n\
output 1st line:\n\
eclipse date,\n\
time of maximum (UT),\n\
core shadow width (negative with total eclipses),\n\
fraction of solar diameter that is eclipsed\n\
Julian day number (6-digit fraction) of maximum\n\
output 2nd line:\n\
start and end times for partial and total phase\n\
output 3rd line:\n\
geographical longitude and latitude of maximum eclipse,\n\
totality duration at that geographical position,\n\
output with -local, see below.\n\
-occult occultation of planet or star by the moon. Use -p to \n\
specify planet (-pf -xfAldebaran for stars) \n\
output format same as with -solecl\n\
-lunecl lunar eclipse\n\
output 1st line:\n\
eclipse date,\n\
time of maximum (UT),\n\
Julian day number (6-digit fraction) of maximum\n\
output 2nd line:\n\
6 contacts for start and end of penumbral, partial, and\n\
total phase\n\
";
static char *infocmd5 = "\
-local only with -solecl or -occult, if the next event of this\n\
kind is wanted for a given geogr. position.\n\
Use -geopos[long,lat,elev] to specify that position.\n\
If -local is not set, the program \n\
searches for the next event anywhere on earth.\n\
output 1st line:\n\
eclipse date,\n\
time of maximum,\n\
fraction of solar diameter that is eclipsed\n\
output 2nd line:\n\
local eclipse duration,\n\
local four contacts,\n\
-hev[type] heliacal events,\n\
type 1 = heliacal rising\n\
type 2 = heliacal setting\n\
type 3 = evening first\n\
type 4 = morning last\n\
type 0 or missing = all four events are listed.\n\
-rise rising and setting of a planet or star.\n\
Use -geopos[long,lat,elev] to specify geographical position.\n\
-metr southern and northern meridian transit of a planet of star\n\
Use -geopos[long,lat,elev] to specify geographical position.\n\
specifications for eclipses:\n\
-total total eclipse (only with -solecl, -lunecl)\n\
-partial partial eclipse (only with -solecl, -lunecl)\n\
-annular annular eclipse (only with -solecl)\n\
-anntot annular-total (hybrid) eclipse (only with -solecl)\n\
-penumbral penumbral lunar eclipse (only with -lunecl)\n\
-central central eclipse (only with -solecl, nonlocal)\n\
-noncentral non-central eclipse (only with -solecl, nonlocal)\n\
specifications for risings and settings:\n\
-norefrac neglect refraction (with option -rise)\n\
-disccenter find rise of disc center (with option -rise)\n\
-hindu hindu version of sunrise (with option -rise)\n\
";
static char *infocmd6 = "\
specifications for heliacal events:\n\
-at[press,temp,rhum,visr]:\n\
pressure in hPa\n\
temperature in degrees Celsius\n\
relative humidity in %\n\
visual range, interpreted as follows:\n\
> 1 : meteorological range in km\n\
1>visr>0 : total atmospheric coefficient (ktot)\n\
= 0 : calculated from press, temp, rhum\n\
Default values are -at1013.25,15,40,0\n\
-obs[age,SN] age of observer and Snellen ratio\n\
Default values are -obs36,1\n\
-opt[age,SN,binocular,magn,diam,transm]\n\
age and SN as with -obs\n\
0 monocular or 1 binocular\n\
telescope magnification\n\
optical aperture in mm\n\
optical transmission\n\
Default values: -opt36,1,1,1,0,0 (naked eye)\n\
backward search:\n\
-bwd\n";
/* characters still available:
bcgijklruvxy
*/
static char *infoplan = "\n\
Planet selection letters:\n\
planetary lists:\n\
d (default) main factors 0123456789mtABCcg\n\
p main factors as above, plus main asteroids DEFGHI\n\
h ficticious factors J..X\n\
a all factors\n\
(the letters above can only appear as a single letter)\n\n\
single planet letters:\n\
0 Sun (character zero)\n\
1 Moon (character 1)\n\
2 Mercury\n\
....\n\
9 Pluto\n\
m mean lunar node\n\
t true lunar node\n\
n nutation\n\
o obliquity of ecliptic\n\
q delta t\n\
y time equation\n\
A mean lunar apogee (Lilith, Black Moon) \n\
B osculating lunar apogee \n\
c intp. lunar apogee \n\
g intp. lunar perigee \n\
C Earth (in heliocentric or barycentric calculation)\n\
dwarf planets, plutoids\n\
F Ceres\n\
9 Pluto\n\
s -xs136199 Eris\n\
s -xs136472 Makemake\n\
s -xs136108 Haumea\n\
some minor planets:\n\
D Chiron\n\
E Pholus\n\
G Pallas \n\
H Juno \n\
I Vesta \n\
s minor planet, with MPC number given in -xs\n\
fixed stars:\n\
f fixed star, with name or number given in -xf option\n\
f -xfSirius Sirius\n\
fictitious objects:\n\
J Cupido \n\
K Hades \n\
L Zeus \n\
M Kronos \n\
N Apollon \n\
O Admetos \n\
P Vulkanus \n\
Q Poseidon \n\
R Isis (Sevin) \n\
S Nibiru (Sitchin) \n\
T Harrington \n\
U Leverrier's Neptune\n\
V Adams' Neptune\n\
W Lowell's Pluto\n\
X Pickering's Pluto\n\
Y Vulcan\n\
Z White Moon\n\
w Waldemath's dark Moon\n\
z hypothetical body, with number given in -xz\n\
e print a line of labels\n\
\n";
/* characters still available
cgjv
*/
static char *infoform = "\n\
Output format SEQ letters:\n\
In the standard setting five columns of coordinates are printed with\n\
the default format PLBRS. You can change the default by providing an\n\
option like -fCCCC where CCCC is your sequence of columns.\n\
The coding of the sequence is like this:\n\
y year\n\
Y year.fraction_of_year\n\
p planet index\n\
P planet name\n\
J absolute juldate\n\
T date formatted like 23.02.1992 \n\
t date formatted like 920223 for 1992 february 23\n\
L longitude in degree ddd mm'ss\"\n\
l longitude decimal\n\
Z longitude ddsignmm'ss\"\n\
S speed in longitude in degree ddd:mm:ss per day\n\
SS speed for all values specified in fmt\n\
s speed longitude decimal (degrees/day)\n\
ss speed for all values specified in fmt\n\
B latitude degree\n\
b latitude decimal\n\
R distance decimal in AU\n\
r distance decimal in AU, Moon in seconds parallax\n\
relative distance (1000=nearest, 0=furthest)\n\
A right ascension in hh:mm:ss\n\
a right ascension hours decimal\n\
D declination degree\n\
d declination decimal\n\
I Azimuth degree\n\
i Azimuth decimal\n\
H Height degree\n\
h Height decimal\n\
K Height (with refraction) degree\n\
k Height (with refraction) decimal\n\
G house position in degrees\n\
g house position in degrees decimal\n\
j house number 1.0 - 12.99999\n\
X x-, y-, and z-coordinates ecliptical\n\
x x-, y-, and z-coordinates equatorial\n\
U unit vector ecliptical\n\
u unit vector equatorial\n\
Q l, b, r, dl, db, dr, a, d, da, dd\n\
n nodes (mean): ascending/descending (Me - Ne); longitude decimal\n\
N nodes (osculating): ascending/descending, longitude; decimal\n\
f apsides (mean): perihel, aphel, second focal point; longitude dec.\n\
F apsides (osc.): perihel, aphel, second focal point; longitude dec.\n\
+ phase angle\n\
- phase\n\
* elongation\n\
/ apparent diameter of disc (without refraction)\n\
= magnitude\n";
static char *infoform2 = "\
v (reserved)\n\
V (reserved)\n\
";
static char *infodate = "\n\
Date entry:\n\
In the interactive mode, when you are asked for a start date,\n\
you can enter data in one of the following formats:\n\
\n\
1.2.1991 three integers separated by a nondigit character for\n\
day month year. Dates are interpreted as Gregorian\n\
after 4.10.1582 and as Julian Calendar before.\n\
Time is always set to midnight.\n\
If the three letters jul are appended to the date,\n\
the Julian calendar is used even after 1582.\n\
If the four letters greg are appended to the date,\n\
the Gregorian calendar is used even before 1582.\n\
\n\
j2400123.67 the letter j followed by a real number, for\n\
the absolute Julian daynumber of the start date.\n\
Fraction .5 indicates midnight, fraction .0\n\
indicates noon, other times of the day can be\n\
chosen accordingly.\n\
\n\
<RETURN> repeat the last entry\n\
\n\
. stop the program\n\
\n\
+20 advance the date by 20 days\n\
\n\
-10 go back in time 10 days\n";
static char *infoexamp = "\n\
\n\
Examples:\n\
\n\
swetest -p2 -b1.12.1900 -n15 -s2\n\
ephemeris of Mercury (-p2) starting on 1 Dec 1900,\n\
15 positions (-n15) in two-day steps (-s2)\n\
\n\
swetest -p2 -b1.12.1900 -n15 -s2 -fTZ -roundsec -g, -head\n\
same, but output format = date and zodiacal position (-fTZ),\n\
separated by comma (-g,) and rounded to seconds (-roundsec),\n\
without header (-head).\n\
\n\
swetest -ps -xs433 -b1.12.1900\n\
position of asteroid 433 Eros (-ps -xs433)\n\
\n\
swetest -pf -xfAldebaran -b1.1.2000\n\
position of fixed star Aldebaran \n\
\n\
swetest -p1 -d0 -b1.12.1900 -n10 -fPTl -head\n\
angular distance of moon (-p1) from sun (-d0) for 10\n\
consecutive days (-n10).\n\
\n\
swetest -p6 -DD -b1.12.1900 -n100 -s5 -fPTZ -head -roundmin\n\
Midpoints between Saturn (-p6) and Chiron (-DD) for 100\n\
consecutive steps (-n100) with 5-day steps (-s5) with\n\
longitude in degree-sign format (-f..Z) rounded to minutes (-roundmin)\n\
\n\
swetest -b5.1.2002 -p -house12.05,49.50,k -ut12:30\n\
Koch houses for a location in Germany at a given date and time\n";
/**************************************************************/
#include "swephexp.h" /* this includes "sweodef.h" */
/*
* programmers warning: It looks much worse than it is!
* Originally swetest.c was a small and simple test program to test
* the main functions of the Swiss Ephemeris and to demonstrate
* its precision.
* It compiles on Unix, on MSDOS and as a non-GUI utility on 16-bit
* and 32-bit windows.
* This portability has forced us into some clumsy constructs, which
* end to hide the actual simplicity of the use of Swiss Ephemeris.
* For example, the mechanism implemented here in swetest.c to find
* the binary ephemeris files overrides the much simpler mechanism
* inside the SwissEph library. This was necessary because we wanted
* swetest.exe to run directly off the CDROM and search with some
* intelligence for ephemeris files already installed on a system.
*/
#if MSDOS
# include <direct.h>
# include <dos.h>
# ifdef _MSC_VER
# include <sys\types.h>
# endif
# include <sys\stat.h>
# include <float.h>
#else
# ifdef MACOS
# include <console.h>
# else
# include <sys/stat.h>
# endif
#endif
#define J2000 2451545.0 /* 2000 January 1.5 */
#define square_sum(x) (x[0]*x[0]+x[1]*x[1]+x[2]*x[2])
#define SEFLG_EPHMASK (SEFLG_JPLEPH|SEFLG_SWIEPH|SEFLG_MOSEPH)
#define BIT_ROUND_SEC 1
#define BIT_ROUND_MIN 2
#define BIT_ZODIAC 4
#define BIT_LZEROES 8
#define BIT_TIME_LZEROES 8
#define BIT_TIME_LMT 16
#define BIT_TIME_LAT 32
#define PLSEL_D "0123456789mtA"
#define PLSEL_P "0123456789mtABCcgDEFGHI"
#define PLSEL_H "JKLMNOPQRSTUVWXYZw"
#define PLSEL_A "0123456789mtABCcgDEFGHIJKLMNOPQRSTUVWXYZw"
#define DIFF_DIFF 'd'
#define DIFF_MIDP 'D'
#define MODE_HOUSE 1
#define MODE_LABEL 2
static char se_pname[AS_MAXCH];
static char *zod_nam[] = {"ar", "ta", "ge", "cn", "le", "vi",
"li", "sc", "sa", "cp", "aq", "pi"};
static char star[AS_MAXCH] = "algol", star2[AS_MAXCH];
static char sastno[AS_MAXCH] = "433";
static char shyp[AS_MAXCH] = "1";
static char *dms(double x, int32 iflag);
static int make_ephemeris_path(int32 iflag, char *argv0, char *ephepath);
static int letter_to_ipl(int letter);
static int print_line(int mode);
static int do_special_event(double tjd, int32 ipl, char *star, int32 special_event, int32 special_mode, double *geopos, double *datm, double *dobs, char *serr) ;
static char *hms_from_tjd(double x);
static void do_printf(char *info);
static char *hms(double x, int32 iflag);
static void remove_whitespace(char *s);
#if MSDOS
static int cut_str_any(char *s, char *cutlist, char *cpos[], int nmax);
#endif
/* globals shared between main() and print_line() */
static char *fmt = "PLBRS";
static char *gap = " ";
static double t, te, tut, jut = 0;
static int jmon, jday, jyear;
static int ipl = SE_SUN, ipldiff = SE_SUN, nhouses = 12;
static char spnam[AS_MAXCH], spnam2[AS_MAXCH], serr[AS_MAXCH];
static char serr_save[AS_MAXCH], serr_warn[AS_MAXCH];
static int gregflag = SE_GREG_CAL;
static int diff_mode = 0;
static AS_BOOL universal_time = FALSE;
static int32 round_flag = 0;
static int32 time_flag = 0;
static AS_BOOL short_output = FALSE;
static int32 special_event = 0;
static int32 special_mode = 0;
static AS_BOOL hel_using_AV = FALSE;
static double x[6], x2[6], xequ[6], xcart[6], xcartq[6], xobl[6], xaz[6], xt[6], hpos, hpos2, hposj, armc, xsv[6];
static int hpos_meth = 0;
static double geopos[10];
static double attr[20], tret[20], datm[4], dobs[6];
static int32 iflag = 0, iflag2; /* external flag: helio, geo... */
static char *hs_nam[] = {"undef",
"Ascendant", "MC", "ARMC", "Vertex"};
static int direction = 1;
static AS_BOOL direction_flag = FALSE;
static int32 helflag = 0;
static double tjd = 2415020.5;
static int32 nstep = 1, istep;
static int32 search_flag = 0;
static char sout[AS_MAXCH];
static int32 whicheph = SEFLG_SWIEPH;
static char *psp;
static int32 norefrac = 0;
static int32 disccenter = 0;
#define SP_LUNAR_ECLIPSE 1
#define SP_SOLAR_ECLIPSE 2
#define SP_OCCULTATION 3
#define SP_RISE_SET 4
#define SP_MERIDIAN_TRANSIT 5
#define SP_HELIACAL 6
# define SP_MODE_HOW 2 /* an option for Lunar */
# define SP_MODE_LOCAL 8 /* an option for Solar */
# define SP_MODE_HOCAL 4096
# define ECL_LUN_PENUMBRAL 1 /* eclipse types for hocal list */
# define ECL_LUN_PARTIAL 2
# define ECL_LUN_TOTAL 3
# define ECL_SOL_PARTIAL 4
# define ECL_SOL_ANNULAR 5
# define ECL_SOL_TOTAL 6
int main(int argc, char *argv[])
{
char sdate_save[AS_MAXCH];
char s1[AS_MAXCH], s2[AS_MAXCH];
char *sp, *sp2;
char *spno;
char *plsel = PLSEL_D;
#if HPUNIX
char hostname[80];
#endif
int i, j, n, iflag_f = -1, iflgt;
int line_count, line_limit = 32000;
double daya;
double top_long = 0.0; /* Greenwich UK */
double top_lat = 51.5;
double top_elev = 0;
AS_BOOL have_geopos = FALSE;
int ihsy = 'p';
AS_BOOL do_houses = FALSE;
char ephepath[AS_MAXCH];
char fname[AS_MAXCH];
char sdate[AS_MAXCH];
char *begindate = NULL;
int32 iflgret;
AS_BOOL with_header = TRUE;
AS_BOOL with_header_always = FALSE;
AS_BOOL do_ayanamsa = FALSE;
int32 sid_mode = SE_SIDM_FAGAN_BRADLEY;
double t2, tstep = 1, thour = 0;
double delt;
datm[0] = 0; datm[1] = 0; datm[2] = 0; datm[3] = 0;
dobs[0] = 0; dobs[1] = 0;
dobs[2] = 0; dobs[3] = 0; dobs[4] = 0; dobs[5] = 0;
/* swe_set_tid_acc(-25.858); * to test delta t output */
serr[0] = serr_save[0] = serr_warn[0] = sdate_save[0] = '\0';
# ifdef MACOS
argc = ccommand(&argv); /* display the arguments window */
# endif
strcpy(ephepath, "");
strcpy(fname, SE_FNAME_DFT);
for (i = 1; i < argc; i++) {
if (strncmp(argv[i], "-ut", 3) == 0) {
universal_time = TRUE;
if (strlen(argv[i]) > 3) {
*s1 = '\0';
strncat(s1, argv[i] + 3, 30);
if ((sp = strchr(s1, ':')) != NULL) {
*sp = '.';
if ((sp = strchr(s1, ':')) != NULL) {
strcpy(s2, sp + 1);
strcpy(sp, s2);
}
}
thour = atof(s1);
/* h.mmss -> decimal */
t = fmod(thour, 1) * 100 + 1e-6;
j = (int) t;
t = fmod(t, 1) * 100 + 1e-6;
thour = (int) thour + j / 60.0 + t / 3600.0;
}
} else if (strncmp(argv[i], "-head", 5) == 0) {
with_header = FALSE;
} else if (strncmp(argv[i], "+head", 5) == 0) {
with_header_always = TRUE;
} else if (strcmp(argv[i], "-j2000") == 0) {
iflag |= SEFLG_J2000;
} else if (strcmp(argv[i], "-icrs") == 0) {
iflag |= SEFLG_ICRS;
} else if (strncmp(argv[i], "-ay", 3) == 0) {
do_ayanamsa = TRUE;
sid_mode = atol(argv[i]+3);
swe_set_sid_mode(sid_mode, 0, 0);
} else if (strncmp(argv[i], "-sidt0", 6) == 0) {
iflag |= SEFLG_SIDEREAL;
sid_mode = atol(argv[i]+6);
if (sid_mode == 0)
sid_mode = SE_SIDM_FAGAN_BRADLEY;
sid_mode |= SE_SIDBIT_ECL_T0;
swe_set_sid_mode(sid_mode, 0, 0);
} else if (strncmp(argv[i], "-sidsp", 6) == 0) {
iflag |= SEFLG_SIDEREAL;
sid_mode = atol(argv[i]+6);
if (sid_mode == 0)
sid_mode = SE_SIDM_FAGAN_BRADLEY;
sid_mode |= SE_SIDBIT_SSY_PLANE;
swe_set_sid_mode(sid_mode, 0, 0);
} else if (strncmp(argv[i], "-sid", 4) == 0) {
iflag |= SEFLG_SIDEREAL;
sid_mode = atol(argv[i]+4);
if (sid_mode > 0)
swe_set_sid_mode(sid_mode, 0, 0);
} else if (strcmp(argv[i], "-jplhora") == 0) {
iflag |= SEFLG_JPLHOR_APPROX;
} else if (strcmp(argv[i], "-jplhor") == 0) {
iflag |= SEFLG_JPLHOR;
} else if (strncmp(argv[i], "-j", 2) == 0) {
begindate = argv[i] + 1;
} else if (strncmp(argv[i], "-ejpl", 5) == 0) {
whicheph = SEFLG_JPLEPH;
if (*(argv[i]+5) != '\0') {
*fname = '\0';
strncat(fname, argv[i]+5, sizeof(fname) - 1);
}
} else if (strncmp(argv[i], "-edir", 5) == 0) {
if (*(argv[i]+5) != '\0') {
*ephepath = '\0';
strncat(ephepath, argv[i]+5, sizeof(ephepath) - 1);
}
} else if (strcmp(argv[i], "-eswe") == 0) {
whicheph = SEFLG_SWIEPH;
} else if (strcmp(argv[i], "-emos") == 0) {
whicheph = SEFLG_MOSEPH;
} else if (strncmp(argv[i], "-helflag", 8) == 0) {
helflag = atoi(argv[i]+8);
if (helflag >= SE_HELFLAG_AV)
hel_using_AV = TRUE;
} else if (strcmp(argv[i], "-hel") == 0) {
iflag |= SEFLG_HELCTR;
} else if (strcmp(argv[i], "-bary") == 0) {
iflag |= SEFLG_BARYCTR;
} else if (strncmp(argv[i], "-house", 6) == 0) {
sout[0] = '\0';
sout[1] = '\0';
sscanf(argv[i] + 6, "%lf,%lf,%c", &top_long, &top_lat, sout);
top_elev = 0;
if (*sout) ihsy = sout[0];
do_houses = TRUE;
have_geopos = TRUE;
} else if (strncmp(argv[i], "-hsy", 4) == 0) {
ihsy = *(argv[i] + 4);
if (ihsy == '\0') ihsy = 'p';
if (strlen(argv[i]) > 5)
hpos_meth = atoi(argv[i] + 5);
have_geopos = TRUE;
} else if (strncmp(argv[i], "-topo", 5) == 0) {
iflag |= SEFLG_TOPOCTR;
sscanf(argv[i] + 5, "%lf,%lf,%lf", &top_long, &top_lat, &top_elev);
have_geopos = TRUE;
} else if (strncmp(argv[i], "-geopos", 7) == 0) {
sscanf(argv[i] + 7, "%lf,%lf,%lf", &top_long, &top_lat, &top_elev);
have_geopos = TRUE;
} else if (strcmp(argv[i], "-true") == 0) {
iflag |= SEFLG_TRUEPOS;
} else if (strcmp(argv[i], "-noaberr") == 0) {
iflag |= SEFLG_NOABERR;
} else if (strcmp(argv[i], "-nodefl") == 0) {
iflag |= SEFLG_NOGDEFL;
} else if (strcmp(argv[i], "-nonut") == 0) {
iflag |= SEFLG_NONUT;
} else if (strcmp(argv[i], "-speed3") == 0) {
iflag |= SEFLG_SPEED3;
} else if (strcmp(argv[i], "-speed") == 0) {
iflag |= SEFLG_SPEED;
} else if (strncmp(argv[i], "-testaa", 7) == 0) {
whicheph = SEFLG_JPLEPH;
strcpy(fname, SE_FNAME_DE200);
if (strcmp(argv[i]+7, "95") == 0)
begindate = "j2449975.5";
if (strcmp(argv[i]+7, "96") == 0)
begindate = "j2450442.5";
if (strcmp(argv[i]+7, "97") == 0)
begindate = "j2450482.5";
fmt = "PADRu";
universal_time = FALSE;
plsel="3";
} else if (strcmp(argv[i], "-lmt") == 0) {
universal_time = TRUE;
time_flag |= BIT_TIME_LMT;
} else if (strcmp(argv[i], "-lat") == 0) {
universal_time = TRUE;
time_flag |= BIT_TIME_LAT;
} else if (strcmp(argv[i], "-lunecl") == 0) {
special_event = SP_LUNAR_ECLIPSE;
} else if (strcmp(argv[i], "-solecl") == 0) {
special_event = SP_SOLAR_ECLIPSE;
have_geopos = TRUE;
} else if (strcmp(argv[i], "-short") == 0) {
short_output = TRUE;
} else if (strcmp(argv[i], "-occult") == 0) {
special_event = SP_OCCULTATION;
have_geopos = TRUE;
} else if (strcmp(argv[i], "-hocal") == 0) {
/* used to create a listing for inclusion in hocal.c source code */
special_mode |= SP_MODE_HOCAL;
} else if (strcmp(argv[i], "-how") == 0) {
special_mode |= SP_MODE_HOW;
} else if (strcmp(argv[i], "-total") == 0) {
search_flag |= SE_ECL_TOTAL|SE_ECL_CENTRAL|SE_ECL_NONCENTRAL;
} else if (strcmp(argv[i], "-annular") == 0) {
search_flag |= SE_ECL_ANNULAR|SE_ECL_CENTRAL|SE_ECL_NONCENTRAL;
} else if (strcmp(argv[i], "-anntot") == 0) {
search_flag |= SE_ECL_ANNULAR_TOTAL|SE_ECL_CENTRAL|SE_ECL_NONCENTRAL;
} else if (strcmp(argv[i], "-partial") == 0) {
search_flag |= SE_ECL_PARTIAL|SE_ECL_CENTRAL|SE_ECL_NONCENTRAL;
} else if (strcmp(argv[i], "-penumbral") == 0) {
search_flag |= SE_ECL_PENUMBRAL;
} else if (strcmp(argv[i], "-noncentral") == 0) {
search_flag &= ~SE_ECL_CENTRAL;
search_flag |= SE_ECL_NONCENTRAL;
} else if (strcmp(argv[i], "-central") == 0) {
search_flag &= ~SE_ECL_NONCENTRAL;
search_flag |= SE_ECL_CENTRAL;
} else if (strcmp(argv[i], "-local") == 0) {
special_mode |= SP_MODE_LOCAL;
} else if (strcmp(argv[i], "-rise") == 0) {
special_event = SP_RISE_SET;
have_geopos = TRUE;
} else if (strcmp(argv[i], "-norefrac") == 0) {
norefrac = 1;
} else if (strcmp(argv[i], "-disccenter") == 0) {
disccenter = 1;
} else if (strcmp(argv[i], "-hindu") == 0) {
norefrac = 1;
disccenter = 1;
} else if (strcmp(argv[i], "-metr") == 0) {
special_event = SP_MERIDIAN_TRANSIT;
have_geopos = TRUE;
} else if (strncmp(argv[i], "-hev", 4) == 0) {
special_event = SP_HELIACAL;
search_flag = 0;
if (strlen(argv[i]) > 4)
search_flag = atoi(argv[i] + 4);
have_geopos = TRUE;
if (strstr(argv[i], "AV")) hel_using_AV = TRUE;
} else if (strncmp(argv[i], "-at", 3) == 0) {
sscanf(argv[i]+3, "%lf,%lf,%lf,%lf", &(datm[0]), &(datm[1]), &(datm[2]), &(datm[3]));
} else if (strncmp(argv[i], "-obs", 4) == 0) {
sscanf(argv[i]+4, "%lf,%lf", &(dobs[0]), &(dobs[1]));
} else if (strncmp(argv[i], "-opt", 4) == 0) {
sscanf(argv[i]+4, "%lf,%lf,%lf,%lf,%lf,%lf", &(dobs[0]), &(dobs[1]), &(dobs[2]), &(dobs[3]), &(dobs[4]), &(dobs[5]));
} else if (strcmp(argv[i], "-bwd") == 0) {
direction = -1;
direction_flag = TRUE;
} else if (strncmp(argv[i], "-p", 2) == 0) {
spno = argv[i]+2;
switch (*spno) {
case 'd':
/*
case '\0':
case ' ':
*/
plsel = PLSEL_D; break;
case 'p': plsel = PLSEL_P; break;
case 'h': plsel = PLSEL_H; break;
case 'a': plsel = PLSEL_A; break;
default: plsel = spno;
}
} else if (strncmp(argv[i], "-xs", 3) == 0) {
/* number of asteroid */
*sastno = '\0';
strncat(sastno, argv[i] + 3, sizeof(sastno) - 1);
} else if (strncmp(argv[i], "-xf", 3) == 0) {
/* name or number of fixed star */
*star = '\0';
strncat(star, argv[i] + 3, sizeof(star) - 1);
} else if (strncmp(argv[i], "-xz", 3) == 0) {
/* number of hypothetical body */
*shyp = '\0';
strncat(shyp, argv[i] + 3, sizeof(shyp) - 1);
} else if (strncmp(argv[i], "-x", 2) == 0) {
/* name or number of fixed star */
*star = '\0';
strncat(star, argv[i] + 2, sizeof(star) - 1);
} else if (strncmp(argv[i], "-n", 2) == 0) {
nstep = atoi(argv[i]+2);
if (nstep == 0)
nstep = 20;
} else if (strncmp(argv[i], "-i", 2) == 0) {
iflag_f = atoi(argv[i]+2);
if (iflag_f & SEFLG_XYZ)
fmt = "PX";
} else if (strncmp(argv[i], "-s", 2) == 0) {
tstep = atof(argv[i]+2);
} else if (strncmp(argv[i], "-b", 2) == 0) {
begindate = argv[i] + 2;
} else if (strncmp(argv[i], "-f", 2) == 0) {
fmt = argv[i] + 2;
} else if (strncmp(argv[i], "-g", 2) == 0) {
gap = argv[i] + 2;
if (*gap == '\0') gap = "\t";
} else if (strncmp(argv[i], "-d", 2) == 0
|| strncmp(argv[i], "-D", 2) == 0) {
diff_mode = *(argv[i] + 1); /* 'd' or 'D' */
sp = argv[i]+2;
ipldiff = letter_to_ipl((int) *sp);
if (ipldiff <0) ipldiff = SE_SUN;
swe_get_planet_name(ipldiff, spnam2);
} else if (strcmp(argv[i], "-roundsec") == 0) {
round_flag |= BIT_ROUND_SEC;
} else if (strcmp(argv[i], "-roundmin") == 0) {
round_flag |= BIT_ROUND_MIN;
} else if (strncmp(argv[i], "-t", 2) == 0) {
if (strlen(argv[i]) > 2) {
*s1 = '\0';
strncat(s1, argv[i] + 2, 30);
if ((sp = strchr(s1, ':')) != NULL) {
*sp = '.';
if ((sp = strchr(s1, ':')) != NULL) {
strcpy(s2, sp + 1);
strcpy(sp, s2);
}
}
thour = atof(s1);
/* h.mmss -> decimal */
t = fmod(thour, 1) * 100 + 1e-6;
j = (int) t;
t = fmod(t, 1) * 100 + 1e-6;
thour = (int) thour + j / 60.0 + t / 3600.0;
}
} else if (strncmp(argv[i], "-h", 2) == 0
|| strncmp(argv[i], "-?", 2) == 0) {
sp = argv[i]+2;
if (*sp == 'c' || *sp == '\0') {
fputs(infocmd0,stdout);
fputs(infocmd1,stdout);
fputs(infocmd2,stdout);
fputs(infocmd3,stdout);
fputs(infocmd4,stdout);
fputs(infocmd5,stdout);
fputs(infocmd6,stdout);
}
if (*sp == 'p' || *sp == '\0')
fputs(infoplan,stdout);
if (*sp == 'f' || *sp == '\0') {
fputs(infoform,stdout);
fputs(infoform2,stdout);
}
if (*sp == 'd' || *sp == '\0')
fputs(infodate,stdout);
if (*sp == 'e' || *sp == '\0')
fputs(infoexamp,stdout);
goto end_main;
} else {
strcpy(sout, "illegal option ");
strncat(sout, argv[i], sizeof(sout) - 50);
strcat(sout, "\n");
fputs(sout,stdout);
exit(1);
}
}
if (special_event == SP_OCCULTATION ||
special_event == SP_RISE_SET ||
special_event == SP_MERIDIAN_TRANSIT ||
special_event == SP_HELIACAL
) {
ipl = letter_to_ipl(*plsel);
if (*plsel == 'f')
ipl = SE_FIXSTAR;
else
*star = '\0';
if (special_event == SP_OCCULTATION && ipl == 1)
ipl = 2; /* no occultation of moon by moon */
}
geopos[0] = top_long;
geopos[1] = top_lat;
geopos[2] = top_elev;
swe_set_topo(top_long, top_lat, top_elev);
#if HPUNIX
gethostname (hostname, 80);
if (strstr(hostname, "as10") != NULL)
line_limit = 1000;
#endif
if (with_header) {
for (i = 0; i < argc; i++) {
fputs(argv[i],stdout);
printf(" ");
}
}
iflag = (iflag & ~SEFLG_EPHMASK) | whicheph;
if (strpbrk(fmt, "SsQ") != NULL)
iflag |= SEFLG_SPEED;
if (*ephepath == '\0') {
if (make_ephemeris_path(iflag, argv[0], ephepath) == ERR) {
iflag = (iflag & ~SEFLG_EPHMASK) | SEFLG_MOSEPH;
whicheph = SEFLG_MOSEPH;
}
}
swe_set_ephe_path(ephepath);
if (whicheph & SEFLG_JPLEPH)
swe_set_jpl_file(fname);
while (TRUE) {
serr[0] = serr_save[0] = serr_warn[0] = '\0';
if (begindate == NULL) {
printf("\nDate ?");
sdate[0] = '\0';
if( !fgets(sdate, AS_MAXCH, stdin) ) goto end_main;
} else {
*sdate = '\0';
strncat(sdate, begindate, AS_MAXCH-1);
begindate = "."; /* to exit afterwards */
}
if (strcmp(sdate, "-bary") == 0) {
iflag = iflag & ~SEFLG_HELCTR;
iflag |= SEFLG_BARYCTR;
*sdate = '\0';
} else if (strcmp(sdate, "-hel") == 0) {
iflag = iflag & ~SEFLG_BARYCTR;
iflag |= SEFLG_HELCTR;
*sdate = '\0';
} else if (strcmp(sdate, "-geo") == 0) {
iflag = iflag & ~SEFLG_BARYCTR;
iflag = iflag & ~SEFLG_HELCTR;
*sdate = '\0';
} else if (strcmp(sdate, "-ejpl") == 0) {
iflag &= ~SEFLG_EPHMASK;
iflag |= SEFLG_JPLEPH;
*sdate = '\0';
} else if (strcmp(sdate, "-eswe") == 0) {
iflag &= ~SEFLG_EPHMASK;
iflag |= SEFLG_SWIEPH;
*sdate = '\0';
} else if (strcmp(sdate, "-emos") == 0) {
iflag &= ~SEFLG_EPHMASK;
iflag |= SEFLG_MOSEPH;
*sdate = '\0';
} else if (strncmp(sdate, "-xs",3) == 0) {
/* number of asteroid */
strcpy(sastno, sdate + 3);
*sdate = '\0';
}
// swe_set_tid_acc((double) (iflag & SEFLG_EPHMASK));
sp = sdate;
if (*sp == '.') {
goto end_main;
} else if (*sp == '\0' || *sp == '\n' || *sp == '\r') {
strcpy (sdate, sdate_save);
} else {
strcpy (sdate_save, sdate);
}
if (*sdate == '\0') {
sprintf(sdate, "j%f", tjd);
}
if (*sp == 'j') { /* it's a day number */
if ((sp2 = strchr(sp, ',')) != NULL)
*sp2 = '.';
sscanf(sp+1,"%lf", &tjd);
if (tjd < 2299160.5)
gregflag = SE_JUL_CAL;
else
gregflag = SE_GREG_CAL;
if (strstr(sp, "jul") != NULL)
gregflag = SE_JUL_CAL;
else if (strstr(sp, "greg") != NULL)
gregflag = SE_GREG_CAL;
swe_revjul(tjd, gregflag, &jyear, &jmon, &jday, &jut);
} else if (*sp == '+') {
n = atoi(sp);
if (n == 0) n = 1;
tjd += n;
swe_revjul(tjd, gregflag, &jyear, &jmon, &jday, &jut);
} else if (*sp == '-') {
n = atoi(sp);
if (n == 0) n = -1;
tjd += n;
swe_revjul(tjd, gregflag, &jyear, &jmon, &jday, &jut);
} else {
if (sscanf (sp, "%d%*c%d%*c%d", &jday,&jmon,&jyear) < 1) exit(1);
if ((int32) jyear * 10000L + (int32) jmon * 100L + (int32) jday < 15821015L)
gregflag = SE_JUL_CAL;
else
gregflag = SE_GREG_CAL;
if (strstr(sp, "jul") != NULL)
gregflag = SE_JUL_CAL;
else if (strstr(sp, "greg") != NULL)
gregflag = SE_GREG_CAL;
jut = 0;
tjd = swe_julday(jyear,jmon,jday,jut,gregflag);
tjd += thour / 24.0;
}
if (special_event > 0) {
do_special_event(tjd, ipl, star, special_event, special_mode, geopos, datm, dobs, serr) ;
swe_close();
return OK;
}
line_count = 0;
for (t = tjd, istep = 1; istep <= nstep; t += tstep, istep++) {
if (t < 2299160.5)
gregflag = SE_JUL_CAL;
else
gregflag = SE_GREG_CAL;
if (strstr(sdate, "jul") != NULL)
gregflag = SE_JUL_CAL;
else if (strstr(sdate, "greg") != NULL)
gregflag = SE_GREG_CAL;
t2 = t;
swe_revjul(t2, gregflag, &jyear, &jmon, &jday, &jut);
if (with_header) {
printf("\ndate (dmy) %d.%d.%d", jday, jmon, jyear);
if (gregflag)
printf(" greg.");
else
printf(" jul.");
t2 = jut;
printf(" % 2d:", (int) t2);
t2 = (t2 - (int32) t2) * 60;
printf("%02d:", (int) t2);
t2 = (t2 - (int32) t2) * 60;
printf("%02d", (int) t2);
if (universal_time)
printf(" UT");
else
printf(" ET");
printf("\t\tversion %s", swe_version(sout));
}
delt = swe_deltat(t);
if (universal_time) {
if (with_header) {
printf("\nUT: %.11f", t);
}
if (with_header) {
printf(" delta t: %f sec", delt * 86400.0);
}
te = t + delt;
tut = t;
} else {
te = t;
tut = t - delt;
}
iflgret = swe_calc(te, SE_ECL_NUT, iflag, xobl, serr);
if (with_header) {
printf("\nET: %.11f", te);
if (iflag & SEFLG_SIDEREAL) {
daya = swe_get_ayanamsa(te);
printf(" ayanamsa = %s", dms(daya, round_flag));
}
if (have_geopos) {
printf("\ngeo. long %f, lat %f, alt %f", geopos[0], geopos[1], geopos[2]);
}
if (iflag_f >=0)
iflag = iflag_f;
if (strchr(plsel, 'o') == NULL) {
printf("\n%-15s %s", "Epsilon (true)", dms(xobl[0],round_flag));
}
if (strchr(plsel, 'n') == NULL) {
fputs("\nNutation ", stdout);
fputs(dms(xobl[2], round_flag), stdout);
fputs(gap, stdout);
fputs(dms(xobl[3], round_flag), stdout);
}
printf("\n");
if (do_houses) {
char *shsy = swe_house_name(ihsy);
if (!universal_time) {
do_houses = FALSE;
printf("option -house requires option -ut for Universal Time\n");
} else {
strcpy(s1, dms(top_long, round_flag));
strcpy(s2, dms(top_lat, round_flag));
printf("Houses system %c (%s) for long=%s, lat=%s\n", ihsy, shsy, s1, s2);
}
}
}
if (with_header && !with_header_always)
with_header = FALSE;
if (do_ayanamsa) {
daya = swe_get_ayanamsa(te);
fputs("Ayanamsa", stdout);
fputs(gap, stdout);
fputs(dms(daya, round_flag), stdout);
fputs("\n", stdout);
/*printf("Ayanamsa%s%s\n", gap, dms(daya, round_flag));*/
continue;
}
if (strchr(plsel, 'e'))
print_line(MODE_LABEL);
for (psp = plsel; *psp != '\0'; psp++) {
if (*psp == 'e') continue;
ipl = letter_to_ipl((int) *psp);
if (ipl == -2) {
printf("illegal parameter -p%s\n", plsel);
exit(1);
}
if (*psp == 'f')
ipl = SE_FIXSTAR;
else if (*psp == 's')
ipl = atoi(sastno) + 10000;
else if (*psp == 'z')
ipl = atoi(shyp) + SE_FICT_OFFSET_1;
if (iflag & SEFLG_HELCTR) {
if (ipl == SE_SUN
|| ipl == SE_MEAN_NODE || ipl == SE_TRUE_NODE
|| ipl == SE_MEAN_APOG || ipl == SE_OSCU_APOG)
continue;
} else if (iflag & SEFLG_BARYCTR) {
if (ipl == SE_MEAN_NODE || ipl == SE_TRUE_NODE
|| ipl == SE_MEAN_APOG || ipl == SE_OSCU_APOG)
continue;
} else /* geocentric */
if (ipl == SE_EARTH)
continue;
/* ecliptic position */
if (iflag_f >=0)
iflag = iflag_f;
if (ipl == SE_FIXSTAR) {
iflgret = swe_fixstar(star, te, iflag, x, serr);
/* magnitude, etc. */
if (iflgret != ERR && strpbrk(fmt, "=") != NULL) {
double mag;
iflgret = swe_fixstar_mag(star, &mag, serr);
attr[4] = mag;
}
strcpy(se_pname, star);
} else {
iflgret = swe_calc(te, ipl, iflag, x, serr);
/* phase, magnitude, etc. */
if (iflgret != ERR && strpbrk(fmt, "+-*/=") != NULL)
iflgret = swe_pheno(te, ipl, iflag, attr, serr);
swe_get_planet_name(ipl, se_pname);
}
if (*psp == 'q') {/* delta t */
x[0] = swe_deltat(te) * 86400;
x[1] = x[2] = x[3] = 0;
strcpy(se_pname, "Delta T");
}
if (*psp == 'o') {/* ecliptic is wanted, remove nutation */
x[2] = x[3] = 0;
strcpy(se_pname, "Ecl. Obl.");
}
if (*psp == 'n') {/* nutation is wanted, remove ecliptic */
x[0] = x[2];
x[1] = x[3];
x[2] = x[3] = 0;
strcpy(se_pname, "Nutation");
}
if (*psp == 'y') {/* time equation */
iflgret = swe_time_equ(tut, &(x[0]), serr);
x[0] *= 86400; /* in seconds */;
x[1] = x[2] = x[3] = 0;
strcpy(se_pname, "Time Equ.");
}
if (iflgret < 0) {
if (strcmp(serr, serr_save) != 0
&& (ipl == SE_SUN || ipl == SE_MOON
|| ipl == SE_MEAN_NODE || ipl == SE_TRUE_NODE
|| ipl == SE_CHIRON || ipl == SE_PHOLUS || ipl == SE_CUPIDO
|| ipl >= SE_AST_OFFSET || ipl == SE_FIXSTAR)) {
fputs("error: ", stdout);
fputs(serr, stdout);
fputs("\n", stdout);
}
strcpy(serr_save, serr);
} else if (*serr != '\0' && *serr_warn == '\0') {
if (strstr(serr, "'seorbel.txt' not found") == NULL)
strcpy(serr_warn, serr);
}
if (diff_mode) {
iflgret = swe_calc(te, ipldiff, iflag, x2, serr);
if (iflgret < 0) {
fputs("error: ", stdout);
fputs(serr, stdout);
fputs("\n", stdout);
}
if (diff_mode == DIFF_DIFF) {
for (i = 1; i < 6; i++)
x[i] -= x2[i];
if ((iflag & SEFLG_RADIANS) == 0)
x[0] = swe_difdeg2n(x[0], x2[0]);
else
x[0] = swe_difrad2n(x[0], x2[0]);
} else { /* DIFF_MIDP */
for (i = 1; i < 6; i++)
x[i] = (x[i] + x2[i]) / 2;
if ((iflag & SEFLG_RADIANS) == 0)
x[0] = swe_deg_midp(x[0], x2[0]);
else
x[0] = swe_rad_midp(x[0], x2[0]);
}
}
/* equator position */
if (strpbrk(fmt, "aADdQ") != NULL) {
iflag2 = iflag | SEFLG_EQUATORIAL;
if (ipl == SE_FIXSTAR)
iflgret = swe_fixstar(star, te, iflag2, xequ, serr);
else
iflgret = swe_calc(te, ipl, iflag2, xequ, serr);
if (diff_mode) {
iflgret = swe_calc(te, ipldiff, iflag2, x2, serr);
if (diff_mode == DIFF_DIFF) {
for (i = 1; i < 6; i++)
xequ[i] -= x2[i];
if ((iflag & SEFLG_RADIANS) == 0)
xequ[0] = swe_difdeg2n(xequ[0], x2[0]);
else
xequ[0] = swe_difrad2n(xequ[0], x2[0]);
} else { /* DIFF_MIDP */
for (i = 1; i < 6; i++)
xequ[i] = (xequ[i] + x2[i]) / 2;
if ((iflag & SEFLG_RADIANS) == 0)
xequ[0] = swe_deg_midp(xequ[0], x2[0]);
else
xequ[0] = swe_rad_midp(xequ[0], x2[0]);
}
}
}
/* azimuth and height */
if (strpbrk(fmt, "IiHhKk") != NULL) {
/* first, get topocentric equatorial positions */
iflgt = whicheph | SEFLG_EQUATORIAL | SEFLG_TOPOCTR;
if (ipl == SE_FIXSTAR)
iflgret = swe_fixstar(star, te, iflgt, xt, serr);
else
iflgret = swe_calc(te, ipl, iflgt, xt, serr);
/* to azimuth/height */
/* atmospheric pressure "0" has the effect that a value
* of 1013.25 mbar is assumed at 0 m above sea level.
* If the altitude of the observer is given (in geopos[2])
* pressure is estimated according to that */
swe_azalt(tut, SE_EQU2HOR, geopos, 0, 10, xt, xaz);
if (diff_mode) {
iflgret = swe_calc(te, ipldiff, iflgt, xt, serr);
swe_azalt(tut, SE_EQU2HOR, geopos, 0, 10, xt, x2);
if (diff_mode == DIFF_DIFF) {
for (i = 1; i < 3; i++)
xaz[i] -= x2[i];
if ((iflag & SEFLG_RADIANS) == 0)
xaz[0] = swe_difdeg2n(xaz[0], x2[0]);
else
xaz[0] = swe_difrad2n(xaz[0], x2[0]);
} else { /* DIFF_MIDP */
for (i = 1; i < 3; i++)
xaz[i] = (xaz[i] + x2[i]) / 2;
if ((iflag & SEFLG_RADIANS) == 0)
xaz[0] = swe_deg_midp(xaz[0], x2[0]);
else
xaz[0] = swe_rad_midp(xaz[0], x2[0]);
}
}
}
/* ecliptic cartesian position */
if (strpbrk(fmt, "XU") != NULL) {
iflag2 = iflag | SEFLG_XYZ;
if (ipl == SE_FIXSTAR)
iflgret = swe_fixstar(star, te, iflag2, xcart, serr);
else
iflgret = swe_calc(te, ipl, iflag2, xcart, serr);
if (diff_mode) {
iflgret = swe_calc(te, ipldiff, iflag2, x2, serr);
if (diff_mode == DIFF_DIFF) {
for (i = 0; i < 6; i++)
xcart[i] -= x2[i];
} else {
xcart[i] = (xcart[i] + x2[i]) / 2;
}
}
}
/* equator cartesian position */
if (strpbrk(fmt, "xu") != NULL) {
iflag2 = iflag | SEFLG_XYZ | SEFLG_EQUATORIAL;
if (ipl == SE_FIXSTAR)
iflgret = swe_fixstar(star, te, iflag2, xcartq, serr);
else
iflgret = swe_calc(te, ipl, iflag2, xcartq, serr);
if (diff_mode) {
iflgret = swe_calc(te, ipldiff, iflag2, x2, serr);
if (diff_mode == DIFF_DIFF) {
for (i = 0; i < 6; i++)
xcartq[i] -= x2[i];
} else {
xcartq[i] = (xcart[i] + x2[i]) / 2;
}
}
}
/* house position */
if (strpbrk(fmt, "gGj") != NULL) {
armc = swe_degnorm(swe_sidtime(tut) * 15 + geopos[0]);
for (i = 0; i < 6; i++)
xsv[i] = x[i];
if (hpos_meth == 1)
xsv[1] = 0;
if (ipl == SE_FIXSTAR)
strcpy(star2, star);
else
*star2 = '\0';
if (hpos_meth >= 2 && tolower(ihsy) == 'g') {
swe_gauquelin_sector(tut, ipl, star2, iflag, hpos_meth, geopos, 0, 0, &hposj, serr);
} else {
hposj = swe_house_pos(armc, geopos[1], xobl[0], ihsy, xsv, serr);
}
if (tolower(ihsy) == 'g')
hpos = (hposj - 1) * 10;
else
hpos = (hposj - 1) * 30;
if (diff_mode) {
for (i = 0; i < 6; i++)
xsv[i] = x2[i];
if (hpos_meth == 1)
xsv[1] = 0;
hpos2 = swe_house_pos(armc, geopos[1], xobl[0], ihsy, xsv, serr);
if (tolower(ihsy) == 'g')
hpos2 = (hpos2 - 1) * 10;
else
hpos2 = (hpos2 - 1) * 30;
if (diff_mode == DIFF_DIFF) {
if ((iflag & SEFLG_RADIANS) == 0)
hpos = swe_difdeg2n(hpos, hpos2);
else
hpos = swe_difrad2n(hpos, hpos2);
} else { /* DIFF_MIDP */
if ((iflag & SEFLG_RADIANS) == 0)
hpos = swe_deg_midp(hpos, hpos2);
else
hpos = swe_rad_midp(hpos, hpos2);
}
}
}
strcpy(spnam, se_pname);
print_line(0);
line_count++;
if (line_count >= line_limit) {
printf("****** line count %d was exceeded\n", line_limit);
break;
}
} /* for psp */
if (do_houses) {
double cusp[100];
int iofs;
if (tolower(ihsy) == 'g')
nhouses = 36;
iofs = nhouses + 1;
iflgret = swe_houses_ex(t,iflag, top_lat, top_long, ihsy, cusp, cusp+iofs);
if (iflgret < 0) {
if (strcmp(serr, serr_save) != 0 ) {
fputs("error: ", stdout);
fputs(serr, stdout);
fputs("\n", stdout);
}
strcpy(serr_save, serr);
} else {
for (ipl = 1; ipl < iofs+4; ipl++) {
x[0] = cusp[ipl];
x[1] = 0; /* latitude */
x[2] = 1.0; /* pseudo radius vector */
if (strpbrk(fmt, "aADdQ") != NULL) {
swe_cotrans(x, xequ, -xobl[0]);
}
print_line(MODE_HOUSE);
line_count++;
}
}
}
if (line_count >= line_limit)
break;
} /* for tjd */
if (*serr_warn != '\0') {
printf("\nwarning: ");
fputs(serr_warn,stdout);
printf("\n");
}
} /* while 1 */
/* close open files and free allocated space */
end_main:
swe_close();
return OK;
}
/*
* The string fmt contains a sequence of format specifiers;
* each character in fmt creates a column, the columns are
* sparated by the gap string.
*/
static int print_line(int mode)
{
char *sp, *sp2;
double t2, ju2 = 0;
double y_frac;
double ar, sinp;
double dret[20];
AS_BOOL is_house = ((mode & MODE_HOUSE) != 0);
AS_BOOL is_label = ((mode & MODE_LABEL) != 0);
int32 iflgret;
for (sp = fmt; *sp != '\0'; sp++) {
if (is_house && strchr("bBsSrRxXuUQnNfF+-*/=", *sp) != NULL) continue;
if (sp != fmt)
fputs(gap,stdout);
switch(*sp) {
case 'y':
if (is_label) { printf("year"); break; }
printf("%d", jyear);
break;
case 'Y':
if (is_label) { printf("year"); break; }
t2 = swe_julday(jyear,1,1,ju2,gregflag);
y_frac = (t - t2) / 365.0;
printf("%.2f", jyear + y_frac);
break;
case 'p':
if (is_label) { printf("obj.nr"); break; }
if (! is_house && diff_mode == DIFF_DIFF) {
printf("%d-%d", ipl, ipldiff);
} else if (! is_house && diff_mode == DIFF_MIDP) {
printf("%d/%d", ipl, ipldiff);
} else {
printf("%d", ipl);
}
break;
case 'P':
if (is_label) { printf("%-15s", "name"); break; }
if (is_house) {
if (ipl <= nhouses) {
printf("house %2d ", ipl);
} else {
printf("%-15s", hs_nam[ipl - nhouses]);
}
} else if (diff_mode == DIFF_DIFF) {
printf("%.3s-%.3s", spnam, spnam2);
} else if (diff_mode == DIFF_MIDP) {
printf("%.3s/%.3s", spnam, spnam2);
} else {
printf("%-15s", spnam);
}
break;
case 'J':
if (is_label) { printf("julday"); break; }
y_frac = (t - floor(t)) * 100;
if (floor(y_frac) != y_frac) {
printf("%.5f", t);
} else {
printf("%.2f", t);
}
break;
case 'T':
if (is_label) { printf("date"); break; }
printf("%02d.%02d.%d", jday, jmon, jyear);
if (jut != 0) {
int h, m, s;
s = (int) (jut * 3600 + 0.5);
h = (int) (s / 3600.0);
m = (int) ((s % 3600) / 60.0);
s %= 60;
printf(" %d:%02d:%02d", h, m, s);
if (universal_time)
printf(" UT");
else
printf(" ET");
}
break;
case 't':
if (is_label) { printf("date"); break; }
printf("%02d%02d%02d", jyear % 100, jmon, jday);
break;
case 'L':
if (is_label) { printf("long."); break; }
if (*psp == 'q' || *psp == 'y') /* delta t or time equation */
goto ldec;
fputs(dms(x[0], round_flag),stdout);
break;
case 'l':
if (is_label) { printf("long"); break; }
ldec:
printf("%# 11.7f", x[0]);
break;
case 'G':
if (is_label) { printf("housPos"); break; }
fputs(dms(hpos, round_flag),stdout);
break;
case 'g':
if (is_label) { printf("housPos"); break; }
printf("%# 11.7f", hpos);
break;
case 'j':
if (is_label) { printf("houseNr"); break; }
printf("%# 11.7f", hposj);
break;
case 'Z':
if (is_label) { printf("long"); break; }
fputs(dms(x[0], round_flag|BIT_ZODIAC),stdout);
break;
case 'S':
case 's':
if (*(sp+1) == 'S' || *(sp+1) == 's' || strpbrk(fmt, "XUxu") != NULL) {
for (sp2 = fmt; *sp2 != '\0'; sp2++) {
if (sp2 != fmt)
fputs(gap,stdout);
switch(*sp2) {
case 'L': /* speed! */
case 'Z': /* speed! */
if (is_label) { printf("lon/day"); break; }
fputs(dms(x[3], round_flag),stdout);
break;
case 'l': /* speed! */
if (is_label) { printf("lon/day"); break; }
printf("%11.7f", x[3]);
break;
case 'B': /* speed! */
if (is_label) { printf("lat/day"); break; }
fputs(dms(x[4], round_flag),stdout);
break;
case 'b': /* speed! */
if (is_label) { printf("lat/day"); break; }
printf("%11.7f", x[4]);
break;
case 'A': /* speed! */
if (is_label) { printf("RA/day"); break; }
fputs(dms(xequ[3]/15, round_flag|SEFLG_EQUATORIAL),stdout);
break;
case 'a': /* speed! */
if (is_label) { printf("RA/day"); break; }
printf("%11.7f", xequ[3]);
break;
case 'D': /* speed! */
if (is_label) { printf("dcl/day"); break; }
fputs(dms(xequ[4], round_flag),stdout);
break;
case 'd': /* speed! */
if (is_label) { printf("dcl/day"); break; }
printf("%11.7f", xequ[4]);
break;
case 'R': /* speed! */
case 'r': /* speed! */
if (is_label) { printf("AU/day"); break; }
printf("%# 14.9f", x[5]);
break;
case 'U': /* speed! */
case 'X': /* speed! */
if (is_label) {
fputs("speed_0", stdout);
fputs(gap, stdout);
fputs("speed_1", stdout);
fputs(gap, stdout);
fputs("speed_2", stdout);
break;
}
if (*sp =='U')
ar = sqrt(square_sum(xcart));
else
ar = 1;
printf("%# 14.9f", xcart[3]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcart[4]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcart[5]/ar);
break;
case 'u': /* speed! */
case 'x': /* speed! */
if (is_label) {
fputs("speed_0", stdout);
fputs(gap, stdout);
fputs("speed_1", stdout);
fputs(gap, stdout);
fputs("speed_2", stdout);
break;
}
if (*sp =='u')
ar = sqrt(square_sum(xcartq));
else
ar = 1;
printf("%# 14.9f", xcartq[3]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcartq[4]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcartq[5]/ar);
break;
default:
break;
}
}
if (*(sp+1) == 'S' || *(sp+1) == 's')
sp++;
} else if (*sp == 'S') {
if (is_label) { printf("deg/day"); break; }
fputs(dms(x[3], round_flag),stdout);
} else {
if (is_label) { printf("deg/day"); break; }
printf("%# 11.7f", x[3]);
}
break;
case 'B':
if (is_label) { printf("lat"); break; }
fputs(dms(x[1], round_flag),stdout);
break;
case 'b':
if (is_label) { printf("lat"); break; }
printf("%# 11.7f", x[1]);
break;
case 'A': /* right ascension */
if (is_label) { printf("RA"); break; }
fputs(dms(xequ[0]/15, round_flag|SEFLG_EQUATORIAL),stdout);
break;
case 'a': /* right ascension */
if (is_label) { printf("RA"); break; }
printf("%# 11.7f", xequ[0]);
break;
case 'D': /* declination */
if (is_label) { printf("decl"); break; }
fputs(dms(xequ[1], round_flag),stdout);
break;
case 'd': /* declination */
if (is_label) { printf("decl"); break; }
printf("%# 11.7f", xequ[1]);
break;
case 'I': /* azimuth */
if (is_label) { printf("azimuth"); break; }
fputs(dms(xaz[0], round_flag),stdout);
break;
case 'i': /* azimuth */
if (is_label) { printf("azimuth"); break; }
printf("%# 11.7f", xaz[0]);
break;
case 'H': /* height */
if (is_label) { printf("height"); break; }
fputs(dms(xaz[1], round_flag),stdout);
break;
case 'h': /* height */
if (is_label) { printf("height"); break; }
printf("%# 11.7f", xaz[1]);
break;
case 'K': /* height (apparent) */
if (is_label) { printf("hgtApp"); break; }
fputs(dms(xaz[2], round_flag),stdout);
break;
case 'k': /* height (apparent) */
if (is_label) { printf("hgtApp"); break; }
printf("%# 11.7f", xaz[2]);
break;
case 'R':
if (is_label) { printf("distAU"); break; }
printf("%# 14.9f", x[2]);
break;
case 'r':
if (is_label) { printf("dist"); break; }
if ( ipl == SE_MOON ) { /* for moon print parallax */
/* geocentric horizontal parallax: */
if (0) {
sinp = 8.794 / x[2]; /* in seconds of arc */
ar = sinp * (1 + sinp * sinp * 3.917402e-12);
/* the factor is 1 / (3600^2 * (180/pi)^2 * 6) */
printf("%# 13.5f\" %# 13.5f'", ar, ar/60.0);
}
swe_pheno(te, ipl, iflag, dret, serr);
printf("%# 13.5f\"", dret[5] * 3600);
} else {
printf("%# 14.9f", x[2]);
}
break;
case 'U':
case 'X':
if (*sp =='U')
ar = sqrt(square_sum(xcart));
else
ar = 1;
printf("%# 14.9f", xcart[0]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcart[1]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcart[2]/ar);
break;
case 'u':
case 'x':
if (is_label) {
fputs("x0", stdout);
fputs(gap, stdout);
fputs("x1", stdout);
fputs(gap, stdout);
fputs("x2", stdout);
break;
}
if (*sp =='u')
ar = sqrt(square_sum(xcartq));
else
ar = 1;
printf("%# 14.9f", xcartq[0]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcartq[1]/ar);
fputs(gap,stdout);
printf("%# 14.9f", xcartq[2]/ar);
break;
case 'Q':
if (is_label) { printf("Q"); break; }
printf("%-15s", spnam);
fputs(dms(x[0], round_flag),stdout);
fputs(dms(x[1], round_flag),stdout);
printf(" %# 14.9f", x[2]);
fputs(dms(x[3], round_flag),stdout);
fputs(dms(x[4], round_flag),stdout);
printf(" %# 14.9f\n", x[5]);
printf(" %s", dms(xequ[0], round_flag));
fputs(dms(xequ[1], round_flag),stdout);
printf(" %s", dms(xequ[3], round_flag));
fputs(dms(xequ[4], round_flag),stdout);
break;
case 'N':
case 'n': {
double xasc[6], xdsc[6];
int imeth = (*sp == tolower(*sp))?SE_NODBIT_MEAN:SE_NODBIT_OSCU;
iflgret = swe_nod_aps(te, ipl, iflag, imeth, xasc, xdsc, NULL, NULL, serr);
if (iflgret >= 0 && (ipl <= SE_NEPTUNE || *sp == 'N') ) {
if (is_label) {
fputs("nodAsc", stdout);
fputs(gap, stdout);
fputs("nodDesc", stdout);
break;
}
printf("%# 11.7f", xasc[0]);
fputs(gap,stdout);
printf("%# 11.7f", xdsc[0]);
}
};
break;
case 'F':
case 'f':
if (! is_house) {
double xfoc[6], xaph[6], xper[6];
int imeth = (*sp == tolower(*sp))?SE_NODBIT_MEAN:SE_NODBIT_OSCU;
iflgret = swe_nod_aps(te, ipl, iflag, imeth, NULL, NULL, xper, xaph, serr);
if (iflgret >= 0 && (ipl <= SE_NEPTUNE || *sp == 'F') ) {
if (is_label) {
fputs("peri", stdout);
fputs(gap, stdout);
fputs("apo", stdout);
break;
}
printf("%# 11.7f", xper[0]);
fputs(gap,stdout);
printf("%# 11.7f", xaph[0]);
}
imeth |= SE_NODBIT_FOPOINT;
iflgret = swe_nod_aps(te, ipl, iflag, imeth, NULL, NULL, xper, xfoc, serr);
if (iflgret >= 0 && (ipl <= SE_NEPTUNE || *sp == 'F') ) {
if (is_label) {
fputs(gap, stdout);
fputs("focus", stdout);
break;
}
fputs(gap,stdout);
printf("%# 11.7f", xfoc[0]);
}
};
break;
case '+':
if (is_house) break;
if (is_label) { printf("phase"); break; }
fputs(dms(attr[0], round_flag),stdout);
break;
case '-':
if (is_label) { printf("phase"); break; }
if (is_house) break;
printf(" %# 14.9f", attr[1]);
break;
case '*':
if (is_label) { printf("elong"); break; }
if (is_house) break;
fputs(dms(attr[2], round_flag),stdout);
break;
case '/':
if (is_label) { printf("diamet"); break; }
if (is_house) break;
fputs(dms(attr[3], round_flag),stdout);
break;
case '=':
if (is_label) { printf("magn"); break; }
if (is_house) break;
printf(" %# 6.1fm", attr[4]);
break;
case 'V': /* human design gates */
case 'v': {
double xhds;
int igate, iline, ihex;
static int hexa[64] = {1, 43, 14, 34, 9, 5, 26, 11, 10, 58, 38, 54, 61, 60, 41, 19, 13, 49, 30, 55, 37, 63, 22, 36, 25, 17, 21, 51, 42, 3, 27, 24, 2, 23, 8, 20, 16, 35, 45, 12, 15, 52, 39, 53, 62, 56, 31, 33, 7, 4, 29, 59, 40, 64, 47, 6, 46, 18, 48, 57, 32, 50, 28, 44};
if (is_label) { printf("hds"); break; }
if (is_house) break;
xhds = swe_degnorm(x[0] - 223.25);
ihex = (int) floor(xhds / 5.625);
iline = ((int) (floor(xhds / 0.9375))) % 6 + 1 ;
igate = hexa[ihex];
printf("%2d.%d", igate, iline);
if (*sp == 'V')
printf(" %2d%%", swe_d2l(100 * fmod(xhds / 0.9375, 1)));
break;
}
} /* switch */
} /* for sp */
printf("\n");
return OK;
}
#define OUTPUT_EXTRA_PRECISION 0
static char *dms(double xv, int32 iflg)
{
int izod;
int32 k, kdeg, kmin, ksec;
char *c = ODEGREE_STRING;
char *sp, s1[50];
static char s[50];
int sgn;
#if MSDOS
if (_isnan(xv))
return "nan";
#else
if (isnan(xv))
return "nan";
#endif
*s = '\0';
if (iflg & SEFLG_EQUATORIAL)
c = "h";
if (xv < 0) {
xv = -xv;
sgn = -1;
} else
sgn = 1;
if (iflg & BIT_ROUND_MIN)
xv = swe_degnorm(xv + 0.5/60);
if (iflg & BIT_ROUND_SEC)
xv = swe_degnorm(xv + 0.5/3600);
if (iflg & BIT_ZODIAC) {
izod = (int) (xv / 30);
xv = fmod(xv, 30);
kdeg = (int32) xv;
sprintf(s, "%2d %s ", kdeg, zod_nam[izod]);
} else {
kdeg = (int32) xv;
sprintf(s, " %3d%s", kdeg, c);
}
xv -= kdeg;
xv *= 60;
kmin = (int32) xv;
if ((iflg & BIT_ZODIAC) && (iflg & BIT_ROUND_MIN)) {
sprintf(s1, "%2d", kmin);
} else {
sprintf(s1, "%2d'", kmin);
}
strcat(s, s1);
if (iflg & BIT_ROUND_MIN)
goto return_dms;
xv -= kmin;
xv *= 60;
ksec = (int32) xv;
if (iflg & BIT_ROUND_SEC) {
sprintf(s1, "%2d\"", ksec);
} else {
sprintf(s1, "%2d", ksec);
}
strcat(s, s1);
if (iflg & BIT_ROUND_SEC)
goto return_dms;
xv -= ksec;
#if OUTPUT_EXTRA_PRECISION
k = (int32) (xv * 100000 + 0.5);
sprintf(s1, ".%05d", k);
#else
k = (int32) (xv * 10000 + 0.5);
sprintf(s1, ".%04d", k);
#endif
strcat(s, s1);
return_dms:;
if (sgn < 0) {
sp = strpbrk(s, "0123456789");
*(sp-1) = '-';
}
if (iflg & BIT_LZEROES) {
while ((sp = strchr(s+2, ' ')) != NULL) *sp = '0';
}
return(s);
}
static int letter_to_ipl(int letter)
{
if (letter >= '0' && letter <= '9')
return letter - '0' + SE_SUN;
if (letter >= 'A' && letter <= 'I')
return letter - 'A' + SE_MEAN_APOG;
if (letter >= 'J' && letter <= 'Z')
return letter - 'J' + SE_CUPIDO;
switch (letter) {
case 'm': return SE_MEAN_NODE;
case 'c': return SE_INTP_APOG;
case 'g': return SE_INTP_PERG;
case 'n':
case 'o': return SE_ECL_NUT;
case 't': return SE_TRUE_NODE;
case 'f': return SE_FIXSTAR;
case 'w': return SE_WALDEMATH;
case 'e': /* swetest: a line of labels */
case 'q': /* swetest: delta t */
case 'y': /* swetest: time equation */
case 's': /* swetest: an asteroid, with number given in -xs[number] */
case 'z': /* swetest: a fictitious body, number given in -xz[number] */
case 'd': /* swetest: default (main) factors 0123456789mtABC */
case 'p': /* swetest: main factors ('d') plus main asteroids DEFGHI */
case 'h': /* swetest: fictitious factors JKLMNOPQRSTUVWXYZw */
case 'a': /* swetest: all factors, like 'p'+'h' */
return -1;
}
return -2;
}
static int32 ut_to_lmt_lat(double t_ut, double *geopos, double *t_ret, char *serr)
{
int32 iflgret = OK;
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
t_ut += geopos[0] / 360.0;
if (time_flag & BIT_TIME_LAT) {
iflgret = swe_lmt_to_lat(t_ut, geopos[0], &t_ut, NULL);
}
}
*t_ret = t_ut;
return iflgret;
}
static int32 call_rise_set(double t_ut, int32 ipl, char *star, int32 whicheph, int32 special_mode, double *geopos, char *serr)
{
int ii;
int32 rsmi = 0;
double tret[10];
double t0, t1;
int32 retc = OK;
swe_set_topo(geopos[0], geopos[1], geopos[2]);
do_printf("\n");
/* loop over days */
for (ii = 0; ii < nstep; ii++, t_ut = tret[1] + 0.1) {
*sout = '\0';
/* swetest -rise
* calculate and print rising and setting */
if (special_event == SP_RISE_SET) {
/* rising */
rsmi = SE_CALC_RISE;
if (norefrac) rsmi |= SE_BIT_NO_REFRACTION;
if (disccenter) rsmi |= SE_BIT_DISC_CENTER;
if (swe_rise_trans(t_ut, ipl, star, whicheph, rsmi, geopos, 1013.25, 10, &(tret[0]), serr) != OK) {
do_printf(serr);
exit(0);
}
/* setting */
rsmi = SE_CALC_SET;
if (norefrac) rsmi |= SE_BIT_NO_REFRACTION;
if (disccenter) rsmi |= SE_BIT_DISC_CENTER;
if (swe_rise_trans(t_ut, ipl, star, whicheph, rsmi, geopos, 1013.25, 10, &(tret[1]), serr) != OK) {
do_printf(serr);
exit(0);
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(tret[0], geopos, &(tret[0]), serr);
retc = ut_to_lmt_lat(tret[1], geopos, &(tret[1]), serr);
}
t0 = 0; t1 = 0;
strcpy(sout, "rise ");
if (tret[0] == 0 || tret[0] > tret[1]) {
strcat(sout, " - ");
} else {
t0 = tret[0];
swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s ", jday, jmon, jyear, hms(jut,BIT_LZEROES));
}
strcat(sout, "set ");
if (tret[1] == 0) {
strcat(sout, " - \n");
} else {
t1 = tret[1];
swe_revjul(tret[1], gregflag, &jyear, &jmon, &jday, &jut);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s ", jday, jmon, jyear, hms(jut,BIT_LZEROES));
}
if (t0 != 0 && t1 != 0) {
t0 = (t1 - t0) * 24;
sprintf(sout + strlen(sout), "dt = %s", hms(t0, BIT_LZEROES));
}
strcat(sout, "\n");
do_printf(sout);
}
/* swetest -metr
* calculate and print transits over meridian (midheaven and lower
* midheaven */
if (special_event == SP_MERIDIAN_TRANSIT) {
/* transit over midheaven */
if (swe_rise_trans(t_ut, ipl, star, whicheph, SE_CALC_MTRANSIT, geopos, 1013.25, 10, &(tret[0]), serr) != OK) {
do_printf(serr);
return ERR;
}
/* transit over lower midheaven */
if (swe_rise_trans(t_ut, ipl, star, whicheph, SE_CALC_ITRANSIT, geopos, 1013.25, 10, &(tret[1]), serr) != OK) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(tret[0], geopos, &(tret[0]), serr);
retc = ut_to_lmt_lat(tret[1], geopos, &(tret[1]), serr);
}
strcpy(sout, "mtransit ");
if (tret[0] == 0) strcat(sout, " - ");
else {
swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s ", jday, jmon, jyear, hms(jut,BIT_LZEROES));
}
strcat(sout, "itransit ");
if (tret[1] == 0) strcat(sout, " - \n");
else {
swe_revjul(tret[1], gregflag, &jyear, &jmon, &jday, &jut);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\n", jday, jmon, jyear, hms(jut,BIT_LZEROES));
}
do_printf(sout);
}
}
return retc;
}
static int32 call_lunar_eclipse(double t_ut, int32 whicheph, int32 special_mode, double *geopos, char *serr)
{
int i, ii, retc, eclflag, ecl_type = 0;
int ihou, imin, isec, isgn;
double dfrc, attr[30], dt;
char s1[AS_MAXCH], sout_short[AS_MAXCH], sfmt[AS_MAXCH];
/* no selective eclipse type set, set all */
if ((search_flag & SE_ECL_ALLTYPES_LUNAR) == 0)
search_flag |= SE_ECL_ALLTYPES_LUNAR;
do_printf("\n");
for (ii = 0; ii < nstep; ii++, t_ut += direction) {
*sout = '\0';
/* swetest -lunecl -how
* type of lunar eclipse and percentage for a given time: */
if (special_mode & SP_MODE_HOW) {
if ((eclflag = swe_lun_eclipse_how(t_ut, whicheph, geopos, attr, serr)) ==
ERR) {
do_printf(serr);
return ERR;
} else {
if (eclflag & SE_ECL_TOTAL) {
ecl_type = ECL_LUN_TOTAL;
strcpy(sfmt, "total lunar eclipse: %f o/o \n");
} else if (eclflag & SE_ECL_PARTIAL) {
ecl_type = ECL_LUN_PARTIAL;
strcpy(sfmt, "partial lunar eclipse: %f o/o \n");
} else if (eclflag & SE_ECL_PENUMBRAL) {
ecl_type = ECL_LUN_PENUMBRAL;
strcpy(sfmt, "penumbral lunar eclipse: %f o/o \n");
} else {
strcpy(sfmt, "no lunar eclipse \n");
}
strcpy(sout, sfmt);
if (strchr(sfmt, '%') != NULL) {
sprintf(sout, sfmt, attr[0]);
}
do_printf(sout);
}
continue;
}
/* swetest -lunecl
* find next lunar eclipse: */
/* locally visible lunar eclipse */
if (special_mode & SP_MODE_LOCAL) {
if ((eclflag = swe_lun_eclipse_when_loc(t_ut, whicheph, geopos,
tret, attr, direction_flag, serr)) == ERR) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
for (i = 0; i < 10; i++) {
if (tret[i] != 0)
retc = ut_to_lmt_lat(tret[i], geopos, &(tret[i]), serr);
}
}
t_ut = tret[0];
if ((eclflag & SE_ECL_TOTAL)) {
strcpy(sout, "total ");
ecl_type = ECL_LUN_TOTAL;
}
if ((eclflag & SE_ECL_PENUMBRAL)) {
strcpy(sout, "penumb. ");
ecl_type = ECL_LUN_PENUMBRAL;
}
if ((eclflag & SE_ECL_PARTIAL)) {
strcpy(sout, "partial ");
ecl_type = ECL_LUN_PARTIAL;
}
strcat(sout, "lunar eclipse\t");
swe_revjul(t_ut, gregflag, &jyear, &jmon, &jday, &jut);
/*if ((eclflag = swe_lun_eclipse_how(t_ut, whicheph, geopos, attr, serr)) == ERR) {
do_printf(serr);
return ERR;
}*/
dt = (tret[3] - tret[2]) * 24 * 60;
sprintf(s1, "%d min %4.2f sec", (int) dt, fmod(dt, 1) * 60);
/* short output:
* date, time of day, umbral magnitude, umbral duration, saros series, member number */
sprintf(sout_short, "%s\t%2d.%2d.%4d\t%s\t%.3f\t%s\t%d\t%d\n", sout, jday, jmon, jyear, hms(jut,0), attr[8],s1, (int) attr[9], (int) attr[10]);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%.4f/%.4f\tsaros %d/%d\t%.6f\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[0],attr[1], (int) attr[9], (int) attr[10], t_ut);
/* second line:
* eclipse times, penumbral, partial, total begin and end */
if (eclflag & SE_ECL_PENUMBBEG_VISIBLE)
sprintf(sout + strlen(sout), " %s ", hms_from_tjd(tret[6]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_PARTBEG_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_TOTBEG_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_TOTEND_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_PARTEND_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_PENUMBEND_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[7]));
else
strcat(sout, " - ");
strcat(sout, "\n");
/* global lunar eclipse */
} else {
if ((eclflag = swe_lun_eclipse_when(t_ut, whicheph, search_flag,
tret, direction_flag, serr)) == ERR) {
do_printf(serr);
return ERR;
}
t_ut = tret[0];
if ((eclflag & SE_ECL_TOTAL)) {
strcpy(sout, "total ");
ecl_type = ECL_LUN_TOTAL;
}
if ((eclflag & SE_ECL_PENUMBRAL)) {
strcpy(sout, "penumb. ");
ecl_type = ECL_LUN_PENUMBRAL;
}
if ((eclflag & SE_ECL_PARTIAL)) {
strcpy(sout, "partial ");
ecl_type = ECL_LUN_PARTIAL;
}
strcat(sout, "lunar eclipse\t");
if ((eclflag = swe_lun_eclipse_how(t_ut, whicheph, geopos, attr, serr)) ==
ERR) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
for (i = 0; i < 10; i++) {
if (tret[i] != 0)
retc = ut_to_lmt_lat(tret[i], geopos, &(tret[i]), serr);
}
}
t_ut = tret[0];
swe_revjul(t_ut, gregflag, &jyear, &jmon, &jday, &jut);
dt = (tret[3] - tret[2]) * 24 * 60;
sprintf(s1, "%d min %4.2f sec", (int) dt, fmod(dt, 1) * 60);
/* short output:
* date, time of day, umbral magnitude, umbral duration, saros series, member number */
sprintf(sout_short, "%s\t%2d.%2d.%4d\t%s\t%.3f\t%s\t%d\t%d\n", sout, jday, jmon, jyear, hms(jut,0), attr[8],s1, (int) attr[9], (int) attr[10]);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%.4f/%.4f\tsaros %d/%d\t%.6f\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[0],attr[1], (int) attr[9], (int) attr[10], t_ut);
/* second line:
* eclipse times, penumbral, partial, total begin and end */
sprintf(sout + strlen(sout), " %s ", hms_from_tjd(tret[6]));
if (tret[2] != 0)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2]));
else
strcat(sout, " - ");
if (tret[4] != 0)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4]));
else
strcat(sout, " - ");
if (tret[5] != 0)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5]));
else
strcat(sout, " - ");
if (tret[3] != 0)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3]));
else
strcat(sout, " - ");
sprintf(sout + strlen(sout), "%s\n", hms_from_tjd(tret[7]));
if (special_mode & SP_MODE_HOCAL) {
swe_split_deg(jut, SE_SPLIT_DEG_ROUND_MIN, &ihou, &imin, &isec, &dfrc, &isgn);
sprintf(sout, "\"%04d %02d %02d %02d.%02d %d\",\n", jyear, jmon, jday, ihou, imin, ecl_type);
}
}
if (short_output)
do_printf(sout_short);
else
do_printf(sout);
}
return OK;
}
static int32 call_solar_eclipse(double t_ut, int32 whicheph, int32 special_mode, double *geopos, char *serr)
{
int i, ii, retc, eclflag, ecl_type = 0;
double dt, tret[30], attr[30], geopos_max[3];
char s1[AS_MAXCH], s2[AS_MAXCH], sout_short[AS_MAXCH];
AS_BOOL has_found = FALSE;
/* no selective eclipse type set, set all */
if ((search_flag & SE_ECL_ALLTYPES_SOLAR) == 0)
search_flag |= SE_ECL_ALLTYPES_SOLAR;
/* for local eclipses: set geographic position of observer */
if (special_mode & SP_MODE_LOCAL)
swe_set_topo(geopos[0], geopos[1], geopos[2]);
do_printf("\n");
for (ii = 0; ii < nstep; ii++, t_ut += direction) {
*sout = '\0';
/* swetest -solecl -local -geopos...
* find next solar eclipse observable from a given geographic position */
if (special_mode & SP_MODE_LOCAL) {
if ((eclflag = swe_sol_eclipse_when_loc(t_ut, whicheph, geopos, tret,
attr, direction_flag, serr)) == ERR) {
do_printf(serr);
return ERR;
} else {
has_found = FALSE;
t_ut = tret[0];
if ((search_flag & SE_ECL_TOTAL) && (eclflag & SE_ECL_TOTAL)) {
strcpy(sout, "total ");
has_found = TRUE;
ecl_type = ECL_SOL_TOTAL;
}
if ((search_flag & SE_ECL_ANNULAR) && (eclflag & SE_ECL_ANNULAR)) {
strcpy(sout, "annular ");
has_found = TRUE;
ecl_type = ECL_SOL_ANNULAR;
}
if ((search_flag & SE_ECL_PARTIAL) && (eclflag & SE_ECL_PARTIAL)) {
strcpy(sout, "partial ");
has_found = TRUE;
ecl_type = ECL_SOL_PARTIAL;
}
if (!has_found) {
ii--;
} else {
swe_calc(t_ut + swe_deltat(t_ut), SE_ECL_NUT, 0, x, serr);
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
for (i = 0; i < 10; i++) {
if (tret[i] != 0)
retc = ut_to_lmt_lat(tret[i], geopos, &(tret[i]), serr);
}
}
t_ut = tret[0];
swe_revjul(t_ut, gregflag, &jyear, &jmon, &jday, &jut);
dt = (tret[3] - tret[2]) * 24 * 60;
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%.4f/%.4f/%.4f\tsaros %d/%d\t%.6f\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[8], attr[0], attr[2], (int) attr[9], (int) attr[10], t_ut);
sprintf(sout + strlen(sout), "\t%d min %4.2f sec\t", (int) dt, fmod(dt, 1) * 60);
if (eclflag & SE_ECL_1ST_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[1]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_2ND_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_3RD_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_4TH_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4]));
else
strcat(sout, " - ");
#if 0
sprintf(sout + strlen(sout), "\t%d min %4.2f sec %s %s %s %s",
(int) dt, fmod(dt, 1) * 60,
strcpy(s1, hms(fmod(tret[1] + 0.5, 1) * 24, BIT_LZEROES)),
strcpy(s3, hms(fmod(tret[2] + 0.5, 1) * 24, BIT_LZEROES)),
strcpy(s4, hms(fmod(tret[3] + 0.5, 1) * 24, BIT_LZEROES)),
strcpy(s2, hms(fmod(tret[4] + 0.5, 1) * 24, BIT_LZEROES)));
#endif
strcat(sout, "\n");
do_printf(sout);
}
}
} /* endif search_local */
/* swetest -solecl
* find next solar eclipse observable from anywhere on earth */
if (!(special_mode & SP_MODE_LOCAL)) {
if ((eclflag = swe_sol_eclipse_when_glob(t_ut, whicheph, search_flag,
tret, direction_flag, serr)) == ERR) {
do_printf(serr);
return ERR;
}
t_ut = tret[0];
if ((eclflag & SE_ECL_TOTAL)) {
strcpy(sout, "total ");
ecl_type = ECL_SOL_TOTAL;
}
if ((eclflag & SE_ECL_ANNULAR)) {
strcpy(sout, "annular ");
ecl_type = ECL_SOL_ANNULAR;
}
if ((eclflag & SE_ECL_ANNULAR_TOTAL)) {
strcpy(sout, "ann-tot ");
ecl_type = ECL_SOL_ANNULAR; /* by Alois: what is this ? */
}
if ((eclflag & SE_ECL_PARTIAL)) {
strcpy(sout, "partial ");
ecl_type = ECL_SOL_PARTIAL;
}
if ((eclflag & SE_ECL_NONCENTRAL) && !(eclflag & SE_ECL_PARTIAL))
strcat(sout, "non-central ");
swe_sol_eclipse_where(t_ut, whicheph, geopos_max, attr, serr);
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
for (i = 0; i < 10; i++) {
if (tret[i] != 0)
retc = ut_to_lmt_lat(tret[i], geopos, &(tret[i]), serr);
}
}
swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
sprintf(sout_short, "%s\t%2d.%2d.%4d\t%s\t%.3f", sout, jday, jmon, jyear, hms(jut,0), attr[8]);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%f km\t%.4f/%.4f/%.4f\tsaros %d/%d\t%.6f\n", jday, jmon, jyear, hms(jut,0), attr[3], attr[8], attr[0], attr[2], (int) attr[9], (int) attr[10], tret[0]);
sprintf(sout + strlen(sout), "\t%s ", hms_from_tjd(tret[2]));
if (tret[4] != 0) {
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4]));
} else {
strcat(sout, " - ");
}
if (tret[5] != 0) {
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5]));
} else {
strcat(sout, " - ");
}
sprintf(sout + strlen(sout), "%s\n", hms_from_tjd(tret[3]));
sprintf(sout + strlen(sout), "\t%s\t%s", strcpy(s1, dms(geopos_max[0], BIT_ROUND_MIN)), strcpy(s2, dms(geopos_max[1], BIT_ROUND_MIN)));
strcat(sout, "\t");
strcat(sout_short, "\t");
if (!(eclflag & SE_ECL_PARTIAL) && !(eclflag & SE_ECL_NONCENTRAL)) {
if ((eclflag = swe_sol_eclipse_when_loc(t_ut - 10, whicheph, geopos_max, tret, attr, 0, serr)) == ERR) {
do_printf(serr);
return ERR;
}
if (fabs(tret[0] - t_ut) > 2)
do_printf("when_loc returns wrong date\n");
dt = (tret[3] - tret[2]) * 24 * 60;
sprintf(s1, "%d min %4.2f sec", (int) dt, fmod(dt, 1) * 60);
strcat(sout, s1);
strcat(sout_short, s1);
}
sprintf(sout_short + strlen(sout_short), "\t%d\t%d", (int) attr[9], (int) attr[10]);
strcat(sout, "\n");
strcat(sout_short, "\n");
if (special_mode & SP_MODE_HOCAL) {
int ihou, imin, isec, isgn;
double dfrc;
swe_split_deg(jut, SE_SPLIT_DEG_ROUND_MIN, &ihou, &imin, &isec, &dfrc, &isgn);
sprintf(sout, "\"%04d %02d %02d %02d.%02d %d\",\n", jyear, jmon, jday, ihou, imin, ecl_type);
}
/*printf("len=%ld\n", strlen(sout));*/
if (short_output)
do_printf(sout_short);
else
do_printf(sout);
}
}
return OK;
}
static int32 call_lunar_occultation(double t_ut, int32 ipl, char *star, int32 whicheph, int32 special_mode, double *geopos, char *serr)
{
int i, ii, ecl_type = 0, eclflag, retc;
double dt, tret[30], attr[30], geopos_max[3];
char s1[AS_MAXCH], s2[AS_MAXCH];
AS_BOOL has_found = FALSE;
/* no selective eclipse type set, set all */
if ((search_flag & SE_ECL_ALLTYPES_SOLAR) == 0)
search_flag |= SE_ECL_ALLTYPES_SOLAR;
/* for local occultations: set geographic position of observer */
if (special_mode & SP_MODE_LOCAL)
swe_set_topo(geopos[0], geopos[1], geopos[2]);
do_printf("\n");
for (ii = 0; ii < nstep; ii++) {
*sout = '\0';
if (special_mode & SP_MODE_LOCAL) {
if ((eclflag = swe_lun_occult_when_loc(t_ut, ipl, star, whicheph, geopos, tret, attr, direction_flag, serr)) == ERR) {
do_printf(serr);
return ERR;
} else {
t_ut = tret[0];
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
for (i = 0; i < 10; i++) {
if (tret[i] != 0)
retc = ut_to_lmt_lat(tret[i], geopos, &(tret[i]), serr);
}
}
has_found = FALSE;
*sout = '\0';
if ((search_flag & SE_ECL_TOTAL) && (eclflag & SE_ECL_TOTAL)) {
strcat(sout, "total");
has_found = TRUE;
ecl_type = ECL_SOL_TOTAL;
}
if ((search_flag & SE_ECL_ANNULAR) && (eclflag & SE_ECL_ANNULAR)) {
strcat(sout, "annular");
has_found = TRUE;
ecl_type = ECL_SOL_ANNULAR;
}
if ((search_flag & SE_ECL_PARTIAL) && (eclflag & SE_ECL_PARTIAL)) {
strcat(sout, "partial");
has_found = TRUE;
ecl_type = ECL_SOL_PARTIAL;
}
if (ipl != SE_SUN) {
if ((eclflag & SE_ECL_OCC_BEG_DAYLIGHT) && (eclflag & SE_ECL_OCC_END_DAYLIGHT))
strcat(sout, "(daytime)"); /* occultation occurs during the day */
else if (eclflag & SE_ECL_OCC_BEG_DAYLIGHT)
strcat(sout, "(sunset) "); /* occultation occurs during the day */
else if (eclflag & SE_ECL_OCC_END_DAYLIGHT)
strcat(sout, "(sunrise)"); /* occultation occurs during the day */
}
while (strlen(sout) < 17)
strcat(sout, " ");
if (!has_found) {
ii--;
} else {
swe_calc_ut(t_ut, SE_ECL_NUT, 0, x, serr);
swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
dt = (tret[3] - tret[2]) * 24 * 60;
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%fo/o\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[0]);
sprintf(sout + strlen(sout), "\t%d min %4.2f sec\t", (int) dt, fmod(dt, 1) * 60);
if (eclflag & SE_ECL_1ST_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[1]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_2ND_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_3RD_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3]));
else
strcat(sout, " - ");
if (eclflag & SE_ECL_4TH_VISIBLE)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4]));
else
strcat(sout, " - ");
#if 0
sprintf(sout + strlen(sout), "\t%d min %4.2f sec %s %s %s %s",
(int) dt, fmod(dt, 1) * 60,
strcpy(s1, hms(fmod(tret[1] + 0.5, 1) * 24, BIT_LZEROES)),
strcpy(s3, hms(fmod(tret[2] + 0.5, 1) * 24, BIT_LZEROES)),
strcpy(s4, hms(fmod(tret[3] + 0.5, 1) * 24, BIT_LZEROES)),
strcpy(s2, hms(fmod(tret[4] + 0.5, 1) * 24, BIT_LZEROES)));
#endif
strcat(sout, "\n");
do_printf(sout);
}
}
} /* endif search_local */
if (!(special_mode & SP_MODE_LOCAL)) {
/* * global search for occultations */
if ((eclflag = swe_lun_occult_when_glob(t_ut, ipl, star, whicheph, search_flag, tret, direction_flag, serr)) == ERR) {
do_printf(serr);
return ERR;
}
if ((eclflag & SE_ECL_TOTAL)) {
strcpy(sout, "total ");
ecl_type = ECL_SOL_TOTAL;
}
if ((eclflag & SE_ECL_ANNULAR)) {
strcpy(sout, "annular ");
ecl_type = ECL_SOL_ANNULAR;
}
if ((eclflag & SE_ECL_ANNULAR_TOTAL)) {
strcpy(sout, "ann-tot ");
ecl_type = ECL_SOL_ANNULAR; /* by Alois: what is this ? */
}
if ((eclflag & SE_ECL_PARTIAL)) {
strcpy(sout, "partial ");
ecl_type = ECL_SOL_PARTIAL;
}
if ((eclflag & SE_ECL_NONCENTRAL) && !(eclflag & SE_ECL_PARTIAL))
strcat(sout, "non-central ");
t_ut = tret[0];
swe_lun_occult_where(t_ut, ipl, star, whicheph, geopos_max, attr, serr);
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
for (i = 0; i < 10; i++) {
if (tret[i] != 0)
retc = ut_to_lmt_lat(tret[i], geopos, &(tret[i]), serr);
}
}
swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%f km\t%f o/o\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[3], attr[0]);
sprintf(sout + strlen(sout), "\t%s ", hms_from_tjd(tret[2]));
if (tret[4] != 0)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4]));
else
strcat(sout, " - ");
if (tret[5] != 0)
sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5]));
else
strcat(sout, " - ");
sprintf(sout + strlen(sout), "%s\n", hms_from_tjd(tret[3]));
sprintf(sout + strlen(sout), "\t%s\t%s", strcpy(s1, dms(geopos_max[0], BIT_ROUND_MIN)), strcpy(s2, dms(geopos_max[1], BIT_ROUND_MIN)));
if (!(eclflag & SE_ECL_PARTIAL) && !(eclflag & SE_ECL_NONCENTRAL)) {
if ((eclflag = swe_lun_occult_when_loc(t_ut - 10, ipl, star, whicheph, geopos_max, tret, attr, 0, serr)) == ERR) {
do_printf(serr);
return ERR;
}
if (fabs(tret[0] - t_ut) > 2)
do_printf("when_loc returns wrong date\n");
dt = (tret[3] - tret[2]) * 24 * 60;
sprintf(sout + strlen(sout), "\t%d min %4.2f sec\t", (int) dt, fmod(dt, 1) * 60);
}
strcat(sout, "\n");
if (special_mode & SP_MODE_HOCAL) {
int ihou, imin, isec, isgn;
double dfrc;
swe_split_deg(jut, SE_SPLIT_DEG_ROUND_MIN, &ihou, &imin, &isec, &dfrc, &isgn);
sprintf(sout, "\"%04d %02d %02d %02d.%02d %d\",\n", jyear, jmon, jday, ihou, imin, ecl_type);
}
do_printf(sout);
}
t_ut += direction;
}
return OK;
}
static void do_print_heliacal(double *dret, int32 event_type, char *obj_name)
{
char *sevtname[] = {"",
"heliacal rising ",
"heliacal setting",
"evening first ",
"morning last ",
"evening rising ",
"morning setting ",};
char *stz = "UT";
char stim0[40], stim1[40], stim2[40];
if (time_flag & BIT_TIME_LMT)
stz = "LMT";
if (time_flag & BIT_TIME_LAT)
stz = "LAT";
*sout = '\0';
swe_revjul(dret[0], gregflag, &jyear, &jmon, &jday, &jut);
if (event_type <= 4) {
if (hel_using_AV) {
strcpy(stim0, hms_from_tjd(dret[0]));
remove_whitespace(stim0);
/* The following line displays only the beginning of visibility. */
sprintf(sout + strlen(sout), "%s %s: %d/%02d/%02d %s %s (%.5f)\n", obj_name, sevtname[event_type], jyear, jmon, jday, stim0, stz, dret[0]);
} else {
/* display the moment of beginning and optimum visibility */
strcpy(stim0, hms_from_tjd(dret[0]));
strcpy(stim1, hms_from_tjd(dret[1]));
strcpy(stim2, hms_from_tjd(dret[2]));
remove_whitespace(stim0);
remove_whitespace(stim1);
remove_whitespace(stim2);
sprintf(sout + strlen(sout), "%s %s: %d/%02d/%02d %s %s (%.5f), opt %s, end %s, dur %.1f min\n", obj_name, sevtname[event_type], jyear, jmon, jday, stim0, stz, dret[0], stim1, stim2, (dret[2] - dret[0]) * 1440);
}
} else {
strcpy(stim0, hms_from_tjd(dret[0]));
remove_whitespace(stim0);
sprintf(sout + strlen(sout), "%s %s: %d/%02d/%02d %s %s (%f)\n", obj_name, sevtname[event_type], jyear, jmon, jday, stim0, stz, dret[0]);
}
do_printf(sout);
}
static int32 call_heliacal_event(double t_ut, int32 ipl, char *star, int32 whicheph, int32 special_mode, double *geopos, double *datm, double *dobs, char *serr)
{
int ii, retc, event_type = 0, retflag;
double dret[40], tsave1, tsave2 = 0;
char obj_name[AS_MAXCH];
helflag |= whicheph;
/* if invalid heliacal event type was required, set 0 for any type */
if (search_flag < 0 || search_flag > 6)
search_flag = 0;
/* optical instruments used: */
if (dobs[3] > 0)
helflag |= SE_HELFLAG_OPTICAL_PARAMS;
if (hel_using_AV)
helflag |= SE_HELFLAG_AV;
if (ipl == SE_FIXSTAR)
strcpy(obj_name, star);
else
swe_get_planet_name(ipl, obj_name);
do_printf("\n");
for (ii = 0; ii < nstep; ii++, t_ut = dret[0] + 1) {
*sout = '\0';
if (search_flag > 0)
event_type = search_flag;
else if (ipl == SE_MOON)
event_type = SE_EVENING_FIRST;
else
event_type = SE_HELIACAL_RISING;
retflag = swe_heliacal_ut(t_ut, geopos, datm, dobs, obj_name, event_type, helflag, dret, serr);
if (retflag == ERR) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(dret[0], geopos, &(dret[0]), serr);
retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
}
do_print_heliacal(dret, event_type, obj_name);
/* list all events within synodic cycle */
if (search_flag == 0) {
if (ipl == SE_VENUS || ipl == SE_MERCURY) {
/* we have heliacal rising (morning first), now find morning last */
event_type = SE_MORNING_LAST;
retflag = swe_heliacal_ut(dret[0], geopos, datm, dobs, obj_name, event_type, helflag, dret, serr);
if (retflag == ERR) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(dret[0], geopos, &(dret[0]), serr);
retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
}
do_print_heliacal(dret, event_type, obj_name);
tsave1 = dret[0];
/* mercury can have several evening appearances without any morning
* appearances in betweeen. We have to find out when the next
* morning appearance is and then find all evening appearances
* that take place before that */
if (ipl == SE_MERCURY) {
event_type = SE_HELIACAL_RISING;
retflag = swe_heliacal_ut(dret[0], geopos, datm, dobs, obj_name, event_type, helflag, dret, serr);
if (retflag == ERR) {
do_printf(serr);
return ERR;
}
tsave2 = dret[0];
}
repeat_mercury:
/* evening first */
event_type = SE_EVENING_FIRST;
retflag = swe_heliacal_ut(tsave1, geopos, datm, dobs, obj_name, event_type, helflag, dret, serr);
if (retflag == ERR) {
do_printf(serr);
return ERR;
}
if (ipl == SE_MERCURY && dret[0] > tsave2)
continue;
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(dret[0], geopos, &(dret[0]), serr);
retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
}
do_print_heliacal(dret, event_type, obj_name);
}
if (ipl == SE_MOON) {
/* morning last */
event_type = SE_MORNING_LAST;
retflag = swe_heliacal_ut(dret[0], geopos, datm, dobs, obj_name, event_type, helflag, dret, serr);
if (retflag == ERR) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(dret[0], geopos, &(dret[0]), serr);
retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
}
do_print_heliacal(dret, event_type, obj_name);
} else {
/* heliacal setting (evening last) */
event_type = SE_HELIACAL_SETTING;
retflag = swe_heliacal_ut(dret[0], geopos, datm, dobs, obj_name, event_type, helflag, dret, serr);
if (retflag == ERR) {
do_printf(serr);
return ERR;
}
if (time_flag & (BIT_TIME_LMT | BIT_TIME_LAT)) {
retc = ut_to_lmt_lat(dret[0], geopos, &(dret[0]), serr);
retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
}
do_print_heliacal(dret, event_type, obj_name);
if (0 && ipl == SE_MERCURY) {
tsave1 = dret[0];
goto repeat_mercury;
}
}
}
}
return OK;
}
static int do_special_event(double tjd, int32 ipl, char *star, int32 special_event, int32 special_mode, double *geopos, double *datm, double *dobs, char *serr)
{
int retc;
/* risings, settings, meridian transits */
if (special_event == SP_RISE_SET ||
special_event == SP_MERIDIAN_TRANSIT)
retc = call_rise_set(tjd, ipl, star, whicheph, special_mode, geopos, serr);
/* lunar eclipses */
if (special_event == SP_LUNAR_ECLIPSE)
retc = call_lunar_eclipse(tjd, whicheph, special_mode, geopos, serr);
/* solar eclipses */
if (special_event == SP_SOLAR_ECLIPSE)
retc = call_solar_eclipse(tjd, whicheph, special_mode, geopos, serr);
/* occultations by the moon */
if (special_event == SP_OCCULTATION)
retc = call_lunar_occultation(tjd, ipl, star, whicheph, special_mode, geopos, serr);
/* heliacal event */
if (special_event == SP_HELIACAL)
retc = call_heliacal_event(tjd, ipl, star, whicheph, special_mode, geopos, datm, dobs, serr);
return retc;
}
static char *hms_from_tjd(double x)
{
static char s[AS_MAXCH];
sprintf(s, "%s ", hms(fmod(x + 1000000.5, 1) * 24, BIT_LZEROES));
return s;
}
static char *hms(double x, int32 iflag)
{
static char s[AS_MAXCH], s2[AS_MAXCH], *sp;
char *c = ODEGREE_STRING;
x += 0.5 / 36000.0; /* round to 0.1 sec */
strcpy(s, dms(x, iflag));
sp = strstr(s, c);
if (sp != NULL) {
*sp = ':';
strcpy(s2, sp + strlen(ODEGREE_STRING));
strcpy(sp + 1, s2);
*(sp + 3) = ':';
*(sp + 8) = '\0';
}
return s;
}
static void do_printf(char *info)
{
#ifdef _WINDOWS
fprintf(fp, info);
#else
fputs(info,stdout);
#endif
}
/* make_ephemeris_path().
* ephemeris path includes
* current working directory
* + program directory
* + default path from swephexp.h on current drive
* + on program drive
* + on drive C:
*/
static int make_ephemeris_path(int32 iflg, char *argv0, char *path)
{
char *sp;
char *dirglue = DIR_GLUE;
size_t pathlen = 0;
/* moshier needs no ephemeris path */
if (iflg & SEFLG_MOSEPH)
return OK;
/* current working directory */
sprintf(path, ".%c", *PATH_SEPARATOR);
/* program directory */
sp = strrchr(argv0, *dirglue);
if (sp != NULL) {
pathlen = sp - argv0;
if (strlen(path) + pathlen < AS_MAXCH-2) {
strncat(path, argv0, pathlen);
sprintf(path + strlen(path), "%c", *PATH_SEPARATOR);
}
}
#if MSDOS
{
char *cpos[20];
char s[2 * AS_MAXCH], *s1 = s + AS_MAXCH;
char *sp[3];
int i, j, np;
strcpy(s1, SE_EPHE_PATH);
np = cut_str_any(s1, PATH_SEPARATOR, cpos, 20);
/*
* default path from swephexp.h
* - current drive
* - program drive
* - drive C
*/
*s = '\0';
/* current working drive */
sp[0] = getcwd(NULL, 0);
if (sp[0] == NULL) {
printf("error in getcwd()\n");
exit(1);
}
if (*sp[0] == 'C')
sp[0] = NULL;
/* program drive */
if (*argv0 != 'C' && (sp[0] == NULL || *sp[0] != *argv0))
sp[1] = argv0;
else
sp[1] = NULL;
/* drive C */
sp[2] = "C";
for (i = 0; i < np; i++) {
strcpy(s, cpos[i]);
if (*s == '.') /* current directory */
continue;
if (s[1] == ':') /* drive already there */
continue;
for (j = 0; j < 3; j++) {
if (sp[j] != NULL && strlen(path) + 2 + strlen(s) < AS_MAXCH-1) {
sprintf(path + strlen(path), "%c:%s%c", *sp[j], s, *PATH_SEPARATOR);
}
}
}
}
#else
if (strlen(path) + strlen(SE_EPHE_PATH) < AS_MAXCH-1)
strcat(path, SE_EPHE_PATH);
#endif
return OK;
}
static void remove_whitespace(char *s)
{
char *sp, s1[AS_MAXCH];
if (s == NULL) return;
for (sp = s; *sp == ' '; sp++)
;
strcpy(s1, sp);
while (*(sp = s1 + strlen(s1) - 1) == ' ')
*sp = '\0';
strcpy(s, s1);
}
#if MSDOS
/**************************************************************
cut the string s at any char in cutlist; put pointers to partial strings
into cpos[0..n-1], return number of partial strings;
if less than nmax fields are found, the first empty pointer is
set to NULL.
More than one character of cutlist in direct sequence count as one
separator only! cut_str_any("word,,,word2",","..) cuts only two parts,
cpos[0] = "word" and cpos[1] = "word2".
If more than nmax fields are found, nmax is returned and the
last field nmax-1 rmains un-cut.
**************************************************************/
static int cut_str_any(char *s, char *cutlist, char *cpos[], int nmax)
{
int n = 1;
cpos [0] = s;
while (*s != '\0') {
if ((strchr(cutlist, (int) *s) != NULL) && n < nmax) {
*s = '\0';
while (*(s + 1) != '\0' && strchr (cutlist, (int) *(s + 1)) != NULL) s++;
cpos[n++] = s + 1;
}
if (*s == '\n' || *s == '\r') { /* treat nl or cr like end of string */
*s = '\0';
break;
}
s++;
}
if (n < nmax) cpos[n] = NULL;
return (n);
} /* cutstr */
#endif
Reference in New Issue View Git Blame Copy Permalink