1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
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
|
+
+
+
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
-
-
+
+
+
+
|
%% CR file for sFMR
%% Just for flow input definition of r, ie must use with flow,r;k_s
in "$MTT_CRS/r/slin.cr";
%%% Sensitivity RS component CR
in "lin.cr";
in "slin.cr";
FOR ALL e_r,f_r,de_r,df_r,Temp,dTemp LET
sRS(flow,6,
e_r,effort,1,
f_r,flow,2,
de_r,effort,3,
df_r,flow,4,
Temp,effort,5,
dTemp,effort,6
FOR ALL gain_cause,r,rs,out_cause,inp,sinp,in_cause,temp,stemp
LET sRS(gain_cause,r,rs,out_cause,4,
inp,flow,1,
temp,effort,2,
sinp,flow,3,
stemp,effort,4
)
= 2*inp*sinp*r/temp
+ (inp^2)*rs/temp
- (inp^2)*r*temps/(temp^2);
%% - effort in
FOR ALL gain_cause,r,rs,out_cause,inp,sinp,in_cause,temp,stemp
LET sRS(gain_cause,r,rs,out_cause,4,
inp,effort,1,
temp,effort,2,
sinp,effort,3,
stemp,effort,4
) = ((e_r*df_r + de_r*f_r)*Temp - e_r*f_r*dTemp)/(Temp^2);
|