-
Notifications
You must be signed in to change notification settings - Fork 1
/
PrintQTIntBCVals.m
46 lines (34 loc) · 1.12 KB
/
PrintQTIntBCVals.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
function PrintQTIntBCVals(IntBC)
QTGlobals
fprintf('Internal BCs:\n');
NetFlow = 0;
for i=1:length(IntBC)
fprintf('\tRegion %i ',i);
TotalFlow = 0;
AveVal = 0;
nb = length(IntBC(i).blocks);
for j = 1:nb
b = IntBC(i).blocks(j);
AveVal = AveVal + Blocks{b}.data(VVAL);
if IntBC(i).mode == 'v'
for i = 1:MAX_NUM_BR
if Blocks{b}.br(br)
TotalFlow = TotalFlow + Blocks{b}.data(VVAL)-...
Blocks{Blocks{b}.br(br)}.data(VVAL);
end
end
elseif IntBC(i).mode == 'g'
TotalFlow = TotalFlow + ...
(Blocks{b}.data(VVAL)-IntBC(i).val(1))*IntBC(i).val(2);
elseif IntBC(i).mode == 'm'
TotalFlow = TotalFlow + ...
(Blocks{b}.data(VVAL)-IntBC(i).val(1))*IntBC(i).val(2) + ...
IntBC(i).val(3);
end
end
NetFlow = NetFlow + TotalFlow;
AveVal = AveVal/nb;
fprintf(' Average Val = %f TotalFlow = %f\n',AveVal,TotalFlow);
end
fprintf('\n\tNetFlow = %f\n',NetFlow);
end