clear
clf
file=input('Enter output file (default = output)','s');
if isempty(file)
   file='output';
end   
load(strcat(file,'\x.out'))
load(strcat(file,'\y.out'))
load(strcat(file,'\umov01'))
load(strcat(file,'\umov02'))
load(strcat(file,'\umov03'))
load(strcat(file,'\umov04'))
load(strcat(file,'\umov05'))
load(strcat(file,'\umov06'))
load(strcat(file,'\umov07'))
load(strcat(file,'\umov08'))
load(strcat(file,'\umov09'))
load(strcat(file,'\umov10'))
load(strcat(file,'\umov11'))
load(strcat(file,'\vmov01'))
load(strcat(file,'\vmov02'))
load(strcat(file,'\vmov03'))
load(strcat(file,'\vmov04'))
load(strcat(file,'\vmov05'))
load(strcat(file,'\vmov06'))
load(strcat(file,'\vmov07'))
load(strcat(file,'\vmov08'))
load(strcat(file,'\vmov09'))
load(strcat(file,'\vmov10'))
load(strcat(file,'\vmov11'))
load(strcat(file,'\tmov01'))
load(strcat(file,'\tmov02'))
load(strcat(file,'\tmov03'))
load(strcat(file,'\tmov04'))
load(strcat(file,'\tmov05'))
load(strcat(file,'\tmov06'))
load(strcat(file,'\tmov07'))
load(strcat(file,'\tmov08'))
load(strcat(file,'\tmov09'))
load(strcat(file,'\tmov10'))
load(strcat(file,'\tmov11'))
h=pcolor(x,y,tmov01);
shading interp
hold on
quiver(x,y,umov01,vmov01)
hold off
M(:,1)=getframe;
h=pcolor(x,y,tmov02);
shading interp
hold on
quiver(x,y,umov02,vmov02)
hold off
M(:,2)=getframe;
h=pcolor(x,y,tmov03);
shading interp
hold on
quiver(x,y,umov03,vmov03)
hold off
M(:,3)=getframe;
h=pcolor(x,y,tmov04);
shading interp
hold on
quiver(x,y,umov04,vmov04)
hold off
M(:,4)=getframe;
h=pcolor(x,y,tmov05);
shading interp
hold on
quiver(x,y,umov05,vmov05)
hold off
M(:,5)=getframe;
h=pcolor(x,y,tmov06);
shading interp
hold on
quiver(x,y,umov06,vmov06)
hold off
M(:,6)=getframe;
h=pcolor(x,y,tmov07);
shading interp
hold on
quiver(x,y,umov07,vmov07)
hold off
M(:,7)=getframe;
h=pcolor(x,y,tmov08);
shading interp
hold on
quiver(x,y,umov08,vmov08)
hold off
M(:,8)=getframe;
h=pcolor(x,y,tmov09);
shading interp
hold on
quiver(x,y,umov09,vmov09)
hold off
M(:,9)=getframe;
h=pcolor(x,y,tmov10);
shading interp
hold on
quiver(x,y,umov10,vmov10)
hold off
M(:,10)=getframe;
h=pcolor(x,y,tmov11);
shading interp
hold on
quiver(x,y,umov11,vmov11)
hold off
M(:,11)=getframe;
movie(M,11,3)
set(gca,'fontweight','bold')
xlabel('X (km)','fontweight','bold')
ylabel('Y (km)','fontweight','bold')
