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