%load amp_data1
%load hdata
h=figure(1);
%set(h,'Position', [79   588   229   342])

set(h,'Position', [200   300   250   400])


a=[0.6793    0.2113
    1.0000   -0.2146];

c1=(10)*[-1:.4:1];
c1=[c1 c1(5:-1:2)];
c1=[c1 c1 c1 c1 c1 c1 c1 c1 c1 c1 c1 c1];
c2a=c1;
c2=(10)*[-1:(1/6):1];
c2=[c2 c2(12:-1:2)];
c2=[c2 c2 c2 c2 c2];
c1=c2;
c2=c2a;
b=a;
a=(1/max(b(:,1)))*a;
rdata=(1/sqrt(2))*[1 sqrt(2) 1 0 1];
tdata=[0 pi/4 2*pi/4 0 pi/4];

   axis image
   axis off
   view(25,3);
   axis((3/4)*[-.5 4.5 -2 2 -8.5 1])
   title('Odd Mode Shapes');

Modd = moviein(120);
j=0;
for i=1:120
   cla
fprintf('%.0f out of %.0f frames for ODD modes...\n',i,120);
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
   hold on;
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'b');
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
   tdata2a=a(1,1)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'b');
   
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k'); 
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'b'); 
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k'); 
   tdata2b=a(1,2)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'b'); 


   tdata2a=-a(1,2)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
   tdata2a=-a(1,2)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'r');
tdata2a=-a(1,2)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'k');
tdata2a=-a(1,2)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a),'r');

   tdata2b=-a(1,1)*c1(i)*(pi/180)+tdata;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k');
   tdata2b=-a(1,1)*c1(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'r');
   tdata2b=-a(1,1)*c1(i)*(pi/180)+tdata+pi;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'k');
   tdata2b=-a(1,1)*c1(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b), 'r');




   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'r');
   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=a(2,1)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+0,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'r');

   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k'); 
   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'b');
   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k');
   tdata2b=a(2,2)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+1,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'b');
   
   tdata2a=-a(2,2)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=-a(2,2)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'r');
   tdata2a=-a(2,2)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'k');
   tdata2a=-a(2,2)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+2,rdata.*cos(tdata2a),rdata.*sin(tdata2a)-3.5,'r');
   
   
   tdata2b=-a(2,1)*c2(i)*(pi/180)+tdata;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k');
   tdata2b=-a(2,1)*c2(i)*(pi/180)+tdata+pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'b');
   tdata2b=-a(2,1)*c2(i)*(pi/180)+tdata+pi;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'k');
   tdata2b=-a(2,1)*c2(i)*(pi/180)+tdata+3*pi/2;
   patch(0*rdata+3,1*rdata.*cos(tdata2b),1*rdata.*sin(tdata2b)-3.5, 'b');

   plot3([0 1],[0 0],[-3.5 -3.5],'k')
   plot3([0 1],[0 0],[0 0],'k')
   plot3(1+[0 1],[0 0],[-3.5 -3.5],'k')
plot3(1+[0 1],[0 0],[0 0],'k')
plot3(2+[0 1],[0 0],[-3.5 -3.5],'k')
plot3(2+[0 1],[0 0],[0 0],'k')

a=-a;
plot3([0 0 1 2 3 3 0],[0 0 0 0 0 0 0],-2+.07*c1(i)*[0 a(1,1) a(1,2) -a(1,2) -a(1,1) 0  0],'b');
plot3([ 0 1 2 3 ],[ 0 0 0 0],-2+.07*c1(i)*[ a(1,1) a(1,2) -a(1,2) -a(1,1) ],'b.');


plot3([0  0 1 2 3 3 0],[0 0 0 0 0 0 0],-3.5-2+.07*c2(i)*[0 a(2,1) a(2,2) -a(2,2) -a(2,1)  0  0],'r');
plot3([0 1 2  3 ],[ 0 0 0 0],-3.5-2+.07*c2(i)*[a(2,1) a(2,2) -a(2,2) -a(2,1) ],'r.');
a=-a;

   
   axis image
   axis off
   view(25,3);
   axis((3/4)*[-.5 4.5 -2 2 -8.5 1])
   j=j+1;
  title('Odd Mode Shapes');


   Modd(:,j) = getframe;
end