Esope source files

   1 : C J3SINT    SOURCE    CHAT      05/01/13    00:47:23     5004           
   2 :       SUBROUTINE J3SINT(XY1,XY2,XY3,XY4,XTOL,XLAM,XGAM,XY)
   3 : C----------------------------------------------------
   4 : C     SEGMENT INTERSECTION
   5 : C
   6 : C     ESCLAVE DE J3COTO
   7 : C
   8 : C     PP 6/97
   9 : C     Pierre Pegon/JRC Ispra
  10 : C----------------------------------------------------
  11 : C
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8(A-H,O-Z)
  14 : -INC CCREEL
  15 :       DIMENSION XY1(2),XY2(2),XY3(2),XY4(2),XY(2)
  16 : C
  17 :       XLAM=XGRAND
  18 :       XGAM=XGRAND
  19 : C
  20 :       DET=(XY2(1)-XY1(1))*(XY3(2)-XY4(2))
  21 :      >   -(XY2(2)-XY1(2))*(XY3(1)-XY4(1))
  22 : C
  23 : C     LES SEGMENTS SONT PARALLELS
  24 : C
  25 :       IF(ABS(DET).LT.XTOL)RETURN
  26 : C
  27 : C     LES SEGMENTS SE COUPENT
  28 : C
  29 :       XLAM=((XY3(1)-XY1(1))*(XY3(2)-XY4(2))
  30 :      >     -(XY3(2)-XY1(2))*(XY3(1)-XY4(1)))/DET
  31 :       XGAM=((XY2(1)-XY1(1))*(XY3(2)-XY1(2))
  32 :      >     -(XY2(2)-XY1(2))*(XY3(1)-XY1(1)))/DET
  33 : C
  34 :       DO IE1=1,2
  35 :         XY(IE1)=XLAM*XY2(IE1)+(1-XLAM)*XY1(IE1)
  36 :       ENDDO
  37 : C
  38 :       RETURN
  39 :       END
  40 :  
  41 :  
  42 :