timour_matlab/mscripts/plottransport.m

clear 
clear path 
whos
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

figure
plot(squeeze(v(Nx/2,33,:)),mdep,'o-');  
hold on
plot(squeeze(v(33,Ny/2,:)),mdep,'x:');
xlabel('transport m/s')
ylabel('depth m')
grid
anchor1=v(Nx/2,Ny/4,Nz/2);
text((anchor1*0.9),mdep(Nz/2),['v(33,58,:) - RADIAL transport']);
anchor2=v(Nx/4,Ny/2,Nz/2);
text((anchor2*0.9),mdep(round(3*Nz/4)),['v(58,33,:) - AXIMUTHAL transport'] );

title([descriptor ' timestep ' num2str(eval(it))])
1	clear
2	clear path
3	whos
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	figure
15	plot(squeeze(v(Nx/2,33,:)),mdep,'o-');
16	hold on
17	plot(squeeze(v(33,Ny/2,:)),mdep,'x:');
18	xlabel('transport m/s')
19	ylabel('depth m')
20	grid
21	anchor1=v(Nx/2,Ny/4,Nz/2);
22	text((anchor1*0.9),mdep(Nz/2),['v(33,58,:) - RADIAL transport']);
23	anchor2=v(Nx/4,Ny/2,Nz/2);
24	text((anchor20.9),mdep(round(3Nz/4)),['v(58,33,:) - AXIMUTHAL transport'] );
25
26	title([descriptor ' timestep ' num2str(eval(it))])