[daip] Could you help me?

[Eric Greisen]

Sang-Sung Lee writes:

 > I have  a problem in running the task VPLOT on AIPS.
 > To check the result from LPCAL, I want to see how the REAL-vs-IMAG plot 
 > change with dopol -1/+1.
 > I could plot those two  STOKES - POLPLOT configurations ;  RR-RR/LL, 
 > LL-LL/RR.
 > But for the others, such as RL-RL/LL(RR) and LR-LR/LL(RR), I always have 
 > the message written in the file.
 > I put in the attached file my input parameter of the task and the 
 > message from its running.
 > 

I have solved the problem and put the correction into the 31DEC04
version of VPLOT.  This will be available beginning later tonight -
say after 9 hours from now - for running a MNJ.  If you do not have
the 31DEC04 version, I send you corrections to make in the form of a
diff file.

Eric Greisen

diff VPLOT.FOR old-vers
39c39
<      *   TVCHN, TVCORN(4), NUMFRQ, DTYPE, KPLOT, LABEL, INVER,
---
>      *   TVCHN, TVCORN(4), NUMFRQ, NUMPLZ, DTYPE, KPLOT, LABEL, INVER,
61c61
<      *   DTYPE, DOTV, DOERRB, MULTI, DOAVG, AVGSCA, AMPPH,
---
>      *   NUMPLZ, DTYPE, DOTV, DOERRB, MULTI, DOAVG, AVGSCA, AMPPH,
78c78
< C;  Copyright (C) 1995-2004
---
> C;  Copyright (C) 1995-2003
195d194
<       IPOLPT = 0
584a584
>       NUMPLZ = CATBLK(KINAX+JLOCS)
967,968c967,968
<             CALL AVGCIF (BUFF1, BCHAN, ECHAN, CIF, BIF, EIF, LOOP, VMIN,
<      *         VMAX, VISO)
---
>             CALL AVGCIF (BUFF1, BCHAN, ECHAN, CIF, BIF, EIF, LOOP, KNCS,
>      *         VMIN, VMAX, VISO)
1309c1309
<      *   AUNITS(9)*8, CHTYPE(17,2)*9, POLLAB*5, STRING*8
---
>      *   AUNITS(9)*8, CHTYPE(17)*8, POLLAB*5, STRING*8
1314,1315c1314
<      *   IAXL2, NGOOD, NNOFIT, NUMVIS, LTYPE, MBASE, NN, ISOLIN, J2,
<      *   LIF, LTY
---
>      *   IAXL2, NGOOD, NNOFIT, NUMVIS, LTYPE, MBASE, NN, ISOLIN, J2, LIF
1333c1332
<       DATA AUNITS /'Janskys ', 'Degrees ',  'Wavlngth', 'Degrees ',
---
>       DATA AUNITS /'Janskys ', 'Degrees ',  'Wavlngth','Degrees ',
1335c1334
<       DATA CHTYPE /'Amplitude', 'Phase   ', 'UV dist ',
---
>       DATA CHTYPE /'Ampltude', 'Phase   ', 'UV dist ',
1339,1344c1338
<      *   'Weight  ', 'Rms     ',
<      *   'RatioAmpl', 'RatioPhas', 'UV dist ',
<      *   'UV pa   ', 'Time    ', 'U       ', 'V      ',
<      *   'W       ', 'RatioReal', 'RatioImag    ', 'GST     ',
<      *   'Time    ', 'Source #', 'Freqid  ', 'Int.time',
<      *   'RatioWeit', 'Rms     '/
---
>      *   'Weight  ', 'Rms     '/
1346,1347d1339
<       LTY = 1
<       IF (POLPLT.GT.0) LTY = 2
1438,1439d1429
<          IF ((CTYP(I,LOCNUM).EQ.AUNITS(1)) .AND. (POLPLT.GT.0)) 
<      *      CTYP(I,LOCNUM) = 'Ratio'
1566,1567c1556,1557
<          TEXT = CHTYPE(TYPEAX(2),LTY)
<          INP = 10
---
>          TEXT = CHTYPE(TYPEAX(2))
>          INP = 9
1569,1570c1559,1560
<             TEXT(INP:) = ' and' // CHTYPE(TYPEAX(3),LTY)(1:8)
<             INP = INP + 13
---
>             TEXT(INP:) = ' and' // CHTYPE(TYPEAX(3))(1:8)
>             INP = INP + 12
1572,1573c1562,1563
<          TEXT(INP:) = ' vs ' // CHTYPE(TYPEAX(1),LTY)(1:8) // ' _for_ '
<          INP = INP + 20
---
>          TEXT(INP:) = ' vs ' // CHTYPE(TYPEAX(1))(1:8) // ' _for_ '
>          INP = INP + 19
1876,1877c1866,1867
<  100     TR = BUFR(1)
<          TI = BUFR(2)
---
>  100     TR = BUFR(1+IPOLPT*INCS)
>          TI = BUFR(2+IPOLPT*INCS)
1944c1934
<  890     XY(I) = BUFR(3)
---
>  890     XY(I) = BUFR(3+IPOLPT*INCS)
3142,3144d3131
< C                                       Modify VIS as appropriate
< C                                       based on POLPLT.
<          IF (POLPLT.GT.0) CALL POLVIS (VIS, POLPLT, IPOLPT)
3197c3184
<      *         VMIN, VMAX, VISO)
---
>      *         KNCS, VMIN, VMAX, VISO)
3729,3730c3716,3717
<       SUBROUTINE AVGCIF (VIS, BCHAN, ECHAN, CIF, BIF, EIF, IBAS, VMIN,
<      *   VMAX, VISOUT)
---
>       SUBROUTINE AVGCIF (VIS, BCHAN, ECHAN, CIF, BIF, EIF, IBAS, JNCS,
>      *   VMIN, VMAX, VISOUT)
3743a3731
> C      JNCS     I         Index for polzn in output array
3750c3738
< C      VISOUT   R(3)      Pseudo-continuum visibility
---
> C      VISOUT   R(*)      Pseudo-continuum visibility
3752c3740
<       INTEGER   BCHAN, ECHAN, CIF, BIF, EIF, IBAS
---
>       INTEGER   BCHAN, ECHAN, CIF, BIF, EIF, IBAS, JNCS
3757c3745,3746
<       INTEGER   LOOPIF, LOOPF, INDEX, INP, FCHAN, SCHAN, FIF, SIF
---
>       INTEGER   LOOPS, LOOPIF, LOOPF, INDEX, INP, OUTDEX, FCHAN, SCHAN,
>      *   FIF, SIF
3772a3762
>       NUMPLZ = 1
3773a3764,3765
>       DO 40 LOOPS = 1,NUMPLZ
>          OUTDEX = 1 + (LOOPS-1)*JNCS
3783c3775
<          INDEX = 1 + IPOLPT*INCS + (LOOPIF-BIF)*INCIF
---
>             INDEX = 1 + (LOOPS-1)*INCS + (LOOPIF-BIF)*INCIF
3824a3817,3819
>             ELSE
>                GO TO 30
>                END IF
3855d3849
<                END IF
3865,3866c3859,3860
<          VISOUT(1) = AMP * COS(FI)
<          VISOUT(2) = AMP * SIN(FI)
---
>             VISOUT(OUTDEX) = AMP * COS(FI)
>             VISOUT(OUTDEX+1) = AMP * SIN(FI)
3868,3869c3862,3863
<          VISOUT(1) = RE
<          VISOUT(2) = IM
---
>             VISOUT(OUTDEX) = RE
>             VISOUT(OUTDEX+1) = IM
3871c3865,3866
<       VISOUT(3) = SUMWT
---
>          VISOUT(OUTDEX+2) = SUMWT
>  40      CONTINUE