Changeset 1786 in Sophya for trunk/ArchTOIPipe/ProcWSophya

Timestamp:

Nov 26, 2001, 4:13:38 PM (24 years ago)

Author:

aubourg

Message:

debug bornes

File:

• trunk/ArchTOIPipe/ProcWSophya/toi2map.cc (modified) (3 diffs)

trunk/ArchTOIPipe/ProcWSophya/toi2map.cc

@@ -3 +3 @@
 //                               Christophe Magneville
 //                               Reza Ansari
-// $Id: toi2map.cc,v 1.15 2001-11-13 16:22:47 aubourg Exp $
+// $Id: toi2map.cc,v 1.16 2001-11-26 15:13:38 aubourg Exp $
 
 #include "machdefs.h"
@@ -162 +162 @@
   // sphere theta entre [0,Pi]   en radians
   double phi=-1.,theta;
-  CoordConvertToStd(mTypCoorIn,coord1,coord2);
+  int crc = CoordConvertToStd(mTypCoorIn,coord1,coord2);
+
+  // On tente de remettre dans l'intervalle si on est hors intervalle
+  if ((mTypCoorIn&TypCoordEq) && (coord1<0))
+    coord1 += 24;
+  if ((mTypCoorIn&TypCoordEq) && (coord1>24))
+    coord1 -= 24;
+  if ((mTypCoorIn&TypCoordEq) && (coord2>90))
+    coord2 -= 180;
+  if ((mTypCoorIn&TypCoordEq) && (coord2<-90))
+    coord2 += 180;
+
+  if ((mTypCoorIn&TypCoordGal) && (coord1<0))
+    coord1 += 360;
+  if ((mTypCoorIn&TypCoordGal) && (coord1>360))
+    coord1 -= 360;
+  if ((mTypCoorIn&TypCoordGal) && (coord2>90)) 
+    coord2 -= 180;
+  if ((mTypCoorIn&TypCoordGal) && (coord2<-90))
+    coord2 += 180;
+
   if(mTypCoorIn&TypCoordEq && mTypCoorOut&TypCoordGal) { // Eq -> Gal
     EqtoGal(mjd,coord1,coord2,&coord1,&coord2);
@@ -175 +195 @@
   }
   theta = (90.-coord2) * M_PI/180.;
+
+  if(phi<0.   || phi>=2*M_PI || theta<0. || theta>=M_PI) {
+    cout << "out of range " << theta << " " << phi << " "
+         << coord1 << " " << coord2 << endl;
+  }
+     
+
   if(phi<0.   || phi>=2*M_PI) continue;
   if(theta<0. || theta>=M_PI) continue;