[daip] Fwd: S2053A and S2053C data
Amy Mioduszewski
amiodusz at nrao.edu
Wed Sep 15 16:43:02 EDT 2010
Hi Svetlana,
You will be thrilled to hear (I am sure) that I managed to reproduce your
problem using a Mac that had 2008 on it (which is pretty much what you have
since you installed in Dec 08), the frozen 2009 version and current 2010 work
fine, so you should update your AIPS and it will work.
Amy
Svetlana Jorstad wrote:
> Thanks, Amy.
> I am sure that I am running 31Dec 2009 installed in Dec 2008.
> I used APCAL without opacities but it did not help.
> I will try another version of AIPS.
> Thank you very much,
> Svetlana
> Amy Mioduszewski wrote:
>> Hi Svetlana,
>>
>> The file size I get is 4.6Gb. But the header I get is identical to
>> yours, i.e., same number of visibilities. Are you sure you are using
>> 31Dec2009 installed in Dec *2008*?
>>
>> I have tried and tried to reproduce your error, but cannot. I checked
>> your AN and WX tables and they look O.K.. It does look like something
>> is going wrong in you opacity fitting, could you try running APCAL
>> without opacities? Other than that I don't know what to tell you
>> except to update your AIPS. I will TASAV the tables I got by just
>> running VLBACALA so you can download and look at them. Go to
>> ftp://ftp.aoc.nrao.edu/pub/staff/amiodusz/svetlana/ and copy
>> svetlana.tables.fits.
>>
>> Thanks!
>>
>> Amy
>>
>> Svetlana Jorstad wrote:
>>> Hi Amy,
>>> I downloaded S2053C data on September 3, the size of the file is
>>> 4.28GB. I use AIPS 2009 installed in December 2008 at MAC OSX 10.5.8.
>>> I can not understand what is going on. I will ask Alan (Marscher) to get
>>> the data
>>> and repeat VLBACALA, he installed AIPS recently.
>>> Thank you very much,
>>> Svetlana
>>> The steps are:
>>> I. FITLD:
>>> localh> FITLD1: Task FITLD (release of 31DEC09) begins
>>> localh> FITLD1: Create S2053 .RAW . 1 (UV) on disk 1
>>> cno 130
>>> localh> FITLD1: Warning: table type OB is of zero length
>>> localh> FITLD1: Image=MULTI (UV) Filename=S2053
>>> .RAW .
>>> 1
>>> localh> FITLD1: Telescope=VLBA Receiver=
>>> localh> FITLD1: Observer=S2053 User #= 107
>>> localh> FITLD1: Observ. date=25-OCT-2009 Map date=13-SEP-2010
>>> localh> FITLD1: # visibilities 1470501 Sort order TB
>>> localh> FITLD1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
>>> localh> FITLD1: SOURCE FREQSEL INTTIM GATEID CORR-ID
>>> localh> FITLD1:
>>> ----------------------------------------------------------------
>>> localh> FITLD1: Type Pixels Coord value at Pixel Coord
>>> incr Rotat
>>> localh> FITLD1: COMPLEX 3 1.0000000E+00 1.00
>>> 1.0000000E+00 0.00
>>> localh> FITLD1: STOKES 4 -1.0000000E+00 1.00
>>> -1.0000000E+00 0.00
>>> localh> FITLD1: FREQ 16 4.3201490E+10 0.56
>>> 5.0000000E+05 0.00
>>> localh> FITLD1: IF 4 1.0000000E+00 1.00
>>> 1.0000000E+00 0.00
>>> localh> FITLD1: RA 1 00 00 00.000 1.00
>>> 0.000000 0.00
>>> localh> FITLD1: DEC 1 00 00 00.000 1.00
>>> 0.000000 0.00
>>> localh> FITLD1:
>>> ----------------------------------------------------------------
>>> localh> FITLD1: Coordinate equinox 2000.00
>>> localh> FITLD1: Maximum version number of extension files of type HI
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type CL
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type AT
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type NX
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type IM
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type CT
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type MC
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type OB
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type GC
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type WX
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type TY
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type PC
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type FQ
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type AN
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type SU
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type CQ
>>> is 1
>>> localh> FITLD1: Maximum version number of extension files of type FG
>>> is 1
>>> localh> FITLD1: Appears to have ended successfully
>>> localh> FITLD1: localhos 31DEC09 TST: Cpu= 53.9 Real= 538
>>> IO= 4413
>>> II. VLBAMCAL
>>>> inp VLBAMCAL
>>> AIPS 1: VLBAMCAL: Procedure to merge redundant calibration data
>>> AIPS 1: Adverbs Values Comments
>>> AIPS 1: ----------------------------------------------------------------
>>> AIPS 1: INNAME 'S2053' Input file name
>>> AIPS 1: INCLASS 'RAW' Input file class
>>> AIPS 1: INSEQ 1 Input file sequence number
>>> AIPS 1: INDISK 1 Disk number for input file
>>> AIPS 1: BADDISK *all 0 List of disks not to be used
>>> AIPS 1: for scratch files.
>>> AIPS 1:
>>> AIPS 1: VLBAMCAL is defined in the VLBAUTIL run file.
>>> localh> TACOP1: Task TACOP (release of 31DEC09) begins
>>> localh> TACOP1: Copied GC file from vol/cno/vers 1 130 1 to 1
>>> 130 2
>>> localh> TACOP1: TACOP / Copied 1 extension files
>>> localh> TACOP1: Appears to have ended successfully
>>> localh> TACOP1: localhost 31DEC09 TST: Cpu= 0.0 Real= 0
>>> localh> TAMRG1: Task TAMRG (release of 31DEC09) begins
>>> localh> TAMRG1: Start sort to scratch version 3
>>> localh> TAMRG1: Start merge to output version 1
>>> localh> TAMRG1: Read 10 rows from GC file version 2
>>> localh> TAMRG1: Wrote 10 rows to GC file version 1
>>> localh> TAMRG1: Appears to have ended successfully
>>> localh> TAMRG1: localhost 31DEC09 TST: Cpu= 0.0 Real= 0
>>> localh> TACOP1: Task TACOP (release of 31DEC09) begins
>>> localh> TACOP1: Copied PC file from vol/cno/vers 1 130 1 to 1
>>> 130 2
>>> localh> TACOP1: TACOP / Copied 1 extension files
>>> localh> TACOP1: Appears to have ended successfully
>>> localh> TACOP1: localhos 31DEC09 TST: Cpu= 0.0 Real= 0
>>> IO= 4
>>> localh> TAMRG1: Task TAMRG (release of 31DEC09) begins
>>> localh> TAMRG1: Start sort to scratch version 3
>>> localh> TAMRG1: Start merge to output version 1
>>> localh> TAMRG1: Read 3574 rows from PC file version 2
>>> localh> TAMRG1: Wrote 3574 rows to PC file version 1
>>> localh> TAMRG1: Appears to have ended successfully
>>> localh> TAMRG1: localhos 31DEC09 TST: Cpu= 0.1 Real= 0
>>> IO= 16
>>> localh> TACOP1: Task TACOP (release of 31DEC09) begins
>>> localh> TACOP1: Copied TY file from vol/cno/vers 1 130 1 to 1
>>> 130 2
>>> localh> TACOP1: TACOP / Copied 1 extension files
>>> localh> TACOP1: Appears to have ended successfully
>>> localh> TACOP1: localhos 31DEC09 TST: Cpu= 0.0 Real= 0
>>> IO= 1
>>> localh> TAMRG1: Task TAMRG (release of 31DEC09) begins
>>> localh> TAMRG1: Start sort to scratch version 3
>>> localh> TAMRG1: Start merge to output version 1
>>> localh> TAMRG1: Read 4968 rows from TY file version 2
>>> localh> TAMRG1: Wrote 4968 rows to TY file version 1
>>> localh> TAMRG1: Appears to have ended successfully
>>> localh> TAMRG1: localhos 31DEC09 TST: Cpu= 0.0 Real= 0
>>> IO= 5
>>> AIPS 1: Resumes
>>> AIPS 1: Extension file type GC version 1 deleted
>>> AIPS 1: Resumes
>>> AIPS 1: Extension file type GC version 2 deleted
>>> AIPS 1: Resumes
>>> AIPS 1: Extension file type PC version 1 deleted
>>> AIPS 1: Resumes
>>> AIPS 1: Extension file type PC version 2 deleted
>>> AIPS 1: Resumes
>>> AIPS 1: Extension file type TY version 1 deleted
>>> AIPS 1: Resumes
>>> AIPS 1: Extension file type TY version 2 deleted
>>> III. VLBACALA
>>>> inp VLBACALA
>>> AIPS 1: VLBACALA: Applies a-priori amplitude and sampling corrections
>>> AIPS 1: Adverbs Values Comments
>>> AIPS 1: ----------------------------------------------------------------
>>> AIPS 1: INNAME 'S2053' Input file name
>>> AIPS 1: INCLASS 'RAW' Input file class
>>> AIPS 1: INSEQ 1 Input file sequence number
>>> AIPS 1: INDISK 1 Disk number for input file
>>> AIPS 1: FREQID -1 Frequency ID to correct
>>> AIPS 1: 0 => 1
>>> AIPS 1: SUBARRAY 0 Subarray to correct. 0 => 1
>>> AIPS 1: DOFIT 1 1 Fit for opacity correction
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: 1 1 1 1
>>> AIPS 1: BADDISK *all 0 Disks to avoid for scratch
>>> AIPS 1:
>>> AIPS 1: VLBACALA is defined in the VLBAUTIL run file.
>>> localh> ACCOR1: Task ACCOR (release of 31DEC09) begins
>>> localh> ACCOR1: Writing to SN table 1
>>> localh> ACCOR1: Appears to have ended successfully
>>> localh> ACCOR1: localhos 31DEC09 TST: Cpu= 5.3 Real= 129
>>> IO= 4372
>>> localh> SNSMO1: Task SNSMO (release of 31DEC09) begins
>>> localh> SNSMO1: Copied SN file from vol/cno/vers 1 130 1 to 1
>>> 130 2
>>> localh> SNSMO1: SNSMTH: Clipping SN table
>>> localh> SNSMO1: SNREF: re-referencing 100.00% of table entries to
>>> antenna 1
>>> localh> SNSMO1: SNSMTH: Smoothing SN table
>>> localh> SNSMO1: Appears to have ended successfully
>>> localh> SNSMO1: localhos 31DEC09 TST: Cpu= 0.5 Real= 3
>>> IO= 71
>>> localh> TACOP1: Task TACOP (release of 31DEC09) begins
>>> localh> TACOP1: Copied SN file from vol/cno/vers 1 130 2 to 1
>>> 130 1
>>> localh> TACOP1: TACOP / Copied 1 extension files
>>> localh> TACOP1: Appears to have ended successfully
>>> localh> TACOP1: localhos 31DEC09 TST: Cpu= 0.0 Real= 1
>>> IO= 2
>>> localh> CLCAL1: Task CLCAL (release of 31DEC09) begins
>>> localh> CLCAL1: Using interpolation mode 2PT
>>> localh> CLCAL1: Processing SN table 1
>>> localh> CLCAL1: SNMRG: Merging SN table
>>> localh> CLCAL1: SNMRG: Write 4537 merged records from 4537 input
>>> records
>>> localh> CLCAL1: SN2CL: Applying SN tables to CL table 1, writing CL
>>> table 2
>>> localh> CLCAL1: Appears to have ended successfully
>>> localh> CLCAL1: localhos 31DEC09 TST: Cpu= 1.4 Real= 14
>>> IO= 318
>>> localh> APCAL1: Task APCAL (release of 31DEC09) begins
>>> localh> APCAL1: TXTWX: Weather data read from WX table 1
>>> localh> APCAL1: DFTAU: INITIAL TAU0 SET TO 0.100
>>> localh> APCAL1: DFTRC:BR RCP Initial Trec set to = 61.40
>>> localh> APCAL1: DFTRC:BR LCP Initial Trec set to = 85.45
>>> localh> APCAL1: DFTRC:FD RCP Initial Trec set to = 60.99
>>> localh> APCAL1: DFTRC:FD LCP Initial Trec set to = 65.90
>>> localh> APCAL1: DFTRC:HN RCP Initial Trec set to = 106.93
>>> localh> APCAL1: DFTRC:HN LCP Initial Trec set to = 98.66
>>> localh> APCAL1: DFTRC:KP RCP Initial Trec set to = 72.36
>>> localh> APCAL1: DFTRC:KP LCP Initial Trec set to = 67.82
>>> localh> APCAL1: DFTRC:LA RCP Initial Trec set to = 65.64
>>> localh> APCAL1: DFTRC:LA LCP Initial Trec set to = 89.13
>>> localh> APCAL1: DFTRC:MK RCP Initial Trec set to = 42.98
>>> localh> APCAL1: DFTRC:MK LCP Initial Trec set to = 56.20
>>> localh> APCAL1: DFTRC:NL RCP Initial Trec set to = 114.53
>>> localh> APCAL1: DFTRC:NL LCP Initial Trec set to = 108.27
>>> localh> APCAL1: DFTRC:OV RCP Initial Trec set to = 70.44
>>> localh> APCAL1: DFTRC:OV LCP Initial Trec set to = 74.91
>>> localh> APCAL1: DFTRC:PT RCP Initial Trec set to = 63.80
>>> localh> APCAL1: DFTRC:PT LCP Initial Trec set to = 63.05
>>> localh> APCAL1: DFTRC:SC RCP Initial Trec set to = 104.14
>>> localh> APCAL1: DFTRC:SC LCP Initial Trec set to = 92.09
>>> localh> APCAL1: Writing SN table 2
>>> localh> APCAL1: BR RCP 1/ 7h 0m Trec (K): 95.51 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: BR LCP 1/ 7h 0m Trec (K): 120.84 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: Plot file version 1 created
>>> localh> APCAL1: FD RCP 1/ 7h 0m Trec (K): 96.11 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: FD LCP 1/ 7h 0m Trec (K): 103.84 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: GFINIS: number records used 84
>>> localh> APCAL1: HN RCP 1/ 7h 0m Trec (K): 133.93 Zen. opac.:
>>> 0.002
>>> localh> APCAL1: HN LCP 1/ 7h 0m Trec (K): 125.56 Zen. opac.:
>>> 0.002
>>> localh> APCAL1: Plot file version 2 created
>>> localh> APCAL1: KP RCP 1/ 7h 0m Trec (K): 97.94 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: KP LCP 1/ 7h 0m Trec (K): 93.17 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: GFINIS: number records used 65
>>> localh> APCAL1: LA RCP 1/ 7h 0m Trec (K): 108.74 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: LA LCP 1/ 7h 0m Trec (K): 133.24 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: Plot file version 3 created
>>> localh> APCAL1: MK RCP 1/ 7h 0m Trec (K): 69.46 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: MK LCP 1/ 7h 0m Trec (K): 80.61 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: GFINIS: number records used 99
>>> localh> APCAL1: NL RCP 1/ 7h 0m Trec (K): 198.99 Zen. opac.:
>>> 0.047
>>> localh> APCAL1: NL LCP 1/ 7h 0m Trec (K): 107.17 Zen. opac.:
>>> 0.171
>>> localh> APCAL1: Plot file version 4 created
>>> localh> APCAL1: OV RCP 1/ 7h 0m Trec (K): 98.18 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: OV LCP 1/ 7h 0m Trec (K): 104.43 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: GFINIS: number records used 72
>>> localh> APCAL1: PT RCP 1/ 7h 0m Trec (K): 90.23 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: PT LCP 1/ 7h 0m Trec (K): 89.17 Zen. opac.:
>>> 0.000
>>> localh> APCAL1: Plot file version 5 created
>>> localh> APCAL1: SC RCP 1/ 7h 0m Trec (K): 44.18 Zen. opac.:
>>> 0.201
>>> localh> APCAL1: SC LCP 1/ 7h 0m Trec (K): 27.91 Zen. opac.:
>>> 0.201
>>> localh> APCAL1: GFINIS: number records used 45
>>> localh> APCAL1: *** At least one source had identically zero flux
>>> localh> APCAL1: *** Usually this means you forgot to run SETJY.
>>> localh> APCAL1: *** This may be important if you have a very
>>> localh> APCAL1: *** bright source (masers?).
>>> localh> APCAL1: Appears to have ended successfully
>>> localh> APCAL1: localhos 31DEC09 TST: Cpu= 10.1 Real= 11
>>> IO= 24
>>> localh> CLCAL1: Task CLCAL (release of 31DEC09) begins
>>> localh> CLCAL1: Using interpolation mode SELF
>>> localh> CLCAL1: Processing SN table 2
>>> localh> CLCAL1: SNMRG: Merging SN table
>>> localh> CLCAL1: SNMRG: Write 1476 merged records from 1476 input
>>> records
>>> localh> CLCAL1: SN2CL: Applying SN tables to CL table 2, writing CL
>>> table 3
>>> localh> CLCAL1: Appears to have ended successfully
>>> localh> CLCAL1: localhos 31DEC09 TST: Cpu= 3.5 Real= 11
>>> IO= 556
>>>
>>>
>>>> Hi Svetlana,
>>>>
>>>> So I can't reproduce it even with your tables (i.e., I can see it in
>>>> your
>>>> tables
>>>> but when I run APCAL it doesn't get rid of everything but MK for
>>>> OJ287).
>>>> I also
>>>> downloaded your data from the archive and don't reproduce it. I
>>>> looked in
>>>> the
>>>> history table of the data you sent me but it is very confused, so could
>>>> you tell
>>>> me the exact steps you took. I just ran VLBACALA and VLBPANG and
>>>> everything was
>>>> fine. I am somewhat at a loss, because you have tried downloading the
>>>> data
>>>> again and both forms of the archive data and tried another computer
>>>> (right?).
>>>> Tell me the exact steps and we will see if there is a clue there. The
>>>> only clue
>>>> I have is when I run VLBACALA I get messages that say:
>>>> APCAL5: BR RCP 1/ 7h 0m Trec (K): 78.14 Zen. opac.: 0.069
>>>> APCAL5: BR LCP 1/ 7h 0m Trec (K): 93.22 Zen. opac.: 0.106
>>>> APCAL5: Plot file version 16 created
>>>> APCAL5: GFINIS: number records used 210
>>>> APCAL5: Plot file version 17 created
>>>> APCAL5: GFINIS: number records used 201
>>>> APCAL5: Plot file version 18 created
>>>> APCAL5: GFINIS: number records used 194
>>>> APCAL5: Plot file version 19 created
>>>> APCAL5: GFINIS: number records used 201
>>>> APCAL5: Plot file version 20 created
>>>> APCAL5: GFINIS: number records used 195
>>>> APCAL5: *** At least one source had identically zero flux
>>>> APCAL5: *** Usually this means you forgot to run SETJY.
>>>> APCAL5: *** This may be important if you have a very
>>>> APCAL5: *** bright source (masers?).
>>>>
>>>> which are different from yours. How old is your AIPS again?
>>>>
>>>> Amy
>>>>
>>>>
>>>> Svetlana Jorstad wrote:
>>>>> Hi Amy,
>>>>> the last time I download data a week ago, B-epoch was good (I
>>>>> reduced it
>>>>> without
>>>>> problems) but A and C epochs have the same problem.
>>>>> cheers,
>>>>> Svetlana
>>>>>> Hi Sveltana,
>>>>>>
>>>>>> When did you download your data and did you download the slightly
>>>>>> processed
>>>>>> version? There was an error (fixed in the archive Jan 2010) with
>>>>>> FITLD
>>>>>> that
>>>>>> messed up the AN table, which could have the symptoms you describe.
>>>>>> And
>>>>>> it
>>>>>> effected your A, B and C data.
>>>>>>
>>>>>> Amy
>>>>>>
>>>>>> Svetlana Jorstad wrote:
>>>>>>> Hi Amy,
>>>>>>> thank you very much for the quick response,
>>>>>>> it does not matter that DOFIT is 1,0. I used 1 for all
>>>>>>> antennas and 0 for all antennas it does not change
>>>>>>> the result.
>>>>>>> >How do you know the data is gone?
>>>>>>> I used to after VLBAPANG (CL table 4) look at phases of a couple
>>>>>>> of sources to check that they behavior correctly using
>>>>>>> SNPLT. I looked at BLLac - there are all baselines and everything
>>>>>>> ok,
>>>>>>> but nothing for OJ287 in S2053C. I looked at amplitudes using CL
>>>>>>> table
>>>>>>> 3, BLLac has all antennas, OJ287 has only MK (very strange).
>>>>>>>
>>>>>>> >Did you do binary install or did you compile AIPS yourself?
>>>>>>> I used binary.
>>>>>>>
>>>>>>> >Could you do a TASAV and write out that file and put it
>>>>>>> somewhere I
>>>>>>> >can download it (make sure the CL table that you produced ).
>>>>>>> You can find the file S2053C.TABS from
>>>>>>> http://www.bu.edu/blazars/filestodownload/.
>>>>>>>
>>>>>>> Really, I was hesitate to ask because the thing is very strange
>>>>>>> that makes me think that I have a stupid mistake somewhere.
>>>>>>> I try to find it but I failed for now. In general neither VLBAMCAL
>>>>>>> or VLBACALA require many parameters where you can make
>>>>>>> mistakes. In addition,
>>>>>>> we have a paper preparing, for which these epochs are very important
>>>>>>> for OJ287 and I need the data urgently.
>>>>>>> Thank you very much,
>>>>>>> Svetlana
>>>>>>>
>>>>>>> Amy Mioduszewski wrote:
>>>>>>>> Hi Svetlana,
>>>>>>>>
>>>>>>>> Paul Dyer downloaded your data and we looked at it and
>>>>>>>> unfortunately
>>>>>>>> cannot reproduce your problem (i.e., for S2053A, OJ287 is still
>>>>>>>> there
>>>>>>>> and looks fine after running VLBACALA). Your APCAL output does
>>>>>>>> look
>>>>>>>> odd, but I suspect that is because you didn't set DOFIT=1,0 rather
>>>>>>>> than DOFIT 1 (which fills all the values in as 1). DOFIT tells
>>>>>>>> APCAL
>>>>>>>> which antennas to fit and by setting DOFIT 1, 0 you are telling
>>>>>>>> it to
>>>>>>>> fit only the first antenna (BR). But that does not solve your
>>>>>>>> problem, because even when we ran VLBACALA with DOFIT 1,0 your
>>>>>>>> sources
>>>>>>>> were still there. We also tried version 21DEC09 and got the same
>>>>>>>> result.
>>>>>>>>
>>>>>>>> So I now have a couple questions for you: How do you know the
>>>>>>>> data is
>>>>>>>> gone? Did you do binary install or did you compile AIPS yourself?
>>>>>>>> Also since this sounds like a table issue and not a data issue
>>>>>>>> could
>>>>>>>> you do a TASAV and write out that file and put it somewhere I can
>>>>>>>> download it (make sure the CL table that you produced that made the
>>>>>>>> data go away is there).
>>>>>>>>
>>>>>>>> Thanks,
>>>>>>>>
>>>>>>>> Amy
>>>>>>>>
>>>>>>>> ------------------------------------------------------------------------
>>>>>>>>
>>>>>>>>> Subject:
>>>>>>>>> S2053A and S2053C data
>>>>>>>>> From:
>>>>>>>>> "Svetlana Jorstad" <jorstad at bu.edu>
>>>>>>>>> Date:
>>>>>>>>> Fri, 3 Sep 2010 00:43:34 -0400 (EDT)
>>>>>>>>> To:
>>>>>>>>> "Data Analysts" <analysts at aoc.nrao.edu>
>>>>>>>>>
>>>>>>>>> To:
>>>>>>>>> "Data Analysts" <analysts at aoc.nrao.edu>
>>>>>>>>>
>>>>>>>>>
>>>>>>>>> It appears that the VLBA data of programs sS2053A and S2053C,
>>>>>>>>> on which I am the PI, are corrupted. The observations
>>>>>>>>> were performed on 2009 October 14 and 25.
>>>>>>>>> The SN and CL tables after running the VLBACALA utility lose
>>>>>>>>> almost all data for sources 10,11,12, and 13,
>>>>>>>>> observed in the second half of the programs, from 8 to 16 (A) and
>>>>>>>>> from 7 to 15 (B) UT. For example,
>>>>>>>>> for source 11(OJ287) S2053A has data only at antennas BR and MK
>>>>>>>>> and S2053C has data only at antenna MK. However, there are no
>>>>>>>> indications
>>>>>>>>> of significant troubles in the log files, although in the VLBA
>>>>>>>> archive
>>>>>>>>> the quality of the data is marked as ERR. On the other hand,
>>>>>>>>> the data for the program S2053B (2009 October 20)
>>>>>>>>> have similar quality but processed more or less normally.
>>>>>>>>> I tried to troubleshoot by doing the following:
>>>>>>>>> 1) I loaded the S2053A epoch twice from the VLBA archive but this
>>>>>>>>> did not change the result.
>>>>>>>>> 2) I processed data on two different computers, PC with LINUX
>>>>>>>> (UBUNTU) with
>>>>>>>>> AIPS version 31DEC2008 and a Mac running OS X 10.5.8 with AIPS
>>>>>>>> version
>>>>>>>>> 31DEC2009. The result is the the same.
>>>>>>>>> 3). For S2053A I performed amplitude calibration by two different
>>>>>>>>> ways: i) - using the TY,GC, and WX tables attached to the file and
>>>>>>>> using
>>>>>>>>> utilities VLBAMCAL and VLBACALA; ii) - deleting these tables and
>>>>>>>>> recreating them in the old manner using tasks VLOG, ANTAB, and
>>>>>>>>> APCAL
>>>>>>>>> and calibration files provided by NRAO. The result is again the
>>>>>>>>> same.
>>>>>>>>> In the case of the first method, although the VLBACALA procedure
>>>>>>>> appears to
>>>>>>>>> have ended successfully, the results of the opacity correction do
>>>>>>>> not look
>>>>>>>>> right:
>>>>>>>>> ===================================================================
>>>>>>>>>
>>>>>>>>>> APCAL1: Task APCAL (release of 31DEC09) begins
>>>>>>>>> localh> APCAL1: TXTWX: Weather data read from WX table 1
>>>>>>>>> localh> APCAL1: DFTAU: INITIAL TAU0 SET TO 0.100
>>>>>>>>> localh> APCAL1: DFTRC:BR RCP Initial Trec set to = 61.40
>>>>>>>>> localh> APCAL1: DFTRC:BR LCP Initial Trec set to = 85.45
>>>>>>>>> localh> APCAL1: DFTRC:FD RCP Initial Trec set to = 60.99
>>>>>>>>> localh> APCAL1: DFTRC:FD LCP Initial Trec set to = 65.90
>>>>>>>>> localh> APCAL1: DFTRC:HN RCP Initial Trec set to = 106.93
>>>>>>>>> localh> APCAL1: DFTRC:HN LCP Initial Trec set to = 98.66
>>>>>>>>> localh> APCAL1: DFTRC:KP RCP Initial Trec set to = 72.36
>>>>>>>>> localh> APCAL1: DFTRC:KP LCP Initial Trec set to = 67.82
>>>>>>>>> localh> APCAL1: DFTRC:LA RCP Initial Trec set to = 65.64
>>>>>>>>> localh> APCAL1: DFTRC:LA LCP Initial Trec set to = 89.13
>>>>>>>>> localh> APCAL1: DFTRC:MK RCP Initial Trec set to = 42.98
>>>>>>>>> localh> APCAL1: DFTRC:MK LCP Initial Trec set to = 56.20
>>>>>>>>> localh> APCAL1: DFTRC:NL RCP Initial Trec set to = 114.53
>>>>>>>>> localh> APCAL1: DFTRC:NL LCP Initial Trec set to = 108.27
>>>>>>>>> localh> APCAL1: DFTRC:OV RCP Initial Trec set to = 70.44
>>>>>>>>> localh> APCAL1: DFTRC:OV LCP Initial Trec set to = 74.91
>>>>>>>>> localh> APCAL1: DFTRC:PT RCP Initial Trec set to = 63.80
>>>>>>>>> localh> APCAL1: DFTRC:PT LCP Initial Trec set to = 63.05
>>>>>>>>> localh> APCAL1: DFTRC:SC RCP Initial Trec set to = 104.14
>>>>>>>>> localh> APCAL1: DFTRC:SC LCP Initial Trec set to = 92.09
>>>>>>>>> localh> APCAL1: Writing SN table 2
>>>>>>>>> localh> APCAL1: BR RCP 1/ 7h 0m Trec (K): 95.51 Zen.
>>>>>>>> opac.: 0.000
>>>>>>>>> localh> APCAL1: BR LCP 1/ 7h 0m Trec (K): 120.84 Zen.
>>>>>>>> opac.: 0.000
>>>>>>>>> localh> APCAL1: Plot file version 11 created
>>>>>>>>> localh> APCAL1: GFINIS: number records used 82
>>>>>>>>> localh> APCAL1: Plot file version 12 created
>>>>>>>>> localh> APCAL1: GFINIS: number records used 59
>>>>>>>>> localh> APCAL1: Plot file version 13 created
>>>>>>>>> localh> APCAL1: GFINIS: number records used 90
>>>>>>>>> localh> APCAL1: Plot file version 14 created
>>>>>>>>> localh> APCAL1: GFINIS: number records used 65
>>>>>>>>> localh> APCAL1: Plot file version 15 created
>>>>>>>>> localh> APCAL1: GFINIS: number records used 42
>>>>>>>>> localh> APCAL1: Appears to have ended successfully
>>>>>>>>> localh> APCAL1: localhos 31DEC09 TST: Cpu= 1.9 Real= 2
>>>>>>>> IO=
>>>>>>>>> 24
>>>>>>>>> localh> CLCAL1: Task CLCAL (release of 31DEC09) begins
>>>>>>>>> localh> CLCAL1: Using interpolation mode SELF
>>>>>>>>> localh> CLCAL1: Processing SN table 2
>>>>>>>>> localh> CLCAL1: SNMRG: Merging SN table
>>>>>>>>> localh> CLCAL1: SNMRG: Write 1476 merged records from 1476
>>>>>>>> input
>>>>>>>>> records
>>>>>>>>> localh> CLCAL1: SN2CL: Applying SN tables to CL table 2, writing
>>>>>>>>> CL
>>>>>>>>> table 3
>>>>>>>>> localh> CLCAL1: Appears to have ended successfully
>>>>>>>>> localh> CLCAL1: localhos 31DEC09 TST: Cpu= 3.5 Real= 9
>>>>>>>> IO=
>>>>>>>>> 556
>>>>>>>>>
>>>>>>>> ========================================================================
>>>>>>>>
>>>>>>>>> In the second way, APCAL ends with an ERROR that indicates a
>>>>>>>>> problem
>>>>>>>>> with the TY table. However, I could not find anything unusual
>>>>>>>>> in the
>>>>>>>>> TY-table, which is attached.
>>>>>>>>> I appreciate any advice you can give me on the possible source of
>>>>>>>>> the
>>>>>>>>> problem and, if needed, a work-around.
>>>>>>>>> Thank you very much,
>>>>>>>>> Svetlana Jorstad
>>>>>>>>>
>>>>>
>>>
>>>
>
More information about the Daip
mailing list