timour_matlab/mscripts/meanps1.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


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


iterations

itstart = input(' Please enter start iteration : ','s') 
itend = input(' Please enter end iteration : ','s') 


sizeit=size(iterations);
for i=1:sizeit(1)
iter(i)=eval(iterations(i,1:10));
end
nitstart=find(iter==eval(itstart))
nitend=find(iter==eval(itend))

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

sum=zeros(Nx,Ny);
counter=0;

for i=nitstart:nitend
psfilename=(['PS.' iterations((i),1:10) ]) ;
ps=rdmds(psfilename,'b');
sum=sum+ps;
counter=counter+1;
end
meanps=sum/counter;

cmax=max(max(meanps));
cmin=min(min(meanps));
V=[cmin cmax];
figure
caxis('manual')

        title='mean surface pressure';
        imagesc(lat,long,meanps');shading flat;axis image;caxis(V);colorbar('vertical');
        set(gca,'ydir','norm')
        text(50,130,descriptor);text(50,-30,num2str(nitend-nitstart+1));text(80,-30,'files');
        text(0,130,title);
        text(0,-10,'from timestep');text(30,-10,num2str(itstart));
        text(70,-10,'to'); text(80,-10,num2str(itend));


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	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
19	iterations
20
21	itstart = input(' Please enter start iteration : ','s')
22	itend = input(' Please enter end iteration : ','s')
23
24
25	sizeit=size(iterations);
26	for i=1:sizeit(1)
27	iter(i)=eval(iterations(i,1:10));
28	end
29	nitstart=find(iter==eval(itstart))
30	nitend=find(iter==eval(itend))
31
32	path = this_path
33	cmdstr=['cd ' path ];
34	eval(cmdstr);
35	path=pwd
36
37	sum=zeros(Nx,Ny);
38	counter=0;
39
40	for i=nitstart:nitend
41	psfilename=(['PS.' iterations((i),1:10) ]) ;
42	ps=rdmds(psfilename,'b');
43	sum=sum+ps;
44	counter=counter+1;
45	end
46	meanps=sum/counter;
47
48	cmax=max(max(meanps));
49	cmin=min(min(meanps));
50	V=[cmin cmax];
51	figure
52	caxis('manual')
53
54	title='mean surface pressure';
55	imagesc(lat,long,meanps');shading flat;axis image;caxis(V);colorbar('vertical');
56	set(gca,'ydir','norm')
57	text(50,130,descriptor);text(50,-30,num2str(nitend-nitstart+1));text(80,-30,'files');
58	text(0,130,title);
59	text(0,-10,'from timestep');text(30,-10,num2str(itstart));
60	text(70,-10,'to'); text(80,-10,num2str(itend));
61
62