[daip] New Ticket - [!NTC-279643]: EVLA HI Spectral Line Calibration in AIPS
John Michael Cannon
do-not-reply at nrao.edu
Thu May 31 13:13:20 EDT 2012
New Ticket: EVLA HI Spectral Line Calibration in AIPS
Hello,
I am using AIPS to calibrate a number of HI spectral line datasets acquired with the EVLA. One of these datasets (program 12A-074) is being worked on by a student of mine; these data contain 16 spectral windows. During our discussions we have been forced to think very carefully about the "standard" treatments in AIPS. Specifically I am referring to Appendix E of the AIPS cookbook and the discussion therein.
I include below my current reduction pipeline, which seems to work well for most simple datasets (single, 1 MHz wide datasets from AC990). A few specific items we have noticed are the following.
1. We request the AIPS UVFITS format from the NRAO archive, meaning that certain tables (perhaps most noticeably SY) are not present. Is it recommended to use the SDMBDF format & then use OBIT?
2. Section E.6 discusses "Calibration with the SysPower" table. Specifically, one might use TYAPL to supply the SY table. I have been operating under the assumption that using REWAY instead of SPLIT to apply calibration and to create single-source datasets does this correctly. Stated differently, I have been assuming that the SY table is not needed and that REWAY does the weighting correctly. Is this the case, or should I revise my strategy?
3. My pipeline does not use RFLAG. I flag on the primary, derive a FRING solution and apply it, and then flag in detail on every source with both TVFLG and SPFLG. I do this on uncalibrated data. I then derive a calibration solution, flag further on the single source dataset(s) if needed, and he use REWAY to split off single source datasets. This differs significantly from the pipeline discussed in Appendix E: initial flagging, calibration, RFLAG, calibration with SysPower table, application of SY table, repeat of calibration steps, then SPLIT or REWAY.
Any and all suggestions or comments on these questions and on the pipeline below are welcome. In the end I need to apply this process uniformly to 30+ datasets, so it is important that all potential mis-steps are identified early on.
Thanks very much,
John
1. UVLOD: load data
default 'uvlod'; bcount 0;
ncount 0;
datain 'VLADATA:AC990.AGC731457.D'; douvcomp -1;
outna 'AGC731457_D'; outcl 'UVDATA'; outseq 0;
outdi 1; doall 1; npoints 0; go
IMHEAD:
256 channels, 1.4183510E+09 Hz at channel 129.00, 3.9062500E+03 Hz/channel
2. INDXR
task 'indxr'; getn 1;
infile ''; prtlev 0; cparm 0,0,0.06666667;
bparm 0,0; go
--> NOTE: using CPARM(3)=4/60, to match 4-second scan intervals
3. LISTR
default 'listr'; getn 1;
optype 'scan'; docrt 132;
go
File = AGC731457_D .UVDATA. 1 Vol = 1 Userid = 9309
Freq = 1.418351008 GHz Ncor = 2 No. vis = 575575
Scan summary listing
Scan Source Qual Calcode Sub Timerange FrqID START VIS END VIS
1 1331+305=3C286 : 0000 1 0/13:25:18 - 0/13:26:10 1 1 4550
2 1331+305=3C286 : 0000 K 1 0/13:26:18 - 0/13:34:10 1 4551 43225
3 J1021+2159 : 0000 D 1 0/13:34:18 - 0/13:40:10 1 43226 71825
4 AGC731457 : 0000 1 0/13:40:18 - 0/13:59:06 1 71826 163150
5 J1021+2159 : 0000 D 1 0/13:59:14 - 0/14:01:06 1 163151 172575
6 AGC731457 : 0000 1 0/14:01:14 - 0/14:20:02 1 172576 263900
7 J1021+2159 : 0000 D 1 0/14:20:10 - 0/14:22:02 1 263901 273325
8 AGC731457 : 0000 1 0/14:22:10 - 0/14:40:58 1 273326 364650
9 J1021+2159 : 0000 D 1 0/14:41:06 - 0/14:42:58 1 364651 374075
10 AGC731457 : 0000 1 0/14:43:06 - 0/15:01:54 1 374076 465400
11 J1021+2159 : 0000 D 1 0/15:02:02 - 0/15:03:54 1 465401 474825
12 AGC731457 : 0000 1 0/15:04:02 - 0/15:22:54 1 474826 566475
13 J1021+2159 : 0000 D 1 0/15:23:02 - 0/15:24:50 1 566476 575575
Source summary
Velocity type = ' ' Definition = ' '
ID Source Qual Calcode RA(2000.0) Dec(2000.0) IFlux QFlux UFlux VFlux No. vis
1 1331+305=3C286 : 0000 13:31:08.2880 30:30:32.959 0.000 0.000 0.000 0.000 4550
2 1331+305=3C286 : 0000 K 13:31:08.2880 30:30:32.959 0.000 0.000 0.000 0.000 38675
3 J1021+2159 : 0000 D 10:21:54.5289 21:59:30.440 0.000 0.000 0.000 0.000 75400
4 AGC731457 : 0000 10:31:55.8000 28:01:33.000 0.000 0.000 0.000 0.000 456950
ID Source Freq(GHz) Velocity(Km/s) Rest freq (GHz)
1 All Sources 1.4184 454.0000 1.4204
Frequency Table summary
FQID IF# Freq(GHz) BW(kHz) Ch.Sep(kHz) Sideband
1 1 1.41835101 1000.0001 3.9063 1
4. PTRAN
task 'prtan'; optype 'scan'; docrt 132; go
Location of VLA antennas
N09 (20)*
N08 (28)*
N07 ( 8)*
N06 (23)*
N05 (10)*
N04 (16)*
( )
N02 (22)*
N01 (13)*
*(24) W01 E01 (21)*
*(17) W02 E02 ( 9)*
*(14) W03 E03 (11)*
*(19) W04 E04 ( 7)*
*( 4) W05 E05 (18)*
*(26) W06 E06 ( 2)*
*( 1) W07 E07 (27)*
*( 5) W08 E08 (12)*
*(15) W09 E09 ( 3)*
( 6)
(25)
* => EVLA ANTENNA
5. VLANT
default 'vlant'; getn 1; go
-No changes needed
- 12 antennas may still get more corrections
6. From log file:
************************************************************************
**** BEGIN EXAMINING DATA NOW ****
****
**** EVLA OPERATOR LOG COMMENTS FOLLOW
****
**** NO DOWNTIME.
****
**** SHAOWING OCCURS AT BEGINNING OF THE SB
****
**** Ant 6: all data, elevation motor problems. Data lost.
****
**** Ant 24: 13:25-15 - 13:51:00 - Cryo issues. Data lost.
****
**** Ant 13: 13:37:00 - 13:52:00 - cryo issues. Data lost.
****
**** Ant 21: 14:14:00 - 15:24:58 - ant in local control. Data lost.
****
**** ANT 1 recently moved, check for baseline errors.
****
************************************************************************
7. Examine PRIMARY CAL data in TVFLG:
-- ASSURE THAT YOU SET 'SWITCH ALL-CH FLAG'!!
**** Primary calibrator, RR:
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1;
dohist -1; sources ''; calcode 'K';
timerang 0; stokes 'RR'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- 10 flags issued
**** Primary calibrator, LL:
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1;
dohist -1; sources ''; calcode 'K';
timerang 0; stokes 'LL'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- 10 flags issued
8. FRING
-Selected 00:13:31:30-00:13:32:30 on 1331+305=3C286
task 'fring'; default fring;
getn 1; timerang 00 13 31 30 00 13 32 30;
calsour '1331+305=3C286'; calcode 'K';
solint 1.05; dparm(9)=1;
--> FINDS 42 good solutions
--> The SN table needs ot be applied to the main CL table created by INDXR:
9. CLCAL
task 'clcal'; default 'clcal';
timerang 0; gainuse 0; gainver 0;
snver 1; inver 1
10. PRTUV: task 'prtuv'; docrt 132; getn 1; go
-Scan length 4.0 seconds
11. Initial TVFLG.
-- ASSURE THAT YOU SET 'SWITCH ALL-CH FLAG'!!
**** Primary calibrator, RR:
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1;
dohist -1; sources ''; calcode 'K';
timerang 0; stokes 'RR'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- No additional flags
**** Primary calibrator, LL:
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1;
dohist -1; sources ''; calcode 'K';
timerang 0; stokes 'LL'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- No additional flags
**** Phase calibrator, RR:
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1;
dohist -1; sources ''; calcode 'D';
timerang 0; stokes 'RR'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- 7 flags entered into data.
**** Phase calibrator, LL:
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1;
dohist -1; sources ''; calcode 'D';
timerang 0; stokes 'LL'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- 7 flags entered into data.
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1; calcode '-cal';
dohist -1; sources '';
timerang 0; stokes 'RR'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 139; echan 159; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- 22 flags entered into data.
tvreset; tvclear;
task 'tvflg'; default tvflg;
getn 1; docat -1; calcode '-cal';
dohist -1; sources '';
timerang 0; stokes 'LL'; selband -1;
selfreq -1; freqid 1; bif 0; eif 0;
bchan 118; echan 138; antennas 0;
baseline 0; uvrange 0; subarray 0;
docalib 2; gainuse 0; blver -1;
flagver 0; doband -1; bpver -1;
smooth 0; dparm(6)=4; baddisk 0;
go
-- 22 flags entered into data.
12. SPFLG - edit line UV data using the TV display and cursor
======>> Primary CAL, RR:
task 'spflg'; default spflg; tvreset; tvclear;
getn 1; sources ''; calcode 'K'; timerang 0;
stokes 'RR'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib 2; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
--> 31 flags entered into data.
======>> Primary CAL, LL:
task 'spflg'; default spflg; tvreset; tvclear;
getn 1; sources ''; calcode 'K'; timerang 0;
stokes 'LL'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib 2; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
--> 10 flags entered into data.
======>> Phase CAL, RR:
task 'spflg'; default spflg; tvreset; tvclear;
getn 1; sources ''; calcode 'd'; timerang 0;
stokes 'RR'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib 2; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
--> 4 flags entered into dataset.
======>> Phase CAL, LL:
task 'spflg'; default spflg; tvreset; tvclear;
getn 1; sources ''; calcode 'D'; timerang 0;
stokes 'LL'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib 2; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
--> 10 flags entered into data.
=====>> Source, RR ----
task 'spflg'; default spflg; tvreset; tvclear;
getn 1; sources ''; calcode '-cal'; timerang 0;
stokes 'RR'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib 2; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
--> 31 flags entered into data.
=====>> Source, LL ---
task 'spflg'; default spflg; tvreset; tvclear;
getn 1; sources ''; calcode '-cal'; timerang 0;
stokes 'LL'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib 2; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
--> 8 flags entered into data.
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
--> FLAGGING IS COMPLETED.
--> Highest FG table version # = 11
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX
13. POSSM; examine bandpasses & identify reliable sections
-- Using antennas 11, 19, 22 (center of array) for BPASS examination
tvreset; tvclear;
task 'possm'; default 'possm'
getn 1; source '';
calcode 'K'; stokes ''; bchan 0; echan 0;
dotv 1; nplots 9; anten 13,24,21;
baseline 0;
docal 1; aparm 1,0; go
14. BPASS
task 'bpass'; default 'bpass';
getn 1; docal 1; solint 0; calcode 'K';
calsour '1331+305=3C286'; ichansel 40,216,1,0;
bpassprm(5)=1; bpassprm(10)=3; go
15. POSSM - examine bandpasses again, after BPASS has executed.
tvreset; tvclear;
task 'possm'; default 'possm'
getn 1; source '';
calcode 'K'; stokes ''; bchan 0; echan 0;
dotv 1; nplots 9;
anten 13,24,21; baseline 0;
docal 1; aparm 1,0; doband 3; go
16. SETJY: MAKE ABSOLUTELY SURE THAT calcode = '' (is blank)
task 'setjy'; default 'setjy';
getn 1; sources '1331+305=3C286';
qual -1; bif 0; eif 0; zerosp 0;
optype 'CALC'; calcode ''; sysvel 0;
restfreq 0; veltyp ''; veldef '';
freqid=1; aparm 0; go
--> SETJY1: '1331+305=3C286 ' IF = 1 FLUX =14.7662 (Jy calcd)
-- CHECK LISTR: task 'LISTR'; default 'listr'; getn 1; optype 'scan'; docrt 132; go
17. CALRD - read in a calibrator FITS image file
task 'calrd';
object '3c286'; band 'L'; go
18. CALIB: first, information on calibrators:
SITUATIONS WHERE SPECIAL RESTRICTIONS ARE NECESSARY FOR FLUX CALIBRATION
The following rules must be carefully followed to ensure proper
flux bootstrapping in the combinations of array scale and band noted
below. For the hybrid configurations (BnA, CnB, DnC) the rule for the
more compact configuration should be adopted (i.e. follow B config
rules for BnA). When specifying inner antennas to be used for the
calibration solution, no antenna on the North arm further out than on
the East or West arms should be used. Finally, it is a good idea to
set WTUV = 0.1 in CALIB to ensure a stable solution.
NUMBER OF
SOURCE BAND UVRANGE CONFIG INNER ANTENNAS NOTES
(cm) (Kl) (per arm)
-----------------------------------------------------------------------------
3C48/3C147 20 0-40 A 7
" B,C,D All
3C286 20 0-18 A 4
" B,C,D All
90-180 A All Reduce flux by 6%
---> UVRANGE for 3c286 is 0-18 for B,C,D arrays
--->> Phase cal = J1021+2159. "S" closure errors in B/C/D, 1.7 Jy at L-band,
UVRANGE = 0-30
task 'calib'; default 'calib'; getn 1;
get2n 2; calcode 'K'; qual -1;
timerang 0; uvrange 0 18;
docal 1; doband 3; solin 0; nmaps 1;
calsour '1331+305=3C286';
ichansel 40,216,1,0;
snver 2; inp;
---> Results:
---->>> Task CALIB (release of 31DEC11) begins
---->>> Using AGC731457 D . UVDATA . 1 Disk= 1 Usid= 9309
---->>> UVGET: Using flag table version 11 to edit data
---->>> Selecting, editing and calibrating the data
---->>> Using Clean Component self-cal source model
---->>> SETGDS: imaging done with reprojected tangent point(s)
---->>> FACSET: source model will be scaled to 14.766 Jy
---->>> FACSET: 16.044003 Jy found from 1424 components
---->>> FACSET: scaling factor set to 9.20354E-01
---->>> QINIT: did a GET of 5120 Kwords, OFF 4504213
---->>> VISDFT: Begin DFT component division
---->>> VISDFT: fields 1 - 1 chns 1 - 177 in 1 CC models
---->>> VISDFT: Model components of type Point
---->>> I Polarization model processed
---->>> VISDFT at visibility number 1
---->>> Model computation is 10 percent complete
---->>> Model computation is 20 percent complete
---->>> Model computation is 30 percent complete
---->>> Model computation is 40 percent complete
---->>> Model computation is 50 percent complete
---->>> Model computation is 60 percent complete
---->>> Model computation is 70 percent complete
---->>> Model computation is 80 percent complete
---->>> Model computation is 90 percent complete
---->>> Model computation is 100 percent complete
---->>> Field 1 used 1424 CCs
---->>> Determining solutions using full vector averaging
---->>> Writing SN table 2
---->>> RPOL, IF= 1 The average gain over these antennas is 8.512E+00
---->>> LPOL, IF= 1 The average gain over these antennas is 8.636E+00
---->>> Found 40 good solutions
---->>> Average closure rms = 0.00028 +- 0.00007
---->>> Fraction of times having data > 2.5 rms from solution
---->>> 1.00000 of the times had 2 - 4 percent outside 2.5 times rms
---->>> Adjusting solutions to a common reference antenna
---->>> Appears to have ended successfully
---->>> localhos 31DEC11 TST: Cpu= 80.4 Real= 82 IO= 503
Header:
AIPS 1: Maximum version number of extension files of type HI is 1
AIPS 1: Maximum version number of extension files of type FQ is 1
AIPS 1: Maximum version number of extension files of type AN is 1
AIPS 1: Maximum version number of extension files of type SU is 1
AIPS 1: Maximum version number of extension files of type WX is 1
AIPS 1: Maximum version number of extension files of type CL is 2
AIPS 1: Maximum version number of extension files of type NX is 1
AIPS 1: Maximum version number of extension files of type FG is 11
AIPS 1: Maximum version number of extension files of type SN is 2
AIPS 1: Maximum version number of extension files of type BP is 1
19. CALIB on phase cals:
clr2name; nmaps 0;
uvrange 0 30;
calsour 'J1021+2159'; calcode 'D';
go
--> Results:
---->>> Task CALIB (release of 31DEC11) begins
---->>> Using AGC731457 D . UVDATA . 1 Disk= 1 Usid= 9309
---->>> UVGET: Using flag table version 11 to edit data
---->>> Selecting, editing and calibrating the data
---->>> Doing cal transfer mode with point model for each source
---->>> This is not self-calibration
---->>> Dividing data by source flux densities
---->>> Determining solutions using full vector averaging
---->>> Writing SN table 2
---->>> RPOL, IF= 1 The average gain over these antennas is 8.466E+00
---->>> LPOL, IF= 1 The average gain over these antennas is 8.769E+00
---->>> Found 188 good solutions
---->>> Average closure rms = 0.00124 +- 0.00006
---->>> Fraction of times having data > 2.5 rms from solution
---->>> 0.25000 of the times had 0 - 2 percent outside 2.5 times rms
---->>> Adjusting solutions to a common reference antenna
---->>> Appears to have ended successfully
---->>> localhos 31DEC11 TST: Cpu= 1.0 Real= 5 IO= 744
20. GETJY. MAKE SURE THAT YOU USE ONLY THE HIGHEST SN TABLE VERSION, AND SET IT EXPLICITLY.
default 'getjy'; getn 1; source 'J1021+2159';
calcode '*'; snver 2;
calsour '1331+305=3C286'; bif 0; eif 0;
inp
--> J1021+2159 : 0 D 1 1.69027 +/- 0.01021 0.02859
--> Calibrator manual: 1.70 Jy; agreement within 0.57%
21. CLCAL
task 'clcal';
default 'clcal'; getn 1; calcode '*';
snver 2; invers snver; go
--> CL table, version #3, written by CLCAL
22. LISTR: print outputs for phase cals. Amplitudes should be roughly constant,
while phases should be close to zero.
Task 'LISTR'; getn 1; Optype 'matx';
Inext 'SN'; Inver 0; Sources '';
Calcode 'D'; Timerang 0; Stokes '';
Selband -1; Selfreq -1; Freqid 1;
Bif 0; Eif 0; Bchan 1; Echan 0;
Antennas 0; Baseline 0; Uvrange 0;
Subarray 0; Docalib 2; Gainuse 0;
Dopol -1; Blver -1; flagver 0;
Doband 3; Bpver 1; Smooth 0;
Dparm 5,1,0; doacor 0; Factor 0;
Docrt -1; baddisk 0;
Outprint 'VLADATA:AGC731457.D.LISTR.CLCAL';
--> Solution: Most phases within +/-5 degrees; a few outliers at
10-12 degrees.
--> If you need to reset calibration:
task 'extdest'; getn 1; inext 'sn'; invers 2; extdest
task 'extdest'; getn 1; inext 'cl'; invers 3; extdest
task 'extdest'; getn 1; inext 'bp'; invers 0; extdest
23. TVFLG ON SOURCE.
tvreset; tvclear;
default 'tvflg'; getn 1; docal 1;
stokes 'RR'; sources 'AGC731457';
doband 3; bchan 40; echan 216;
bif 0; eif 0; nchav echan-bchan+1;
calcode '-cal'; dparm(6)=4; go
tvreset; tvclear;
default 'tvflg'; getn 1; docal 1;
stokes 'LL'; sources 'AGC731457';
doband 3; bchan 50; echan 216;
bif 0; eif 0; nchav echan-bchan+1;
calcode '-cal'; dparm(6)=4; go
24. REWAY - replaces SPLIT, applies calibration and attempts to get the
data weighting correct.
task 'reway'; default 'reway'; getn 1;
docal 1; doband 3;
source 'AGC749741';
aparm 11,44,12,0,10,4;
calcode '-cal'; bpver 0;
outname 'agc731457';
outclass 'split';
-- Note on aparm: to use a rolling buffer of 11 times separated by no
more than 44 seconds and then smoothed further with a Gaussian 12
seconds in FWHM. Data are flagged if the rms is more than 4 times
the variance away from the mean averaged over all baselines, IFs,
and polarizations. Flagging on the variance of the rms from the
mean on a baseline basis is essentially turned off by the 10.
-- This is the suggestion from APPENDIX E of the AIPS COOKBOOK; it has changed.
25. SPFLG on final split data:
=====>> Source, RR
task 'spflg'; default spflg; tvreset; tvclear;
getn 3; sources '''; calcode ''; timerang 0;
stokes 'RR'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib -1; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
=====>> Source, LL
task 'spflg'; default spflg; tvreset; tvclear;
getn 3; sources '''; calcode ''; timerang 0;
stokes 'LL'; selband -1; selfreq -1;
freqid -1; bif 0; eif 0; bchan 1; echan 0;
antenn 0; baseli 0; docat -1;
uvrange 0; subarray 0; docalib -1; gainuse 0;
dopol -1; blver -1; flagver 0; outfgver 0;
doband -1; bpver -1; smooth 0;
dparm 0; antennas 0;
baddisk 0; dparm(6)=4; go
Ticket Details
===================
Ticket ID: NTC-279643
Department: AIPS Data Processing
Priority: Default
Status: Open
Link: https://help.nrao.edu/staff/index.php?_m=tickets&_a=viewticket&ticketid=1831
More information about the Daip
mailing list