timour_matlab/mscripts/plotpsb.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 ps
global iterations


%param_file_name =   ...
%   input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
%feval(param_file_name) ;

%iterations

%it = input(' Please enter iteration : ','s') 

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

%psfilename=(['PS.' it ]);
%ps=rdmds(psfilename,'b');

%tfilename=(['T.' it ]);
%t=rdmds(tfilename,'b');

sum=zeros(Nx,Ny);

for k=Nz:-1:1,

sum=sum-dz*0.0002*(t(:,:,k)-20);
% sum=sum-9.8*dz*0.0002*(t(:,:,k)-20);
% sum=sum+999.8*9.8*dz*0.0002*(t(:,:,k)-20);
end
psb=sum+ps;

cmax=max(max(ps));
cmin=min(min(ps));
V=[cmin cmax];
V1=[0 -min(min(sum))];

caxis('manual')

subplot(1,2,1)
        title='bottom pressure';
        imagesc(lat,long,psb');shading flat;axis image;caxis(V);colorbar('horizontal');
        set(gca,'ydir','norm')
%        text(150,140,descriptor);
        text(0,140,title);
        text(0,-40,'timestep');text(100,-40,num2str(it));

subplot(1,2,2)
        title='surface pressure';
        imagesc(lat,long,ps');shading flat;axis image;caxis(V);colorbar('horizontal');
        set(gca,'ydir','norm')
        text(0,140,title);

figure
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 ps
11	global iterations
12
13
14	%param_file_name = ...
15	% input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
16	%feval(param_file_name) ;
17
18	%iterations
19
20	%it = input(' Please enter iteration : ','s')
21
22	%path = this_path
23	%cmdstr=['cd ' path ];
24	%eval(cmdstr);
25	%path=pwd
26
27	%psfilename=(['PS.' it ]);
28	%ps=rdmds(psfilename,'b');
29
30	%tfilename=(['T.' it ]);
31	%t=rdmds(tfilename,'b');
32
33	sum=zeros(Nx,Ny);
34
35	for k=Nz:-1:1,
36
37	sum=sum-dz0.0002(t(:,:,k)-20);
38	% sum=sum-9.8dz0.0002*(t(:,:,k)-20);
39	% sum=sum+999.89.8dz0.0002(t(:,:,k)-20);
40	end
41	psb=sum+ps;
42
43	cmax=max(max(ps));
44	cmin=min(min(ps));
45	V=[cmin cmax];
46	V1=[0 -min(min(sum))];
47
48	caxis('manual')
49
50	subplot(1,2,1)
51	title='bottom pressure';
52	imagesc(lat,long,psb');shading flat;axis image;caxis(V);colorbar('horizontal');
53	set(gca,'ydir','norm')
54	% text(150,140,descriptor);
55	text(0,140,title);
56	text(0,-40,'timestep');text(100,-40,num2str(it));
57
58	subplot(1,2,2)
59	title='surface pressure';
60	imagesc(lat,long,ps');shading flat;axis image;caxis(V);colorbar('horizontal');
61	set(gca,'ydir','norm')
62	text(0,140,title);
63
64	figure
65	caxis('manual')
66
67	title='baroclinic pressure';
68	imagesc(lat,long,-sum');shading flat;axis image;caxis(V1);colorbar('vertical');
69	set(gca,'ydir','norm')
70	% text(80,130,descriptor);
71	text(0,130,title);
72	text(0,-20,'timestep');text(50,-20,num2str(it));
73