public:meetings:glowbi2008:tutorial_monday

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision		 Previous revision
public:meetings:glowbi2008:tutorial_monday [2008-11-11 18:20] – James M Andersonpublic: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.  The student will be introduced to interferometry data, visibilities, the //(u,v)// plane, calibration, and imaging.   This tutorial session is intended to get people familiar with the basic concepts of radio interferometry.  The student will be introduced to interferometry data, visibilities, the //(u,v)// plane, calibration, and imaging.  
 +
 +
 +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.mpifr-bonn.mpg.de/outgoing/anderson/GD_1|GD_1]] and 
 +[[ftp://ftp.mpifr-bonn.mpg.de/outgoing/anderson/GD_2|GD_2]].
 +
 +While you are there, grab the new set of commands from
 +[[ftp://ftp.mpifr-bonn.mpg.de/outgoing/anderson/Sun.txt|Sun.txt]].
  
 ====== 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\\ 
-dohist=1\\ 
-docrt=132\\ 
-dotv=1\\ 
 task 'fillm'\\ task 'fillm'\\
-default\\ 
 datain = 'PWD:GD_\\ datain = 'PWD:GD_\\
 nfiles=0\\ nfiles=0\\
 ncount=2\\ ncount=2\\
 vlamode='S '\\ vlamode='S '\\
 +band 'l'\\
 doweight=10\\ doweight=10\\
 doconcat=true\\ doconcat=true\\
Line 70: Line 83:
 cparm(4)=28\\ cparm(4)=28\\
 cparm(8)=10./60\\ cparm(8)=10./60\\
 +dparm 0\\
 +bparm 0\\
 +timer 0\\
 +calcode ' '\\
 +clron\\
 +outdisk 1\\
 go\\ go\\
 '' ''
Line 75: 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
 '' ''
- 
  
 <file> <file>
 AIPS 1: Catalog on disk  1 AIPS 1: Catalog on disk  1
 AIPS 1:  Cat Usid Mapname      Class   Seq  Pt     Last access      Stat AIPS 1:  Cat Usid Mapname      Class   Seq  Pt     Last access      Stat
-AIPS 1:    1  100 19810926    .C BAND.    1 UV 11-NOV-2008 18:44:06 +AIPS 1:    1  100 19810926    .L BAND.    1 UV 14-NOV-2008 21:35:31
-AIPS 1:    2  100 19810926    .L BAND.    1 UV 11-NOV-2008 18:46:44 +
-AIPS 1:    3  100 19810926    .C BAND.    2 UV 11-NOV-2008 18:40:21 +
-AIPS 1:    4  100 19810926    .L BAND.    2 UV 11-NOV-2008 18:40:21+
 </file> </file>
  
  
 '' ''
-default listr\\ +getn 1\\ 
-task 'listr'\\ +imhe
-indi 1\\ +
-getn 2\\ +
-optype 'scan'\\ +
-docrt=132\\ +
-go\\+
 '' ''
 +
 +
  
 <file> <file>
-AIPS 1: Got(1)   disk= 1  user= 100   type=UV   19810926.L BAND.1 +AIPS 1: Image=MULTI     (UV)         Filename=19810926    .L BAND.   1 
- vlb054    LISTR(31DEC08)    100     11-NOV-2008  18:46:44    Page    1 +AIPS 1: Telescope=VLA                Receiver=VLA 
-File = 19810926    .L BAND.   1 Vol = 1  Userid =  100 +AIPS 1: Observer=GD                  User #=  100 
-Freq =  1.446150006 GHz   Ncor =  4   No. vis =    278465 +AIPS 1: Observ. date=26-SEP-1981     Map date=14-NOV-2008 
-Scan summary listing +AIPS 1: # visibilities    556884     Sort order  TB 
-  +AIPS 1: Rand axes: UU-L-SIN  VV-L-SIN  WW-L-SIN  BASELINE  TIME1 
-Scan      Source      Qual  Calcode Sub         Timerange          FrqID   START VIS  END VIS +AIPS 1:            SOURCE  FREQSEL 
-   1 1148-001        : 0000  C        1  0/13:53:25 -   0/13:54:35     1       1       2195 +AIPS 1: ---------------------------------------------------------------- 
-   2 SUN             : 0000           1  0/13:59:15 -   0/13:59:35     1    2196       3150 +AIPS 1: Type    Pixels   Coord value     at Pixel     Coord incr   Rotat 
-   3 1148-001        : 0000  C        1  0/14:18:05 -   0/14:19:35     1    3151       5809 +AIPS 1: COMPLEX      3   1.0000000E+00       1.00  1.0000000E+00    0.00 
-   4 SUN             : 0000           1  0/14:20:25 -   0/14:24:35     1    5810      14662 +AIPS 1: STOKES       4  -1.0000000E+00       1.00 -1.0000000E+00    0.00 
-   5 1148-001        : 0000  C        1  0/14:43:05 -   0/14:44:35     1   14663      17320 +AIPS 1: FREQ         1   1.4461500E+09       1.00  1.2500000E+07    0.00 
-   6 SUN             : 0000           1  0/14:45:25 -   0/14:49:35     1   17321      25867 +AIPS 1: IF           1   1.0000000E+00       1.00  1.0000000E+00    0.00 
-   7 1148-001        : 0000  C        1  0/15:07:55 -   0/15:09:25     1   25868      28503 +AIPS 1: RA           1    00 00 00.000       1.00       3600.000    0.00 
-   8 SUN             : 0000           1  0/15:10:15 -   0/15:14:25     1   28504      37058 +AIPS 1: DEC          1    00 00 00.000       1.00       3600.000    0.00 
-   9 1148-001        : 0000  C        1  0/15:32:55 -   0/15:34:25     1   37059      39716 +AIPS 1: ---------------------------------------------------------------- 
-  10 SUN             : 0000           1  0/15:35:15 -   0/15:39:25     1   39717      48338 +AIPS 1: Coordinate equinox    0.00 
-  11 1148-001        : 0000  C        1  0/15:57:55 -   0/15:59:25     1   48339      50720 +AIPS 1: Maximum version number of extension files of type HI is   1 
-  12 SUN             : 0000           1  0/16:00:15 -   0/16:04:15     1   50721      58994 +AIPS 1: Maximum version number of extension files of type AN is   1 
-  13 1148-001        : 0000  C        1  0/16:22:45 -   0/16:24:15     1   58995      61584 +AIPS 1: Maximum version number of extension files of type NX is   1 
-  14 SUN             : 0000           1  0/16:25:05 -   0/16:29:15     1   61585      70334 +AIPS 1: Maximum version number of extension files of type SU is   1 
-  15 1148-001        : 0000  C        1  0/16:47:45 -   0/16:49:15     1   70335      72949 +AIPS 1: Maximum version number of extension files of type FQ is   1 
-  16 SUN             : 0000           1  0/16:50:05 -   0/16:54:15     1   72950      81664 +AIPS 1: Maximum version number of extension files of type CL is   1 
-  17 1148-001        : 0000  C        1  0/17:12:35 -   0/17:12:35     1   81665      81664 +AIPS 1: Maximum version number of extension files of type TY is   1 
-  18 1148-001        : 0000  C        1  0/17:12:45 -   0/17:14:05     1   81665      84070 +AIPS 1: Maximum version number of extension files of type WX is   1 
-  19 SUN             : 0000           1  0/17:14:55 -   0/17:19:05     1   84071      92937 +AIPS 1: Maximum version number of extension files of type OF is   1 
-  20 1148-001        : 0000  C        1  0/17:37:35 -   0/17:39:05     1   92938      95507 +AIPS 1: Keyword = 'CORRMODE'  value = '        ' 
-  21 SUN             : 0000           1  0/17:39:45 -   0/17:44:25     1   95508     105079 +AIPS 1: Keyword = 'VLAIFS  '  value = 'AC      ' 
-  22 1148-001        : 0000  C        1  0/18:02:25 -   0/18:03:55     1  105080     107535 +AIPS 1: Keyword = 'CORRCOEF'  value =           -1
-  23 SUN             : 0000           1  0/18:04:45 -   0/18:08:55     1  107536     116401 +
-  24 1148-001        : 0000  C        1  0/18:27:25 -   0/18:28:55     1  116402     119012 +
-  25 SUN             : 0000           1  0/18:29:45 -   0/18:33:55     1  119013     127842 +
-  26 1148-001        : 0000  C        1  0/18:52:35 -   0/18:53:55     1  127843     129825 +
-  27 SUN             : 0000           1  0/18:54:45 -   0/18:59:05     1  129826     137705 +
-  28 1148-001        : 0000  C        1  0/19:17:25 -   0/19:18:45     1  137706     139972 +
-  29 SUN             : 0000           1  0/19:19:35 -   0/19:23:45     1  139973     148118 +
-  30 1148-001        : 0000  C        1  0/19:42:15 -   0/19:43:45     1  148119     150467 +
-  31 SUN             : 0000           1  0/19:44:35 -   0/19:53:05     1  150468     168022 +
-  32 1148-001        : 0000  C        1  0/20:07:05 -   0/20:08:35     1  168023     170342 +
-  33 SUN             : 0000           1  0/20:09:25 -   0/20:13:35     1  170343     178889 +
-  34 1148-001        : 0000  C        1  0/20:32:05 -   0/20:33:35     1  178890     181422 +
-  35 SUN             : 0000           1  0/20:34:25 -   0/20:38:35     1  181423     190047 +
-  36 1148-001        : 0000  C        1  0/20:57:05 -   0/20:58:35     1  190048     192547 +
-  37 SUN             : 0000           1  0/20:59:25 -   0/21:03:35     1  192548     201271 +
-  38 1148-001        : 0000  C        1  0/21:21:55 -   0/21:23:25     1  201272     203709 +
-  39 SUN             : 0000           1  0/21:24:15 -   0/21:28:25     1  203710     212708 +
-  40 1148-001        : 0000  C        1  0/21:47:05 -   0/21:48:25     1  212709     215193 +
-  41 SUN             : 0000           1  0/21:49:15 -   0/21:53:25     1  215194     224112 +
-  42 1148-001        : 0000  C        1  0/22:11:45 -   0/22:13:15     1  224113     226323 +
-  43 SUN             : 0000           1  0/22:14:05 -   0/22:18:15     1  226324     234634 +
-  44 1148-001        : 0000  C        1  0/22:36:55 -   0/22:38:15     1  234635     236916 +
-  45 SUN             : 0000           1  0/22:39:05 -   0/22:43:15     1  236917     245740 +
-  46 1148-001        : 0000  C        1  0/23:01:55 -   0/23:03:15     1  245741     248122 +
-  47 SUN             : 0000           1  0/23:04:05 -   0/23:08:05     1  248123     256112 +
-  48 1148-001        : 0000  C        1  0/23:26:45 -   0/23:28:05     1  256113     258141 +
-  49 SUN             : 0000           1  0/23:28:55 -   0/23:33:05     1  258142     265941 +
-  50 1148-001        : 0000  C        1  0/23:51:45 -   0/23:53:05     1  265942     267662 +
-  51 SUN             : 0000           1  0/23:54:05 -   0/23:58:05     1  267663     275388 +
-  52 SUN             : 0000           1  1/00:07:25 -   1/00:07:55     1  275389     276069 +
-  53 3C286           : 0000  B        1  1/00:21:35 -   1/00:22:55     1  276070     278465 +
-  +
-Source summary +
-Velocity type = '        '    Definition = '        ' +
-  +
-  ID Source           Qual  Calcode RA(   0.0)     Dec(   0.0)   IFlux   QFlux   UFlux   VFlux  No. vis +
-   1 1148-001        : 0000   C     11:48:10.1300 -00:07:13.300   0.000   0.000   0.000   0.000   59936 +
-   2 SUN             : 0000         12:10:30.1272 -01:08:22.648   0.000   0.000   0.000   0.000  216133 +
-   3 3C286           : 0000   B     13:28:49.6570  30:45:58.640   0.000   0.000   0.000   0.000    2396 +
-  +
-  ID Source            Freq(GHz) Velocity(Km/s) Rest freq (GHz) +
-   1 All Sources          1.4462         0.0000          0.0000 +
-  +
-Frequency Table summary +
-FQID IF#      Freq(GHz)      BW(kHz)   Ch.Sep(kHz)  Sideband +
-   1   1       1.44615001   12500.0010  12500.0010      1 +
-AIPS 1: Resumes+
 </file> </file>
 +
 +===== LISTR --- scan listing =====
  
  
 '' ''
-getn 4\\+task 'listr'\\ 
 +indi 1\\ 
 +getn 1\\ 
 +optype 'scan'\\ 
 +docrt=132\\ 
 +flagver 0\\ 
 +sources ' '\\ 
 +stokes ' '\\ 
 +docalib 0\\ 
 +gainuse 1\\ 
 +dopol -1\\ 
 +dparm 0\\
 go\\ go\\
 '' ''
  
 <file> <file>
-AIPS 1: Got(1)   disk= 1  user= 100   type=UV   19810926.L BAND.2 + vlb054    LISTR(31DEC10)    100     21-OCT-2010  17:17:44    Page    1 
- vlb054    LISTR(31DEC08)    100     11-NOV-2008  18:49:17    Page    1 +File = 19810926    .L BAND.   1 Vol = 1  Userid =  100 
-File = 19810926    .L BAND.   2 Vol = 1  Userid =  100 +Freq =  1.446150006 GHz   Ncor =  4   No. vis =    421247
-Freq =  1.446150006 GHz   Ncor =  4   No. vis =    207072+
 Scan summary listing Scan summary listing
    
 Scan      Source      Qual  Calcode Sub         Timerange          FrqID   START VIS  END VIS Scan      Source      Qual  Calcode Sub         Timerange          FrqID   START VIS  END VIS
-   1 SUN             : 0000           1  0/14:05:15 -   0/14:09:35     1       1       8935 +   1 1148-001        : 0000  C        1  0/13:53:25 -   0/13:53:25     1       1          0 
-   2 SUN             : 0000           1  0/14:30:05 -   0/14:34:35     1    8936      17885 +   2 1148-001        : 0000  C        1  0/13:53:35 -   0/13:54:35     1       1       1113 
-   3 SUN             : 0000           1  0/14:55:05 -   0/14:59:25     1   17886      26458 +   3 SUN             : 0000           1  0/13:59:15 -   0/13:59:35     1    1114       1575 
-   4 SUN             : 0000           1  0/15:19:55 -   0/15:24:25     1   26459      35408 +   4 SUN             : 0000           1  0/14:05:25 -   0/14:09:35     1    1576       5952 
-   5 SUN             : 0000           1  0/15:44:55 -   0/15:49:25     1   35409      44108 +   5 1148-001        : 0000  C        1  0/14:18:25 -   0/14:19:35     1    5953       7738 
-   6 SUN             : 0000           1  0/16:09:45 -   0/16:14:15     1   44109      52834 +   6 SUN             : 0000           1  0/14:20:25 -   0/14:24:35     1    7739      12262 
-   7 SUN             : 0000           1  0/16:34:45 -   0/16:39:15     1   52835      61669 +   7 SUN             : 0000           1  0/14:30:25 -   0/14:34:35     1   12263      17489 
-   8 SUN             : 0000           1  0/16:59:45 -   0/17:04:05     1   61670      70226 +   8 1148-001        : 0000  C        1  0/14:43:25 -   0/14:44:35     1   17490      19470 
-   9 SUN             : 0000           1  0/17:24:35 -   0/17:29:05     1   70227      79180 +   9 SUN             : 0000           1  0/14:45:25 -   0/14:49:35     1   19471      25116 
-  10 SUN             : 0000           1  0/17:50:15 -   0/17:54:05     1   79181      86877 +  10 SUN             : 0000           1  0/14:55:25 -   0/14:59:25     1   25117      30808 
-  11 SUN             : 0000           1  0/18:14:25 -   0/18:18:55     1   86878      95785 +  11 1148-001        : 0000  C        1  0/15:08:15 -   0/15:09:25     1   30809      32950 
-  12 SUN             : 0000           1  0/18:39:35 -   0/18:44:05     1   95786     104689 +  12 SUN             : 0000           1  0/15:10:15 -   0/15:14:25     1   32951      39133 
-  13 SUN             : 0000           1  0/19:04:25 -   0/19:08:45     1  104690     112638 +  13 SUN             : 0000           1  0/15:20:15 -   0/15:24:25     1   39134      45646 
-  14 SUN             : 0000           1  0/19:29:15 -   0/19:33:45     1  112639     120854 +  14 1148-001        : 0000  C        1  0/15:33:15 -   0/15:34:25     1   45647      47809 
-  15 SUN             : 0000           1  0/19:54:25 -   0/19:58:45     1  120855     128861 +  15 SUN             : 0000           1  0/15:35:15 -   0/15:39:25     1   47810      54604 
-  16 SUN             : 0000           1  0/20:19:05 -   0/20:23:35     1  128862     137749 +  16 SUN             : 0000           1  0/15:45:15 -   0/15:49:25     1   54605      61507 
-  17 SUN             : 0000           1  0/20:44:05 -   0/20:48:35     1  137750     146600 +  17 1148-001        : 0000  C        1  0/15:58:15 -   0/15:59:25     1   61508      63889 
-  18 SUN             : 0000           1  0/21:08:55 -   0/21:13:25     1  146601     155488 +  18 SUN             : 0000           1  0/16:00:15 -   0/16:04:15     1   63890      70412 
-  19 SUN             : 0000           1  0/21:33:55 -   0/21:38:25     1  155489     164306 +  19 SUN             : 0000           1  0/16:10:05 -   0/16:14:15     1   70413      77315 
-  20 SUN             : 0000           1  0/21:58:55 -   0/22:03:25     1  164307     173005 +  20 1148-001        : 0000  C        1  0/16:23:05 -   0/16:24:15     1   77316      79905 
-  21 SUN             : 0000           1  0/22:23:55 -   0/22:28:15     1  173006     181479 +  21 SUN             : 0000           1  0/16:25:05 -   0/16:29:15     1   79906      88655 
-  22 SUN             : 0000           1  0/22:48:55 -   0/22:53:25     1  181480     190299 +  22 SUN             : 0000           1  0/16:35:05 -   0/16:39:15     1   88656      97490 
-  23 SUN             : 0000           1  0/23:13:45 -   0/23:18:05     1  190300     199102 +  23 1148-001        : 0000  C        1  0/16:48:05 -   0/16:49:15     1   97491     100105 
-  24 SUN             : 0000           1  0/23:38:45 -   0/23:43:05     1  199103     207072 +  24 SUN             : 0000           1  0/16:50:05 -   0/16:54:15     1  100106     108820 
-  25 SUN             : 0000           1  1/00:01:15 -   1/00:01:15     1  207073     207072+  25 SUN             : 0000           1  0/17:00:05 -   0/17:04:05     1  108821     117377 
 +  26 1148-001        : 0000  C        1  0/17:12:55 -   0/17:14:05     1  117378     119783 
 +  27 SUN             : 0000           1  0/17:14:55 -   0/17:19:05     1  119784     128650 
 +  28 SUN             : 0000           1  0/17:24:55 -   0/17:29:05     1  128651     137604 
 +  29 1148-001        : 0000  C        1  0/17:37:55 -   0/17:39:05     1  137605     140174 
 +  30 SUN             : 0000           1  0/17:39:45 -   0/17:44:25     1  140175     149746 
 +  31 SUN             : 0000           1  0/17:50:25 -   0/17:54:05     1  149747     157443 
 +  32 1148-001        : 0000  C        1  0/18:02:45 -   0/18:03:55     1  157444     159899 
 +  33 SUN             : 0000           1  0/18:04:45 -   0/18:08:55     1  159900     168765 
 +  34 SUN             : 0000           1  0/18:14:45 -   0/18:18:55     1  168766     177673 
 +  35 1148-001        : 0000  C        1  0/18:27:45 -   0/18:28:55     1  177674     180284 
 +  36 SUN             : 0000           1  0/18:29:45 -   0/18:33:55     1  180285     189114 
 +  37 SUN             : 0000           1  0/18:39:45 -   0/18:44:05     1  189115     198018 
 +  38 1148-001        : 0000  C        1  0/18:52:45 -   0/18:53:55     1  198019     200001 
 +  39 SUN             : 0000           1  0/18:54:45 -   0/18:59:05     1  200002     207881 
 +  40 SUN             : 0000           1  0/19:04:35 -   0/19:08:45     1  207882     215830 
 +  41 1148-001        : 0000  C        1  0/19:17:35 -   0/19:18:45     1  215831     218097 
 +  42 SUN             : 0000           1  0/19:19:35 -   0/19:23:45     1  218098     226243 
 +  43 SUN             : 0000           1  0/19:29:35 -   0/19:33:45     1  226244     234459 
 +  44 1148-001        : 0000  C        1  0/19:42:35 -   0/19:43:45     1  234460     236808 
 +  45 SUN             : 0000           1  0/19:44:35 -   0/19:53:05     1  236809     254363 
 +  46 SUN             : 0000           1  0/19:54:35 -   0/19:58:45     1  254364     262370 
 +  47 1148-001        : 0000  C        1  0/20:07:25 -   0/20:08:35     1  262371     264690 
 +  48 SUN             : 0000           1  0/20:09:25 -   0/20:13:35     1  264691     273237 
 +  49 SUN             : 0000           1  0/20:19:25 -   0/20:23:35     1  273238     282125 
 +  50 1148-001        : 0000  C        1  0/20:32:25 -   0/20:33:35     1  282126     284658 
 +  51 SUN             : 0000           1  0/20:34:25 -   0/20:38:35     1  284659     293283 
 +  52 SUN             : 0000           1  0/20:44:25 -   0/20:48:35     1  293284     302134 
 +  53 1148-001        : 0000  C        1  0/20:57:25 -   0/20:58:35     1  302135     304634 
 +  54 SUN             : 0000           1  0/20:59:25 -   0/21:03:35     1  304635     313358 
 +  55 SUN             : 0000           1  0/21:09:15 -   0/21:13:25     1  313359     322246 
 +  56 1148-001        : 0000  C        1  0/21:22:15 -   0/21:23:25     1  322247     324222 
 +  57 SUN             : 0000           1  0/21:24:15 -   0/21:28:25     1  324223     333221 
 +  58 SUN             : 0000           1  0/21:34:15 -   0/21:38:25     1  333222     341392 
 +  59 1148-001        : 0000  C        1  0/21:47:15 -   0/21:48:25     1  341393     343406 
 +  60 SUN             : 0000           1  0/21:49:15 -   0/21:53:25     1  343407     350399 
 +  61 SUN             : 0000           1  0/21:59:15 -   0/22:03:25     1  350400     357289 
 +  62 1148-001        : 0000  C        1  0/22:12:05 -   0/22:13:15     1  357290     359048 
 +  63 SUN             : 0000           1  0/22:14:05 -   0/22:18:15     1  359049     365542 
 +  64 SUN             : 0000           1  0/22:24:05 -   0/22:28:15     1  365543     372227 
 +  65 1148-001        : 0000  C        1  0/22:37:05 -   0/22:38:15     1  372228     373947 
 +  66 SUN             : 0000           1  0/22:39:05 -   0/22:43:15     1  373948     380876 
 +  67 SUN             : 0000           1  0/22:49:05 -   0/22:53:25     1  380877     387324 
 +  68 1148-001        : 0000  C        1  0/23:02:05 -   0/23:03:15     1  387325     388927 
 +  69 SUN             : 0000           1  0/23:04:05 -   0/23:08:05     1  388928     394642 
 +  70 SUN             : 0000           1  0/23:13:55 -   0/23:18:05     1  394643     401050 
 +  71 1148-001        : 0000  C        1  0/23:26:55 -   0/23:28:05     1  401051     402397 
 +  72 SUN             : 0000           1  0/23:28:55 -   0/23:33:05     1  402398     407527 
 +  73 SUN             : 0000           1  0/23:38:55 -   0/23:43:05     1  407528     412709 
 +  74 1148-001        : 0000  C        1  0/23:51:55 -   0/23:53:05     1  412710     413612 
 +  75 SUN             : 0000           1  0/23:54:05 -   0/23:58:05     1  413613     418616 
 +  76 SUN             : 0000           1  1/00:07:35 -   1/00:07:55     1  418617     419028 
 +  77 3C286           : 0000  B        1  1/00:21:45 -   1/00:22:55     1  419029     421247
    
 Source summary Source summary
Line 228: Line 256:
    
   ID Source           Qual  Calcode RA(   0.0)     Dec(   0.0)   IFlux   QFlux   UFlux   VFlux  No. vis   ID Source           Qual  Calcode RA(   0.0)     Dec(   0.0)   IFlux   QFlux   UFlux   VFlux  No. vis
-   1 SUN             : 0000         12:10:30.1272 -01:08:22.648   0.000   0.000   0.000   0.000  207072+   1 1148-001        : 0000   C     11:48:10.1300 -00:07:13.300   0.000   0.000   0.000   0.000   52089 
 +   2 SUN             : 0000         12:10:30.1272 -01:08:22.648   0.000   0.000   0.000   0.000  366939 
 +   3 3C286           : 0000   B     13:28:49.6570  30:45:58.640   0.000   0.000   0.000   0.000    2219
    
   ID Source            Freq(GHz) Velocity(Km/s) Rest freq (GHz)   ID Source            Freq(GHz) Velocity(Km/s) Rest freq (GHz)
Line 245: Line 275:
  
 '' ''
-getn 2\\ 
 go prtan\\ go prtan\\
 '' ''
  
 <file> <file>
- vlb054    PRTAN(31DEC08)    100     11-NOV-2008  18:53:09    Page    1+ vlb054    PRTAN(31DEC08)    100     14-NOV-2008  21:45:10    Page    1
 File=19810926    .L BAND.   1     An.ver=   1     Vol= 1     User=  100 File=19810926    .L BAND.   1     An.ver=   1     Vol= 1     User=  100
 Array= VLA          Freq=  1446.150006 MHz     Ref.date= 26-SEP-1981 Array= VLA          Freq=  1446.150006 MHz     Ref.date= 26-SEP-1981
Line 402: Line 431:
 AIPS 1: Resumes AIPS 1: Resumes
 </file> </file>
 +
 +The really important part for you, the data reducer, is the antenna layout at the bottom.  This tells you in a nice graphical form where each antenna is located.  This helps you to figure out which baselines are small, and which baselines are long.  This is also useful for figuring out other things related to antenna position and baseline direction.
 +
 +
 +
 +===== UVPLT =====
 +
 +Let's have a look at the information about the visibilities.
 +
 +''
 +task 'uvplt'\\
 +source '3c286',' '\\
 +stokes 'half'\\
 +calcode ' '\\
 +uvrange 0\\
 +antennas 0\\
 +basel 0\\
 +xinc 1\\
 +aparm 0\\
 +bparm 0\\
 +doweight 1\\
 +refant 3\\
 +do3col 1\\
 +dotv 1\\
 +tvinit\\
 +go\\
 +''
 +
 +
 +{{:public:meetings:glowbi2008:fig1.jpg|}}
 +
 +
 +
 +''
 +source '1148-001', ' '\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig2.jpg|}}
 +
 +
 +
 +Warning: the next step will take a long time.  You are free to skip this, and just look at the picture.
 +
 +''
 +source 'sun',' '\\
 +go\\
 +''
 +
 +{{:public:meetings:glowbi2008:fig3.jpg|}}
 +
 +
 +''
 +source '3c286',' '\\
 +bparm 0 2\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig4.jpg|}}
 +
 +
 +
 +''
 +source '1148-001', ' '\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig5.jpg|}}
 +
 +
 +
 +''
 +source '3c286',' '\\
 +bparm 11 1\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig6.jpg|}}
 +
 +
 +
 +
 +''
 +source '1148-001', ' '\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig7.jpg|}}
 +
 +
 +
 +''
 +source '3c286',' '\\
 +bparm 6 7\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig8.jpg|}}
 +
 +
 +
 +''
 +source '1148-001', ' '\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig9.jpg|}}
 +
 +
 +''
 +source ' '\\
 +bparm 11 15\\
 +xinc 50\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig10.jpg|}}
 +
 +
 +''
 +bparm 11 18\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig11.jpg|}}
 +
 +
 +
 +===== LISTR =====
 +
 +Now let's look at the raw visibility numbers.  
 +
 +
 +
 +''
 +xinc 1\\
 +task 'listr'\\
 +optype 'list'\\
 +source '3c286', ' '\\
 +stokes ' '\\
 +bif 0\\
 +eif 0\\
 +flagver 0\\
 +dparm 0\\
 +antenna 3,0\\
 +go\\
 +''
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +<file>
 + vlb054    LISTR(31DEC08)    100     14-NOV-2008  22:05:45    Page    1
 +File = 19810926    .L BAND.   1 Vol = 1  Userid =  100   Channels =  1-   1 IF = 1
 +Freq=  1.446150006 GHz   Ncor=  1   No. vis=      3159
 +Stokes = RR   Subarray =   1
 + 
 + 
 +Source=3C286           : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = B   , Freq =  1.446150006 GHz
 +Amplitudes, 1000 =   1.000 Jy, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   1/00:21:35  293 259 258 282  45 360 431 333     563 326 247 264 256 469 287 127 283  84  56 455  12 344 307 313 321 347
 +Amplitudes, 1000 =  10.000 Jy, averging type = Vector
 +   1/00:21:45  195 180 182 168  15 212 234 203     247 135 160 177 169 187 174 156 178 129  66 205   1 185 186 189 158 158
 +   1/00:21:55  206 196 194 176  67 223 246 214     252 139 169 192 183 193 183 185 190 188 205 214   1 192 199 199 165 168
 +   1/00:22:05  209 197 197 179 165 223 246 217     251 144 172 193 185 194 186 186 193 190 235 216   1 193 201 202 167 169
 +   1/00:22:15  205 194 196 178 198 223 244 214     246 140 170 189 182 193 183 186 190 190 235 212   0 193 198 199 164 169
 +   1/00:22:25  206 192 197 177 196 222 243 214     242 139 165 189 180 193 183 184 190 186 231 212   1 191 197 199 164 168
 +   1/00:22:35  207 195 196 178 198 223 246 215     240 140 166 192 182 192 184 185 190 190 233 214   0 192 199 199 165 168
 +   1/00:22:45  206 192 192 177 199 222 243 213     232 141 168 191 183 190 183 185 191 187 233 214   1 191 197 198 164 169
 +   1/00:22:55  206 193 195 179 199 223 243 213     230 140 170 190 182 193 183 185 190 186 233 215   1 191 197 199 164 169
 +AIPS 1: Resumes
 +</file>
 +
 +
 +
 +
 +
 +===== Flagging =====
 +
 +
 +Just enter these commands.  We can talk later about why they
 +needed to be applied.
 +
 +
 +''
 +task 'uvflg'\\
 +antenna 23,0\\
 +basel 0\\
 +timer 0\\
 +aparm 0\\
 +opcode 'flag'\\
 +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 'listr'\\
 +optype 'list'\\
 +source '1148-001', ' '\\
 +stokes ' '\\
 +calcode ' '\\
 +bif 0\\
 +eif 0\\
 +timer 0\\
 +flagver 0\\
 +dparm 0\\
 +antenna 3,0\\
 +go\\
 +''
 +
 +
 +
 +<file>
 + vlb054    LISTR(31DEC08)    100     14-NOV-2008  22:35:27    Page    1
 +File = 19810926    .L BAND.   1 Vol = 1  Userid =  100   Channels =  1-   1 IF = 1
 +Freq=  1.446150006 GHz   Ncor=  1   No. vis=     84072
 +Stokes = RR   Subarray =   1
 +Applying flag table   1
 + 
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Amplitudes, 1000 =  10.000 Jy, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/13:53:55   57 126  58  53      86 123  73     195  32  57  13  23  66  30  36  27  75  30  89      24  20  28  45 103
 +   0/13:54:05   67 134  66  24      65 155  47     213  26  66  28  28  92  30  53  35  60  26  63      28  33  55  77 103
 +   0/13:54:15   43 124  77  29  43  26 144  44     184  36  77  40  29  96  19  39  29  34  31  34      26  34  76  62 134
 +   0/13:54:25    6 104  91  62  40  68  96  73     177  59  93  30  30  87  35   4  26  41  45  78      32  31  74  14 139
 +   0/13:54:35   28  84  94  70  21  93  69  80     191  56 104  19  31  61  34  50  31  68  42  94      31  24  41  46  99
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Amplitudes, 1000 =  10.000 Jy, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/14:18:25   37  65  58  31     101 107  83     230  68  57  35  34  81  33  90  40  56  18  32      32  45  69  32  47
 +   0/14:18:35    4  65  55  22      78  85  80     191  67  58  28  28  78  26  58  28  58  33  50      29  39  28  17  39
 +   0/14:18:45   36  62  56  30  17  32  45  47     128  48  57  28  32  61  32  48  24  31  44  48      34  26  18  29  17
 +   0/14:18:55   54  54  70  40  24  39  15  18     194  27  61  30  36  28  33  80  34   6  45  27      34  25  52  36  35
 +   0/14:19:05   48  27  87  54  30  90  70  24     176  21  56  29  47  14  39 114  42  29  50  42      43  44  63  24  47
 +   0/14:19:15   38   7 105  40  33 117 116  67     113  59  63  42  43  51  44 105  45  70  51  62      44  58  66  29  39
 +   0/14:19:25   23  48 112  35  30  75 117  92     198  80  50  46  44  84  46  57  38  74  59  63      47  48  21  50  30
 +   0/14:19:35   58  72  99  36  26  39  70  66     123  60  31  38  40  86  34  35  41  42  49  31      42  27  33  46  35
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Amplitudes, 1000 =  10.000 Jy, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/14:43:25   32  32  95  40      47  27  55     237  14  83  30  33  58  31  33  38  40  20  35      34  35  13  28  35
 +   0/14:43:35   34  56  91  30      44  48  60     173  29  80  34  33  46  32  29  36  38  35  37      34  33  38  29  37
 +   0/14:43:45   27  59 103  14  27  21  51  54     130  51  66  25  25  72  24  23  22  14  27  26      21  23  48  23  26
 +   0/14:43:55   22  49 112  29  24  18  35  23     182  49  43  19  26  95  23  30  20   9  30  28      27  22  32  20  21
 +   0/14:44:05   39  59 141  41  28  42  36   2     163  35  32  32  33 112  31  32  37  35  41  35      34  35  29  30  30
 +   0/14:44:15   36  49 157  37  23  50  28  29     179  16  17  31  32  98  32  31  38  44  36  34      29  33  18  28  35
 +   0/14:44:25   27  28 159  25  22  37  38  54     200  17  17  30  29  71  28  26  27  38  33  31      27  28  36  28  31
 +   0/14:44:35   28   1 150  20  31  23  47  60     242  39  22  29  30  42  28  29  27  22  39  34      30  29  51  28  27
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Amplitudes, 1000 =  10.000 Jy, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/15:08:15   17  35  52  25      33  42  35     148  14  29  28  29  49  29  31  30  27  15  38      29  32  35  25  33
 +   0/15:08:25   29  38  64  22      27  33  32     146  14  38  27  25  52  25  31  33  18  29  30      27  22  23  22  35
 +   0/15:08:35   33  25  70  27  23  25  27  25     135  25  35  21  25  49  21  29  26  17  30  25      24  19  18  18  29
 +   0/15:08:45   50  35  91  34  35  42  44  43     185  32  52  37  34  49  34  34  30  38  45  40      36  39  43  33  33
 +   0/15:08:55   30  28  84  30  32  37  43  34     188  26  53  32  33  26  32  26  29  42  37  40      32  36  41  27  26
 +   0/15:09:05   21  27  72  28  31  37  46  35     185  26  56  31  32   9  31  27  37  36  37  39      32  30  35  29  35
 +</file>
 +
 +
 +''
 +dparm 1,0\\
 +go\\
 +''
 +
 +
 +<file>
 + vlb054    LISTR(31DEC08)    100     14-NOV-2008  22:36:24    Page    1
 +File = 19810926    .L BAND.   1 Vol = 1  Userid =  100   Channels =  1-   1 IF = 1
 +Freq=  1.446150006 GHz   Ncor=  1   No. vis=     84072
 +Stokes = RR   Subarray =   1
 +Applying flag table   1
 + 
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Phase, 1000 = 1000.00 degrees, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/13:53:55 -154 176 177-145      82  59 -83     -35-118 -77 101  41   3  59  63 114  69 -97  80      76  88 -57  67 121
 +   0/13:54:05  163 138-130 -83     141 119 -21     -20 -30 -50  52  43  52  78  97 122 122-106 147      78  88  46 166-137
 +   0/13:54:15  127 100 -91  65-134-101 172  78    -173  44 -26  79  44  90  51 137 129-158-129 -79      72 102  95-124 -55
 +   0/13:54:25  118  61 -61 125 -91  13-128 145     172 105  -3  99  45 126  51-125 122 -51-119  21      73 111 140   2  12
 +   0/13:54:35 -125   9 -32 168 -49  68 -31-164     161 144  17  64  39 167  75  31 115  12-106  79      77  97-174 142  90
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Phase, 1000 = 1000.00 degrees, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/14:18:25  135 174  69 167      66 121-173      51 109  26  75  39  58  61 136 118  27-117  30      73  88 140-157  -8
 +   0/14:18:35  -72 153 124 160     123 168-128      93 145  32  67  42  93  49-154 126  72-128  56      67 111-158 177  13
 +   0/14:18:45 -131 132 173 146-116-171-153 -91     147-177  42  68  36 129  50 -56 105 107-120  83      73 112  25 160 -43
 +   0/14:18:55 -165 103-132 149-123 -19  -9 -18       6-122  55  73  40 167  61  23 101-175-110  89      76  81  66-175 -41
 +   0/14:19:05  162  87 -88 160-120  46  81 143     -37  33  72  73  37  23  53  78 112  -3-118  46      73  83 112-159 -16
 +   0/14:19:15  141 -86 -66 164-116  89 113 174    -114  69  87  66  40  22  49 120 122  26-126  54      73  98 143 165  -7
 +   0/14:19:25 -137-143 -42 155-114 141 151-151     118 111  94  69  35  50  56 167 112  65-120  78      73 111 155 166 -18
 +   0/14:19:35 -139-164 -20 142-130 -91-161-118     166 144 107  72  37  79  55 -66  99 109-114  88      78  99  37-168 -38
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Phase, 1000 = 1000.00 degrees, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/14:43:25 -174-121 103 162      65 109 166      41  37  70  66  32-169  49  86  97  42-121  64      64  91  89-178 -29
 +   0/14:43:35 -169-148 150 173      81 104-159     138  70  88  65  32 -85  48  76 109  66-129  60      67  92  71 178 -20
 +   0/14:43:45 -165-171-168 150-127  97 120-123     116 107 114  64  28 -20  54  69 109  87-132  61      56  94  90 170 -18
 +   0/14:43:55 -170 163-127 133-123  58 147 -84      90 142 139  68  35  32  53  74  91   6-133  61      58  85 121 175 -25
 +   0/14:44:05 -177 139 -92 148-112  53 149 -32      72 177 164  69  35  70  47  82  97  19-127  64      64  90 134-174 -30
 +   0/14:44:15  180 123 -62 162-122  68 116 147      62-116-143  65  32 105  54  81 109  35-124  64      66  90  87 175 -30
 +   0/14:44:25 -176 101 -37 168-136  85 105 177      51  15 -85  66  30 145  50  76 112  59-132  57      64  88  72 171 -20
 +   0/14:44:35 -162 105  -7 146-133  79 117-152      41  68 -46  68  30-163  47  71  96  67-132  61      59  87  97 173 -24
 + 
 +Source=1148-001        : 0000, Stokes=RR  , IF=  1, Chans=   1-   1
 +Flux =  0.0000 Jy, Calcode = C   , Freq =  1.446150006 GHz
 +Phase, 1000 = 1000.00 degrees, averging type = Vector
 + 
 +Baselines      1 3 2 3 3 4 3 5 3 6 3 7 3 8 3 9 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329
 +   0/15:08:15 -152 172-105 148      73 125-161     -32 127 -46  69  28  39  48  65  87  38-129  62      62  93 108 179 -40
 +   0/15:08:25 -136 152 -86 139      76 126-157     -27  89 -36  66  28  63  53  80  95  41-139  66      61  96 114 175 -30
 +   0/15:08:35 -164 140 -63 149-130  64 135-148     -15 107 -11  67  22  87  50  78 114   9-134  60      67  80  93 177 -17
 +   0/15:08:45  178 144 -39 151-126  66 128-148      -3 119   9  68  28 106  51  90 108  12-127  58      65  82  91-177 -16
 +   0/15:08:55  167 150 -18 151-124  70 124-149       9 132  23  69  27 121  48  70  87  27-127  57      62  87 101-179 -33
 +   0/15:09:05 -160 162   3 151-129  72 126-157      21 135  33  65  28  74  51  48  93  40-134  67      61  96 111 177 -36
 +</file>
 +
 +====== Calibration ======
 +
 +
 +===== SETJY =====
 +
 +''
 +task 'setjy'\\
 +source '3c286',' '\\
 +optype 'calc'\\
 +aparm 0\\
 +aparm(2)=4\\
 +zerosp 0\\
 +go\\
 +''
 +
 +<file>
 +vlb054> SETJY1: Task SETJY  (release of 31DEC10) begins
 +vlb054> SETJY1: **WARNING: OPCODE=CALC AND FREQID =   -1
 +vlb054> SETJY1:      FREQID WILL BE RESET TO 1, CHECK YOUR RESULTS CAREFULLY
 +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: '3C286           ' IF =  1 FLUX =14.6428 (Jy calcd)
 +vlb054> SETJY1: / Using (1995.2) VLA or Reynolds (1934-638) coefficients
 +vlb054> SETJY1: Appears to have ended successfully
 +vlb054> SETJY1: vlb054       31DEC10 TST: Cpu=       0.0  Real=       0
 +</file>
 +
 +===== CALRD =====
 +
 +''
 +task 'calrd'\\
 +object '3c286'\\
 +band 'l'\\
 +go\\
 +''
 +
 +
 +<file>
 +vlb054> CALRD1: Task CALRD  (release of 31DEC08) begins
 +vlb054> CALRD1: Reading disk file AIPSTARS:3C286_L.MODEL
 +vlb054> CALRD1: Create 3C286_L     .MODEL .   1 (MA)  on disk  1  cno    2
 +vlb054> CALRD1: Appears to have ended successfully
 +vlb054> CALRD1: vlb054       31DEC08 TST: Cpu=       0.0  Real=       0
 +</file>
 +
 +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\\
 +''
 +
 +{{:public:meetings:glowbi2008:fig14.jpg|}}
 +
 +
 +===== UVFIX =====
 +
 +The sky positions (RA and Dec) were entered as coordinates of date, because the Sun is a moving object.  However, they were entered in an odd fashion which means that the current software has trouble dealing with the calibration image for proper calibration.  Run the fix program.
 +
 +
 +''
 +task 'uvfix'\\
 +getn 1\\
 +clron\\
 +shift 0\\
 +uvfixprm 0\\
 +go\\
 +''
 +
 +UVFIX change in April means we have to regenerate NX table\\
 +
 +''
 +task 'indxr'\\
 +getn 3\\
 +cparm 0\\
 +cparm(3)= -1\\
 +bparm 0\\
 +go\\
 +''
 +
 +
 +
 +===== CALIB =====
 +
 +First, we do the amplitude calibrator.
 +
 +
 +''
 +task 'calib'\\
 +getn 3\\
 +calsour '3c286',' '\\
 +uvrange 0\\
 +antennas 0\\
 +timer 0\\
 +refant 24\\
 +weightit 1\\
 +in2di 1\\
 +get2n 2\\
 +ncomp 0\\
 +solmode 'A&P'\\
 +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 '1148-001',' '\\
 +clr2n\\
 +uvrange 1.5,0\\
 +solint 0\\
 +go\\
 +clrmsg\\
 +''
 +
 +===== SNPLT =====
 +
 +Look at the amplitude calibrator results.
 +
 +''
 +task 'snplt'\\
 +source ' '\\
 +inver 1\\
 +inext 'sn'\\
 +stokes ' '\\
 +optype 'phas'\\
 +opcode 'alsi'\\
 +xinc 1\\
 +nplots 9\\
 +antenna 0\\
 +timer 0\\
 +tvinit\\
 +dotv 1\\
 +go\\
 +optype 'amp'\\
 +go\\
 +''
 +
 +Look at the phase calibrator results.
 +
 +
 +
 +''
 +inver 2\\
 +optype 'phas'\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig12.jpg|}}
 +
 +
 +
 +
 +''
 +optype 'amp'\\
 +go\\
 +''
 +
 +{{:public:meetings:glowbi2008:fig13.jpg|}}
 +
 +
 +===== GETJY =====
 +
 +Now we need to transfer the amplitude calibration information from the amplitude calibrator to the phase calibrator.
 +
 +''
 +task 'getjy'\\
 +sources '1148-001',' '\\
 +calsour '3c286',' '\\
 +calcode ' '\\
 +bif 0\\
 +eif 0\\
 +antenna 0\\
 +timer 0\\
 +snver 0\\
 +go\\
 +''
 +
 +
 +<file>
 +vlb054> GETJY1: Task GETJY  (release of 31DEC08) begins
 +vlb054> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy)
 +vlb054> GETJY1: 1148-001        :  0   C     1    2.90731 +/-   0.04192
 +vlb054> GETJY1: Appears to have ended successfully
 +vlb054> GETJY1: vlb054       31DEC08 TST: Cpu=       0.0  Real=       0
 +</file>
 +
 +2.9 Jy agrees relatively well with the rough value given by the VLA calibrator list.
 +
 +<file>
 +1150-003   J2000  A 11h50m43.870761s -00d23'54.204900"  Aug01         
 +1148-001   B1950  A 11h48m10.124900s -00d07'13.164000"
 +-----------------------------------------------------
 +BAND        A B C D    FLUX(Jy)    UVMIN(kL)  UVMAX(kL)
 +=====================================================
 + 20cm    L  P P P P       2.80                       visplot
 +  6cm    C  P P P P       1.92
 +3.7cm    X  P P P P       1.25                       visplot
 +  2cm    U  P P P P       1.40
 +1.3cm    K  S S S S       0.63                       visplot
 +0.7cm    Q  W W W W       0.65''
 +</file>
 +
 +===== Redo the CALIB calibration for the phase calibrator =====
 +
 +
 +''
 +inext 'sn'\\
 +inver 2\\
 +extdest\\
 +tget calib\\
 +go\\
 +clrmsg\\
 +''
 +
 +===== Apply the calibration to all sources, both amplitude and phase =====
 +
 +
 +''
 +task 'clcal'\\
 +getn 3\\
 +sources ' '\\
 +calsour '3c286','1148-001',' '\\
 +calcode ' '\\
 +opcode 'cali'\\
 +gainver 1\\
 +gainuse 2\\
 +refant 24\\
 +bparm 0\\
 +interpol '2pt'\\
 +timer 0\\
 +antenna 0\\
 +doblank 1\\
 +go\\
 +''
 +
 +
 +===== Check the calibration results =====
 +
 +''
 +task 'snplt'\\
 +source ' '\\
 +inver 2\\
 +inext 'cl'\\
 +stokes ' '\\
 +optype 'phas'\\
 +opcode 'alsi'\\
 +xinc 1\\
 +nplots 9\\
 +antenna 0\\
 +timer 0\\
 +tvinit\\
 +dotv 1\\
 +go\\
 +''
 +
 +{{:public:meetings:glowbi2008:fig15.jpg|}}
 +
 +
 +''
 +optype 'amp'\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig16.jpg|}}
 +
 +
 +
 +====== Test Imaging ======
 +
 +===== Amplitude Calibrator =====
 +
 +
 +Ok, let's start with a few test images.  First, we will image the amplitude calibrator.  Because of the funny epoch conversion stuff, it is not quite at the center where it is expected.  This is not a problem.
 +
 +
 +
 +''
 +task 'imagr'\\
 +getn 3\\
 +source '3c286',' '\\
 +uvrange 0\\
 +docalib 1\\
 +gainuse 2\\
 +stokes 'i'\\
 +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.
 +
 +{{:public:meetings:glowbi2008:fig17.jpg|}}
 +
 +
 +And this is the point spread function of the actual observations (called the dirty beam).  It looks like the amplitude calibrator image before you started cleaning.  Why?
 +
 +
 +''
 +getn 5\\
 +tvinit\\
 +tvlod\\
 +tvfiddle\\
 +''
 +
 +{{:public:meetings:glowbi2008:fig18.jpg|}}
 +
 +===== Phase Calibrator =====
 +
 +Ok, clean up after ourselves
 +
 +
 +''
 +getn 5;zap\\
 +getn 6;zap\\
 +''
 +
 +
 +Now do the imaging and light cleaning
 +
 +
 +
 +''
 +tget imagr\\
 +source '1148-001',' '\\
 +rashift 0\\
 +decshift 0\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig19.jpg|}}
 +
 +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\\
 +''
 +{{:public:meetings:glowbi2008:fig2.jpg|}}
 +
 +====== Imaging the Sun ======
 +
 +
 +Ok, to do the best job, we need to separate the Sun data from the rest of the dataset.  First, cleanup after the phase calibrator imaging.
 +
 +''
 +getn 5;zap\\
 +getn 6;zap\\
 +''
 +
 +
 +
 +===== SPLIT =====
 +
 +
 +
 +''
 +task 'split'\\
 +getn 3\\
 +source 'sun',' '\\
 +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.
 +
 +{{:public:meetings:glowbi2008:fig21.jpg|}}
 +{{:public:meetings:glowbi2008:fig22.jpg|}}
 +{{:public:meetings:glowbi2008:fig23.jpg|}}
 +
 +
 +
 +===== Self-Calibration =====
 +
 +Ok, let's do some self-calibration.  This is just the same as "normal" calibration, except that you are saying to yourself that you are less confident that you know what the source is supposed to look like.
 +
 +We run exactly the same calibration software, it exactly the same manner.  Because the phase calibrator was resolved at short spacings, and because it is quite difficult for us to clean the extended emission, I limit the //(u,v)// coverage.  Note that this also is only doing phase calibration, not amplitude calibration.
 +
 +
 +''
 +tget calib\\
 +getn 4\\
 +get2n 7\\
 +solint 30/60\\
 +solmode 'p!a'\\
 +docalib 0\\
 +uvrange 0.5,0\\
 +calsour ' '\\
 +inver 0\\
 +snver 0\\
 +go\\
 +clrmsg\\
 +''
 +
 +==== Check the results ====
 +
 +
 +
 +''
 +task 'snplt'\\
 +inext 'sn'\\
 +inver 1\\
 +stokes ' '\\
 +optype 'phas'\\
 +tvinit\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig24.jpg|}}
 +
 +
 +Note some of the really large phase deviations!  The phase calibration was not perfect.
 +
 +===== Image again, using the new calibration =====
 +
 +
 +
 +
 +
 +''
 +getn 5;zap\\
 +getn 6;zap\\
 +getn 7;zap\\
 +tget imagr\\
 +docalib 1\\
 +gainuse 1\\
 +go\\
 +''
 +{{:public:meetings:glowbi2008:fig25.jpg|}}
 +
 +This one is a bit better. 
 +
 +
 +
 +===== More Self-Calibration =====
 +But we can do more.  Now let's do both amplitude and phase calibration.  First, we do a phase-only self-calibration on short timescales, to get rid of the phase jitter.
 +
 +==== Phase Self-Cal ====
 +
 +
 +
 +''
 +tget calib\\
 +getn 4\\
 +get2n 7\\
 +solint 20/60\\
 +solmode 'p!a'\\
 +docalib 0\\
 +uvrange 0.5,0\\
 +calsour ' '\\
 +inver 0\\
 +snver 0\\
 +go\\
 +clrmsg\\
 +''
 +==== Check ====
 +
 +
 +
 +''
 +task 'snplt'\\
 +inext 'sn'\\
 +inver 2\\
 +stokes ' '\\
 +optype 'phas'\\
 +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 'A&P'\\
 +cparm(2) 1\\
 +go\\
 +clrmsg\\
 +''
 +
 +
 +==== Check ====
 +
 +
 +''
 +task 'snplt'\\
 +inext 'sn'\\
 +inver 1\\
 +stokes ' '\\
 +optype 'phas'\\
 +tvinit\\
 +go\\
 +optype 'amp'\\
 +go\\
 +''
 +
 +
 +==== Cleanup ====
 +
 +Clean up the files.  But save the initial "CALIB" uv dataset!  It is the one with the current amplitude calibration information.
 +
 +
 +''
 +getn 6;zap\\
 +getn 7;zap\\
 +getn 8;zap\\
 +''
 +
 +
 +===== Final Imaging =====
 +
 +
 +
 +
 +
 +''
 +tget imagr\\
 +getn 5\\
 +docalib 1\\
 +gainuse 1\\
 +''
 +
 +{{:public:meetings:glowbi2008:fig26.jpg|}}
 +
 +
 +
 +
 +====== Discussion ======
 +
 +What have you learned?
 +
 +Does your image look like this image?
 +{{http://images.nrao.edu/images/SUN_VR16C_lo.jpg}}
 +
 +What is different?
 +
 +Have a look at http://images.nrao.edu/8 and see what additional information was used to generate the NRAO web image.  How would such additional information change your image?
  
  
  
  • Last modified: 2008-11-11 18:20
  • by James M Anderson