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