Calibration Checks ------------------------------------------------ CHECK BANDPASS SOLUTIONS ------------------------------------------------ Check Bandpass: 1) POSSM *** *** Check bandpass solutions by plotting the BP table. *** *** Smooth functions are okay. Juergen would be worried if the *** phases wrap around a lot. (quick wrapping) *** DEFAULT POSSM INDISK 5; QUAL -1; CALCODE ''; SELBAND -1; SELFREQ -1; FREQID 1; UVRANG 0, 0; TIMERANG 0; STOKES ''; BIF 0; EIF 0; BCHAN 1;ECHAN 0; SUBARRAY 0; ANTENNAS 0; BASELINE 0; DOCALIB -1; GAINUSE 2; DOPOL -1; BLVER -1; DOBAND -1; BPVER 1; SMOOTH 0; SHIFT 0, 0; CODETYPE ''; POLPLOT ''; BPARM 0; OUTFILE ''; LTYPE 3; FACTOR 0; XYRATIO 0; BADDISK 0; GRCHAN 0; FLAGVER 1 APARM 0, 1, 0.7, 1.3, -180, 180, 0, 2, 0, 0 $ Plot BP, with amp/ph ranges SOURCE flux1, flux2, '' $ POSSM doesn't work with $ source '' for some reason! SOLINT -1 $ Separate plots for each scan NPLOTS 9 $ 9 plots per page BPARM 0 DOTV 1 FREQID 1; BPVER 0; getn (+2) $ *.LINCOP 2) POSSM *** *** Apply BP table to secondary calibrator and vector-average all *** data. *** *** Your secondary calibrator should be quite flat in both *** amplitude and phase as a function of frequency, with no slope. *** If some channels are off, note which ones those are and keep an *** eye out for interference or other bad data. *** If there are large errors consider running BPASS on the *** secondary calibrator and using that to correct the galaxy. Note *** that this will be somewhat painful since AIPS does not allow *** incremental BP tables -- unlike SN or CL tables. Sigh. *** *** *** We can try anten=vla;baseline=vla;aparm 0 for scalar *** averaging of VLA antennas. This should be flat. *** *** NOTE: To plot versus frequency, use aparm(7) = 1 DEFAULT POSSM INDISK 5; QUAL -1; CALCODE ''; SELBAND -1; SELFREQ -1; FREQID 1; UVRANG 0, 0; TIMERANG 0; STOKES ''; BIF 0; EIF 0; BCHAN 1;ECHAN 0; SUBARRAY 0; ANTENNAS 0; BASELINE 0; DOCALIB -1; GAINUSE 2; DOPOL -1; BLVER -1; DOBAND -1; BPVER 1; SMOOTH 0; SHIFT 0, 0; CODETYPE ''; POLPLOT ''; BPARM 0; OUTFILE ''; LTYPE 3; FACTOR 0; XYRATIO 0; BADDISK 0; GRCHAN 0; FLAGVER 0 APARM 0 $ Plot data SOLINT 0 $ average all time NPLOTS 0 $ average all baselines APARM 0 APARM(1) 1 $ vector average $ APARM(1) 0 $ ams: scalar average anten=vla baseline=vla $ ams: we can only trust vla DOCAL 1 ; GAINUSE 3 ; DOBAND 1 $ average all BP entries DOTV 1 FREQID 1; BPVER 0 $ for single-FREQID data sets source=phase,'' $ Secondary (phase) calibrator getn (+2) $ *.LINCOP ------------------------------------------------ CHECK CALIB SOLUTIONS ------------------------------------------------ 3) SNPLT - phase DEFAULT SNPLT; INDISK 5; INEXT 'CL'; INVERS 0; SOURCES ''; QUAL -1; TIMERANG 0; STOKES ''; SELBAND -1; SELFREQ -1; FREQID 1; SUBARRAY 0; BIF 0; EIF 0; ANTENNAS 0; PIXRANGE 0; NPLOTS 4; XINC 1; OPTYPE 'PHAS'; OPCODE 'ALSI'; DO3COL 1; BCOUNT 1; XAXIS 0; SYMBOL 5; FACTOR 2; CUTOFF 0; LTYPE 3; DOTV 1; GRCHAN 0; pixrange -180, 180; getn (+3); $ *.LINCOP 4) SNPLT - amplitude DEFAULT SNPLT; INDISK 5; INEXT 'CL'; INVERS 0; SOURCES ''; QUAL -1; TIMERANG 0; STOKES ''; SELBAND -1; SELFREQ -1; FREQID 1; SUBARRAY 0; BIF 0; EIF 0; ANTENNAS 0; PIXRANGE 0; NPLOTS 4; XINC 1; OPTYPE 'AMP'; OPCODE 'ALSI'; DO3COL 1; BCOUNT 1; XAXIS 0; SYMBOL 5; FACTOR 2; CUTOFF 0; LTYPE 3; DOTV 1; GRCHAN 0; pixrange 2 6; getn (+3); $ *.LINCOP ------------------------------------------------ CHECK CALIBRATOR AMP AND PHASE ------------------------------------------------ 5a) UVPLT *** *** Check the amplitude/phase vs. uv-distance for all calibrators. *** Amplitude should match the results of SETJY/GETJY. If there are *** obvious outliers which are not expected due to source *** structure, go back and flag those (and possibly re-run CALIB, *** etc.) *** *** *** If you see problems in a plot with amp/phase versus uv distance, *** try plotting versus time. *** to do this, set bparm(1) = 11 *** *** for galactic hi, you have to tacop the cl table to the *** newch0 to accurately check *** DEFAULT UVPLT DOCAL 1 ; GAINUSE 5 $ Apply calibration FLAGVER 0 $ set this to the latest FG version DOTV 1 DO3COL 1 BPARM 0 $ amp vs. uv-distance CALCODE '*' getn (+3) $ *.NEWCH0.1 $ do this for each NEWCH0 file 5b) UVPLT *** *** Check phase vs. uv-distance *** DEFAULT UVPLT DOCAL 1 ; GAINUSE 0 $ Apply calibration FLAGVER 0 $ set this to the latest FG version DOTV 1 DO3COL 1 BPARM 0,2 $ phase vs. uv-distance CALCODE '*' source phase,'' getn (+3) $ *.NEWCH0.1 $ do this for each NEWCH0 file 6a) UVPLT *** *** Real versus imaginary *** DEFAULT UVPLT DOCAL 1 ; GAINUSE 0 $ Apply calibration FLAGVER 0 $ set this to the latest FG version DOTV 1 DO3COL 1 BPARM 9,10 $ real vs. imaginary $ no arcs! looks like a donut. CALCODE '*' source flux1,'' getn (+3) $ *.NEWCH0.1 $ do this for each NEWCH0 file 6b) UVPLT DEFAULT UVPLT DOCAL 1 ; GAINUSE 0 $ Apply calibration FLAGVER 0 $ set this to the latest FG version DOTV 1 DO3COL 1 BPARM 0,2 $ phase vs. uv-distance CALCODE '*' source phase,'' getn (+3) $ *.NEWCH0.1 $ do this for each NEWCH0 file 6a) TVFLG DEFAULT TVFLG INDISK 5; DOCAT -1; IN2SEQ 0; IN2DISK 0; DOHIST -1; SOURCES ''; CALCODE ''; TIMERANG 0; STOKES ''; SELBAND -1; SELFREQ -1; FREQID 1; BIF 0; EIF 0; BCHAN 1; ECHAN 0; NCHAV 1; CHINC 1; ANTENNAS 0; BASELINE 0; UVRANGE 0, 0; SUBARRAY 0; DOPOL -1; BLVER -1; FLAGVER 1; OUTFGVER 1; DOBAND -1; BPVER -1; SMOOTH 0; DOCALIB 1; GAINUSE 0 $ Apply calibrations DPARM(3) = 1 $ Show baselines twice, treat all antennas the $ same. this shows baseline 27-1 as well as $ 1-27 DPARM(6) = 10 $ set to calibrator's integ. time in seconds $ (use prtuv to get this) source flux1, flux2,'' getn (+3) $ *.NEWCH0 6b) TVFLG (also check phase cal) source phase,'' go tvflg also time - just phase cal ------------------------------------------------ CHECK CALIBRATOR AMP AND PHASE ------------------------------------------------ 7a) ANBPL *** *** Checks data weights. Data with very high weights (factor fo *** 5-10 more than normal) should be flagged with UVFLG. *** *** Occasionally, Adrienne has seen slopes in the weights and isn't *** sure what to do about them. She thinks they may be caused by *** slopes in antenna gains over time. *** DEFAULT ANBPL INDISK 5; QUAL -1; SELBAND -1; SELFREQ -1; FREQID -1; BCHAN 1; ECHAN 0; BIF 0; EIF 0; AVGIF 0; XINC 1; UVRANGE 0 0; TIMERANG 0; STOKES ''; ANTENNAS 0; BASELINE 0; DOPOL -1; BPVER -1; DOBAND -1; BPVER -1; SMOOTH 0; SUBARRAY 0; SOLINT 0; SYMBOL 0; LTYPE 3; DOTV 1; GRCHAN 0; OUTPRINT ''; NDIG 0; DOCALIB 1;GAINUSE 0 FLAGVER 0 BPARM 2,17,0 $ Plot antenna-based weight vs. time NPLOTS 9; DOTV 1 DOCRT 132 $ Print as well as plotting weights -- useful for $ finding exact times of bad weights $ Note you can also use outprint to send to a file. OPCODE 'alsi' $ Plot all IFs together DO3COL 1 $ ...using different colors getn (+3) $ *.NEWCH0