Page 454 - Advanced Mine Ventilation

Appendix B: Ventilation Network Analyzer in Fortran IV 423

P(1, J)¼0.
64 P(2, J)¼1.
W(1)¼NP
BETA(1)¼0.
IF(NPe7)66, 65, 65
65 N1¼6
GO TO 67
66 N1¼NPe1
67 DO 73 1¼1, N1
K¼Iþ1
Z(I)¼0.
DO 68 J¼1, NP
68 Z(l)¼Z(I)þFX(J)*P(K, J)
S(I)¼Z(I)/W(I)
DSQeDSQ-S(I)*S(I)*W(I)
SIGMA2(I)eDSQ/FLOAT(NPeI)
IF(SIGMA2(I)e1.De6)78, 78, 69
69 IF(IeN1)70, 74, 74
70 ALPHA(I)¼0.
DO 71 J¼1, NP
71 ALPHA(I)¼ALPHA(I) þX (J)*P(K, J)*P(K, J)
ALPHA(I)¼ALPHA(I)/W(I)
W(K)¼0.
DO 72 J¼1, NP
P(Kþ1, J)¼(X(J)eALPHA(I))*P(K, J)eBETA(I)*P(I, J)
72 W(K)¼W(K)þP(Kþ1, J)*P(Kþ1, J)
73 BETA(K)¼W(K)/W(I)
74 SMALL¼SIGMA2(1)
I ¼ 1
DO 77 J¼2, N1
IF(SMALLeSIGMA2(J))77, 77, 76
76 I¼J
SMALL¼SIGMA2(J)
77 CONTINUE
78 PRINT 517, SIGMA2(I)
N1 ¼ I
N¼N1e1
IF (N)780, 780, 79
780 ND(L) ¼ 1
C(L, 1)¼ S(l)
C(L, 2)¼0.
GO TO 840
79 ND(L)¼N