From dshepher@aoc.nrao.edu Mon Mar 5 10:11:08 2001 Date: Fri, 2 Mar 2001 17:45:35 -0700 (MST) From: Debra Shepherd To: Steven Myers Cc: Ed Fomalont , Greg Taylor , John Benson , John Hibbard Subject: aips script for qband 3c84 data (fwd) Steve, this message and the next one are detailed scripts in aips and aips++ which reduce the same data in as similar a fashion as possible. Although the aips++ script is not ready for outside users, it might be useful for our testing group to have the two scripts available for comparison. Perhaps you could put them on the web site? Thanks, Debra ---------- Forwarded message ---------- Date: Fri, 2 Mar 2001 17:36:19 -0700 (MST) From: Debra Shepherd To: Tim Cornwell , George Moellenbrock Cc: Debra Shepherd Subject: aips script for qband 3c84 data 22dec00 VLA Q and K band tests on a bright quasar program number DSTST Q band cont VLA A array synth beam Q: 0.05'' Theta(LAS) : 1.3'' Primary beam: 60'' ** Visibility Tape Information ** Tape # File # Time of Final Record: N13920 1 restore 0 uc tape; inp mount; intape 1; remhost ''; tape; mount task 'fillm' vlaobs 'DSTST'; doall 1; ncount 1; nfiles 0; doweight 1; cparm 0 dparm 0; douvcomp -1; doconcat -1; outdi 2; outna '' cparm(6) 1; inp go fillm - cparm(6) specifies subarray... subarray 1 = data - default averaging time (no averaging). ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** ** spy dismount uc AIPS 1: Catalog on disk 2 AIPS 1: 59 2160 20001222 .K BAND. 1 UV 24-JAN-2001 11:48:42 AIPS 1: 60 2160 20001222 .Q BAND. 1 UV 24-JAN-2001 11:51:02 src = '''' (3c84) gcal = '0319+415''' (3c84) 5s scans, no fast switching, just sit and stare at 3c84 for ~1 hr at Q and ~1 hr at K ******* ******* ******* K & Q-band continuum in subarray 1 ******* ******* ******* ******************************************************************* getn 59 task 'listr' indi 2; getn 59; tput calib; optyp 'scan'; inext ''; factor 1; bif 1 eif 1; docal -1; docrt 0; sourc ''; outpri 'klistr.scan';inp go listr - output in klistr.scan: task 'listr' indi 2; getn 60; tput calib; optyp 'scan'; inext ''; factor 1; bif 1 eif 1; docal -1; docrt 0; sourc ''; outpri 'qlistr.scan';inp go listr - output in qlistr.scan: task 'prtan' ncount 0; invers 1; docrt 0; outpri 'prtan.out'; inp go prtan - output in prtan.out - Non Q-band rxs: - totally out of the array: 5 & 29, - 15 & 19 had servo problems, probably delete all *********************** - see if weights were filled during FILLM: There are several routines that allow you to examine the weights: 1. PRTUV to look at the weights directly. 2. UVPLT to plot the weights versus elevation, time, etc... with and without DOCAL = 2 to examine the effect. 3. VPLOT to plot the weights for each baseline on separate pages. task 'prtuv' sour '0319+415'''; xinc 1; cparm 0; dparm 0; docrt 132; outpr ''; inp - OK, weights were filled. - see if we need to quack scans and get gen'l idea of data quality: task 'listr' indi 2; getn 60 optyp 'list'; anten 8,0; basel 0; sour ''; calcode ''; stokes 'rr' bif 1; docrt 132; outpr ''; docal -1; dparm 0; inp go listr - Don't quack for now. - AN 17 has very low gains at Q at first, then gains creep up... - AN 8 looks like a reasonable refant indi 2; getn 59; inp go listr - AN 17 gains at K good at first, then then drop suddenly. - AN 10 has some occasional drop outs -- edit carefully - AN 8 looks like a reasonable refant ******************************************************************* - plot time vs elevation for srcs: task 'uvplt' indi 2; getn 60; anten 1 0; xinc 10; stokes 'rr' bparm 5 15 0; refant 8; dotv -1; inp go uvplt plver 1; outfi 'q.el.tm.ps' go lwpla - src rising fast at Q band indi 2; getn 59; inp go uvplt plver 1; outfi 'k.el.tm.ps' go lwpla - src roughly uniform elevation at K band **** Didn't observe flux cal. **** manually insert flux density into stokes I of SU table: 0319+415 J2000 A 03h19m48.1601s 41d30'42.106" 3C84 0316+413 B1950 A 03h16m29.5673s 41d19'51.919" NGC1275 ----------------------------------------------------- BAND A B C D FLUX(Jy) UVMIN(kL) UVMAX(kL) ===================================================== 90cm P S X X X 8 13 20cm L P P X X 23.9 12 6cm C P P P P 23.3 3.7cm X P P P P 21.70 visplot 2cm U P P P P 20.70 visplot 1.3cm K X S S S 16.4 0.7cm Q X S S S 9.00 1800 visplot task 'setjy' indi 2; getn 60; sour '0319+415''' opty ''; zerosp 9.0 0 0 0 0 freqid = 1; bif 1; eif 2; aparm = 0; xinc 1; inp go setjy - SETJY1: SOURCES = '0319+415 ' / Included - SETJY1: ZEROSP = 9.000 0.000 0.000 0.000 /I,Q,U,V fluxes task 'prtab' indi 2; getn 60; docrt 132; inex 'su'; invers 1; inp go prtab task 'setjy' indi 2; getn 59; sour '0319+415''' opty ''; zerosp 16.4 0 0 0 0 freqid = 1; bif 1; eif 2; aparm = 0; xinc 1; inp go setjy - SETJY1: SOURCES = '0319+415 ' / Included - SETJY1: ZEROSP = 16.400 0.000 0.000 0.000 /I,Q,U,V fluxes task 'prtab' indi 2; getn 59; docrt 132; inex 'su'; invers 1; inp go prtab task 'tvflg' indi 2; getn 59; sour '0319+415''' stokes ''; bif 1; eif 2; docat 0; anten 0; calcode '' docal -1; flagver 0; dparm(6) 5; inp go tvflg - rr1 rr2 ll2 ll1 - K band flagged: - 8 partial dropouts on AN 10 - 2nd 1/2 of AN 17 - A few dropouts on AN 20 indi 2; getn 60; inp go tvflg - rr1 rr2 ll2 ll1 - Q band flagged: - 1st 1/2 of AN 17 - occasional dropouts on 10,14,22 are flagged online - calib solns show some bad 1st scans. quack first 5s of each scan: task 'quack' indi 2; getn 59; sour '0319+415'''; timer 0; antennas 0; flagver 1 opcode 'beg'; reason 'bad start of scan'; aparm 0, 5/60,0; inp go quack - QUACK1: Flagged 0.33 min. of data in 2 scans indi 2; getn 60; inp go quack - QUACK1: Flagged 0.33 min. of data in 2 scans * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * solve for phase and amp solns * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * src = '0319+415''' (3c84) gcal = '0319+415''' task 'calib' indi 2; getn 59 refant 8; soltype 'l1'; solmod 'a&p'; solint 5/60; flagver 0; docal 2 cparm 0 0 20 20 1 0; aparm 0; aparm(6) 3; outna ''; outdi 2; timerang 0 calsour '0319+415'''; snver 1; inp go calib - one solution for each 5 sec scan to see rx stability: - snr = 500 - 1500 - CALIB1: Found 46388 good solutions - CALIB1: Failed on 12 solutions - CALIB1: 16 solutions had insufficient data tvall task 'snplt' xinc 1; bif 2; eif 2; stokes 'rrll'; sour ''; timerang 0 opty 'amp'; opco 'alif'; inex 'sn'; nplots 6; dotv 1 invers 1; inp go snplt - ants 10 & 20 still having minor 5s dropouts I didn't find with tvflg - ant 27 amps decrease slowly with time by 1/3 then increase - *many* ants show encoder oscillation problems dotv -1; go plver 2; outfi 'k.gain1.ps' go lwpla plver 3; outfi 'k.gain2.ps' go lwpla plver 4; outfi 'k.gain3.ps' go lwpla plver 5; outfi 'k.gain4.ps' go lwpla plver 6; outfi 'k.gain5.ps' go lwpla plver 7; outfi 'k.gain6.ps' go lwpla plver 8; outfi 'k.gain7.ps' go lwpla plver 9; outfi 'k.gain8.ps' go lwpla plver 10; outfi 'k.gain9.ps' go lwpla opty 'phas'; dotv 1; go - phase tracks OK. task 'calib' indi 2; getn 60 refant 8; soltype 'l1'; solmod 'a&p'; solint 5/60; flagver 0; docal 2 cparm 0 0 20 20 1 0; aparm 0; aparm(6) 3; outna ''; outdi 2; timerang 0 calsour '0319+415'''; snver 1; inp go calib - one solution for each 5 sec scan to see rx stability: - snr = 100 - 300, adequate - CALIB1: Found 43500 good solutions - CALIB1: 8 solutions had insufficient data task 'snplt' indi 2; getn 60; xinc 1; bif 2; eif 2; stokes 'rrll'; sour ''; timerang 0 opty 'amp'; opco 'alif'; inex 'sn'; nplots 6; dotv 1 invers 1; inp go snplt - *many* ants show encoder oscillation problems dotv -1; go plver 2; outfi 'q.gain1.ps' go lwpla plver 3; outfi 'q.gain2.ps' go lwpla plver 4; outfi 'q.gain3.ps' go lwpla plver 5; outfi 'q.gain4.ps' go lwpla plver 6; outfi 'q.gain5.ps' go lwpla plver 7; outfi 'q.gain6.ps' go lwpla plver 8; outfi 'q.gain7.ps' go lwpla plver 9; outfi 'q.gain8.ps' go lwpla opty 'phas'; dotv 1; go - phases track OK, no problems - don't need to do getjy, since I setjy the exact value. - create CL2 table: CL1 has amps = 1 and phase= 0 for all time stamps. take the SN gains (amp & phase) and interpolate between cal points to make a table with all time stamps that can be applied to the src. Use linear interpolation, no smoothing (2PT). Combine all SN tables into a CL table (CALI). task 'clcal' indi 2 getn 59 interp '2pt'; refant 8; qual -1; opco 'cali'; gainver 1; gainuse 2 snver 1; calsour '0319+415'''; sour '0319+415'''; inp go clcal - interpolate on 3c84 - CLCAL1: SN2CL: Applying SN tables to CL table 1, writing CL table 2 indi 2 getn 60; inp go clcal - CLCAL1: SN2CL: Applying SN tables to CL table 1, writing CL table 2 - look at antenna-based gain solutions: task 'listr' indi 2 getn 59 sour ''; opty 'gain'; inex 'cl'; inver 2; docal 2; gainuse 2; flagver 1 docrt 132; outpr ''; dparm 0; inp go listr - looks reasonable indi 2 getn 60; inp go listr - looks reasonable for Q band - look at scalar scan averages of amp and phase The matrix average amplitudes for the calibrators in this listing should be very close to the values that were entered with SETJY (or which were derived by GETJY) and all the phases in all rows and columns should be close to zero. indi 2 getn 59 opty 'matx'; dparm 5 1 0; stokes ''; sour '0319+415'''; inp go listr - K residual phases generally less than 6 deg - AMP Matrix average = 1.64174E+01 sigma = 8.67675E-03 - PHA Matrix average = -3.07788E-01 sigma = 1.35761E-01 indi 2 getn 60; inp opty 'matx'; dparm 5 1 0; stokes ''; sour '0319+415'''; inp go listr - Q residual phases generally less than 15 deg - AMP Matrix average = 9.08537E+00 sigma = 1.78546E-02 - PHA Matrix average = 3.12542E-01 sigma = 2.07259E-01 - split out calibrated source data: task 'split' indi 2; getn 59; sour '0319+415'''; stokes ''; bif 1; eif 2; docal 2; gainuse 2; flagver 1; outcl 'ksplit'; outdi 2; subarray 0 douvcomp -1; aparm 0; inp go split - SPLIT1: Create 0319+415 .KSPLIT. 1 (UV) on disk 2 cno 61 indi 2; getn 60; sour '0319+415'''; stokes ''; bif 1; eif 2; docal 2; gainuse 2; flagver 1; outcl 'qsplit'; outdi 2; subarray 0 douvcomp -1; aparm 0; inp go split - SPLIT1: Create 0319+415 .QSPLIT. 1 (UV) on disk 2 cno 62 uc AIPS 1: Catalog on disk 2 AIPS 1: 61 2160 0319+415 .KSPLIT. 1 UV 24-JAN-2001 16:54:26 AIPS 1: 62 2160 0319+415 .QSPLIT. 1 UV 24-JAN-2001 16:55:25 - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - image details: VLA A array synth beam Q: 0.05'' K: 0.1'' Primary beam: 60'' 120'' - Q band theory rms approx. 0.11 mJy/beam = 0.00011 Jy/bm - K band theory rms approx. 0.14 mJy/beam = 0.00014 Jy/bm - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - Cannot get anywhere close to theoretical rms with straight clean. Dynamic range limited. Best I can do without selfcal: tvall Q task 'imagr' indi 2; getn 62; docal -1; sour '0319+415'''; outna '3c84q'; outdi 2 cell = 0.01; imsize = 1024; minpatch 51; gain 0.1; clbox 0; cparm 0; dotv 1; xtype 5; ytype 5; uvwtfn 'na'; flux 0.00011; niter 1000000; zerosp 0; inp go imagr - cell = 0.01; imsize = 1024 niter 1000000: - IMAGR1: Field 1 fit FWHM = 62.419 x 43.653 Milliarcsec, PA= -32.0 - IMAGR1: Reached min. Clean flux density= 109.961 MicroJy iter= 823316 - IMAGR1: Total Cleaned flux density = 8.921 Jy 823316 comps getn 65 tvbox blc 500 500; trc 525 525; imstat - AIPS 1: Maximum= 8.7809E+00 at 512 513 1 1 1 1 1 tvbox blc 467.00 461.00; trc 557.00 501.00; imstat - AIPS 1: Mean= 9.2047E-03 rms= 2.3293E-02 JY/BEAM over 3731. pixels tvall K task 'imagr' indi 2; getn 61; docal -1; sour '0319+415'''; outna '3c84k'; outdi 2 cell = 0.02; imsize = 1024; minpatch 51; gain 0.1; clbox 0; cparm 0; dotv 1; xtype 5; ytype 5; uvwtfn 'na'; flux 0.00014; niter 100000; zerosp 0; inp go imagr - IMAGR1: Field 1 fit FWHM = 110.222 x 96.442 Milliarcsec, PA= -24.3 - IMAGR1: Reached Iter. limit, Max resid=-632.473 MicroJy iter= 100000 - IMAGR1: Total Cleaned flux density = 16.385 Jy 100000 comps getn 67 tvbox blc 500 500; trc 525 525; imstat - AIPS 1: Maximum= 1.6261E+01 at 512 513 1 1 1 1 1 tvbox blc 467.00 461.00; trc 557.00 501.00; imstat - AIPS 1: Mean= 6.8843E-04 rms= 2.2115E-02 JY/BEAM over 3731. pixels tvinit task 'kntr' indi 2; getn xx; dogrey -1; blc 950 950; trc 1100 1100; clev xx; levs -3,3,4,5,6,7,8,9,11,13,15,17,19; dotv -1; inp go kntr plver 1; outfi 'xx.ps'; go lwpla