@swisseph/node 1.0.1 → 1.0.2

package/libswe/swetest.c

@@ -0,0 +1,4079 @@
+/* 
+  swetest.c	A test program
+   
+  Authors: Dieter Koch and Alois Treindl, Astrodienst Zuerich
+
+**************************************************************/
+
+/* Copyright (C) 1997 - 2021 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 Affero General Public License (AGPL)
+  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 AGPL 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 AGPL or a compatible license.
+  See https://www.gnu.org/licenses/agpl-3.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:                                                                                   \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\
+        -glp  report file location of library\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[:MM[:SS]]  input time (as Ephemeris Time)\n\
+        -ut     input date is Universal Time (UT1)\n\
+	-utHH[:MM[:SS]] input time (as Universal Time)\n\
+	-utcHH[:MM[:SS]] input time (as Universal Time Coordinated UTC)\n\
+		H,M,S can have one or two digits. Their limits are unchecked.\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\
+        -dhX    differential ephemeris: print differential ephemeris between\n\
+                heliocentric body X and each body in list given by -p\n\
+                example: -p8 -dh8 -ftl -n36600 -b1.1.1500 -s5 prints the longitude\n\
+                distance between geocentric and heliocentric Neptune (planet 8)\n\
+                for 500 year starting at 1 Jan 1500.\n\
+		Using this option mostly makes sense for a single planet\n\
+		to find out how much its geocentric and heliocentric positions can differ\n\
+		over extended periods of time\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 timesteps; 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\
+                If an 'y' is appended, the time step is in years instead of days, \n\
+                for example -s10y for a time step of 10 years.\n\
+                If an 'mo' is appended, the time step is in months instead of days, \n\
+                for example -s3mo for a time step of 3 months.\n\
+                If an 'm' is appended, the time step is in minutes instead of days, \n\
+                for example -s15m for a time step of 15 minutes.\n\
+                If an 's' is appended, the time step is in seconds instead of days, \n\
+                for example -s1s for a time step of 1 second.\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\
+        -hor	list data for multiple planets 'horizontally' in same line.\n\
+		all columns of -fSEQ are repeated except time colums tTJyY.\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, Vertex,\n\
+		Equatorial Ascendant, co-Ascendant as defined by Walter Koch, \n\
+		co-Ascendant as defined by Michael Munkasey, and Polar Ascendant. \n\
+		Houses can only be computed if option -ut is given.\n\
+                   A  equal\n\
+                   B  Alcabitius\n\
+                   C  Campanus\n\
+                   D  equal / MC\n\
+                   E  equal = A\n\
+                   F  Carter poli-equatorial\n\
+                   G  36 Gauquelin sectors\n\
+                   H  horizon / azimuth\n\
+                   I  Sunshine\n\
+                   i  Sunshine alternative\n\
+                   K  Koch\n\
+                   L  Pullen S-delta\n\
+                   M  Morinus\n\
+                   N  Whole sign, Aries = 1st house\n\
+                   O  Porphyry\n\
+                   P  Placidus\n\
+                   Q  Pullen S-ratio\n\
+                   R  Regiomontanus\n\
+                   S  Sripati\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\
+		 The use of lower case letters is deprecated. They will have a\n\
+		 different meaning in future releases of Swiss Ephemeris.\n\
+        -hsy[hsys]	\n\
+		 house system to be used (for house positions of planets)\n\
+		 for long, lat, hsys, see -house\n\
+		 The use of lower case letters is deprecated. They will have a\n\
+		 different meaning in future releases of Swiss Ephemeris.\n\
+";
+static char *infocmd3 = "\
+        -geopos[long,lat,elev]	\n\
+		Geographic position. Can be used for azimuth and altitude\n\
+                or house cusps 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\
+		For topocentric planet positions please user the parameter -topo\n\
+     sidereal astrology:\n\
+	-ay..   ayanamsha, with number of method, e.g. ay0 for Fagan/Bradley\n\
+	-sid..    sidereal, with number of method (see below)\n\
+	-sidt0..  dito, but planets are projected on the ecliptic plane of the\n\
+	          reference date of the ayanamsha (more info in general documentation\n\
+		  www.astro.com/swisseph/swisseph.htm)\n\
+	-sidsp..  dito, but planets are projected on the solar system plane.\n\
+		  (see www.astro.com/swisseph/swisseph.htm)\n\
+        -sidudef[jd,ay0,...]  sidereal, with user defined ayanamsha; \n\
+	          jd=julian day number in TT/ET\n\
+	          ay0=initial value of ayanamsha, \n\
+		  ...=optional parameters, comma-sparated:\n\
+		  'jdisut': ayanamsha reference date is UT\n\
+		  'eclt0':  project on ecliptic of reference date (like -sidt0..)\n\
+		  'ssyplane':  project on solar system plane (like -sidsp..)\n\
+		  e.g. '-sidudef2452163.8333333,25.0,jdisut': ayanamsha is 25.0° on JD 2452163.8333333 UT\n\
+           number of ayanamsha method:\n\
+	   0 for Fagan/Bradley\n\
+	   1 for Lahiri\n\
+	   2 for De Luce\n\
+	   3 for Raman\n\
+	   4 for Usha/Shashi\n\
+	   5 for Krishnamurti\n\
+	   6 for Djwhal Khul\n\
+	   7 for Yukteshwar\n\
+	   8 for J.N. Bhasin\n\
+	   9 for Babylonian/Kugler 1\n\
+	   10 for Babylonian/Kugler 2\n\
+	   11 for Babylonian/Kugler 3\n\
+	   12 for Babylonian/Huber\n\
+	   13 for Babylonian/Eta Piscium\n\
+	   14 for Babylonian/Aldebaran = 15 Tau\n\
+	   15 for Hipparchos\n\
+	   16 for Sassanian\n\
+	   17 for Galact. Center = 0 Sag\n\
+	   18 for J2000\n\
+	   19 for J1900\n\
+	   20 for B1950\n\
+	   21 for Suryasiddhanta\n\
+	   22 for Suryasiddhanta, mean Sun\n\
+	   23 for Aryabhata\n\
+	   24 for Aryabhata, mean Sun\n\
+	   25 for SS Revati\n\
+	   26 for SS Citra\n\
+	   27 for True Citra\n\
+	   28 for True Revati\n\
+	   29 for True Pushya (PVRN Rao)\n\
+	   30 for Galactic (Gil Brand)\n\
+	   31 for Galactic Equator (IAU1958)\n\
+	   32 for Galactic Equator\n\
+	   33 for Galactic Equator mid-Mula\n\
+	   34 for Skydram (Mardyks)\n\
+	   35 for True Mula (Chandra Hari)\n\
+	   36 Dhruva/Gal.Center/Mula (Wilhelm)\n\
+	   37 Aryabhata 522\n\
+	   38 Babylonian/Britton\n\
+   	   39 Vedic/Sheoran\n\
+	   40 Cochrane (Gal.Center = 0 Cap)\n\
+	   41 Galactic Equator (Fiorenza)\n\
+	   42 Vettius Valens\n\
+	   43 Lahiri 1940\n\
+	   44 Lahiri VP285 (1980)\n\
+	   45 Krishnamurti VP291\n\
+	   46 Lahiri ICRC\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\
+";
+static char *infocmd4 = "\
+        -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 more precise \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\
+	-ep		  use extra precision in output for some data\n\
+	-dms              use dms instead of fractions, at some places\n\
+	-lim		  print ephemeris file range\n\
+     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\
+		Greenwich is used 0.00,51.50,0\n\
+        -pc...  compute planetocentric positions\n\
+                to specify the central body, use the internal object number\n\
+		of Swiss Ephemeris, e.g. 3 for Venus, 4 for Mars, \n\
+        -pc3 	Venus-centric \n\
+        -pc4 	Mars-centric \n\
+        -pc5 	Jupiter-centric (barycenter)\n\
+	-pc9599 Jupiter-centric (center of body)\n\
+	-pc9699 Saturn-centric (center of body)\n\
+		For asteroids use MPC number + 10000, e.g.\n\
+	-pc10433 Eros-centric (Eros = 433 + 10000)\n\
+     orbital elements:\n\
+        -orbel  compute osculating orbital elements relative to the\n\
+	        mean ecliptic J2000. (Note, all values, including time of\n\
+		pericenter vary considerably depending on the date for which the\n\
+		osculating ellipse is calculated\n\
+\n\
+     special events:\n\
+        -solecl solar eclipse\n\
+                output 1st line:\n\
+                  eclipse date,\n\
+                  time of maximum (UT):\n\
+		    geocentric angle between centre of Sun and Moon reaches minimum.\n\
+                  core shadow width (negative with total eclipses),\n\
+		  eclipse magnitudes:\n\
+		    1. NASA method (= 2. with partial ecl. and \n\
+		       ratio lunar/solar diameter with total and annular ecl.)\n\
+		    2. fraction of solar diameter covered by moon;\n\
+		       if the value is > 1, it means that Moon covers more than\n\
+		       just the solar disk\n\
+		    3. fraction of solar disc covered by moon (obscuration)\n\
+		       with total and annular eclipses it is the ratio of\n\
+		       the sizes of the solar disk and the lunar disk.\n\
+		  Saros series and eclipse number\n\
+		  Julian day number (6-digit fraction) of maximum\n\
+                output 2nd line:\n\
+                  start and end times for partial and total phases\n\
+		  delta t in sec\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, with the following differences:\n\
+		  Magnitude is defined like no. 2. with solar eclipses.\n\
+		  There are no saros series.\n\
+";
+static char *infocmd5 = "\
+        -lunecl lunar eclipse\n\
+                output 1st line:\n\
+                  eclipse date,\n\
+                  time of maximum (UT),\n\
+                  eclipse magnitudes: umbral and penumbral\n\
+		    method as method 2 with solar eclipses\n\
+		  Saros series and eclipse number \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\
+		  delta t in sec\n\
+                output 3rd line:\n\
+                  geographic position where the Moon is in zenith at maximum eclipse\n\
+        -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\
+                  eclipse magnitudes, as with global solar eclipse function \n\
+		    (with occultations: only diameter method, see solar eclipses, method 2)\n\
+		  Saros series and eclipse number (with solar eclipses only)\n\
+		  Julian day number (6-digit fraction) of maximum\n\
+                output 2nd line:\n\
+                  local eclipse duration for totality (zero with partial occultations)\n\
+                  local four contacts,\n\
+		  delta t in sec\n\
+		Occultations with the remark \"(daytime)\" cannot be observed because\n\
+		they are taking place by daylight. Occultations with the remark\n\
+		\"(sunrise)\" or \"(sunset)\" can be observed only partly because part\n\
+		of them takes place in daylight.\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\
+";
+static char *infocmd6 = "\
+     specifications for risings and settings:\n\
+        -norefrac   neglect refraction (with option -rise)\n\
+        -disccenter find rise of disc center (with option -rise)\n\
+        -discbottom find rise of disc bottom (with option -rise)\n\
+	-hindu      hindu version of sunrise (with option -rise)\n\
+     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:
+  ijklruv
+ */
+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 body numbers/letters:\n\
+        0 Sun (character zero)\n\
+        1 Moon (character 1)\n\
+        2 Mercury\n\
+        3 Venus\n\
+        4 Mars\n\
+        5 Jupiter\n\
+        6 Saturn\n\
+        7 Uranus\n\
+        8 Neptune\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\
+	b ayanamsha\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\
+        For planets Jupiter to Pluto the center of body (COB) can be\n\
+        calculated using the additional parameter -cob\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\
+     some planetary moons and center of body of a planet:\n\
+        v with moon number given in -xv:\n\
+        v -xv9501 Io/Jupiter:\n\
+        v -xv9599 Jupiter, center of body (COB):\n\
+        v -xv94.. Mars moons:\n\
+        v -xv95.. Jupiter moons and COB:\n\
+        v -xv96.. Saturn moons and COB:\n\
+        v -xv97.. Uranus moons and COB:\n\
+        v -xv98.. Neptune moons and COB:\n\
+        v -xv99.. Pluto moons and COB:\n\
+          The numbers of the moons are given here: \n\
+	  https://www.astro.com/ftp/swisseph/ephe/sat/plmolist.txt\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\
+     sidereal time:\n\
+        x sidereal time\n\
+        e print a line of labels\n\
+          \n";
+/* characters still available 
+   CcEeMmOoqWwz
+*/
+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\
+        W distance decimal in light years\n\
+        w distance decimal in km\n\
+        q relative distance (1000=nearest, 0=furthest)\n\
+        A right ascension in hh:mm:ss\n\
+        a right ascension hours decimal\n\
+	m Meridian distance \n\
+	z Zenith distance \n\
+        D declination degree\n\
+        d declination decimal\n\
+        I azimuth degree\n\
+        i azimuth decimal\n\
+        H altitude degree\n\
+        h altitude decimal\n\
+        K altitude (with refraction) degree\n\
+        k altitude (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): perihelion, aphelion, second focal point; longitude dec.\n\
+	F apsides (osc.): perihelion, aphelion, 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 (0 h).\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\
+\n\
+    swetest -b1.1.2016  -g -fTlbR -p0123456789Dmte -hor -n366 -roundsec\n\
+	tabular ephemeris (all planets Sun - Pluto, Chiron, mean node, true node)\n\
+	in one horizontal row, tab-separated, for 366 days. For each planet\n\
+	list longitude, latitude and geocentric distance.\n";
+/**************************************************************/
+
+#include "swephexp.h" 	/* this includes  "sweodef.h" */
+#include "swephlib.h"
+#include "sweph.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
+#ifdef __MINGW32__
+#  include <sys/stat.h>
+#else
+#  include <sys\stat.h>
+#endif
+#  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 BIT_ALLOW_361    64
+
+#define PLSEL_D  "0123456789mtA"
+#define PLSEL_P  "0123456789mtABCcgDEFGHI"
+#define PLSEL_H  "JKLMNOPQRSTUVWXYZw"
+#define PLSEL_A  "0123456789mtABCcgDEFGHIJKLMNOPQRSTUVWXYZw"
+
+#define DIFF_DIFF	'd'
+#define DIFF_GEOHEL	'h'
+#define DIFF_MIDP	'D'
+#define MODE_HOUSE	1
+#define MODE_LABEL	2
+#define MODE_AYANAMSA	4
+
+#define SEARCH_RANGE_LUNAR_CYCLES 20000
+
+#define LEN_SOUT    1000 // length of output string variable
+#define SIND(x) sin((x) * DEGTORAD)
+#define COSD(x) cos((x) * DEGTORAD)
+#define ACOSD(x)        (acos((x)) * RADTODEG)
+
+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(char *argv0, char *ephepath);
+static int letter_to_ipl(int letter);
+static int print_line(int mode, AS_BOOL is_first, int sid_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 int32 orbital_elements(double tjd_et, int32 ipl, int32 iflag, 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
+static int32 call_swe_fixstar(char *star, double te, int32 iflag, double *x, char *serr);
+static void jd_to_time_string(double jut, char *stimeout);
+static char *our_strcpy(char *to, char *from);
+
+/* 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 int iplctr = SE_SUN;
+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 AS_BOOL gregflag_auto = TRUE;
+static int diff_mode = 0;
+static AS_BOOL use_dms = FALSE;
+static AS_BOOL universal_time = FALSE;
+static AS_BOOL universal_time_utc = FALSE;
+static int32 round_flag = 0;
+static int32 time_flag = 0;
+static AS_BOOL short_output = FALSE;
+static AS_BOOL list_hor = FALSE;
+static int32 special_event = 0;
+static int32 special_mode = 0;
+static AS_BOOL do_orbital_elements = FALSE;
+static AS_BOOL hel_using_AV = FALSE;
+static AS_BOOL with_header = TRUE;
+static AS_BOOL with_chart_link = 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", "equat. Asc.","co-Asc. W.Koch", "co-Asc Munkasey", "Polar Asc."};
+static int direction = 1;
+static AS_BOOL direction_flag = FALSE;
+static AS_BOOL step_in_minutes = FALSE;
+static AS_BOOL step_in_seconds = FALSE;
+static AS_BOOL step_in_years = FALSE;
+static AS_BOOL step_in_months = FALSE;
+static int32 helflag = 0;
+static double tjd = 2415020.5;
+static int32 nstep = 1, istep;
+static int32 search_flag = 0;
+static char sout[LEN_SOUT];
+static int32 whicheph = SEFLG_SWIEPH;   
+static char *psp;
+static int32 norefrac = 0;
+static int32 disccenter = 0;
+static int32 discbottom = 0;
+static int32 hindu = 0;
+/* for test of old models only */
+static char *astro_models;
+static int do_set_astro_models = FALSE;
+static char smod[2000];
+static AS_BOOL inut = FALSE; /* for Astrodienst internal feature */
+static AS_BOOL have_gap_parameter = FALSE; 
+static AS_BOOL use_swe_fixstar2 = FALSE; 
+static AS_BOOL output_extra_prec = FALSE;
+static AS_BOOL show_file_limit = FALSE;
+
+#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 = 36525; // days in a century
+  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';
+  int year_start = 0, mon_start = 1, day_start = 1;
+  AS_BOOL do_houses = FALSE;
+  char ephepath[AS_MAXCH];
+  char fname[AS_MAXCH];
+  char sdate[AS_MAXCH];
+  char *begindate = NULL;
+  char stimein[AS_MAXCH];
+  char stimeout[AS_MAXCH];
+  int32 iflgret;
+  AS_BOOL is_first = TRUE;
+  AS_BOOL with_glp = FALSE;
+  AS_BOOL with_header_always = FALSE;
+  AS_BOOL do_ayanamsa = FALSE;
+  AS_BOOL do_planeto_centric = FALSE;
+  double aya_t0 = 0, aya_val0 = 0;
+  AS_BOOL no_speed = FALSE;
+  int32 sid_mode = SE_SIDM_FAGAN_BRADLEY;
+  double t2, tstep = 1, thour = 0;
+  double delt;
+  double tid_acc = 0;
+  datm[0] = 1013.25; datm[1] = 15; datm[2] = 40; datm[3] = 0;
+  dobs[0] = 0; dobs[1] = 0;
+  dobs[2] = 0; dobs[3] = 0; dobs[4] = 0; dobs[5] = 0;
+  serr[0] = serr_save[0] = serr_warn[0] = sdate_save[0] = '\0';
+# ifdef MACOS
+  argc = ccommand(&argv); /* display the arguments window */    
+# endif
+  *stimein = '\0';
+  strcpy(ephepath, "");
+  strcpy(fname, SE_FNAME_DFT);
+  for (i = 1; i < argc; i++) {
+    if (strncmp(argv[i], "-utc", 4) == 0) {
+      universal_time = TRUE;
+      universal_time_utc = TRUE;
+      if (strlen(argv[i]) > 4) {
+	strncpy(stimein, argv[i] + 4, 30);
+	stimein[30] = '\0';
+      }
+    } else if (strncmp(argv[i], "-ut", 3) == 0) {
+      universal_time = TRUE;
+      if (strlen(argv[i]) > 3) {
+	strncpy(stimein, argv[i] + 3, 30);
+	stimein[30] = '\0';
+      }
+    } else if (strncmp(argv[i], "-glp", 4) == 0) {
+      with_glp = TRUE;
+    } else if (strncmp(argv[i], "-hor", 4) == 0) {
+      list_hor = TRUE;
+    } 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 (strcmp(argv[i], "-cob") == 0) {
+      iflag |= SEFLG_CENTER_BODY;
+    } 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;
+    } else if (strncmp(argv[i], "-sidudef", 8) == 0) {
+      iflag |= SEFLG_SIDEREAL;
+      sid_mode = SE_SIDM_USER;
+      strcpy(s1, argv[i] + 8);
+      aya_t0 = atof(s1);
+      if ((sp = strchr(s1, ',')) != NULL) {
+        aya_val0 = atof(sp+1);
+      }
+      if (strstr(sp, "jdisut") != NULL) {
+        sid_mode |= SE_SIDBIT_USER_UT;
+      }
+      /*swe_set_sid_mode(sid_mode, 0, 0);*/
+    } else if (strncmp(argv[i], "-sidbit", 7) == 0) {
+      sid_mode |= atoi(argv[i]+7);
+    } 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], "-tpm") == 0) {
+      iflag |= SEFLG_TEST_PLMOON;
+    } 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') {
+        strncpy(fname, argv[i]+5, AS_MAXCH - 1);
+        fname[AS_MAXCH-1] = '\0';
+      }
+    } else if (strncmp(argv[i], "-edir", 5) == 0) {
+      if (*(argv[i]+5) != '\0') {
+        strncpy(ephepath, argv[i]+5, AS_MAXCH - 1);
+	ephepath[AS_MAXCH-1] = '\0';
+      }
+    } 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';
+      sp = argv[i] + 6;
+      if (*sp == '[') sp++;
+      sscanf(sp, "%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;
+      sp = argv[i] + 5;
+      if (*sp == '[') sp++;
+      sscanf(sp, "%lf,%lf,%lf", &top_long, &top_lat, &top_elev);
+      have_geopos = TRUE;
+    } else if (strncmp(argv[i], "-geopos", 7) == 0) {
+      sp = argv[i] + 7;
+      if (*sp == '[') sp++;
+      sscanf(sp, "%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 (strcmp(argv[i], "-nospeed") == 0) {
+      no_speed = TRUE;
+    } 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 (strncmp(argv[i], "-lmt", 4) == 0) {
+      universal_time = TRUE;
+      time_flag |= BIT_TIME_LMT;
+      if (strlen(argv[i]) > 4) {
+	strncpy(stimein, argv[i] + 4, 30);
+	stimein[30] = '\0';
+      }
+    } else if (strcmp(argv[i], "-lat") == 0) {
+      universal_time = TRUE;
+      time_flag |= BIT_TIME_LAT;
+    } else if (strcmp(argv[i], "-lim") == 0) {
+      show_file_limit = TRUE;
+    } else if (strcmp(argv[i], "-clink") == 0) {
+      with_chart_link = TRUE;
+    } 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], "-ep") == 0) {
+      output_extra_prec = 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;
+    } else if (strcmp(argv[i], "-annular") == 0) {
+      search_flag |= SE_ECL_ANNULAR;
+    } else if (strcmp(argv[i], "-anntot") == 0) {
+      search_flag |= SE_ECL_ANNULAR_TOTAL;
+    } else if (strcmp(argv[i], "-partial") == 0) {
+      search_flag |= SE_ECL_PARTIAL;
+    } 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) {
+      hindu = 1;
+      norefrac = 1;
+      disccenter = 1;
+    } else if (strcmp(argv[i], "-discbottom") == 0) {
+      discbottom = 1;
+    } else if (strcmp(argv[i], "-metr") == 0) {
+      special_event = SP_MERIDIAN_TRANSIT;
+      have_geopos = TRUE;
+    /* undocumented test feature */
+    } else if (strncmp(argv[i], "-amod",5) == 0) {
+      astro_models = argv[i] + 5;
+      do_set_astro_models = TRUE;
+    /* undocumented test feature */
+    } else if (strncmp(argv[i], "-tidacc",7) == 0) {
+      tid_acc = atof(argv[i] + 7);
+    } else if (strncmp(argv[i], "-hev", 4) == 0) {
+      special_event = SP_HELIACAL;
+      search_flag = 0;
+      sp = argv[i] + 4;
+      if (*sp == '[') sp++;
+      if (strlen(sp) > 0)
+	search_flag = atoi(sp);
+      have_geopos = TRUE;
+      if (strstr(argv[i], "AV")) hel_using_AV = TRUE;
+    } else if (strncmp(argv[i], "-at", 3) == 0) {
+      sp = argv[i]+3;
+      if (*sp == '[') sp++;
+      j = 0;
+      while (j < 4 && sp != NULL) {
+        datm[j] = atof(sp);
+	sp = strchr(sp, ',');
+	if (sp != NULL) sp += 1;
+        j++;
+      }
+    } else if (strncmp(argv[i], "-obs", 4) == 0) {
+      sp = argv[i] + 4;
+      if (*sp == '[') sp++;
+      sscanf(sp, "%lf,%lf", &(dobs[0]), &(dobs[1]));
+    } else if (strncmp(argv[i], "-opt", 4) == 0) {
+      sp = argv[i] + 4;
+      if (*sp == '[') sp++;
+      sscanf(sp, "%lf,%lf,%lf,%lf,%lf,%lf", &(dobs[0]), &(dobs[1]), &(dobs[2]), &(dobs[3]), &(dobs[4]), &(dobs[5]));
+    } else if (strcmp(argv[i], "-orbel") == 0) {
+      do_orbital_elements = TRUE;
+    } else if (strcmp(argv[i], "-bwd") == 0) {
+      direction = -1;
+      direction_flag = TRUE;
+    } else if (strncmp(argv[i], "-pc", 3) == 0) {
+      iplctr = atoi(argv[i]+3);
+      do_planeto_centric = 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 */
+      strncpy(sastno, argv[i] + 3, AS_MAXCH - 1);
+      sastno[AS_MAXCH-1] = '\0';
+    } else if (strncmp(argv[i], "-xv", 3) == 0) {
+      /* number of planetary moon */
+      strncpy(sastno, argv[i] + 3, AS_MAXCH - 1);
+      sastno[AS_MAXCH-1] = '\0';
+    } else if (strncmp(argv[i], "-xf", 3) == 0) {
+      /* name or number of fixed star */
+      strncpy(star, argv[i] + 3, AS_MAXCH - 1);
+      star[AS_MAXCH-1] = '\0';
+    } else if (strncmp(argv[i], "-xz", 3) == 0) {
+      /* number of hypothetical body */
+      strncpy(shyp, argv[i] + 3, AS_MAXCH - 1);
+      shyp[AS_MAXCH-1] = '\0';
+    } else if (strncmp(argv[i], "-x", 2) == 0) {
+      /* name or number of fixed star */
+      strncpy(star, argv[i] + 2, AS_MAXCH - 1);
+      star[AS_MAXCH-1] = '\0';
+    } else if (strcmp(argv[i], "-nut") == 0) {
+      inut = TRUE;
+    } 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 (strcmp(argv[i], "-swefixstar2") == 0) {
+      use_swe_fixstar2 = TRUE; 
+    } else if (strncmp(argv[i], "-s", 2) == 0) {
+      tstep = atof(argv[i]+2);
+      if (*(argv[i] + strlen(argv[i]) -1) == 'm')
+        step_in_minutes = TRUE;
+      if (*(argv[i] + strlen(argv[i]) -1) == 's')
+        step_in_seconds = TRUE;
+      if (*(argv[i] + strlen(argv[i]) -1) == 'y')
+        step_in_years = TRUE;
+      if (*(argv[i] + strlen(argv[i]) -1) == 'o') {
+        step_in_minutes = FALSE;
+        step_in_months = TRUE;
+      }
+    } 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;
+      have_gap_parameter = TRUE;
+      if (*gap == '\0') gap = "\t";
+    } else if (strcmp(argv[i], "-dms") == 0) {
+      use_dms = TRUE;
+    } 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;
+      if (*(argv[i] + 2) == 'h') {
+        sp++;
+        diff_mode = 'h';   // diff helio to geo
+      }
+      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], "-timeout", 8) == 0) {
+      swe_set_timeout(atoi(argv[i]) + 8);*/
+    } else if (strncmp(argv[i], "-t", 2) == 0) {
+      if (strlen(argv[i]) > 2) {
+	strncat(stimein, argv[i] + 2, 30);
+	stimein[30] = '\0';
+      }
+    } else if (strncmp(argv[i], "-h", 2) == 0
+      || strncmp(argv[i], "-?", 2) == 0) {
+      sp = argv[i]+2;
+      if (*sp == 'c' || *sp == '\0') {
+        //char si0[strlen(infocmd0)+1]; // Microsoft Visual Studio does not like this
+        char si0[2000];  
+	swe_version(sout);
+	strcpy(si0, infocmd0);
+        sp2 = strstr(si0, "Version:");
+	if (sp2 != NULL && strlen(sp2) > 10 + strlen(sout)) 
+	  strcpy(sp2 + 9, sout);
+        fputs(si0,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], 100);
+      sout[100] = '\0';
+      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 {
+      if (*plsel == 's')
+	ipl = atoi(sastno) + SE_AST_OFFSET;
+      *star = '\0';
+    }
+    if (special_event == SP_OCCULTATION && ipl == 1)
+      ipl = 2; /* no occultation of moon by moon */
+  }
+  if (*stimein != '\0') {
+    t = 0;
+    if ((sp = strchr(stimein, ':')) != NULL) {
+      if ((sp2 = strchr(sp + 1, ':')) != NULL) {
+        t += atof(sp2 + 1) / 60.0;
+      }
+      t += atoi(sp + 1);
+      t /= 60.0;
+    }
+    if (atoi(stimein) < 0)
+      t = -t;
+    t += atoi(stimein);
+    //t += 0.0000000001;
+    thour = t;
+  }
+#if HPUNIX
+  gethostname (hostname, 80);
+  if (strstr(hostname, "as10") != NULL) 
+    line_limit = 1000;
+#endif
+#if MSDOS
+  SetConsoleOutputCP(65001);	// set console to utf-8,
+  				// works only from Windows Vista upwards, not on XP.
+#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_SPEED3) && !no_speed) 
+    iflag |= SEFLG_SPEED;
+  if (*ephepath == '\0') {
+    if (make_ephemeris_path(argv[0], ephepath) == ERR) {
+      iflag = (iflag & ~SEFLG_EPHMASK) | SEFLG_MOSEPH;
+      whicheph = SEFLG_MOSEPH;
+    }
+  }
+  if (whicheph != SEFLG_MOSEPH) 
+    swe_set_ephe_path(ephepath);
+  if (whicheph & SEFLG_JPLEPH)
+    swe_set_jpl_file(fname);
+  /* the following is only a test feature */
+  if (do_set_astro_models) {
+    swe_set_astro_models(astro_models, iflag); /* secret test feature for dieter */
+    swe_get_astro_models(astro_models, smod, iflag);
+  }
+#if 1
+  if (inut) /* Astrodienst internal feature */
+    swe_set_interpolate_nut(TRUE);
+#endif
+  if ((iflag & SEFLG_SIDEREAL) || do_ayanamsa) {
+    if (sid_mode & SE_SIDM_USER) 
+      swe_set_sid_mode(sid_mode, aya_t0, aya_val0);
+    else
+      swe_set_sid_mode(sid_mode, 0, 0);
+  }
+  geopos[0] = top_long;
+  geopos[1] = top_lat;
+  geopos[2] = top_elev;
+  swe_set_topo(top_long, top_lat, top_elev);
+  if (tid_acc != 0) 
+    swe_set_tid_acc(tid_acc); 
+  serr[0] = serr_save[0] = serr_warn[0] = '\0';
+  while (TRUE) {
+    if (begindate == NULL) {
+      printf("\nDate ?");
+      sdate[0] = '\0';
+      if( !fgets(sdate, AS_MAXCH, stdin) ) goto end_main;
+    } else {
+      strncpy(sdate, begindate, AS_MAXCH-1);
+      sdate[AS_MAXCH-1] = '\0';
+      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';
+    }
+    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;
+	gregflag_auto = FALSE;
+      } else if (strstr(sp, "greg") != NULL) {
+        gregflag = SE_GREG_CAL;
+	gregflag_auto = FALSE;
+      }
+      swe_revjul(tjd, gregflag, &jyear, &jmon, &jday, &jut);
+      year_start = jyear;
+      mon_start = jmon;
+      day_start = jday;
+    } 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);
+      year_start = jyear;
+      mon_start = jmon;
+      day_start = jday;
+      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;
+	gregflag_auto = FALSE;
+      } else if (strstr(sp, "greg") != NULL) {
+        gregflag = SE_GREG_CAL;
+	gregflag_auto = FALSE;
+      }
+      jut = 0;
+      if (universal_time_utc) {
+	int ih = 0, im = 0;
+	double ds = 0.0;
+	if (*stimein != '\0') {
+	  sscanf(stimein, "%d:%d:%lf", &ih, &im, &ds);
+	}
+        if (swe_utc_to_jd(jyear,jmon,jday, ih, im, ds, gregflag, tret, serr) == ERR) {
+	  printf(" error in swe_utc_to_jd(): %s\n", serr);
+	  exit(-1);
+	}
+	tjd = tret[1];
+      } else {
+	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 (step_in_minutes) 
+        t = tjd + (istep -1) * tstep / 1440;
+      if (step_in_seconds) 
+        t = tjd + (istep -1) * tstep / 86400;
+      if (step_in_years) {
+	t = swe_julday(year_start + (istep - 1) * (int) tstep, mon_start, day_start, jut, gregflag);
+      }
+      if (step_in_months) {
+	jmon = mon_start + (istep - 1) * (int) tstep;
+	jyear = year_start + (int) ((jmon - 1) / 12);
+	jmon = ((jmon - 1) % 12) + 1;
+	t = swe_julday(jyear, jmon, day_start, jut, gregflag);
+      }
+      if (gregflag_auto) {
+	if (t < 2299160.5)
+	  gregflag = SE_JUL_CAL;
+	else
+	  gregflag = SE_GREG_CAL;
+      }
+      // must repeat because gregflag may have changed
+      if (step_in_years) {
+	t = swe_julday(year_start + (istep - 1) * (int) tstep, mon_start, day_start, jut, gregflag);
+      }
+      if (step_in_months) {
+	jmon = mon_start + (istep - 1) * (int) tstep;
+	jyear = year_start + (int) ((jmon - 1) / 12);
+	jmon = ((jmon - 1) % 12) + 1;
+	t = swe_julday(jyear, jmon, day_start, jut, gregflag);
+      }
+      delt = swe_deltat_ex(t, iflag, serr);
+      if (!universal_time) {
+	delt = swe_deltat_ex(t - delt, iflag, serr);
+      }
+      t2 = t;
+      // output line: 
+      // "date (dmy) 4.6.2017 greg.   2:07:00 TT		version 2.07.02"
+      swe_revjul(t2, gregflag, &jyear, &jmon, &jday, &jut);
+      if (with_header) {
+#ifndef NO_SWE_GLP      // -DNO_SWE_GLP to suppress this function
+	if (with_glp) {
+	  swe_get_library_path(sout);
+	  printf("\npath: %s", sout);
+	}
+#endif
+        printf("\ndate (dmy) %d.%d.%04d", jday, jmon, jyear);
+        if (gregflag)
+          printf(" greg.");
+        else
+          printf(" jul.");
+	jd_to_time_string(jut, stimeout);
+	printf(stimeout);
+        if (universal_time) {
+	  if (time_flag & BIT_TIME_LMT)
+	    printf(" LMT");
+	  else
+	    printf(" UT");
+        } else {
+          printf(" TT");
+	}
+	printf("\t\tversion %s", swe_version(sout));
+      }
+      if (universal_time) {
+	// "LMT: 2457908.588194444"
+	if (time_flag & BIT_TIME_LMT) {
+	  if (with_header) {
+	    printf("\nLMT: %.9f", t);
+	    t -= geopos[0] / 15.0 / 24.0;   
+	  }
+	}
+	// "UT:  2457908.565972222     delta t: 68.761612 sec"
+        if (with_header) {
+          printf("\nUT:  %.9f", t);
+          printf("     delta t: %f sec", delt * 86400.0);
+        }
+        te = t + delt;
+        tut = t;
+      } else {
+        te = t;
+        tut = t - delt;
+	// "UT:  2457908.565972222     delta t: 68.761612 sec"
+        if (with_header) {
+          printf("\nUT:  %.9f", tut);
+          printf("     delta t: %f sec", delt * 86400.0);
+        }
+      }
+      iflgret = swe_calc(te, SE_ECL_NUT, iflag, xobl, serr);
+      if (with_header) {
+	// "TT:  2457908.566768074
+        printf("\nTT:  %.9f", te);
+	// "ayanamsa =   24° 5'51.6509 (Lahiri)"
+	if (iflag & SEFLG_SIDEREAL) {
+	  if (swe_get_ayanamsa_ex(te, iflag, &daya, serr) == ERR) {
+	    printf("   error in swe_get_ayanamsa_ex(): %s\n", serr);
+	    exit(1);
+	  }
+	  printf("   ayanamsa = %s (%s)", dms(daya, round_flag), swe_get_ayanamsa_name(sid_mode));
+	}
+	// "geo. long 8.000000, lat 47.000000, alt 0.000000"
+	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) {
+	  if (iflag & (SEFLG_NONUT | SEFLG_SIDEREAL)) {
+	    printf("\n%-15s %s", "Epsilon (m)", dms(xobl[0],round_flag)); 
+	  } else {
+	    printf("\n%-15s %s%s", "Epsilon (t/m)", dms(xobl[0],round_flag), gap);
+	    printf("%s", dms(xobl[1],round_flag)); 
+	  }
+	}
+	if (strchr(plsel, 'n') == NULL && !(iflag & (SEFLG_NONUT | SEFLG_SIDEREAL))) {
+	  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) {
+	  const 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) {
+	if (swe_get_ayanamsa_ex(te, iflag, &daya, serr) == ERR) {
+	  printf("   error in swe_get_ayanamsa_ex(): %s\n", serr);
+	  exit(1);
+	}
+	x[0] = daya;
+        print_line(MODE_AYANAMSA, TRUE, sid_mode);
+	continue;
+      }
+      if (t == tjd && strchr(plsel, 'e')) {
+	if (list_hor) {
+	  is_first = TRUE;
+	  for (psp = plsel; *psp != '\0'; psp++) {
+	    if (*psp == 'e') continue;
+	    ipl = letter_to_ipl((int) *psp);
+	    *spnam = '\0';
+	    if (ipl >= SE_SUN && ipl <= SE_VESTA)
+	      swe_get_planet_name(ipl, spnam);
+	    print_line(MODE_LABEL, is_first, 0);
+	    is_first = FALSE;
+	  }
+	  printf("\n");
+	} else {
+	  print_line(MODE_LABEL, TRUE, 0);
+        }
+      }
+      is_first = TRUE;
+      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')      // fixed star
+          ipl = SE_FIXSTAR;
+        else if (*psp == 's') // asteroid
+          ipl = atoi(sastno) + 10000;
+        else if (*psp == 'v') // planetary moon
+          ipl = atoi(sastno);
+        else if (*psp == 'z') // fictitious object
+          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 && !do_orbital_elements)
+            continue;
+	}
+        /* ecliptic position */
+	if (iflag_f >=0)
+	  iflag = iflag_f;
+        if (ipl == SE_FIXSTAR) {
+          iflgret = call_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 if (do_planeto_centric) {
+	  iflgret = swe_calc_pctr(te, ipl, iplctr, iflag, x, serr);
+	  swe_get_planet_name(ipl, se_pname);
+        } 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 (show_file_limit && ipl > SE_AST_OFFSET) {
+	    const char *fnam;
+	    char sbeg[40], send[40];
+	    double tfstart, tfend;
+            int denum;
+	    fnam = swe_get_current_file_data(3, &tfstart, &tfend, &denum);
+	    if (fnam != NULL) {
+	      swe_revjul(tfstart, gregflag, &jyear, &jmon, &jday, &jut);
+	      sprintf(sbeg, "%d.%02d.%04d", jday, jmon, jyear);
+	      swe_revjul(tfend, gregflag, &jyear, &jmon, &jday, &jut);
+	      sprintf(send, "%d.%02d.%04d", jday, jmon, jyear);
+	      printf("range %s: %.1lf = %s to %.1lf = %s de=%d\n", fnam, tfstart, sbeg, tfend, send, denum);
+	      show_file_limit = FALSE;
+	    }
+	  }
+        }
+	if (*psp == 'q') {/* delta t */
+	  x[0] = swe_deltat_ex(tut, iflag, serr) * 86400;
+	  x[1] = x[2] = x[3] = 0;
+	  x[1] = x[0] / 3600.0; // to hours
+	  strcpy(se_pname, "Delta T");
+	}
+	if (*psp == 'x') {/* sidereal time */
+	  x[0] = swe_degnorm(swe_sidtime(tut) * 15 + geopos[0]);
+	  x[1] = x[2] = x[3] = 0;
+	  strcpy(se_pname, "Sidereal Time");
+	}
+	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 (*psp == 'b') {/* ayanamsha */
+	  if (swe_get_ayanamsa_ex(te, iflag, &(x[0]), serr) == ERR) {
+	    printf("   error in swe_get_ayanamsa_ex(): %s\n", serr);
+	    iflgret = -1;
+	  }
+	  x[1] = 0;
+	  strcpy(se_pname, "Ayanamsha");
+	}
+        if (iflgret < 0) { 
+          if (strcmp(serr, serr_save) != 0 
+            && (ipl == SE_SUN || ipl == SE_MOON || ipl <= SE_PLUTO
+                || ipl == SE_MEAN_NODE || ipl == SE_TRUE_NODE
+                || ipl == SE_CERES || ipl == SE_PALLAS || ipl == SE_JUNO || ipl == SE_VESTA 
+                || ipl == SE_CHIRON || ipl == SE_PHOLUS || ipl == SE_CUPIDO 
+		|| ipl >= SE_PLMOON_OFFSET
+                || ipl >= SE_AST_OFFSET || ipl == SE_FIXSTAR
+		|| *psp == 'y')) {
+            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 (diff_mode == DIFF_GEOHEL)
+	    iflgret = swe_calc(te, ipldiff, iflag|SEFLG_HELCTR, x2, serr);
+          if (iflgret < 0) { 
+            fputs("error: ", stdout);
+            fputs(serr, stdout);
+            fputs("\n", stdout);
+          }
+	  if (diff_mode == DIFF_DIFF || diff_mode == DIFF_GEOHEL) {
+	    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, "aADdQmz") != NULL) {
+          iflag2 = iflag | SEFLG_EQUATORIAL;
+          if (ipl == SE_FIXSTAR) {
+            iflgret = call_swe_fixstar(star, te, iflag2, xequ, serr);
+	  } else if (do_planeto_centric) {
+	    iflgret = swe_calc_pctr(te, ipl, iplctr, 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 || diff_mode == DIFF_GEOHEL) {
+	      if (diff_mode == DIFF_GEOHEL)
+		iflgret = swe_calc(te, ipldiff, iflag2|SEFLG_HELCTR, x2, serr);
+	      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 = call_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, datm[0], datm[1], xt, xaz);
+          if (diff_mode) {
+            iflgret = swe_calc(te, ipldiff, iflgt, xt, serr);
+            swe_azalt(tut, SE_EQU2HOR, geopos, datm[0], datm[1], xt, x2);
+	    if (diff_mode == DIFF_DIFF || diff_mode == DIFF_GEOHEL) {
+	      if (diff_mode == DIFF_GEOHEL) {	// makes little sense for a heliocentric 
+		iflgret = swe_calc(te, ipldiff, iflgt|SEFLG_HELCTR, xt, serr);
+		swe_azalt(tut, SE_EQU2HOR, geopos, datm[0], datm[1], xt, x2);
+	      }
+	      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 = call_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 || diff_mode == DIFF_GEOHEL) {
+	      if (diff_mode == DIFF_GEOHEL)
+		iflgret = swe_calc(te, ipldiff, iflag2|SEFLG_HELCTR, x2, serr);
+	      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 = call_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 || diff_mode == DIFF_GEOHEL) {
+	      if (diff_mode == DIFF_GEOHEL)
+		iflgret = swe_calc(te, ipldiff, iflag2|SEFLG_HELCTR, x2, serr);
+	      for (i = 0; i < 6; i++) 
+		xcartq[i] -= x2[i];
+	    } else {
+		xcartq[i] = (xcart[i] + x2[i]) / 2;
+	    }
+          }
+        }
+        /* house position */
+        if (strpbrk(fmt, "gGjzm") != 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 && toupper(ihsy) == 'G') {
+	    swe_gauquelin_sector(tut, ipl, star2, iflag, hpos_meth, geopos, 0, 0, &hposj, serr);
+	  } else {
+	    double cusp[100];
+	    if (ihsy == 'i' || ihsy == 'I') {
+	      iflgret = swe_houses_ex(t,iflag, top_lat, top_long, ihsy, cusp, cusp+13);
+	    }
+	    hposj = swe_house_pos(armc, geopos[1], xobl[0], ihsy, xsv, serr);
+	  }
+	  if (toupper(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 (toupper(ihsy) == 'G')
+	      hpos2 = (hpos2 - 1) * 10;
+            else
+	      hpos2 = (hpos2 - 1) * 30;
+	    if (diff_mode == DIFF_DIFF || diff_mode == DIFF_GEOHEL) {
+	      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, is_first, 0);
+	is_first = FALSE;
+        if (! list_hor) line_count++;
+	if (do_orbital_elements) {
+	  orbital_elements(te, ipl, iflag, serr);
+	  continue;
+	}
+        if (line_count >= line_limit) {
+          printf("****** line count %d was exceeded\n", line_limit);
+          break;
+        }
+      }         /* for psp */
+      if (list_hor) {
+	printf("\n");
+        line_count++;
+      }
+      if (do_houses) {
+	double cusp[37];
+	double cusp_speed[37];
+	double ascmc[10];
+	double ascmc_speed[10];
+	int iofs; 
+	if (toupper(ihsy) == 'G') // Gauquelin has 36 cusps
+	  nhouses = 36;
+	iofs = nhouses + 1;
+	iflgret = swe_houses_ex2(t,iflag, top_lat, top_long, ihsy, cusp, ascmc, cusp_speed, ascmc_speed, serr);
+	// when swe_houses_ex() fails (e.g. with Placidus, Gauquelin, Makranski),
+	// it always returns Porphyry cusps instead
+        if (iflgret < 0) { 
+	  const char *shsy = swe_house_name(ihsy);
+	  sprintf(serr, "House method %s failed, Porphyry calculated instead", shsy); 
+          if (strcmp(serr, serr_save) != 0 ) {
+            fputs("error: ", stdout);
+            fputs(serr, stdout);
+            fputs("\n", stdout);
+          }
+          strcpy(serr_save, serr);
+	  ihsy = 'O';
+	  nhouses = 12; // instead of 36 with 'G'
+	  iofs = nhouses + 1;
+	} 
+	is_first = TRUE;
+	for (ipl = 1; ipl < iofs+8; ipl++) {
+	  x[0] = cusp[ipl];
+	  if (ipl >= iofs) {
+	    x[0] =  ascmc[ipl - iofs];
+	    x[3] =  ascmc_speed[ipl - iofs];
+	  } else {
+	    x[3] =  cusp_speed[ipl];
+	  }
+	  x[1] = 0;	/* latitude */
+	  x[2] = 1.0;	/* pseudo radius vector */
+	  if (ipl == iofs+2) { /* armc is already equatorial! */
+	    xequ[0] = x[0];
+	    xequ[1] = x[1];
+	    xequ[2] = x[2];
+	  } else if (strpbrk(fmt, "aADdQ") != NULL) {
+	    swe_cotrans(x, xequ, -xobl[0]);
+	  }
+	  if (strpbrk(fmt, "IiHhKk") != NULL) {
+	    double gpos[3];
+	    gpos[0] = top_long;
+	    gpos[1] = top_lat;
+	    gpos[2] = 0;
+	    swe_azalt(t, SE_ECL2HOR, gpos, datm[0], datm[1], x, xaz);
+	  }
+	  if (strpbrk(fmt, "gGj") != NULL) {
+	    hposj = swe_house_pos(armc, geopos[1], xobl[0], ihsy, x, serr);
+	    if (toupper(ihsy) == 'G')
+	      hpos = (hposj - 1) * 10;
+	    else
+	      hpos = (hposj - 1) * 30;
+	  }
+	  print_line(MODE_HOUSE, is_first, 0);
+	  is_first = FALSE;
+	  if (! list_hor) line_count++;
+	}
+	if (list_hor) {
+	  printf("\n");
+	  line_count++;
+	}
+      }
+      if (line_count >= line_limit) {
+	printf("****** line count %d was exceeded\n", 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:
+  if (do_set_astro_models) {
+    printf(smod);
+  }
+  swe_close();
+  return  OK;
+}
+
+static int32 call_swe_fixstar(char *star, double te, int32 iflag, double *x, char *serr)
+{
+  if (use_swe_fixstar2) 
+    return swe_fixstar2(star, te, iflag, x, serr);
+  else
+    return swe_fixstar(star, te, iflag, x, serr);
+}
+
+/* This function calculates the geocentric relative distance of a planet,
+ * where the closest position has value 1000, and remotest position has 
+ * value 0.
+ * The value is returned as an integer. The algorithm does not allow 
+ * much higher accuracy.
+ *
+ * With the Moon we measure the distance relative to the maximum and minimum
+ * found between 12000 BCE and 16000 CE.
+ * If the distance value were given relative to the momentary osculating 
+ * ellipse, then the apogee would always have the value 1000 and the perigee
+ * the value 0. It is certainly more interesting to know how much it is
+ * relative to a greater time range.
+ */
+static int32 get_geocentric_relative_distance(double tjd_et, int32 ipl, int32 iflag, char *serr)
+{
+  int32 iflagi = (iflag & (SEFLG_EPHMASK | SEFLG_HELCTR | SEFLG_BARYCTR));
+  int32 retval;
+  double ar = 0;
+  double xx[6];
+  double dmax, dmin, dtrue;
+  if ((0) && ipl == SE_MOON) {
+    dmax = 0.002718774; // jd = 283030.8
+    dmin = 0.002381834; // jd = -1006731.3
+    if ((retval = swe_calc(tjd_et, SE_MOON, iflagi | SEFLG_J2000 | SEFLG_TRUEPOS, xx, serr)) == ERR)
+      return 0;
+    dtrue = xx[2];
+  } else {
+    if (swe_orbit_max_min_true_distance(tjd_et, ipl, iflagi, &dmax, &dmin, &dtrue, serr) == ERR)
+      return 0;
+  }
+  if (dmax - dmin == 0) {
+    ar = 0;
+  } else {
+    ar = (1 - (dtrue - dmin) / (dmax - dmin)) * 1000.0;
+    ar += 0.5; // rounding
+  }
+  return (int32) ar;
+}
+
+/*
+ * The string fmt contains a sequence of format specifiers;
+ * each character in fmt creates a column, the columns are
+ * sparated by the gap string.
+ * Time columns tTJyY are only printed, if is_first is TRUE,
+ * so that they are not repeated in list_hor (horizontal list) mode.
+ * In list_hor mode, no newline is printed.
+ */
+static int print_line(int mode, AS_BOOL is_first, int sid_mode)
+{
+  char *sp, *sp2;
+  double t2, ju2 = 0;
+  double y_frac;
+  double ar, sinp;
+  double dret[20];
+  char slon[40];
+  char pnam[30];
+  AS_BOOL is_house = ((mode & MODE_HOUSE) != 0);
+  AS_BOOL is_label = ((mode & MODE_LABEL) != 0);
+  AS_BOOL is_ayana = ((mode & MODE_AYANAMSA) != 0);
+  int32 iflgret, dar;
+  // build planet name column, just in case
+  if (is_house) {
+    if (ipl <= nhouses) {
+      sprintf(pnam, "house %2d       ", ipl);
+    } else {
+      sprintf(pnam, "%-15s", hs_nam[ipl - nhouses]);
+    }
+  } else if (diff_mode == DIFF_DIFF) {
+    sprintf(pnam, "%.3s-%.3s", spnam, spnam2);
+  } else if (diff_mode == DIFF_GEOHEL) {
+    sprintf(pnam, "%.3s-%.3sHel", spnam, spnam2);
+  } else if (diff_mode == DIFF_MIDP) {
+    sprintf(pnam, "%.3s/%.3s", spnam, spnam2);
+  } else {
+    sprintf(pnam, "%-15.15s", spnam);
+  }
+  if (list_hor && strchr(fmt, 'P') == NULL) {
+    sprintf(slon, "%.8s %s", pnam, "long.");
+  } else {  
+    sprintf(slon, "%-14s", "long.");
+  }
+  for (sp = fmt; *sp != '\0'; sp++) {
+    // if (is_house && ipl <= nhouses && strchr("bBsSrRxXuUQnNfFj+-*/=", *sp) != NULL) continue;
+    if (is_house && strchr("bBrRxXuUQnNfFj+-*/=", *sp) != NULL) continue;
+    if (is_ayana && strchr("bBsSrRxXuUQnNfFj+-*/=", *sp) != NULL) continue;
+    if (sp != fmt)
+      fputs(gap,stdout);
+    if (sp == fmt && list_hor && !is_first && strchr("yYJtT", *sp) == NULL)
+      fputs(gap,stdout);	
+    switch(*sp) {
+    case 'y':
+        if (list_hor && ! is_first) {
+	  break;
+	}
+        if (is_label) { printf("year"); break; }
+	printf("%d", jyear);
+	break;
+    case 'Y':
+        if (list_hor && ! is_first) {
+	  break;
+	}
+        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_GEOHEL) {
+	  printf("%d-%dhel", 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 (is_ayana) {
+	  // printf("Ayanamsha       ");
+	  printf("Ayanamsha %s ", swe_get_ayanamsa_name(sid_mode));
+	} else if (diff_mode == DIFF_DIFF || diff_mode == DIFF_GEOHEL) {
+	  printf("%.3s-%.3s", spnam, spnam2);
+	} else if (diff_mode == DIFF_MIDP) {
+	  printf("%.3s/%.3s", spnam, spnam2);
+	} else {
+	  printf("%-15s", spnam);
+	}
+	break;
+    case 'J':
+        if (list_hor && ! is_first) {
+	  break;
+	}
+        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 (list_hor && ! is_first) {
+	  break;
+	}
+        if (is_label) { printf("date    "); break; }
+	printf("%02d.%02d.%04d", jday, jmon, jyear);
+	if (gregflag == SE_JUL_CAL) printf("j");
+	if (jut != 0 || step_in_minutes || step_in_seconds ) {
+	  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(" TT");
+	}
+	break;
+    case 't':
+        if (list_hor && ! is_first) {
+	  break;
+	}
+        if (is_label) { printf("date"); break; }
+	printf("%02d%02d%02d", jyear % 100, jmon, jday);
+	break;
+    case 'L':
+        if (is_label) { printf(slon); break; }
+	if (psp != NULL && (*psp == 'q' || *psp == 'y')) { /* delta t or time equation */
+	  printf("%# 11.7f", x[0]);
+	  printf("s");
+	  break;
+	}
+	fputs(dms(x[0], round_flag),stdout);
+	break;
+    case 'l':
+        if (is_label) { printf(slon); break; }
+	if (output_extra_prec)
+	  printf("%# 11.11f", x[0]);
+	else
+	  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(slon); 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; }
+		if (output_extra_prec)
+		  printf("%# 11.9f", x[3]);
+		else
+		  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; }
+		if (output_extra_prec)
+		  printf("%# 11.9f", x[4]);
+		else
+		  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; }
+		if (output_extra_prec)
+		  printf("%# 11.9f", xequ[3]);
+		else
+		  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; }
+		if (output_extra_prec)
+		  printf("%# 11.9f", xequ[4]);
+		else
+		  printf("%# 11.7f", xequ[4]);
+		break;
+	      case 'R':   /* speed! */
+	      case 'r':   /* speed! */
+		if (is_label) { printf("AU/day"); break; }
+		if (output_extra_prec)
+		  printf("%# 16.14f", x[5]);
+		else
+		  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') {
+	  int flag = round_flag;
+	  if (is_house) flag |= BIT_ALLOW_361;	// speed of houses can be > 360
+	  if (is_label) { printf("deg/day"); break; }
+	  fputs(dms(x[3], flag),stdout);
+	} else {
+	  if (is_label) { printf("deg/day"); break; }
+	  if (output_extra_prec)
+	    printf("%# 11.17f", x[3]);
+	  else
+	    printf("%# 11.7f", x[3]);
+	}
+	break;
+    case 'B':
+	if (is_label) { printf("lat.    "); break; }
+	if (*psp == 'q') { /* delta t */
+	  printf("%# 11.7f", x[1]);
+	  printf("h");
+	  break;
+	}
+	fputs(dms(x[1], round_flag),stdout);
+	break;
+    case 'b':
+	if (is_label) { printf("lat.    "); break; }
+	if (output_extra_prec)
+	  printf("%# 11.11f", x[1]);
+	else
+	  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; }
+	if (output_extra_prec)
+	  printf("%# 11.11f", xequ[0]);
+	else
+	  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; }
+	if (output_extra_prec)
+	  printf("%# 11.11f", xequ[1]);
+	else
+	  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; }
+	if (output_extra_prec)
+	  printf("%# 16.14f", x[2]);
+	else
+	  printf("%# 14.9f", x[2]);
+	break;
+    case 'W':
+	if (is_label) { printf("distLY   "); break; }
+	printf("%# 14.9f", x[2] * SE_AUNIT_TO_LIGHTYEAR);
+	break;
+    case 'w':
+	if (is_label) { printf("distkm   "); break; }
+	printf("%# 14.9f", x[2] * SE_AUNIT_TO_KM);
+	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 'q':
+	if (is_label) { printf("reldist"); break; }
+        dar = get_geocentric_relative_distance(te, ipl, iflag, serr);
+	printf("% 5d", dar);
+	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;
+	if (output_extra_prec) {
+	  printf("%# .17f", xcartq[0]/ar);
+	  fputs(gap,stdout);
+	  printf("%# .17f", xcartq[1]/ar);
+	  fputs(gap,stdout);
+	  printf("%# .17f", xcartq[2]/ar);
+	} else {
+	  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; 
+	  }
+	  if (use_dms)
+	    fputs(dms(xasc[0], round_flag|BIT_ZODIAC),stdout);
+	  else
+	    printf("%# 11.7f", xasc[0]);
+	  fputs(gap,stdout);
+	  if (use_dms)
+	    fputs(dms(xdsc[0], round_flag|BIT_ZODIAC),stdout);
+	  else
+	    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;
+//	fprintf(stderr, "c=%c\n", *sp);
+	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);
+	    fputs(gap, stdout);
+	    fputs("focus", 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') ) {
+	  fputs(gap,stdout);
+	  printf("%# 11.7f", xfoc[0]);
+	}
+      };
+      break;
+    case '+':
+	if (is_house) break;
+        if (is_label) { printf("phase"); break; }
+	if (strchr(fmt, 'l') != NULL) {	// if decimal longitude is present, do phae angle also decimal
+	  printf("%# 11.7f", attr[0]);
+	} else {
+	  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;
+	if (strchr(fmt, 'l') != NULL) {	// if decimal longitude is present, do elongation also decimal
+	  printf("%# 11.7f", attr[2]);
+	} else {
+	  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.3fm", 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;
+      }
+    case 'm': {	// Meridian distance
+	if (is_label) { printf("MD      "); break; }
+	double md = swe_difdeg2n(xequ[0], armc);
+	if (md < 0) md = -md;
+	if (output_extra_prec)
+	  printf("%# 11.11f", md);
+	else
+	  printf("%# 11.7f", md);
+    	break;
+      }
+    case 'z': {	// Zenith distance
+	if (is_label) { printf("ZD      "); break; }
+	swe_azalt(tut, SE_EQU2HOR, geopos, datm[0], datm[1], xequ, xaz);
+	double zd = 90 - xaz[1];
+	if (output_extra_prec)
+	  printf("%# 11.11f", zd);
+	else
+	  printf("%# 11.7f", zd);
+    	break;
+      }
+    }     /* switch */
+  }       /* for sp */
+  if (! list_hor)
+    printf("\n");
+  return OK;
+}
+
+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
+  if (xv >= 360 && !(iflg & BIT_ALLOW_361))
+    xv = 0;
+  *s = '\0';
+  if (iflg & SEFLG_EQUATORIAL)
+    c = "h";
+  if (xv < 0) {
+    xv = -xv;
+    sgn = -1;
+  } else {
+    sgn = 1;
+  }
+  if (iflg & BIT_ROUND_MIN) {
+    if (!(iflg & BIT_ALLOW_361))
+      xv = swe_degnorm(xv + 0.5/60);
+  } else if (iflg & BIT_ROUND_SEC) {
+    if (!(iflg & BIT_ALLOW_361))
+      xv = swe_degnorm(xv + 0.5/3600);
+  } else {
+  /* rounding 0.9999999999 to 1 */
+    if (output_extra_prec)
+      xv += (xv < 0 ? -1 : 1 ) * 0.000000005 / 3600.0;
+    else
+      xv += (xv < 0 ? -1 : 1 ) * 0.00005 / 3600.0;
+  }
+  if (iflg & BIT_ZODIAC) {
+    izod = (int) (xv / 30); 
+    if (izod == 12) izod = 0;
+    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_prec) {
+    k = (int32) (xv * 100000000);
+    sprintf(s1, ".%08d", k);
+  } else {
+    k = (int32) (xv * 10000);
+    sprintf(s1, ".%04d", k);
+  }
+  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 'x': /* swetest: sidereal time */
+  case 'b': /* swetest: ayanamsha */
+  case 's': /* swetest: an asteroid, with number given in -xs[number] */
+  case 'v': /* swetest: a planetary moon, with number given in -xv[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, serr);
+    }
+  }
+  *t_ret = t_ut;
+  return iflgret;
+}
+
+static int32 orbital_elements(double tjd_et, int32 ipl, int32 iflag, char *serr) 
+{
+  int32 retval;
+  double dret[20], jut;
+  int32 jyear, jmon, jday;
+  char sdateperi[20];
+  retval = swe_get_orbital_elements(tjd_et, ipl, iflag, dret, serr);
+  if (retval == ERR) {
+    printf("%s\n", serr);
+    return ERR;
+  } else {
+    swe_revjul(dret[14], gregflag, &jyear, &jmon, &jday, &jut);
+    sprintf(sdateperi, "%2d.%02d.%04d,%s", jday, jmon, jyear, hms(jut,BIT_LZEROES));
+    printf("semiaxis         \t%f\neccentricity     \t%f\ninclination      \t%f\nasc. node       \t%f\narg. pericenter  \t%f\npericenter       \t%f\n", dret[0], dret[1], dret[2], dret[3], dret[4], dret[5]);
+    printf("mean longitude   \t%f\nmean anomaly     \t%f\necc. anomaly     \t%f\ntrue anomaly     \t%f\n", dret[9], dret[6], dret[8], dret[7]);
+    printf("time pericenter  \t%f %s\ndist. pericenter \t%f\ndist. apocenter  \t%f\n", dret[14], sdateperi, dret[15], dret[16]);
+    printf("mean daily motion\t%f\nsid. period (y)  \t%f\ntrop. period (y) \t%f\nsynodic cycle (d)\t%f\n", dret[11], dret[10], dret[12], dret[13]);
+  }
+  return OK;
+}
+
+static void insert_gap_string_for_tabs(char *sout, char *gap)
+{
+  char *sp;
+  char s[LEN_SOUT];
+  if (!have_gap_parameter)
+    return;
+  if (*gap == '\t')
+    return;
+  while((sp = strchr(sout, '\t')) != NULL && strlen(sout) + strlen(gap) < LEN_SOUT) {
+    strcpy(s, sp + 1);
+    strcpy(sp, gap);
+    strcat(sp, s);
+  }
+}
+
+static int32 print_rise_set_line(double trise, double tset, double *geopos, char *serr) {
+  double t0;
+  int retc = OK;
+  *sout = '\0';
+  if (trise != 0) retc = ut_to_lmt_lat(trise, geopos, &(trise), serr);
+  if (tset != 0) retc = ut_to_lmt_lat(tset, geopos, &(tset), serr);
+  strcpy(sout, "rise     ");
+  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+  if (trise == 0) {
+    strcat(sout, "         -\t           -    ");
+  } else {
+    swe_revjul(trise, gregflag, &jyear, &jmon, &jday, &jut);
+    sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s    ", jday, jmon, jyear, hms(jut,BIT_LZEROES));
+  }
+  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+  strcat(sout, "set      ");
+  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+  if (tset == 0 ) {
+  strcat(sout, "         -\t           -    ");
+  } else {
+    swe_revjul(tset, gregflag, &jyear, &jmon, &jday, &jut);
+    sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s    ", jday, jmon, jyear, hms(jut,BIT_LZEROES));
+  }
+  if (trise != 0 && tset != 0) {
+    if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+    sprintf(sout + strlen(sout), "dt =");
+    if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+    t0 = (tset - trise) * 24;
+    sprintf(sout + strlen(sout), "%s", hms(t0, BIT_LZEROES));
+  }
+  strcat(sout, "\n");
+  if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+  do_printf(sout);
+  return retc;
+}
+
+static int32 call_rise_set(double t_ut, int32 ipl, char *star, int32 whicheph, double *geopos, char *serr)
+{
+  int ii, rval, loop_count;
+  int32 rsmi = 0;
+  double dayfrac = 0.0001;
+  double tret[10], trise, tset, tnext, tret1sv = 0;
+  AS_BOOL do_rise, do_set;
+  AS_BOOL last_was_empty = FALSE;
+  int32 retc = OK;
+  int rsmior = 0;
+  if (norefrac) rsmior |= SE_BIT_NO_REFRACTION;
+  if (disccenter) rsmior |= SE_BIT_DISC_CENTER;
+  if (discbottom) rsmior |= SE_BIT_DISC_BOTTOM;
+  if (hindu) rsmior |= SE_BIT_HINDU_RISING;
+  if (fabs(geopos[1]) < 60 && ipl >= SE_SUN && ipl <= SE_PLUTO) 
+    dayfrac = 0.01;
+  swe_set_topo(geopos[0], geopos[1], geopos[2]); 
+  // "geo. long 8.000000, lat 47.000000, alt 0.000000"
+  if (with_header)
+    printf("\ngeo. long %f, lat %f, alt %f", geopos[0], geopos[1], geopos[2]);
+  do_printf("\n");
+  tnext = t_ut;
+  // the code is designed for looping with -nxxx over many days, during which
+  // the object might become circumpolar, or never rise at all.
+  while (special_event == SP_RISE_SET && tnext < t_ut + nstep) {
+    // the following 'if' avoids unnecessary calculations for circumpolar 
+    // objects. even without it, the output would be correct, but 
+    // could be considerably slower.
+    if (last_was_empty && (star == NULL || *star == '\0')) {
+      rval = swe_calc_ut(tnext, ipl, whicheph|SEFLG_EQUATORIAL, tret, serr);
+      if (rval >= 0 ) {
+        double edist = geopos[1] + tret[1];
+        double edist2 = geopos[1] - tret[1];
+        if ((edist - 2 > 90 || edist + 2 < -90) 
+        || (edist2 - 2 > 90 || edist2 + 2 < -90)) {
+	  tnext += 1;
+	  continue;
+	}
+      }
+    }
+    /* rising */
+    rsmi = SE_CALC_RISE | rsmior;
+    rval= swe_rise_trans(tnext, ipl, star, whicheph, rsmi, geopos, datm[0], datm[1], &trise, serr);
+    if (rval == ERR) {
+      do_printf(serr);
+      exit(0);
+    } 
+    do_rise = (rval == OK);
+    /* setting */
+    rsmi = SE_CALC_SET | rsmior;
+    do_set = FALSE;
+    loop_count = 0;
+    //tnext = trise; // dieter 14-feb-17
+    while (! do_set && loop_count < 2) {
+      rval = swe_rise_trans(tnext, ipl, star, whicheph, rsmi, geopos, datm[0], datm[1], &tset, serr);
+      if (rval == ERR) {
+	do_printf(serr);
+	exit(0);
+      } 
+      do_set = (rval == OK);
+      if (!do_set && do_rise ) {
+        tnext = trise;
+      }
+      loop_count++;
+    }
+    if (do_rise && do_set && trise > tset) {
+      do_rise = FALSE;	// ignore rises happening before setting
+      trise = 0;	// we hope that exact time 0 never happens, is highly unlikely.
+    }
+    if (do_rise && do_set) {
+      rval = print_rise_set_line(trise, tset, geopos, serr);
+      last_was_empty = FALSE;
+      tnext = tset + dayfrac;
+    } else if (do_rise && !do_set) {
+      rval = print_rise_set_line(trise, 0, geopos, serr);
+      last_was_empty = FALSE;
+      tnext = trise + dayfrac;
+    } else if (do_set && ! do_rise) {
+      tnext = tset + dayfrac;
+      rval = print_rise_set_line(0, tset, geopos, serr);
+      last_was_empty = FALSE;
+    } else { // neither rise nor set 
+      // for sequences of days without rise or set, the line '-   -' is printed only once.
+      if (! last_was_empty) rval = print_rise_set_line(0, 0, geopos, serr);
+      tnext += 1;
+      last_was_empty = TRUE;
+    }
+    if (rval == ERR) {
+      do_printf(serr);
+      exit(0);
+    }
+    if (nstep == 1) break;
+  }
+  /* swetest -metr
+   * calculate and print transits over meridian (midheaven and lower
+   * midheaven */
+  if (special_event == SP_MERIDIAN_TRANSIT) {
+  /* loop over days */
+    for (ii = 0; ii < nstep; ii++, t_ut = tret1sv + 0.001) {
+      /* transit over midheaven */
+      if (swe_rise_trans(t_ut, ipl, star, whicheph, SE_CALC_MTRANSIT, geopos, datm[0], datm[1], &(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, datm[0], datm[1], &(tret[1]), serr) != OK) {
+	do_printf(serr);
+	return ERR;
+      } 
+      tret1sv = tret[1];
+      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 (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[0] == 0 || tret[0] > tret[1]) strcat(sout, "         -\t           -    ");
+      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));
+      }
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      strcat(sout, "itransit ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[1] == 0) strcat(sout, "         -\t           -    \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));
+      }
+      if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+      do_printf(sout);
+    }
+  }
+  return retc;
+}
+
+static char* get_gregjul(int gregflag, int year)
+{
+  if (gregflag == SE_JUL_CAL) return "jul";
+  if (year < 1700) return "greg";
+  return "";
+}
+
+// print lon and lat string in minute precision
+static void format_lon_lat(char *slon, char *slat, double lon, double lat)
+{
+  int roundflag, ideg, imin, isec, isgn;
+  double dsecfr;
+  char c;
+  roundflag = SE_SPLIT_DEG_ROUND_SEC;
+  swe_split_deg(lon, roundflag, &ideg, &imin, &isec, &dsecfr, &isgn);
+  c = (lon < 0) ? 'w' : 'e';
+  sprintf(slon, "%d%c%02d%02d", abs(ideg), c, imin, isec);
+  swe_split_deg(lat, roundflag, &ideg, &imin, &isec, &dsecfr, &isgn);
+  c = (lat < 0) ? 's' : 'n';
+  sprintf(slat, "%d%c%02d%02d", abs(ideg), c, imin, isec);
+}
+
+static int32 call_lunar_eclipse(double t_ut, int32 whicheph, int32 special_mode, double *geopos, char *serr)
+{
+  int i, ii, retc = OK, eclflag, ecl_type = 0;
+  int rval, ihou, imin, isec, isgn;
+  double dfrc, attr[30], dt, xx[6], geopos_max[3];
+  char s1[AS_MAXCH], s2[AS_MAXCH], sout_short[2 * LEN_SOUT], sfmt[AS_MAXCH], *styp = "none", *sgj;
+  char slon[8], slat[8], saros[40];
+  geopos_max[0] = 0; geopos_max[1] = 0;
+  /* no selective eclipse type set, set all */
+  if (with_chart_link) do_printf("<pre>");
+  if ((search_flag & SE_ECL_ALLTYPES_LUNAR) == 0)
+    search_flag |= SE_ECL_ALLTYPES_LUNAR;
+  // "geo. long 8.000000, lat 47.000000, alt 0.000000"
+  if (special_mode & SP_MODE_LOCAL) {
+    if (with_header)
+      printf("\ngeo. long %f, lat %f, alt %f", geopos[0], geopos[1], geopos[2]);
+  }
+  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);
+	    if (retc == 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");
+      swe_revjul(t_ut, gregflag, &jyear, &jmon, &jday, &jut);
+      sgj = get_gregjul(gregflag, jyear);
+      /*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(saros, "%d/%d", (int) attr[9], (int) attr[10]);
+      sprintf(sout_short, "%s\t%2d.%2d.%4d%s\t%s\t%.3f\t%s\t%s\n", sout, jday, jmon, jyear, sgj, hms(jut,0), attr[8],s1, saros);
+      sprintf(sout + strlen(sout), "%2d.%02d.%04d%s\t%s\t%.4f/%.4f\tsaros %s\t%.6f\n", jday, jmon, jyear, sgj, hms(jut,BIT_LZEROES), attr[0],attr[1], saros, t_ut);
+      /* second line:
+       * eclipse times, penumbral, partial, total begin and end */
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (eclflag & SE_ECL_PENUMBBEG_VISIBLE)
+	sprintf(sout + strlen(sout), "  %s ", hms_from_tjd(tret[6])); 
+      else
+	strcat(sout, "      -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (eclflag & SE_ECL_PARTBEG_VISIBLE)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2])); 
+      else
+	strcat(sout, "    -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (eclflag & SE_ECL_TOTBEG_VISIBLE)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4])); 
+      else
+	strcat(sout, "    -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (eclflag & SE_ECL_TOTEND_VISIBLE)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5])); 
+      else
+	strcat(sout, "    -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (eclflag & SE_ECL_PARTEND_VISIBLE)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3])); 
+      else
+	strcat(sout, "    -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (eclflag & SE_ECL_PENUMBEND_VISIBLE)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[7])); 
+      else
+	strcat(sout, "    -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "dt=%.1f", swe_deltat_ex(tret[0], whicheph, serr) * 86400.0);
+      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)) {
+	styp = "Total";
+        strcpy(sout, "total ");
+	ecl_type = ECL_LUN_TOTAL;
+      }
+      if ((eclflag & SE_ECL_PENUMBRAL)) {
+        styp = "Penumbral";
+        strcpy(sout, "penumb. ");
+	ecl_type = ECL_LUN_PENUMBRAL;
+      }
+      if ((eclflag & SE_ECL_PARTIAL)) {
+        styp = "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);
+	    if (retc == ERR) {
+	      do_printf(serr);
+	      return ERR;
+	    }
+	  }
+	}
+      }
+      t_ut = tret[0];
+      rval = swe_calc_ut(t_ut, SE_MOON, whicheph|SEFLG_EQUATORIAL, xx, s1);
+      if (rval < 0) 
+        strcat(s1, "\n");
+      do_printf(s1);
+      swe_revjul(t_ut, gregflag, &jyear, &jmon, &jday, &jut);
+      geopos_max[0] = swe_degnorm(xx[0]- swe_sidtime(t_ut) * 15);
+      if (geopos_max[0] > 180) geopos_max[0] -= 360;
+      geopos_max[1] = xx[1];
+      sgj = get_gregjul(gregflag, jyear);
+      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(saros, "%d/%d", (int) attr[9], (int) attr[10]);
+      sprintf(sout_short, "%s\t%2d.%2d.%4d%s\t%s\t%.3f\t%s\t%s\n", sout, jday, jmon, jyear, sgj, hms(jut,0), attr[8],s1, saros);
+      //sprintf(sout + strlen(sout), "%2d.%02d.%04d%s\t%s\t%.4f/%.4f\tsaros %s\t%.6f\tdt=%.2f\n", jday, jmon, jyear, sgj, hms(jut,BIT_LZEROES), attr[0],attr[1], saros, t_ut, swe_deltat_ex(t_ut, whicheph, serr) * 86400);
+      sprintf(sout + strlen(sout), "%2d.%02d.%04d%s\t%s\t%.4f/%.4f\tsaros %s\t%.6f\n", jday, jmon, jyear, sgj, hms(jut,BIT_LZEROES), attr[0],attr[1], saros, t_ut);
+      /* second line:
+       * eclipse times, penumbral, partial, total begin and end */
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "  %s ", hms_from_tjd(tret[6])); 
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[2] != 0)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2])); 
+      else
+	strcat(sout, "   -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[4] != 0)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4])); 
+      else
+	strcat(sout, "   -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[5] != 0)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5])); 
+      else
+	strcat(sout, "   -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[3] != 0)
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3])); 
+      else
+	strcat(sout, "   -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "%s", hms_from_tjd(tret[7])); 
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "dt=%.1f", swe_deltat_ex(tret[0], whicheph, serr) * 86400.0);
+      strcat(sout, "\n");
+      if (special_mode & SP_MODE_HOCAL) {
+	swe_split_deg(jut, SE_SPLIT_DEG_ROUND_MIN, &ihou, &imin, &isec, &dfrc, &isgn);
+	sprintf(sout, "\"%04d%s %02d %02d %02d.%02d %d\",\n", jyear, sgj, jmon, jday, ihou, imin, ecl_type);
+      } 
+      sprintf(sout + strlen(sout), "\t%s\t%s\n", strcpy(s1, dms(geopos_max[0], BIT_ROUND_SEC)), strcpy(s2, dms(geopos_max[1], BIT_ROUND_SEC)));
+    }
+    //dt = (tret[7] - tret[6]) * 24 * 60;
+    //sprintf(sout + strlen(sout), "\t%d min %4.2f sec\n", (int) dt, fmod(dt, 1) * 60);
+    if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+    if (short_output) {
+      do_printf(sout_short);
+    } else {
+      do_printf(sout);
+    }
+    if (with_chart_link) {
+      char snat[2 * AS_MAXCH];
+      char stim[80];
+      int iflg = 0; 
+      char cal = gregflag ? 'g' : 'j';
+      strcpy(stim, hms(jut,BIT_LZEROES));
+      format_lon_lat(slon, slat, geopos_max[0], geopos_max[1]);
+      while (*stim == ' ') our_strcpy(stim, stim + 1);
+      if (*stim == '0') our_strcpy(stim, stim + 1);
+      sprintf(snat, "Lunar Eclipse %s,%s,e,%d,%d,%d,%s,h0e,%cnu,%d,Moon Zenith location,,%s,%s,u,0,0,0", saros, styp, jday, jmon, jyear, stim, cal, iflg, slon, slat);
+      sprintf(sout, "<a href='https://www.astro.com/cgi/chart.cgi?muasp=1;nhor=1;act=chmnat;nd1=%s;rs=1;iseclipse=1' target='eclipse'>chart link</a>\n\n", snat);
+      do_printf(sout);
+    }
+  }
+  if (with_chart_link) do_printf("</pre>\n");
+  return OK;
+}
+
+static int32 call_solar_eclipse(double t_ut, int32 whicheph, int32 special_mode, double *geopos, char *serr)
+{
+  int i, ii, retc = OK, eclflag, ecl_type = 0;
+  double dt, tret[30], attr[30], geopos_max[3];
+  char slon[8], slat[8], saros[20];
+  char s1[AS_MAXCH], s2[AS_MAXCH], sout_short[AS_MAXCH + LEN_SOUT], *styp = "none", *sgj;
+  AS_BOOL has_found = FALSE;
+  /* no selective eclipse type set, set all */
+  if (with_chart_link) do_printf("<pre>");
+  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]); 
+    // "geo. long 8.000000, lat 47.000000, alt 0.000000"
+    if (with_header)
+      printf("\ngeo. long %f, lat %f, alt %f", 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 (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	if (!has_found) {
+	  ii--;
+	} else {
+	  swe_calc(t_ut + swe_deltat_ex(t_ut, whicheph, serr), 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);
+		if (retc == ERR) {
+		  do_printf(serr);
+		  return ERR;
+		}
+	      }
+	    }
+	  }
+	  t_ut = tret[0];
+	  swe_revjul(t_ut, gregflag, &jyear, &jmon, &jday, &jut);
+	  dt = (tret[3] - tret[2]) * 24 * 60;
+	  sgj = get_gregjul(gregflag, jyear);
+	  sprintf(saros, "%d/%d", (int) attr[9], (int) attr[10]);
+	  sprintf(sout + strlen(sout), "%2d.%02d.%04d%s\t%s\t%.4f/%.4f/%.4f\tsaros %s\t%.6f\n", jday, jmon, jyear, sgj,hms(jut,BIT_LZEROES), attr[8], attr[0], attr[2], saros, 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 (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	  if (eclflag & SE_ECL_2ND_VISIBLE) {
+	    sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2])); 
+	  } else {
+	    strcat(sout, "   -         ");
+	  }
+	  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	  if (eclflag & SE_ECL_3RD_VISIBLE) {
+	    sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3])); 
+	  } else {
+	    strcat(sout, "   -         ");
+	  }
+	  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	  if (eclflag & SE_ECL_4TH_VISIBLE) {
+	    sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4])); 
+	  } else {
+	    strcat(sout, "   -         ");
+	  }
+	  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+#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
+	  sprintf(sout + strlen(sout), "dt=%.1f", swe_deltat_ex(tret[0], whicheph, serr) * 86400.0);
+	  strcat(sout, "\n");
+	  if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+	  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)) {
+	styp = "Total";
+        strcpy(sout, "total");
+	ecl_type = ECL_SOL_TOTAL;
+      }
+      if ((eclflag & SE_ECL_ANNULAR)) {
+	styp = "Annular";
+        strcpy(sout, "annular");
+	ecl_type = ECL_SOL_ANNULAR;
+      }
+      if ((eclflag & SE_ECL_ANNULAR_TOTAL)) {
+	styp = "Annular-Total";
+        strcpy(sout, "ann-tot");
+	ecl_type = ECL_SOL_ANNULAR;	/* by Alois: what is this ? */
+      }
+      if ((eclflag & SE_ECL_PARTIAL)) {
+	styp = "Partial";
+        strcpy(sout, "partial");
+	ecl_type = ECL_SOL_PARTIAL;
+      }
+      if ((eclflag & SE_ECL_NONCENTRAL) && !(eclflag & SE_ECL_PARTIAL))
+        strcat(sout, " non-central");
+      sprintf(sout + strlen(sout), " solar\t");
+      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);
+	    if (retc == ERR) {
+	      do_printf(serr);
+	      return ERR;
+	    }
+	  }
+	}
+      }
+      swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
+      sgj = get_gregjul(gregflag, jyear);
+      sprintf(saros, "%d/%d", (int) attr[9], (int) attr[10]);
+      sprintf(sout_short, "%s\t%2d.%2d.%4d%s\t%s\t%.3f", sout, jday, jmon, jyear, sgj, hms(jut,0), attr[8]);
+      sprintf(sout + strlen(sout), "%2d.%02d.%04d%s\t%s\t%f km\t%.4f/%.4f/%.4f\tsaros %s\t%.6f\n", jday, jmon, jyear, sgj, hms(jut,0), attr[3], attr[8], attr[0], attr[2], saros, tret[0]);
+      sprintf(sout + strlen(sout), "\t%s ", hms_from_tjd(tret[2])); 
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[4] != 0) {
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4])); 
+      } else {
+	strcat(sout, "   -         ");
+      }
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[5] != 0) {
+	sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5])); 
+      } else {
+	strcat(sout, "   -         ");
+      }
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "%s", hms_from_tjd(tret[3])); 
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "dt=%.1f", swe_deltat_ex(tret[0], whicheph, serr) * 86400.0);
+      strcat(sout, "\n");
+      sprintf(sout + strlen(sout), "\t%s\t%s", strcpy(s1, dms(geopos_max[0], BIT_ROUND_SEC)), strcpy(s2, dms(geopos_max[1], BIT_ROUND_SEC)));
+      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%s %02d %02d %02d.%02d %d\",\n", jyear, sgj, jmon, jday, ihou, imin, ecl_type);
+      } 
+      /*printf("len=%ld\n", strlen(sout));*/
+      if (short_output) {
+	do_printf(sout_short);
+      } else {
+	if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+	do_printf(sout);
+      }
+      if (with_chart_link) {
+	char snat[2 * AS_MAXCH];
+	char stim[80];
+	int iflg = 0; // NAT_IFLG_UNKNOWN_TIME;
+	char cal = gregflag ? 'g' : 'j';
+	format_lon_lat(slon, slat, geopos_max[0], geopos_max[1]);
+	strcpy(stim, hms(jut,BIT_LZEROES));
+	while (*stim == ' ') our_strcpy(stim, stim + 1);
+	if (*stim == '0') our_strcpy(stim, stim + 1);
+	sprintf(snat, "Solar Eclipse %s,%s,e,%d,%d,%d,%s,h0e,%cnu,%d,Location of Maximum,,%s,%s,u,0,0,0", saros, styp, jday, jmon, jyear, stim, cal, iflg, slon, slat);
+	sprintf(sout, "<a href='https://www.astro.com/cgi/chart.cgi?muasp=1;nhor=1;act=chmnat;nd1=%s;rs=1;iseclipse=1;topo=1' target='eclipse'>chart link</a>\n\n", snat);
+	do_printf(sout);
+      }
+    }
+  }
+  if (with_chart_link) do_printf("</pre>\n");
+  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 = OK;
+  double dt, tret[30], attr[30], geopos_max[3];
+  char s1[AS_MAXCH], s2[AS_MAXCH];
+  AS_BOOL has_found = FALSE;
+  int nloops = 0;
+  /* 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]); 
+    if (with_header)
+      printf("\ngeo. long %f, lat %f, alt %f", geopos[0], geopos[1], geopos[2]);
+  }
+  do_printf("\n");
+  for (ii = 0; ii < nstep; ii++) {
+    *sout = '\0';
+    nloops++;
+    if (nloops > SEARCH_RANGE_LUNAR_CYCLES) {
+      sprintf(serr, "event search ended after %d lunar cycles at jd=%f\n", SEARCH_RANGE_LUNAR_CYCLES, t_ut);
+      do_printf(serr);
+      return ERR;
+    }
+    if (special_mode & SP_MODE_LOCAL) {
+      /* * local search for occultation, test one lunar cycle only (SE_ECL_ONE_TRY) */
+      if (ipl != SE_SUN) {
+        search_flag &= ~(SE_ECL_ANNULAR|SE_ECL_ANNULAR_TOTAL);
+	if (search_flag == 0)
+	  search_flag = SE_ECL_ALLTYPES_SOLAR;
+      }
+      if ((eclflag = swe_lun_occult_when_loc(t_ut, ipl, star, whicheph, geopos, tret, attr, direction_flag|SE_ECL_ONE_TRY, serr)) == ERR) {
+        do_printf(serr);
+        return ERR;
+      } else if (eclflag == 0) {  /* event not found, try next conjunction */
+        t_ut = tret[0] + direction * 10;  /* try again with start date increased by 10 */
+	ii--;
+      } 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);
+	      if (retc == ERR) {
+		do_printf(serr);
+		return ERR;
+	      }
+	    }
+	  }
+	}
+	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 */
+	}
+	if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	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%f\t%.6f\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[0], tret[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 (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	  if (eclflag & SE_ECL_2ND_VISIBLE)
+	    sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[2])); 
+	  else
+	    strcat(sout, "   -         ");
+	  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	  if (eclflag & SE_ECL_3RD_VISIBLE)
+	    sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[3])); 
+	  else
+	    strcat(sout, "   -         ");
+	  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+	  if (eclflag & SE_ECL_4TH_VISIBLE)
+	    sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4])); 
+	  else
+	    strcat(sout, "   -         ");
+	  if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+#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
+	  sprintf(sout + strlen(sout), "dt=%.1f", swe_deltat_ex(tret[0], whicheph, serr) * 86400.0);
+	  strcat(sout, "\n");
+	  if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+	  do_printf(sout);
+	}
+      }
+    }   /* endif search_local */
+    if (!(special_mode & SP_MODE_LOCAL)) {
+    /* * global search for occultations, test one lunar cycle only (SE_ECL_ONE_TRY) */
+      if ((eclflag = swe_lun_occult_when_glob(t_ut, ipl, star, whicheph, search_flag, tret, direction_flag|SE_ECL_ONE_TRY, serr)) == ERR) {
+        do_printf(serr);
+        return ERR;
+      } 
+      if (eclflag == 0) { /* no occltation was found at next conjunction, try next conjunction */
+	t_ut = tret[0] + direction;
+	ii--;
+	continue;
+      }
+      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);
+      /* for (i = 0; i < 8; i++) {
+        printf("attr[%d]=%.17f\n", i, attr[i]);
+      } */
+      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);
+	    if (retc == ERR) {
+	      do_printf(serr);
+	      return ERR;
+	    }
+	  }
+	}
+      }
+      swe_revjul(tret[0], gregflag, &jyear, &jmon, &jday, &jut);
+      sprintf(sout + strlen(sout), "%2d.%02d.%04d\t%s\t%f km\t%f\t%.6f\n", jday, jmon, jyear, hms(jut,BIT_LZEROES), attr[3], attr[0], tret[0]);
+      sprintf(sout + strlen(sout), "\t%s ", hms_from_tjd(tret[2])); 
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[4] != 0)
+	  sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[4])); 
+      else
+          strcat(sout, "   -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      if (tret[5] != 0)
+	  sprintf(sout + strlen(sout), "%s ", hms_from_tjd(tret[5])); 
+      else
+          strcat(sout, "   -         ");
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "%s", hms_from_tjd(tret[3])); 
+      if (have_gap_parameter) sprintf(sout + strlen(sout), "\t");
+      sprintf(sout + strlen(sout), "dt=%.1f", swe_deltat_ex(tret[0], whicheph, serr) * 86400.0);
+      strcat(sout, "\n");
+      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", (int) dt, fmod(dt, 1) * 60);
+      } 
+      strcat(sout, "\n");
+      if (have_gap_parameter) insert_gap_string_for_tabs(sout, gap);
+      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, double *geopos, double *datm, double *dobs, char *serr)
+{
+  int ii, retc = OK, 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);
+  }
+  if (with_header) {
+    printf("\ngeo. long %f, lat %f, alt %f", geopos[0], geopos[1], geopos[2]);
+    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);
+      if (retc != ERR) retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
+      if (retc != ERR) retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
+      if (retc == ERR) {
+	do_printf(serr);
+	return ERR;
+      }
+    }
+    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);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
+	  if (retc == ERR) {
+	    do_printf(serr);
+	    return ERR;
+	  }
+	}
+	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);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
+	  if (retc == ERR) {
+	    do_printf(serr);
+	    return ERR;
+	  }
+	}
+	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);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
+	  if (retc == ERR) {
+	    do_printf(serr);
+	    return ERR;
+	  }
+	}
+	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);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[1], geopos, &(dret[1]), serr);
+	  if (retc != ERR) retc = ut_to_lmt_lat(dret[2], geopos, &(dret[2]), serr);
+	  if (retc == ERR) {
+	    do_printf(serr);
+	    return ERR;
+	  }
+	}
+	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 = 0;
+  /* risings, settings, meridian transits */
+  if (special_event == SP_RISE_SET || 
+      special_event == SP_MERIDIAN_TRANSIT)
+    retc = call_rise_set(tjd, ipl, star, whicheph, 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, geopos, datm, dobs, serr);
+  return retc;
+}
+
+static char *hms_from_tjd(double tjd)
+{
+  static char s[AS_MAXCH];
+  double x;
+  /* tjd may be negative, 0h corresponds to day number 9999999.5 */ 
+  x = fmod(tjd, 1);  /* may be negative ! */
+  x = fmod(x + 1.5, 1); /* is positive day fraction */
+  sprintf(s, "%s ", hms(x * 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(char *argv0, char *path)
+{
+  char *sp;
+  char *dirglue = DIR_GLUE;
+  size_t pathlen = 0; 
+  /* 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);
+}
+
+static void jd_to_time_string(double jut, char *stimeout)
+{
+  double t2;
+  t2 = jut + 0.5 / 3600000.0; // rounding to millisec
+  sprintf(stimeout, "  % 2d:", (int) t2); // hour
+  t2 = (t2 - (int32) t2) * 60;
+  sprintf(stimeout + strlen(stimeout), "%02d:", (int) t2);  // min
+  t2 = (t2 - (int32) t2) * 60;
+  sprintf(stimeout + strlen(stimeout), "%02d", (int) t2); // sec
+  t2 = (t2 - (int32) t2) * 1000;
+  if ((int32) t2 > 0) {
+    sprintf(stimeout + strlen(stimeout), ".%03d", (int) t2); // millisec, if > 0
+  }
+}
+
+#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
+
+static char *our_strcpy(char *to, char *from)
+{
+  char *sp, s[AS_MAXCH];
+  if (*from == '\0') {
+    *to = '\0';
+    return to;
+  }
+  if (strlen(from) < AS_MAXCH) {
+    strcpy(s, from);
+    strcpy(to, s);
+  } else {
+    sp = strdup(from);
+    if (sp == NULL) {
+      strcpy(to, from);
+    } else {
+      strcpy(to, sp);
+      free(sp);
+    }
+  }
+  return to;
+}