Page 449 - Advanced Mine Ventilation

418 Appendix B: Ventilation Network Analyzer in Fortran IV

C
C
C
IMPLICIT REAL*8(AeH, OeZ)
REAL*8 NVP(550), LENGTH(550)
REAL*8 CDB(550), CDJ(400), QDWALL(550), QDJUNC(400)
REAL*8 TLVMXB(550), TLVMXJ(400)
INTEGER BRANCH
INTEGER*4 JCTN(400), NBOUT(550), JOUT(400), NBIN(550), JIN(400),
1MB(550), MJ(400)
C
DIMENSION NA(4200), JC(999), BRANCH(55 0), R(550), Q(550), RR(550)
1 OUT(550), J1(400), J2(400), ME(350), SUMNVP(350), ND(40), DD(80)
2 X(20), FX(20), W(6), Z(6), S(6), ALPHA(6), BETA(6), SIGMA2(6),
C(40, 6), 3P(7, 20), FQ(40), HH(550), TITLE(10), HEIGHT(550),
WIDTH(550)
DIMENSION SUBTIT(10)
C
EQUIVALENCE (NA, JC), (BRANCH, SUMNVP), (BRANCH(351), C), (Q, RR, OUT
1 NVP, X), Q(21), FX), CQ(41), W), (Q(51), Z), (Q(61), S), (Q(71),
SIGMA2)
2 (Q(81), ALPHA), (Q(91), BETA), (Q(101), P)
C
C READ IDENTIFICATION OF THE PROBLEM
10 READ 500, (TITLE(I), I ¼ 1, 10) Rel
C
C READ PARAMETERS
READ 501, NB, NJ, MAXJ, NF, MAXIT, NBWNVP, NFIXB, E Re2
IF(NB)11, 11, 12
11 STOP
12 PRINT 610, (TITLE(I), I¼1, 10)
NM¼NBeNJþl
NFBPF¼NFIXBþNF
PRINT 502, NB, NJ, MAXJ, NM, NF, MAXIT, NBWNVP, NFIXB, E Pel
C
C READ BRANCH CHARACTERISTICS
DO 15 I¼1, NB
READ 503, J1(I), J2(I), R(I) , FF, HEIGHT(I), WIDTH( I), LENGTH( I)
BRANCH(I)¼I
IF(R(I)) 13, 13, 14
13 T¼HEIGHT(I)*WIDTH(I)
R(I) ¼ FF*(HEIGHT(I)þWIDTH(I))*LENGTH(I)/(2.6*T*T*T )