[daip] problem with fring when including source structure
Juan-Carlos Algaba-Marcos
algaba at kasi.re.kr
Wed Jan 28 07:14:32 EST 2015
Dear AIPS guru,
I am trying to fringe fitting a uv dataset with a source model for phase reference methods. The source has actually been splited from the very same uvdata and imaged in difmap, The second fringe fitting is done using the CL table prior the initial fringe fitting.
When I run FRING, I get the following error message:
localh> FRING1: Task FRING (release of 31DEC12) begins
localh> FRING1: UVGET: doing no flagging this time
localh> FRING1: Set integration time to 0.819200 seconds, hope that is ok
localh> FRING1: Selecting and calibrating the data
localh> FRING1: SETGDS: IMAGE 1 EPOCH 2000.0 UV EPOCH -1.0
localh> FRING1: USE UVFIX ON THE UV DATA OR EPOSW ON THE IMAGE
localh> FRING1: Purports to die of UNNATURAL causes
Apparently, uvdata and image do not have the same epoch. But if I look at their headers, it seems that the coordinate equinox is the same (J2000):
>imh
AIPS 1: Image=MULTI (UV) Filename=IMOGABA9F043.UVDATA. 1
AIPS 1: Telescope=VLBA Receiver=VLBA
AIPS 1: Observer=N13SL01 User #= 251
AIPS 1: Observ. date=19-NOV-2013 Map date=26-JAN-2015
AIPS 1: # visibilities 249176 Sort order TB
AIPS 1: Rand axes: UU-L-SIN VV-L-SIN WW-L-SIN TIME1 BASELINE
AIPS 1: SOURCE FREQSEL INTTIM GATEID CORR-ID
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: COMPLEX 3 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: STOKES 1 -2.0000000E+00 1.00 -1.0000000E+00 0.00
AIPS 1: FREQ 256 4.3400000E+10 0.52 6.2500000E+04 0.00
AIPS 1: IF 4 1.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: RA 1 00 00 00.000 1.00 0.000000 0.00
AIPS 1: DEC 1 00 00 00.000 1.00 0.000000 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Maximum version number of extension files of type FQ is 1
AIPS 1: Maximum version number of extension files of type NX is 1
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type AT is 1
AIPS 1: Maximum version number of extension files of type CT is 1
AIPS 1: Maximum version number of extension files of type AN is 1
AIPS 1: Maximum version number of extension files of type BP is 1
AIPS 1: Maximum version number of extension files of type CL is 5
AIPS 1: Maximum version number of extension files of type CQ is 1
AIPS 1: Maximum version number of extension files of type GC is 1
AIPS 1: Maximum version number of extension files of type IM is 1
AIPS 1: Maximum version number of extension files of type MC is 1
AIPS 1: Maximum version number of extension files of type SN is 4
AIPS 1: Maximum version number of extension files of type SU is 1
AIPS 1: Maximum version number of extension files of type TY is 1
AIPS 1: Keyword = 'OLDRFQ ' value = 2.17000000D+10
>
>im2h
AIPS 1: Header for second input image or data set
AIPS 1: Image=0827+243 (MA) Filename=0827+243 .ICC043. 1
AIPS 1: Telescope=VLBA Receiver=VLBA
AIPS 1: Observer=N13SL01 User #= 251
AIPS 1: Observ. date=19-NOV-2013 Map date=28-JAN-2015
AIPS 1: Minimum=-1.19117908E-02 Maximum= 7.42246687E-01 JY/BEAM
AIPS 1: ----------------------------------------------------------------
AIPS 1: Type Pixels Coord value at Pixel Coord incr Rotat
AIPS 1: RA---SIN 512 08 33 49.342 256.00 -0.000250 0.00
AIPS 1: DEC--SIN 512 24 00 41.714 257.00 0.000250 -0.24
AIPS 1: FREQ 1 4.3432000E+10 1.00 6.4000000E+07 0.00
AIPS 1: STOKES 1 -2.0000000E+00 1.00 1.0000000E+00 0.00
AIPS 1: ----------------------------------------------------------------
AIPS 1: Coordinate equinox 2000.00
AIPS 1: Map type=NORMAL Number of iterations= 1
AIPS 1: Conv size= 0.00473 X 0.00162 Position angle= 58.26
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type CC is 1
>
In any case, I have tried the suggestion and I have run EPOSW a couple of times on the image (to get B1950 and get back J2000). I still get the same error message. Alternatively, I copied the uvdata and run UVFIX on it, this time being able to run FRING, giving very poor solution results (>60% solution failed), whereas simple fringe fitting without including the source structure finds almost all solutions.
This leads me to think that the coordinates are not properly set-up in either the uvdata or the image. But this is quite puzzling, as image has been done from the dataset and, if one coordinate shall be wrong, so the other too.
Below are the inputs for FRING. I have tried this with a couple of (point-like) sources, either including clean components or gaussian fit as model in the image CC table. My version of AIPS is 31DEC12. Do you have any suggestions about this issue?
Thank you very much,
J. C. Algaba
AIPS 1: FRING: Task to fringe fit data
AIPS 1: Adverbs Values Comments
AIPS 1: ----------------------------------------------------------------
AIPS 1: Input uv data.
AIPS 1: INNAME 'IMOGABA9F043' UV file name (name)
AIPS 1: INCLASS 'UVDATA' UV file name (class)
AIPS 1: INSEQ 1 UV file name (seq. #)
AIPS 1: INDISK 1 UV file disk drive #
AIPS 1: Data selection (multisource):
AIPS 1: CALSOUR '0827+243' Calibrator sources
AIPS 1: *rest ' '
AIPS 1: QUAL -1 Calibrator 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: TIMERANG *all 0 Time range to use.
AIPS 1: BCHAN 1 Lowest channel number 0=>all
AIPS 1: ECHAN 0 Highest channel number
AIPS 1: CHINC 1 Spectral channel increment
AIPS 1: (data are averaged)
AIPS 1: ANTENNAS *all 0 Antennas to select. 0=all
AIPS 1: DOFIT *all 0 Subset of ANTENNAS list
AIPS 1: for which solns are desired.
AIPS 1: SUBARRAY 0 Subarray, 0=>all
AIPS 1: UVRANGE 0 0 Range of uv distance for full
AIPS 1: weight
AIPS 1: WTUV 0 Weight outside UVRANGE 0=0.
AIPS 1: WEIGHTIT 0 Modify data weights function
AIPS 1:
AIPS 1: Cal. info for input:
AIPS 1: DOCALIB 1 > 0 calibrate data & weights
AIPS 1: > 99 do NOT calibrate weights
AIPS 1: GAINUSE 4 CL table to apply.
AIPS 1: DOPOL -1 If >0 correct polarization.
AIPS 1: PDVER 0 PD table to apply (DOPOL>0)
AIPS 1: BLVER -1 BL table to apply.
AIPS 1: FLAGVER 0 Flag table version
AIPS 1: DOBAND 0 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:
AIPS 1: CLEAN map (optional)
AIPS 1: IN2NAME '0827+243' Cleaned map name (name)
AIPS 1: IN2CLASS 'ICC043' Cleaned map name (class)
AIPS 1: IN2SEQ 1 Cleaned map name (seq. #)
AIPS 1: IN2DISK 1 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 0 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: SMODEL *all 0 Source model, 1=flux,2=x,3=y
AIPS 1: See HELP SMODEL for models.
AIPS 1:
AIPS 1: Output uv data file.
AIPS 1: OUTNAME ' ' UV file name (name)
AIPS 1: OUTCLASS ' ' UV file name (class)
AIPS 1: OUTSEQ 0 UV file name (seq. #)
AIPS 1: OUTDISK 0 UV file disk drive #
AIPS 1:
AIPS 1: Solution control adverbs:
AIPS 1: REFANT 1 Reference antenna
AIPS 1: SEARCH *all 0 Prioritized reference antenna
AIPS 1: list - supplements REFANT
AIPS 1: - but only if APARM(9)>0
AIPS 1: SOLINT 0.5 Solution interval (min)
AIPS 1: 0 => 10 min
AIPS 1: SOLSUB 0 Solution subinterval
AIPS 1: SOLMIN 0 Min solution interval
AIPS 1:
AIPS 1: APARM 2 0 General parameters
AIPS 1: 0 1 1=min. no. antennas
AIPS 1: 1 2 2 > 0 => data divided
AIPS 1: 5 *rest 0 3 > 0 => avg. RR,LL
AIPS 1: 4 > 0 => avg. freq. in IFs
AIPS 1: 5 = 1 => combine all IFs
AIPS 1: = 2 => also MB delay
AIPS 1: = 3 => combine IFs in
AIPS 1: halves
AIPS 1: = 4 => combine IFs in
AIPS 1: quarters
AIPS 1: 6=print level, 1=some
AIPS 1: 7=SNR cutoff (0=>5)
AIPS 1: 8=max. ant. # (no AN)
AIPS 1: 9 > 0 => do exhaustive
AIPS 1: baseline search
AIPS 1: DPARM 0 100 Delay-rate parameters
AIPS 1: 100 *rest 0 1=no. bl combo. (def=3)
AIPS 1: 2=delay win (nsec), if <0
AIPS 1: no delay search done
AIPS 1: 3=rate win (mHz)
AIPS 1: 4=int. time (sec)
AIPS 1: 0 => min. found in data
AIPS 1: 5 >0 => don't do ls. soln
AIPS 1: 6 >0 => don't avg. in freq
AIPS 1: 7 >0 => don't rereference
AIPS 1: phase
AIPS 1: 8 > 0 => activate zero'ing
AIPS 1: options
AIPS 1: 9 > 0 => do not fit rate
AIPS 1: SNVER 5 Output SN table, 0=>new table
AIPS 1: ANTWT *all 0 Ant. weights (0=>1.0)
AIPS 1: BADDISK *all 0 Disk no. not to use for
AIPS 1: scratch files.
>
>go
AIPS 1: Resumes
AIPS 1: RETURN CODE 5 RECEIVED: STOPPING
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://listmgr.nrao.edu/pipermail/daip/attachments/20150128/efa09451/attachment-0001.html>
More information about the Daip
mailing list