C SURF3P    SOURCE    CHAT      05/01/13    03:27:49     5004           
      SUBROUTINE SURF3P(P1,P2,P3,S)
C=============================================================
C                                                            =
C     INPUT:                                                 =
C     P1,P2 ET P3 SONT TROIS POINTS DANS UN PLAN             =
C                                                            =
C     OUTPUT:                                                =
C      S : SURFACE DU TRIANGLE                               =
C                                                            =
C     G. L. JIANG                                            =
C                                                            =
C=============================================================
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION P1(*),P2(*),P3(*)
      P=P2(1)-P1(1)
      Q=P2(2)-P1(2)
      R=P2(3)-P1(3)
      A=P3(1)-P1(1)
      B=P3(2)-P1(2)
      C=P3(3)-P1(3)
      C1=C*Q-B*R
      C2=A*R-C*P
      C3=B*P-A*Q
      S=0.5D0*SQRT(C1*C1+C2*C2+C3*C3)
      END