MITgcm_contrib/gag_subduction_scripts/get_mldvel.m

function [ustar,vstar,wstar] = get_mldvel(u,v,w,delZ,h)
%function [ustar,vstar,wstar] = get_mldvel(u,v,w,delZ,h)
%
% Get velocity at a surface h = h(x,y).
% Velocity remains on the C-grid with depths "depth".
%
% Started: D. Jamous, 1996, FORTRAN diags.
%
% Translated: G. Gebbie, MIT-WHOI, November 2003.

 [nx,ny,nz]=size(u);

 ustar = zeros(nx,ny);
 vstar = zeros(nx,ny);
 wstar = zeros(nx,ny);

 zbot = cumsum(delZ);

 depth = get_depths(delZ);

 for i=2:nx-1
   for j=2:ny-1
     ustar(i,j) = interp1( depth, sq(u(i,j,:)),(h(i,j)+h(i-1,j))./2,'linear');
     vstar(i,j) = interp1( depth, sq(v(i,j,:)),(h(i,j)+h(i,j-1))./2,'linear');
   end
 end
 for i=1:nx-1
   for j=1:ny-1
     wstar(i,j) = interp1( sq(zbot(1:23)), sq(w(i,j,:)), h(i,j), 'linear');
   end
 end

 ustar(isnan(ustar))= 0;
 vstar(isnan(vstar))= 0;
 wstar(isnan(wstar))= 0;
1	function [ustar,vstar,wstar] = get_mldvel(u,v,w,delZ,h)
2	%function [ustar,vstar,wstar] = get_mldvel(u,v,w,delZ,h)
3	%
4	% Get velocity at a surface h = h(x,y).
5	% Velocity remains on the C-grid with depths "depth".
6	%
7	% Started: D. Jamous, 1996, FORTRAN diags.
8	%
9	% Translated: G. Gebbie, MIT-WHOI, November 2003.
10
11	[nx,ny,nz]=size(u);
12
13	ustar = zeros(nx,ny);
14	vstar = zeros(nx,ny);
15	wstar = zeros(nx,ny);
16
17	zbot = cumsum(delZ);
18
19	depth = get_depths(delZ);
20
21	for i=2:nx-1
22	for j=2:ny-1
23	ustar(i,j) = interp1( depth, sq(u(i,j,:)),(h(i,j)+h(i-1,j))./2,'linear');
24	vstar(i,j) = interp1( depth, sq(v(i,j,:)),(h(i,j)+h(i,j-1))./2,'linear');
25	end
26	end
27	for i=1:nx-1
28	for j=1:ny-1
29	wstar(i,j) = interp1( sq(zbot(1:23)), sq(w(i,j,:)), h(i,j), 'linear');
30	end
31	end
32
33	ustar(isnan(ustar))= 0;
34	vstar(isnan(vstar))= 0;
35	wstar(isnan(wstar))= 0;