[daip] PCAL and ORI- error

[Andy Biggs]

Hi Eric. Thanks - I've rerun PCAL and SPLIT and both work fine now.

Andy

On 12 Mar 2012, at 23:45, Eric Greisen wrote:

> Andy Biggs wrote:
>> Hi there. I'm reducing some 22-GHz VLBA + EB data and need to calibrate the D-terms. I can use LPCAL and that seems to work quite nicely, but the leakage is quite high on EB, OV and MK - the last two agree very nicely with values published in VLBA scientific memo 30 so I definitely believe them. For LCP of OV they approach 15% so I don't think I should be using the linearised model. I've therefore tried PCAL as that offers the 'ORI-' option - this is supposed to be better when dealing with high levels of leakage. Using OPTY='RAPR' I get similar results to LPCAL so I think I've at least got this task working correctly. If I try OPTY='ORI-' it seems to work fine, but any attempt to then use the results via e.g. SPLIT, ends in the following error message:
>>   1    5   12-MAR-2012  20:43:39     SPLIT     Task SPLIT  (release of 31DEC12) begins
>>   1    2   12-MAR-2012  20:43:39     SPLIT     You are using a non-standard program
>>   1    2   12-MAR-2012  20:43:39     SPLIT     Doing subarray  1
>>   1    3   12-MAR-2012  20:43:39     SPLIT     UVGET: Using flag table version   3 to edit data
>>   1    2   12-MAR-2012  20:43:59     SPLIT     Create 0218+357    .SPLIT .   1 (UV)  on disk  3  cno  202
>>       pc016987   AIPS (31DEC12)     218               12-MAR-2012   20:52:07               Page    2
>> Pops  Prior    Date        Time       Task          Messages for user  218
>>    1    3   12-MAR-2012  20:43:59     SPLIT     Applying CL Table version    13
>>   1    7   12-MAR-2012  20:43:59     SPLIT     ZCLOSE: DOES NOT PERFORM LINE PRINTER CLOSES AS OF 15JAN91
>>   1    6   12-MAR-2012  20:43:59     SPLIT     TABINI: I/O ERROR FROM ISTAB ON AN
>>   1    6   12-MAR-2012  20:43:59     SPLIT     ANINI: ERROR   3 FROM TABINI OPCODE = READ
>>   1    7   12-MAR-2012  20:43:59     SPLIT     ANINI: ERROR   3 INITIALIZING ANTENNA TABLE FOR READ
>>   1    8   12-MAR-2012  20:43:59     SPLIT     PDRGET: ERROR  3 OPEN-FOR-WRITE AN FILE
>>   1    4   12-MAR-2012  20:43:59     SPLIT     Destroyed  1 extension files of type NX
>>   1    4   12-MAR-2012  20:44:00     SPLIT     Destroyed UV image file: catno=    202 disk= 3
>>   1    6   12-MAR-2012  20:44:00     SPLIT     PROBLEM WITH SOURCE: 0218+357
>>   1    9   12-MAR-2012  20:44:00     SPLIT     FAILED TO SELECT VALID DATA
>>   1    3   12-MAR-2012  20:44:00     SPLIT     Purports to die of UNNATURAL causes
>>   1    5   12-MAR-2012  20:44:00     SPLIT     pc016987 31DEC12 TST: Cpu=      3.9  Real=     21  IO=         9
>> I do not get this message if I run PCAL using the RAPR option - that seems to work fine. I'm running an absolutely up-to-date version of AIPS 31DEC12.
>> By the way, I'm also confused about what DOPOL=2 does. This applies the second-order corrections from the linear model, but that name implies that this solution doesn't have second-order terms. Or is this referring to the correction of the parallel hand visibilities?
>> 
> I have found an error where a 512-word buffer should have been and instead a local 16 word variable was used in initializing the pol cal for ORI-.  DOPOL=2,3 is checked only for the special VLBI mode not offered by PCAL.  That mode can be applied only to cross-hands (DOPOL=1 or 6) or to all 4 with DOPOL=2,3,7,8.  ORI is applied to all 4 like it or not.
> 
> I will correct this code now and re-link lots of things.  Can you try it tomorrow?
> 
> Thanks
> 
> Eric