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:47] – 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 439: Line 613:
  
  
 +
 +
 +
 +===== 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\\
 +''
  
  
Line 502: Line 726:
 </file> </file>
  
 +
 +''
 +dparm 1,0\\
 +go\\
 +''
  
  
Line 563: Line 792:
    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    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> </file>
 +
 +====== Calibration ======
 +
  
 ===== SETJY ===== ===== SETJY =====
Line 570: 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 607: Line 840:
 </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:47
  • by James M Anderson