[daip] SNPLT Phase-cal plotting bug

[Mark Kettenis]

The plots made by SNPLT from a PC table are not quite what the labels
say they are.  If you're plotting a PC table that contains more than a
single tone per IF, and are plotting more than a single IF, the plot
labelled "Phase tone 1, IF 2" is actually for "Phase tone 2, IF 1".
The attached patch fixes that for PHAS and AMP plots.  A similar fix
might be needed for 'RATE' plots.

Thanks,

Mark


Index: SNPLT.FOR
===================================================================
RCS file: /home/aipssrc/cvs/31DEC04/Y/PGM/NOTST/SNPLT.FOR,v
retrieving revision 1.11
diff -u -p -r1.11 SNPLT.FOR
--- SNPLT.FOR	23 Nov 2004 04:47:13 -0000	1.11
+++ SNPLT.FOR	19 Jul 2005 14:36:41 -0000
@@ -1386,20 +1386,20 @@ C                                       
  100  IF (ICODE.EQ.1) THEN
          DO 110 IIS = 1,MUMPOL
             LP = IIS - MUMPOL
-            IP1 = REKOL(IIS) - 1
-            JP1 = IMKOL(IIS) - 1
+            IP1 = REKOL(IIS) - NTONE
+            JP1 = IMKOL(IIS) - NTONE
             IF (ITPLOT.GE.4) THEN
-               IP2 = REKOL(IIS+2) - 1
-               JP2 = IMKOL(IIS+2) - 1
+               IP2 = REKOL(IIS+2) - NTONE
+               JP2 = IMKOL(IIS+2) - NTONE
             ELSE
-               IP2 = REKOL(2) - 1
-               JP2 = IMKOL(2) - 1
+               IP2 = REKOL(2) - NTONE
+               JP2 = IMKOL(2) - NTONE
                END IF
             DO 105 IIF = 1,MUMIF
-               IP1 = IP1 + 1
-               JP1 = JP1 + 1
-               IP2 =  IP2 + 1
-               JP2 = JP2 + 1
+               IP1 = IP1 + NTONE
+               JP1 = JP1 + NTONE
+               IP2 =  IP2 + NTONE
+               JP2 = JP2 + NTONE
                LP = LP + MUMPOL
                IF ((GNREC(IP1).NE.FBLANK) .AND. (GNREC(JP1).NE.FBLANK))
      *            THEN
@@ -1425,20 +1425,20 @@ C                                       
  120  ELSE IF (ICODE.EQ.2) THEN
          DO 130 IIS = 1,MUMPOL
             LP = IIS - MUMPOL
-            IP1 = REKOL(IIS) - 1
-            JP1 = IMKOL(IIS) - 1
+            IP1 = REKOL(IIS) - NTONE
+            JP1 = IMKOL(IIS) - NTONE
             IF (ITPLOT.GE.4) THEN
-               IP2 = REKOL(IIS+2) - 1
-               JP2 = IMKOL(IIS+2) - 1
+               IP2 = REKOL(IIS+2) - NTONE
+               JP2 = IMKOL(IIS+2) - NTONE
             ELSE
-               IP2 = REKOL(2) - 1
-               JP2 = IMKOL(2) - 1
+               IP2 = REKOL(2) - NTONE
+               JP2 = IMKOL(2) - NTONE
                END IF
             DO 125 IIF = 1,MUMIF
-               IP1 = IP1 + 1
-               JP1 = JP1 + 1
-               IP2 =  IP2 + 1
-               JP2 = JP2 + 1
+               IP1 = IP1 + NTONE
+               JP1 = JP1 + NTONE
+               IP2 =  IP2 + NTONE
+               JP2 = JP2 + NTONE
                LP = LP + MUMPOL
                IF ((GNREC(IP1).NE.FBLANK) .AND. (GNREC(JP1).NE.FBLANK))
      *            THEN