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