timour_matlab/mscripts/working.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) ;
path   = this_path 
cmdstr=['cd ' path ];
eval(cmdstr);
path=pwd

% create circular sampling rings


Nn=Nx/2;

rad=[1:1:(Nn-1)];
ring=ones(Nx,Ny,Nz);
weight=zeros(1,Nn-1);
radialmean=zeros(1,Nn-1);
XCenter = (Nx+1)/3 ;
YCenter = (Ny+1)/2 ;


for k=1:Nn-1
  for i=1:Nx
    for j=1:Ny
      radius(i,j) = sqrt((i-XCenter)^2+(j-YCenter)^2);
    end
  end
  radius(find( round(radius) > (k+1) ))=radius(find (round(radius) > (k+1) ))*0;
  radius(find( round(radius) < (k)   ))=radius(find (round(radius) < (k)   ))*0;
  radius(find(radius>0))=1.0;
  ring(:,:,k)=radius;
  weight(:,k)=length(find(ring(:,:,k) > 0));
end

sizeit=size(iterations);
theta20=zeros(sizeit(1),1);
figure;hold on
for i=3:sizeit(1)-1
tfilename=(['T.' iterations((i),1:10) ]);
t=rdmds(tfilename,'b');

tring=t(:,:,2).*ring(:,:,41);
meantee=sum(sum(tring))/weight(:,41);

tval(i,:)=meantee-20;
centralt=tval;
plot(eval(iterations((i),1:10)),centralt(i,:),'.');drawnow
tincol(i)=sum((centralt(i,:)).*dz);

%time(i)=(eval(iterations(i,1:10))-su_its)*delt;
time(i)=eval(iterations(i,1:10));


end


tcol=tval;
title(descriptor)
xlabel('timestep')
ylabel('central temperature')
grid

computertime=time;

%save centralt centralt
%save iterations iterations
%save tcol tcol
%save tincol tincol
%save mdep mdep 
save working tcol tincol mdep computertime
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	path = this_path
19	cmdstr=['cd ' path ];
20	eval(cmdstr);
21	path=pwd
22
23	% create circular sampling rings
24
25
26	Nn=Nx/2;
27
28	rad=[1:1:(Nn-1)];
29	ring=ones(Nx,Ny,Nz);
30	weight=zeros(1,Nn-1);
31	radialmean=zeros(1,Nn-1);
32	XCenter = (Nx+1)/3 ;
33	YCenter = (Ny+1)/2 ;
34
35
36	for k=1:Nn-1
37	for i=1:Nx
38	for j=1:Ny
39	radius(i,j) = sqrt((i-XCenter)^2+(j-YCenter)^2);
40	end
41	end
42	radius(find( round(radius) > (k+1) ))=radius(find (round(radius) > (k+1) ))*0;
43	radius(find( round(radius) < (k) ))=radius(find (round(radius) < (k) ))*0;
44	radius(find(radius>0))=1.0;
45	ring(:,:,k)=radius;
46	weight(:,k)=length(find(ring(:,:,k) > 0));
47	end
48
49	sizeit=size(iterations);
50	theta20=zeros(sizeit(1),1);
51	figure;hold on
52	for i=3:sizeit(1)-1
53	tfilename=(['T.' iterations((i),1:10) ]);
54	t=rdmds(tfilename,'b');
55
56	tring=t(:,:,2).*ring(:,:,41);
57	meantee=sum(sum(tring))/weight(:,41);
58
59	tval(i,:)=meantee-20;
60	centralt=tval;
61	plot(eval(iterations((i),1:10)),centralt(i,:),'.');drawnow
62	tincol(i)=sum((centralt(i,:)).*dz);
63
64	%time(i)=(eval(iterations(i,1:10))-su_its)*delt;
65	time(i)=eval(iterations(i,1:10));
66
67
68	end
69
70
71
72	tcol=tval;
73	title(descriptor)
74	xlabel('timestep')
75	ylabel('central temperature')
76	grid
77
78	computertime=time;
79
80	%save centralt centralt
81	%save iterations iterations
82	%save tcol tcol
83	%save tincol tincol
84	%save mdep mdep
85	save working tcol tincol mdep computertime