% script plotbas(file1)
% Esta funcion toma un estado b´asico y calcula su diagram de velocidades, 
%
%
load datos.fluidos.mat
x = cos(pi.*(0:n)./n);
y = cos(pi.*(0:m)./m);
% Plot
 
% Dibujamos la temperatura

% Reescalamos

b=[0.01;0.02];
A=[1 1; -1 1];
c=A\b;
xp=polyval(c,x);

b=[0;long];

c=A\b;
yp=polyval(c,y);

subplot(2,1,1)
pcolor(xp,yp,t0);shading interp;colorbar ;axis image
ps=get(gca,'Position');
set(gca,'fontsi',16)
xlabel('x','fonts',16);ylabel('y','fonts',16)

%hold on
%contourf(xp,yp,T',20)
%hold off

%title (' Temperature and velocity field for dif=',num2str(dif),'and c=',num2str(c));
%Dibujo quiver

subplot(2,1,2)
H=quiver(xp,yp,ur0,uz0,2);axis image
axis([min(xp) max(xp), yp(1) yp(end)])
pfin=get(gca,'Position');
ps(2)=pfin(2);
set(gca,'Position',ps)
% pcolor(x,y,w');shading interp,colorbar
set(gca,'fontsi',16)
xlabel('x','fonts',16);ylabel('y','fonts',16)