clear 
clear path 
whos
global Nx Ny Nz
global lat long dz dm mdep
global delt_su su_its t_su delt
global descriptor this_path
global f deltaf Q beta r_expt r_heat H
global time rots it
global g Cp rho_bar alpha
global u v t w 
global iterations

for iii=1:4,
getuvwt;

%i=find(t<20.1);
%t(i)=NaN;
t(1,:,:)=NaN;
t(Nx,:,:)=NaN;
t(:,:,3)=30;
t(:,:,Nz)=NaN;

        colormap gray
        cmap=flipud(colormap);
        colormap(cmap);


cmin=20;
cmax=30;
V=[cmin cmax];
% figure
caxis('manual')

subplot(4,2,2*iii-1)
%        title='in-situ temperature at z=5';
        imagesc(lat,long,squeeze(t(:,:,5))');shading flat;axis image;caxis(V);colorbar
        set(gca,'ydir','norm','XTickLabel','|||','YTickLabel',' |||')
	xlabel('x','FontSize',15)
	ylabel('y','FontSize',15,'Rotation',[0])

	if iii==1
	text(-30,110,'a)','FontSize',15)
	end

	if iii==2
	text(-30,110,'b)','FontSize',15)
	end

	if iii==3
	text(-30,110,'c)','FontSize',15)
	end

	if iii==4
	text(-30,110,'d)','FontSize',15)
	end

%        text(0,-50,descriptor);
%        text(0,140,title);
%        text(0,-30,'timestep');text(80,-30,num2str(it))
subplot(4,2,2*iii)
        [c,h]=contourf(flipud(squeeze(t(:,Ny/2,3:30))'),7);shading flat;caxis(V);box on;
%	set(h(:),'LineWidth',[0.001])
%        imagesc(lat,long,squeeze(t(:,Ny/2,3:30))');shading flat;axis image;caxis(V);
        set(gca,'DataAspectRatio',[4,1,1],'XTickLabel','|||','YTickLabel',' |||')
	xlabel('x','FontSize',15)
	ylabel('z','FontSize',15,'Rotation',[0])
%        text(0,-2,'in-situ east-west temperature at mid-tank')

end