# # Template Checklist and Worksheet for reduction of VSOP Survey data. # (http://halca.vsop.isas.ac.jp/survey/template_checklist.txt) # # Version 0.3 - J. Lovell 14-Apr-1999 # # NOTE: Copy this file to your data reduction outputs directory (e.g # /something/vs05a_jones) and name it "checklist.txt" # # Comments are given on lines beginning with a # and questions are # enclosed in {}'s. # Please enter responses after the colons. # Multiple-line responses are OK. # Yes/No questions may be answered with "yes", "no", "Y", "N", "y" or "n" # An example of a completed Checklist is given in # http://halca.vsop.isas.ac.jp/survey/checklist_eg.txt # # # Do not delete or alter anything past this line #--------------------------------------------------------------------------- # Section 1 : General Information #--------------------------------------------------------------------------- {Experiment code} : VS05q {Your family name} : Lovell {Your given name} : Jim {The date you prepared this document (format is day-month-year, e.g 23-04-1999)} : 27-04-1999 {The name of the survey source that was observed} : J1316-33 {The name of the calibration source(s) observed} : 1334-127 {The year the experiment took place (e.g. 1998)} : 1998 {The starting day number of the experiment} : 20 {List the GRTs that observed} : Hart, Shanghai, Mopra, Kashima, Usuda {List the tracking stations that took part} : NZ, TZ {Which correlator was used to correlate the data? (Mitaka/Penticton/Socorro)} : Penticton {AIPS Version used (e.g 15OCT98)} : 15OCT98 {Difmap Version used (e.g 2.3c)} : 2.3c {Was vsop_difmap used? (yes/no)} : yes {Abstract. Write a brief summary of the data reduction and its outcomes} : Fringes to HALCA were found during only part of the NZ pass (about one hour in total). No fringes in IF-2 at Shanghai and very low amplitudes at Kashima, suggesting a serious problem. Kashima data were eventually flagged. Analysis of the calibrator source showed Shanghai gain to be underestimated. Multiplying Shanghai gains by a factor of ~1.7 resolved this problem. Usuda gains were also found to be low by a factor of ~1.6. An image and modelfit of the survey source was made. EBF: 07-OCT-2005: Some editing of data core is 0.3 +/- 0.1 #--------------------------------------------------------------------------- # Section 2 : Pre-AIPS Checklist #--------------------------------------------------------------------------- {PI letter retrieved (yes/no)} : yes {Data quality reports and plots retrieved? (yes/no)} : yes {Pre-launch 5GHz survey data retrieved? (yes/no)} : yes {Light-curve data retrieved? (yes/no)} : yes # Input files for INDXR, ANTAB and UVFLG {Has a INDXR input file (INDXR.TXT) been created? (yes/no)} : n {Has a ANTAB input file (ANTAB.TXT) been created? (yes/no)} : y {Has a UVFLG input file (UVFLG.TXT) been created? (yes/no)} : n {Any comments on Pre-AIPS file preparation?} : No INDXR.TXT needed as there were no gain jumps during this experiment. #--------------------------------------------------------------------------- # Section 3 : AIPS Output File Checklist # # Names of expected output file names are contained in square brackets. #--------------------------------------------------------------------------- {PRTAN [PRTAN.LIS] (yes/no)} : y {LISTR (optyp = 'SCAN') [LISTR.LIS] (yes/no)} : y {DTSUM (aparm = 0) [DTSUM_0.LIS] (yes/no)} : y {DTSUM (aparm = 1) [DTSUM_1.LIS] (yes/no)} : y {UVPLT (u vs v) [UVPLT_UV.PS] (yes/no)} : y {UVPLT ((u,v) distance vs time) [UVPLT_UVD.PS] (yes/no)} : y {UVPLT ((u,v) position angle vs time). [UVPLT_UVPA.PS] (yes/no)} : y {VPLOT antenna weights vs time (after adjustment) [WEIGHTS.PS] (yes/no)} : n {ACCOR SNPLT (amplitude), IF 1 [SN_ACCOR_1.PS] (yes/no)} : n {IF 2 [SN_ACCOR_2.PS] (yes/no)} : n {APCAL SNPLT (amplitude), IF 1 [SN_APCAL_1.PS] (yes/no)} : y {IF 2 [SN_APCAL_2.PS] (yes/no)} : y # FRING SNPLT (delay, rate, SNR), as many as necessary to document the # clear discovery of fringes at all antennas, numbered with "nn" for # both IFs. Do NOT include the final SN table as it is plotted later. # # Format is SN_FRINGnn_DEL_i.PS, SN_FRINGnn_RAT_i.PS, # SN_FRINGnn_SNR_i.PS for delay, rate and SNR plots respectively, # where "i" is the IF number (1 or 2). {Number of times preliminary FRING solutions were plotted (i.e maximum value of nn)} : 2 {Post-FRING CALIB SNPLT, Phase: IF 1 [SN_CALIB_PHA_1] (yes/no)} : y {IF 2 [SN_CALIB_PHA_2] (yes/no)} : y {Delay: IF 1 [SN_CALIB_DEL_1] (yes/no)} : y {IF 2 [SN_CALIB_DEL_2] (yes/no)} : y {Rate: IF 1 [SN_CALIB_RAT_1] (yes/no)} : y {IF 2 [SN_CALIB_RAT_2] (yes/no)} : y {SNR: IF 1 [SN_CALIB_SNR_1] (yes/no)} : y {IF 2 [SN_CALIB_SNR_2] (yes/no)} : y {Post-calibration VPLOT (amp and phase, with channel averaging) IF 1 [VPLOT_AP_1.PS] (yes/no)} : n {IF 2 [VPLOT_AP_1.PS] (yes/no)} : n {Post-calibration POSSM (amplitude and phase, at least one 5 to 10 min integration spectrum per tracking pass, indexed with nn = 01, 02, 03, ...), IF 1 [POSSMnn_AP_1.PS] (yes/no)} : n {IF 2 [POSSMnn_AP_2.PS] (yes/no)} : n {Output TASAV FITS file (containing all 'multisource' cal tables) [VSxxx__TASAV.FITS] (yes/no)} : y {Calibrated and spectrally averaged FITS data files (one for each source), [__AIPS.FITS] (yes/no)} : y #--------------------------------------------------------------------------- # Section 4 : Supplementary files # # Were any files created, other than those listed, to document the data # reduction? If so, please list them below using as many (File, Description) # groups as required. #--------------------------------------------------------------------------- #--------------------------------------------------------- {File} : POSSM01_AP.PS {Description} : Post-calibration spectra (IF 1 and 2) at 0/21:00 #--------------------------------------------------------- {File} : POSSM02_AP.PS {Description} : Post-calibration spectra (IF 1 and 2) at 1/00:20 #--------------------------------------------------------- #--------------------------------------------------------------------------- # Section 5 : AIPS data-reduction notes # # The format is as follows: # # Task : the name of the task that is run or other operation performed # ADVERBs : a list of the relevant non-default adverbs # Input CL : the input calibration table used (if any) in running the task # Input SN : the SN table that was used (if any) # Output SN : the SN table produced (if any) # Output CL : the CL table produced or modified (if any) # Comments: any comments on the results of the operation # # Use as many of the (Task, ADVERBs, Comments, Input CL, Input SN, # Output SN, Output CL) groups as required. #--------------------------------------------------------------------------- #--------------------------------------------------------- {Task} : FITLD {ADVERBs} : douvcomp = 1 {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Loaded previously calibrated data for re-analysis #--------------------------------------------------------- {Task} : EXTDEST {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Removed previous CL, SN, TY, GC, NX tables to return data to an uncalibrated state #--------------------------------------------------------- {Task} : MSORT {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Not really necessary, but did it anyway. #--------------------------------------------------------- {Task} : INDXR {ADVERBs} : cparm(3) 0.25 {Input CL} : {Input SN} : {Output SN} : {Output CL} : 1 {Comments} : #--------------------------------------------------------- {Task} : UVFLG {ADVERBs} : antennas = 2,0 bif 2; eif 2 opcode 'flag' reason 'Shanghai IF 2 bad' {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Flagged shaghai IF 2 as it is known to contain invalid data. #--------------------------------------------------------- {Task} : POSSM {ADVERBs} : docalib -1 {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : These plots and DQA show Mopra to be best reference antenna #--------------------------------------------------------- {Task} : ANTAB {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Gain curve information from VSOP calibration web page. Nominal Tsys used for Kashima, the rest from log files. Limited data available from Mopra and Usuda. #--------------------------------------------------------- {Task} : APCAL {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : 1 {Output CL} : {Comments} : #--------------------------------------------------------- {Task} : CLCAL {ADVERBs} : refant 4 {Input CL} : 1 {Input SN} : 1 {Output SN} : {Output CL} : 2 {Comments} : Amplitude calibration solutions applied #--------------------------------------------------------- {Task} : TACOP {ADVERBs} : inext 'cl' {Input CL} : 2 {Input SN} : {Output SN} : {Output CL} : 3 {Comments} : copy CL2 to CL3 to protect CL2 in case of a mistake in CLCOR #--------------------------------------------------------- {Task} : CLCOR {ADVERBs} : opcode 'gain' clcorprm=0.906,0 gainver 3 {Input CL} : 3 {Input SN} : {Output SN} : {Output CL} : 3 {Comments} : Penticton gain correction #--------------------------------------------------------- {Task} : CLCOR {ADVERBs} : opcode 'sbdl' antennas = 1,0 clcorprm -480,-480 antennas = 2,0 clcorprm 1480,1480 antennas = 3,0 clcorprm 1845,1845 antennas = 5,0 clcorprm -1700,-1700 antennas = 6,0 clcorprm -1888,-1888 {Input CL} : 3 {Input SN} : {Output SN} : {Output CL} : 3 {Comments} : Known delay offsets from DQA report applied before fringe-fitting. #--------------------------------------------------------- {Task} : KRING {ADVERBs} : docalib = 1 gainuse = 3 dofit = 0 refant = 0 soltyp 'nols' search = 6, 4, 1 solmode = 'SBD' bparm(5) = 0.5 cparm = 0, 0, 2, 1000, 100, -1, -1, -1, 0, 0 antwt = 0 solint = 6 prtlev 1 {Input CL} : 3 {Input SN} : {Output SN} : 2 {Output CL} : {Comments} : fringes found at most times on all baselines except to HALCA on the second tracking pass, despite the projected baselines crossing some ground baselines. #--------------------------------------------------------- {Task} : KRING {ADVERBs} : docalib = 1 gainuse = 3 dofit = 0 refant = 0 soltyp 'nols' search = 6, 4, 1 solmode = 'SBD' bparm(5) = 0.5 cparm = 0, 0, 2, 2500, 1000, -1, -1, -1, 0, 0 antwt = 0 solint = 12 prtlev 1 timerang 0,23,30,00,1,1,30,00 {Input CL} : 3 {Input SN} : {Output SN} : 3, 4 {Output CL} : {Comments} : concentrating on second tracking pass. Increased solint to 12min but still no convincing fringes. Return to FRING and 6 min solint with smaller windows for final fringe-fit. #--------------------------------------------------------- {Task} : FRING {ADVERBs} : docalib = 1 gainuse = 3 refant = 4 search = 6, 1 dparm = 3, 200, 100, 0.5 solint = 6 aparm = 0, 0, 0, 0, 0, 1, 0.1, 0, 1, 0 {Input CL} : 3 {Input SN} : {Output SN} : 5 {Output CL} : {Comments} : Final fringe-fit. #--------------------------------------------------------- {Task} : TBOUT {ADVERBs} : outfile 'survey:sn5.txt' docrt = 1000 {Input CL} : {Input SN} : 5 {Output SN} : {Output CL} : {Comments} : Wrote fringe fit SN table out for editing in snedit.pl #--------------------------------------------------------- {Task} : snedit.pl {ADVERBs} : output file SN6.TXT {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Used to flag bad solutions. Much easier than SNSMO #--------------------------------------------------------- {Task} : TBIN {ADVERBs} : infile 'survey:sn6.txt' {Input CL} : {Input SN} : {Output SN} : 6 {Output CL} : {Comments} : Edited SN table read in #--------------------------------------------------------- {Task} : CLCAL {ADVERBs} : refant 4 {Input CL} : 3 {Input SN} : 6 {Output SN} : {Output CL} : 4 {Comments} : Amplitude calibration solutions applied #--------------------------------------------------------- {Task} : POSSM {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Plot 10min at two times: 0/21:00 to 21:10 and 1/00:20 to 1/00:30. Shows big phase variation near edge of bandpass for Shanghai and smaller variations for Usuda. Chans 30 to 220 OK #--------------------------------------------------------- {Task} : SPLIT {ADVERBs} : aparm = 2, 0.5, 0 chansel 30, 220, 1 flagver 1 docalib 1 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : {Comments} : #--------------------------------------------------------- {Task} : FITTP {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Data written out for analysis in Difmap #--------------------------------------------------------- {Task} : CLCOR {ADVERBs} : opcode 'gain' antennas 2,0 clcorprm 1.8,0 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : 4 {Comments} : Gain correction for Shanghai after Difmap examination. #--------------------------------------------------------- {Task} : SPLIT {ADVERBs} : aparm = 2, 0.5, 0 chansel 30, 220, 1 flagver 1 docalib 1 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : {Comments} : #--------------------------------------------------------- {Task} : FITTP {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Data written out for analysis in Difmap again #--------------------------------------------------------- {Task} : CLCOR {ADVERBs} : opcode 'gain' antennas 6,0 clcorprm 1.7,0 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : 4 {Comments} : Gain correction for Usuda after second Difmap examination. #--------------------------------------------------------- {Task} : SPLIT {ADVERBs} : aparm = 2, 0.5, 0 chansel 30, 220, 1 flagver 1 docalib 1 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : {Comments} : #--------------------------------------------------------- {Task} : FITTP {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Data written out for analysis in Difmap again! #--------------------------------------------------------- {Task} : CLCOR {ADVERBs} : opcode 'gain' antennas 5,0 clcorprm 32.5,0 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : 4 {Comments} : Gain correction for Kashima after third Difmap examination. #--------------------------------------------------------- {Task} : SPLIT {ADVERBs} : aparm = 2, 0.5, 0 chansel 30, 220, 1 flagver 1 docalib 1 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : {Comments} : #--------------------------------------------------------- {Task} : FITTP {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Data written out for analysis in Difmap again! #--------------------------------------------------------- {Task} : CLCOR {ADVERBs} : opcode 'gain' antennas 2,6,0 clcorprm 0.92,0 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : 4 {Comments} : Gain correction for Usuda and Shanghai after fourth Difmap examination. #--------------------------------------------------------- {Task} : SPLIT {ADVERBs} : aparm = 2, 0.5, 0 chansel 30, 220, 1 flagver 1 docalib 1 {Input CL} : 4 {Input SN} : {Output SN} : {Output CL} : {Comments} : #--------------------------------------------------------- {Task} : FITTP {ADVERBs} : {Input CL} : {Input SN} : {Output SN} : {Output CL} : {Comments} : Data written out for analysis in Difmap again! #--------------------------------------------------------- #--------------------------------------------------------------------------- # Section 6 : Difmap data reduction notes #--------------------------------------------------------------------------- {Phase Selfcal applied? (yes/no)} : y {Global Amplitude Selfcal applied? (yes/no)} : y {Amplitude Selfcal applied on non-global timescales? (yes/no)} : n {Modelfit made and saved using the vs_final_mod macro? (yes/no)} : y {Image made and saved using the vs_final_img macro? (yes/no)} : y {Comments on Difmap processing. Include the global amplitude offsets that were applied if gscal was used} : After initial AIPS extraction, 1334-127 data loaded and examined. Only Hart, Mopra and Shanghai observed this source. By comparison with VLBA data (and accounting for variability), Shanghai gain may be down by a factor of about 1.8. First look at J1316-33 supports this and indicates Usuda may be low too. CLCOR in AIPS used to apply 1.8 Gain correction. Corrected 1334-127 data looks OK now. A simple single circular Gaussian was modelfitted and saved with vs_final_mod. The corrected J1316 data looks a little odd. Ududa gain may be out by about a factor of 1.7. CLCOR used to apply this and data exported again. J1316 data now looks OK. Flux on short Usuda-Shang baseline ~1.5 Jy maybe a bit high (expect about 1.2 Jy) but will see how it goes. Also, Kashima data has very low amplitudes! Assuming Mopra and Usuda gains are now OK, the Mop-Kashima (0/22:30) and Mop-Usuda (0/21:00) baselines are almost the same, so can guess Kashima correction = 1.3 Jy / 0.04 Jy = 32.5. Corrected in AIPS, but resulting Kashima data is very noisy - bad pointing perhaps. Decided to flag Kashima. After best possible image made with just cleam/selfcal, gscal was run. The result of this flagged all of IF2 (because no amplitude closure) and in IF1, Shanghai and Usuda were re-scaled by 0.91 and 0.93 respectively while VSOP and Mopra were increased by 1.04 and 1.07 respectively (no solution for Hart as no closure). This is most likely an over-compensation of Usuda and Shanghai gains, so applying just those corrections (0.92 each) in AIPS with CLCOR. Looked at 1334-127 data again, fitted simple circular gaussian to data and saved using vs_final_mod. J1316 imaged again for last time. Used selfflag to stop flagging of unusable data after phase-selfcal because Mopra-Hart baseline provides constraints that would otherwise be lost with no significant loss of image quality as a result. Saved with vs_final_img. Model fitting: Started with two circular Gaussians. vs_wplot suggests a VSOP weight increase of about 100. This is set using selfant vsop,false,100. Modelfit converged after about 70 iterations with uvstat(rms) = 0.173 Jy. Changed secondary to elliptical Gaussian and re-converged to get uvstat(rms) = 0.171 Jy, so no significant improvement and circulars Gaussians will do. saved with vs_final_mod. EBF: 08-OCT-2005: Some editing. See uv data in J1616_FOMALONT_AIPS.FITS Core 0.3 +/- 0.1 #--------------------------------------------------------------------------- # Do not delete or alter anything past this line {VSOP Survey Checklist and Worksheet version} : 0.3 #---------------------------------------------------------------------------