Google
 

Trailing-Edge - PDP-10 Archives - decuslib10-06 - 43,50402/model4.pgm
There are 2 other files named model4.pgm in the archive. Click here to see a list.
      SUBROUTINE MODEL4 (NPROB, B, FY, NOB, NC, X, NVARX, NOBMAX, NCMAX,07M40010
     1 KTOU)                                                            07M10011
C                                                                       07M40020
C --- 5 OXIDES( - MIN VALUES) VS D AND A.  B, G AND H ARE COMMON.       07M40030
C                                                                       07M40040
      IMPLICIT REAL*8(A-H,O-Z)                                          07M40050
      DIMENSION  B(NCMAX), FY(NOBMAX), X(NVARX,NOBMAX)                  07M40060
C                                                                       07M40070
      RB  = 1.0/B(1)                                                    07M40080
      DO 10  J=1,NOB                                                    07M40090
      X1 = (X(2,J) - 21.0) / 10.66                                      07M40100
      X2 = (X(3,J) - 3.48) / 10.66                                      07M40110
      X3 = (X(4,J) - 1.18) / 10.66                                      07M40120
      X4 = (X(5,J) - 61.6) / 10.66                                      07M40130
      X5 = ( 97.92 - (X(2,J) + X(3,J) + X(4,J) + X(5,J))) / 10.66       07M40140
      X6 = X1*X3                                                        07M40150
      X7 = X1*X4                                                        07M40160
      X8 = X2*X4                                                        07M40170
      X9 = X3*X4                                                        07M40180
      X10 = X1*X2                                                       07M40190
      D   = B( 2)*X1 + B( 3)*X2 + B( 4)*X3 + B( 5)*X4 + B( 6)*X5        07M40200
     1               - B( 7)*X6 + B( 8)*X8                              07M40210
      A   = B( 9)*X1 + B(10)*X2 + B(11)*X3 + B(12)*X4 + B(13)*X5        07M40220
     1               + B(14)*X10- B(15)*X7                              07M40230
      C = B(16)*D - B(17)                                               07M40240
      F = (A/C)**RB                                                     07M40250
      FY(J) = D - C + ( A / ( (X(1,J) + F )**B(1) )    )                07M40260
   10 CONTINUE                                                          07M40270
      RETURN                                                            07M40280
      END                                                               07M40290