Isotau.m


close all;
clear all;
load res.tau;
t=res(1,1);
ni=res(1,2);
nj=res(1,3);
for j=1:nj
for i=1:ni
X(i,j)=res(1+j+(i-1)*nj,1);
Y(i,j)=res(1+j+(i-1)*nj,2);
TAU(i,j)=res(1+j+(i-1)*nj,3);
end
end
% echelle conservee
% axis equal tight
% ajustage de la longueur des axes
axes('Position',[0.075 0.3 0.9 0.30]);
v=[0.2 0.5 0.8]
%[C,h]=contour(Y,X,TAU,1,'r');
[C,h]=contour(X,Y,TAU,v);
%set(gca,'XLim',[0.0 3.],'YLim',[0.0 0.5])
set(gca,'XLim',[0.0 3.],'YLim',[0.0 0.5])
%text(2.5,2.1,['t = ',num2str(t)],...
%'HorizontalAlignment','Center')
text(5.,0.12,['t = ',num2str(t)],...
'HorizontalAlignment','Center')

Return