Page 456 - Advanced Mine Ventilation

Appendix B: Ventilation Network Analyzer in Fortran IV 425

C
C BEGIN ITERATION
IB¼NFIXBþ1
DO 108 IT¼1, MAXIT
IF(NFIXB)820, 820, 810
810 JE¼ME(NFIX8)
GO TO 830
820 JE¼0
830 L¼0
SUMD¼0.
DO 107 I¼IB, NM
IM¼IeNFIXB
JS¼JEþ1
JE¼ME(I)
SUMH¼ eSUMNVP(I)
SUMDH¼0.
DNF¼0.
IF(NFBPFeI)98, 94, 94
94 N¼ND(IM)
FQ(IM)¼C(IM, Nþ1)*Q(I)
IF(Ne1)97, 97, 95
95 J¼N
DO 96 IJ¼2, N
FQ(IM)¼(FQ(IM)þC(IM, J))*Q(I)
DHF¼(DHFþFLOAT(J)*C(IM, Jþ1))*Q(I)
96 J¼Je1
97 FQ(IM)¼FQ(IM)þC(IM, 1)
SUMH¼SUMHeFQ(IM)
DHF¼DHFþC(IM, 2)
98 DO 101 J¼JS, JE
K¼IABS(NA(J))
DH¼R(K)*DABS(Q(K))
H¼DH*Q(K)
SUMDH¼SUMDHþDH
IF(NA(J))100, 99, 99
99 SUMH¼SUMHþH
GO TO 101
100 SUMH¼ SUMHeH
101 CONTINUE
SUMDH¼SUMDHþSUMDHeDHF
IF(DABS(SUMDH)el.De20)106, 106, 102
102 D¼eSUMH/SUMDH
DO 105 J¼JS, JE