C                                                                       20170000
C***********************************************************************
C                 COPYRIGHT ATLANTIC RICHFIELD COMPANY 1991            *
C***********************************************************************
CTITLESADM3DM - SOLVE FOR 3 SIMULTANEOUS EQUATIONS                      21330000
CA                                                                      21340000
CA   DESIGNER     JAMES SUN                                             21350000
CA   AUTHOR       JAMES SUN                                             21360000
CA   LANGUAGE     FORTRAN 77                                            21370000
CA   SYSTEM       CRAY/IBM                                              00100000
CA   WRITTEN      12/02/87                                              21390000
C    REVISED      12/17/91   JJC   MODIFIED TO MEET EDP STANDARDS       21400000
C                                                                       21410000
C
C    PURPOSE: TO SOLVE FOR                                              21420000
C                A(1)*X+B(1)*Y+Z=D(1)                                   21430000
C                A(2)*X+B(2)*Y+Z=D(2)                                   21440000
C                A(3)*X+B(3)*Y+Z=D(3)                                   21450000
CA                                                                      21460000
CA CALLING PROCEDURE:
CA     SUBROUTINE SADM3DM(A,B,D,X,Y,Z)
CA
C   CALLING ARGUMENTS
CA
CA   INPUT  A        : COEFF MATRIX ARRAY                            R4 21480000
CA   INPUT  B        : COEFF MATRIX ARRAY                            R4 21490000
CA   INPUT  D        : COEFF MATRIX ARRAY                            R4 21500000
CA   OUTPUT X        : SOLVED X COEFFICIENT                          R4 21510000
CA   OUTPUT Y        : SOLVED Y COEFFICIENT                          R4 21520000
CA   OUTPUT Z        : SOLVED Z COEFFICIENT                          R4 21530000
C                                                                       21540000
      SUBROUTINE SADM3DM(A,B,D,X,Y,Z)                                   21550000
C                                                                       21570000
C
      IMPLICIT INTEGER (A-Z)                                            21560008
C                                                                       21570000
C
      REAL  A
      REAL  A2A1
      REAL  A3A1
      REAL  B
      REAL  B2B1
      REAL  B3B1
      REAL  D
      REAL  D2D1
      REAL  D3D1
      REAL  TEMP1
      REAL  X
      REAL  Y
      REAL  Z
C                                                                       21580000
C
      DIMENSION A(3),B(3),D(3)                                          21590000
C                                                                       21600000
C
      A3A1=A(3)-A(1)                                                    21610000
      A2A1=A(2)-A(1)                                                    21620000
      D3D1=D(3)-D(1)                                                    21630000
      D2D1=D(2)-D(1)                                                    21640000
      B3B1=B(3)-B(1)                                                    21650000
      B2B1=B(2)-B(1)                                                    21660000
C                                                                       21670000
C
      IF(A3A1.EQ.0.  .AND. A2A1.EQ.0.) THEN                             21680000
        X=0.                                                            21690000
        Y=D3D1/B3B1                                                     21700000
        Z=D(1)-B(1)*Y                                                   21710000
      ELSE IF(B3B1.EQ.0.  .AND. B2B1.EQ.0.) THEN                        21720000
        Y=0.                                                            21730000
        X=D3D1/A3A1                                                     21740000
        Z=D(1)-A(1)*X                                                   21750000
      ELSE                                                              21760000
        TEMP1 = A3A1*B2B1 - A2A1 * B3B1                                 21770000
        X=(D3D1*B2B1-D2D1*B3B1)/(TEMP1)                                 21780000
        Y=-1.0 * ((D3D1*A2A1-D2D1*A3A1)/(TEMP1))                        21790000
        Z=D(1)-A(1)*X-B(1)*Y                                            21800000
      END IF                                                            21810000
      RETURN                                                            21820000
      END                                                               21830000