forked from csherwood-usgs/fm_flocmod
-
Notifications
You must be signed in to change notification settings - Fork 0
/
fm_comp_fsd.m
67 lines (51 loc) · 1.43 KB
/
fm_comp_fsd.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
% fm_comp_fsd
% This processes NNin and returns NNout
NNout = zeros(size(NNin));
tmp_g1=0.0;
tmp_g3=0.0;
tmp_g4=0.0;
tmp_l1=0.0;
tmp_l3=0.0;
tmp_l4=0.0;
f_g1_tmp=zeros(nv_mud,nv_mud,nv_mud);
f_l1_tmp=zeros(nv_mud,nv_mud);
if (l_COLLFRAG)
fm_collfrag
end
for iv1=1:nv_mud
% for iv2=1:nv_mud
% for iv3=1:nv_mud
%if (l_ASH)
f_g1_tmp(:,:,iv1)=f_g1_tmp(:,:,iv1)+l_ASH*f_g1_sh(:,:,iv1)*Gval;
%end
%if (l_ADS)
f_g1_tmp(:,:,iv1)=f_g1_tmp(:,:,iv1)+l_ADS*f_g1_ds(:,:,iv1);
%end
tmp_g1=tmp_g1+(NNin'*(f_g1_tmp(:,:,iv1))*NNin);
%if (l_COLLFRAG)
tmp_g4=tmp_g4+l_COLLFRAG*(NNin'*(f_g4(:,:,iv1)*Gval)*NNin);
%end
% end
tmp_g3=tmp_g3+f_g3(:,iv1)'*NNin*Gval^1.5;
%if (l_ASH)
f_l1_tmp(:,iv1)=f_l1_tmp(:,iv1)+l_ASH*f_l1_sh(:,iv1)*Gval;
%end
%if (l_ADS)
f_l1_tmp(:,iv1)=f_l1_tmp(:,iv1)+l_ADS*f_l1_ds(:,iv1)*Gval;
%end
tmp_l1=tmp_l1+(f_l1_tmp(:,iv1))'*NNin;
%if (l_COLLFRAG)
tmp_l4=tmp_l4+l_COLLFRAG*(f_l4(:,iv1)*Gval)'*NNin;
%end
% end
tmp_l1=tmp_l1*NNin(iv1);
tmp_l4=tmp_l4*NNin(iv1);
tmp_l3=f_l3(iv1)*Gval^1.50*NNin(iv1);
NNout(iv1)=NNin(iv1)+f_dt*(tmp_g1+tmp_g3+tmp_g4-(tmp_l1+tmp_l3+tmp_l4));
tmp_g1=0.0;
tmp_g3=0.0;
tmp_g4=0.0;
tmp_l1=0.0;
tmp_l3=0.0;
tmp_l4=0.0;
end