timour_matlab/mscripts/multi.m

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

sizeit=size(iterations);
  for i=1:round(sizeit(1)/16)
figurenumber=i;
figurenumber
      figure(i)
title(descriptor)
for j=1:16;
figurenumber
subplotnumber=j;
subplotnumber
((i-1)*16)+j
   it = iterations(((figurenumber-1)*16)+j,1:10);
time   = (it-su_its)*delt ;
rots = time*f/(4*pi);
   path   = this_path 
   cmdstr=['cd ' path ];
   eval(cmdstr);
   path=pwd
   % T (oC)
   tfilename=(['T.' it ]);
   t=rdmeta(tfilename,'b');

       subplot(4,4,j)
      polyplot
      drawnow
set(gca,'XTick',[])           
set(gca,'YTick',[])
title(it)

   end
end
1	clear
2	clear path
3
4	global Nx Ny Nz
5	global lat long dz dm mdep
6	global delt_su su_its t_su delt
7	global descriptor this_path
8	global f deltaf Q beta r_expt r_heat H
9	global time rots it
10	global g Cp rho_bar alpha
11	global u v t w
12	global iterations
13
14
15	param_file_name = ...
16	input(' Please enter the name of the m-file with the parameters for this run : ','s') ;
17	feval(param_file_name)
18
19	sizeit=size(iterations);
20	for i=1:round(sizeit(1)/16)
21	figurenumber=i;
22	figurenumber
23	figure(i)
24	title(descriptor)
25	for j=1:16;
26	figurenumber
27	subplotnumber=j;
28	subplotnumber
29	((i-1)*16)+j
30	it = iterations(((figurenumber-1)*16)+j,1:10);
31	time = (it-su_its)*delt ;
32	rots = timef/(4pi);
33	path = this_path
34	cmdstr=['cd ' path ];
35	eval(cmdstr);
36	path=pwd
37	% T (oC)
38	tfilename=(['T.' it ]);
39	t=rdmeta(tfilename,'b');
40
41	subplot(4,4,j)
42	polyplot
43	drawnow
44	set(gca,'XTick',[])
45	set(gca,'YTick',[])
46	title(it)
47
48	end
49	end