# CASA 4.3.1 stable prefix = '5B_Carray' msfileorig = '10C-119.sb8318174.eb9574290.56017.104540011576.ms' msfile= prefix+'_h.ms' dummy='' phasecal='' fluxcal='' polleakcal='' source='' myrefant='' mySPW='*' fluxcalname='3C286' #either 3C286 or 3C48 for our data # Run listobs # See output in oldmydatareduction.txt # Split galaxy, gain cal, flux+pol cal, leakage cal # Run listobs again since field names will have changed split(vis=msfileorig,outputvis='n4244_cl.ms',field='6,7,12,13',spw='', correlation='',datacolumn='all') # Can delete msfileorig now listobs(vis='n4244_cl.ms') 2015-05-19 17:01:22 INFO listobs ########################################## 2015-05-19 17:01:22 INFO listobs ##### Begin Task: listobs ##### 2015-05-19 17:01:22 INFO listobs listobs(vis="n4244_cl.ms",selectdata=True,spw="",field="",antenna="", 2015-05-19 17:01:22 INFO listobs uvrange="",timerange="",correlation="",scan="",intent="", 2015-05-19 17:01:22 INFO listobs feed="",array="",observation="",verbose=True,listfile="", 2015-05-19 17:01:22 INFO listobs listunfl=False,cachesize=50,overwrite=False) 2015-05-19 17:01:22 INFO listobs ================================================================================ 2015-05-19 17:01:22 INFO listobs MeasurementSet Name: /home/rjr/changes/n4244/carraylband/n4244_cl.ms MS Version 2 2015-05-19 17:01:22 INFO listobs ================================================================================ 2015-05-19 17:01:22 INFO listobs Observer: Dr. Judith Irwin Project: T.B.D. 2015-05-19 17:01:22 INFO listobs Observation: EVLA 2015-05-19 17:01:24 INFO listobs Data records: 6188832 Total elapsed time = 21390 seconds 2015-05-19 17:01:24 INFO listobs Observed from 31-Mar-2012/04:18:20.0 to 31-Mar-2012/10:14:50.0 (UTC) 2015-05-19 17:01:40 INFO listobs 2015-05-19 17:01:40 INFO listobs ObservationID = 0 ArrayID = 0 2015-05-19 17:01:40 INFO listobs Date Timerange (UTC) Scan FldId FieldName nRows SpwIds Average Interval(s) ScanIntent 2015-05-19 17:01:40 INFO listobs 31-Mar-2012/04:18:20.0 - 04:23:10.0 17 0 J1227+3635 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_AMPLI#UNSPECIFIED,CALIBRATE_PHASE#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 04:23:20.0 - 04:30:30.0 18 1 N4244 482976 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 04:30:40.0 - 04:40:00.0 19 1 N4244 628992 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 04:40:10.0 - 04:45:00.0 20 0 J1227+3635 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_AMPLI#UNSPECIFIED,CALIBRATE_PHASE#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 05:37:00.0 - 05:46:50.0 29 2 1331+305=3C286 662688 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_BANDPASS#UNSPECIFIED,CALIBRATE_FLUX#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED,UNSPECIFIED#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 05:47:00.0 - 05:56:50.0 30 3 J1407+2827 662688 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED,UNSPECIFIED#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 07:17:50.0 - 07:22:50.0 44 0 J1227+3635 336960 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_AMPLI#UNSPECIFIED,CALIBRATE_PHASE#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 07:23:00.0 - 07:27:50.0 45 1 N4244 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 07:28:00.0 - 07:36:50.0 46 1 N4244 595296 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 07:37:00.0 - 07:41:20.0 47 0 J1227+3635 292032 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_AMPLI#UNSPECIFIED,CALIBRATE_PHASE#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 09:51:20.0 - 09:56:10.0 69 0 J1227+3635 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_AMPLI#UNSPECIFIED,CALIBRATE_PHASE#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 09:56:20.0 - 10:01:10.0 70 1 N4244 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 10:01:20.0 - 10:10:10.0 71 1 N4244 595296 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs 10:10:20.0 - 10:14:50.0 72 0 J1227+3635 303264 [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [CALIBRATE_AMPLI#UNSPECIFIED,CALIBRATE_PHASE#UNSPECIFIED,OBSERVE_TARGET#UNSPECIFIED] 2015-05-19 17:01:40 INFO listobs (nRows = Total number of rows per scan) 2015-05-19 17:01:41 INFO listobs Fields: 4 2015-05-19 17:01:41 INFO listobs ID Code Name RA Decl Epoch SrcId nRows 2015-05-19 17:01:41 INFO listobs 0 D J1227+3635 12:27:58.725445 +36.35.11.82734 J2000 0 1909440 2015-05-19 17:01:41 INFO listobs 1 NONE N4244 12:17:29.660000 +37.48.25.60001 J2000 1 2954016 2015-05-19 17:01:41 INFO listobs 2 N 1331+305=3C286 13:31:08.287984 +30.30.32.95886 J2000 2 662688 2015-05-19 17:01:41 INFO listobs 3 H J1407+2827 14:07:00.394410 +28.27.14.68997 J2000 3 662688 2015-05-19 17:01:42 INFO listobs Spectral Windows: (32 unique spectral windows and 1 unique polarization setups) 2015-05-19 17:01:42 INFO listobs SpwID Name #Chans Frame Ch0(MHz) ChanWid(kHz) TotBW(kHz) CtrFreq(MHz) BBC Num Corrs 2015-05-19 17:01:42 INFO listobs 0 Subband:0 64 TOPO 1247.000 250.000 16000.0 1254.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 1 Subband:1 64 TOPO 1263.000 250.000 16000.0 1270.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 2 Subband:2 64 TOPO 1279.000 250.000 16000.0 1286.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 3 Subband:3 64 TOPO 1295.000 250.000 16000.0 1302.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 4 Subband:4 64 TOPO 1311.000 250.000 16000.0 1318.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 5 Subband:5 64 TOPO 1327.000 250.000 16000.0 1334.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 6 Subband:6 64 TOPO 1343.000 250.000 16000.0 1350.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 7 Subband:7 64 TOPO 1359.000 250.000 16000.0 1366.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 8 Subband:8 64 TOPO 1375.000 250.000 16000.0 1382.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 9 Subband:9 64 TOPO 1391.000 250.000 16000.0 1398.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 10 Subband:10 64 TOPO 1407.000 250.000 16000.0 1414.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 11 Subband:11 64 TOPO 1423.000 250.000 16000.0 1430.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 12 Subband:12 64 TOPO 1439.000 250.000 16000.0 1446.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 13 Subband:13 64 TOPO 1455.000 250.000 16000.0 1462.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 14 Subband:14 64 TOPO 1471.000 250.000 16000.0 1478.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 15 Subband:15 64 TOPO 1487.000 250.000 16000.0 1494.8750 12 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 16 Subband:0 64 TOPO 1647.000 250.000 16000.0 1654.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 17 Subband:1 64 TOPO 1663.000 250.000 16000.0 1670.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 18 Subband:2 64 TOPO 1679.000 250.000 16000.0 1686.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 19 Subband:3 64 TOPO 1695.000 250.000 16000.0 1702.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 20 Subband:4 64 TOPO 1711.000 250.000 16000.0 1718.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 21 Subband:5 64 TOPO 1727.000 250.000 16000.0 1734.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 22 Subband:6 64 TOPO 1743.000 250.000 16000.0 1750.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 23 Subband:7 64 TOPO 1759.000 250.000 16000.0 1766.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 24 Subband:8 64 TOPO 1775.000 250.000 16000.0 1782.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 25 Subband:9 64 TOPO 1791.000 250.000 16000.0 1798.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 26 Subband:10 64 TOPO 1807.000 250.000 16000.0 1814.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 27 Subband:11 64 TOPO 1823.000 250.000 16000.0 1830.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 28 Subband:12 64 TOPO 1839.000 250.000 16000.0 1846.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 29 Subband:13 64 TOPO 1855.000 250.000 16000.0 1862.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 30 Subband:14 64 TOPO 1871.000 250.000 16000.0 1878.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs 31 Subband:15 64 TOPO 1887.000 250.000 16000.0 1894.8750 15 RR RL LR LL 2015-05-19 17:01:42 INFO listobs Sources: 128 2015-05-19 17:01:42 INFO listobs ID Name SpwId RestFreq(MHz) SysVel(km/s) 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 0 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 1 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 2 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 3 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 4 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 5 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 6 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 7 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 8 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 9 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 10 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 11 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 12 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 13 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 14 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 15 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 16 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 17 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 18 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 19 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 20 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 21 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 22 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 23 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 24 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 25 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 26 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 27 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 28 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 29 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 30 - - 2015-05-19 17:01:42 INFO listobs 0 J1227+3635 31 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 0 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 1 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 2 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 3 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 4 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 5 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 6 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 7 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 8 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 9 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 10 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 11 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 12 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 13 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 14 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 15 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 16 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 17 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 18 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 19 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 20 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 21 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 22 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 23 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 24 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 25 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 26 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 27 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 28 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 29 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 30 - - 2015-05-19 17:01:42 INFO listobs 1 N4244 31 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 0 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 1 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 2 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 3 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 4 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 5 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 6 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 7 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 8 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 9 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 10 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 11 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 12 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 13 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 14 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 15 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 16 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 17 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 18 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 19 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 20 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 21 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 22 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 23 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 24 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 25 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 26 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 27 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 28 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 29 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 30 - - 2015-05-19 17:01:42 INFO listobs 2 1331+305=3C286 31 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 0 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 1 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 2 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 3 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 4 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 5 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 6 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 7 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 8 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 9 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 10 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 11 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 12 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 13 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 14 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 15 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 16 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 17 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 18 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 19 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 20 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 21 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 22 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 23 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 24 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 25 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 26 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 27 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 28 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 29 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 30 - - 2015-05-19 17:01:42 INFO listobs 3 J1407+2827 31 - - 2015-05-19 17:01:45 INFO listobs Antennas: 27: 2015-05-19 17:01:45 INFO listobs ID Name Station Diam. Long. Lat. Offset from array center (m) ITRF Geocentric coordinates (m) 2015-05-19 17:01:45 INFO listobs East North Elevation x y z 2015-05-19 17:01:45 INFO listobs 0 ea01 N08 25.0 m -107.37.07.5 +33.54.15.8 -68.9092 433.1903 -5.0594 -1601147.943700 -5041733.835900 3555235.957600 2015-05-19 17:01:45 INFO listobs 1 ea02 E06 25.0 m -107.36.55.6 +33.53.57.7 236.9148 -126.3352 -2.4529 -1600951.577000 -5042125.906000 3554773.008900 2015-05-19 17:01:45 INFO listobs 2 ea03 E18 25.0 m -107.35.57.2 +33.53.35.1 1732.7282 -823.9896 3.1267 -1599644.861100 -5042953.662300 3554197.033200 2015-05-19 17:01:45 INFO listobs 3 ea04 W14 25.0 m -107.37.46.9 +33.53.38.9 -1078.9030 -708.5656 1.1470 -1602304.909100 -5042039.849500 3554291.739800 2015-05-19 17:01:45 INFO listobs 4 ea05 W08 25.0 m -107.37.21.6 +33.53.53.0 -432.1178 -272.1466 -1.5043 -1601614.092200 -5042001.652900 3554652.510900 2015-05-19 17:01:45 INFO listobs 5 ea06 N10 25.0 m -107.37.08.2 +33.54.22.4 -86.6491 636.0548 -6.1307 -1601130.329000 -5041619.765300 3555403.729400 2015-05-19 17:01:45 INFO listobs 6 ea07 E04 25.0 m -107.37.00.8 +33.53.59.7 102.8106 -63.7722 -2.6333 -1601068.788100 -5042051.920300 3554824.836500 2015-05-19 17:01:45 INFO listobs 7 ea08 N12 25.0 m -107.37.09.0 +33.54.30.0 -107.1583 870.2666 -7.3175 -1601110.028900 -5041488.093200 3555597.449200 2015-05-19 17:01:45 INFO listobs 8 ea09 E02 25.0 m -107.37.04.4 +33.54.01.1 9.8170 -20.4321 -2.7728 -1601150.069400 -5042000.625300 3554860.731500 2015-05-19 17:01:45 INFO listobs 9 ea10 N14 25.0 m -107.37.09.9 +33.54.38.5 -130.2432 1134.2226 -8.4603 -1601087.168300 -5041339.854000 3555815.873700 2015-05-19 17:01:45 INFO listobs 10 ea11 W16 25.0 m -107.37.57.4 +33.53.33.0 -1348.7206 -890.6160 1.2834 -1602592.857600 -5042054.977300 3554140.699900 2015-05-19 17:01:45 INFO listobs 11 ea12 E08 25.0 m -107.36.48.9 +33.53.55.1 407.8285 -206.0065 -3.2272 -1600801.926000 -5042219.366500 3554706.448200 2015-05-19 17:01:45 INFO listobs 12 ea13 W09 25.0 m -107.37.25.2 +33.53.51.0 -521.9458 -332.7614 -1.1878 -1601710.021500 -5042006.925500 3554602.375100 2015-05-19 17:01:45 INFO listobs 13 ea14 W12 25.0 m -107.37.37.4 +33.53.44.2 -835.3706 -544.2307 0.5891 -1602044.903400 -5042025.823600 3554427.836900 2015-05-19 17:01:45 INFO listobs 14 ea15 W18 25.0 m -107.38.08.9 +33.53.26.5 -1643.9568 -1089.8836 2.3855 -1602908.188700 -5042072.296200 3553975.903500 2015-05-19 17:01:45 INFO listobs 15 ea16 N04 25.0 m -107.37.06.5 +33.54.06.1 -42.6378 132.8529 -3.5476 -1601173.991500 -5041902.650000 3554987.526200 2015-05-19 17:01:45 INFO listobs 16 ea17 E12 25.0 m -107.36.31.7 +33.53.48.5 848.6927 -411.6117 -2.7480 -1600416.525000 -5042462.430200 3554536.056100 2015-05-19 17:01:45 INFO listobs 17 ea18 E14 25.0 m -107.36.21.3 +33.53.44.5 1114.9846 -535.8413 -1.5280 -1600183.959300 -5042609.986600 3554433.619800 2015-05-19 17:01:45 INFO listobs 18 ea19 W02 25.0 m -107.37.07.5 +33.54.00.9 -67.9749 -26.5393 -2.7144 -1601225.258200 -5041980.372390 3554855.695000 2015-05-19 17:01:45 INFO listobs 19 ea20 N18 25.0 m -107.37.12.0 +33.54.58.3 -183.8798 1746.4967 -10.3485 -1601034.401400 -5040996.541900 3556322.934300 2015-05-19 17:01:45 INFO listobs 20 ea21 E10 25.0 m -107.36.40.9 +33.53.52.0 612.4523 -301.4708 -3.5406 -1600622.920100 -5042331.781200 3554627.035200 2015-05-19 17:01:45 INFO listobs 21 ea22 N02 25.0 m -107.37.06.2 +33.54.03.5 -35.6384 53.1810 -3.1429 -1601180.872480 -5041947.442400 3554921.624300 2015-05-19 17:01:45 INFO listobs 22 ea23 N06 25.0 m -107.37.06.9 +33.54.10.3 -54.0788 263.8798 -4.2322 -1601162.602700 -5041828.989900 3555095.895300 2015-05-19 17:01:45 INFO listobs 23 ea25 W04 25.0 m -107.37.10.8 +33.53.59.1 -152.8575 -83.8015 -2.4633 -1601315.889600 -5041985.317370 3554808.306800 2015-05-19 17:01:45 INFO listobs 24 ea26 W06 25.0 m -107.37.15.6 +33.53.56.4 -275.8441 -166.7260 -2.0605 -1601447.209000 -5041992.486500 3554739.702600 2015-05-19 17:01:45 INFO listobs 25 ea27 E16 25.0 m -107.36.09.8 +33.53.40.0 1410.0320 -673.4683 -0.7927 -1599926.109000 -5042772.965300 3554319.791200 2015-05-19 17:01:45 INFO listobs 26 ea28 N16 25.0 m -107.37.10.9 +33.54.48.0 -155.8610 1426.6357 -9.3817 -1601061.967000 -5041175.882700 3556058.031600 2015-05-19 17:01:45 INFO listobs ##### End Task: listobs ##### #Automatically find antenna position updates: gencal(vis='n4244_cl.ms',caltable ...: ='n4244_cl.antpos',caltype='antpos') #Result is small offsets for ea01, as log.ps said. #Create cal table that calibrates for antenna shape deformity vs. elevation. Apply in all later cals: gencal(vis='n4244_cl.ms',caltable='n4244_cl.gaincurve',caltype='gc') #Generate ionospheric correction table to be used in calibration: from recipes import tec_maps #Generate "TEC surface image". Also plot of TEC vs time at zenith. Do screenshot to save it. This is tec.png #This died the first time because it was looking for a file 'igsg0910.12i' which wasn't there. However, in my #directory there was a .Z version of this file. Upon uncompressing it and rerunning the program, it worked. tec_maps.create(vis='n4244_cl.ms',doplot=T,imname='tec.im') #Generate the ionospheric correction caltable, to be used as a 'gaintable' parameter during next calibration. gencal(vis='n4244_cl.ms',caltable='tec.cal',caltype='tecim',infile='tec.im.IGS_TEC.im') # 'tec.cal' should now be applied as a prior calibration in subsequent calibration solves using parameter gaintable # and blank for gainfield, and applied in the final applycal #Delay calibration. Use 3c286, but there is heavy interference so it needs to be flagged first, so follow Sec 4.10 #Split 3c286 from rest: split(vis='n4244_cl.ms',outputvis='3c286_h.ms',field=2,datacolumn='data') #Hanning smooth the result: hanningsmooth(vis='3c286_h.ms') #Flag bad data from 3c286 hanning-smoothed file. Use Arpad's email, but skip some, e.g. ea10 looks fine, plus more, e.g. ea01,ea05,ea23 should be completely flagged. #Flag 5 end channels as he did, some others slightly modified. flagdata commands not included in this copy: #Use ea22 as refant right from the start, not ea19 as in first calibration last year. #Create delay table on hanning smoothed 3c286 ms according to Sec 4.8. Now including tec.cal in gaintable, but not sure it matters for delays. gaincal(vis='3c286_h.ms',caltable='3c286_h.K0',gaintable=['n4244_cl.antpos','n4244_cl.gaincurve','tec.cal'],field='',spw='',gaintype='K',refant='ea22',combine='scan',solint='60s',minsnr=5,selectdata=T, timerange='05:40:00~05:42:00') #Plot delays vs. antenna for poln LL and RR in turn. Delays should be a few ms or less. They all are in this case. plotcal(caltable='3c286_h.K0',xaxis='time',yaxis='delay',poln='L',spw='',subplot=221,iteration='antenna') #Apply delay cal to original ms file applycal(vis='n4244_cl.ms',field='',gaintable=['3c286_h.K0','n4244_cl.antpos'],interp='nearest',calwt=False) #Split delay corrected data into new ms split(vis='n4244_cl.ms',outputvis='n4244_cl_delaycaled.ms',datacolumn='corrected') #Hanning smooth this new ms, into another new ms. hanningsmooth(vis='n4244_cl_delaycaled.ms',outputvis='n4244_cl_delay_h.ms',datacolumn='data') #Delete n4244_cl_delaycaled.ms # # #DATA EDITING - commands not included here # # # # # STOKES I CALIBRATION # # # Set flux density for 3C286 setjy(vis='n4244_cl_delay_h.ms',field='2',standard='Perley-Butler 2010',modimage='/home/rjr/changes/n4244/carraylband/3C286_L.im',scalebychan=True) 2015-06-05 20:48:33 INFO setjy ########################################## 2015-06-05 20:48:33 INFO setjy ##### Begin Task: setjy ##### 2015-06-05 20:48:33 INFO setjy setjy(vis="n4244_cl_delay_h.ms",field="2",spw="",selectdata=False,timerange="", 2015-06-05 20:48:33 INFO setjy scan="",intent="",observation="",scalebychan=True,standard="Perley-Butler 2010", 2015-06-05 20:48:33 INFO setjy model="",modimage="/home/rjr/changes/n4244/carraylband/3C286_L.im",listmodels=False,fluxdensity=-1,spix=0.0, 2015-06-05 20:48:33 INFO setjy reffreq="1GHz",polindex=[],polangle=[],rotmeas=0.0,fluxdict={}, 2015-06-05 20:48:33 INFO setjy useephemdir=False,interpolation="nearest",usescratch=False) 2015-06-05 20:48:35 INFO Imager Opening MeasurementSet /home/rjr/changes/n4244/carraylband/n4244_cl_delay_h.ms 2015-06-05 20:48:43 INFO setjy The modimage parameter is deprecated please use model instead 2015-06-05 20:48:47 INFO imager Using channel dependent flux densities 2015-06-05 20:48:50 INFO imager Selected 662688 out of 6188832 rows. 2015-06-05 20:48:50 INFO imager 1331+305=3C286 (fld ind 2) spw 0 [I=15.662, Q=0, U=0, V=0] Jy @ 1.247e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:50 INFO imager 1331+305=3C286 (fld ind 2) spw 1 [I=15.573, Q=0, U=0, V=0] Jy @ 1.263e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:50 INFO imager 1331+305=3C286 (fld ind 2) spw 2 [I=15.485, Q=0, U=0, V=0] Jy @ 1.279e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:50 INFO imager 1331+305=3C286 (fld ind 2) spw 3 [I=15.398, Q=0, U=0, V=0] Jy @ 1.295e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:50 INFO imager 1331+305=3C286 (fld ind 2) spw 4 [I=15.312, Q=0, U=0, V=0] Jy @ 1.311e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:51 INFO imager 1331+305=3C286 (fld ind 2) spw 5 [I=15.228, Q=0, U=0, V=0] Jy @ 1.327e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:51 INFO imager 1331+305=3C286 (fld ind 2) spw 6 [I=15.144, Q=0, U=0, V=0] Jy @ 1.343e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:51 INFO imager 1331+305=3C286 (fld ind 2) spw 7 [I=15.062, Q=0, U=0, V=0] Jy @ 1.359e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:51 INFO imager 1331+305=3C286 (fld ind 2) spw 8 [I=14.981, Q=0, U=0, V=0] Jy @ 1.375e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:51 INFO imager 1331+305=3C286 (fld ind 2) spw 9 [I=14.901, Q=0, U=0, V=0] Jy @ 1.391e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:51 INFO imager 1331+305=3C286 (fld ind 2) spw 10 [I=14.822, Q=0, U=0, V=0] Jy @ 1.407e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:52 INFO imager 1331+305=3C286 (fld ind 2) spw 11 [I=14.744, Q=0, U=0, V=0] Jy @ 1.423e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:52 INFO imager 1331+305=3C286 (fld ind 2) spw 12 [I=14.667, Q=0, U=0, V=0] Jy @ 1.439e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:52 INFO imager 1331+305=3C286 (fld ind 2) spw 13 [I=14.59, Q=0, U=0, V=0] Jy @ 1.455e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:52 INFO imager 1331+305=3C286 (fld ind 2) spw 14 [I=14.515, Q=0, U=0, V=0] Jy @ 1.471e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:52 INFO imager 1331+305=3C286 (fld ind 2) spw 15 [I=14.441, Q=0, U=0, V=0] Jy @ 1.487e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:52 INFO imager 1331+305=3C286 (fld ind 2) spw 16 [I=13.747, Q=0, U=0, V=0] Jy @ 1.647e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:53 INFO imager 1331+305=3C286 (fld ind 2) spw 17 [I=13.682, Q=0, U=0, V=0] Jy @ 1.663e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:53 INFO imager 1331+305=3C286 (fld ind 2) spw 18 [I=13.618, Q=0, U=0, V=0] Jy @ 1.679e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:53 INFO imager 1331+305=3C286 (fld ind 2) spw 19 [I=13.555, Q=0, U=0, V=0] Jy @ 1.695e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:53 INFO imager 1331+305=3C286 (fld ind 2) spw 20 [I=13.492, Q=0, U=0, V=0] Jy @ 1.711e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:53 INFO imager 1331+305=3C286 (fld ind 2) spw 21 [I=13.43, Q=0, U=0, V=0] Jy @ 1.727e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:53 INFO imager 1331+305=3C286 (fld ind 2) spw 22 [I=13.368, Q=0, U=0, V=0] Jy @ 1.743e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:54 INFO imager 1331+305=3C286 (fld ind 2) spw 23 [I=13.308, Q=0, U=0, V=0] Jy @ 1.759e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:54 INFO imager 1331+305=3C286 (fld ind 2) spw 24 [I=13.248, Q=0, U=0, V=0] Jy @ 1.775e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:54 INFO imager 1331+305=3C286 (fld ind 2) spw 25 [I=13.188, Q=0, U=0, V=0] Jy @ 1.791e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:54 INFO imager 1331+305=3C286 (fld ind 2) spw 26 [I=13.129, Q=0, U=0, V=0] Jy @ 1.807e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:54 INFO imager 1331+305=3C286 (fld ind 2) spw 27 [I=13.071, Q=0, U=0, V=0] Jy @ 1.823e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:54 INFO imager 1331+305=3C286 (fld ind 2) spw 28 [I=13.014, Q=0, U=0, V=0] Jy @ 1.839e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:55 INFO imager 1331+305=3C286 (fld ind 2) spw 29 [I=12.957, Q=0, U=0, V=0] Jy @ 1.855e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:55 INFO imager 1331+305=3C286 (fld ind 2) spw 30 [I=12.901, Q=0, U=0, V=0] Jy @ 1.871e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:55 INFO imager 1331+305=3C286 (fld ind 2) spw 31 [I=12.845, Q=0, U=0, V=0] Jy @ 1.887e+09Hz, (Perley-Butler 2010) 2015-06-05 20:48:59 INFO imager Using model image /home/rjr/changes/n4244/carraylband/3C286_L.im 2015-06-05 20:48:59 INFO imager Scaling spw(s) [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31]'s model image by channel to I = 15.663, 14.05, 12.7912 Jy @(1.24683e+09, 1.57483e+09, 1.90258e+09)Hz for visibility prediction (a few representative values are shown). 2015-06-05 20:48:59 INFO imager The model image's reference pixel is 0.000254726 arcsec from 1331+305=3C286's phase center. 2015-06-05 20:48:59 INFO imager Will clear any existing model with matching field=1331+305=3C286 and spw=* 2015-06-05 20:48:59 INFO Clearing model records in MS header for selected fields. 2015-06-05 20:48:59 INFO 1331+305=3C286 (id = 2) not found. 2015-06-05 20:49:00 INFO imager Selected 20709 out of 6188832 rows. 2015-06-05 20:49:04 INFO imager Selected 20709 out of 6188832 rows. etc. 2015-06-05 20:50:18 INFO setjy ##### End Task: setjy ##### 2015-06-05 20:50:18 INFO setjy ########################################## #Initial gain calibration for passband calibration, but now with refant 22 and applying tec.cal on the fly gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G0',field='2',refant='ea22',spw='0~31:15~50',calmode='p',solint='int',minsnr=5,solnorm=T,gaintable=['tec.cal']) #Rerun gaincal. Stick with G0. By default, append=F and it will overwrite the existing caltable. gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G0',field='2',refant='ea22',spw='0~31:15~50',calmode='p',solint='int',minsnr=5,solnorm=T,gaintable=['tec.cal']) #Reports all attempted calibrations succeeded #Plot gain phase vs time for all SPWs for each antenna in turn. plotcal(caltable='n4244_cl_delay_h.G0',xaxis='time',yaxis='phase',poln='L',field='2',spw='0~31',iteration='antenna') #All look smooth with time. Only one scan. plotcal(caltable='n4244_cl_delay_h.G0',xaxis='time',yaxis='phase',poln='R',field='2',spw='0~31',iteration='antenna') #Okay ######################## #Bandpass cal. refant 22, tec.cal, G0 gain cal applied OTF. gaincurve parameter no longer exists. This was the elevation gain correction curve. Presumably #applied automatically now. The casa webpages sort of imply that we should have been doing this for VLA data anyway when the parameter existed. bandpass(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.B0',field='2',spw='',refant='ea22',solnorm=False,combine='scan',solint='inf',bandtype='B',gaintable=['n4244_cl_delay_h.G0','tec.cal']) #Check solns plotcal(caltable='n4244_cl_delay_h.B0',poln='L',xaxis='chan',yaxis='amp',field='2',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.B0',poln='R',xaxis='chan',yaxis='amp',field='2',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.B0',poln='L',xaxis='chan',yaxis='phase',field='2',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.B0',poln='R',xaxis='chan',yaxis='phase',field='2',subplot=221,spw='',iteration='antenna') #Seem ok. I don't think I see the problem from last time: "For about half the antennas, spw13, L and R, chs 30-45 solns show a phase variation of up to 10 deg. #Can see in plotms. Exists for all times during scan. Unclear if it exists for all tracks." ################################ #Gain calibration. Use refant 22, tec.cal OTF, create version G1, not G2 as last time. gaincurve doesn't exist anymore #Phases first. First for fluxcal. Do solns for 10 sec integrations for some reason, apparently because these phase solns will only get used for calibrating the amp # of the flux cal later. gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G1',solint='int',field='2',spw='0~31:5~58',refant='ea22',gaintype='G',calmode='p',solnorm=F,gaintable=['n4244_cl_delay_h.B0','tec.cal']) #For other cals, solint=inf and combine='' apparently means do a solution for each scan. Append=T means appending onto previous table gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G1',solint='inf',field='0,3',spw='0~31:5~58',refant='ea22',gaintype='G',calmode='p',solnorm=F,append=T,gaintable=['n4244_cl_delay_h.B0','tec.cal']) #Plot phase solns for each cal and poln. Does something funny with polcal (field 2) when xaxis=time. With solint=inf, there is only one sol per scan, so there should be #only one sol time plotted, but it splits the antennae up into three times, even though it doesn't plot a time axis. With xaxis=scan, it looks the same as the flux cal plot #plotted the same way. May be a bug in how it handles time when scan averaging sols. Data look fine in ms, if much noisier than flux cal. plotcal(caltable='n4244_cl_delay_h.G1',poln='R',xaxis='time',yaxis='phase',field='2',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.G1',poln='L',xaxis='time',yaxis='phase',field='2',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.G1',poln='R',xaxis='time',yaxis='phase',field='3',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.G1',poln='L',xaxis='time',yaxis='phase',field='3',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.G1',poln='R',xaxis='time',yaxis='phase',field='0',subplot=221,spw='',iteration='antenna') plotcal(caltable='n4244_cl_delay_h.G1',poln='L',xaxis='time',yaxis='phase',field='0',subplot=221,spw='',iteration='antenna') #Now amp cal on a per scan basis for each cal in turn. Gainfield param means use 3c286 from passband gaintable and the phases for the field for the #source being amp calibrated from the gain table G1. Applying B0, G1, tec.cal OTF, refant 22, creating G2, not G3 as before. gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G2',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','tec.cal'],gainfield=['2','2',''],field='2',refant='ea22',solnorm=F,spw='0~31:5~58',solint='inf',gaintype='G',calmode='a') gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G2',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','tec.cal'],gainfield=['2','0',''],field='0',refant='ea22',solnorm=F,spw='0~31:5~58',solint='inf',gaintype='G',calmode='a',append=T) gaincal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G2',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','tec.cal'],gainfield=['2','3',''],field='3',refant='ea22',solnorm=F,spw='0~31:5~58',solint='inf',gaintype='G',calmode='a',append=T) # Find other cal fluxes from 3C286. For 1227, we get 1.8-2.1 Jy. It is 2.1 at 20cm according to cal manual. For 1407, we get 0.8-1.3 Jy. Cal manual has 0.8 Jy. fluxscale(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.G2',fluxtable='n4244_cl_delay_h.F1',reference=['2'],transfer=['0,3']) #Resulting fluxes output 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=0 (freq=1.247e+09 Hz) is: 2.11091 +/- 0.0232638 (SNR = 90.738, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=1 is: INSUFFICIENT DATA 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=2 is: INSUFFICIENT DATA 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=3 (freq=1.295e+09 Hz) is: 2.08742 +/- 0.0172723 (SNR = 120.854, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=4 (freq=1.311e+09 Hz) is: 2.0742 +/- 0.0104032 (SNR = 199.381, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=5 (freq=1.327e+09 Hz) is: 2.0705 +/- 0.0109704 (SNR = 188.735, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=6 (freq=1.343e+09 Hz) is: 2.06161 +/- 0.0111965 (SNR = 184.13, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=7 (freq=1.359e+09 Hz) is: 2.05142 +/- 0.00331967 (SNR = 617.96, N = 42) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=8 (freq=1.375e+09 Hz) is: 2.0469 +/- 0.00323993 (SNR = 631.773, N = 42) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=9 (freq=1.391e+09 Hz) is: 2.05101 +/- 0.0108496 (SNR = 189.039, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=10 (freq=1.407e+09 Hz) is: 2.04319 +/- 0.011283 (SNR = 181.086, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=11 (freq=1.423e+09 Hz) is: 2.03323 +/- 0.0108927 (SNR = 186.66, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=12 (freq=1.439e+09 Hz) is: 2.02506 +/- 0.0118269 (SNR = 171.225, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=13 (freq=1.455e+09 Hz) is: 2.02364 +/- 0.0122409 (SNR = 165.317, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=14 (freq=1.471e+09 Hz) is: 2.01679 +/- 0.0126728 (SNR = 159.143, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=15 (freq=1.487e+09 Hz) is: 2.01186 +/- 0.0127855 (SNR = 157.355, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=16 (freq=1.647e+09 Hz) is: 1.95225 +/- 0.00891864 (SNR = 218.896, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=17 (freq=1.663e+09 Hz) is: 1.95592 +/- 0.00890262 (SNR = 219.702, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=18 is: INSUFFICIENT DATA 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=19 (freq=1.695e+09 Hz) is: 1.93297 +/- 0.00964557 (SNR = 200.399, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=20 (freq=1.711e+09 Hz) is: 1.92113 +/- 0.0107628 (SNR = 178.498, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=21 (freq=1.727e+09 Hz) is: 1.91503 +/- 0.011588 (SNR = 165.259, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=22 (freq=1.743e+09 Hz) is: 1.90746 +/- 0.0114643 (SNR = 166.382, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=23 (freq=1.759e+09 Hz) is: 1.89744 +/- 0.00995411 (SNR = 190.619, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=24 (freq=1.775e+09 Hz) is: 1.89032 +/- 0.00875055 (SNR = 216.023, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=25 (freq=1.791e+09 Hz) is: 1.87989 +/- 0.00879656 (SNR = 213.707, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=26 (freq=1.807e+09 Hz) is: 1.87486 +/- 0.00903426 (SNR = 207.528, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=27 (freq=1.823e+09 Hz) is: 1.8628 +/- 0.00936346 (SNR = 198.944, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=28 (freq=1.839e+09 Hz) is: 1.85756 +/- 0.00916382 (SNR = 202.706, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=29 (freq=1.855e+09 Hz) is: 1.85269 +/- 0.00905527 (SNR = 204.598, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=30 (freq=1.871e+09 Hz) is: 1.84522 +/- 0.00718614 (SNR = 256.774, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1227+3635 in SpW=31 (freq=1.887e+09 Hz) is: 1.8339 +/- 0.00928353 (SNR = 197.544, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=0 (freq=1.247e+09 Hz) is: 0.745819 +/- 0.00331061 (SNR = 225.281, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=1 is: INSUFFICIENT DATA 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=2 is: INSUFFICIENT DATA 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=3 (freq=1.295e+09 Hz) is: 0.78149 +/- 0.00344599 (SNR = 226.783, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=4 (freq=1.311e+09 Hz) is: 0.794661 +/- 0.00334695 (SNR = 237.428, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=5 (freq=1.327e+09 Hz) is: 0.812757 +/- 0.0035138 (SNR = 231.305, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=6 (freq=1.343e+09 Hz) is: 0.827035 +/- 0.00350097 (SNR = 236.23, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=7 (freq=1.359e+09 Hz) is: 0.838588 +/- 0.00240574 (SNR = 348.577, N = 42) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=8 (freq=1.375e+09 Hz) is: 0.858833 +/- 0.00288527 (SNR = 297.661, N = 42) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=9 (freq=1.391e+09 Hz) is: 0.872482 +/- 0.00346463 (SNR = 251.825, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=10 (freq=1.407e+09 Hz) is: 0.890895 +/- 0.00363012 (SNR = 245.418, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=11 (freq=1.423e+09 Hz) is: 0.905105 +/- 0.00366366 (SNR = 247.05, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=12 (freq=1.439e+09 Hz) is: 0.916393 +/- 0.00338063 (SNR = 271.072, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=13 (freq=1.455e+09 Hz) is: 0.931102 +/- 0.00309768 (SNR = 300.58, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=14 (freq=1.471e+09 Hz) is: 0.942962 +/- 0.00308483 (SNR = 305.677, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=15 (freq=1.487e+09 Hz) is: 0.95551 +/- 0.00362046 (SNR = 263.919, N = 44) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=16 (freq=1.647e+09 Hz) is: 1.10524 +/- 0.00249279 (SNR = 443.373, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=17 (freq=1.663e+09 Hz) is: 1.11739 +/- 0.00263747 (SNR = 423.66, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=18 is: INSUFFICIENT DATA 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=19 (freq=1.695e+09 Hz) is: 1.14445 +/- 0.00268405 (SNR = 426.388, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=20 (freq=1.711e+09 Hz) is: 1.15816 +/- 0.00279482 (SNR = 414.394, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=21 (freq=1.727e+09 Hz) is: 1.17194 +/- 0.00285317 (SNR = 410.751, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=22 (freq=1.743e+09 Hz) is: 1.18816 +/- 0.00311174 (SNR = 381.832, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=23 (freq=1.759e+09 Hz) is: 1.2023 +/- 0.00317813 (SNR = 378.303, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=24 (freq=1.775e+09 Hz) is: 1.21832 +/- 0.0032735 (SNR = 372.176, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=25 (freq=1.791e+09 Hz) is: 1.23006 +/- 0.00333594 (SNR = 368.729, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=26 (freq=1.807e+09 Hz) is: 1.24335 +/- 0.00332227 (SNR = 374.246, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=27 (freq=1.823e+09 Hz) is: 1.25561 +/- 0.00305208 (SNR = 411.395, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=28 (freq=1.839e+09 Hz) is: 1.27147 +/- 0.00295938 (SNR = 429.639, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=29 (freq=1.855e+09 Hz) is: 1.28419 +/- 0.00315076 (SNR = 407.581, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=30 (freq=1.871e+09 Hz) is: 1.29938 +/- 0.0031697 (SNR = 409.936, N = 46) 2015-06-08 17:17:05 INFO fluxscale Flux density for J1407+2827 in SpW=31 (freq=1.887e+09 Hz) is: 1.31095 +/- 0.0033746 (SNR = 388.475, N = 46) 2015-06-08 17:17:05 INFO fluxscale Fitted spectrum for J1227+3635 with fitorder=1: Flux density = 1.96397 +/- 0.00183961 (freq=1.56968 GHz) spidx=-0.318786 +/- 0.00698435 2015-06-08 17:17:05 INFO fluxscale Fitted spectrum for J1407+2827 with fitorder=1: Flux density = 1.02742 +/- 0.00099646 (freq=1.56968 GHz) spidx=1.36342 +/- 0.00757383 #Apply calibration, using B0 table for passband, G1 for phase and F1 for fluxes (all flux calibration is now contained in F1, #so no need to use G2 for fluxes, although I'm not sure why you couldn't). Apply first to 3C286, then gain cal, pol cal and source (2,0,3,1). #Don't apply tec.cal yet (although I'm not sure it matters), keep applying it on the fly in polarization calibration steps, then apply in final applycal. applycal(vis='n4244_cl_delay_h.ms',field='2',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1'],gainfield=['','2','2'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') applycal(vis='n4244_cl_delay_h.ms',field='0',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1'],gainfield=['','0','0'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') applycal(vis='n4244_cl_delay_h.ms',field='3',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1'],gainfield=['','3','3'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') applycal(vis='n4244_cl_delay_h.ms',field='1',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1'],gainfield=['','0','0'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') #Flag galaxy as before, before splitting to save time later in polcal. commands not included here. #Split n4244 into own file split(vis='n4244_cl_delay_h.ms',outputvis='n4244clcaled.ms',datacolumn='corrected',field='1') listobs(vis='n4244clcaled.ms') 2015-06-12 22:23:37 INFO listobs Data records: 2954016 Total elapsed time = 20810 seconds 2015-06-12 22:23:37 INFO listobs Observed from 31-Mar-2012/04:23:20.0 to 31-Mar-2012/10:10:10.0 (UTC) 2015-06-12 22:23:44 INFO listobs 2015-06-12 22:23:44 INFO listobs ObservationID = 0 ArrayID = 0 2015-06-12 22:23:44 INFO listobs Date Timerange (UTC) Scan FldId FieldName nRows SpwIds Average Interval(s) ScanIntent 2015-06-12 22:23:44 INFO listobs 31-Mar-2012/04:23:20.0 - 04:30:30.0 18 0 N4244 482976 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-06-12 22:23:44 INFO listobs 04:30:40.0 - 04:40:00.0 19 0 N4244 628992 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-06-12 22:23:44 INFO listobs 07:23:00.0 - 07:27:50.0 45 0 N4244 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-06-12 22:23:44 INFO listobs 07:28:00.0 - 07:36:50.0 46 0 N4244 595296 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-06-12 22:23:44 INFO listobs 09:56:20.0 - 10:01:10.0 70 0 N4244 325728 [0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-06-12 22:23:44 INFO listobs 10:01:20.0 - 10:10:10.0 71 0 N4244 595296 [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31] [10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10, 10] [OBSERVE_TARGET#UNSPECIFIED] 2015-06-12 22:23:44 INFO listobs (nRows = Total number of rows per scan) 2015-06-12 22:23:45 INFO listobs Fields: 1 2015-06-12 22:23:45 INFO listobs ID Code Name RA Decl Epoch SrcId nRows 2015-06-12 22:23:45 INFO listobs 0 NONE N4244 12:17:29.660000 +37.48.25.60001 J2000 0 2954016 2015-06-12 22:23:46 INFO listobs Spectral Windows: (32 unique spectral windows and 1 unique polarization setups) 2015-06-12 22:23:46 INFO listobs SpwID Name #Chans Frame Ch0(MHz) ChanWid(kHz) TotBW(kHz) CtrFreq(MHz) BBC Num Corrs 2015-06-12 22:23:46 INFO listobs 0 Subband:0 64 TOPO 1247.000 250.000 16000.0 1254.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 1 Subband:1 64 TOPO 1263.000 250.000 16000.0 1270.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 2 Subband:2 64 TOPO 1279.000 250.000 16000.0 1286.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 3 Subband:3 64 TOPO 1295.000 250.000 16000.0 1302.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 4 Subband:4 64 TOPO 1311.000 250.000 16000.0 1318.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 5 Subband:5 64 TOPO 1327.000 250.000 16000.0 1334.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 6 Subband:6 64 TOPO 1343.000 250.000 16000.0 1350.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 7 Subband:7 64 TOPO 1359.000 250.000 16000.0 1366.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 8 Subband:8 64 TOPO 1375.000 250.000 16000.0 1382.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 9 Subband:9 64 TOPO 1391.000 250.000 16000.0 1398.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 10 Subband:10 64 TOPO 1407.000 250.000 16000.0 1414.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 11 Subband:11 64 TOPO 1423.000 250.000 16000.0 1430.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 12 Subband:12 64 TOPO 1439.000 250.000 16000.0 1446.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 13 Subband:13 64 TOPO 1455.000 250.000 16000.0 1462.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 14 Subband:14 64 TOPO 1471.000 250.000 16000.0 1478.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 15 Subband:15 64 TOPO 1487.000 250.000 16000.0 1494.8750 12 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 16 Subband:0 64 TOPO 1647.000 250.000 16000.0 1654.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 17 Subband:1 64 TOPO 1663.000 250.000 16000.0 1670.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 18 Subband:2 64 TOPO 1679.000 250.000 16000.0 1686.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 19 Subband:3 64 TOPO 1695.000 250.000 16000.0 1702.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 20 Subband:4 64 TOPO 1711.000 250.000 16000.0 1718.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 21 Subband:5 64 TOPO 1727.000 250.000 16000.0 1734.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 22 Subband:6 64 TOPO 1743.000 250.000 16000.0 1750.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 23 Subband:7 64 TOPO 1759.000 250.000 16000.0 1766.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 24 Subband:8 64 TOPO 1775.000 250.000 16000.0 1782.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 25 Subband:9 64 TOPO 1791.000 250.000 16000.0 1798.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 26 Subband:10 64 TOPO 1807.000 250.000 16000.0 1814.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 27 Subband:11 64 TOPO 1823.000 250.000 16000.0 1830.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 28 Subband:12 64 TOPO 1839.000 250.000 16000.0 1846.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 29 Subband:13 64 TOPO 1855.000 250.000 16000.0 1862.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 30 Subband:14 64 TOPO 1871.000 250.000 16000.0 1878.8750 15 RR RL LR LL 2015-06-12 22:23:46 INFO listobs 31 Subband:15 64 TOPO 1887.000 250.000 16000.0 1894.8750 15 RR RL LR LL #Leaving out mapping and cleaning of Stokes I ######################################################################################################################################################################## #Go on to Pol Cal. Back to unsplit ms. Did no flagging after splitting, and splitting included all galaxy data, so straight to 8.2, flagging RL and LR of pol and flux cals. #If either LR or RL is bad, should flag both, I guess. #1331 LR: ea19 about 3 times as high as all others. RL too Needs flagging in all sources. Commands not included. #Sec 8.4 Set Q and U values for flux cal 3C286, which is used for pol PA and flux. execfile('PolCalSet_Lband.py') #Delay calibration for cross-hand pols. Not sure about the gaintables to apply OTF here. Do these now incorporate antpos and elev gain curve ones? gaincal(vis='n4244_cl_delay_h.ms', solint='60s',caltable='n4244_cl_delay_h.KCROSS',field='2', calmode='p',spw='*:20~45',refant='ea22',gaintype='KCROSS',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','tec.cal'],selectdata=T, parang=T, timerange='05:40:00~05:42:00') #Solve for leakage terms using zero-pol calibrator, correcting for gains, bandpass, cross-term delays, tec, OTF. polcal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.D1',field='3',spw='',refant='ea22',poltype='Df',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.G2','n4244_cl_delay_h.KCROSS','tec.cal'],gainfield=['','3','3','','3']) #Plot cal solutions table plotcal(caltable='n4244_cl_delay_h.D1',xaxis='freq',yaxis='amp',spw='',field='',iteration='antenna') #Leakage a few %, except for SPW 0 ch 35-36, and SPW 29, ch 44-47 on most ants, so flag them all on all sources, LR and RL only flagdata(vis='n4244_cl_delay_h.ms',correlation='LR,RL',spw='0:35~36,29:44~47',field='',antenna='',scan='',timerange='') #Redo polcal and plot polcal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.D1',field='3',spw='',refant='ea22',poltype='Df',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.G2','n4244_cl_delay_h.KCROSS','tec.cal'],gainfield=['','3','3','','3']) plotcal(caltable='n4244_cl_delay_h.D1',xaxis='freq',yaxis='amp',spw='',field='',iteration='antenna') #ea13 has about 25% leakage on pol=0 (this means R, plot just R using poln='R' in plotcal). ea22 (refant) has 0 leakage for pol=0 (R), 2% for pol=1 (L) #Onto pol angle cal with 3C286. Refant still 22, apply gains, bandpass, cross-term delays, tec.cal OTF. #At first couldn't find soln for SPW 0. As per cookbook, it then fails on all other SPWs. SPW 0 only has a few non-zero channels in cals, and is flagged out for source #anyway, so flag it in LR, RL for all sources. #Redo polcal starting with SPW 3, which is now the first SPW with pol data. polcal(vis='n4244_cl_delay_h.ms',caltable='n4244_cl_delay_h.X1',field='2',refant='ea22',spw='3~31',poltype='Xf',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.G2','n4244_cl_delay_h.D1','n4244_cl_delay_h.KCROSS','tec.cal'],gainfield=['', '2','2','','',''],smodel=[1.0,0.407,0.914,0.0]) plotcal(caltable='n4244_cl_delay_h.X1',xaxis='freq',yaxis='phase',spw='',field='',iteration='antenna') #Now apply pol angle (X) and leakage calibrations (D) to all sources. Reminder: 3C286, gain cal, pol cal and source are (2,0,3,1). #The first time I did this I forgot to apply tec.cal. I then copied the .ms to a new copy, n4244_cl_delay_h2.ms, and proceeded, hence #the change in the ms name below. #Flux calibrator applycal(vis='n4244_cl_delay_h2.ms',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D1','n4244_cl_delay_h.X1','tec.cal'],parang=True,field='2',gainfield=['','2','2','','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Zero pol cal applycal(vis='n4244_cl_delay_h2.ms',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D1','n4244_cl_delay_h.X1','tec.cal'],parang=True,field='3',gainfield=['','3','3','','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Gain cal applycal(vis='n4244_cl_delay_h2.ms',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D1','n4244_cl_delay_h.X1','tec.cal'],parang=True,field='0',gainfield=['','0','0','','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Source applycal(vis='n4244_cl_delay_h2.ms',gaintable=['n4244_cl_delay_h.B0','n4244_cl_delay_h.G1','n4244_cl_delay_h.F1','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D1','n4244_cl_delay_h.X1','tec.cal'],parang=True,field='1',gainfield=['','0','0','','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Split off source again split(vis='n4244_cl_delay_h2.ms',outputvis='n4244_cl_delay_h2_IQUVcaled.ms',antenna='',datacolumn='corrected',field='1',scan='',spw='') ################################################################################################################################################################### ################################################################################################################################################################### ################################################################################################################################################################### #Redo of polcal and mapping without tec.cal at all. Create new copy of n4244_cl_delay_h2.ms called n4244_cl_delay_h3.ms. Should have all flagging done. Redo pol angle and intensity #cal. All other cal files the same. #Initial gain calibration for passband calibration, but now with refant 22 gaincal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G3',field='2',refant='ea22',spw='0~31:15~50',calmode='p',solint='int',minsnr=5,solnorm=T) #Reports all attempted calibrations succeeded ######################## #Bandpass cal. refant 22, tec.cal, G3 gain cal applied OTF. gaincurve parameter no longer exists. This was the elevation gain correction curve. Presumably #applied automatically now. The casa webpages sort of imply that we should have been doing this for VLA data anyway when the parameter existed. bandpass(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.B1',field='2',spw='',refant='ea22',solnorm=False,combine='scan',solint='inf',bandtype='B',gaintable=['n4244_cl_delay_h.G3']) ################################ #Gain calibration. Use refant 22, create version G4, not G3 as last time. gaincurve doesn't exist anymore #Phases first. First for fluxcal. Do solns for 10 sec integrations for some reason, apparently because these phase solns will only get used for calibrating the amp # of the flux cal later. gaincal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G4',solint='int',field='2',spw='0~31:5~58',refant='ea22',gaintype='G',calmode='p',solnorm=F,gaintable=['n4244_cl_delay_h.B1']) #For other cals, solint=inf and combine='' apparently means do a solution for each scan. Append=T means appending onto previous table gaincal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G4',solint='inf',field='0,3',spw='0~31:5~58',refant='ea22',gaintype='G',calmode='p',solnorm=F,append=T,gaintable=['n4244_cl_delay_h.B1']) #Now amp cal on a per scan basis for each cal in turn. Gainfield param means use 3c286 from passband gaintable and the phases for the field for the #source being amp calibrated from the gain table G4. Applying B1, G4, OTF, refant 22, creating G5. gaincal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G5',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4'],gainfield=['2','2',''],field='2',refant='ea22',solnorm=F,spw='0~31:5~58',solint='inf',gaintype='G',calmode='a') gaincal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G5',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4'],gainfield=['2','0',''],field='0',refant='ea22',solnorm=F,spw='0~31:5~58',solint='inf',gaintype='G',calmode='a',append=T) gaincal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G5',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4'],gainfield=['2','3',''],field='3',refant='ea22',solnorm=F,spw='0~31:5~58',solint='inf',gaintype='G',calmode='a',append=T) # Find other cal fluxes from 3C286. For 1227, we get 1.8-2.1 Jy. It is 2.1 at 20cm according to cal manual. For 1407, we get 0.8-1.3 Jy. Cal manual has 0.8 Jy. fluxscale(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.G5',fluxtable='n4244_cl_delay_h.F2',reference=['2'],transfer=['0,3']) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=3 (freq=1.295e+09 Hz) is: 2.08647 +/- 0.0068149 (SNR = 306.164, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=4 (freq=1.311e+09 Hz) is: 2.07252 +/- 0.0035668 (SNR = 581.06, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=5 (freq=1.327e+09 Hz) is: 2.06838 +/- 0.0032996 (SNR = 626.857, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=6 (freq=1.343e+09 Hz) is: 2.06145 +/- 0.00312041 (SNR = 660.635, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=7 (freq=1.359e+09 Hz) is: 2.05015 +/- 0.00330129 (SNR = 621.016, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=8 (freq=1.375e+09 Hz) is: 2.04736 +/- 0.00320794 (SNR = 638.216, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=9 (freq=1.391e+09 Hz) is: 2.05032 +/- 0.00305226 (SNR = 671.739, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=10 (freq=1.407e+09 Hz) is: 2.03794 +/- 0.00326591 (SNR = 624.004, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=11 (freq=1.423e+09 Hz) is: 2.03285 +/- 0.00300165 (SNR = 677.244, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=12 (freq=1.439e+09 Hz) is: 2.02264 +/- 0.00331652 (SNR = 609.867, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=13 (freq=1.455e+09 Hz) is: 2.02357 +/- 0.003429 (SNR = 590.135, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=14 (freq=1.471e+09 Hz) is: 2.01302 +/- 0.00357394 (SNR = 563.249, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=15 (freq=1.487e+09 Hz) is: 2.01056 +/- 0.00434846 (SNR = 462.362, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=16 (freq=1.647e+09 Hz) is: 1.95034 +/- 0.0020876 (SNR = 934.249, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=17 (freq=1.663e+09 Hz) is: 1.95456 +/- 0.00205443 (SNR = 951.389, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=18 is: INSUFFICIENT DATA 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=19 (freq=1.695e+09 Hz) is: 1.93203 +/- 0.00237006 (SNR = 815.181, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=20 (freq=1.711e+09 Hz) is: 1.92189 +/- 0.00174893 (SNR = 1098.89, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=21 (freq=1.727e+09 Hz) is: 1.91279 +/- 0.00166655 (SNR = 1147.75, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=22 (freq=1.743e+09 Hz) is: 1.90735 +/- 0.0018544 (SNR = 1028.56, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=23 (freq=1.759e+09 Hz) is: 1.89663 +/- 0.00201785 (SNR = 939.926, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=24 (freq=1.775e+09 Hz) is: 1.8895 +/- 0.00228282 (SNR = 827.702, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=25 (freq=1.791e+09 Hz) is: 1.87947 +/- 0.00278504 (SNR = 674.844, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=26 (freq=1.807e+09 Hz) is: 1.87491 +/- 0.00249057 (SNR = 752.804, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=27 (freq=1.823e+09 Hz) is: 1.86244 +/- 0.00263568 (SNR = 706.625, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=28 (freq=1.839e+09 Hz) is: 1.85738 +/- 0.00252576 (SNR = 735.373, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=29 (freq=1.855e+09 Hz) is: 1.85152 +/- 0.00251786 (SNR = 735.353, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=30 (freq=1.871e+09 Hz) is: 1.84432 +/- 0.00259067 (SNR = 711.91, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1227+3635 in SpW=31 (freq=1.887e+09 Hz) is: 1.83392 +/- 0.00275262 (SNR = 666.246, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=0 is: INSUFFICIENT DATA 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=1 is: INSUFFICIENT DATA 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=2 is: INSUFFICIENT DATA 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=3 (freq=1.295e+09 Hz) is: 0.781161 +/- 0.00290901 (SNR = 268.531, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=4 (freq=1.311e+09 Hz) is: 0.796415 +/- 0.00228537 (SNR = 348.484, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=5 (freq=1.327e+09 Hz) is: 0.813593 +/- 0.00215825 (SNR = 376.969, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=6 (freq=1.343e+09 Hz) is: 0.825898 +/- 0.00234179 (SNR = 352.678, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=7 (freq=1.359e+09 Hz) is: 0.838976 +/- 0.00237412 (SNR = 353.384, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=8 (freq=1.375e+09 Hz) is: 0.858414 +/- 0.00300977 (SNR = 285.21, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=9 (freq=1.391e+09 Hz) is: 0.872754 +/- 0.00290581 (SNR = 300.348, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=10 (freq=1.407e+09 Hz) is: 0.891103 +/- 0.00291792 (SNR = 305.39, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=11 (freq=1.423e+09 Hz) is: 0.90785 +/- 0.00283524 (SNR = 320.202, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=12 (freq=1.439e+09 Hz) is: 0.91586 +/- 0.00218548 (SNR = 419.065, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=13 (freq=1.455e+09 Hz) is: 0.93013 +/- 0.00200335 (SNR = 464.288, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=14 (freq=1.471e+09 Hz) is: 0.94385 +/- 0.00240319 (SNR = 392.749, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=15 (freq=1.487e+09 Hz) is: 0.957518 +/- 0.00303328 (SNR = 315.67, N = 42) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=16 (freq=1.647e+09 Hz) is: 1.10631 +/- 0.00214951 (SNR = 514.681, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=17 (freq=1.663e+09 Hz) is: 1.11997 +/- 0.00227167 (SNR = 493.015, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=18 is: INSUFFICIENT DATA 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=19 (freq=1.695e+09 Hz) is: 1.14539 +/- 0.00202478 (SNR = 565.684, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=20 (freq=1.711e+09 Hz) is: 1.15828 +/- 0.0020037 (SNR = 578.071, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=21 (freq=1.727e+09 Hz) is: 1.17301 +/- 0.00196566 (SNR = 596.751, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=22 (freq=1.743e+09 Hz) is: 1.18927 +/- 0.00217682 (SNR = 546.334, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=23 (freq=1.759e+09 Hz) is: 1.20426 +/- 0.00223462 (SNR = 538.911, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=24 (freq=1.775e+09 Hz) is: 1.21957 +/- 0.00230952 (SNR = 528.063, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=25 (freq=1.791e+09 Hz) is: 1.23063 +/- 0.00256705 (SNR = 479.396, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=26 (freq=1.807e+09 Hz) is: 1.24448 +/- 0.00249185 (SNR = 499.42, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=27 (freq=1.823e+09 Hz) is: 1.25597 +/- 0.00241451 (SNR = 520.175, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=28 (freq=1.839e+09 Hz) is: 1.27184 +/- 0.00234057 (SNR = 543.391, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=29 (freq=1.855e+09 Hz) is: 1.28385 +/- 0.00236521 (SNR = 542.805, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=30 (freq=1.871e+09 Hz) is: 1.29972 +/- 0.00242767 (SNR = 535.377, N = 44) 2015-10-07 19:18:17 INFO fluxscale Flux density for J1407+2827 in SpW=31 (freq=1.887e+09 Hz) is: 1.31243 +/- 0.00264505 (SNR = 496.184, N = 44) 2015-10-07 19:18:17 INFO fluxscale Fitted spectrum for J1227+3635 with fitorder=1: Flux density = 1.9631 +/- 0.00216575 (freq=1.58264 GHz) spidx=-0.325417 +/- 0.00940658 2015-10-07 19:18:17 INFO fluxscale Fitted spectrum for J1407+2827 with fitorder=1: Flux density = 1.03987 +/- 0.00104437 (freq=1.58264 GHz) spidx=1.36215 +/- 0.00820404 #Apply calibration, using B1 table for passband, G4 for phase and F2 for fluxes (all flux calibration is now contained in F2, #so no need to use G5 for fluxes, although I'm not sure why you couldn't). Apply first to 3C286, then gain cal, pol cal and source (2,0,3,1). applycal(vis='n4244_cl_delay_h3.ms',field='2',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2'],gainfield=['','2','2'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') applycal(vis='n4244_cl_delay_h3.ms',field='0',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2'],gainfield=['','0','0'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') applycal(vis='n4244_cl_delay_h3.ms',field='3',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2'],gainfield=['','3','3'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') applycal(vis='n4244_cl_delay_h3.ms',field='1',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2'],gainfield=['','0','0'],interp=['nearest','','nearest',''],parang=F,calwt=F,spw='0~31:5~58') ######################################################################################################################################################################## #Go on to Pol Cal. Back to unsplit ms. Not repeating delay cal for cross-hand pols, assuming that is not affected by tec. #Solve for leakage terms using zero-pol calibrator, correcting for gains, bandpass, cross-term delays, OTF. polcal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.D3',field='3',spw='',refant='ea22',poltype='Df',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.G5','n4244_cl_delay_h.KCROSS'],gainfield=['','3','3','','3']) #Onto pol angle cal with 3C286. Refant still 22, apply gains, bandpass, cross-term delays, OTF. polcal(vis='n4244_cl_delay_h3.ms',caltable='n4244_cl_delay_h.X3',field='2',refant='ea22',spw='3~31',poltype='Xf',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.G5','n4244_cl_delay_h.D3','n4244_cl_delay_h.KCROSS'],gainfield=['', '2','2','','',''],smodel=[1.0,0.407,0.914,0.0]) #Now apply pol angle (X) and leakage calibrations (D) to all sources. Reminder: 3C286, gain cal, pol cal and source are (2,0,3,1). #Flux calibrator applycal(vis='n4244_cl_delay_h3.ms',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D3','n4244_cl_delay_h.X3'],parang=True,field='2',gainfield=['','2','2','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Zero pol cal applycal(vis='n4244_cl_delay_h3.ms',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D3','n4244_cl_delay_h.X3'],parang=True,field='3',gainfield=['','3','3','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Gain cal applycal(vis='n4244_cl_delay_h3.ms',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D3','n4244_cl_delay_h.X3'],parang=True,field='0',gainfield=['','0','0','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Source applycal(vis='n4244_cl_delay_h3.ms',gaintable=['n4244_cl_delay_h.B1','n4244_cl_delay_h.G4','n4244_cl_delay_h.F2','n4244_cl_delay_h.KCROSS','n4244_cl_delay_h.D3','n4244_cl_delay_h.X3'],parang=True,field='1',gainfield=['','0','0','','',''], spw='',interp=['','','nearest','nearest','',''],calwt=F) #Split off source again split(vis='n4244_cl_delay_h3.ms',outputvis='n4244_cl_delay_h3_IQUVcaled.ms',antenna='',datacolumn='corrected',field='1',scan='',spw='')