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