Page 457 - Advanced Mine Ventilation
P. 457

426                            Appendix B: Ventilation Network Analyzer in Fortran IV

            K¼IABS(NA(J))
            IF(NA(J))104, 103, 103
            103 Q(K)¼Q(K)þD
            GO TO 105
            104 Q(K) ¼ Q(K)eD
            105 CONTINUE
            SUMD¼ SUMDþDABS (D)
            IF(DABS(D)eE)107, 107, 106
            106 Lel
            107 CONTINUE
            IF(L)109, 109, 108
            108 CONTINUE
            109 PRINT 519, IT, SUMD Pe4
            PRINT 520
            GO TO 110
            109 PRINT 519, IT, SUMD Pe5
            C
            C COMPUTE HEAD LOSS
            C PRINT SOLUTION
            110 PRINT 521
            DO 111 I¼1, NB
            HH(I)¼R(I)*DABS(Q(I))*Q(I )
            111 PRINT 522, I, J1(I ), J2(I), R(I), Q(I), HH(I) Pe6
            C
            C COMPUTE THE PRESURE TO BE ADJUSTED IN A FIXED QUANTITY BRANCH
            C -- R(I)*DABS(Q(I) )*Q(I ) IS ZERO IF R(I) OF A FIXED QUANTITY BRANCH IS
            C -- NOT GIVEN
            IF(NFIXB)880, 880, 850
            850 PRINT 600
            JE¼0
            DO 870 I¼1, NFIXB
            JS¼JEþ1
            JE¼ME(I)
            SUMH¼eSUMNVP(I)
            DO 860 J¼JS, JE
            K¼IABS(NA(J))
            IF(NA(J)) 856, 855, 855
            855 SUMH¼SUMHþHH(K)
            GO TO 860
            586 SUMH¼SUMHeHH(K)
            860 CONTINUE
            PRINT 522, I, J1(I), J2(I), R(I), Q(I), SUMH
            C CALCULATE THE SIZE OF THE REGULATOR USING MURGUE'S FORMULA
   452   453   454   455   456   457   458   459   460   461   462