public:meetings:glowbi2008:tutorial_wednesday

Differences

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

Link to this comparison view

public:meetings:glowbi2008:tutorial_wednesday [2010-11-16 17:13] – James M Andersonpublic:meetings:glowbi2008:tutorial_wednesday [2017-03-08 15:27] (current) – external edit 127.0.0.1
Line 8: Line 8:
 ===== Web Image of 3C129 ===== ===== Web Image of 3C129 =====
  
-{{http://rsd-www.nrl.navy.mil/7210/7213/backgrnd/3c129_pub.jpg}}+{{https://www.nrl.navy.mil/rsd/sites/www.nrl.navy.mil.rsd/files/Images/radiogalaxy.jpg}}
  
 Wide-field 90 cm image of the large head-tail radio galaxy 3C129 (right). The bright source at the left is 3C129.1. From Lane et al., AJ, 123, 2985, 2002. Wide-field 90 cm image of the large head-tail radio galaxy 3C129 (right). The bright source at the left is 3C129.1. From Lane et al., AJ, 123, 2985, 2002.
Line 553: Line 553:
 ===== Flag some bad data.  Just run the commands. ===== ===== Flag some bad data.  Just run the commands. =====
  
 +
 +''
 +task 'quack'\\
 +reason 'slow start'\\
 +sources ' '\\
 +calcode ' '\\
 +stokes ' '\\
 +timer 0\\
 +antenna 17 0\\
 +flagver 1\\
 +opcode 'beg'\\
 +basel 0\\
 +bif 1\\
 +eif 1\\
 +aparm 0\\
 +aparm(2) 72.5/60\\
 +go\\
 +\\
 +\\
 +reason 'quack'\\
 +sources ' '\\
 +calcode ' '\\
 +stokes ' '\\
 +timer 0\\
 +antenna 0\\
 +flagver 1\\
 +opcode 'beg'\\
 +basel 0\\
 +bif 2\\
 +eif 2\\
 +aparm 0\\
 +aparm(2) 22.5/60\\
 +go\\
 +''
  
 '' ''
Line 558: Line 592:
 default\\ default\\
 dohist 1\\ dohist 1\\
-getn 4\\ +getn 3\\ 
-antenna 1,0\\+antenna 14,0\\
 basel 0\\ basel 0\\
-timer 0 13 45 20 0 15 0 0\\+timer 0 15 10 50 0 15 55 0\\
 stokes ' '\\ stokes ' '\\
 +bif 2\\
 +eif 2\\
 opcode 'flag'\\ opcode 'flag'\\
 reason 'bad ant'\\ reason 'bad ant'\\
 outfgver 1\\ outfgver 1\\
 go\\ go\\
-antenna 14,0\\ +timer 0 13 25 30 0 13 27 40\\
-timer 0 16 50 50 0 16 51 0\\+
 go\\ go\\
-timer 0+timer 0 13 59 0 0 14 2 20\\ 
 +go\\ 
 +timer 0\\ 
 +bif 0\\ 
 +eif 0\\ 
 +antenna 0\\
 '' ''
- 
  
 ===== Start the calibration ===== ===== Start the calibration =====
Line 582: Line 621:
 task 'setjy'\\ task 'setjy'\\
 default\\ default\\
-getn 4\\+getn 3\\
 source '0518+165',' '\\ source '0518+165',' '\\
 bif 0\\ bif 0\\
Line 600: Line 639:
 vlb054> SETJY1: / Flux calculated using known spectrum vlb054> SETJY1: / Flux calculated using known spectrum
 vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs
-vlb054> SETJY1: '0518+165        ' IF =  1 FLUX = 3.7742 (Jy calcd) +vlb054> SETJY1: '0518+165        ' IF =  1 FLUX = 2.4983 (Jy calcd) 
-vlb054> SETJY1: '0518+165        ' IF =  2 FLUX = 3.9464 (Jy calcd)+vlb054> SETJY1: '0518+165        ' IF =  2 FLUX = 2.4166 (Jy calcd)
 vlb054> SETJY1: / Using (1995.2) 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
Line 611: Line 650:
 eif 1\\ eif 1\\
 optype ' '\\ optype ' '\\
-zerosp 3.7742, 0.380, -0.153, 0.0\\+zerosp 2.4983, 0.3092, -0.0628, 0.0\\
 go\\ go\\
 bif 2\\ bif 2\\
 eif 2\\ eif 2\\
-zerosp 3.9464, 0.397, -0.161, 0.0\\+zerosp 2.4166, 0.2990, -0.0607, 0.0\\
 go\\ go\\
 '' ''
Line 631: Line 670:
  
 <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: A source model for this calibrator may be available +
-vlb054> SETJY1: Use the verb CALDIR to see if there is one+
 vlb054> SETJY1: / Flux calculated using known spectrum vlb054> SETJY1: / Flux calculated using known spectrum
 vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs
-vlb054> SETJY1: '0134+329        ' IF =  1 FLUX = 5.4891 (Jy calcd) +vlb054> SETJY1: '0134+329        ' IF =  1 FLUX = 3.2045 (Jy calcd) 
-vlb054> SETJY1: '0134+329        ' IF =  2 FLUX = 5.8249 (Jy calcd)+vlb054> SETJY1: '0134+329        ' IF =  2 FLUX = 3.0714 (Jy calcd)
 vlb054> SETJY1: / Using (1995.2) 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>
 +
 +''
 +bif 1\\
 +eif 1\\
 +optype ' '\\
 +zerosp 3.2045, 0.0772, -0.1722, 0.0\\
 +go\\
 +bif 2\\
 +eif 2\\
 +zerosp 3.0714, 0.0740, -0.1650, 0.0\\
 +go\\
 +\\
 +bif 0\\
 +eif 0\\
 +''
 +
  
  
Line 710: Line 763:
  
 '' ''
-task 'indxr'+task 'indxr'\\
 getn 12\\ getn 12\\
 cparm 0\\ cparm 0\\
Line 782: Line 835:
 do3col 1\\ do3col 1\\
 dotv 1\\ dotv 1\\
 +xinc 1\\
 nplots 9\\ nplots 9\\
 tvinit\\ tvinit\\
Line 819: Line 873:
 <file> <file>
 vlb054> GETJY1: Task GETJY  (release of 31DEC10) begins vlb054> GETJY1: Task GETJY  (release of 31DEC10) begins
-vlb054> GETJY1: Calibrator robust averaging used    200 of    200 gain samples+vlb054> GETJY1: Calibrator robust averaging used    800 of    800 gain samples
 vlb054> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy) vlb054> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy)
-vlb054> GETJY1: 0420+417        :  3   B     1    1.46712 +/-   0.00294  0.00294 +vlb054> GETJY1: 0420+417        :  2   B     1    1.41467 +/-   0.00432  0.00431 
-vlb054> GETJY1:                              2    1.47148 +/-   0.00234  0.00344+vlb054> GETJY1:                              2    1.40641 +/-   0.00455  0.00455
 vlb054> GETJY1:    Source:Qual    CALCODE    used   total     bad used  tot  bad vlb054> GETJY1:    Source:Qual    CALCODE    used   total     bad used  tot  bad
-vlb054> GETJY1: 0420+417        :  3 B        900     900       0   97  100    0+vlb054> GETJY1: 0420+417        :  2 B        904     904       0  100  100    0
 vlb054> GETJY1: Appears to have ended successfully vlb054> GETJY1: Appears to have ended successfully
-vlb054> GETJY1: vlb054       31DEC10 TST: Cpu=       0.0  Real=       0+vlb054> GETJY1: vlb054 31DEC10 TST: Cpu=      0.0  Real=      0  IO=         1
 </file> </file>
  
Line 864: Line 918:
 ===== Have a look at the calibration ===== ===== Have a look at the calibration =====
  
-==== Plots ====+
  
 ==== Plots ==== ==== Plots ====
Line 982: Line 1036:
 <file> <file>
 vlb054> PCAL 1: Calibration source   1   0420+417 vlb054> PCAL 1: Calibration source   1   0420+417
-vlb054> PCAL 1:   Q+iU=(  -0.01034,   0.03759) +/- (  0.000106,  0.000106) Jy +vlb054> PCAL 1:   Q+iU=(  -0.02115,  -0.06575) +/- (  0.000088,  0.000088) Jy 
-vlb054> PCAL 1:   Pol. inten. =   0.03899 +/-  0.000150 Jy +vlb054> PCAL 1:   Pol. inten. =   0.06906 +/-  0.000124 Jy 
-vlb054> PCAL 1:   Pol. angle  =  52.69 +/-  0.078 deg +vlb054> PCAL 1:   Pol. angle  = -53.92 +/-  0.036 deg 
-vlb054> PCAL 1: 0420+417         I,Q,U,V=   1.4671 -0.01034  0.03759  0.00000 Jy+vlb054> PCAL 1: 0420+417         I,Q,U,V=   1.4147 -0.02115 -0.06575  0.00000 Jy
 vlb054> PCAL 1: Calibration source   1   0420+417 vlb054> PCAL 1: Calibration source   1   0420+417
-vlb054> PCAL 1:   Q+iU=(  -0.03422,   0.00814) +/- (  0.000113,  0.000113) Jy +vlb054> PCAL 1:   Q+iU=(   0.03641,   0.05888) +/- (  0.000116,  0.000116) Jy 
-vlb054> PCAL 1:   Pol. inten. =   0.03518 +/-  0.000160 Jy +vlb054> PCAL 1:   Pol. inten. =   0.06922 +/-  0.000163 Jy 
-vlb054> PCAL 1:   Pol. angle  =  83.31 +/-  0.092 deg +vlb054> PCAL 1:   Pol. angle  =  29.14 +/-  0.048 deg 
-vlb054> PCAL 1: 0420+417         I,Q,U,V=   1.4715 -0.03422  0.00814  0.00000 Jy+vlb054> PCAL 1: 0420+417         I,Q,U,V=   1.4064  0.03641  0.05888  0.00000 Jy
 </file> </file>
  
Line 1002: Line 1056:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0,45\\+uvrange 0\\
 docalib 1\\ docalib 1\\
 gainuse 2\\ gainuse 2\\
Line 1011: Line 1065:
  
 <file> <file>
-Ampscalar average of matrix  =   -30.20(  0.091) sigma =  1.307=  1.160 +Ampscalar average of matrix  =    90.43(  0.020) sigma =  0.496 
-Ampscalar average of matrix  =    33.92(  0.103) sigma =  1.487=  1.389+Ampscalar average of matrix  =  -102.71(  0.065) sigma =  1.583
 </file> </file>
  
Line 1024: Line 1078:
 AIPS 1: Adverbs     Values                 Comments AIPS 1: Adverbs     Values                 Comments
 AIPS 1: ---------------------------------------------------------------- AIPS 1: ----------------------------------------------------------------
-AIPS 1: CLCORPRM    -30.2029     33.9241   R-L phase difference returned+AIPS 1: CLCORPRM     90.4283   -102.7064   R-L phase difference returned
 AIPS 1:            *rest 0 AIPS 1:            *rest 0
 </file> </file>
Line 1044: Line 1098:
 http://www.vla.nrao.edu/astro/calib/polar/1999/ http://www.vla.nrao.edu/astro/calib/polar/1999/
  
-Correct number is -21.71+Correct number is -22.98
  
  
Line 1050: Line 1104:
  
 '' ''
-for i=1 to 2;clcorp(i) = -21.71 - clcorp(i);end\\+for i=1 to 2;clcorp(i) = -22.98 - clcorp(i);end\\
 '' ''
  
Line 1077: Line 1131:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0,45\\+uvrange 0\\
 docalib 1\\ docalib 1\\
 gainuse 3\\ gainuse 3\\
Line 1095: Line 1149:
 AIPS 1: Adverbs     Values                 Comments AIPS 1: Adverbs     Values                 Comments
 AIPS 1: ---------------------------------------------------------------- AIPS 1: ----------------------------------------------------------------
-AIPS 1: CLCORPRM    -21.71      -21.71     R-L phase difference returned+AIPS 1: CLCORPRM    -22.98      -22.98     R-L phase difference returned
 AIPS 1:            *rest 0 AIPS 1:            *rest 0
 </file> </file>
  
  
 +
 +==== Save the tables again ====
 +''
 +task 'tasav'\\
 +clro\\
 +outdi 1\\
 +go\\
 +''
  
  
  
-############### 
  
 ===== Now start with the smaller array data.  Load it in ===== ===== Now start with the smaller array data.  Load it in =====
  
-(Actually, in the 2010 Nov 16 version of these instructions, it is already loaded in slot 7.  You can skip this step if you want.)+(Actually, in the 2010 Nov 16 version of these instructions, it is already loaded in slot 9.  You can skip this step if you want.)
  
 ==== FILLM ==== ==== FILLM ====
Line 1151: Line 1212:
 </file> </file>
  
-==== Makre sure we have the right one ====+==== Make sure we have the right one ====
  
 '' ''
-getn 6\\+getn 9\\
 imhe\\ imhe\\
 '' ''
Line 1166: Line 1227:
 task 'listr'\\ task 'listr'\\
 default\\ default\\
-getn 7\\+getn 9\\
 optype 'scan'\\ optype 'scan'\\
 doweight 1\\ doweight 1\\
Line 1173: Line 1234:
  
 <file> <file>
- vlb054    LISTR(31DEC10)    100     16-NOV-2010  16:10:46    Page    1 + vlb054    LISTR(31DEC10)    100     16-NOV-2010  21:19:54    Page    1 
-File = 19940725    .C BAND.   4 Vol = 1  Userid =  100 +File = 19940725    .X BAND.   4 Vol = 1  Userid =  100 
-Freq =  7.814900000 GHz   Ncor =  4   No. vis =    313007+Freq =  8.514900000 GHz   Ncor =  4   No. vis =    287494
 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 0134+329        : 0001  C        1101/03:50:05 - 101/03:50:05     1       1        351 +   1 0134+329        : 0002  C        1101/04:01:46 - 101/04:03:55     1       1       3635 
-   2 0134+329        : 0001  C        1101/03:50:15 - 101/04:01:35     1     352      24492 +   2 0134+329        : 0002  C        1101/04:04:15 - 101/04:05:05     1    3636       5560 
-   3 0056-001        : 0001  C        1101/04:09:15 - 101/04:11:15     1   24493      24525 +   3 0056-001        : 0002  C        1101/04:11:25 - 101/04:13:45     1    5561       7731 
-   4 0053-016        : 0001           1101/04:14:05 - 101/04:22:05     1   24526      41158 +   4 0053-016        : 0002           1101/04:22:15 - 101/04:30:35     1    7732      23734 
-   5 0053-015        : 0001           1101/04:30:55 - 101/04:38:55     1   41159      57845 +   5 0053-015        : 0002           1101/04:39:04 - 101/04:47:15     1   23735      39425 
-   6 0056-001        : 0001  C        1101/04:47:35 - 101/04:49:36     1   57846      61842 +   6 0056-001        : 0002  C        1101/04:49:45 - 101/04:52:05     1   39426      43741 
-   7 0053-016        : 0001           1101/04:52:25 - 101/05:00:25     1   61843      78492 +   7 0053-016        : 0002           1101/05:00:35 - 101/05:08:55     1   43742      59757 
-   8 0053-015        : 0001           1101/05:09:15 - 101/05:17:15     1   78493      95142 +   8 0053-015        : 0002           1101/05:17:25 - 101/05:25:35     1   59758      74973 
-   9 0056-001        : 0001  C        1101/05:25:55 - 101/05:27:55     1   95143      99139 +   9 0056-001        : 0002  C        1101/05:28:05 - 101/05:30:25     1   74974      79289 
-  10 0053-016        : 0001           1101/05:30:45 - 101/05:38:45     1   99140     115772 +  10 0053-016        : 0002           1101/05:38:55 - 101/05:47:15     1   79290      95250 
-  11 0053-015        : 0001           1101/05:47:35 - 101/05:55:35     1  115773     132440 +  11 0053-015        : 0002           1101/05:55:45 - 101/06:03:45     1   95251     110616 
-  12 0056-001        : 0001  C        1101/06:04:05 - 101/06:06:14     1  132441     136762 +  12 0056-001        : 0002  C        1101/06:06:25 - 101/06:08:45     1  110617     114932 
-  13 0053-016        : 0001           1101/06:09:05 - 101/06:17:05     1  136763     153395 +  13 0053-016        : 0002           1101/06:17:15 - 101/06:25:36     1  114933     130948 
-  14 0053-015        : 0001           1101/06:25:55 - 101/06:33:55     1  153396     170063 +  14 0053-015        : 0002           1101/06:34:05 - 101/06:42:15     1  130949     146625 
-  15 0056-001        : 0001  C        1101/06:42:35 - 101/06:44:35     1  170064     174077 +  15 0056-001        : 0002  C        1101/06:44:45 - 101/06:47:05     1  146626     150928 
-  16 0053-016        : 0001           1101/06:47:25 - 101/06:55:25     1  174078     190710 +  16 0053-016        : 0002           1101/06:55:34 - 101/07:03:55     1  150929     166906 
-  17 0053-015        : 0001           1101/07:04:15 - 101/07:12:15     1  190711     207378 +  17 0053-015        : 0002           1101/07:12:25 - 101/07:20:36     1  166907     182583 
-  18 0056-001        : 0001  C        1101/07:20:55 - 101/07:22:55     1  207379     211375 +  18 0056-001        : 0002  C        1101/07:23:05 - 101/07:25:24     1  182584     186886 
-  19 0420+417        : 0001  B        1101/07:29:35 - 101/07:31:26     1  211376     214487 +  19 0420+417        : 0002  B        1101/07:31:35 - 101/07:33:26     1  186887     190226 
-  20 0518+165        : 0001  C        1101/07:36:55 - 101/07:38:55     1  214488     218426 +  20 0518+165        : 0002  C        1101/07:39:04 - 101/07:42:25     1  190227     196417 
-  21 3C129           : 0001           1101/07:47:35 - 101/08:02:35     1  218427     249671 +  21 3C129           : 0002           1101/08:02:45 - 101/08:17:55     1  196418     225758 
-  22 0420+417        : 0001  B        1101/08:23:55 - 101/08:25:25     1  249672     252804 +  22 0420+417        : 0002  B        1101/08:25:34 - 101/08:27:25     1  225759     228387 
-  23 3C129           : 0001           1101/08:29:46 - 101/08:44:45     1  252805     283300 +  23 0420+417        : 0002  B        1101/08:47:25 - 101/08:49:15     1  228388     231675 
-  24 0420+417        : 0001  B        1101/08:45:45 - 101/08:47:15     1  283301     286784 +  24 3C129           : 0002           1101/08:51:55 - 101/09:06:55     1  231676     260851 
-  25 0420+417        : 0001  B        1101/09:13:26 - 101/09:14:55     1  286785     289484 +  25 0420+417        : 0002  B        1101/09:15:04 - 101/09:16:55     1  260852     264196 
-  26 3C129           : 0001           1101/09:19:45 - 101/09:29:45     1  289485     309848 +  26 3C129           : 0002           1101/09:29:56 - 101/09:40:15     1  264197     284173 
-  27 0420+417        : 0001  B        1101/09:45:55 - 101/09:47:25     1  309849     313007+  27 0420+417        : 0002  B        1101/09:47:34 - 101/09:49:25     1  284174     287494
 Type Q to stop, just hit RETURN to continue                                                                                     Type Q to stop, just hit RETURN to continue                                                                                    
  
- vlb054    LISTR(31DEC10)    100     16-NOV-2010  16:11:06    Page    2 + vlb054    LISTR(31DEC10)    100     16-NOV-2010  21:19:56    Page    2 
-File = 19940725    .C BAND.   4 Vol = 1  Userid =  100+File = 19940725    .X BAND.   4 Vol = 1  Userid =  100
    
 Source summary Source summary
Line 1215: Line 1276:
    
   ID Source           Qual  Calcode RA(1950.0)     Dec(1950.0)   IFlux   QFlux   UFlux   VFlux  No. vis   ID Source           Qual  Calcode RA(1950.0)     Dec(1950.0)   IFlux   QFlux   UFlux   VFlux  No. vis
-   1 0134+329        : 0001   C     01:34:49.8320  32:54:20.520   0.000   0.000   0.000   0.000   24492 +   1 0134+329        : 0002   C     01:34:49.8320  32:54:20.520   0.000   0.000   0.000   0.000    5560 
-   2 0056-001        : 0001   C     00:56:31.7550 -00:09:18.750   0.000   0.000   0.000   0.000   20360 +             IF( 2)                                               0.000   0.000   0.000   0.000 
-   3 0053-016        : 0001         00:53:28.0000 -01:36:45.000   0.000   0.000   0.000   0.000   83182 +   2 0056-001        : 0002   C     00:56:31.7550 -00:09:18.750   0.000   0.000   0.000   0.000   23725 
-   4 0053-015        : 0001         00:53:52.2350 -01:31:57.790   0.000   0.000   0.000   0.000   83341 +             IF( 2)                                               0.000   0.000   0.000   0.000 
-   5 0420+417        : 0001   B     04:20:27.9370  41:43:08.045   0.000   0.000   0.000   0.000   15588 +   3 0053-016        : 0002         00:53:28.0000 -01:36:45.000   0.000   0.000   0.000   0.000   79974 
-   6 0518+165        : 0001   C     05:18:16.5320  16:35:26.900   0.000   0.000   0.000   0.000    3939 +             IF( 2)                                               0.000   0.000   0.000   0.000 
-   7 3C129           : 0001         04:45:31.6950  44:55:19.950   0.000   0.000   0.000   0.000   82105+   4 0053-015        : 0002         00:53:52.2350 -01:31:57.790   0.000   0.000   0.000   0.000   77627 
 +             IF( 2)                                               0.000   0.000   0.000   0.000 
 +   5 0420+417        : 0002   B     04:20:27.9370  41:43:08.045   0.000   0.000   0.000   0.000   15923 
 +             IF( 2)                                               0.000   0.000   0.000   0.000 
 +   6 0518+165        : 0002   C     05:18:16.5320  16:35:26.900   0.000   0.000   0.000   0.000    6191 
 +             IF( 2)                                               0.000   0.000   0.000   0.000 
 +   7 3C129           : 0002         04:45:31.6950  44:55:19.950   0.000   0.000   0.000   0.000   78494 
 +             IF( 2)                                               0.000   0.000   0.000   0.000
    
   ID Source            Freq(GHz) Velocity(Km/s) Rest freq (GHz)   ID Source            Freq(GHz) Velocity(Km/s) Rest freq (GHz)
-   1 All Sources          7.8149         0.0000          0.0000+   1 All Sources          8.5149         0.0000          0.0000 
 +     IF(  2)              8.8851         0.0000          0.0000
    
 Frequency Table summary Frequency Table summary
 FQID IF#      Freq(GHz)      BW(kHz)   Ch.Sep(kHz)  Sideband FQID IF#      Freq(GHz)      BW(kHz)   Ch.Sep(kHz)  Sideband
-   1   1       7.81490000   50000.0039  50000.0039      1+   1   1       8.51490000   50000.0039  50000.0039      1 
 +       2       8.88510000   50000.0039  50000.0039      1
 </file> </file>
  
Line 1290: Line 1360:
 bparm 0 1\\ bparm 0 1\\
 aparm 0\\ aparm 0\\
-refant 4\\+refant 3\\
 antennas 0\\ antennas 0\\
 do3col 1\\ do3col 1\\
Line 1355: Line 1425:
 ===== Flag some bad data.  Just run the commands and think about this later. ===== ===== Flag some bad data.  Just run the commands and think about this later. =====
  
 +
 +''
 +task 'quack'\\
 +reason 'quack'\\
 +sources ' '\\
 +calcode ' '\\
 +stokes ' '\\
 +timer 0\\
 +antenna 0\\
 +flagver 1\\
 +opcode 'beg'\\
 +basel 0\\
 +bif 0\\
 +eif 0\\
 +aparm 0\\
 +aparm(2) 22.5/60\\
 +go\\
 +''
  
 '' ''
Line 1360: Line 1448:
 default\\ default\\
 dohist 1\\ dohist 1\\
-getn 7\\ +getn 9\\ 
-antenna 4 0\\+antenna 0\\
 basel 0\\ basel 0\\
-timer 101 4 0 50 101 4 5 0\\+timer 101 4 3 50 101 4 4 20\\
 stokes ' '\\ stokes ' '\\
 +bif 0\\
 +eif 0\\
 opcode 'flag'\\ opcode 'flag'\\
-reason 'bad ant\\+reason 'dropout'\\
 outfgver 1\\ outfgver 1\\
 go\\ go\\
 +timer 0\\
 +bif 0\\
 +eif 0\\
 +antenna 0\\
 '' ''
- 
- 
  
  
Line 1383: Line 1475:
 task 'setjy'\\ task 'setjy'\\
 default\\ default\\
-getn 6\\+getn 9\\
 source '0518+165',' '\\ source '0518+165',' '\\
 bif 0\\ bif 0\\
Line 1401: Line 1493:
 vlb054> SETJY1: / Flux calculated using known spectrum vlb054> SETJY1: / Flux calculated using known spectrum
 vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs
-vlb054> SETJY1: '0518+165        ' IF =  1 FLUX = 3.7742 (Jy calcd) +vlb054> SETJY1: '0518+165        ' IF =  1 FLUX = 2.4983 (Jy calcd) 
-vlb054> SETJY1: '0518+165        ' IF =  2 FLUX = 3.9464 (Jy calcd)+vlb054> SETJY1: '0518+165        ' IF =  2 FLUX = 2.4166 (Jy calcd)
 vlb054> SETJY1: / Using (1995.2) 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
Line 1412: Line 1504:
 eif 1\\ eif 1\\
 optype ' '\\ optype ' '\\
-zerosp 3.7742, 0.380, -0.153, 0.0\\+zerosp 2.4983, 0.3092, -0.0628, 0.0\\
 go\\ go\\
 bif 2\\ bif 2\\
 eif 2\\ eif 2\\
-zerosp 3.9464, 0.397, -0.161, 0.0\\+zerosp 2.4166, 0.2990, -0.0607, 0.0\\
 go\\ go\\
 '' ''
Line 1432: Line 1524:
  
 <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: A source model for this calibrator may be available +
-vlb054> SETJY1: Use the verb CALDIR to see if there is one+
 vlb054> SETJY1: / Flux calculated using known spectrum vlb054> SETJY1: / Flux calculated using known spectrum
 vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs vlb054> SETJY1: BIF =  1 EIF =  2 /Range of IFs
-vlb054> SETJY1: '0134+329        ' IF =  1 FLUX = 5.4891 (Jy calcd) +vlb054> SETJY1: '0134+329        ' IF =  1 FLUX = 3.2045 (Jy calcd) 
-vlb054> SETJY1: '0134+329        ' IF =  2 FLUX = 5.8249 (Jy calcd)+vlb054> SETJY1: '0134+329        ' IF =  2 FLUX = 3.0714 (Jy calcd)
 vlb054> SETJY1: / Using (1995.2) 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>
  
 +
 +''
 +bif 1\\
 +eif 1\\
 +optype ' '\\
 +zerosp 3.2045, 0.0772, -0.1722, 0.0
 +go\\
 +bif 2\\
 +eif 2\\
 +zerosp 3.0714, 0.0740, -0.1650, 0.0
 +go\\
 +\\
 +bif 0\\
 +eif 0\\
 +''
 +
 +
 +
 +
 +==== Convert the (u,v) data to the epoch of the calibration images ====
 +
 +
 +Note: it is possible to change the epoch of the calibration image headers to B1950, but there would remain a small error in the rotation of the coordinate system.  The proper fix is to convert the (u,v) coordinates or to rotate the images.
 +
 +
 +''
 +task 'uvfix'\\
 +clron\\
 +outdisk 1\\
 +shift 0\\
 +uvfixprm 0\\
 +go\\
 +''
 +
 +Now fix the index table.
 +
 +''
 +pcat
 +''
 +
 +''
 +task 'indxr'\\
 +getn 16\\
 +cparm 0\\
 +cparm(3)= -1\\
 +bparm 0\\
 +go\\
 +cparm 0\\
 +''
  
 ==== Run the calibration for 3C138 ==== ==== Run the calibration for 3C138 ====
Line 1455: Line 1594:
 task 'calib'\\ task 'calib'\\
 default\\ default\\
-getn 13\\+getn 16\\
 calsour '0518+165',' '\\ calsour '0518+165',' '\\
 uvrange 0\\ uvrange 0\\
 antenna 0\\ antenna 0\\
-refant 4\\+refant 3\\
 minamper 10\\ minamper 10\\
 minphser 10\\ minphser 10\\
 weightit 1\\ weightit 1\\
-get2n 6\\+get2n 10\\
 ncomp 0\\ ncomp 0\\
 solmode 'A&P'\\ solmode 'A&P'\\
 aparm(6) 2\\ aparm(6) 2\\
-solint 2\\+solint 1\\
 solsub 2\\ solsub 2\\
 cparm(3) 10\\ cparm(3) 10\\
Line 1474: Line 1613:
 clrmsg\\ clrmsg\\
 '' ''
 +
 +==== Calibrate 3C48 ====
 +
 +''
 +get2n 11\\
 +calsour '0134+329',' '\\
 +solint 1\\
 +uvrange 0\\
 +go\\
 +clrmsg\\
 +''
 +
  
 ==== Calibrate the phase calibrator ==== ==== Calibrate the phase calibrator ====
 +
 +In this step we will try to get an amplitude calibration of the phase calibrator.  Only the first scan seems to be good in amplitude, so set the time range for just it.
 +
  
 '' ''
 clr2n\\ clr2n\\
 calsour '0420+417', ' '\\ calsour '0420+417', ' '\\
-solint 0\\+solint 40./60\\
 uvrange 0\\ uvrange 0\\
 +timer 101 0 0 0 101 8 0 0\\
 go\\ go\\
 clrmsg\\ clrmsg\\
 '' ''
 +
 ==== Check the results ==== ==== Check the results ====
  
Line 1507: Line 1663:
 go\\ go\\
 inver 2\\ inver 2\\
 +optype 'amp'\\
 +go\\
 +optype 'phas'\\
 +go\\
 +inver 3\\
 optype 'amp'\\ optype 'amp'\\
 go\\ go\\
Line 1520: Line 1681:
 task 'getjy'\\ task 'getjy'\\
 sources '0420+417', ' '\\ sources '0420+417', ' '\\
-calsour '0518+165',' '\\+calsour '0518+165','0134+329',' '\\
 calcode ' '\\ calcode ' '\\
 bif 0\\ bif 0\\
Line 1532: Line 1693:
  
 <file> <file>
-vlb054> GETJY1: Task GETJY  (release of 31DEC08) begins\\+vlb054> GETJY1: Task GETJY  (release of 31DEC10) begins 
 +vlb054> GETJY1: Calibrator robust averaging used    818 of    818 gain samples
 vlb054> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy) vlb054> GETJY1:    Source:Qual      CALCODE IF  Flux (Jy)
-vlb054> GETJY1: 0420+417        :  3   B     1    1.33893 +/-   0.00641 +vlb054> GETJY1: 0420+417        :  2   B     1    1.50919 +/-   0.00369  0.00369 
-vlb054> GETJY1:                              2    1.34714 +/-   0.00525+vlb054> GETJY1:                              2    1.50224 +/-   0.00401  0.00401 
 +vlb054> GETJY1:    Source:Qual    CALCODE    used   total     bad used  tot  bad 
 +vlb054> GETJY1: 0420+417        :  2 B        416     416       0  104  104    0
 vlb054> GETJY1: Appears to have ended successfully vlb054> GETJY1: Appears to have ended successfully
-vlb054> GETJY1: vlb054       31DEC08 TST: Cpu=       0.0  Real=       0+vlb054> GETJY1: vlb054       31DEC10 TST: Cpu=       0.0  Real=       0
 </file> </file>
  
Line 1543: Line 1707:
 ==== Redo the phase calibrator calibration ==== ==== Redo the phase calibrator calibration ====
  
-Now that we know the flux density of the phase calibrator, run through its calibration again, just to be sure AIPS doesn't have some information left around from calibrating before with an unknown flux density.+Now that we know the flux density of the phase calibrator, run through its calibration again, just to be sure AIPS doesn't have some information left around from calibrating before with an unknown flux density.  This time, run through the entire time range.
  
  
Line 1549: Line 1713:
 '' ''
 inext 'sn'\\ inext 'sn'\\
-inver 2\\+inver 3\\
 extdest\\ extdest\\
 inver 0\\ inver 0\\
 tget calib\\ tget calib\\
 +timer 0\\
 go\\ go\\
 clrmsg\\ clrmsg\\
Line 1562: Line 1727:
 '' ''
 task 'clcal'\\ task 'clcal'\\
-sources '0518+165','0420+417', '3C129', ' '\\ +sources '0518+165','0420+417', '0134+329','3C129', ' '\\ 
-calsour '0518+165','0420+417', ' '\\+calsour '0518+165','0420+417', '0134+329',' '\\
 opcode 'cali'\\ opcode 'cali'\\
 gainver 1\\ gainver 1\\
 gainuse 2\\ gainuse 2\\
-refant 4\\+refant 3\\
 interp '2pt'\\ interp '2pt'\\
 doblank 1\\ doblank 1\\
 go\\ go\\
 '' ''
 +
 ===== Check the results ===== ===== Check the results =====
  
Line 1616: Line 1782:
 task 'tasav'\\ task 'tasav'\\
 clro\\ clro\\
 +outdi 1\\
 go\\ go\\
 '' ''
 +
 ==== Run the polarization calibration on the phase calibrator ==== ==== Run the polarization calibration on the phase calibrator ====
  
Line 1627: Line 1795:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0\\+uvrange 20,0\\
 bif 0\\ bif 0\\
 eif 0\\ eif 0\\
Line 1637: Line 1805:
 soltype 'appr'\\ soltype 'appr'\\
 prtlev 1\\ prtlev 1\\
-refant 4\\+refant 3\\
 bparm 0\\ bparm 0\\
 cparm 0\\ cparm 0\\
Line 1646: Line 1814:
  
 <file> <file>
-vlb054> PCAL 1: Calibration source   1 +vlb054> PCAL 1: Calibration source   1   0420+417 
-vlb054> PCAL 1:   Q+iU=(   0.00444,  -0.04706) +/- (  0.000273,  0.000273) Jy +vlb054> PCAL 1:   Q+iU=(  -0.01935,   0.08692) +/- (  0.000536,  0.000536) Jy 
-vlb054> PCAL 1:   Pol. inten. =   0.04727 +/-  0.000386 Jy, angle = -42.30 +/- +vlb054> PCAL 1:   Pol. inten. =   0.08905 +/-  0.000758 Jy 
-vlb054>         0.165 deg +vlb054> PCAL 1:   Pol. angle  =  51.28 +/-  0.172 deg 
-vlb054> PCAL 1: 0420+417         I,Q,U,V =   1.3389  0.00444 -0.04706  0.00000 J +vlb054> PCAL 1: 0420+417         I,Q,U,V=   1.5092 -0.01935  0.08692  0.00000 Jy 
-vlb054>         y +vlb054> PCAL 1: Calibration source   1   0420+417 
-vlb054> PCAL 1: Calibration source   1 +vlb054> PCAL 1:   Q+iU=(   0.08929,  -0.01579) +/- (  0.000699,  0.000699) Jy 
-vlb054> PCAL 1:   Q+iU=(   0.02002,   0.03901) +/- (  0.000277,  0.000277) Jy +vlb054> PCAL 1:   Pol. inten. =   0.09068 +/-  0.000988 Jy 
-vlb054> PCAL 1:   Pol. inten. =   0.04385 +/-  0.000392 Jy, angle =  31.42 +/- +vlb054> PCAL 1:   Pol. angle  =  -5.02 +/-  0.221 deg 
-vlb054>         0.181 deg +vlb054> PCAL 1: 0420+417         I,Q,U,V=   1.5022  0.08929 -0.01579  0.00000 Jy
-vlb054> PCAL 1: 0420+417         I,Q,U,V =   1.3471  0.02002  0.03901  0.00000 J +
-vlb054>         y+
 </file> </file>
  
Line 1669: Line 1835:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0,45\\+uvrange 0\\
 docalib 1\\ docalib 1\\
 gainuse 2\\ gainuse 2\\
Line 1686: Line 1852:
 AIPS 1: Adverbs     Values                 Comments AIPS 1: Adverbs     Values                 Comments
 AIPS 1: ---------------------------------------------------------------- AIPS 1: ----------------------------------------------------------------
-AIPS 1: CLCORPRM    132.6554    -74.984    R-L phase difference returned+AIPS 1: CLCORPRM    -55.2632   -171.0893   R-L phase difference returned
 AIPS 1:            *rest 0 AIPS 1:            *rest 0
 </file> </file>
Line 1695: Line 1861:
 task 'tasav'\\ task 'tasav'\\
 clro\\ clro\\
 +outdi 1\\
 go\\ go\\
 '' ''
Line 1704: Line 1871:
 http://www.vla.nrao.edu/astro/calib/polar/1999/ http://www.vla.nrao.edu/astro/calib/polar/1999/
  
-Correct number is -21.71+Correct number is -22.98
  
  
Line 1710: Line 1877:
  
 '' ''
-for i=1 to 2;clcorp(i) = -21.71 - clcorp(i);end+for i=1 to 2;clcorp(i) = -22.98 - clcorp(i);end
 '' ''
 +
 ==== Apply the corrections ==== ==== Apply the corrections ====
  
Line 1737: Line 1905:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0,45\\+uvrange 0\\
 docalib 1\\ docalib 1\\
 gainuse 3\\ gainuse 3\\
Line 1766: Line 1934:
  
 ########################## ##########################
 +
 ===== Work on just the target data now ===== ===== Work on just the target data now =====
  
Line 1776: Line 1945:
 '' ''
 task 'split'\\ task 'split'\\
-getn 4\\ +getn 12\\ 
-source '3C129', ' '\\+source '0518+165','0134+329','0420+417','3C129', ' '\\
 timer 0\\ timer 0\\
 bif 0\\ bif 0\\
Line 1788: Line 1957:
 aparm 0\\ aparm 0\\
 go\\ go\\
-getn 13\\+getn 16\\
 go\\ go\\
 clrmsg\\ clrmsg\\
Line 1831: Line 2000:
 '' ''
 task 'dbcon'\\ task 'dbcon'\\
-getn 16\\ +getn 20\\ 
-get2n 17\\+get2n 26\\
 go\\ go\\
 pcat\\ pcat\\
Line 1839: Line 2008:
  
 <file> <file>
-AIPS 1:   18  101 3C129       .DBCON .    1 UV 13-NOV-2008 18:28:31+AIPS 1:   27  100 3C129       .DBCON .    1 UV 16-NOV-2010 22:13:56
 </file> </file>
 +
 ===== Imaging ===== ===== Imaging =====
  
Line 1851: Line 2021:
 '' ''
 task 'imagr'\\ task 'imagr'\\
-getn 18\\+default\\ 
 +getn 27\\
 stokes 'i'\\ stokes 'i'\\
 dopol 0\\ dopol 0\\
-cellsize 0.25,0.25\\+cellsize 0.15,0.15\\
 imsize 2048,2048\\ imsize 2048,2048\\
 clr2n\\ clr2n\\
Line 1862: Line 2033:
 niter 20000\\ niter 20000\\
 overlap 1\\ overlap 1\\
 +dotv 1\\
 go\\ go\\
 clrmsg\\ clrmsg\\
Line 1874: Line 2046:
 clrmsg\\ clrmsg\\
 '' ''
 +
 +
  
 ==== Check where the files went ==== ==== Check where the files went ====
  • Last modified: 2010-11-16 17:13
  • by James M Anderson