Télécharger j3poin.eso

Retour à la liste

Numérotation des lignes :

  1. C J3POIN SOURCE CHAT 05/01/13 00:47:17 5004
  2. SUBROUTINE J3POIN(WORK,NPTO,INOEUD,XY,DENS1)
  3. C----------------------------------------------------
  4. C INSERSION D'UN POINT
  5. C
  6. C ESCLAVE DE J3COTO
  7. C
  8. C CODE IST(1,I): 0 point non traite
  9. C 1 est sur le segment IST(2,I)
  10. C 2 est sur les segments IST(2,I) et IST(3,I)
  11. C -1 est a l'interieur
  12. C -2 est a l'exterieur
  13. C
  14. C PP 6/97
  15. C Pierre Pegon/JRC Ispra
  16. C----------------------------------------------------
  17. IMPLICIT INTEGER(I-N)
  18. IMPLICIT REAL*8(A-H,O-Z)
  19. DIMENSION XY(2)
  20. C
  21. SEGMENT WORK
  22. REAL*8 XYC(2,NPTO)
  23. INTEGER IST(3,NPTO)
  24. REAL*8 DENS(NPTO)
  25. INTEGER JUN
  26. ENDSEGMENT
  27. C
  28. SEGADJ,WORK
  29. C
  30. C ON SHIFTE LE CONTENU DE WORK DE INOEUD->NPTO-1 A INOEUD+1->NPTO
  31. C
  32. IF(NPTO.GT.INOEUD)THEN
  33. DO IE1=NPTO,INOEUD+1,-1
  34. DO IE2=1,2
  35. XYC(IE2,IE1)=XYC(IE2,IE1-1)
  36. ENDDO
  37. DO IE2=1,3
  38. IST(IE2,IE1)=IST(IE2,IE1-1)
  39. ENDDO
  40. DENS(IE1)=DENS(IE1-1)
  41. ENDDO
  42. ENDIF
  43. C
  44. C ON INCERE LE NOUVEAU NOEUD
  45. C
  46. IST(1,INOEUD)=0
  47. IST(2,INOEUD)=0
  48. IST(3,INOEUD)=0
  49. DO IE1=1,2
  50. XYC(IE1,INOEUD)=XY(IE1)
  51. ENDDO
  52. DENS(INOEUD)=DENS1
  53. C
  54. RETURN
  55. END
  56.  
  57.  
  58.  

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