[daip] Re: complex clean
Uwe Bach
bach at to.astro.it
Fri Apr 29 04:34:25 EDT 2005
Hi Krisztina,
I'm not sure were your problem comes from, but the files you get are not
complete (okay, that was 10 minutes ago, you can read the possible
solution further down, but I don't want to rewrite everything!). The
output you sent stops after the first half of the original CXPOLN
procedure. CXPOLN is properly finished when this message appears:
"AIPS 1: 'NORMALIZATION FACTOR =' 1.002638"
Here the complete output which should appear in the normal AIPS command
window (in brackets are my comments):
AIPS 1: Resumes (First run of UVPOL in IMAGE mode)
AIPS 1: Resumes (First run of MX)
AIPS 1: Resumes (Second run of MX)
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= 264 disk= 7
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= 267 disk= 7
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= 270 disk= 7
AIPS 1: Destroyed 1 extension files of type HI
AIPS 1: Destroyed MA image file: catno= 265 disk= 7
AIPS 1: Destroyed 1 extension files of type HI
AIPS 1: Destroyed MA image file: catno= 268 disk= 7
AIPS 1: Resumes (Second run of UVPOL in BEAM mode)
AIPS 1: Resumes (First run of MX)
AIPS 1: Resumes (Second run of MX)
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= 264 disk= 7
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= 268 disk= 7
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= 272 disk= 7
AIPS 1: Destroyed 1 extension files of type HI
AIPS 1: Destroyed MA image file: catno= 265 disk= 7
AIPS 1: Destroyed 1 extension files of type HI
AIPS 1: Destroyed MA image file: catno= 270 disk= 7
AIPS 1: Rename TEST_A2 .QIM001. 2 (MA) on disk 7 cno 267
AIPS 1: To TEST_A2 .QBM001. 1 (MA) on disk 7 cno 267
AIPS 1: Rename TEST_A2 .UIM001. 2 (MA) on disk 7 cno 271
AIPS 1: To TEST_A2 .UBM001. 1 (MA) on disk 7 cno 271
AIPS 1: 'NORMALIZATION FACTOR =' 1.002638
I ran CXPOLN for my VSOP data and it works fine. I also had a look at
the CXPOLN version at the MPI and it is the same. Did you perform the
run CXPOLN this year
Because the old version before Jan 2005 was corrupt. You should try to
load it again.
There is also no DOPOL parameter for MX so you need to apply your
d-terms first with SPLIT to your uv-data. Otherwise they will not be
used. Just SPLIT your single source uv-data file gain with SPLIT using
DOPOL = 1 and I would also average the IFs aparm(1)=3.
Another possibility to avoid that some parameters in UVPOL or MX are set
to some strange values is to type
default UVPOL
defualt MX
before running CXPOLN.
Yuppiii!!! Here it is. I just looked at your output again and found that
it aborts with
> AIPS 1: I WON'T DESTROY ON DEFAULT NAME INNAME = ' '
> AIPS 1: NO DESTROY
So the error is that the inname in the ZAP procedure was not set. CXPOLN
does that from the outname with inn=outn. Therefore you have to specify
the outname in CXPOLN. I guess you left this blank. You also need to set
the OUTSEQ number.
Try this and tell me if it works. If not then just send me the data file
and I will try it myself.
Cheers
Uwe
Krisztina Eva Gabanyi wrote:
> 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
>
--
_________________________________________
Uwe Bach
INAF - Osservatorio Astronomico di Torino
Via Osservatorio 20
10025 Pino Torinese (TO), ITALY
http://www.to.astro.it/blazars/bach
phone: +39 011 8101902
fax: +39 011 8101930
More information about the Daip
mailing list