[daip] PCAL claims no model provided - wrong.
R. Craig Walker
cwalker at nrao.edu
Thu Apr 14 17:25:57 EDT 2016
After almost exactly 3 years from when I tried using PCAL to get frequency
dependent D-terms and had to report problems, I am finally trying it
again. At that time, other distractions prevented me from testing the fix
that was provided. This time I am having a problem that is probably
simpler but prevents the program from finishing. This problem has
appeared since my use in 2013.
The data are BW098A, a VLBA 43 GHz observation of M87 from Jan. 2013.
Looking back at email history, I reported problems to daip on April 12,
2013. The problem was that, after split, the data were all zero. A FITS
file from April 16, 2013 is at
/home/noatak2/datafiles/M87_transient/BW098A/BW098A_43.MULTI.1
That will likely suffer from the problem reported here.
At the time of the original discussion, I provided that file so you may
still have a copy, although at 13 GB it was big enough to likely get
purged.
The current issue is that PCAL dies claiming
"POLARIZATION MODEL NOT PROVIDED"
This is after giving details of the apparently successful loads of the I,
Q, U, and V images of the calibrator.
I have tried using exactly the same inputs as used in 2013 on the same
data set. At that time, the program ran. It just did not give good
calibrated data downstream after split (all zero).
I can't see why it doesn't think I have a model.
I am running on a MAC under OSX 10.11 with an aips updated a few minutes ago.
Here are the messages as PCAL ran:
>go pcal
PCAL 1: Task PCAL (release of 31DEC16) begins
AIPS 1: Resumes
>PCAL 1: UVGET: Using flag table version 4 to edit data
PCAL 1: TOO LITTLE DATA FOR SPECTRAL FIT OR SINGULAR
PCAL 1: Processing IF number 1 channel 1
PCAL 1: PCLCOP: Copying data for IF 1 to work file
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 1
PCAL 1: Processing IF number 1 channel 2
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 2
PCAL 1: Processing IF number 1 channel 3
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 3
PCAL 1: Processing IF number 1 channel 4
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 4
PCAL 1: Processing IF number 1 channel 5
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 5
PCAL 1: Processing IF number 1 channel 6
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 6
PCAL 1: Processing IF number 1 channel 7
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 7
PCAL 1: Processing IF number 1 channel 8
PCAL 1: PCLSEL: NO DATA SELECTED FOR IF 1 CHAN 8
PCAL 1: Processing IF number 1 channel 9
PCAL 1: Using Clean Component source model
PCAL 1: NO FLUX FOUND FOR CALIBRATOR SOURCE IN SU TABLE
PCAL 1: MODEL WILL USE THE FLUX CONTAINED IN IT ONLY
PCAL 1: FACSET: 3.145688 Jy found from 3000 components
PCAL 1: QINIT: did a GET of 5120 Kwords, OFF 4832271
PCAL 1: ALGMEM field 1 Ipol gridded model subtraction chns 1- 1
PCAL 1: ALGMEM: at visibility number 1
PCAL 1: ALGMEM: 512 - 0 cells, with 16740 Pts
PCAL 1: Using Clean Component source model
PCAL 1: NO FLUX FOUND FOR CALIBRATOR SOURCE IN SU TABLE
PCAL 1: MODEL WILL USE THE FLUX CONTAINED IN IT ONLY
PCAL 1: FACSET: -0.023715 Jy found from 3000 components
PCAL 1: ALGMEM field 1 Qpol gridded model subtraction chns 1- 1
PCAL 1: ALGMEM: at visibility number 1
PCAL 1: ALGMEM: 512 - 0 cells, with 16740 Pts
PCAL 1: Using Clean Component source model
PCAL 1: NO FLUX FOUND FOR CALIBRATOR SOURCE IN SU TABLE
PCAL 1: MODEL WILL USE THE FLUX CONTAINED IN IT ONLY
PCAL 1: FACSET: -0.008601 Jy found from 3000 components
PCAL 1: ALGMEM field 1 Upol gridded model subtraction chns 1- 1
PCAL 1: ALGMEM: at visibility number 1
PCAL 1: ALGMEM: 512 - 0 cells, with 16740 Pts
PCAL 1: Using Clean Component source model
PCAL 1: NO FLUX FOUND FOR CALIBRATOR SOURCE IN SU TABLE
PCAL 1: MODEL WILL USE THE FLUX CONTAINED IN IT ONLY
PCAL 1: FACSET: 0.001776 Jy found from 3000 components
PCAL 1: ALGMEM field 1 Vpol gridded model subtraction chns 1- 1
PCAL 1: ALGMEM: at visibility number 1
PCAL 1: ALGMEM: 512 - 0 cells, with 16740 Pts
PCAL 1: POLARIZATION MODEL NOT PROVIDED
PCAL 1: Purports to die of UNNATURAL causes
PCAL 1: localhos 31DEC16 TST: Cpu= 2.5 Real= 17 IO= 1106
I presume that the initial no data complaints are because of ICHANSEL,
which specifies starting with channel 9 (this is a bit of a harsh way
to deal with a normal input). I tried with ichansel=0. The initial
no data complaints were gone, but what followed was the same.
I presume the no flux in the SU table complaints are benign given that
the images are available.
My notes from when I tried this is in 2013 indicate it was running slowly
because it was doing all channels. This time, it quit.
The inputs are:
>inp pcal
AIPS 1: PCAL Task to compute polarization corrections
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INNAME 'BW098A_43' Input UV file name (name)
AIPS 1: INCLASS 'MULTI' Input UV file name (class)
AIPS 1: INSEQ 1 Input UV file name (seq. #)
AIPS 1: INDISK 5 Input UV file disk unit #
AIPS 1: Data selection (multisource):
AIPS 1: CALSOUR 'OJ287' Sources to calibrate with
AIPS 1: *rest ' '
AIPS 1: QUAL -1 Calibrator qualifier -1=>all
AIPS 1: CALCODE ' ' Calibrator code ' '=>all
AIPS 1: TIMERANG *all 0 Time range to use.
AIPS 1: SELBAND -1 Bandwidth to select (kHz)
AIPS 1: SELFREQ -1 Frequency to select (MHz)
AIPS 1: FREQID -1 Freq. ID to select.
AIPS 1: BIF 0 Lowest IF number 0=>all
AIPS 1: EIF 0 Highest IF number 0=>all
AIPS 1: ANTENNAS *all 0 Antennas to solve for.
AIPS 1: UVRANGE 0 0 UV range in kilolamdba
AIPS 1: SUBARRAY 0 Subarray, 0=>all
AIPS 1: Cal. info for input:
AIPS 1: DOCALIB 2 > 0 calibrate data & weights
AIPS 1: > 99 do NOT calibrate weights
AIPS 1: GAINUSE 19 CAL table to apply.
AIPS 1: CLEAN map (optional)
AIPS 1: BLVER -1 BL table to apply.
AIPS 1: FLAGVER 4 Flag table version
AIPS 1: DOBAND 1 If >0 apply bandpass cal.
AIPS 1: Method used depends on value
AIPS 1: of DOBAND (see HELP file).
AIPS 1: BPVER 3 Bandpass table version
AIPS 1: SMOOTH 6 8 Smoothing function. See
AIPS 1: 16 HELP SMOOTH for details.
AIPS 1: ICHANSEL 9 56 Array of start and stop chan
AIPS 1: 1 *rest 0 numbers, plus a channel
AIPS 1: increment and IF to be used
AIPS 1: to select channels to sum to
AIPS 1: find the polarization.
AIPS 1: 0 => all channels
AIPS 1: IN2NAME 'OJ287_98A_43' Cleaned map name (name)
AIPS 1: IN2CLASS 'ICL001' Cleaned map name (class)
AIPS 1: IN2SEQ 71 Cleaned map name (seq. #)
AIPS 1: IN2DISK 5 Cleaned map disk unit #
AIPS 1: INVERS 1 CC file version #.
AIPS 1: NCOMP *all 0 # comps to use for model.
AIPS 1: 1 value per field
AIPS 1: FLUX 0 Lowest CC component used.
AIPS 1: NMAPS 1 No. Clean map files
AIPS 1: CMETHOD ' ' Modeling method:
AIPS 1: 'DFT','GRID',' '
AIPS 1: CMODEL ' ' Model type: 'COMP','IMAG'
AIPS 1: 'SUBI' (see HELP re images)
AIPS 1: DOMODEL 1 > 0 => use model, do not
AIPS 1: fit source Q and U
AIPS 1: = 2 use CP table model
AIPS 1: PMODEL *all 0 Source poln. model
AIPS 1: SPECPARM *all 0 Spectral index: I Q U V for
AIPS 1: each CALSOUR
AIPS 1: DOSCALE 0 >= 0 use spectral index
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
#
#
AIPS 1: = 2 solve for curvature
AIPS 1: SOLINT 1 Soln. interval (min) 0=>5.
AIPS 1: SOLTYPE 'RAPR' Solution type:
AIPS 1: 'ORI-', 'APPR', 'RAPR'
AIPS 1: SPECTRAL 1 > 0 do spectral PCAL
AIPS 1: <= 0 do continuum PCAL
AIPS 1: INTPARM *all 0 Smoothing parameters for
AIPS 1: spectral mode
AIPS 1: PRTLEV 2 Print statistics 0=>none
AIPS 1: 1 = some, 2 = lots. Use 1.
AIPS 1: REFANT 5 Reference antenna, 0->pick
AIPS 1: BPARM *all 0 Task enrichment parameters
AIPS 1: for SOLTYPE 'ORI-' only:
AIPS 1: (1) if > 0 use default
AIPS 1: initial feed parameters.
AIPS 1: (2) > 0 -> no error calc.
AIPS 1: (3) if > 0 then fit for R-L
AIPS 1: phase difference
AIPS 1: (4) initial R-L phase
AIPS 1: (5) >0 solve for Vpol.
AIPS 1: (6) >0 fix ref. ori. 1
AIPS 1: (7) >0 fix ref. ori. 2
AIPS 1: (8) >0 fix all orientations
AIPS 1: (9) >0 fix all ellipticities
AIPS 1: (10) >0 fix source poln.
AIPS 1: CPARM *all 0 Task enrichment parameters
AIPS 1: (1) >0 => average in IF
AIPS 1: and find common solution.
AIPS 1: (2) >0 => update source
AIPS 1: table with pol. solution
AIPS 1: (3) >0 => do NOT interpolate
AIPS 1: over flagged channels
AIPS 1: (4) CP table version if
AIPS 1: DOMODEL = 2
AIPS 1: (5) >0 Use linear model for
AIPS 1: linear polarization data
AIPS 1: (7) >0 => use initial guess
AIPS 1: from AN or PD tables
AIPS 1: (8) >0 => max. no. iter
AIPS 1: (9) >0 => conv. tolerance
AIPS 1: (10) >0 => conv. tol.
AIPS 1: BADDISK *all 0 Disk no. not to use for
AIPS 1: scratch files.
Am I missing something - is something set wrong? If not, there is a
problem.
The OJ287 images load to the TV ok for whatever that is worth (probably no
un-applied format changes).
Cheers,
Craig
------------------------------------------------------------------
R. Craig Walker Scientist Emeritus
1305 Vista Dr. Array Operations Center
Socorro NM 87801 USA National Radio Astronomy Observatory
cwalker at nrao.edu P.O. Box O
Phone 575 835 3972 Socorro, NM 87801 USA
575 835 7247
------------------------------------------------------------------
More information about the Daip
mailing list