Changeset 2551 in Sophya for trunk/SophyaExt/XephemAstroLib/testxephem.c

Timestamp:

Jun 15, 2004, 6:54:12 PM (21 years ago)

Author:

cmv

Message:

nouvelle version de xephem/libastro (3.6) cmv 15/6/04

File:

• trunk/SophyaExt/XephemAstroLib/testxephem.c (modified) (3 diffs)

trunk/SophyaExt/XephemAstroLib/testxephem.c

@@ -4 +4 @@
 #include <math.h>
 #include <unistd.h>
-#include "P_.h"
+/* include "P_.h" */
 #include "astro.h"
 
@@ -11 +11 @@
 int main(int narg,char** arg)
 {
-  int year=1989,month=7,day=21; /* UTC*/
+ int year=1989,month=7,day=21; /* UTC*/
  double hour=12.5;
  int npnum=8, pnum[8]={MERCURY,VENUS,MARS,JUPITER,SATURN,URANUS,NEPTUNE,PLUTO};
  int i;
- double mjd,sunecl,sunecb,sundist,geodist,geoecl,geoecb,diamang,mag,msp,mdp;
+ double mymjd,sunecl,sunecb,sundist,geodist,geoecl,geoecb,diamang,mag,msp,mdp;
 
- cal_mjd(month,day+hour/24.,year,&mjd);
- printf("date %d/%d/%d %gh  -> mjd=%f\n",day,month,year,hour,mjd);
+ cal_mjd(month,day+hour/24.,year,&mymjd);
+ printf("date %d/%d/%d %gh  -> mjd=%f\n",day,month,year,hour,mymjd);
 
  /* Soleil */
- sunpos(mjd,&geoecl,&geodist,&geoecb);
+ sunpos(mymjd,&geoecl,&geodist,&geoecb);
  printf("Sun: geo ecl=%g ecb=%g dist=%g\n",geoecl*R2D,geoecb*R2D,geodist);
 
  /* Lune */
- moon(mjd,&geoecl,&geoecb,&geodist,&msp,&mdp);
+ moon(mymjd,&geoecl,&geoecb,&geodist,&msp,&mdp);
  printf("Moon: geo ecl=%g ecb=%g dist=%g\n",geoecl*R2D,geoecb*R2D,geodist);
 
@@ -31 +31 @@
  printf("--- Planets\n");
  for(i=0;i<npnum;i++) {
-   plans(mjd,pnum[i],&sunecl,&sunecb,&sundist,&geodist,&geoecl,&geoecb,&diamang,&mag);
+   plans(mymjd,pnum[i],&sunecl,&sunecb,&sundist,&geodist,&geoecl,&geoecb,&diamang,&mag);
    printf("pnum=%d: sun ecl=%g ecb=%g dist=%g,",pnum[i],sunecl*R2D,sunecb*R2D,sundist);
    printf(" geo ecl=%g ecb=%g dist=%g,",geoecl*R2D,geoecb*R2D,geodist);