Differences
This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision | ||
public:meetings:glowbi2008:tutorial_monday [2008-11-11 17:50] – more James M Anderson | public:meetings:glowbi2008:tutorial_monday [2017-03-08 15:27] (current) – external edit 127.0.0.1 | ||
---|---|---|---|
Line 5: | Line 5: | ||
This tutorial session is intended to get people familiar with the basic concepts of radio interferometry. | This tutorial session is intended to get people familiar with the basic concepts of radio interferometry. | ||
+ | |||
+ | |||
+ | History: | ||
+ | 2008 Nov ?? Initial version | ||
+ | 2010 NOv 14 Update for 31DEC10 AIPS | ||
====== Step 1 --- Figure Out What to Observe (Reduce) ====== | ====== Step 1 --- Figure Out What to Observe (Reduce) ====== | ||
Line 29: | Line 34: | ||
'' | '' | ||
+ | You should be able to download them from | ||
+ | [[ftp:// | ||
+ | [[ftp:// | ||
+ | |||
+ | While you are there, grab the new set of commands from | ||
+ | [[ftp:// | ||
====== Step 3 --- Start up AIPS ====== | ====== Step 3 --- Start up AIPS ====== | ||
Line 44: | Line 55: | ||
+ | '' | ||
+ | dowait=true\\ | ||
+ | dohist=1\\ | ||
+ | docrt=132\\ | ||
+ | dotv=1\\ | ||
+ | '' | ||
====== Step 4 --- Initial look at data ====== | ====== Step 4 --- Initial look at data ====== | ||
Line 54: | Line 71: | ||
'' | '' | ||
- | dowait=true\\ | ||
task ' | task ' | ||
- | default\\ | ||
datain = ' | datain = ' | ||
nfiles=0\\ | nfiles=0\\ | ||
ncount=2\\ | ncount=2\\ | ||
vlamode=' | vlamode=' | ||
+ | band ' | ||
doweight=10\\ | doweight=10\\ | ||
doconcat=true\\ | doconcat=true\\ | ||
Line 67: | Line 83: | ||
cparm(4)=28\\ | cparm(4)=28\\ | ||
cparm(8)=10./ | cparm(8)=10./ | ||
+ | dparm 0\\ | ||
+ | bparm 0\\ | ||
+ | timer 0\\ | ||
+ | calcode ' '\\ | ||
+ | clron\\ | ||
+ | outdisk 1\\ | ||
go\\ | go\\ | ||
'' | '' | ||
Line 72: | Line 94: | ||
This sets us up to read the data, sets the VLA mode to Solar, tells AIPS not to change the source if the position appears to be moving (which the Sun does), sets the shadowing limit to 28 meters, and sets the CL table interval to 10 seconds. | This sets us up to read the data, sets the VLA mode to Solar, tells AIPS not to change the source if the position appears to be moving (which the Sun does), sets the shadowing limit to 28 meters, and sets the CL table interval to 10 seconds. | ||
+ | ===== Header information ===== | ||
- | ===== LISTR --- scan listing ===== | ||
- | |||
- | First, let's check out which file we need. | ||
'' | '' | ||
- | pcat\\ | + | indisk 1\\ |
+ | pcat | ||
'' | '' | ||
- | |||
< | < | ||
AIPS 1: Catalog on disk 1 | AIPS 1: Catalog on disk 1 | ||
AIPS 1: Cat Usid Mapname | AIPS 1: Cat Usid Mapname | ||
- | AIPS 1: 1 100 19810926 | + | AIPS 1: 1 100 19810926 |
- | AIPS 1: 2 100 19810926 | + | |
- | AIPS 1: 3 100 19810926 | + | |
- | AIPS 1: | + | |
</ | </ | ||
'' | '' | ||
- | default listr\\ | + | getn 1\\ |
+ | imhe | ||
+ | '' | ||
+ | |||
+ | |||
+ | |||
+ | < | ||
+ | AIPS 1: Image=MULTI | ||
+ | AIPS 1: Telescope=VLA | ||
+ | AIPS 1: Observer=GD | ||
+ | AIPS 1: Observ. date=26-SEP-1981 | ||
+ | AIPS 1: # visibilities | ||
+ | AIPS 1: Rand axes: UU-L-SIN | ||
+ | AIPS 1: SOURCE | ||
+ | AIPS 1: ---------------------------------------------------------------- | ||
+ | AIPS 1: Type Pixels | ||
+ | AIPS 1: COMPLEX | ||
+ | AIPS 1: STOKES | ||
+ | AIPS 1: FREQ | ||
+ | AIPS 1: IF | ||
+ | AIPS 1: RA | ||
+ | AIPS 1: DEC 1 00 00 00.000 | ||
+ | AIPS 1: ---------------------------------------------------------------- | ||
+ | AIPS 1: Coordinate equinox | ||
+ | AIPS 1: Maximum version number of extension files of type HI 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 NX is 1 | ||
+ | AIPS 1: Maximum version number of extension files of type SU is 1 | ||
+ | AIPS 1: Maximum version number of extension files of type FQ is 1 | ||
+ | AIPS 1: Maximum version number of extension files of type CL is 1 | ||
+ | AIPS 1: Maximum version number of extension files of type TY is 1 | ||
+ | AIPS 1: Maximum version number of extension files of type WX is 1 | ||
+ | AIPS 1: Maximum version number of extension files of type OF is 1 | ||
+ | AIPS 1: Keyword = ' | ||
+ | AIPS 1: Keyword = ' | ||
+ | AIPS 1: Keyword = ' | ||
+ | </ | ||
+ | |||
+ | ===== LISTR --- scan listing ===== | ||
+ | |||
+ | |||
+ | '' | ||
task ' | task ' | ||
indi 1\\ | indi 1\\ | ||
- | getn 2\\ | + | getn 1\\ |
optype ' | optype ' | ||
docrt=132\\ | docrt=132\\ | ||
+ | flagver 0\\ | ||
+ | sources ' '\\ | ||
+ | stokes ' '\\ | ||
+ | docalib 0\\ | ||
+ | gainuse 1\\ | ||
+ | dopol -1\\ | ||
+ | dparm 0\\ | ||
go\\ | go\\ | ||
'' | '' | ||
< | < | ||
- | AIPS 1: Got(1) | + | |
- | vlb054 | + | |
File = 19810926 | File = 19810926 | ||
- | Freq = 1.446150006 GHz Ncor = 4 No. vis = | + | Freq = 1.446150006 GHz Ncor = 4 No. vis = |
Scan summary listing | Scan summary listing | ||
Scan Source | Scan Source | ||
- | 1 1148-001 | + | 1 1148-001 |
- | 2 SUN : 0000 | + | 2 1148-001 |
- | 3 1148-001 | + | 3 SUN : 0000 |
- | 4 SUN : 0000 | + | 4 SUN : 0000 |
- | 5 1148-001 | + | |
- | 6 SUN : 0000 | + | 6 SUN : 0000 |
- | 7 1148-001 | + | 7 SUN : 0000 |
- | 8 SUN : 0000 | + | |
- | 9 1148-001 | + | 9 SUN : 0000 |
- | | + | 10 SUN : 0000 |
- | | + | 11 1148-001 |
- | | + | |
- | | + | 13 SUN : 0000 |
- | | + | 14 1148-001 |
- | 15 1148-001 | + | |
- | | + | |
- | | + | 17 1148-001 |
- | | + | |
- | | + | |
- | | + | 20 1148-001 |
- | | + | |
- | | + | |
- | | + | 23 1148-001 |
- | | + | |
- | | + | |
- | | + | |
- | | + | |
- | | + | 28 SUN : 0000 |
- | | + | |
- | | + | |
- | | + | 31 SUN : 0000 |
- | | + | |
- | | + | |
- | | + | |
- | | + | 35 1148-001 |
- | | + | |
- | | + | |
- | | + | 38 1148-001 |
- | | + | |
- | | + | |
- | | + | 41 1148-001 |
- | | + | |
- | | + | |
- | | + | 44 1148-001 |
- | | + | |
- | | + | |
- | | + | 47 1148-001 |
- | | + | |
- | | + | |
- | | + | 50 1148-001 |
- | | + | |
- | | + | |
- | | + | 53 1148-001 |
+ | | ||
+ | | ||
+ | 56 1148-001 | ||
+ | | ||
+ | | ||
+ | 59 1148-001 | ||
+ | | ||
+ | | ||
+ | 62 1148-001 | ||
+ | | ||
+ | | ||
+ | 65 1148-001 | ||
+ | | ||
+ | | ||
+ | 68 1148-001 | ||
+ | | ||
+ | | ||
+ | 71 1148-001 | ||
+ | | ||
+ | | ||
+ | 74 1148-001 | ||
+ | | ||
+ | | ||
+ | | ||
Source summary | Source summary | ||
Line 168: | Line 256: | ||
ID Source | ID Source | ||
- | 1 1148-001 | + | 1 1148-001 |
- | 2 SUN : 0000 | + | 2 SUN : 0000 |
- | 3 3C286 : 0000 | + | 3 3C286 : 0000 |
ID Source | ID Source | ||
Line 181: | Line 269: | ||
</ | </ | ||
+ | |||
+ | ===== PRTAN ===== | ||
+ | |||
+ | Print the Antenna positions --- useful for thinking about calibration. | ||
'' | '' | ||
- | getn 4\\ | + | go prtan\\ |
- | go\\ | + | |
'' | '' | ||
< | < | ||
- | AIPS 1: Got(1) | + | |
- | vlb054 | + | File=19810926 |
- | File = 19810926 | + | Array= VLA |
- | Freq = | + | |
- | Scan summary listing | + | |
- | Scan Source | + | Array reference position in meters (Earth centered) |
- | 1 SUN : 0000 | + | Array BX= -1601185.36500 |
- | 2 SUN : 0000 | + | Polar X = 0.00000 Polar Y = 0.00000 arcsec |
- | | + | Earth rotation rate = 360.9856449713 degrees |
- | 4 SUN : 0000 | + | GST at UT=0 = 364.7139688925 degrees |
- | 5 SUN : 0000 | + | UT1-UTC= 0.0000000 |
- | 6 SUN : 0000 | + | Solutions not yet determined for a particular FREQID |
- | 7 SUN : 0000 | + | |
- | 8 SUN : 0000 | + | |
- | 9 SUN : 0000 | + | |
- | 10 SUN : 0000 | + | |
- | 11 SUN : 0000 | + | |
- | 12 SUN : 0000 | + | |
- | 13 SUN : 0000 | + | |
- | 14 SUN : 0000 | + | |
- | 15 SUN : 0000 | + | |
- | 16 SUN : 0000 | + | |
- | 17 SUN : 0000 | + | |
- | 18 SUN : 0000 | + | |
- | 19 SUN : 0000 | + | |
- | 20 SUN : 0000 | + | |
- | 21 SUN : 0000 | + | |
- | 22 SUN : 0000 | + | |
- | 23 SUN : 0000 | + | |
- | 24 SUN : 0000 | + | |
- | 25 SUN : 0000 | + | |
- | Source summary | + | Ant 1 = VLA: |
- | Velocity type = ' | + | Mount=ALAZ Axis offset= |
+ | Feed polarization type = | ||
- | | + | Ant |
- | 1 SUN : 0000 | + | Mount=ALAZ |
+ | Feed polarization type = R L | ||
- | | + | Ant |
- | 1 All Sources | + | Mount=ALAZ |
+ | Feed polarization type = R L | ||
- | Frequency Table summary | + | Ant 4 = VLA: |
- | FQID IF# Freq(GHz) | + | Mount=ALAZ |
- | 1 | + | Feed polarization type = R L |
+ | |||
+ | Ant 5 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 6 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 10 = VLA:_OUT BX= 0.0000 BY= 0.0000 BZ= 0.0000 | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 11 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 12 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 13 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 14 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 15 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 16 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 17 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 18 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 19 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 20 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 21 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 22 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 23 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 24 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 25 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 26 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 27 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 28 = VLA: | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Ant 29 = VPT:_OUT BX= 0.0000 BY= 0.0000 BZ= 0.0000 | ||
+ | Mount=ALAZ | ||
+ | Feed polarization type = R L | ||
+ | |||
+ | Location Of VLA Antennas | ||
+ | |||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | ( 8) W1 | ||
+ | ( 6) W2 | ||
+ | ( 7) W3 | ||
+ | (28) W4 | ||
+ | (22) W5 | ||
+ | (27) W6 | ||
+ | (14) W7 | ||
+ | (17) W8 | ||
+ | (15) W9 | ||
+ | VLA:_OUT (10) | ||
+ | VPT:_OUT (29) | ||
AIPS 1: Resumes | AIPS 1: Resumes | ||
</ | </ | ||
+ | |||
+ | The really important part for you, the data reducer, is the antenna layout at the bottom. | ||
+ | |||
+ | |||
+ | |||
+ | ===== UVPLT ===== | ||
+ | |||
+ | Let's have a look at the information about the visibilities. | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | source ' | ||
+ | stokes ' | ||
+ | calcode ' '\\ | ||
+ | uvrange 0\\ | ||
+ | antennas 0\\ | ||
+ | basel 0\\ | ||
+ | xinc 1\\ | ||
+ | aparm 0\\ | ||
+ | bparm 0\\ | ||
+ | doweight 1\\ | ||
+ | refant 3\\ | ||
+ | do3col 1\\ | ||
+ | dotv 1\\ | ||
+ | tvinit\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | Warning: the next step will take a long time. You are free to skip this, and just look at the picture. | ||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | bparm 0 2\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | bparm 11 1\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | bparm 6 7\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | source ' | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | '' | ||
+ | source ' '\\ | ||
+ | bparm 11 15\\ | ||
+ | xinc 50\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | '' | ||
+ | bparm 11 18\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | ===== LISTR ===== | ||
+ | |||
+ | Now let's look at the raw visibility numbers. | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | xinc 1\\ | ||
+ | task ' | ||
+ | optype ' | ||
+ | source ' | ||
+ | stokes ' '\\ | ||
+ | bif 0\\ | ||
+ | eif 0\\ | ||
+ | flagver 0\\ | ||
+ | dparm 0\\ | ||
+ | antenna 3,0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | < | ||
+ | | ||
+ | File = 19810926 | ||
+ | Freq= 1.446150006 GHz | ||
+ | Stokes = RR | ||
+ | |||
+ | |||
+ | Source=3C286 | ||
+ | Flux = 0.0000 Jy, Calcode = B , Freq = 1.446150006 GHz | ||
+ | Amplitudes, 1000 = 1.000 Jy, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | Amplitudes, 1000 = 10.000 Jy, averging type = Vector | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | AIPS 1: Resumes | ||
+ | </ | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | ===== Flagging ===== | ||
+ | |||
+ | |||
+ | Just enter these commands. | ||
+ | needed to be applied. | ||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | antenna 23,0\\ | ||
+ | basel 0\\ | ||
+ | timer 0\\ | ||
+ | aparm 0\\ | ||
+ | opcode ' | ||
+ | reason 'bad ant'\\ | ||
+ | dohist 1\\ | ||
+ | go\\ | ||
+ | antenna 26,0\\ | ||
+ | timer 0 18 29 0 0 20 7 20\\ | ||
+ | go\\ | ||
+ | timer 0 0 0 0 0 15 8 10\\ | ||
+ | go\\ | ||
+ | timer 0 21 48 30 0 22 37 10\\ | ||
+ | go\\ | ||
+ | antenna 25,0\\ | ||
+ | timer 0 22 38 20 1 0 21 50\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | ===== More LISTR, now on the phase calibrator ===== | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | optype ' | ||
+ | source ' | ||
+ | stokes ' '\\ | ||
+ | calcode ' '\\ | ||
+ | bif 0\\ | ||
+ | eif 0\\ | ||
+ | timer 0\\ | ||
+ | flagver 0\\ | ||
+ | dparm 0\\ | ||
+ | antenna 3,0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | |||
+ | < | ||
+ | | ||
+ | File = 19810926 | ||
+ | Freq= 1.446150006 GHz | ||
+ | Stokes = RR | ||
+ | Applying flag table 1 | ||
+ | |||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Amplitudes, 1000 = 10.000 Jy, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Amplitudes, 1000 = 10.000 Jy, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Amplitudes, 1000 = 10.000 Jy, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Amplitudes, 1000 = 10.000 Jy, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | </ | ||
+ | |||
+ | |||
+ | '' | ||
+ | dparm 1,0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | < | ||
+ | | ||
+ | File = 19810926 | ||
+ | Freq= 1.446150006 GHz | ||
+ | Stokes = RR | ||
+ | Applying flag table 1 | ||
+ | |||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Phase, 1000 = 1000.00 degrees, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Phase, 1000 = 1000.00 degrees, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Phase, 1000 = 1000.00 degrees, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | |||
+ | Source=1148-001 | ||
+ | Flux = 0.0000 Jy, Calcode = C , Freq = 1.446150006 GHz | ||
+ | Phase, 1000 = 1000.00 degrees, averging type = Vector | ||
+ | |||
+ | Baselines | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | | ||
+ | </ | ||
+ | |||
+ | ====== Calibration ====== | ||
+ | |||
+ | |||
+ | ===== SETJY ===== | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | source ' | ||
+ | optype ' | ||
+ | aparm 0\\ | ||
+ | aparm(2)=4\\ | ||
+ | zerosp 0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | < | ||
+ | vlb054> SETJY1: Task SETJY (release of 31DEC10) begins | ||
+ | vlb054> SETJY1: **WARNING: OPCODE=CALC AND FREQID = -1 | ||
+ | vlb054> SETJY1: | ||
+ | vlb054> SETJY1: A source model for this calibrator is available | ||
+ | vlb054> SETJY1: Consult the help file for CALRD for assistance | ||
+ | vlb054> SETJY1: / Flux calculated using known spectrum | ||
+ | vlb054> SETJY1: BIF = 1 EIF = 1 /Range of IFs | ||
+ | vlb054> SETJY1: ' | ||
+ | vlb054> SETJY1: / Using (1995.2) VLA or Reynolds (1934-638) coefficients | ||
+ | vlb054> SETJY1: Appears to have ended successfully | ||
+ | vlb054> SETJY1: vlb054 | ||
+ | </ | ||
+ | |||
+ | ===== CALRD ===== | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | object ' | ||
+ | band ' | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | < | ||
+ | vlb054> CALRD1: Task CALRD (release of 31DEC08) begins | ||
+ | vlb054> CALRD1: Reading disk file AIPSTARS: | ||
+ | vlb054> CALRD1: Create 3C286_L | ||
+ | vlb054> CALRD1: Appears to have ended successfully | ||
+ | vlb054> CALRD1: vlb054 | ||
+ | </ | ||
+ | |||
+ | So it looks like the calibration data image was put into catalog entry number 2. | ||
+ | |||
+ | Let's have a look at it. | ||
+ | |||
+ | |||
+ | '' | ||
+ | getn 2\\ | ||
+ | tvinit\\ | ||
+ | tvlod\\ | ||
+ | tvfiddle\\ | ||
+ | '' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | ===== UVFIX ===== | ||
+ | |||
+ | The sky positions (RA and Dec) were entered as coordinates of date, because the Sun is a moving object. | ||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | getn 1\\ | ||
+ | clron\\ | ||
+ | shift 0\\ | ||
+ | uvfixprm 0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | UVFIX change in April means we have to regenerate NX table\\ | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | getn 3\\ | ||
+ | cparm 0\\ | ||
+ | cparm(3)= -1\\ | ||
+ | bparm 0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | |||
+ | ===== CALIB ===== | ||
+ | |||
+ | First, we do the amplitude calibrator. | ||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | getn 3\\ | ||
+ | calsour ' | ||
+ | uvrange 0\\ | ||
+ | antennas 0\\ | ||
+ | timer 0\\ | ||
+ | refant 24\\ | ||
+ | weightit 1\\ | ||
+ | in2di 1\\ | ||
+ | get2n 2\\ | ||
+ | ncomp 0\\ | ||
+ | solmode ' | ||
+ | aparm 0\\ | ||
+ | aparm(6) 2\\ | ||
+ | minamper 10\\ | ||
+ | minphser 10\\ | ||
+ | docalib 1\\ | ||
+ | gainuse 1\\ | ||
+ | solint 30\\ | ||
+ | solsub 2\\ | ||
+ | cparm 0\\ | ||
+ | cparm(3) 10\\ | ||
+ | cparm(4) 10\\ | ||
+ | soltype ' '\\ | ||
+ | clrmsg\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | Next, the phase calibrator. | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | calsour ' | ||
+ | clr2n\\ | ||
+ | uvrange 1.5,0\\ | ||
+ | solint 0\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | |||
+ | ===== SNPLT ===== | ||
+ | |||
+ | Look at the amplitude calibrator results. | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | source ' '\\ | ||
+ | inver 1\\ | ||
+ | inext ' | ||
+ | stokes ' '\\ | ||
+ | optype ' | ||
+ | opcode ' | ||
+ | xinc 1\\ | ||
+ | nplots 9\\ | ||
+ | antenna 0\\ | ||
+ | timer 0\\ | ||
+ | tvinit\\ | ||
+ | dotv 1\\ | ||
+ | go\\ | ||
+ | optype ' | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | Look at the phase calibrator results. | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | inver 2\\ | ||
+ | optype ' | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | optype ' | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | ===== GETJY ===== | ||
+ | |||
+ | Now we need to transfer the amplitude calibration information from the amplitude calibrator to the phase calibrator. | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | sources ' | ||
+ | calsour ' | ||
+ | calcode ' '\\ | ||
+ | bif 0\\ | ||
+ | eif 0\\ | ||
+ | antenna 0\\ | ||
+ | timer 0\\ | ||
+ | snver 0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | < | ||
+ | vlb054> GETJY1: Task GETJY (release of 31DEC08) begins | ||
+ | vlb054> GETJY1: | ||
+ | vlb054> GETJY1: 1148-001 | ||
+ | vlb054> GETJY1: Appears to have ended successfully | ||
+ | vlb054> GETJY1: vlb054 | ||
+ | </ | ||
+ | |||
+ | 2.9 Jy agrees relatively well with the rough value given by the VLA calibrator list. | ||
+ | |||
+ | < | ||
+ | 1150-003 | ||
+ | 1148-001 | ||
+ | ----------------------------------------------------- | ||
+ | BAND A B C D FLUX(Jy) | ||
+ | ===================================================== | ||
+ | | ||
+ | 6cm C P P P P 1.92 | ||
+ | 3.7cm X P P P P | ||
+ | 2cm U P P P P 1.40 | ||
+ | 1.3cm K S S S S | ||
+ | 0.7cm Q W W W W | ||
+ | </ | ||
+ | |||
+ | ===== Redo the CALIB calibration for the phase calibrator ===== | ||
+ | |||
+ | |||
+ | '' | ||
+ | inext ' | ||
+ | inver 2\\ | ||
+ | extdest\\ | ||
+ | tget calib\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | |||
+ | ===== Apply the calibration to all sources, both amplitude and phase ===== | ||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | getn 3\\ | ||
+ | sources ' '\\ | ||
+ | calsour ' | ||
+ | calcode ' '\\ | ||
+ | opcode ' | ||
+ | gainver 1\\ | ||
+ | gainuse 2\\ | ||
+ | refant 24\\ | ||
+ | bparm 0\\ | ||
+ | interpol ' | ||
+ | timer 0\\ | ||
+ | antenna 0\\ | ||
+ | doblank 1\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | ===== Check the calibration results ===== | ||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | source ' '\\ | ||
+ | inver 2\\ | ||
+ | inext ' | ||
+ | stokes ' '\\ | ||
+ | optype ' | ||
+ | opcode ' | ||
+ | xinc 1\\ | ||
+ | nplots 9\\ | ||
+ | antenna 0\\ | ||
+ | timer 0\\ | ||
+ | tvinit\\ | ||
+ | dotv 1\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | '' | ||
+ | optype ' | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | ====== Test Imaging ====== | ||
+ | |||
+ | ===== Amplitude Calibrator ===== | ||
+ | |||
+ | |||
+ | Ok, let's start with a few test images. | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | getn 3\\ | ||
+ | source ' | ||
+ | uvrange 0\\ | ||
+ | docalib 1\\ | ||
+ | gainuse 2\\ | ||
+ | stokes ' | ||
+ | cellsize 10,10\\ | ||
+ | imsize 512,512\\ | ||
+ | antenna 0\\ | ||
+ | niter 10000\\ | ||
+ | overlap 1\\ | ||
+ | DO3DIMAG 1\\ | ||
+ | rashift 1800.00,0\\ | ||
+ | decshift -920.00,0\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | After you have cleaned it a bit, this is what it looks like. | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | And this is the point spread function of the actual observations (called the dirty beam). | ||
+ | |||
+ | |||
+ | '' | ||
+ | getn 5\\ | ||
+ | tvinit\\ | ||
+ | tvlod\\ | ||
+ | tvfiddle\\ | ||
+ | '' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | ===== Phase Calibrator ===== | ||
+ | |||
+ | Ok, clean up after ourselves | ||
+ | |||
+ | |||
+ | '' | ||
+ | getn 5;zap\\ | ||
+ | getn 6;zap\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | Now do the imaging and light cleaning | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | tget imagr\\ | ||
+ | source ' | ||
+ | rashift 0\\ | ||
+ | decshift 0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | Note that in this image, there is a lot of large-scale junk. Why? | ||
+ | |||
+ | |||
+ | And here is the dirty beam. Notice how the dirty beam looks like the initial image of the phase calibrator before you began to clean it. Why? | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | getn 5\\ | ||
+ | tvinit\\ | ||
+ | tvlod\\ | ||
+ | tvfiddle\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | ====== Imaging the Sun ====== | ||
+ | |||
+ | |||
+ | Ok, to do the best job, we need to separate the Sun data from the rest of the dataset. | ||
+ | |||
+ | '' | ||
+ | getn 5;zap\\ | ||
+ | getn 6;zap\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | |||
+ | ===== SPLIT ===== | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | getn 3\\ | ||
+ | source ' | ||
+ | timer 0\\ | ||
+ | stokes ' '\\ | ||
+ | aparm 0\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | ===== IMAGR ===== | ||
+ | |||
+ | |||
+ | Now do an initial imaging and cleaning step. | ||
+ | |||
+ | |||
+ | '' | ||
+ | tget imagr\\ | ||
+ | getn 4\\ | ||
+ | source ' '\\ | ||
+ | docalib 0\\ | ||
+ | gainuse 1\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | Here are some of the stages of my cleaning. | ||
+ | |||
+ | {{: | ||
+ | {{: | ||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | ===== Self-Calibration ===== | ||
+ | |||
+ | Ok, let's do some self-calibration. | ||
+ | |||
+ | We run exactly the same calibration software, it exactly the same manner. | ||
+ | |||
+ | |||
+ | '' | ||
+ | tget calib\\ | ||
+ | getn 4\\ | ||
+ | get2n 7\\ | ||
+ | solint 30/60\\ | ||
+ | solmode ' | ||
+ | docalib 0\\ | ||
+ | uvrange 0.5,0\\ | ||
+ | calsour ' '\\ | ||
+ | inver 0\\ | ||
+ | snver 0\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | |||
+ | ==== Check the results ==== | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | inext ' | ||
+ | inver 1\\ | ||
+ | stokes ' '\\ | ||
+ | optype ' | ||
+ | tvinit\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | |||
+ | Note some of the really large phase deviations! | ||
+ | |||
+ | ===== Image again, using the new calibration ===== | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | getn 5;zap\\ | ||
+ | getn 6;zap\\ | ||
+ | getn 7;zap\\ | ||
+ | tget imagr\\ | ||
+ | docalib 1\\ | ||
+ | gainuse 1\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | {{: | ||
+ | |||
+ | This one is a bit better. | ||
+ | |||
+ | |||
+ | |||
+ | ===== More Self-Calibration ===== | ||
+ | But we can do more. Now let's do both amplitude and phase calibration. | ||
+ | |||
+ | ==== Phase Self-Cal ==== | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | tget calib\\ | ||
+ | getn 4\\ | ||
+ | get2n 7\\ | ||
+ | solint 20/60\\ | ||
+ | solmode ' | ||
+ | docalib 0\\ | ||
+ | uvrange 0.5,0\\ | ||
+ | calsour ' '\\ | ||
+ | inver 0\\ | ||
+ | snver 0\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | ==== Check ==== | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | inext ' | ||
+ | inver 2\\ | ||
+ | stokes ' '\\ | ||
+ | optype ' | ||
+ | tvinit\\ | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | ==== Amplitude Self-Cal ==== | ||
+ | |||
+ | Doing both phase and amplitude calibration requires far more from the data S/N. It is best to have a longer time interval for amplitude calibration. | ||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | tget calib\\ | ||
+ | getn 5\\ | ||
+ | solint 20\\ | ||
+ | solmode ' | ||
+ | cparm(2) 1\\ | ||
+ | go\\ | ||
+ | clrmsg\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | ==== Check ==== | ||
+ | |||
+ | |||
+ | '' | ||
+ | task ' | ||
+ | inext ' | ||
+ | inver 1\\ | ||
+ | stokes ' '\\ | ||
+ | optype ' | ||
+ | tvinit\\ | ||
+ | go\\ | ||
+ | optype ' | ||
+ | go\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | ==== Cleanup ==== | ||
+ | |||
+ | Clean up the files. | ||
+ | |||
+ | |||
+ | '' | ||
+ | getn 6;zap\\ | ||
+ | getn 7;zap\\ | ||
+ | getn 8;zap\\ | ||
+ | '' | ||
+ | |||
+ | |||
+ | ===== Final Imaging ===== | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | |||
+ | '' | ||
+ | tget imagr\\ | ||
+ | getn 5\\ | ||
+ | docalib 1\\ | ||
+ | gainuse 1\\ | ||
+ | '' | ||
+ | |||
+ | {{: | ||
+ | |||
+ | |||
+ | |||
+ | |||
+ | ====== Discussion ====== | ||
+ | |||
+ | What have you learned? | ||
+ | |||
+ | Does your image look like this image? | ||
+ | {{http:// | ||
+ | |||
+ | What is different? | ||
+ | |||
+ | Have a look at http:// | ||