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

figure
subplot(2,1,1)
plot(squeeze(w(Nx/2,Ny/2,1:Nz)),'o-');
%set(gca,'YLim',[mdep(length(mdep)) 0]);
xlabel('depth , gridpoints')
ylabel('w m/s')
grid

figure

cmin=-1.e-4;
cmax=1.e-4,
V=[cmin cmax];

title='in-situ vertical velocity at the 4rd level from the top';
        imagesc(lat,long,squeeze(w(:,:,4))');shading flat;axis image;caxis(V);colorbar('vertical');
        set(gca,'ydir','norm')
        text(0,-30,descriptor);
        text(0,130,title);
        text(0,-20,'timestep');text(50,-20,num2str(it))


%anchor1=w(Nx/2,Ny/2,Nz/2);
%text((anchor1*0.9),mdep(Nz/2),num2str(anchor1));
%anchor2=eval(it);
%text((anchor1*0.9),mdep(round(3*Nz/4)),['timestep  ' num2str(anchor2)] );

%subplot(2,1,2)
%for i=1:length(mdep)
%hold on;plot(squeeze(w(:,Nx/2,i)),':');
%end
%plot(squeeze(w(:,Nx/2,2)),'--');
%plot(squeeze(w(:,Nx/2,Nz/2)),'-');
%set(gca,'XLim',[0 Nx]);
%grid on
%ylabel('w m/s')
%xlabel('distance across the domain (gridpoints)')
%
%subplot(2,1,1)
%title(descriptor)