Page 455 - Advanced Mine Ventilation

424 Appendix B: Ventilation Network Analyzer in Fortran IV

DO 81 I¼1, N1
DO 80 J¼1, I
80 P(I, J)¼0.
81 P(I, Jþl)¼l.
DO 82 J¼1, N
82 P(l, Jþ2)¼eP(l, Jþl)*ALPHA(J)eP(l, J)*BETA(J)
DO 83 I¼2, N
DO 83 J¼ I, N
83 P(I, Jþ2) ¼ P(Ie1, Jþ1)eP(I, Jþ1)*ALPHA(J)eP(I, J)*BETA(J)
DO 84 I¼1, N1
C(L, I)¼0.
DO 84 J¼I, Nl
84 C(L, I)¼C(L, I)þP(I, Jþ1)*S(J)
840 PRINT 518, (C(L, I), I¼l , N1)
GO TO 87
85 ND(L)¼ND(KEY)
N1¼ND(KEY)þ1
DD(JJ)¼DD(KEYþNFIXB)
DO 86 J¼l, Nl
86 C(L, J)¼C(KEY, J)
87 CONTINUE
287 IE ¼ NFBPF
IF(NFIXB)89, 89, 288
288 READ 515, (DD(I), I¼1, NFIXB)
GO TO 89
C
C SET INITIAL VALUES
88 DD(I)¼100.
IE¼1
89 DO 90 I¼1, NB
Q(I) ¼ 0.
90 R(I)¼R(I)*l.De04
JE¼0
DO 93 I¼1, IE
JS¼JEþ1
JE¼ME(I)
DO 93 J¼JS, JE
K¼1ABS(NA(J))
IF(NA(J))92, 91, 9l
91 Q(K) ¼ Q(K)þDD(I)
GO TO 93
92 Q(K)¼Q(K)eDD(I)
93 CONTINUE