timour_matlab/mscripts/plotpb.m

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));

1	edhill	1.1	clear path
2
3			global Nx Ny Nz
4			global lat long dz dm mdep
5			global delt_su su_its t_su delt
6			global descriptor this_path
7			global f deltaf Q beta r_expt r_heat H
8			global time rots it
9			global g Cp rho_bar alpha
10			global u v t w
11			global iterations
12
13
14			path = this_path
15			cmdstr=['cd ' path ];
16			eval(cmdstr);
17			path=pwd
18
19
20			sum=zeros(Nx,Ny);
21
22			for k=Nz:-1:2,
23			sum=sum+9.8dz0.0002*(t(:,:,k)-20);
24			end
25
26			cmax=0.0008;
27			cmin=0.0;
28			V1=[cmin cmax];
29
30			caxis('manual')
31
32			caxis('manual')
33
34			title='baroclinic pressure';
35			imagesc(lat,long,sum');shading flat;axis image;caxis(V1);colorbar('vertical');
36			set(gca,'ydir','norm')
37			text(80,130,descriptor);
38			text(0,130,title);
39			text(0,-20,'timestep');text(50,-20,num2str(it));
40