Sources Esope | Cast3M

coupec
C COUPEC    SOURCE    CHAT      05/01/12    22:24:27     5004
      SUBROUTINE COUPEC(XC,XT,X2,X3,XN,TOL, DIS)
C
C     Esclave de COUPE
C     Distance du point XC au plan defini par les points XT, X2 et X3
C              et la normale XN
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
C
      DIMENSION XC(3),XT(3),X2(3),X3(3),XN(3)
C
      DIS=SQRT((XC(1)-XT(1))**2+(XC(2)-XT(2))**2+(XC(3)-XT(3))**2)
        IF(DIS.GT.TOL)THEN
      DIS=((XC(1)-XT(1))*XN(1)+(XC(2)-XT(2))*XN(2)+(XC(3)-XT(3))*XN(3))
        ELSE
      DIS=SQRT((XC(1)-X2(1))**2+(XC(2)-X2(2))**2+(XC(3)-X2(3))**2)
          IF(DIS.GT.TOL)THEN
      DIS=((XC(1)-X2(1))*XN(1)+(XC(2)-X2(2))*XN(2)+(XC(3)-X2(3))*XN(3))
          ELSE
      DIS=SQRT((XC(1)-X3(1))**2+(XC(2)-X3(2))**2+(XC(3)-X3(3))**2)
      DIS=((XC(1)-X3(1))*XN(1)+(XC(2)-X3(2))*XN(2)+(XC(3)-X3(3))*XN(3))
          ENDIF
        ENDIF
C
      RETURN
      END