From 06d01d959a8c62d4781f19985e775017a9f00540 Mon Sep 17 00:00:00 2001
From: William Joye <wjoye@cfa.harvard.edu>
Date: Fri, 10 May 2019 14:34:01 -0400
Subject: check for AST WCS first

---
 tksao/frame/fitsimage.C | 28 +++++++++++++++-------------
 1 file changed, 15 insertions(+), 13 deletions(-)

diff --git a/tksao/frame/fitsimage.C b/tksao/frame/fitsimage.C
index 73a9e1d..dfbcb7f 100644
--- a/tksao/frame/fitsimage.C
+++ b/tksao/frame/fitsimage.C
@@ -3105,19 +3105,21 @@ void FitsImage::scanWCS(FitsHead* hd)
   astBegin;
   int nframes = astGetI(ast_, "Nframe");
 
-  // wcs_[]
-  // since we have ast_
-  wcs_[0] =1; 
-  wcsNaxes_[0] = astGetI(ast_,"Naxes");
-
-  // fill out wcs_ array
-  for (int kk=0; kk<nframes; kk++) {
-    AstFrameSet* ff = (AstFrameSet*)astGetFrame(ast_,kk+1);
-    const char* id = astGetC(ff, "Ident");
-    if (id && *id) {
-      int jj = (*id == ' ') ? 0 : *id-'@';
-      wcs_[jj] = 1;
-      wcsNaxes_[jj] = astGetI(ff,"Naxes");
+  // special case for AST native encoding
+  if (findKeyword("BEGAST_A")) {
+    wcs_[0] =1; 
+    wcsNaxes_[0] = astGetI(ast_,"Naxes");
+  }
+  else {
+    // fill out wcs_ array
+    for (int kk=0; kk<nframes; kk++) {
+      AstFrameSet* ff = (AstFrameSet*)astGetFrame(ast_,kk+1);
+      const char* id = astGetC(ff, "Ident");
+      if (id && *id) {
+	int jj = (*id == ' ') ? 0 : *id-'@';
+	wcs_[jj] = 1;
+	wcsNaxes_[jj] = astGetI(ff,"Naxes");
+      }
     }
   }
   
-- 
cgit v0.12