active new WCS code

15 files changed, 70 insertions, 66 deletions

diff --git a/tksao/frame/base.C b/tksao/frame/base.C
index 83911e7..4f8482c 100644
--- a/tksao/frame/base.C
+++ b/tksao/frame/base.C
@@ -314,7 +314,7 @@ void Base::alignWCS(Coord::CoordSystem sys, Coord::SkyFrame sky)
 	       &wcsOrientation, &wcsOrientationMatrix, &wcsRotation);
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void Base::alignWCS(FitsImage* ptr, Coord::CoordSystem sys)
 {
   if (!wcsAlign_ || !ptr || !context->cfits || !hasWCS(wcsSystem_)) {
@@ -390,7 +390,7 @@ void Base::calcAlignWCS(FitsImage* fits1,
   }
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void Base::calcAlignWCS(FitsImage* fits1, FitsImage* fits2, 
 			Coord::CoordSystem sys1, Coord::CoordSystem sys2, 
 			Coord::SkyFrame sky,
@@ -467,7 +467,7 @@ void Base::calcAlignWCS(FitsImage* fits1, FitsImage* fits2,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 Matrix Base::calcAlignWCS(FitsImage* fits1, FitsImage* fits2, 
 			  Coord::CoordSystem sys1, Coord::CoordSystem sys2, 
 			  Coord::SkyFrame sky)
@@ -1611,7 +1611,9 @@ void Base::updateBin(const Matrix& mx)
 
   alignWCS();
   updateColorScale();
-  update(MATRIX); // because we have changed zoom
+  // because we have changed zoom
+  //   now because repeated bin cmds will mult markers/contours relative
+  updateNow(MATRIX);
 
   // update markers call backs
   // wait til matrices have been updated so that any dialogs will print
@@ -1634,7 +1636,9 @@ void Base::updateBlock(const Vector& vv)
 
   alignWCS();
   updateColorScale();
-  update(MATRIX);
+  // because we have changed zoom
+  //   now because repeated bin cmds will mult markers/contours relative
+  updateNow(MATRIX);
 
   // update markers call backs
   // wait til matrices have been updated so that any dialogs will print
diff --git a/tksao/frame/base.h b/tksao/frame/base.h
index 26ed2aa..7899cee 100644
--- a/tksao/frame/base.h
+++ b/tksao/frame/base.h
@@ -528,7 +528,7 @@ public:
   virtual FrameType frameType() =0;
   void calcAlignWCS(FitsImage*, Coord::CoordSystem, Coord::SkyFrame,
 		    Coord::Orientation*, Matrix*, double*);
-#ifndef NEWWCS
+#ifdef OLDWCS
   void calcAlignWCS(FitsImage*, FitsImage*, Coord::CoordSystem, Coord::CoordSystem, Coord::SkyFrame,
 		    Coord::Orientation*, Matrix*, double*, Vector*);
 #endif
diff --git a/tksao/frame/coord.h b/tksao/frame/coord.h
index 6176d23..8f28480 100644
--- a/tksao/frame/coord.h
+++ b/tksao/frame/coord.h
@@ -19,7 +19,7 @@ class Coord {
   enum InternalSystem {WINDOW, CANVAS, WIDGET, USER, REF, 
 		       PANNER, MAGNIFIER, PS};
 
-#ifndef NEWWCS
+#ifdef OLDWCS
   enum CoordSystem {DATA, IMAGE, PHYSICAL, AMPLIFIER, DETECTOR, WCS, 
 		    WCSA, WCSB, WCSC, WCSD, WCSE, WCSF, WCSG, WCSH, WCSI,
 		    WCSJ, WCSK, WCSL, WCSM, WCSN, WCSO, WCSP, WCSQ, WCSR,
diff --git a/tksao/frame/fitsimage.C b/tksao/frame/fitsimage.C
index 27af278..3823e84 100644
--- a/tksao/frame/fitsimage.C
+++ b/tksao/frame/fitsimage.C
@@ -107,7 +107,7 @@ FitsImage::FitsImage(Context* cx, Tcl_Interp* pp)
   refToImage3d = dataToImage3d;
 
   manageWCS_ =1;
-#ifndef NEWWCS
+#ifdef OLDWCS
   wcs_ =NULL;
   ast_ =NULL;
   wcsx_ =NULL;
@@ -175,7 +175,7 @@ FitsImage::~FitsImage()
       delete analysisdata_;
   }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
   if (wcs_) {
     for (int ii=0; ii<MULTWCSA; ii++)
       if (manageWCS_ && wcs_[ii])
@@ -1070,7 +1070,7 @@ void FitsImage::iisSetFileName(const char* fn)
   iisFileName = dupstr(fn);
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::initWCS(FitsHead* hd, FitsHead* prim)
 {
   if (wcs_) {
@@ -1327,7 +1327,7 @@ void FitsImage::initWCS(FitsHead* hd, FitsHead* prim)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::initWCS0(const Vector& pix)
 {
   FitsHead* hd =NULL;
@@ -1388,7 +1388,7 @@ void FitsImage::initWCS0(const Vector& pix)
   // new fitschan
   AstFitsChan* chan = astFitsChan(NULL, NULL, "");
   if (!astOK || chan == AST__NULL)
-    return NULL;
+    return;
 
   // no warning messages
   astClear(chan,"Warnings");
@@ -1485,7 +1485,7 @@ void FitsImage::initWCS0(const Vector& pix)
   // do we have anything?
   if (!astOK || frameSet == AST__NULL || 
       strncmp(astGetC(frameSet,"Class"), "FrameSet", 8))
-    return NULL;
+    return;
 
   astShow(frameSet);
 
@@ -1494,7 +1494,7 @@ void FitsImage::initWCS0(const Vector& pix)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::initWCSPhysical()
 {
   // now see if we have a 'physical' wcs, if so, set LTMV keywords
@@ -1573,7 +1573,7 @@ void FitsImage::load()
   data_ = analysisdata_;
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::match(const char* xxname1, const char* yyname1,
 		      Coord::CoordSystem sys1, Coord::SkyFrame sky1,
 		      const char* xxname2, const char* yyname2,
@@ -2481,7 +2481,7 @@ void FitsImage::resetWCS()
 	  image_->primary() && image_->inherit() ? image_->primary() : NULL);
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::resetWCS0()
 {
   int ii = Coord::WCS0-Coord::WCS;
@@ -3050,7 +3050,7 @@ void FitsImage::updatePS(Matrix3d ps)
 
 // WCS
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 Vector FitsImage::getWCScdelt(Coord::CoordSystem sys)
 {
   if (hasWCS(sys)) {
@@ -3080,7 +3080,7 @@ Vector FitsImage::getWCScdelt(Coord::CoordSystem sys)
 }
 #endif
 
-#ifdef NEWWCS
+#ifndef OLDWCS
 double FitsImage::getWCSPixelSize(Coord::CoordSystem sys)
 {
   if (!hasWCS(sys))
@@ -3121,7 +3121,7 @@ double FitsImage::getWCSPixelArea(Coord::CoordSystem sys)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 Coord::Orientation FitsImage::getWCSOrientation(Coord::CoordSystem sys,
 						Coord::SkyFrame sky)
 {
@@ -3180,7 +3180,7 @@ Coord::Orientation FitsImage::getWCSOrientation(Coord::CoordSystem sys,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double FitsImage::getWCSRotation(Coord::CoordSystem sys, Coord::SkyFrame sky)
 {
   if (hasWCS(sys)) {
@@ -3250,7 +3250,7 @@ double FitsImage::getWCSRotation(Coord::CoordSystem sys, Coord::SkyFrame sky)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 const char* FitsImage::getWCSName(Coord::CoordSystem sys) 
 {
   return (wcs_ && wcs_[sys-Coord::WCS]) ? 
@@ -3266,7 +3266,7 @@ const char* FitsImage::getWCSName(Coord::CoordSystem sys)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 Vector FitsImage::pix2wcs(const Vector& in, Coord::CoordSystem sys,
 			  Coord::SkyFrame sky)
 {
@@ -3301,7 +3301,7 @@ Vector FitsImage::pix2wcs(const Vector& in, Coord::CoordSystem sys,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 char* FitsImage::pix2wcs(const Vector& in, Coord::CoordSystem sys,
 			 Coord::SkyFrame sky, Coord::SkyFormat format,
 			 char* lbuf)
@@ -3430,7 +3430,7 @@ char* FitsImage::pix2wcs(const Vector& in, Coord::CoordSystem sys,
 }
 #endif
 
-#ifdef NEWWCS
+#ifndef OLDWCS
 Vector3d FitsImage::pix2wcs(const Vector3d& in, Coord::CoordSystem sys,
 			    Coord::SkyFrame sky)
 {
@@ -3512,7 +3512,7 @@ char* FitsImage::pix2wcs(const Vector3d& in, Coord::CoordSystem sys,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 Vector FitsImage::wcs2pix(const Vector& vv, Coord::CoordSystem sys,
 			  Coord::SkyFrame sky)
 {
@@ -3566,7 +3566,7 @@ Vector3d FitsImage::wcs2pix(const Vector3d& vv, Coord::CoordSystem sys,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double FitsImage::getWCSDist(const Vector& vv1, const Vector& vv2,
 			     Coord::CoordSystem sys)
 {
@@ -3596,7 +3596,7 @@ double FitsImage::getWCSDist(const Vector& vv1, const Vector& vv2,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 int FitsImage::hasWCS(Coord::CoordSystem sys)
 {
   int ss = sys-Coord::WCS;
@@ -3673,7 +3673,7 @@ int FitsImage::hasWCSCel(Coord::CoordSystem sys)
 
 // WCSX
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 
 int FitsImage::hasWCS3D(Coord::CoordSystem sys)
 {
@@ -3714,7 +3714,7 @@ int FitsImage::hasWCS3D(Coord::CoordSystem sys)
 
 // WCS/AST support
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::wcsShow(WorldCoor* ww)
 {
   if (!ww)
@@ -3980,7 +3980,7 @@ void FitsImage::wcsHPXInit()
 
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::astinit0(int ss, FitsHead* hd, FitsHead* prim)
 {
   if (!wcs_[ss]) {
@@ -4013,7 +4013,7 @@ int FitsImage::checkWCS(Vector3d& vv)
 	  fabs(vv[2]) < FLT_MAX ) ? 1 : 0;
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::setWCSFormat(AstFrameSet* aa, int id, const char* format)
 {
   // is it already set?
@@ -4049,7 +4049,7 @@ void FitsImage::setWCSFormat(int id, const char* format)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FitsImage::setWCSSkyFrame(AstFrameSet* ast, Coord::SkyFrame sky)
 {
   // is sky frame
@@ -4194,7 +4194,7 @@ void FitsImage::setWCSSkyFrame(Coord::CoordSystem sys, Coord::SkyFrame sky)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 int FitsImage::wcsIsASkyFrame(AstFrameSet* ast)
 {
   astClearStatus;
@@ -4205,7 +4205,7 @@ int FitsImage::wcsIsASkyFrame(AstFrameSet* ast)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 Vector FitsImage::wcsTran(AstFrameSet* ast, const Vector& in, int forward)
 {
   double xout, yout;
@@ -4433,7 +4433,7 @@ Vector3d FitsImage::wcsTran(const Vector3d& in, int forward)
 
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double FitsImage::wcsDistance(AstFrameSet* ast, const Vector& vv1,
 			      const Vector& vv2)
 {
@@ -4484,7 +4484,7 @@ double FitsImage::wcsDistance(const Vector& vv1, const Vector& vv2)
 
 #endif
 
-#ifdef NEWWCS
+#ifndef OLDWCS
 double FitsImage::wcsAngle(const Vector& vv1, const Vector& vv2,
 			   const Vector& vv3)
 {
@@ -4580,7 +4580,7 @@ double FitsImage::wcsAxAngle(const Vector& vv1, const Vector& vv2)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 AstFrameSet* FitsImage::fits2ast(FitsHead* hd) 
 {
   // we may have an error, just reset
@@ -4809,7 +4809,7 @@ AstFrameSet* FitsImage::fits2ast(FitsHead* hd)
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 AstFrameSet* FitsImage::buildast(int ss, FitsHead* hd, FitsHead* prim) 
 {
   if (DebugAST)
diff --git a/tksao/frame/fitsimage.h b/tksao/frame/fitsimage.h
index 29f9013..f1fdeda 100644
--- a/tksao/frame/fitsimage.h
+++ b/tksao/frame/fitsimage.h
@@ -111,7 +111,7 @@ class FitsImage {
   int address[FTY_MAXAXES];
 
   int manageWCS_;
-#ifndef NEWWCS
+#ifdef OLDWCS
   WorldCoor** wcs_;    // wcs list
   WCSx** wcsx_;        // xth Axis WCS
 #else
@@ -128,7 +128,7 @@ class FitsImage {
   Matrix wcsToRef_;          // iraf/wcs matrix
 
  public:
-#ifndef NEWWCS
+#ifdef OLDWCS
   AstFrameSet** ast_;  // ast frameset;
 #else
   AstFrameSet* ast_;  // ast frameset;
@@ -151,7 +151,7 @@ class FitsImage {
   void initHPX();
 
   void initWCS(FitsHead*, FitsHead*);
-#ifndef NEWWCS
+#ifdef OLDWCS
   void initWCSPhysical();
   void wcsShow(WorldCoor*);
   void astinit(int, FitsHead*, FitsHead*);
@@ -385,7 +385,7 @@ class FitsImage {
 
   char* pix2wcs(const Vector&, Coord::CoordSystem, Coord::SkyFrame, Coord::SkyFormat, char*);
 
-#ifndef NEWWCS
+#ifdef OLDWCS
   WCSx** wcsx() {return wcsx_;}
   double pix2wcsx(double, Coord::CoordSystem);
   double wcs2pixx(double, Coord::CoordSystem);
@@ -416,7 +416,7 @@ class FitsImage {
   double getWCSDist(const Vector&, const Vector&, Coord::CoordSystem);
   const char* getWCSName(Coord::CoordSystem);
 
-#ifndef NEWWCS
+#ifdef OLDWCS
   WorldCoor* getWCS(Coord::CoordSystem sys) 
   {return (wcs_ && wcs_[sys-Coord::WCS]) ? wcs_[sys-Coord::WCS] : NULL;}
   Vector getWCScdelt(Coord::CoordSystem);
@@ -470,7 +470,7 @@ class FitsImage {
   Vector mapFromRef(const Vector&, Coord::CoordSystem, Coord::SkyFrame =Coord::FK5);
   void mapFromRef(const Vector&, Coord::CoordSystem, Coord::SkyFrame, Coord::SkyFormat, char*);
   Vector mapToRef(const Vector&, Coord::CoordSystem, Coord::SkyFrame =Coord::FK5);
-#ifdef NEWWCS
+#ifndef OLDWCS
   Vector3d mapFromRef(const Vector3d&, Coord::CoordSystem, Coord::SkyFrame =Coord::FK5);
   Vector3d mapToRef(const Vector3d&, Coord::CoordSystem, Coord::SkyFrame =Coord::FK5);
 #endif
diff --git a/tksao/frame/fitsmap.C b/tksao/frame/fitsmap.C
index c8ab918..97f714e 100644
--- a/tksao/frame/fitsmap.C
+++ b/tksao/frame/fitsmap.C
@@ -59,7 +59,7 @@ Vector FitsImage::mapToRef(const Vector& vv, Coord::CoordSystem in,
   return Vector();
 }
 
-#ifdef NEWWCS
+#ifndef OLDWCS
 Vector3d FitsImage::mapFromRef(const Vector3d& vv, Coord::CoordSystem out,
 			       Coord::SkyFrame sky)
 {
@@ -143,7 +143,7 @@ void FitsImage::listFromRef(ostream& str, const Vector& vv,
 
 // Map Length
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double FitsImage::mapLenFromRef(double dd, Coord::CoordSystem sys,
 				Coord::DistFormat dist)
 {
@@ -247,7 +247,7 @@ Vector FitsImage::mapLenFromRef(const Vector& vv, Coord::CoordSystem sys,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double FitsImage::mapLenToRef(double dd, Coord::CoordSystem sys, 
 			      Coord::DistFormat dist)
 {
@@ -516,7 +516,7 @@ void FitsImage::listDistFromRef(ostream& str,
 
 // 3D
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double FitsImage::mapFromImage3d(double dd, Coord::CoordSystem sys)
 {
   if (sys >= Coord::WCS)
diff --git a/tksao/frame/frame3dbase.C b/tksao/frame/frame3dbase.C
index 248b495..9f9dd2f 100644
--- a/tksao/frame/frame3dbase.C
+++ b/tksao/frame/frame3dbase.C
@@ -152,7 +152,7 @@ void Frame3dBase::getInfoCmd(const Vector& vv, Coord::InternalSystem ref,
   getInfoClearValue(var);
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void Frame3dBase::getInfoWCS(char* var, Vector3d& rr, FitsImage* ptr, 
 			     FitsImage* sptr)
 {
@@ -256,7 +256,7 @@ void Frame3dBase::getInfoWCS(char* var, Vector3d& rr, FitsImage* ptr,
 }
 #endif
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void Frame3dBase::coordToTclArray(FitsImage* ptr, const Vector3d& vv, 
 				  Coord::CoordSystem out,
 				  const char* var, const char* base)
@@ -1010,7 +1010,7 @@ void Frame3dBase::updatePannerMatrices()
   Base::updatePannerMatrices();
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void Frame3dBase::updatePanner()
 {
   // do this first
diff --git a/tksao/frame/framebase.C b/tksao/frame/framebase.C
index e5a9fdb..d4be1d7 100644
--- a/tksao/frame/framebase.C
+++ b/tksao/frame/framebase.C
@@ -284,7 +284,7 @@ void FrameBase::updateBin(const Matrix& mx)
   Base::updateBin(mx);
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FrameBase::updatePanner()
 {
   Base::updatePanner();
diff --git a/tksao/frame/framergb.C b/tksao/frame/framergb.C
index cfd1ddf..b197f05 100644
--- a/tksao/frame/framergb.C
+++ b/tksao/frame/framergb.C
@@ -90,7 +90,7 @@ void FrameRGB::alignWCS(Coord::CoordSystem sys, Coord::SkyFrame sky)
   updateRGBMatrices();
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 void FrameRGB::alignWCS(FitsImage* ptr, Coord::CoordSystem sys)
 {
   if (!wcsAlign_ || !(keyContext->fits) || !ptr || 
diff --git a/tksao/frame/frblt.C b/tksao/frame/frblt.C
index c488408..acfe79f 100644
--- a/tksao/frame/frblt.C
+++ b/tksao/frame/frblt.C
@@ -372,7 +372,7 @@ int Base::markerAnalysisRadial(Marker* pp, double** x, double** y, double** e,
   int unit =0;
   double xaxis =1;
   if (ptr->hasWCS(sys)) {
-#ifndef NEWWCS
+#ifdef OLDWCS
     Vector cdelt= ptr->getWCScdelt(sys);
     double ll = fabs(cdelt[0]);
 #else
@@ -389,7 +389,7 @@ int Base::markerAnalysisRadial(Marker* pp, double** x, double** y, double** e,
     }
   }
     
-#ifndef NEWWCS
+#ifdef OLDWCS
   Vector cdelt= ptr->getWCScdelt(sys);
   double aa = fabs(cdelt[0]*cdelt[1]);
 #else
@@ -495,7 +495,7 @@ int Base::markerAnalysisPanda(Marker* pp, double** x, double** y, double** e,
   int unit =0;
   double xaxis =1;
   if (ptr->hasWCS(sys)) {
-#ifndef NEWWCS
+#ifdef OLDWCS
     Vector cdelt= ptr->getWCScdelt(sys);
     double ll = fabs(cdelt[0]);
 #else
@@ -512,7 +512,7 @@ int Base::markerAnalysisPanda(Marker* pp, double** x, double** y, double** e,
     }
   }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
   Vector cdelt= ptr->getWCScdelt(sys);
   double aa = fabs(cdelt[0]*cdelt[1]);
 #else
@@ -832,7 +832,7 @@ int Base::markerAnalysisStats1(Marker* pp,FitsImage* ptr, ostream& str,
     return 0;
   default: 
     {
-#ifndef NEWWCS
+#ifdef OLDWCS
       Vector cdelt= ptr->getWCScdelt(sys);
       double ll = fabs(cdelt[0]);
 #else
@@ -878,7 +878,7 @@ void Base::markerAnalysisStats2(FitsImage* ptr, ostream& str,
   case 1:
     {
       // Cel WCS
-#ifndef NEWWCS
+#ifdef OLDWCS
       Vector cdelt= ptr->getWCScdelt(sys);
       double aa = fabs(cdelt[0]*cdelt[1]);
 #else
@@ -890,7 +890,7 @@ void Base::markerAnalysisStats2(FitsImage* ptr, ostream& str,
   case 2:
     {
       // Linear WCS
-#ifndef NEWWCS
+#ifdef OLDWCS
       Vector cdelt= ptr->getWCScdelt(sys);
       double aa = fabs(cdelt[0]*cdelt[1]);
 #else
diff --git a/tksao/frame/frmap.C b/tksao/frame/frmap.C
index eb9f565..cf1047c 100644
--- a/tksao/frame/frmap.C
+++ b/tksao/frame/frmap.C
@@ -64,7 +64,7 @@ double Base::mapAngleToRef(double angle, Coord::CoordSystem sys,
   return zeroTWOPI(rr);
 }
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 double Base::mapLenToRef(double d, Coord::InternalSystem sys)
 {
   Vector r = mapLenToRef(Vector(d,0),sys);
diff --git a/tksao/frame/frsave.C b/tksao/frame/frsave.C
index 770c7f2..b14a9b1 100644
--- a/tksao/frame/frsave.C
+++ b/tksao/frame/frsave.C
@@ -439,7 +439,7 @@ void FrameBase::saveFitsResampleKeyword(OutFitsStream& str, FitsHead& dst)
   // WCS
   if (currentContext->fits->hasWCS(Coord::WCS)) {
 
-#ifndef NEWWCS
+#ifdef OLDWCS
     WorldCoor* wcs = currentContext->fits->getWCS(Coord::WCS);
 
     // abort if this is a DSS, ZPN, TNX
diff --git a/tksao/frame/grid25d.C b/tksao/frame/grid25d.C
index 6b46d1a..f3a698f 100644
--- a/tksao/frame/grid25d.C
+++ b/tksao/frame/grid25d.C
@@ -58,7 +58,7 @@ int Grid25d::doit(RenderMode rm)
   default:
     {
       // set desired skyformat
-#ifndef NEWWCS
+#ifdef OLDWCS
       AstFrameSet* ast = (AstFrameSet*)astCopy(fits->ast_[system_-Coord::WCS]);
       fits->setWCSSkyFrame(ast, sky_);
 #else
diff --git a/tksao/frame/grid2d.C b/tksao/frame/grid2d.C
index 6709b0c..b247648 100644
--- a/tksao/frame/grid2d.C
+++ b/tksao/frame/grid2d.C
@@ -59,7 +59,7 @@ int Grid2d::doit(RenderMode rm)
     {
 
       // set desired skyformat
-#ifndef NEWWCS
+#ifdef OLDWCS
       AstFrameSet* ast = (AstFrameSet*)astCopy(fits->ast_[system_-Coord::WCS]);
       fits->setWCSSkyFrame(ast, sky_);
 #else
diff --git a/tksao/frame/grid3d.C b/tksao/frame/grid3d.C
index 37a770e..09c3557 100644
--- a/tksao/frame/grid3d.C
+++ b/tksao/frame/grid3d.C
@@ -22,7 +22,7 @@ Grid3d::Grid3d(Widget* p, Coord::CoordSystem sys, Coord::SkyFrame sky,
 Grid3d::~Grid3d()
 {}
 
-#ifndef NEWWCS
+#ifdef OLDWCS
 static FitsImage* foobar;
 
 void bar(AstMapping* that, int npoint, int ncoord_in, const double* ptr_in[],
@@ -81,7 +81,7 @@ int Grid3d::doit(RenderMode rm)
   default:
     {
 
-#ifndef NEWWCS
+#ifdef OLDWCS
       foobar = fits;
       AstFrameSet* ast = (AstFrameSet*)astCopy(fits->ast_[system_-Coord::WCS]);
       fits->setWCSSkyFrame(ast, sky_);