Changeset 1515 in Sophya for trunk/SophyaExt/XAstroPack/xastropack.cc

Timestamp:

Jun 12, 2001, 5:36:47 PM (24 years ago)

Author:

cmv

Message:

intro enum TypAstroCoord cmv 12/6/01

File:

• trunk/SophyaExt/XAstroPack/xastropack.cc (modified) (3 diffs)

trunk/SophyaExt/XAstroPack/xastropack.cc

@@ -7 +7 @@
 //        mois [1,12] (mn)
 //        annee       (yr)
-//        universal time [0,23[ (utc)
-//        Greenwich mean siderial [0,23[ (gst)
-//        Greenwich mean siderial at 0h UT [0,23[ (gst0)
+//        universal time [0,24[ (utc)
+//        Greenwich mean siderial [0,24[ (gst)
+//        Greenwich mean siderial at 0h UT [0,24[ (gst0)
 // EQUATORIALE: ascension droite en heures [0,24[ (ra)
 //              declinaison en degres [-90,90]    (dec)
@@ -15 +15 @@
 // GALACTIQUE: longitude en degres [0,360[  (glng)
 //             latitude en degres [-90,90]  (glat)
-// HORIZONTAL: altitude en degres [-90,90] (alt)
-//             azimuth en degres [0,360[   (az)
-//               (angle round to the east from north+)
-// ECLIPTIQUE: latitude ecliptique en degres [-90,90] (eclat)
-//             lontitude ecliptique en degres [0,360[ (eclng)
-//               (angle round counter clockwise from the vernal equinoxe)
-// GEOGRAPHIE: latitude en degres [-90,90] (north>0) (geolat)
-//             longitude en degres ]-180,180]        (geolng)
-//               (angle + vers l'ouest, - vers l'est)
+// HORIZONTAL: azimuth en degres [0,360[   (az)
+//                 (angle round to the east from north+)
+//             altitude en degres [-90,90] (alt)
+// ECLIPTIQUE: lontitude ecliptique en degres [0,360[ (eclng)
+//                 (angle round counter clockwise from the vernal equinoxe)
+//             latitude ecliptique en degres [-90,90] (eclat)
+// GEOGRAPHIE: longitude en degres ]-180,180]        (geolng)
+//                 (angle + vers l'ouest, - vers l'est)
+//             latitude en degres [-90,90] (north>0) (geolat)
 
 double TrueJDfrMJD(double mjd)
@@ -309 +309 @@
   PlanetPos(mjd,SATURN,ecl,ecb,diamang);
 }
+
+/* Given a coordinate type "typ", convert to standard for astropack */
+int  CoordConvertToStd(TypAstroCoord typ,double& coord1,double& coord2)
+// Return : 0 = OK
+//          1 = Type de coordonnees non connu
+//          2 = Mauvais range pour coord1
+//          4 = Mauvais range pour coord2
+//          6 = Mauvais range pour coord1 et coord2
+{
+  int rc = 0;
+
+  // ---- Equatoriales alpha,delta
+  //    - standard = [0,24[ , [-90,90]
+  if(typ&TypCoordEq) {
+    if(typ&TypCoordDD) {
+      coord1 = coord1 / 180. * 12.;
+    } else if(typ&TypCoordRR) {
+      coord1 = coord1 / PI * 12.;
+      coord2 = coord2 / PI * 180.;
+    }
+    if(coord1==24.) coord1 = 0.;
+    if(coord1<0.   || coord1>=24.) rc+= 2;
+    if(coord2<-90. || coord2>90. ) rc+= 4;
+
+  // ---- Galactiques gLong, gLat
+  // ---- Horizontales azimuth,altitude
+  // ---- Ecliptiques EclLong,EclLat
+  //    - standard = [0,360[ , [-90,90]
+  } else if( typ&TypCoordGal || typ&TypCoordHor || typ&TypCoordEcl) {
+    if(typ&TypCoordHD) {
+      coord1 = coord1 / 12. * 180.;
+    } else if(typ&TypCoordRR) {
+      coord1 = coord1 / PI * 180.;
+      coord2 = coord2 / PI * 180.;
+    }
+    if(coord1==360.) coord1 = 0.;
+    if(coord1<0.   || coord1>=360.) rc+= 2;
+    if(coord2<-90. || coord2>90. )  rc+= 4;
+
+  } else {          // Coordonnees non-connues
+    rc= 1;
+  }
+
+  return rc;
+}