[daip] POSSM dies abruptly if STOKES 'HALF' and CL entries blank for LCP
Laura Spitler
lspitler at aoc.nrao.edu
Wed Jul 23 13:24:42 EDT 2003
Dear DAIP,
consider a VLBA data set, 8 IFs, full polarization, with many blank entries in
the CL table (as is not unusual for such data! ;). In particular, *all* of
Stokes LL is blank in the CL table for FD-HN.
We run POSSM with STOKES 'RR'; FD-HN is (correctly) happily plotted.
With Stokes 'LL', FD-HN is (correctly) skipped over, since there are
no data
BUT, with Stokes 'HALF', FD-HN is nicely plotted for RR; but POSSM
dies (without aborting!) when trying to do LL, with the messages
...
spitle> POSSM2: POSSUV: on baseline(s) FD - HN
spitle> POSSM2: Covering timerange 00/20:46:33 - 00/20:47:28
spitle> POSSM2: Will plot 8 channels * 4 IFs polarization RR
spitle> POSSM2: FILLPL: Amplitude mean: 2.2049E-02 rms: 4.3033E-02
spitle> POSSM2: FILLPL: Phase mean: 7.7973E+00 rms: 5.6373E+01
spitle> POSSM2: Doing: FD - HN IF: 1
spitle> POSSM2: -------------------------
spitle> POSSM2: Waiting 30 seconds:
spitle> POSSM2: Hit TV button A to pause indefinitely. Hit button
spitle> POSSM2: B or C to continue sooner, button D to stop plotting
spitle> POSSM2: Continuing
spitle> POSSM2: Will plot 8 channels * 4 IFs polarization LL
spitle> POSSM2: FILLPL: Amplitude mean: 0.0000E+00 rms: 0.0000E+00
spitle> POSSM2: FILLPL: Phase mean: 0.0000E+00 rms: 0.0000E+00
spitle> POSSM2: TICINC: TIC MARK ALGORITHM FAILED! CONTINUING.
appearing in the MSGSRV, and the lovely error message
I/O recursion: I/O started while already doing I/O
apparent state: unit 86549631
appearing in the main AIPS (typing) window. SPY shows that POSSM no longer
exists, although the file still has READ status noted in UCAT. Note that using
the identical inputs, with either docal=-1, or (docal=2; gainuse 1), gives
perfectly normal (no error) results.
Inputs to POSSM are given below. This is AIPS TST (31DEC03) on spitler
(Linux box). The data are available on spitler, user no. 2442; save/get file for
offending POSSM inputs is badpossm.
Cheers,
Michael & Laura
AIPS 2: POSSM Task to plot total and cross-power spectra.
AIPS 2: Adverbs Values Comments
AIPS 2: ----------------------------------------------------------------
AIPS 2: USERID 0 User number
AIPS 2: INNAME 'BM111A ' Input UV file name (name)
AIPS 2: INCLASS 'FQ-1 ' Input UV file name (class)
AIPS 2: INSEQ 1 Input UV file name (seq. #)
AIPS 2: INDISK 1 Input UV file disk unit #
AIPS 2: SOURCES 'J1907+0127 ' Source list
AIPS 2: *rest ' '
AIPS 2: QUAL -1 Source qualifier -1=>all
AIPS 2: CALCODE ' ' Calibrator code ' '=>all
AIPS 2: SELBAND -1 Bandwidth to select (kHz)
AIPS 2: SELFREQ -1 Frequency to select (MHz)
AIPS 2: FREQID 1 Freq. ID to select.
AIPS 2: UVRANG 0 0 UV range to be plotted
AIPS 2: TIMERANG *all 0 Time range to be plotted
AIPS 2: STOKES 'HALF' Stokes type(s) to select.
AIPS 2: BIF 0 Lowest IF number 0=>all
AIPS 2: EIF 0 Highest IF number 0=>all
AIPS 2: BCHAN 1 Lowest channel number 0=>all
AIPS 2: ECHAN 0 Highest channel number 0=>all
AIPS 2: SUBARRAY 0 Subarray, 0=>1
AIPS 2: ANTENNAS 2 *rest 0 Antennas to select
AIPS 2: BASELINE *all 0 Baselines with ANTENNAS
AIPS 2: DOCALIB 2 If >0 calibrate data
AIPS 2: = 2 calibrate weights
AIPS 2: GAINUSE 6 CL (or SN) table to apply
AIPS 2: DOPOL -1 If >0 correct polarization.
AIPS 2: FLAGVER 1 Flag table version
AIPS 2: DOBAND -1 If >0 apply bandpass cal.
AIPS 2: Method used depends on value
AIPS 2: of DOBAND (see HELP file).
AIPS 2: BPVER -1 Bandpass table version
AIPS 2: SMOOTH *all 0 Smoothing function. See
AIPS 2: HELP SMOOTH for details.
AIPS 2: SHIFT 0 0 Position shift:
AIPS 2: RA, Dec (arcsec)
AIPS 2: 0 => no shift
AIPS 2: APARM 1 1 Control information:
AIPS 2: 0 0 1: = 0 => scalar average
AIPS 2: -180 180 > 0 => vector average
AIPS 2: 0 0 2: = 0 => self-scale
AIPS 2: 1 0 > 0 => fixed scale
AIPS 2: (use APARM(3-6))
AIPS 2: 3: min. amplitude
AIPS 2: 4: max. amplitude
AIPS 2: 5: min. phase (degrees)
AIPS 2: 6: max. phase
AIPS 2: 7: x-axis labelling
AIPS 2: = 0 = > in channels.
AIPS 2: = 1 = > in Hz (or secs
AIPS 2: if corr. fn)
AIPS 2: = 2 = > in m/s
AIPS 2: 8: = 0 => plot cross power
AIPS 2: = 1 => plot total power
AIPS 2: = 2 => plot BP table
AIPS 2: = 3 => plot ACF
AIPS 2: = 4 => plot XCF
AIPS 2: = 5 => plot cross power of
AIPS 2: individual antennas
AIPS 2: for cross polaris.
AIPS 2: (STOKES=RL, or LR)
AIPS 2: 9: > 0 => plot several IF's
AIPS 2: and/or polarizations
AIPS 2: together as though one
AIPS 2: long spectrum (see HELP)
AIPS 2: 10: > 0 => reverse direction
AIPS 2: of plotted spectrum, so
AIPS 2: velocity increases to
AIPS 2: right.
AIPS 2: CODETYPE ' ' 'A&P ', 'AMP ', 'PHAS',
AIPS 2: 'R&I ', 'REAL', 'IMAG'
AIPS 2: other => 'A&P '
AIPS 2: POLPLOT ' ' Option to display various
AIPS 2: combinations of polzns to
AIPS 2: plot: 'RL/RR', 'RL/LL',
AIPS 2: 'LR/RR', 'LR/LL', 'RR/LL'
AIPS 2: 'LL/RR'; other = don't use
AIPS 2: this option.
AIPS 2: SOLINT -1 If SOLINT > 0 then it enables
AIPS 2: the user to make multiple
AIPS 2: plots per pass of POSSM.
AIPS 2: It defines the averaging time
AIPS 2: for each individual plot.
AIPS 2: Task will start at TIMERANG
AIPS 2: and make a plot for every
AIPS 2: SOLINT minutes. If SOLINT
AIPS 2: = -1 will do the same but
AIPS 2: will do scan averages if NX
AIPS 2: table is present.
AIPS 2: NPLOTS 1 Number of plots per page
AIPS 2: BPARM *all 0 More control information:
AIPS 2: 1: If = 1 divide by 'channel
AIPS 2: 0' before plotting data.
AIPS 2: 0 => do not divide.
AIPS 2: 2: Start chn. of 'channel 0'
AIPS 2: (0 => determined by POSSM)
AIPS 2: 3: Stop chn. of 'channel 0'
AIPS 2: (0 => determined by POSSM)
AIPS 2: 4: ignore spectrum when ampl.
AIPS 2: channel 0 < BPARM(4) Jy
AIPS 2: 5-9: unused
AIPS 2: 10: =1 => don't write header
AIPS 2: info when writing to outfile
AIPS 2: useful for appending several
AIPS 2: spectra into a single outfile
AIPS 2: [see EXPLAIN POSSM]
AIPS 2: OUTFILE ' '
AIPS 2: Filename in which to write
AIPS 2: spectrum. Default = ' ' =
AIPS 2: do not write spectrum. The
AIPS 2: file is written only if
AIPS 2: NPLOTS = 0
AIPS 2: LTYPE 0 Type of labeling: 1 border,
AIPS 2: 2 no ticks, 3 - 6 standard,
AIPS 2: 7 - 10 only tick labels
AIPS 2: <0 -> no date/time
AIPS 2: BADDISK *all 0 Disks to avoid for scratch
AIPS 2: DOTV 1 > 0 Do plot on the TV, else
AIPS 2: make a plot file
AIPS 2: GRCHAN 0 Graphics channel 0 => 1.
More information about the Daip
mailing list