[daip] complex clean

Krisztina Eva Gabanyi gabanyik at mpifr-bonn.mpg.de
Thu Apr 28 11:54:25 EDT 2005


Hi Uwe, (Dear AIPS friend,)

i had problmes with polarization data reduction and it seems that i have to 
bother you, the specialist with questions :-)

I tried to do this complex clean thing, and i faced with some problems. 
I consulted with Simone and she told me, that it seems that CXPOLN does not 
run properly till the end. It does not destroy some files wich should be 
destroyed and does not create some  which should be created.

After running cxpoln, i got the following output files:
CXPOLN TEMP .QBM001.    1 
CXPOLN TEMP .QIM001.    1 
CXPOLN TEMP .UBM001.    1 
CXPOLN TEMP .UIM001.    1

This is, what aips writes in the message window:
pc186.> UVPOL1: Task UVPOL  (release of 31DEC05) begins
pc186.> UVPOL1: Create CXPOLN TEMP .IMAGE . 666 (UV)  on disk  1  cno  443
pc186.> UVPOL1: Copied FQ file from vol/cno/vers  1  436   1 to  1  443   1
pc186.> UVPOL1: Copied AN file from vol/cno/vers  1  436   1 to  1  443   1
pc186.> UVPOL1: Appears to have ended successfully
pc186.> UVPOL1: pc186        31DEC05 TST: Cpu=       0.2  Real=       0
pc186.> MX   1: Task MX     (release of 31DEC05) begins
pc186.> MX   1: MXIN  : Sort ='TB' not 'X*', Using Un-sorted gridding
pc186.> MX   1: GRDFLT: X and Y convolution type = SPHEROIDAL
pc186.> MX   1: GRDFLT: X and Y parms =  3.0000  1.0000
pc186.> MX   1: GRDFLT: convolution function sampled every 1/100 of a cell
pc186.> MX   1: Create CXPOLN TEMP .QBM001.   1 (MA)  on disk  1  cno  554
pc186.> MX   1: Create CXPOLN TEMP .QIM001.   1 (MA)  on disk  1  cno  555
pc186.> MX   1: Create CXPOLN MX   .TEMP  . 666 (UV)  on disk  1  cno  556
pc186.> MX   1: Copied AN file from vol/cno/vers  1  443   1 to  1  556   1
pc186.> MX   1: Using all     41426 Visibilities
pc186.> MX   1: Got data     Time= 17:25:00 CPU time=      0.30
pc186.> MX   1: UVTBUN: Weighting grid =  1024 X  1024, Box =   0
pc186.> MX   1: UVGRTB: All    89 Rows in  5120k AP (  519 Rows Max)
pc186.> MX   1: UVGRTB: Max U Baseline   161787712. lambda (=   82 cells)
pc186.> MX   1: UVGRTB: Ave    2 Channels; 4.800000E+09 to 4.823938E+09 Hz
pc186.> MX   1: Field   1 Sum of gridding weights =  9.15950E+09
pc186.> MX   1: Fit Gaussian FWHM =  947.703 x  751.850 Microarcsec, PA=  
-40.5
pc186.> MX   1: Beam made    Time= 17:25:01 CPU time=      0.97
pc186.> MX   1: Field    1 min =  -11.39 MilliJy, max =     22.74 MilliJy
pc186.> MX   1: Maps made    Time= 17:25:02 CPU time=      1.57
pc186.> MX   1: MXDISP: Display range =   -11.389    22.744 MilliJy
pc186.> MX   1: ***** Channel    1 done ****: 17:25:02 CPU time=      1.75 
*****
pc186.> MX   1: Field    1 min =  -11.39 MilliJy, max =     22.74 MilliJy
pc186.> MX   1: Appears to have ended successfully
pc186.> MX   1: pc186        31DEC05 TST: Cpu=       1.9  Real=       2
pc186.> MX   1: Task MX     (release of 31DEC05) begins
pc186.> MX   1: MXIN  : Sort ='TB' not 'X*', Using Un-sorted gridding
pc186.> MX   1: GRDFLT: X and Y convolution type = SPHEROIDAL
pc186.> MX   1: GRDFLT: X and Y parms =  3.0000  1.0000
pc186.> MX   1: GRDFLT: convolution function sampled every 1/100 of a cell
pc186.> MX   1: Create CXPOLN TEMP .UBM001.   1 (MA)  on disk  1  cno  557
pc186.> MX   1: Create CXPOLN TEMP .UIM001.   1 (MA)  on disk  1  cno  558
pc186.> MX   1: Create CXPOLN MX   .TEMP  . 667 (UV)  on disk  1  cno  559
pc186.> MX   1: Copied AN file from vol/cno/vers  1  443   1 to  1  559   1
pc186.> MX   1: Using all     41426 Visibilities
pc186.> MX   1: Got data     Time= 17:25:03 CPU time=      0.30
pc186.> MX   1: UVTBUN: Weighting grid =  1024 X  1024, Box =   0
pc186.> MX   1: UVGRTB: All    89 Rows in  5120k AP (  519 Rows Max)
pc186.> MX   1: UVGRTB: Max U Baseline   161787712. lambda (=   82 cells)
pc186.> MX   1: UVGRTB: Ave    2 Channels; 4.800000E+09 to 4.823938E+09 Hz
pc186.> MX   1: Field   1 Sum of gridding weights =  9.15950E+09
pc186.> MX   1: Fit Gaussian FWHM =  947.703 x  751.850 Microarcsec, PA=  
-40.5
pc186.> MX   1: Beam made    Time= 17:25:03 CPU time=      0.96
pc186.> MX   1: Field    1 min =   -8.25 MilliJy, max =     14.24 MilliJy
pc186.> MX   1: Maps made    Time= 17:25:06 CPU time=      1.61
pc186.> MX   1: MXDISP: Display range =    -8.252    14.237 MilliJy
pc186.> MX   1: ***** Channel    1 done ****: 17:25:06 CPU time=      1.81 
*****
pc186.> MX   1: Field    1 min =   -8.25 MilliJy, max =     14.24 MilliJy
pc186.> MX   1: Appears to have ended successfully
pc186.> MX   1: pc186        31DEC05 TST: Cpu=       1.9  Real=       4

And this is which is in the "normal" aips terminal:
AIPS 1: Resumes
AIPS 1: Resumes
AIPS 1: Resumes
AIPS 1: Destroyed  1 extension files of type AN
AIPS 1: Destroyed  1 extension files of type FQ
AIPS 1: Destroyed  1 extension files of type HI
AIPS 1: Destroyed UV image file: catno=    443 disk= 1
AIPS 1: Destroyed  1 extension files of type HI
AIPS 1: Destroyed  1 extension files of type FQ
AIPS 1: Destroyed  1 extension files of type AN
AIPS 1: Destroyed UV image file: catno=    556 disk= 1
AIPS 1: Destroyed  1 extension files of type HI
AIPS 1: Destroyed  1 extension files of type FQ
AIPS 1: Destroyed  1 extension files of type AN
AIPS 1: Destroyed UV image file: catno=    559 disk= 1
AIPS 1: I WON'T DESTROY ON DEFAULT NAME INNAME = '            '
AIPS 1: NO DESTROY

Do you have any idea why it does not delete temporary(looking) files? 
When i used these files in CXCLN, i managed to create QCL and UCL files, and 
with COMB create the polarized intensity and polarization angle maps. They 
looked "nice", but compared to the polarization map obtained via "normal 
clean" (IMAGR)  (on the same datasets), the complex method gave only 70% of 
the polarization.

Maybe the midnight job changed something?

Cheers,

Krisztina Gabanyi

-- 
Krisztina Eva Gabanyi
Max-Planck-Institut fuer Radioastronomie
Auf dem Huegel 69
D-53121 Bonn
Germany

Tel.: (+49)228 525 473
e-mail: gabanyik at mpifr-bonn.mpg.de
Homepage: http://www.mpifr-bonn.mpg.de/staff/kgabanyi




More information about the Daip mailing list