Télécharger la plaquette Cast3M


fr en

Télécharger j3poin.eso

Retour à la liste

Numérotation des lignes : 

j3poin

  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