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-14 21:08] – update 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 79: Line 98:
  
 '' ''
 +indisk 1\\
 pcat pcat
 '' ''
Line 90: Line 110:
  
 '' ''
-getn 2\\+getn 1\\
 imhe imhe
 '' ''
Line 132: Line 152:
  
 '' ''
-default listr\\ 
 task 'listr'\\ task 'listr'\\
 indi 1\\ indi 1\\
-getn 2\\+getn 1\\
 optype 'scan'\\ optype 'scan'\\
 docrt=132\\ docrt=132\\
 +flagver 0\\
 +sources ' '\\
 +stokes ' '\\
 +docalib 0\\
 +gainuse 1\\
 +dopol -1\\
 +dparm 0\\
 go\\ go\\
 '' ''
  
 <file> <file>
- vlb054    LISTR(31DEC08)    100     14-NOV-2008  21:42:46    Page    1+ vlb054    LISTR(31DEC10)    100     21-OCT-2010  17:17:44    Page    1
 File = 19810926    .L BAND.   1 Vol = 1  Userid =  100 File = 19810926    .L BAND.   1 Vol = 1  Userid =  100
-Freq =  1.446150006 GHz   Ncor =  4   No. vis =    556884+Freq =  1.446150006 GHz   Ncor =  4   No. vis =    421247
 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 1148-001        : 0000  C        1  0/13:53:25 -   0/13:53:25     1       1        351 +   1 1148-001        : 0000  C        1  0/13:53:25 -   0/13:53:25     1       1          0 
-   2 1148-001        : 0000  C        1  0/13:53:35 -   0/13:54:35     1     352       2808 +   2 1148-001        : 0000  C        1  0/13:53:35 -   0/13:54:35     1       1       1113 
-   3 SUN             : 0000           1  0/13:59:05 -   0/13:59:35     1    2809       4212 +   3 SUN             : 0000           1  0/13:59:15 -   0/13:59:35     1    1114       1575 
-   4 SUN             : 0000           1  0/14:05:15 -   0/14:09:35     1    4213      13689 +   4 SUN             : 0000           1  0/14:05:25 -   0/14:09:35     1    1576       5952 
-   5 1148-001        : 0000  C        1  0/14:18:05 -   0/14:19:35     1   13690      17199 +   5 1148-001        : 0000  C        1  0/14:18:25 -   0/14:19:35     1    5953       7738 
-   6 SUN             : 0000           1  0/14:20:15 -   0/14:24:35     1   17200      26676 +   6 SUN             : 0000           1  0/14:20:25 -   0/14:24:35     1    7739      12262 
-   7 SUN             : 0000           1  0/14:30:05 -   0/14:34:35     1   26677      36504 +   7 SUN             : 0000           1  0/14:30:25 -   0/14:34:35     1   12263      17489 
-   8 1148-001        : 0000  C        1  0/14:43:05 -   0/14:44:35     1   36505      40014 +   8 1148-001        : 0000  C        1  0/14:43:25 -   0/14:44:35     1   17490      19470 
-   9 SUN             : 0000           1  0/14:45:15 -   0/14:49:35     1   40015      49491 +   9 SUN             : 0000           1  0/14:45:25 -   0/14:49:35     1   19471      25116 
-  10 SUN             : 0000           1  0/14:55:05 -   0/14:59:25     1   49492      58968 +  10 SUN             : 0000           1  0/14:55:25 -   0/14:59:25     1   25117      30808 
-  11 1148-001        : 0000  C        1  0/15:07:55 -   0/15:09:25     1   58969      62478 +  11 1148-001        : 0000  C        1  0/15:08:15 -   0/15:09:25     1   30809      32950 
-  12 SUN             : 0000           1  0/15:10:05 -   0/15:14:25     1   62479      71955 +  12 SUN             : 0000           1  0/15:10:15 -   0/15:14:25     1   32951      39133 
-  13 SUN             : 0000           1  0/15:19:55 -   0/15:24:25     1   71956      81783 +  13 SUN             : 0000           1  0/15:20:15 -   0/15:24:25     1   39134      45646 
-  14 1148-001        : 0000  C        1  0/15:32:55 -   0/15:34:25     1   81784      85293 +  14 1148-001        : 0000  C        1  0/15:33:15 -   0/15:34:25     1   45647      47809 
-  15 SUN             : 0000           1  0/15:35:05 -   0/15:39:25     1   85294      94770 +  15 SUN             : 0000           1  0/15:35:15 -   0/15:39:25     1   47810      54604 
-  16 SUN             : 0000           1  0/15:44:55 -   0/15:49:25     1   94771     104598 +  16 SUN             : 0000           1  0/15:45:15 -   0/15:49:25     1   54605      61507 
-  17 1148-001        : 0000  C        1  0/15:57:55 -   0/15:59:25     1  104599     108108 +  17 1148-001        : 0000  C        1  0/15:58:15 -   0/15:59:25     1   61508      63889 
-  18 SUN             : 0000           1  0/16:00:05 -   0/16:04:15     1  108109     117234 +  18 SUN             : 0000           1  0/16:00:15 -   0/16:04:15     1   63890      70412 
-  19 SUN             : 0000           1  0/16:09:45 -   0/16:14:15     1  117235     127062 +  19 SUN             : 0000           1  0/16:10:05 -   0/16:14:15     1   70413      77315 
-  20 1148-001        : 0000  C        1  0/16:22:45 -   0/16:24:15     1  127063     130572 +  20 1148-001        : 0000  C        1  0/16:23:05 -   0/16:24:15     1   77316      79905 
-  21 SUN             : 0000           1  0/16:24:55 -   0/16:29:15     1  130573     140049 +  21 SUN             : 0000           1  0/16:25:05 -   0/16:29:15     1   79906      88655 
-  22 SUN             : 0000           1  0/16:34:45 -   0/16:39:15     1  140050     149877 +  22 SUN             : 0000           1  0/16:35:05 -   0/16:39:15     1   88656      97490 
-  23 1148-001        : 0000  C        1  0/16:47:45 -   0/16:49:15     1  149878     153387 +  23 1148-001        : 0000  C        1  0/16:48:05 -   0/16:49:15     1   97491     100105 
-  24 SUN             : 0000           1  0/16:49:55 -   0/16:54:15     1  153388     162864 +  24 SUN             : 0000           1  0/16:50:05 -   0/16:54:15     1  100106     108820 
-  25 SUN             : 0000           1  0/16:59:45 -   0/17:04:05     1  162865     172341 +  25 SUN             : 0000           1  0/17:00:05 -   0/17:04:05     1  108821     117377 
-  26 1148-001        : 0000  C        1  0/17:12:35 -   0/17:14:05     1  172342     175851 +  26 1148-001        : 0000  C        1  0/17:12:55 -   0/17:14:05     1  117378     119783 
-  27 SUN             : 0000           1  0/17:14:45 -   0/17:19:05     1  175852     185328 +  27 SUN             : 0000           1  0/17:14:55 -   0/17:19:05     1  119784     128650 
-  28 SUN             : 0000           1  0/17:24:35 -   0/17:29:05     1  185329     195156 +  28 SUN             : 0000           1  0/17:24:55 -   0/17:29:05     1  128651     137604 
-  29 1148-001        : 0000  C        1  0/17:37:35 -   0/17:39:05     1  195157     198666 +  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  198667     208845 +  30 SUN             : 0000           1  0/17:39:45 -   0/17:44:25     1  140175     149746 
-  31 SUN             : 0000           1  0/17:50:15 -   0/17:54:05     1  208846     217269 +  31 SUN             : 0000           1  0/17:50:25 -   0/17:54:05     1  149747     157443 
-  32 1148-001        : 0000  C        1  0/18:02:25 -   0/18:03:55     1  217270     220779 +  32 1148-001        : 0000  C        1  0/18:02:45 -   0/18:03:55     1  157444     159899 
-  33 SUN             : 0000           1  0/18:04:35 -   0/18:08:55     1  220780     230256 +  33 SUN             : 0000           1  0/18:04:45 -   0/18:08:55     1  159900     168765 
-  34 SUN             : 0000           1  0/18:14:25 -   0/18:18:55     1  230257     240084 +  34 SUN             : 0000           1  0/18:14:45 -   0/18:18:55     1  168766     177673 
-  35 1148-001        : 0000  C        1  0/18:27:25 -   0/18:28:55     1  240085     243594 +  35 1148-001        : 0000  C        1  0/18:27:45 -   0/18:28:55     1  177674     180284 
-  36 SUN             : 0000           1  0/18:29:35 -   0/18:33:55     1  243595     253071 +  36 SUN             : 0000           1  0/18:29:45 -   0/18:33:55     1  180285     189114 
-  37 SUN             : 0000           1  0/18:39:35 -   0/18:44:05     1  253072     262899 +  37 SUN             : 0000           1  0/18:39:45 -   0/18:44:05     1  189115     198018 
-  38 1148-001        : 0000  C        1  0/18:52:35 -   0/18:53:55     1  262900     266058 +  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  266059     275535 +  39 SUN             : 0000           1  0/18:54:45 -   0/18:59:05     1  200002     207881 
-  40 SUN             : 0000           1  0/19:04:25 -   0/19:08:45     1  275536     285012 +  40 SUN             : 0000           1  0/19:04:35 -   0/19:08:45     1  207882     215830 
-  41 1148-001        : 0000  C        1  0/19:17:25 -   0/19:18:45     1  285013     288171 +  41 1148-001        : 0000  C        1  0/19:17:35 -   0/19:18:45     1  215831     218097 
-  42 SUN             : 0000           1  0/19:19:25 -   0/19:23:45     1  288172     297648 +  42 SUN             : 0000           1  0/19:19:35 -   0/19:23:45     1  218098     226243 
-  43 SUN             : 0000           1  0/19:29:15 -   0/19:33:45     1  297649     307476 +  43 SUN             : 0000           1  0/19:29:35 -   0/19:33:45     1  226244     234459 
-  44 1148-001        : 0000  C        1  0/19:42:15 -   0/19:43:45     1  307477     310986 +  44 1148-001        : 0000  C        1  0/19:42:35 -   0/19:43:45     1  234460     236808 
-  45 SUN             : 0000           1  0/19:44:25 -   0/19:53:05     1  310987     329589 +  45 SUN             : 0000           1  0/19:44:35 -   0/19:53:05     1  236809     254363 
-  46 SUN             : 0000           1  0/19:54:25 -   0/19:58:45     1  329590     339066 +  46 SUN             : 0000           1  0/19:54:35 -   0/19:58:45     1  254364     262370 
-  47 1148-001        : 0000  C        1  0/20:07:05 -   0/20:08:35     1  339067     342576 +  47 1148-001        : 0000  C        1  0/20:07:25 -   0/20:08:35     1  262371     264690 
-  48 SUN             : 0000           1  0/20:09:15 -   0/20:13:35     1  342577     352053 +  48 SUN             : 0000           1  0/20:09:25 -   0/20:13:35     1  264691     273237 
-  49 SUN             : 0000           1  0/20:19:05 -   0/20:23:35     1  352054     361881 +  49 SUN             : 0000           1  0/20:19:25 -   0/20:23:35     1  273238     282125 
-  50 1148-001        : 0000  C        1  0/20:32:05 -   0/20:33:35     1  361882     365391 +  50 1148-001        : 0000  C        1  0/20:32:25 -   0/20:33:35     1  282126     284658 
-  51 SUN             : 0000           1  0/20:34:15 -   0/20:38:35     1  365392     374868 +  51 SUN             : 0000           1  0/20:34:25 -   0/20:38:35     1  284659     293283 
-  52 SUN             : 0000           1  0/20:44:05 -   0/20:48:35     1  374869     384696 +  52 SUN             : 0000           1  0/20:44:25 -   0/20:48:35     1  293284     302134 
-  53 1148-001        : 0000  C        1  0/20:57:05 -   0/20:58:35     1  384697     388206 +  53 1148-001        : 0000  C        1  0/20:57:25 -   0/20:58:35     1  302135     304634 
-  54 SUN             : 0000           1  0/20:59:15 -   0/21:03:35     1  388207     397657 +  54 SUN             : 0000           1  0/20:59:25 -   0/21:03:35     1  304635     313358 
-  55 SUN             : 0000           1  0/21:08:55 -   0/21:13:25     1  397658     407485 +  55 SUN             : 0000           1  0/21:09:15 -   0/21:13:25     1  313359     322246 
-  56 1148-001        : 0000  C        1  0/21:21:55 -   0/21:23:25     1  407486     410969 +  56 1148-001        : 0000  C        1  0/21:22:15 -   0/21:23:25     1  322247     324222 
-  57 SUN             : 0000           1  0/21:24:05 -   0/21:28:25     1  410970     420446 +  57 SUN             : 0000           1  0/21:24:15 -   0/21:28:25     1  324223     333221 
-  58 SUN             : 0000           1  0/21:33:55 -   0/21:38:25     1  420447     430274 +  58 SUN             : 0000           1  0/21:34:15 -   0/21:38:25     1  333222     341392 
-  59 1148-001        : 0000  C        1  0/21:47:05 -   0/21:48:25     1  430275     433433 +  59 1148-001        : 0000  C        1  0/21:47:15 -   0/21:48:25     1  341393     343406 
-  60 SUN             : 0000           1  0/21:49:05 -   0/21:53:25     1  433434     442858 +  60 SUN             : 0000           1  0/21:49:15 -   0/21:53:25     1  343407     350399 
-  61 SUN             : 0000           1  0/21:58:55 -   0/22:03:25     1  442859     452686 +  61 SUN             : 0000           1  0/21:59:15 -   0/22:03:25     1  350400     357289 
-  62 1148-001        : 0000  C        1  0/22:11:45 -   0/22:13:15     1  452687     456196 +  62 1148-001        : 0000  C        1  0/22:12:05 -   0/22:13:15     1  357290     359048 
-  63 SUN             : 0000           1  0/22:13:55 -   0/22:18:15     1  456197     465647 +  63 SUN             : 0000           1  0/22:14:05 -   0/22:18:15     1  359049     365542 
-  64 SUN             : 0000           1  0/22:23:55 -   0/22:28:15     1  465648     475124 +  64 SUN             : 0000           1  0/22:24:05 -   0/22:28:15     1  365543     372227 
-  65 1148-001        : 0000  C        1  0/22:36:55 -   0/22:38:15     1  475125     478283 +  65 1148-001        : 0000  C        1  0/22:37:05 -   0/22:38:15     1  372228     373947 
-  66 SUN             : 0000           1  0/22:38:55 -   0/22:43:15     1  478284     487734 +  66 SUN             : 0000           1  0/22:39:05 -   0/22:43:15     1  373948     380876 
-  67 SUN             : 0000           1  0/22:48:55 -   0/22:53:25     1  487735     497562 +  67 SUN             : 0000           1  0/22:49:05 -   0/22:53:25     1  380877     387324 
-  68 1148-001        : 0000  C        1  0/23:01:55 -   0/23:03:15     1  497563     500721 +  68 1148-001        : 0000  C        1  0/23:02:05 -   0/23:03:15     1  387325     388927 
-  69 SUN             : 0000           1  0/23:03:55 -   0/23:08:05     1  500722     509847 +  69 SUN             : 0000           1  0/23:04:05 -   0/23:08:05     1  388928     394642 
-  70 SUN             : 0000           1  0/23:13:45 -   0/23:18:05     1  509848     519324 +  70 SUN             : 0000           1  0/23:13:55 -   0/23:18:05     1  394643     401050 
-  71 1148-001        : 0000  C        1  0/23:26:45 -   0/23:28:05     1  519325     522249 +  71 1148-001        : 0000  C        1  0/23:26:55 -   0/23:28:05     1  401051     402397 
-  72 SUN             : 0000           1  0/23:28:45 -   0/23:33:05     1  522250     531700 +  72 SUN             : 0000           1  0/23:28:55 -   0/23:33:05     1  402398     407527 
-  73 SUN             : 0000           1  0/23:38:45 -   0/23:43:05     1  531701     540450 +  73 SUN             : 0000           1  0/23:38:55 -   0/23:43:05     1  407528     412709 
-  74 1148-001        : 0000  C        1  0/23:51:45 -   0/23:53:05     1  540451     543350 +  74 1148-001        : 0000  C        1  0/23:51:55 -   0/23:53:05     1  412710     413612 
-  75 SUN             : 0000           1  0/23:53:45 -   0/23:58:05     1  543351     552125 +  75 SUN             : 0000           1  0/23:54:05 -   0/23:58:05     1  413613     418616 
-  76 SUN             : 0000           1  1/00:01:15 -   1/00:01:15     1  552126     552450 +  76 SUN             : 0000           1  1/00:07:35 -   1/00:07:55     1  418617     419028 
-  77 SUN             : 0000           1  1/00:07:25 -   1/00:07:55     1  552451     553725 +  77 3C286           : 0000  B        1  1/00:21:45 -   1/00:22:55     1  419029     421247
-  78 3C286           : 0000  B        1  1/00:21:35 -   1/00:22:55     1  553726     556884+
    
 Source summary Source summary
Line 231: 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 1148-001        : 0000   C     11:48:10.1300 -00:07:13.300   0.000   0.000   0.000   0.000   84072 +   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  469653 +   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    3159+   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 408: Line 433:
  
 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. 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\\
 +''
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
 +
  
  
Line 437: Line 611:
 AIPS 1: Resumes AIPS 1: Resumes
 </file> </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 ===== ===== SETJY =====
Line 444: Line 802:
 source '3c286',' '\\ source '3c286',' '\\
 optype 'calc'\\ optype 'calc'\\
-aparm(2)=3\\+aparm 0\\ 
 +aparm(2)=4\\ 
 +zerosp 0\\
 go\\ go\\
 '' ''
- 
  
 <file> <file>
-vlb054> SETJY1: Task SETJY  (release of 31DEC08) begins+vlb054> SETJY1: Task SETJY  (release of 31DEC10) begins
 vlb054> SETJY1: **WARNING: OPCODE=CALC AND FREQID =   -1 vlb054> SETJY1: **WARNING: OPCODE=CALC AND FREQID =   -1
 vlb054> SETJY1:      FREQID WILL BE RESET TO 1, CHECK YOUR RESULTS CAREFULLY vlb054> SETJY1:      FREQID WILL BE RESET TO 1, CHECK YOUR RESULTS CAREFULLY
-vlb054> SETJY1: A source model for this calibrator may be available +vlb054> SETJY1: A source model for this calibrator is available 
-vlb054> SETJY1: Use the verb CALDIR to see if there is one+vlb054> SETJY1: Consult the help file for CALRD for assistance
 vlb054> SETJY1: / Flux calculated using known spectrum vlb054> SETJY1: / Flux calculated using known spectrum
 vlb054> SETJY1: BIF =  1 EIF =  1 /Range of IFs vlb054> SETJY1: BIF =  1 EIF =  1 /Range of IFs
-vlb054> SETJY1: '3C286           ' IF =  1 FLUX =14.5979 (Jy calcd) +vlb054> SETJY1: '3C286           ' IF =  1 FLUX =14.6428 (Jy calcd) 
-vlb054> SETJY1: / Using (1990) VLA or Reynolds (1934-638) coefficients+vlb054> SETJY1: / Using (1995.2) VLA or Reynolds (1934-638) coefficients
 vlb054> SETJY1: Appears to have ended successfully vlb054> SETJY1: Appears to have ended successfully
-vlb054> SETJY1: vlb054       31DEC08 TST: Cpu=       0.0  Real=       0+vlb054> SETJY1: vlb054       31DEC10 TST: Cpu=       0.0  Real=       0
 </file> </file>
  
Line 476: Line 835:
 vlb054> CALRD1: Task CALRD  (release of 31DEC08) begins vlb054> CALRD1: Task CALRD  (release of 31DEC08) begins
 vlb054> CALRD1: Reading disk file AIPSTARS:3C286_L.MODEL vlb054> CALRD1: Reading disk file AIPSTARS:3C286_L.MODEL
-vlb054> CALRD1: Create 3C286_L     .MODEL .   1 (MA)  on disk  1  cno    5+vlb054> CALRD1: Create 3C286_L     .MODEL .   1 (MA)  on disk  1  cno    2
 vlb054> CALRD1: Appears to have ended successfully vlb054> CALRD1: Appears to have ended successfully
 vlb054> CALRD1: vlb054       31DEC08 TST: Cpu=       0.0  Real=       0 vlb054> CALRD1: vlb054       31DEC08 TST: Cpu=       0.0  Real=       0
 </file> </file>
  
-So it looks like the calibration data image was put into catalog entry number 5.+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 ===== ===== CALIB =====
 +
 +First, we do the amplitude calibrator.
 +
  
 '' ''
 task 'calib'\\ task 'calib'\\
-getn 2\\+getn 3\\
 calsour '3c286',' '\\ calsour '3c286',' '\\
 uvrange 0\\ uvrange 0\\
 antennas 0\\ antennas 0\\
 +timer 0\\
 refant 24\\ refant 24\\
 weightit 1\\ weightit 1\\
 in2di 1\\ in2di 1\\
-get2n 5\\+get2n 2\\
 ncomp 0\\ ncomp 0\\
 solmode 'A&P'\\ solmode 'A&P'\\
 +aparm 0\\
 aparm(6) 2\\ aparm(6) 2\\
 minamper 10\\ minamper 10\\
 minphser 10\\ minphser 10\\
 docalib 1\\ docalib 1\\
-solint 1\\+gainuse 1\\ 
 +solint 30\\
 solsub 2\\ solsub 2\\
 +cparm 0\\
 cparm(3) 10\\ cparm(3) 10\\
 cparm(4) 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-14 21:08
  • by James M Anderson