polpack: Correct indexing of WCS frames in polimage

Starlink · Sep 28, 2012 · 3bb0366 · 3bb0366
1 parent 1f87f07
commit 3bb0366
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/applications/polpack/polimage.f b/applications/polpack/polimage.f
@@ -177,6 +177,8 @@ SUBROUTINE POLIMAGE( STATUS )
 *        Correct binning of pixel coords.
 *     13-JUL-2009 (DSB):
 *        Changed IPVAR array from DOUBLE PRECISION to REAL.
+*     28-SEP-2012 (DSB):
+*        Correct indexing of WCS Frames when removing the AXIS Frame.
 *     {enter_further_changes_here}
 
 *  Bugs:
@@ -669,11 +671,15 @@ SUBROUTINE POLIMAGE( STATUS )
 
 *  Remove the AXIS Frame in the output NDF default FrameSet, so that any
 *  AXIS Frame in the catalogue WCS FrameSet will be retained.
-            DO I = 1, AST_GETI( IWCSO, 'NFRAME', STATUS )
+            I = 1
+            DO WHILE( I .LE.  AST_GETI( IWCSO, 'NFRAME', STATUS ) .AND.
+     :                STATUS .EQ. SAI__OK )
                FRM = AST_GETFRAME( IWCSO, I, STATUS )
                IF( AST_GETC( FRM, 'DOMAIN', STATUS ) .EQ.
      :                       'AXIS' ) THEN
                   CALL AST_REMOVEFRAME( IWCSO, I, STATUS )
+               ELSE
+                  I = I + 1
                END IF
             END DO