clear path 

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


path   = this_path 
cmdstr=['cd ' path ];
eval(cmdstr);
path=pwd


sum=zeros(Nx,Ny);

for k=Nz:-1:2,
sum=sum+9.8*dz*0.0002*(t(:,:,k)-20);
end

cmax=0.0008;
cmin=0.0;
V1=[cmin cmax];

caxis('manual')

caxis('manual')

        title='baroclinic pressure';
        imagesc(lat,long,sum');shading flat;axis image;caxis(V1);colorbar('vertical');
        set(gca,'ydir','norm')
        text(80,130,descriptor);
        text(0,130,title);
        text(0,-20,'timestep');text(50,-20,num2str(it));