Trailing-Edge
-
PDP-10 Archives
-
decuslib20-04
-
decus/20-0101/model5.for
There is 1 other file named model5.for in the archive. Click here to see a list.
SUBROUTINE MODEL5 (NPROB, B, FY, NOB, NC, X, NVARX, NOBMAX, NCMAX,08M50010
1 KTOU) 08M10011
C OXIDES MEASURED IN MOL PERCENTS. D' = D1 TO D5 08M50020
C --- 5 OXIDES( - MIN VALUES) VS D AND A. B, G AND H ARE COMMON. 08M50030
C 08M50040
IMPLICIT REAL*8(A-H,O-Z) 08M50050
DIMENSION B(NCMAX), FY(NOBMAX), X(NVARX,NOBMAX) 08M50060
RB = 1.0/B(1) 08M50070
DO 10 J=1,NOB 08M50080
X1 = (X(2,J) - 20.9) / 7.43 08M50090
X2 = (X(3,J) - 2.03) / 7.43 08M50100
X3 = (X(4,J) - 0.44) / 7.43 08M50110
X4 = (X(5,J) - 66.1) / 7.43 08M50120
X5 = (X(6,J) - 3.10) / 7.43 08M50130
X12 = X1*X2 08M50140
X13 = X1*X3 08M50150
X14 = X1*X4 08M50160
X24 = X2*X4 08M50170
X23 = X2*X3 08M50180
X34 = X3*X4 08M50190
A = B( 4)*X1 - B( 5)*X2 + B( 6)*X3 + B( 7)*X4 + B( 8)*X5 08M50200
1 + B( 9)*X12 - B(10)*X13 - B(11)*X14 - B(12)*X23 + B(13)*X24 08M50210
2 - B(14)*X34 08M50220
D = B(15)*X1 + B(16)*X2 + B(17)*X3 + B(18)*X4 + B(19)*X5 08M50230
C = B( 2)*D - B( 3) 08M50240
F = (A/C)**RB 08M50250
FY(J) = D - C + ( A / ( (X(1,J) + F )**B(1) ) ) 08M50260
10 CONTINUE 08M50270
RETURN 08M50280
END 08M50290