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