h=plotdriftcontours(ring) load S25ERRseed1.mat
0001 function h=plotContours(ring,indices,colorplotstr) 0002 %h=plotdriftcontours(ring) 0003 %load S25ERRseed1.mat 0004 0005 [dat,prm]=atx(ring); 0006 sig=cat(3,dat.beam66); 0007 0008 sigxx=squeeze(sig(1,1,:)); 0009 sigxy=squeeze(sig(1,3,:)); 0010 sigyy=squeeze(sig(3,3,:)); 0011 0012 %dr278=findcells(ring,'FamName','DR_278'); 0013 0014 hold on 0015 for j=indices 0016 [x,y]=atmakeXYProjectionEllipse(sigxx(j),sigyy(j),sigxy(j)); 0017 plot(x,y,colorplotstr) 0018 end