[daip] Re: aips task pbeam
Leonia Kogan
lkogan at nrao.edu
Wed Jun 8 10:59:19 EDT 2005
Juan,
You'll find the two file at the attachment:
1. colswitch.f - the fortran codes to switch the first two collumns of
the hol. file
Just for the case if you do not know how to compile the file:
1. type f77 colswitch.f to compile
2. type a.out to run
3. do not forget to change the input and output file names in
colswitch.f
2. COL_SWITCH the output file of the colswitch.f
I used the input file as the shorted (to the first antenna only)
your file (as I remember
h1LL01)
I ran PBEAM with INFILE = COL_SWITCH IN2FILE =''
It works as I expected.
TRy it and let me know the result please
Leonia
Juan Uson wrote:
>
> Leonia,
>
> On further thought, it seems to me that the order of the rows is
> determined by the direction that the array was scanning. However, I
> still think that I could do a reordering in linux if that is needed.
>
> I could try to pre-average the data but I don't know off-hand how to
> do this in UVHOL. I'll see if I can figure it out.
>
> Meanwhile, I'll run UVHOL one antenna at a time. I have found that I
> can do this using all 4 reference antennas and PBEAM will work on the
> output of UVHOL and actually give me an answer that seems reasonable.
> So maybe this is the best way for me to proceed right now.
>
> Once again, thank you for taking the time to look into this and for
> your comments.
>
> Juan
>
>
> On Jun 7, 2005, at 7:54 PM, Leonia Kogan wrote:
>
>> Juan,
>>
>> I have looked PBEAM with debugger and found:
>>
>> 1.PBEAM aborts with your original infile because there is too much
>> blocks (segments)
>> with different ref antenna (100). the relevant array has the
>> dimension 50.
>> When I increased it to 200, PBEAM stop aborting but spend too much
>> time to finish.
>>
>> 2. I made your INFILE shorter (1 segment for each ref antenna
>> 4,14,7,25. So one antenna=1). It works and finishes, but still I
>> think you'll have a problem with
>> your INFILE structure.
>>
>> PBEAM pay a lot of attention in recognizing identical X for the given
>> Y and vise versa.
>> PBEAM consider that Y column (the second column) changes faster than
>> the first one.
>> So there is many Y's for the given X and then another set of Y's for
>> another X..
>>
>> Your INFILE has the opposite structure. As a result there are just a
>> few point at each plots.
>>
>> So pay attention on this creating the file by UVHOL if it is possible.
>>
>> Leonia
>>
>> Juan Uson wrote:
>>
>>> Hello again.
>>>
>>> Your note on the unreasonable values causing it to blow up was a
>>> great help! I also had a brief chat with Rick which helped me
>>> understand the options. I have now run UVHOL only for one antenna
>>> (with one reference and with the 4 reference antennas) and used the
>>> output of that to run PBEAM and it works in this mode.
>>>
>>> So it seems that I might have one or more of the antennas giving me
>>> some bad values. I can run the procedure one antenna at a time to
>>> see which one makes it fail...
>>>
>>> In any case, it seems that I can make some progress this way without
>>> you wasting time. Thank you very much for your help,
>>>
>>> Juan
>>>
>>>
>>> On Jun 7, 2005, at 4:43 PM, Leonia Kogan wrote:
>>>
>>>> Juan,
>>>>
>>>> I have tried PBEAM with my INFILE and it works.
>>>> Using your file it aborts with division on zero or sqrt of negative
>>>> numbers.
>>>>
>>>> So I need to look the PBEAM with debugger using your INFILE
>>>>
>>>> Leonia
>>>>
>>>>
>>>>
>>>> Juan Uson wrote:
>>>>
>>>>>
>>>>> Hi Leonia,
>>>>>
>>>>> Thank you for your quick reply.
>>>>>
>>>>> As far as I can see, the format of the files seems OK. I have
>>>>> copied them to my account at the AOC. They are in the directory
>>>>>
>>>>> /users/juson/holo
>>>>>
>>>>> which I have accessed via slogin to zaurak. I trust that you can
>>>>> access them from your account.
>>>>>
>>>>> Cheers,
>>>>>
>>>>> Juan
>>>>>
>>>>>
>>>>> On Jun 7, 2005, at 1:03 PM, Leonia Kogan wrote:
>>>>>
>>>>>> Juan,
>>>>>>
>>>>>> Look please the INFILE (IN2FILE) which you use as input for PBEAM.
>>>>>> Compare the format of the file(s) with those described at the
>>>>>> PBEAM help.
>>>>>> Send the file to me.
>>>>>>
>>>>>> Leonia
>>>>>>
>>>>>>
>>>>>> Juan Uson wrote:
>>>>>>
>>>>>>>
>>>>>>> Hi Rick and Leonia,
>>>>>>>
>>>>>>> I am trying to analyze holography data. I have calibrated these
>>>>>>> and next I ran 'UVHOL'
>>>>>>> which produces output that seems reasonable.
>>>>>>>
>>>>>>> However, 'PBEAM' crashes. I have tried this on several
>>>>>>> different aips installations using
>>>>>>> tst, new and old. For example, old (31DEC03) complains that:
>>>>>>>
>>>>>>> >go
>>>>>>> PBEAM1: Task PBEAM (release of 31DEC03) begins
>>>>>>> PBEAM1: ZTXOP2: using translated file name =
>>>>>>> PBEAM1: ZTXOP2: /users/juson/holo/hlLL01
>>>>>>> AIPS 1: Resumes
>>>>>>> >PBEAM1: ZTXOP2: using translated file name =
>>>>>>> PBEAM1: ZTXOP2: /users/juson/holo/hlRR01
>>>>>>> PBEAM1: ZTXOP2: using translated file name =
>>>>>>> PBEAM1: ZTXOP2: /users/juson/holo/hlLL01
>>>>>>> PBEAM1: ZTXOP2: using translated file name =
>>>>>>> PBEAM1: ZTXOP2: /users/juson/holo/hlRR01
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> ZMSGOP: FILE DA01:MSD000000.000; NOT FOUND
>>>>>>> ZTTOP2: error translating /dev/pts/0
>>>>>>> ERROR 2 ASSIGNING OUTPUT TO /dev/pts/0
>>>>>>> ZMSGOP: ZTTOPN RETURNS ERROR = 1 FOR LUN = 6
>>>>>>> I/O recursion: I/O started while already doing I/O
>>>>>>> apparent state: internal I/O
>>>>>>> last format: (Reference antennas ,30I3)
>>>>>>> lately writing sequential formatted internal IO
>>>>>>>
>>>>>>> and it crashes at this point.
>>>>>>>
>>>>>>> Is there some trick that I am missing, some initialization that
>>>>>>> is needed?
>>>>>>>
>>>>>>> Cheers,
>>>>>>>
>>>>>>> Juan
>>>>>>
>>>>>>
>>>>>>
>>>>>>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: COL_SWITCH
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20050608/b2f08592/attachment.ksh>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: colswitch.f
Type: text/x-fortran
Size: 893 bytes
Desc: not available
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20050608/b2f08592/attachment.bin>
More information about the Daip
mailing list