Télécharger j3sint.eso

Retour à la liste

Numérotation des lignes :

  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.  
  43.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales