[daip] AIPS bug POSSM_JAN05
Leonid Petrov
Leonid.Petrov at lpetrov.net
Thu Jan 6 14:12:43 EST 2005
Dear AIPS-maintainer,
BUG REPORT
==========
Synopsis: possm works incorrectly in an attempt to plot spectrum of continuum
VLBI data
Why bug: 1) The spectrum looks very weird. Apparently indexing is messed up.
2) AIPS, version name "31DEC04", version date: 20040602 works
correctly;
AIPS, version name "31DEC04", version date: 20050106
and AIPS, version name "31DEC05", version date: 20050106 work
incorrectly using exactly the same input.
AIPS version: version name "31DEC04", version date: 20050106 and
version name "31DEC05", version date: 20050106
How to reproduce:
-----------------
1) Download the fits file
ftp://cddis1.gsfc.nasa.gov/vlbi/raw/fits/1995_07_15_bb023c_fits.tar.bz2
2) untar and use FITLD for reading it into AIPS. The list of input parameters
can be found in the attachment 1.
3) Use possm. The list of input parameters can be found in the attachment 2.
4) Compare that what you see with
a) http://gemini.gsfc.nasa.gov/pet/discussion/bugs/possm_20050101.gif
(erroneous plot, generated by AIPS name 31DEC04, AIPS date 20050105)
b) http://gemini.gsfc.nasa.gov/pet/discussion/bugs/possm_20040602.gif
(Correct plot, generated by AIPS name 31DEC04, AIPS date 20040602)
For the case I've put my .AIPSRC file in the attachment 3.
Sincerely,
Leonid Petrov
06-JAN-2005 14:05:26
------------------------------------------------------------------------------
Leonid Petrov http://gemini.gsfc.nasa.gov/petrov.html
Leonid.Petrov at lpetrov.net
Tel.: +1(301)-614-6096 Fax: +1(301)-614-6099
Code 926, NVI, Inc./NASA Goddard Space Flight Center, Greenbelt, MD 20771 USA
--------------------------- Attachment 1: --------------------------------------
>inp
AIPS 1: FITLD: Task to store an image or UV data from a FITS tape
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: INTAPE 1 Input tape drive # (0 => 1)
AIPS 1: NFILES 0 # of files to advance on tape
AIPS 1: INFILE '/vlbi/AIPS/incoming/19 Disk file name
AIPS 1: 95_07_15_bb023c_s.fit
AIPS 1: '
AIPS 1: OUTNAME 'VCS1_3S ' File name (name)
AIPS 1: OUTCLASS 'UVDATA' File name (class)
AIPS 1: OUTSEQ 0 File name (seq. #)
AIPS 1: 0 => highest unique number
AIPS 1: => matching (on VLBA)
AIPS 1: -1 => FITS tape value
AIPS 1: OUTDISK 1 Disk drive # (0 => any)
AIPS 1: OPTYPE ' ' Type of data to load,
AIPS 1: ' ' => all types
AIPS 1: 'UV' => UV data
AIPS 1: 'IM' => images
AIPS 1: NCOUNT 0 Number of files to load.
AIPS 1: DOTABLE 1 True (1.0) means load tables
AIPS 1: for images.
AIPS 1: DOUVCOMP -1 >0 => compressed data (FITS)
AIPS 1: DOCONCAT -1 >0 -> if VLBA or uv table,
AIPS 1: append data to existing
AIPS 1: files, or if no appropriate
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
AIPS 1: files exist create a new file
AIPS 1: and append all data to that
AIPS 1: file.
AIPS 1: *****************************
AIPS 1: Following adverbs are useful
AIPS 1: only when reading VLBA
AIPS 1: distribution tapes.
AIPS 1: *****************************
AIPS 1: CLINT 0.5 CL entry interval (min)
AIPS 1: 0 => 1 minute.
AIPS 1: SOURCES *all ' ' Source list to accept.
AIPS 1: QUAL -1 Source qualifier -1=>all
AIPS 1: TIMERANG *all 0 Timerange selected
AIPS 1: BCHAN 1 Lowest spectral channel
AIPS 1: number to select in each IF.
AIPS 1: 0=>1
AIPS 1: ECHAN 0 Highest spectral channel
AIPS 1: number to select in each IF.
AIPS 1: 0=>highest
AIPS 1: BIF 0 Lowest IF number 0=>1
AIPS 1: EIF 0 Highest IF number 0=>all
AIPS 1: DIGICOR 1 Controls the application of
AIPS 1: the VLBA correlator's digital
AIPS 1: correction.
AIPS 1: 0 or 1 => perform corrections
AIPS 1: 2 => perform cross-power
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
AIPS 1: corrections, do total-power
AIPS 1: only if zero padding used.
AIPS 1: -1 => do NOT perform
AIPS 1: corrections.
AIPS 1: SEE HELP for more details
AIPS 1: SELBAND -1 Bandwidth to select (kHz)
AIPS 1: SELFREQ -1 Frequency to select (MHz).
AIPS 1: See HELP.
AIPS 1: FQTOL -1 Frequency tolerance assigned
AIPS 1: to SELFREQ, FITLD will select
AIPS 1: data with freq. of SELFREQ
AIPS 1: +/- FQTOL.
AIPS 1: The unit of FQTOL is kHz.
AIPS 1: <=0 => 10 kHz
AIPS 1: WTTHRESH 0.7 Flagging threshold based on
AIPS 1: weights. See HELP.
AIPS 1: 0 => no flagging
AIPS 1: 1 => will flag ALL data
AIPS 1: OPCODE ' ' Specify if wish to keep VT
AIPS 1: (VLBA Tape Statistics) table.
AIPS 1: ' ' => do not copy VT table
AIPS 1: 'VT' => keep VT table
AIPS 1: NPIECE 0 Maximum uv table piece to
AIPS 1: load (ignored for tape unless
AIPS 1: NCOUNT = 1)
--------------------------- Attachment 2: --------------------------------------
AIPS 1: POSSM Task to plot total and cross-power spectra.
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: USERID 0 User number
AIPS 1: INNAME 'VCS1_3S ' Input UV file name (name)
AIPS 1: INCLASS 'UVDATA' Input UV file name (class)
AIPS 1: INSEQ 1 Input UV file name (seq. #)
AIPS 1: INDISK 1 Input UV file disk unit #
AIPS 1: SOURCES *all ' ' Source list
AIPS 1: QUAL -1 Source qualifier -1=>all
AIPS 1: CALCODE ' ' Calibrator code ' '=>all
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: UVRANG 0 0 UV range to be plotted
AIPS 1: TIMERANG 0 23 Time range to be plotted
AIPS 1: 46 9 0 23
AIPS 1: 47 58
AIPS 1: STOKES ' ' Stokes type(s) to select.
AIPS 1: BIF 0 Lowest IF number 0=>all
AIPS 1: EIF 0 Highest IF number 0=>all
AIPS 1: BCHAN 1 Lowest channel number 0=>all
AIPS 1: ECHAN 0 Highest channel number 0=>all
AIPS 1: SUBARRAY 0 Subarray, 0=>1
AIPS 1: ANTENNAS *all 0 Antennas to select
AIPS 1: BASELINE *all 0 Baselines with ANTENNAS
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
AIPS 1: DOCALIB -1 If >0 calibrate data
AIPS 1: = 2 calibrate weights
AIPS 1: GAINUSE 0 CL (or SN) table to apply
AIPS 1: DOPOL -1 If >0 correct polarization.
AIPS 1: BLVER -1 BL table to apply.
AIPS 1: FLAGVER 0 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 -1 Bandpass table version
AIPS 1: SMOOTH *all 0 Smoothing function. See
AIPS 1: HELP SMOOTH for details.
AIPS 1: SHIFT 0 0 Position shift:
AIPS 1: RA, Dec (arcsec)
AIPS 1: 0 => no shift
AIPS 1: APARM 1 1 Control information:
AIPS 1: 0 0 1: = 0 => scalar average
AIPS 1: -180 180 > 0 => vector average
AIPS 1: 0 0 2: = 0 => self-scale
AIPS 1: 1 0 > 0 => fixed scale
AIPS 1: (use APARM(3-6))
AIPS 1: 3: min. amplitude
AIPS 1: 4: max. amplitude
AIPS 1: 5: min. phase (degrees)
AIPS 1: 6: max. phase
AIPS 1: 7: x-axis labelling
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
AIPS 1: = 0 = > in channels.
AIPS 1: = 1 = > in Hz (or secs
AIPS 1: if corr. fn)
AIPS 1: = 2 = > in m/s
AIPS 1: 8: = 0 => plot cross power
AIPS 1: = 1 => plot total power
AIPS 1: = 2 => plot BP table
AIPS 1: = 3 => plot ACF
AIPS 1: = 4 => plot XCF
AIPS 1: = 5 => plot cross power of
AIPS 1: individual antennas
AIPS 1: for cross polaris.
AIPS 1: (STOKES=RL, or LR)
AIPS 1: 9: > 0 => plot several IF's
AIPS 1: and/or polarizations
AIPS 1: together as though one
AIPS 1: long spectrum (see HELP)
AIPS 1: 10: > 0 => reverse direction
AIPS 1: of plotted spectrum, so
AIPS 1: velocity increases to
AIPS 1: right.
AIPS 1: CODETYPE ' ' 'A&P ', 'AMP ', 'PHAS',
AIPS 1: 'R&I ', 'REAL', 'IMAG'
AIPS 1: other => 'A&P '
AIPS 1: POLPLOT ' ' Option to display various
AIPS 1: combinations of polzns to
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
AIPS 1: plot: 'RL/RR', 'RL/LL',
AIPS 1: 'LR/RR', 'LR/LL', 'RR/LL'
AIPS 1: 'LL/RR'; other = don't use
AIPS 1: this option.
AIPS 1: SOLINT 0 If SOLINT > 0 then it enables
AIPS 1: the user to make multiple
AIPS 1: plots per pass of POSSM.
AIPS 1: It defines the averaging time
AIPS 1: for each individual plot.
AIPS 1: Task will start at TIMERANG
AIPS 1: and make a plot for every
AIPS 1: SOLINT minutes. If SOLINT
AIPS 1: = -1 will do the same but
AIPS 1: will do scan averages if NX
AIPS 1: table is present.
AIPS 1: NPLOTS 9 Number of plots per page
AIPS 1: BPARM *all 0 More control information:
AIPS 1: 1: If = 1 divide by 'channel
AIPS 1: 0' before plotting data.
AIPS 1: 0 => do not divide.
AIPS 1: 2: Start chn. of 'channel 0'
AIPS 1: (0 => determined by POSSM)
AIPS 1: 3: Stop chn. of 'channel 0'
AIPS 1: (0 => determined by POSSM)
AIPS 1: 4: ignore spectrum when ampl.
AIPS 1: channel 0 < BPARM(4) Jy
AIPS 1: ** press RETURN for more, enter Q or next line to quit print **
AIPS 1: 5-9: unused
AIPS 1: 10: =1 => don't write header
AIPS 1: info when writing to outfile
AIPS 1: useful for appending several
AIPS 1: spectra into a single outfile
AIPS 1: [see EXPLAIN POSSM]
AIPS 1: OUTFILE ' '
AIPS 1: Filename in which to write
AIPS 1: spectrum. Default = ' ' =
AIPS 1: do not write spectrum. The
AIPS 1: file is written only if
AIPS 1: NPLOTS = 0
AIPS 1: LTYPE 3 Type of labeling: 1 border,
AIPS 1: 2 no ticks, 3 - 6 standard,
AIPS 1: 7 - 10 only tick labels
AIPS 1: <0 -> no date/time
AIPS 1: BADDISK *all 0 Disks to avoid for scratch
AIPS 1: DOTV 1 > 0 Do plot on the TV, else
AIPS 1: make a plot file
AIPS 1: GRCHAN 0 Graphics channel 0 => 1.
-------------------------- Attachemt 3: ---------------------------------------
# AIPSRC created by install.pl on Thu Jan 6 10:49:01 2005
#
LAST_VERSION = 31DEC04
AIPS_ROOT = /opt/aips_20050101
GROUP = users
GROUPWRITE = YES
INSTALL_FROM_CD =
INSTALL_FROM_NET =
ARCH = LINUX
SITE = GSFC
THISHOST = LACERTA
HOSTS = "LACERTA"
PRINTERS = ""
DATA_AREAS = "/vlbi/AIPS/LACERTA_1: /p3/data/AIPS/LACERTA_1:"
TAPE_DRIVES = ""
TPHOSTS = "127.0.0.1"
APSIZE = 5
PRINTMED = A
FORT = /usr/local/bin/g77
FORTOPT = -c -fno-automatic -fno-globals -Wno-globals -malign-double -Wimplicit -Wall
CCOM = /usr/local/bin/gcc
CCOMOPT = -c -O3 -fomit-frame-pointer -funroll-loops -I/opt/aips_20050101/31DEC04/INC -D_FILE_OFFSET_BITS=64 -DHAVE_LINUX_GLIBC
LINK = /usr/local/bin/g77
READLINE = /usr/lib/libreadline.a
DOTWOLIB = NO
DOWNLOADED = 2005.01.06:15:48.01
UNPACKED = Unpacked in /opt/aips_20050101 at Thu Jan 6 10:48:13 2005
#---- End of AIPSRC file
More information about the Daip
mailing list