AIPS version = /home/auto/aips/31DEC07/LINUX/SYSTEM/VLA_PIPE/aips Prep version = Dec 7 2007 /home/vlbacald-a/vlbacald/bin/AVLAPREP.001 ------------- Pops Prior Date Time Task Messages for user23917 1 0 01-JUN-2010 08:31:25 AIPS version'ABIN' 1 0 01-JUN-2010 08:31:25 AIPS run avlaprep 1 2 01-JUN-2010 08:31:25 AIPS Found in Version=ABIN 1 2 01-JUN-2010 08:31:25 AIPS $ Lorant Sjouwerman, versions 1.. 2005, for NRAO 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS $$ General naming utilities - default to data disk 1 (safest of all) 1 2 01-JUN-2010 08:31:25 AIPS procedure pipedef1 1 2 01-JUN-2010 08:31:25 AIPS string*12 pnam 1 2 01-JUN-2010 08:31:25 AIPS string*6 pcls 1 2 01-JUN-2010 08:31:25 AIPS scalar pseq,pdsk 1 2 01-JUN-2010 08:31:25 AIPS pdsk=1 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS pipedef1 1 2 01-JUN-2010 08:31:25 AIPS procedure namepipe 1 2 01-JUN-2010 08:31:25 AIPS pnam=inname;pcls=inclass;pseq=inseq;indisk=pdsk 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS procedure pipename 1 2 01-JUN-2010 08:31:25 AIPS inname=pnam;inclass=pcls;inseq=pseq;indisk=pdsk 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure checkcat 1 2 01-JUN-2010 08:31:25 AIPS $LOS - count and make sure catalog starts at one (and not more than 10 files) 1 2 01-JUN-2010 08:31:25 AIPS scalar numcat 1 2 01-JUN-2010 08:31:25 AIPS clrname;indisk=pdsk;chkname;numcat=1-error 1 2 01-JUN-2010 08:31:25 AIPS $ for later: find out how to deal with CH0/LINE data in MOVE 1 2 01-JUN-2010 08:31:25 AIPS if(numcat>0)then;recat;end;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS return numcat 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure gx 1 2 01-JUN-2010 08:31:25 AIPS $ inputs 1 2 01-JUN-2010 08:31:25 AIPS go;wait 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure pre_fillm 1 2 01-JUN-2010 08:31:25 AIPS $LOS - sets VLA archive defaults for fillm 1 2 01-JUN-2010 08:31:25 AIPS default fillm;outdisk=pdsk;doall=1;douvcomp=-1;cparm(8)=0.05;bparm(10)=0.75 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure fname 1 2 01-JUN-2010 08:31:25 AIPS return inname!!'.'!!inclass!!'.'!!char(inseq) 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure shortname 1 2 01-JUN-2010 08:31:25 AIPS scalar idx,nsrc,ll,ii,spc 1 2 01-JUN-2010 08:31:25 AIPS string*12 newstrng 1 2 01-JUN-2010 08:31:25 AIPS keyword'NUM ROW';inext'SU';getthead;nsrc=keyvalue(1);keyword'';keyvalue 0 1 2 01-JUN-2010 08:31:25 AIPS for idx=1:nsrc;pixxy=idx,2,1;tabget;ll=length(keystrng);newstrng'' 1 2 01-JUN-2010 08:31:25 AIPS if(ll>12)then;spc=0;ll=1 1 2 01-JUN-2010 08:31:25 AIPS while(spc=0&ll<17);substr(newstrng,ll,ll)=substr(keystrng,ll,ll) 1 2 01-JUN-2010 08:31:25 AIPS if(substr(keystrng,ll,ll+1)=' ')then 1 2 01-JUN-2010 08:31:25 AIPS if(spc=0)then;spc=1;substr(newstrng,ll,ll)='@';end 1 2 01-JUN-2010 08:31:25 AIPS end;ll=ll+1;clrtemp;end;ii=ll 1 2 01-JUN-2010 08:31:25 AIPS while(ll<17) 1 2 01-JUN-2010 08:31:25 AIPS if(substr(keystrng,ll,ll)<>' ')then 1 2 01-JUN-2010 08:31:25 AIPS substr(newstrng,ii,ii)=substr(keystrng,ll,ll);ii=ii+1 1 2 01-JUN-2010 08:31:25 AIPS end;ll=ll+1;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS end;ll=length(newstrng);nsrc=min(12,ll) $ for safety during tests 1 2 01-JUN-2010 08:31:25 AIPS if(ll>12)then;keystrng=newstrng;ll=ll+1;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS for ii=1:nsrc;substr(newstrng,ii,ii)=substr(keystrng,ll-ii,ll-ii);end 1 2 01-JUN-2010 08:31:25 AIPS end;task'TABED';clroname;inext'SU';optype'REPL';keyword'';keyvalue 0 1 2 01-JUN-2010 08:31:25 AIPS aparm 2 1;bcount=idx;ecount=idx;keystrng=newstrng;go;wait;end;end;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure del 1 2 01-JUN-2010 08:31:25 AIPS error=1;zap 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure numants 1 2 01-JUN-2010 08:31:25 AIPS scalar numarow,numbout,numbant,numbase,idx;string*16 nameant 1 2 01-JUN-2010 08:31:25 AIPS numbout=0;keyword'NUM ROW';inext'AN';invers=0;getthead;numarow=keyvalue(1) 1 2 01-JUN-2010 08:31:25 AIPS for numbant=1:numarow 1 2 01-JUN-2010 08:31:25 AIPS pixxy=numbant,1,0;tabget;nameant=keystrng;idx=1;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS while(substr(nameant,idx,idx)<>':');idx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(substr(nameant,idx+1,idx+3)='OUT')then;numbout=numbout+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(substr(nameant,idx+2,idx+4)='OUT')then;numbout=numbout+1;end 1 2 01-JUN-2010 08:31:25 AIPS end;numbase=(numarow-numbout);clrtemp 1 2 01-JUN-2010 08:31:25 AIPS $ type char(numarow-numbout)!!' ANTENNAS',char(numbase)!!' BASELINES' 1 2 01-JUN-2010 08:31:25 AIPS return numbase 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure skipdata 1 2 01-JUN-2010 08:31:25 AIPS $LOS - delete unwanted data sets - currently all but continuum X and C bands 1 2 01-JUN-2010 08:31:25 AIPS $ if(inclass='LINE')then;del;inclass'CH 0';del;type'ZAP LINE :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS if((inclass='CH 0')!(inclass='LINE'))then;del;type'ZAP LINE :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS if(inclass='4 BAND')then;del;type'ZAP 4-BAND :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS if(inclass='P BAND')then;del;type'ZAP P-BAND :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS $ if(inclass='L BAND')then;del;type'ZAP L-BAND :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS $ if(inclass='U BAND')then;del;type'ZAP U-BAND :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS $ if(inclass='K BAND')then;del;type'ZAP K-BAND :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS $ if(inclass='Q BAND')then;del;type'ZAP Q-BAND :'fname;end;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS if(error<1)then $ file not zapped, check for single dish vlbi; #vis=0 1 2 01-JUN-2010 08:31:25 AIPS task'TABGET';userid 0;inext'NX';invers=0;keyvalue=0;keystrng'' 1 2 01-JUN-2010 08:31:25 AIPS keyword'NUM ROW';getthead;pixxy=keyvalue(1),6,0;tabget;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS if(keyvalue(1)=0)then;del;type'ZAP S-DISH :'fname;end 1 2 01-JUN-2010 08:31:25 AIPS if(error<1)then;idx=numants $ file not zapped, check for # ants (<9) 1 2 01-JUN-2010 08:31:25 AIPS if(idx<9)then;del;type'ZAP MINI-ARRAY ('!!char(idx)!!') :'!!fname 1 2 01-JUN-2010 08:31:25 AIPS end;clrtemp;end;end 1 2 01-JUN-2010 08:31:25 AIPS $also zap if solar/planetary/pulsar mode - su table proper motion?? 1 2 01-JUN-2010 08:31:25 AIPS $ and what about PT-link? antenna VPT is not OUT, later also EVLA names 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure checkids 1 2 01-JUN-2010 08:31:25 AIPS $LOS - account for fast switching source names - check on pos, qual, calco 1 2 01-JUN-2010 08:31:25 AIPS scalar ra1a,ra1b,ra2a,ra2b,dc1a,dc1b,dc2a,dc2b,n,m,l,idx 1 2 01-JUN-2010 08:31:25 AIPS scalar q1,q2,k,x,y1,y2 1 2 01-JUN-2010 08:31:25 AIPS string*1 co1,co2 1 2 01-JUN-2010 08:31:25 AIPS string*16 n1,n2 1 2 01-JUN-2010 08:31:25 AIPS array alist(100),blist(100) 1 2 01-JUN-2010 08:31:25 AIPS inext'SU';invers 0;keyword'NUM ROW';keyvalue 0;keystrng'';getthead 1 2 01-JUN-2010 08:31:25 AIPS n=keyvalue(1);idx=1;aparm 0;bparm 0;alist 0;blist 0 1 2 01-JUN-2010 08:31:25 AIPS if (n > 1) then 1 2 01-JUN-2010 08:31:25 AIPS for m=1:(n-1);clrtemp 1 2 01-JUN-2010 08:31:25 AIPS pixxy m,11;inext'SU';tabget;ra1a keyvalue(1);ra1b keyvalue(2) 1 2 01-JUN-2010 08:31:25 AIPS for l=(m+1):n;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS pixxy l,11;tabget;ra2a keyvalue(1);ra2b keyvalue(2) 1 2 01-JUN-2010 08:31:25 AIPS if (abs(ra1a-ra2a)=0) then 1 2 01-JUN-2010 08:31:25 AIPS if ( abs((ra1b-ra2b)*3600*1000) < 3 ) then $ RA within 3 mas! 1 2 01-JUN-2010 08:31:25 AIPS pixxy m,12;tabget;dc1a keyvalue(1);dc1b keyvalue(2) 1 2 01-JUN-2010 08:31:25 AIPS pixxy l,12;tabget;dc2a keyvalue(1);dc2b keyvalue(2) 1 2 01-JUN-2010 08:31:25 AIPS if (abs(dc1a-dc2a)=0) then 1 2 01-JUN-2010 08:31:25 AIPS if ( abs((dc1b-dc2b)*3600*1000) < 3 ) then $ DC within 3 mas! 1 2 01-JUN-2010 08:31:25 AIPS pixxy m,3;tabget;q1 keyvalue(1);pixxy l,3;tabget;q2 keyvalue(1) 1 2 01-JUN-2010 08:31:25 AIPS $if one has zero vis then next checks not needed and declare/fix same source 1 2 01-JUN-2010 08:31:25 AIPS inext'NX';keyword'NUM ROW';getthead;k=keyvalue(1);y1=0;y2=0 1 2 01-JUN-2010 08:31:25 AIPS for x=1:k;pixxy x,3;tabget 1 2 01-JUN-2010 08:31:25 AIPS if(keyval(1)=m)then;y1=y1+1;else; 1 2 01-JUN-2010 08:31:25 AIPS if(keyval(1)=l)then;y2=y2+1;end;end 1 2 01-JUN-2010 08:31:25 AIPS end $ if ((y1=0)!(y2=0)) then one has 0 visibilities: do dsorc 1 2 01-JUN-2010 08:31:25 AIPS if (q1=q2)!((y1=0)!(y2=0)) then ;inext'SU' $ same calcode? 1 2 01-JUN-2010 08:31:25 AIPS pixxy m,4;tabget;co1 keystrng;pixxy l,4;tabget;co2 keystrng 1 2 01-JUN-2010 08:31:25 AIPS if (co1=co2)!((y1=0)!(y2=0)) then 1 2 01-JUN-2010 08:31:25 AIPS pixxy m,1;tabget;q1 keyvalue(1) 1 2 01-JUN-2010 08:31:25 AIPS pixxy l,1;tabget;q2 keyvalue(1) 1 2 01-JUN-2010 08:31:25 AIPS pixxy m,2;tabget;n1 keystrng;pixxy l,2;tabget;n2 keystrng 1 2 01-JUN-2010 08:31:25 AIPS if (idx>100) then; 1 2 01-JUN-2010 08:31:25 AIPS type 'TOO MANY SOURCES TO DO DSORC - DO BY HAND!' 1 2 01-JUN-2010 08:31:25 AIPS type 'NEXT =.'!!char(idx)!!' DO:'!!char(l)!!'.&.'!!char(m) 1 2 01-JUN-2010 08:31:25 AIPS else 1 2 01-JUN-2010 08:31:25 AIPS if (length(n1)1) then;task'DSORC';n=0 1 2 01-JUN-2010 08:31:25 AIPS outname inname;outclass inclass;outseq inseq;outdisk indisk 1 2 01-JUN-2010 08:31:25 AIPS while (n100)then;type'TOO MANY SOURCES FOR DSORC';error 1;end;clrtemp 1 2 01-JUN-2010 08:31:25 AIPS end 1 2 01-JUN-2010 08:31:25 AIPS finish 1 2 01-JUN-2010 08:31:25 AIPS 1 2 01-JUN-2010 08:31:25 AIPS procedure checkcal 1 2 01-JUN-2010 08:31:25 AIPS $LOS Currently (Sep 22, 2005) can/will do 3C286 and 3C48 at X and C 1 2 01-JUN-2010 08:31:25 AIPS $LOS (and for now skip 3C286/3C48/3C138/3C147 Q,K,U and 3C48 L) 1 2 01-JUN-2010 08:31:25 AIPS scalar idx,jdx 1 2 01-JUN-2010 08:31:25 AIPS string*16 sname 1 2 01-JUN-2010 08:31:25 AIPS task'TABGET';userid 0;inext'SU';invers=0;keyvalue=0;keystrng'';error=1 1 2 01-JUN-2010 08:31:25 AIPS keyword'NUM ROW';getthead;idx=keyvalue(1);clrtemp 1 2 01-JUN-2010 08:31:25 AIPS for jdx=1:idx 1 2 01-JUN-2010 08:31:25 AIPS pixxy=jdx,2,0;tabget;sname=substr(keystrng,1,length(keystrng)) 1 2 01-JUN-2010 08:31:25 AIPS if(sname='3C286')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='3C48')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='3C138')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='3C147')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='3C295')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(error>0)then $ try B1950 name 1 2 01-JUN-2010 08:31:25 AIPS if(sname='1328+307')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='0134+329')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='0518+165')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='0538+498')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='1409+524')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(error>0)then $ try J2000 1 2 01-JUN-2010 08:31:25 AIPS if(sname='1331+305')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='0137+331')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='0521+166')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='0542+498')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(sname='1411+522')then;error=-1;jdx=idx+1;end 1 2 01-JUN-2010 08:31:25 AIPS if(error>0)then $ try J2000 name 1 2 01-JUN-2010 08:31:25 AIPS if(sname='J1331+3030')then;error=-1;jdx=idx+1;en