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=0.75;
cmax=1.75,
V=[cmin cmax];
figure
caxis('manual')

subplot(1,2,1)
        title='in-situ time mean qgpv at z=2';
        imagesc(lat,long,squeeze(meanq(:,:,2))');shading flat;axis image;caxis(V);colorbar('horizontal');
        text(0,-30,descriptor);
        text(0,140,title);
        text(0,-20,'timestep');text(40,-20,int2str(eval(itstart)));
        text(40,-12,'to');text(40,-4,int2str(eval(itend)));
subplot(1,2,2)
        contour(flipud(squeeze(meanq(:,Ny/2,:))'),100);caxis(V);colorbar;grid
        text(0,-2,'in-situ cross-tank time mean qgpv')