 News |
FAQ
|
| Search
|
Home
|
| Getting Started | Documentation | Glish | Learn More | Programming | Contact Us |
|
|
|
Here we show an example run of imager that makes an image, uses it to calculate a calibration table, and then prints out the fit acheived for each time range.
aips++
AIPS++ will use .glishrc file in home directory
AIPSPATH is /aips++/daily sun4sol_gnu socorro tarzan
*********************** Welcome to AIPS++ ***************************
Copyright (C) 1996,1997,1999,2000 Associated Universities, Inc. Washington DC, USA.
defaultlogger available for use, can abbreviate to dl
Server started: misc
Compiled at AIPS++ version: 07.453.00 Thu 1997/01/02 01:30:10 GMT
defaultservers available for user
defaultutility ready for use, can abbreviate to du
defaultplotter ready for use, can abbreviate to dp
table system ready for use
defaultdisplay ready for use, can abbreviate to dd
Starting executable : imager
DO server for imager
defaultimager ready for use, can abbreviate to di or img
AIPS++ packages available and ready for use: aips and synthesis
Initializing help system....
help system initialized
help() for help
Glish version 2.5.0.13.
- di.initialize("3C273XC1")
F
- di.inputs()
# Inputs:
data := [format='uvfits', file='3C273XC1.FITS']
MS := 'MS'
images := [model='model', derivedmodel='model', restored='restored',
residual='residual', psf='psf']
cell := [ 0, 0 ]
npix := [ 128, 128 ]
npol := 1
nchan := 1
wt := [type='uniform', rmode='norm', noise=1.0, robust=0.0]
clean := [niter=1000, threshold=0, gain=0.1, mask='']
nnls := [niter=0, tolerance=1e-06, fluxmask='', datamask='']
beam := [ 0, 0, 0 ]
model := [zero=[ F, F, F, T ], constrain=F, threshold=0, mask='']
mask := [interactive=F, threshold=0]
flags := [threshold=0]
T := [t=0, solve=F, phaseonly=T, initialtable='', table='tcal']
G := [t=0, solve=F, phaseonly=F, initialtable='', table='gcal']
D := [t=0, solve=F, phaseonly=F, initialtable='', table='dcal']
P := [t=0]
M := [t=0, solve=F, initialtable='', table='mcal']
T
- di.init_ms()
Filling MS MS from uvfits file 3C273XC1.FITS
uvfitsfiller: Version 2 - UV Fits to MS filler
Reading and writing 7669 visibility groups
Found binary table of type AIPS FQ following data
Found binary table of type AIPS AN following data
Found binary table of type AIPS BL following data
Table skipped (duplicate or unrecognized)
Done
T
- di.cell:=[0.7,0.7]
- di.npol:=4
- di.T.solve:=T
- di.T.t:=60.0
- di.inputs()
# Inputs:
data := [format='uvfits', file='3C273XC1.FITS']
MS := 'MS'
images := [model='model', derivedmodel='model', restored='restored',
residual='residual', psf='psf']
cell := [ 0.7, 0.7 ]
npix := [ 128, 128 ]
npol := 4
nchan := 1
wt := [type='uniform', rmode='norm', noise=1.0, robust=0.0]
clean := [niter=1000, threshold=0, gain=0.1, mask='']
nnls := [niter=0, tolerance=1e-06, fluxmask='', datamask='']
beam := [ 0, 0, 0 ]
model := [zero=[ F, F, F, T ], constrain=F, threshold=0, mask='']
mask := [interactive=F, threshold=0]
flags := [threshold=0]
T := [t=60, solve=T, phaseonly=T, initialtable='', table='tcal']
G := [t=0, solve=F, phaseonly=F, initialtable='', table='gcal']
B := [t=0, solve=F, phaseonly=F, initialtable='', table='bcal']
D := [t=0, solve=F, phaseonly=F, initialtable='', table='dcal']
P := [t=0]
M := [t=0, solve=F, initialtable='', table='mcal']
T
- di.deconvolve()
Solving for model image
Making template image model
setting image size
making Stokes I,Q,U,V image
status:
MeasurementSet = MS
Image size = 128 by 128
Number of polarizations = 4
Number of channels = 1
Cell size = 0.7 , 0.7
Shift = 0 , 0
constructing and writing image model
Initializing calibration
Initializing solveable atmospheric gain/transmission (T-matrix)
For interval of 60 seconds, found 8 slots
Setting mode phase for T
Storing T matrix in table tcal
Applying calibration
Hogbom Clean without mask
solving for image
making initial residual image
Processing channel 0
Iteration 0 peak is [29.7735, -0.882242, 0.428236, -0.000704129] at 64,64
Iteration 100 peak is [0.178827, -0.0181872, 0.0147749, 0.000874124] at 82,42
Iteration 200 peak is [0.0158964, 0.0225238, 0.0519298, -0.00174803] at 75,58
Iteration 300 peak is [0.00231173, 0.0022082, -0.0501237, -0.001432] at 55,48
Iteration 400 peak is [0.000466952, 0.0199586, 0.0380751, -0.000104814] at 37,38
Iteration 500 peak is [-0.0190946, -0.0473223, -0.0321825, -0.00015517] at 62,61
Iteration 600 peak is [-0.0134061, -0.00958348, 0.046338, -0.000698718] at 59,70
Iteration 700 peak is [0.00182738, -0.0197189, -0.021323, 0.000256499] at 54,94
Iteration 800 peak is [-0.000505099, 0.00379298, 0.0284649, 0.00100731] at 83,72
Iteration 900 peak is [0.00213773, -0.00242627, 0.0241637, 0.000277366] at 64,53
Final iteration 1000 peak is [0.0075263, 0.00320095, -0.0189989, 7.1419e-05] at 64,64
Found Sky Model [user: 19.77 system: 0.76 real: 22.9]
Making restored and residual images restored residual
Fitted Beam: [3.7406, 3.26098, 32.7928]
T
- di.display()
imager: Displaying restored
T
- di.solve_cal()
Applying calibration
Solving for calibration in MS MS
Setting mode phase for T
Solving for T
T Jones Slot 0: 1989/06/27/01:03:20 to 01:04:10
T Jones Initial fit per unit weight = 1.54623 Jy, sum of weights = 3588
T Jones Final fit per unit weight = 1.54575 Jy, after 1 iterations
T Jones Slot 1: 1989/06/27/01:04:20 to 01:04:20
T Jones Initial fit per unit weight = 1.53308 Jy, sum of weights = 876
T Jones Final fit per unit weight = 1.53244 Jy, after 2 iterations
T Jones Slot 2: 1989/06/27/03:31:10 to 03:31:20
T Jones Initial fit per unit weight = 0.723268 Jy, sum of weights = 2028
T Jones Final fit per unit weight = 0.721122 Jy, after 2 iterations
T Jones Slot 3: 1989/06/27/03:31:30 to 03:32:10
T Jones Initial fit per unit weight = 0.710631 Jy, sum of weights = 5196
T Jones Final fit per unit weight = 0.708408 Jy, after 2 iterations
T Jones Slot 4: 1989/06/27/03:32:20 to 03:32:20
T Jones Initial fit per unit weight = 0.707964 Jy, sum of weights = 1048
T Jones Final fit per unit weight = 0.705133 Jy, after 2 iterations
T Jones Slot 5: 1989/06/27/04:43:20 to 04:43:20
T Jones Initial fit per unit weight = 0.889494 Jy, sum of weights = 1092
T Jones Final fit per unit weight = 0.88845 Jy, after 2 iterations
T Jones Slot 6: 1989/06/27/04:43:30 to 04:44:10
T Jones Initial fit per unit weight = 0.896663 Jy, sum of weights = 4236
T Jones Final fit per unit weight = 0.895928 Jy, after 2 iterations
T Jones Slot 7: 1989/06/27/04:44:20 to 04:44:30
T Jones Initial fit per unit weight = 0.864949 Jy, sum of weights = 2068
T Jones Final fit per unit weight = 0.864008 Jy, after 2 iterations
Storing T matrix in table tcal
T
- gaintable(di.T.table).fit()
successful readonly open of table tcal : 6 columns, 8 rows
[1.54575 1.53244 0.721122 0.708408 0.705133 0.88845 0.895928 0.864008]