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

% i=find(t<20);
% t(i)=20;

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

%subplot(2,1,1)
        title='in-situ temperature at z=3';
        imagesc(lat,long,squeeze(t(1:Nx-1,1:Ny-1,3))');shading flat;axis off;box on;
%        imagesc(lat,long,squeeze(t(:,:,5))');shading flat;axis image;caxis(V);colorbar('vertical');
        set(gca,'ydir','norm')
        text(0,-50,descriptor);
        text(0,140,title);
        text(0,-30,'timestep');text(80,-30,num2str(it))
%subplot(2,1,2)
figure
        contour(flipud(squeeze(t(1:Nx-2,Ny/4,1:Nz-1))'),50);grid;axis off;
%        contour(flipud(squeeze(t(1:Nx-1,Ny/2,3:Nz))'),50);caxis(V);colorbar;grid
%        text(0,-2,'in-situ east-west temperature ')
%figure
%        contour(flipud(squeeze(t(Ny/4,1:Ny-1,1:Nz))'),50);colorbar;grid
%        text(0,-2,'in-situ north-south temperature')