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