Isovit.m


close all;
clear all;
load uvw.res
t=uvw(1,1);
ni=uvw(1,2);
nj=uvw(1,3);
nk=uvw(1,4);
for j=1:nj
for i=1:ni
X(i,j)=uvw(1+j+(i-1)*nj,1);
Y(i,j)=uvw(1+j+(i-1)*nj,2);
U(i,j)=uvw(1+j+(i-1)*nj,3)*10.;
V(i,j)=uvw(1+j+(i-1)*nj,4)*10.;
end
end
%
%axis equal tight
% ajustage de la longueur des axes
axes('Position',[0.075 0.3 0.9 0.30]);
hold on quiver(X,Y,U,V,0.5);
hold off


Return