MITgcm_contrib/gag_subduction_scripts/diag_sann.m

function [sann,stot,latind,wmstar] = diag_sann(umean,vmean,wmean,maxmld,delZ,dxc,dyc,raw,mask);
%function [sann,stot,latind,wmstar] = diag_sann(umean,vmean,wmean,maxmld,delZ,dxc,dyc,raw,mask)
%
% Diagnose annual subduction rate of Marshall et al 1993.
% S_ann = -w_H - u_H . del H, [m/yr]
%
% Also, diagnose S_tot, total subduction estimated from
%  annual subduction rate.
% S_tot = \int S_ann dt dA, [Sv]
%
% intermediate terms of calculation:
% latind = u_H . del H = lateral induction
% wmstar = w_H = vertical velocity at h = maxmld.
%
% mask = 2D mask for calculation of subrate.
%
% Started: D. Jamous 1996, Fortran diagnostics.
% Updated: G. Gebbie, 2003, MIT-WHOI for Matlab.

 %% map the mean velocity onto the maxmld surface.
 [umstar,vmstar,wmstar] = get_mldvel(umean,vmean,wmean,delZ,maxmld);

 %% compute mean lateral induction.
 [latind] = diag_induction(umstar,vmstar,maxmld,dxc,dyc);

 sann = -wmstar - latind;

 sann = sann .*86400 .*365; %convert to meters/year.
 
 sanntmp = sann;
 sanntmp(isnan(sanntmp))=0;
 stot=sum(sum(sanntmp.*raw.*mask))./(86400)./365 
 
 return
 
 
 
1	gebbie	1.1	function [sann,stot,latind,wmstar] = diag_sann(umean,vmean,wmean,maxmld,delZ,dxc,dyc,raw,mask);
2			%function [sann,stot,latind,wmstar] = diag_sann(umean,vmean,wmean,maxmld,delZ,dxc,dyc,raw,mask)
3			%
4			% Diagnose annual subduction rate of Marshall et al 1993.
5			% S_ann = -w_H - u_H . del H, [m/yr]
6			%
7			% Also, diagnose S_tot, total subduction estimated from
8			% annual subduction rate.
9			% S_tot = \int S_ann dt dA, [Sv]
10			%
11			% intermediate terms of calculation:
12			% latind = u_H . del H = lateral induction
13			% wmstar = w_H = vertical velocity at h = maxmld.
14			%
15			% mask = 2D mask for calculation of subrate.
16			%
17			% Started: D. Jamous 1996, Fortran diagnostics.
18			% Updated: G. Gebbie, 2003, MIT-WHOI for Matlab.
19
20			%% map the mean velocity onto the maxmld surface.
21			[umstar,vmstar,wmstar] = get_mldvel(umean,vmean,wmean,delZ,maxmld);
22
23			%% compute mean lateral induction.
24			[latind] = diag_induction(umstar,vmstar,maxmld,dxc,dyc);
25
26			sann = -wmstar - latind;
27
28			sann = sann .86400 .365; %convert to meters/year.
29
30			sanntmp = sann;
31			sanntmp(isnan(sanntmp))=0;
32			stot=sum(sum(sanntmp.raw.mask))./(86400)./365
33
34			return
35
36
37