Telescope notes: - BR: - FD: - HN: - KP: - LA: - MK: - NL: - OV: - PT: - SC: * ** Get files from archive, and then if necessary, run *rename_archive.pl * *setenv U $PWD *export U=$PWD */project/ariel/AIPS-31DEC05/START_AIPS.old *137 * *task 'snplt'; default; tput snplt *task 'possm'; default; tput possm *task 'fring'; default; tput fring *task 'snsmo'; default; tput snsmo *task 'split'; default; tput split *task 'prtab'; default; tput prtab *task 'calib'; default; tput calib *task 'clcal'; default; tput clcal *task 'clcor'; default; tput clcor *task 'bpass'; default; tput bpass *task 'vlog'; default; tput vlog *task 'uvavg'; default; tput uvavg *task 'setjy'; default; tput setjy *task 'indxr'; default; tput indxr *task 'lwpla'; default; tput lwpla *task 'tacop'; default; tput tacop * prtask '';clrmsg task 'fitld' INTAPE = -1; NCOUNT = 1; OUTNAME = 'VLBACODE'; outcl 'UV2CM'; outdisk 1; DOUVCOMP = FALSE ;CLINT = 0.5; outcl 'uv2cm' infil 'U:f ncount 100 infil 'U:VLBACODE.FITS **************** **************** vlbaload; indisk 1; getn 1; vlbafix *delete other frequency files ********** * For 512 Mbps experiments, may need to vlbaload 2 files separately, * then combine them with VBGLU ***************** ***************** ******************* recat task 'edita'; default;indisk 1; getn 1; inext 'ty'; crowded 1;freqid 0;inver 1;flagver 1;go docrt false ;outprint 'U:VLBASUMM.TXT vlbasumm task 'snplt'; inext 'ty'; pixr 0; timer 0;anten 0 inv 1; nplot 8; optyp 'tsys'; source '';dotv false; stokes ''; runwait('snplt'); plver 1; inv 20; aspmm 0; outfil 'U:TSYS.PS runwait('lwpla'); inext 'pl'; inv -1;extdest REFANT = antnum('REFANTENNA') optyp ''; aparm 0; clrst; SOURCE(1) = 'src01';zerosp flux01,0; runwait('setjy') SOURCE(1) = 'src02';zerosp flux02,0; runwait('setjy') SOURCE(1) = 'src03';zerosp flux03,0; runwait('setjy') SOURCE(1) = 'src04';zerosp flux04,0; runwait('setjy') SOURCE(1) = 'src05';zerosp flux05,0; runwait('setjy') SOURCE(1) = 'src06';zerosp flux06,0; runwait('setjy') SOURCE(1) = 'src07';zerosp flux07,0; runwait('setjy') SOURCE(1) = 'src08';zerosp flux08,0; runwait('setjy') SOURCE(1) = 'src09';zerosp flux09,0; runwait('setjy') SOURCE(1) = 'src10';zerosp flux10,0; runwait('setjy') SOURCE(1) = 'src11';zerosp flux11,0; runwait('setjy') SOURCE(1) = 'src12';zerosp flux12,0; runwait('setjy') SOURCE(1) = 'src13';zerosp flux13,0; runwait('setjy') SOURCE(1) = 'src14';zerosp flux14,0; runwait('setjy') SOURCE(1) = 'src15';zerosp flux15,0; runwait('setjy') SOURCE(1) = 'src16';zerosp flux16,0; runwait('setjy') SOURCE(1) = 'src17';zerosp flux17,0; runwait('setjy') SOURCE(1) = 'src18';zerosp flux18,0; runwait('setjy') SOURCE(1) = 'src19';zerosp flux19,0; runwait('setjy') SOURCE(1) = 'src20';zerosp flux20,0; runwait('setjy') SOURCE(1) = 'src21';zerosp flux21,0; runwait('setjy') SOURCE(1) = 'src22';zerosp flux22,0; runwait('setjy') SOURCE(1) = 'src23';zerosp flux23,0; runwait('setjy') SOURCE(1) = 'src24';zerosp flux24,0; runwait('setjy') SOURCE(1) = 'src25';zerosp flux25,0; runwait('setjy') SOURCE(1) = 'src26';zerosp flux16,0; runwait('setjy') SOURCE(1) = 'src27';zerosp flux17,0; runwait('setjy') SOURCE(1) = 'src28';zerosp flux18,0; runwait('setjy') SOURCE(1) = 'src29';zerosp flux19,0; runwait('setjy') SOURCE(1) = 'src30';zerosp flux20,0; runwait('setjy') SOURCE(1) = 'src31';zerosp flux21,0; runwait('setjy') SOURCE(1) = 'src32';zerosp flux22,0; runwait('setjy') docrt -1;inext 'su';inv 1;outprint 'U:SUTABLE.TXT';box 2 5;runwait('prtab') timer 0; solint = 2.0; RUNWAIT ('ACCOR') SAMPTYPE = 'MWF'; DOBLANK = -1; SMOTYPE = 'AMPL' CPARM = 0.5, 0, 0, 0, 0, 0.02; INVERS = MAXTAB('SN') OUTVERS = INVERS + 1; RUNWAIT ('SNSMO') INEXT = 'SN'; EXTDEST; VBA_SN = INVERS INVERS = VBA_SN + 1; OUTVERS = VBA_SN; OUTNAME = INNAME INEXT = 'SN';keyw ''; OUTCLASS = INCLASS; OUTSEQ = INSEQ; OUTDISK = INDISK RUNWAIT ('TACOP'); EXTDEST SOUCODE = ' '; CALSOUR = ' '; QUAL = -1; CALCODE = ' ' OPCODE = 'CALI'; INTERPOL = '2PT'; INTPARM = 0; CUTOFF = 0 SMOTYPE = ''; SNVER = VBA_SN; DOBLANK = -1;source ''; GAINVER = MAXTAB('CL'); GAINUSE = GAINVER + 1; RUNWAIT ('CLCAL') task 'snplt'; inext 'sn'; dotv -1;nplot 8; source ''; inv 0; stokes = ''; optyp 'amp'; opcode ''; runwait('snplt'); plver = 1; inver = maxtab('pl'); outfil 'U:ACCOR.PS' ; runwait('lwpla'); inext 'pl'; inv -1;extdest infil ''; SOUR = ' ';timer 0; STOKES = ' '; TYVER = MAXTAB('TY'); GCVER = MAXTAB('GC') SNVER = MAXTAB('SN') + 1; OPCODE = 'GRID'; APARM = 0; SOLINT = 0 TAU0 = 0; DOFIT = 1; PRTLEV = 0; INVERS = 1 DOTV = FALSE; LTYPE = 3; APARM(5) = 4; TREC = 0; TREC=73,61, 47,53, 67,64, 64,64, 52,59, 44,52, 56,55, 64,66, 53,53, 54,62 dofit=1,1,1,1,1,1,1,1,1,1 RUNWAIT ('APCAL') task 'lwpla'; plver = 1;inver =maxtab('pl'); outfil 'U:OPACITY.PS' ; runwait('lwpla'); inext 'pl'; inv -1;extdest task 'snedt'; inext 'sn'; solint 0.5; inver=maxtab('sn');timer 0; anten 0; expert 0;dotwo -1;go ; SOUCODE = ' '; CALSOUR = ' '; QUAL = -1; CALCODE = ' ' OPCODE = 'CALI'; INTERPOL = 'SELF'; INTPARM = 0; CUTOFF = 0 SMOTYPE = ' '; SNVER = MAXTAB('SN'); GAINVER = MAXTAB('CL') GAINUSE = GAINVER + 1; RUNWAIT ('CLCAL') task 'snplt'; inext 'cl'; dotv -1;nplot 8; inv 0; stokes = ''; optyp 'amp'; opcode ''; runwait('snplt'); plver = 1; inver = maxtab('pl'); outfil 'U:APCAL.PS' ; runwait('lwpla'); inext 'pl'; inv -1;extdest vlbapang *********************************************************************** *********************************************************************** *********************************************************************** *********************************************************************** *********************************************************************** ** If TY or GC tables incorrect, append vlba_gains.key to end of cal file, ** then run VLOG and ANTAB. fg docal 1; flagver 1; gainuse 0; bif 0; eif 0; anten 0; solint 0 aparm 1,1,0,0,-180,180,0,0,1,0; nplot 9; dotv 1;codetyp 'p' doband -1;bparm 0;baseline 2 0; stokes 'll'; ****** identify a scan (#) on a strong calibrator that has all antennas present TIMER = scantime(##) go possm; *********************** *********************** *********************** *********************** opcode ''; anten 0;gainuse 0;calsour '' gainuse 0;calsour '' vlbapcor ************ *** If vlbapcor doesn't work, try a different scan, or use procedure at ** the bottom of this script. ******************************* tget possm;aparm 1,1,0,0,-180,180,0,0,1,0; nplot 9; dotv -1;codetyp 'a&p' baseline 2 0; stokes 'll'; runwait('possm') tget lwpla; outfil 'U:POSSM_2.PS' ;plver = 1; inv = maxtab('pl'); runwait('lwpla'); inext 'pl'; inv -1; extdest tget possm; aparm 1,1,0,0,-180,180,0,0,1,0; nplot 9;codetyp 'a&p' baseline 2 0; anten 0; doband -1; solint 0; stokes 'rl'; dotv -1; runwait('possm'); tget lwpla; plver= 1; inver=maxtab('pl'); outfil 'U:POSSM_CROSS.PS' ; runwait('lwpla'); inext 'pl'; inv -1; extdest docal 2; bpassprm 0; bpassprm(10) 2; bpassprm(11) 2; bpver 0;clr3nam gainuse 0;outver 0; bif 0; eif 0;solint -1; smodel 0; doband -1 runwait('bpass') task 'possm'; aparm 1,1,0,9,-180,180, 0,0,1,0; nplot 9;doband 1; bpver 0; dotv -1; baseline 2 0; codetyp 'a&p'; solint 0 stokes 'll'; runwait('possm') task 'lwpla';plver= 1; inver=maxtab('pl'); outfil 'U:POSSM_BPASS.PS' ; runwait('lwpla'); inext 'pl'; inv -1; extdest tget possm; aparm 1, 1, 0.4, 1.1, 0, 0, 0, 2, 1, 0 doband 1; bpver 0; timer 0; codetyp 'amp'; nplot 4; dotv -1 baseline 0; stokes ''; runwait('possm') tget lwpla;plver= 1; inver=maxtab('pl'); outfil 'U:BPASS.PS' ; runwait('lwpla'); inext 'pl'; inv -1; extdest indisk 2; for i = 1 to 200;getn i;zap;end innam 'VLBACODE'; indisk 1;getn 1; outdisk 2 outcl 'split'; dopol -1; gainuse 0; doband 1; bpver 0; docal 2; aparm 2 2 0; douvcomp false; timer 0; inv 0; stokes '' outseq 2; sources ''; runwait('split') task 'calib';default;indisk 2;refant=antnum('REFANTENNA'); outdisk 2; solint 2/60; smodel 1,0;aparm 0; aparm(7) 1e-2; aparm(1) 3; cparm 0; anten 0; docal -1; outseq 2; cmethod 'dft';solmode 'p'; outcl 'calib'; doband -1; snver 0; for i=1 to NSRC.0;getn i; inseq 2; outn inn;runwait('calib');end dowait true task 'fittp' incl 'split';inseq 2;indisk 2 outfil 'U:bsrc01.u.epoch.uvf_raw innam 'src01';go outfil 'U:bsrc02.u.epoch.uvf_raw innam 'src02';go outfil 'U:bsrc03.u.epoch.uvf_raw innam 'src03';go outfil 'U:bsrc04.u.epoch.uvf_raw innam 'src04';go outfil 'U:bsrc05.u.epoch.uvf_raw innam 'src05';go outfil 'U:bsrc06.u.epoch.uvf_raw innam 'src06';go outfil 'U:bsrc07.u.epoch.uvf_raw innam 'src07';go outfil 'U:bsrc08.u.epoch.uvf_raw innam 'src08';go outfil 'U:bsrc09.u.epoch.uvf_raw innam 'src09';go outfil 'U:bsrc10.u.epoch.uvf_raw innam 'src10';go outfil 'U:bsrc11.u.epoch.uvf_raw innam 'src11';go outfil 'U:bsrc12.u.epoch.uvf_raw innam 'src12';go outfil 'U:bsrc13.u.epoch.uvf_raw innam 'src13';go outfil 'U:bsrc14.u.epoch.uvf_raw innam 'src14';go outfil 'U:bsrc15.u.epoch.uvf_raw innam 'src15';go outfil 'U:bsrc16.u.epoch.uvf_raw innam 'src16';go outfil 'U:bsrc17.u.epoch.uvf_raw innam 'src17';go outfil 'U:bsrc18.u.epoch.uvf_raw innam 'src18';go outfil 'U:bsrc19.u.epoch.uvf_raw innam 'src19';go outfil 'U:bsrc20.u.epoch.uvf_raw innam 'src20';go outfil 'U:bsrc21.u.epoch.uvf_raw innam 'src21';go outfil 'U:bsrc22.u.epoch.uvf_raw innam 'src22';go outfil 'U:bsrc23.u.epoch.uvf_raw innam 'src23';go outfil 'U:bsrc24.u.epoch.uvf_raw innam 'src24';go outfil 'U:bsrc25.u.epoch.uvf_raw innam 'src25';go outfil 'U:bsrc26.u.epoch.uvf_raw innam 'src26';go outfil 'U:bsrc27.u.epoch.uvf_raw innam 'src27';go outfil 'U:bsrc28.u.epoch.uvf_raw innam 'src28';go outfil 'U:bsrc29.u.epoch.uvf_raw innam 'src29';go outfil 'U:bsrc30.u.epoch.uvf_raw innam 'src30';go outfil 'U:bsrc31.u.epoch.uvf_raw innam 'src31';go outfil 'U:bsrc32.u.epoch.uvf_raw innam 'src32';go incl 'calib'; outfil 'U:bsrc01.u.epoch.uvf_raw_edt innam 'src01';go outfil 'U:bsrc02.u.epoch.uvf_raw_edt innam 'src02';go outfil 'U:bsrc03.u.epoch.uvf_raw_edt innam 'src03';go outfil 'U:bsrc04.u.epoch.uvf_raw_edt innam 'src04';go outfil 'U:bsrc05.u.epoch.uvf_raw_edt innam 'src05';go outfil 'U:bsrc06.u.epoch.uvf_raw_edt innam 'src06';go outfil 'U:bsrc07.u.epoch.uvf_raw_edt innam 'src07';go outfil 'U:bsrc08.u.epoch.uvf_raw_edt innam 'src08';go outfil 'U:bsrc09.u.epoch.uvf_raw_edt innam 'src09';go outfil 'U:bsrc10.u.epoch.uvf_raw_edt innam 'src10';go outfil 'U:bsrc11.u.epoch.uvf_raw_edt innam 'src11';go outfil 'U:bsrc12.u.epoch.uvf_raw_edt innam 'src12';go outfil 'U:bsrc13.u.epoch.uvf_raw_edt innam 'src13';go outfil 'U:bsrc14.u.epoch.uvf_raw_edt innam 'src14';go outfil 'U:bsrc15.u.epoch.uvf_raw_edt innam 'src15';go outfil 'U:bsrc16.u.epoch.uvf_raw_edt innam 'src16';go outfil 'U:bsrc17.u.epoch.uvf_raw_edt innam 'src17';go outfil 'U:bsrc18.u.epoch.uvf_raw_edt innam 'src18';go outfil 'U:bsrc19.u.epoch.uvf_raw_edt innam 'src19';go outfil 'U:bsrc20.u.epoch.uvf_raw_edt innam 'src20';go outfil 'U:bsrc21.u.epoch.uvf_raw_edt innam 'src21';go outfil 'U:bsrc22.u.epoch.uvf_raw_edt innam 'src22';go outfil 'U:bsrc23.u.epoch.uvf_raw_edt innam 'src23';go outfil 'U:bsrc24.u.epoch.uvf_raw_edt innam 'src24';go outfil 'U:bsrc25.u.epoch.uvf_raw_edt innam 'src25';go outfil 'U:bsrc26.u.epoch.uvf_raw_edt innam 'src26';go outfil 'U:bsrc27.u.epoch.uvf_raw_edt innam 'src27';go outfil 'U:bsrc28.u.epoch.uvf_raw_edt innam 'src28';go outfil 'U:bsrc29.u.epoch.uvf_raw_edt innam 'src29';go outfil 'U:bsrc30.u.epoch.uvf_raw_edt innam 'src30';go outfil 'U:bsrc31.u.epoch.uvf_raw_edt innam 'src31';go outfil 'U:bsrc32.u.epoch.uvf_raw_edt innam 'src32';go * *================END OF INITIAL AIPS CALIBRATION PROCEDURES========= ***************************************************************************** *********************************************************************** *********************************************************************** *********************************************************************** * load files into difmap, use getwin to get clean windows, and flag files ** and build up best clean model, with medt, csm, acal, then use msave filename to save. insert into VLBAsumfiles (obscode, EPOCH,freq) values ('VLBACODE','epoch','15.4'); cd $U/DTERMS fg quit fg * ** Delete SPLIT and CALIB AIPS files indisk 2 for i = 1 to 200;getn i;zap;end * ** import CMAP and uvf_raw_edt files back to aips dowait true tget fitld douvcomp -1;ncount 0 outcl 'EDIT';outseq 1; outdisk 2 infil 'U:bsrc01.u.epoch.uvf_raw_edt outnam 'src01';go infil 'U:bsrc02.u.epoch.uvf_raw_edt outnam 'src02';go infil 'U:bsrc03.u.epoch.uvf_raw_edt outnam 'src03';go infil 'U:bsrc04.u.epoch.uvf_raw_edt outnam 'src04';go infil 'U:bsrc05.u.epoch.uvf_raw_edt outnam 'src05';go infil 'U:bsrc06.u.epoch.uvf_raw_edt outnam 'src06';go infil 'U:bsrc07.u.epoch.uvf_raw_edt outnam 'src07';go infil 'U:bsrc08.u.epoch.uvf_raw_edt outnam 'src08';go infil 'U:bsrc09.u.epoch.uvf_raw_edt outnam 'src09';go infil 'U:bsrc10.u.epoch.uvf_raw_edt outnam 'src10';go infil 'U:bsrc11.u.epoch.uvf_raw_edt outnam 'src11';go infil 'U:bsrc12.u.epoch.uvf_raw_edt outnam 'src12';go infil 'U:bsrc13.u.epoch.uvf_raw_edt outnam 'src13';go infil 'U:bsrc14.u.epoch.uvf_raw_edt outnam 'src14';go infil 'U:bsrc15.u.epoch.uvf_raw_edt outnam 'src15';go infil 'U:bsrc16.u.epoch.uvf_raw_edt outnam 'src16';go infil 'U:bsrc17.u.epoch.uvf_raw_edt outnam 'src17';go infil 'U:bsrc18.u.epoch.uvf_raw_edt outnam 'src18';go infil 'U:bsrc19.u.epoch.uvf_raw_edt outnam 'src19';go infil 'U:bsrc20.u.epoch.uvf_raw_edt outnam 'src20';go infil 'U:bsrc21.u.epoch.uvf_raw_edt outnam 'src21';go infil 'U:bsrc22.u.epoch.uvf_raw_edt outnam 'src22';go infil 'U:bsrc23.u.epoch.uvf_raw_edt outnam 'src23';go infil 'U:bsrc24.u.epoch.uvf_raw_edt outnam 'src24';go infil 'U:bsrc25.u.epoch.uvf_raw_edt outnam 'src25';go infil 'U:bsrc26.u.epoch.uvf_raw_edt outnam 'src26';go infil 'U:bsrc27.u.epoch.uvf_raw_edt outnam 'src27';go infil 'U:bsrc28.u.epoch.uvf_raw_edt outnam 'src28';go infil 'U:bsrc29.u.epoch.uvf_raw_edt outnam 'src29';go infil 'U:bsrc30.u.epoch.uvf_raw_edt outnam 'src30';go infil 'U:bsrc31.u.epoch.uvf_raw_edt outnam 'src31';go infil 'U:bsrc32.u.epoch.uvf_raw_edt outnam 'src32';go outcl 'CMAP';outseq 1;optype '' infil 'U:bsrc01.u.epoch.CMAP outnam 'src01';go infil 'U:bsrc02.u.epoch.CMAP outnam 'src02';go infil 'U:bsrc03.u.epoch.CMAP outnam 'src03';go infil 'U:bsrc04.u.epoch.CMAP outnam 'src04';go infil 'U:bsrc05.u.epoch.CMAP outnam 'src05';go infil 'U:bsrc06.u.epoch.CMAP outnam 'src06';go infil 'U:bsrc07.u.epoch.CMAP outnam 'src07';go infil 'U:bsrc08.u.epoch.CMAP outnam 'src08';go infil 'U:bsrc09.u.epoch.CMAP outnam 'src09';go infil 'U:bsrc10.u.epoch.CMAP outnam 'src10';go infil 'U:bsrc11.u.epoch.CMAP outnam 'src11';go infil 'U:bsrc12.u.epoch.CMAP outnam 'src12';go infil 'U:bsrc13.u.epoch.CMAP outnam 'src13';go infil 'U:bsrc14.u.epoch.CMAP outnam 'src14';go infil 'U:bsrc15.u.epoch.CMAP outnam 'src15';go infil 'U:bsrc16.u.epoch.CMAP outnam 'src16';go infil 'U:bsrc17.u.epoch.CMAP outnam 'src17';go infil 'U:bsrc18.u.epoch.CMAP outnam 'src18';go infil 'U:bsrc19.u.epoch.CMAP outnam 'src19';go infil 'U:bsrc20.u.epoch.CMAP outnam 'src20';go infil 'U:bsrc21.u.epoch.CMAP outnam 'src21';go infil 'U:bsrc22.u.epoch.CMAP outnam 'src22';go infil 'U:bsrc23.u.epoch.CMAP outnam 'src23';go infil 'U:bsrc24.u.epoch.CMAP outnam 'src24';go infil 'U:bsrc25.u.epoch.CMAP outnam 'src25';go infil 'U:bsrc26.u.epoch.CMAP outnam 'src26';go infil 'U:bsrc27.u.epoch.CMAP outnam 'src27';go infil 'U:bsrc28.u.epoch.CMAP outnam 'src28';go infil 'U:bsrc29.u.epoch.CMAP outnam 'src29';go infil 'U:bsrc30.u.epoch.CMAP outnam 'src30';go infil 'U:bsrc31.u.epoch.CMAP outnam 'src31';go infil 'U:bsrc32.u.epoch.CMAP outnam 'src32';go * ** run CALIB on EDIT files using CMAPs, 30 sec A&P -> create APCAL.1 task 'calib';default; outcl 'apcal';indisk 2;in2disk 2;outdisk 2; in2cl 'cmap'; inseq 1; inver 1 cmeth 'dft';refant=antnum('REFANTENNA');solint 0.5;aparm 4,0;solmode 'a&p' in2seq 1; for i =1 to NSRC.0;getn i;in2nam= innam; outn=innam; go calib;end **run CCEDT on CMAPs task 'ccedt'; default;cparm 0,0,4, 0;inv 1; outv 2;indisk 2 for i = (NSRC.0+1) to (NSRC.0*2); getn i; go ccedt; end ******************run LPCAL on APCALs using CMAPs task 'lpcal';default; cparm 0; solint 0; in2ver 2;indisk 2;in2disk 2 in2clas 'cmap' for i = (2*NSRC.0+1) to 3*NSRC.0;getn i;in2nam = innam; go lpcal;end ************ run TBOUT on APCALs to write AN tables to disk task 'tbout';default; docrt 1200; inext 'an'; inclas 'apcal' outfil 'U:DTERMS/src01.AN innam 'src01';go outfil 'U:DTERMS/src02.AN innam 'src02';go outfil 'U:DTERMS/src03.AN innam 'src03';go outfil 'U:DTERMS/src04.AN innam 'src04';go outfil 'U:DTERMS/src05.AN innam 'src05';go outfil 'U:DTERMS/src06.AN innam 'src06';go outfil 'U:DTERMS/src07.AN innam 'src07';go outfil 'U:DTERMS/src08.AN innam 'src08';go outfil 'U:DTERMS/src09.AN innam 'src09';go outfil 'U:DTERMS/src10.AN innam 'src10';go outfil 'U:DTERMS/src11.AN innam 'src11';go outfil 'U:DTERMS/src12.AN innam 'src12';go outfil 'U:DTERMS/src13.AN innam 'src13';go outfil 'U:DTERMS/src14.AN innam 'src14';go outfil 'U:DTERMS/src15.AN innam 'src15';go outfil 'U:DTERMS/src16.AN innam 'src16';go outfil 'U:DTERMS/src17.AN innam 'src17';go outfil 'U:DTERMS/src18.AN innam 'src18';go outfil 'U:DTERMS/src19.AN innam 'src19';go outfil 'U:DTERMS/src20.AN innam 'src20';go outfil 'U:DTERMS/src21.AN innam 'src21';go outfil 'U:DTERMS/src22.AN innam 'src22';go outfil 'U:DTERMS/src23.AN innam 'src23';go outfil 'U:DTERMS/src24.AN innam 'src24';go outfil 'U:DTERMS/src25.AN innam 'src25';go outfil 'U:DTERMS/src26.AN innam 'src26';go outfil 'U:DTERMS/src27.AN innam 'src27';go outfil 'U:DTERMS/src28.AN innam 'src28';go outfil 'U:DTERMS/src29.AN innam 'src29';go outfil 'U:DTERMS/src30.AN innam 'src30';go outfil 'U:DTERMS/src31.AN innam 'src31';go outfil 'U:DTERMS/src32.AN innam 'src32';go * *********************************************************************** *********************************************************************** *********************************************************************** ********************************************************************** ********************************************************************** ** Flag out cross hands if necessary *task 'uvflg'; stokes 'cros'; infil '';aparm 0; *anten 1 2 3 4 5 6 7 8 9 10 *for i = 2*NSRC.0+1 to 3*NSRC.0; getn i;runwait('uvflg');end * ************ Go back to Linux prompt ********** *cd $U/DTERMS *process_AN.pl *.AN **** move bad dterm source AN files to OUT/DTERMS/NOTUSED directory ** If expt is missing an antenna then you need to edit the DTERMS file *mv DTERMS FINAL.DTERMS *cp $U/DTERMS/FINAL.DTERMS ~/a/MOJAVE/evpa/epoch.u.dterms.REFANTENNA.-49.0 *scp $U/DTERMS/FINAL.DTERMS mlister@titania:/project/prospero/mlister/MOJAVE/evpa/epoch.u.dterms.REFANTENNA.-49.0 *cd /project/prospero/mlister/MOJAVE/evpa * *new_process_EVPA.pl *.u.* * scp /project/prospero/mlister/MOJAVE/evpa/epoch* mlister@polaris.cv.nrao.edu:/home/www.cv.nrao.edu/content/2cmVLBA/dterms *** establishes EVPA correction and adds ref antenna to mysql database ************ Go back to AIPS * cd $U fg * **extdest AN 1 of all APCAL files indisk 2 for i = 1+(2*NSRC.0) to 3*NSRC.0; getn i;inext 'an';inv 1;extdest;end infil 'U:DTERMS/FINAL.DTERMS';outdisk 2 for i = 1+(2*NSRC.0) to 3*NSRC.0;geton i;go tbin;end task 'split';default; docal -1; dopol 3; outcl 'PSPLIT'; doband -1;douvcomp -1; outdisk 2;indisk 2 for i = 1+(2*NSRC.0) to 3*NSRC.0;getn i;go split; end task 'calib';default; outcl 'PINAL'; inver 1;indisk 2;outdisk 2; in2disk 2;cmeth 'dft' in2cl 'cmap'; solmode 'a&p';solint 0.5;refant=antnum('REFANTENNA');aparm 4,0 for i = 1+3*NSRC.0 to 4*NSRC.0;getn i;outn innam;in2nam innam; go calib;end task 'multi'; default; outcl 'multi' indisk 2; outdisk 2;aparm 0; for i = 1+(4*NSRC.0) to 5*NSRC.0;getn i;outnam innam; outseq 1;inseq 1;go multi;end *****clcorprm is twice the evpa correction task 'clcor';default; indisk 2; OPCODE='PHAS' ;stokes 'l' *vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv clcorprm -98,-98,-98, -98 *^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ for i = 1+(5*NSRC.0) to 6*NSRC.0;getn i;go clcor;end task 'indxr';default;cparm 0; indisk 2; for i = 1+(5*NSRC.0) to 6*NSRC.0;getn i;go indxr;end task 'split';default; outcl 'uvf';dopol = -1 gainuse 2;docal 1;douvcomp -1; indisk 2;outdisk 2 for i = 1+(5*NSRC.0) to 6*NSRC.0;getn i;go split;end task 'fittp';default; incl 'uvf';inseq 1;indisk 2 outfil 'U:bsrc01.u.epoch.uvf innam 'src01';go outfil 'U:bsrc02.u.epoch.uvf innam 'src02';go outfil 'U:bsrc03.u.epoch.uvf innam 'src03';go outfil 'U:bsrc04.u.epoch.uvf innam 'src04';go outfil 'U:bsrc05.u.epoch.uvf innam 'src05';go outfil 'U:bsrc06.u.epoch.uvf innam 'src06';go outfil 'U:bsrc07.u.epoch.uvf innam 'src07';go outfil 'U:bsrc08.u.epoch.uvf innam 'src08';go outfil 'U:bsrc09.u.epoch.uvf innam 'src09';go outfil 'U:bsrc10.u.epoch.uvf innam 'src10';go outfil 'U:bsrc11.u.epoch.uvf innam 'src11';go outfil 'U:bsrc12.u.epoch.uvf innam 'src12';go outfil 'U:bsrc13.u.epoch.uvf innam 'src13';go outfil 'U:bsrc14.u.epoch.uvf innam 'src14';go outfil 'U:bsrc15.u.epoch.uvf innam 'src15';go outfil 'U:bsrc16.u.epoch.uvf innam 'src16';go outfil 'U:bsrc17.u.epoch.uvf innam 'src17';go outfil 'U:bsrc18.u.epoch.uvf innam 'src18';go outfil 'U:bsrc19.u.epoch.uvf innam 'src19';go outfil 'U:bsrc20.u.epoch.uvf innam 'src20';go outfil 'U:bsrc21.u.epoch.uvf innam 'src21';go outfil 'U:bsrc22.u.epoch.uvf innam 'src22';go outfil 'U:bsrc23.u.epoch.uvf innam 'src23';go outfil 'U:bsrc24.u.epoch.uvf innam 'src24';go outfil 'U:bsrc25.u.epoch.uvf innam 'src25';go outfil 'U:bsrc26.u.epoch.uvf innam 'src26';go outfil 'U:bsrc27.u.epoch.uvf innam 'src27';go outfil 'U:bsrc28.u.epoch.uvf innam 'src28';go outfil 'U:bsrc29.u.epoch.uvf innam 'src29';go outfil 'U:bsrc30.u.epoch.uvf innam 'src30';go outfil 'U:bsrc31.u.epoch.uvf innam 'src31';go outfil 'U:bsrc32.u.epoch.uvf innam 'src32';go ****************************************************** *********************************************************************** ***** exit to linux prompt *********************************************************************** ****************************************************** mv -f $U/bsrc01* $U/bsrc01/epoch mv -f $U/bsrc02* $U/bsrc02/epoch mv -f $U/bsrc03* $U/bsrc03/epoch mv -f $U/bsrc04* $U/bsrc04/epoch mv -f $U/bsrc05* $U/bsrc05/epoch mv -f $U/bsrc06* $U/bsrc06/epoch mv -f $U/bsrc07* $U/bsrc07/epoch mv -f $U/bsrc08* $U/bsrc08/epoch mv -f $U/bsrc09* $U/bsrc09/epoch mv -f $U/bsrc10* $U/bsrc10/epoch mv -f $U/bsrc11* $U/bsrc11/epoch mv -f $U/bsrc12* $U/bsrc12/epoch mv -f $U/bsrc13* $U/bsrc13/epoch mv -f $U/bsrc14* $U/bsrc14/epoch mv -f $U/bsrc15* $U/bsrc15/epoch mv -f $U/bsrc16* $U/bsrc16/epoch mv -f $U/bsrc17* $U/bsrc17/epoch mv -f $U/bsrc18* $U/bsrc18/epoch mv -f $U/bsrc19* $U/bsrc19/epoch mv -f $U/bsrc20* $U/bsrc20/epoch mv -f $U/bsrc21* $U/bsrc21/epoch mv -f $U/bsrc22* $U/bsrc22/epoch mv -f $U/bsrc23* $U/bsrc23/epoch mv -f $U/bsrc24* $U/bsrc24/epoch mv -f $U/bsrc25* $U/bsrc25/epoch mv -f $U/bsrc26* $U/bsrc26/epoch mv -f $U/bsrc27* $U/bsrc27/epoch mv -f $U/bsrc28* $U/bsrc28/epoch mv -f $U/bsrc29* $U/bsrc29/epoch mv -f $U/bsrc30* $U/bsrc30/epoch mv -f $U/bsrc31* $U/bsrc31/epoch mv -f $U/bsrc32* $U/bsrc32/epoch * cd $U make_clean_sources.pl *** Now clean the uvf files to make FITS files ************ difmap @clean_sources_n quit difmap @clean_sources_t quit ******************* insert into evpacorr (obscode,EPOCH,corr,refant,frequency) values ('VLBACODE','epoch','-49','REFANTENNA','15.4'); **** Now produce ps images *********** *cd $U *blmap.pl * ****(change frange, cbase, pbase, gap parameters as necessary for each source) * *./mmap_script * *** Now produce tapered images: *blmap.pl *./mmap_script * * difmap @radplot_sources convert2jpg.pl chmod ogu+x make_jpg ./make_jpg ****************; *******************; update VLBAsumfiles set comment='in archive' where freq like '15%' and obscode='VLBACODE'; *********************; *********************; polstats.pl * *** Go back to AIPS and output TASAV file to disk fg indisk 1;outdisk 1; outfil 'U:VLBACODE_U.TASAV getn 1; geton 1; outcl 'tasav'; * * go tasav;zap;getn 2;go fittp;zap;exit * mkdir $U/AIPS mv *TXT $U/AIPS mv *.PS $U/AIPS mv *TASAV $U/AIPS gzip $U/AIPS/* gzip $U/*/*/*CMAP gzip $U/*/*/*CMOD gzip $U/*/*/*t.mod gzip $U/*/*/*ps gzip $U/*/*/*.fits gzip $U/*/*/*uvf_raw_edt gzip $U/*/*/*uvf_raw rm -rf $U/difma* rm -rf $U/win rm -rf $U/outwin rm -rf $U/*MKIII rm -rf $U/*FLAG rm -rf $U/*SCAN rm -rf $U/*PCAL rm -rf $U/*TSYS rm -rf $U/*WX rm -rf $U/*FITS rm -rf $U/*fits rm -rf f? rm -rf f?? export POLARIS=mlister@polaris.cv.nrao.edu:/home/www.cv.nrao.edu/content/2cmVLBA/data setenv POLARIS mlister@polaris.cv.nrao.edu:/home/www.cv.nrao.edu/content/2cmVLBA/data setenv TITANIA /project/titania/CV_archive/data cp -r $U/bsrc01 $TITANIA cp -r $U/bsrc02 $TITANIA cp -r $U/bsrc03 $TITANIA cp -r $U/bsrc04 $TITANIA cp -r $U/bsrc05 $TITANIA cp -r $U/bsrc06 $TITANIA cp -r $U/bsrc07 $TITANIA cp -r $U/bsrc08 $TITANIA cp -r $U/bsrc09 $TITANIA cp -r $U/bsrc10 $TITANIA cp -r $U/bsrc11 $TITANIA cp -r $U/bsrc12 $TITANIA cp -r $U/bsrc13 $TITANIA cp -r $U/bsrc14 $TITANIA cp -r $U/bsrc15 $TITANIA cp -r $U/bsrc16 $TITANIA cp -r $U/bsrc17 $TITANIA cp -r $U/bsrc18 $TITANIA cp -r $U/bsrc19 $TITANIA cp -r $U/bsrc20 $TITANIA cp -r $U/bsrc21 $TITANIA cp -r $U/bsrc22 $TITANIA cp -r $U/bsrc23 $TITANIA cp -r $U/bsrc24 $TITANIA cp -r $U/bsrc25 $TITANIA cp -r $U/bsrc26 $TITANIA cp -r $U/bsrc27 $TITANIA cp -r $U/bsrc28 $TITANIA cp -r $U/bsrc29 $TITANIA cp -r $U/bsrc30 $TITANIA cp -r $U/bsrc31 $TITANIA cp -r $U/bsrc32 $TITANIA **************************** **************************** **************************** tar -cvf tar.tar 0* 1* 2* scp -pr $U/tar.tar $POLARIS * polaris cd $WWW/data tar -xvf tar.tar rm -rf tar.tar ** If non-MOJAVE EPOCH, do the following: chmod o+r */*/*epoch* exit * rm -rf $U/tar.tar ********************************************************************** ********************************************************************** ********************************************************************** ********************************************************************** ********************************************************************** ********************************************************************** * * *** burn cdrom of $U directory using xcdroast * cd /project/prospero/mlister/public_html/MOJAVE/ nrao_gen.pl mv -f /project/prospero/mlister/public_html/MOJAVE/sourcepages/* /project/astro/www/MOJAVE/sourcepages * *** update the data.html page on purdue web site using mozilla -edit /project/astro/www/MOJAVE/data.html * *** delete files from AIPS * ------------------------------------------------------------------ ***** END OF CALIBRATION/ARCHIVING ********* ------------------------------------------------------------------- * ------------------------------------------------------------------- SN CL 0 1 vlbaload 1 2 ACCOR 2 APCAL 3 3 SNEDT of APCAL 4 TACOP (of CL 2) 5 vlbapang 4 6 vlbapcor 7 clcor * Misc files - found in directory VLBACODE ---------------------------------------- clean_sources = difmap script for creating final images clean_stokes = script for cleaning all 4 stokes parameters -- called by clean_sources final_clean = script called by clean_stokes difmap.log = logfile of creating deep cleaned images final.images.dat = datafile with statistics about final images * * Data files - found in directory VLBACODE/SOURCENAME ----------------------------------------------- sourcename.CMAP = initial FITS image used for AIPS selfcals sourcename.CMOD = initial DIFMAP model file sourcename.SPLIT = initial split *sourcename.CALIB = split.1 with 2 second phase self-cal applied (file not kept) sourcename.EDIT = calib.1 that has been edited, uvaver 30,true *sourcename.APCAL = edit.1 with 30 sec A&P selfcal (file not kept) *sourcename.PSPLIT = apcal.1 with dterms applied (file not kept) sourcename.uvf = final calibrated file, with EVPAs correctly rotated sourcename.IMAP.gz = gzipped deep cleaned I image sourcename.imod = modelfile for above image sourcename.QMAP.gz = gzipped deep cleaned Q image sourcename.qmod = modelfile for above image sourcename.UMAP.gz = gzipped deep cleaned U image sourcename.umod = modelfile for above image sourcename.ps = multi-image plot of I, P, & Fractional P sourcename.wf.ps = wide-field plot of I difmap.log = log file(s) from difmap mmap.opt = options file for mmap plotting script * *AIPS files that are not saved * DTERM files - found in directory VLBACODE/DTERMS ----------------------------------------------- sourcename.AN = DTERM solutions from LPCAL FINAL.DTERMS = final resulting DTERM table from process_AN.pl script (renamed from "DTERMS" which is output by process_AN.pl) * Output files - found in directory VLBACODE/AIPS ----------------------------------------------- ACCOR.PS = Amplitude solutions from ACCOR APCAL.PS = Amplitude solutions from APCAL BPASS.PS = Bandpass solutions from BPASS DELAYS.PS = Delay solutions from FRING OPACITY.PS = Opacity plots from APCAL PCOR.PS = Delay corrections from PCOR POSSM_1.PS = POSSM plot prior to PCOR POSSM_2.PS = POSSM plot with PCOR corrections applied POSSM_BPASS.PS = POSSM plot with BPASS corrections applied POSSM_CROSS.PS = POSSM plot with CRSFRING solutions applied RATES.PS = Rate solutions from FRING TSYS.PS = Plot of TY table SUTABLE.TXT = SU table with SETJY values VLBASUMM.TXT = VLBA summary file with scan times, antenna numbers exptcode.TASAV.gz = AIPS tables, AN table 2 contains D-terms exptcode.SCRIPT = script file * *** If vlbapcor doesn't work, try a different scan, otherwise use following: * outcl 'TASAV'; outn innam; outdi indi;runwait('tasav'); task 'vlog'; infil 'U:vlbacodecal.vlba fqtol -1;prtlev 0;outfil 'U:VLBACODE'; runwait('vlog'); * clrst; inext 'pc'; inv 1; extdest; task 'pclod'; outver 0; infil 'U:VLBACODE.PCAL'; runwait('pclod'); * task 'fring' clr2nam;dofit 0 calsour 'src01','' timer=scantime(##) docal 1; gainuse 0;refant=antnum('REFANTENNA') solint 0; aparm 0; dparm 0; aparm(1) 2 aparm(6) 1 dparm(1) 1 dparm(8) 1 snver 0; anten 0;gainu 0 runwait('fring') * task 'clcal'; timer 0 calsour 'src01','' source ''; interp '';snver=maxtab('SN'); refant=antnum('REFANTENNA') gainver 0; gainuse 0;opcode '' runwait('clcal'); * TIMER=0 SOLINT 2 REFANT = antnum('FD') BADDISK 0 calsour '' dparm(4) 2 source(1) ='src01' source(2) ='src02' source(3) ='src03' source(4) ='src04' source(5) ='src05' source(6) ='src06' SOURCE(7) = 'src07' SOURCE(8) = 'src08' SOURCE(9) = 'src09' SOURCE(10) = 'src10' SOURCE(11) = 'src11' SOURCE(12) = 'src12' SOURCE(13) = 'src13' SOURCE(14) = 'src14' SOURCE(15) = 'src15' SOURCE(16) = 'src16' SOURCE(17) = 'src17' SOURCE(18) = 'src18' SOURCE(19) = 'src19' SOURCE(20) = 'src20' SOURCE(21) = 'src21' SOURCE(22) = 'src22' SOURCE(23) = 'src23' SOURCE(24) = 'src24' SOURCE(25) = 'src25' docrt -1;inext 'su';inv 1;outprint 'U:SUTABLE.TXT';box 2 5;runwait('prtab') timer 0; solint = 2.0; RUNWAIT ('ACCOR') SAMPTYPE = 'MWF'; DOBLANK = -1; SMOTYPE = 'AMPL' CPARM = 0.5, 0, 0, 0, 0, 0.02; INVERS = MAXTAB('SN') OUTVERS = INVERS + 1; RUNWAIT ('SNSMO') INEXT = 'SN'; EXTDEST; VBA_SN = INVERS vlbafrng