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 20:27] – 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 1448: Line 1448:
 default\\ default\\
 dohist 1\\ dohist 1\\
-getn 3\\ +getn 9\\ 
-antenna 14,0\\+antenna 0\\
 basel 0\\ basel 0\\
-timer 0 15 10 50 0 15 55 0\\+timer 101 4 3 50 101 4 4 20\\
 stokes ' '\\ stokes ' '\\
-bif 2\\ +bif 0\\ 
-eif 2\\+eif 0\\
 opcode 'flag'\\ opcode 'flag'\\
-reason 'bad ant'\\+reason 'dropout'\\
 outfgver 1\\ outfgver 1\\
-go\\ 
-timer 0 13 25 30 0 13 27 40\\ 
-go\\ 
-timer 0 13 59 0 0 14 2 20\\ 
 go\\ go\\
 timer 0\\ timer 0\\
Line 1479: Line 1475:
 task 'setjy'\\ task 'setjy'\\
 default\\ default\\
-getn 6\\+getn 9\\
 source '0518+165',' '\\ source '0518+165',' '\\
 bif 0\\ bif 0\\
Line 1497: 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 1508: 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 1528: 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 1551: 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 1570: 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 1603: Line 1663:
 go\\ go\\
 inver 2\\ inver 2\\
 +optype 'amp'\\
 +go\\
 +optype 'phas'\\
 +go\\
 +inver 3\\
 optype 'amp'\\ optype 'amp'\\
 go\\ go\\
Line 1616: 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 1628: 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 1639: 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 1645: 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 1658: 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 1712: 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 1723: Line 1795:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0\\+uvrange 20,0\\
 bif 0\\ bif 0\\
 eif 0\\ eif 0\\
Line 1733: Line 1805:
 soltype 'appr'\\ soltype 'appr'\\
 prtlev 1\\ prtlev 1\\
-refant 4\\+refant 3\\
 bparm 0\\ bparm 0\\
 cparm 0\\ cparm 0\\
Line 1742: 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 1765: Line 1835:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0,45\\+uvrange 0\\
 docalib 1\\ docalib 1\\
 gainuse 2\\ gainuse 2\\
Line 1782: 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 1791: Line 1861:
 task 'tasav'\\ task 'tasav'\\
 clro\\ clro\\
 +outdi 1\\
 go\\ go\\
 '' ''
Line 1800: 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 1806: 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 1833: Line 1905:
 timer 0\\ timer 0\\
 antenna 0\\ antenna 0\\
-uvrange 0,45\\+uvrange 0\\
 docalib 1\\ docalib 1\\
 gainuse 3\\ gainuse 3\\
Line 1862: Line 1934:
  
 ########################## ##########################
 +
 ===== Work on just the target data now ===== ===== Work on just the target data now =====
  
Line 1872: 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 1884: Line 1957:
 aparm 0\\ aparm 0\\
 go\\ go\\
-getn 13\\+getn 16\\
 go\\ go\\
 clrmsg\\ clrmsg\\
Line 1927: Line 2000:
 '' ''
 task 'dbcon'\\ task 'dbcon'\\
-getn 16\\ +getn 20\\ 
-get2n 17\\+get2n 26\\
 go\\ go\\
 pcat\\ pcat\\
Line 1935: 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 1947: 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 1958: Line 2033:
 niter 20000\\ niter 20000\\
 overlap 1\\ overlap 1\\
 +dotv 1\\
 go\\ go\\
 clrmsg\\ clrmsg\\
Line 1970: Line 2046:
 clrmsg\\ clrmsg\\
 '' ''
 +
 +
  
 ==== Check where the files went ==== ==== Check where the files went ====
  • Last modified: 2010-11-16 20:27
  • by James M Anderson