MITgcm_contrib/gag_subduction_scripts/diag_induction.m

function [induction,gradx,grady] = diag_induction(ustar,vstar,h,dxc,dyc);
%function [induction,gradx,grady] = diag_induction(ustar,vstar,h,dxc,dyc)
%
% Diagnose lateral induction u_h . grad h 
%
% G. Gebbie, 2003. 

 [nx,ny] = size(ustar);

 gradx(2:nx,:) = (h(2:nx,:)  - h(1:nx-1,:));
 grady(:,2:ny) =  h(:,2:ny)  - h(:,1:ny-1);

 gradx = gradx ./ dxc;
 grady = grady ./ dyc;

 udelh = ustar .* gradx;
 vdelh = vstar .* grady;

%% now move udelh from U points to H points, in order to match up with W*.
%% involves an average.
 udelh2 = (udelh(2:nx,:)+udelh(1:nx-1,:))./2;
 vdelh2 = (vdelh(:,2:ny)+vdelh(:,1:ny-1))./2;

 udelh2(nx,:) = 0;
 vdelh2(:,ny)=0;

induction = udelh2 + vdelh2;

1	function [induction,gradx,grady] = diag_induction(ustar,vstar,h,dxc,dyc);
2	%function [induction,gradx,grady] = diag_induction(ustar,vstar,h,dxc,dyc)
3	%
4	% Diagnose lateral induction u_h . grad h
5	%
6	% G. Gebbie, 2003.
7
8	[nx,ny] = size(ustar);
9
10	gradx(2:nx,:) = (h(2:nx,:) - h(1:nx-1,:));
11	grady(:,2:ny) = h(:,2:ny) - h(:,1:ny-1);
12
13	gradx = gradx ./ dxc;
14	grady = grady ./ dyc;
15
16	udelh = ustar .* gradx;
17	vdelh = vstar .* grady;
18
19	%% now move udelh from U points to H points, in order to match up with W*.
20	%% involves an average.
21	udelh2 = (udelh(2:nx,:)+udelh(1:nx-1,:))./2;
22	vdelh2 = (vdelh(:,2:ny)+vdelh(:,1:ny-1))./2;
23
24	udelh2(nx,:) = 0;
25	vdelh2(:,ny)=0;
26
27	induction = udelh2 + vdelh2;
28