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

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

t0=max(max(t(:,:,3)));
[I,J]=find(t(:,:,3)==t0);
I
J

subplot(1,2,1)
        title='in-situ temperature at z=3';
        imagesc(lat,long,squeeze(t(:,:,3))');shading flat;axis image;caxis(V);colorbar('horizontal');
        set(gca,'ydir','norm')
        text(0,-50,descriptor);
        text(0,140,title);
        text(0,-30,'timestep');text(90,-30,num2str(it))
subplot(1,2,2)
        contour(flipud(squeeze(t(:,J,:))'),50);caxis(V);colorbar;grid
        text(0,-2,'in-situ east-west temperature through the warmest point')
figure
        contour(flipud(squeeze(t(I,:,:))'),50);caxis(V);colorbar;grid
        text(0,-2,'in-situ north-south temperature through the warmest point');text(100,-2, num2str(Ny/2-4))

% figure
%        contour(flipud(squeeze(t(:,Ny/2-8,:))'),50);caxis(V);colorbar;grid
%        text(0,-2,'in-situ east-west temperature at');text(100,-2, num2str(Ny/2-8))
%
%figure
%        contour(flipud(squeeze(t(:,Ny/2-12,:))'),50);caxis(V);colorbar;grid
%        text(0,-2,'in-situ east-west temperature at');text(100,-2, num2str(Ny/2-12))