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