C J3POIN    SOURCE    CHAT      05/01/13    00:47:17     5004
      SUBROUTINE J3POIN(WORK,NPTO,INOEUD,XY,DENS1)
C----------------------------------------------------
C     INSERSION D'UN POINT
C
C     ESCLAVE DE J3COTO
C
C     CODE IST(1,I): 0 point non traite
C                    1 est sur le segment IST(2,I)
C                    2 est sur les segments IST(2,I) et IST(3,I)
C                   -1 est a l'interieur
C                   -2 est a l'exterieur
C
C     PP 6/97
C     Pierre Pegon/JRC Ispra
C----------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8(A-H,O-Z)
      DIMENSION XY(2)
C
      SEGMENT WORK
        REAL*8 XYC(2,NPTO)
        INTEGER IST(3,NPTO)
        REAL*8 DENS(NPTO)
        INTEGER JUN
      ENDSEGMENT
C
      SEGADJ,WORK
C
C     ON SHIFTE LE CONTENU DE WORK DE INOEUD->NPTO-1 A INOEUD+1->NPTO
C
      IF(NPTO.GT.INOEUD)THEN
        DO IE1=NPTO,INOEUD+1,-1
          DO IE2=1,2
            XYC(IE2,IE1)=XYC(IE2,IE1-1)
          ENDDO
          DO IE2=1,3
            IST(IE2,IE1)=IST(IE2,IE1-1)
          ENDDO
          DENS(IE1)=DENS(IE1-1)
        ENDDO
      ENDIF
C
C     ON INCERE LE NOUVEAU NOEUD
C
      IST(1,INOEUD)=0
      IST(2,INOEUD)=0
      IST(3,INOEUD)=0
      DO IE1=1,2
        XYC(IE1,INOEUD)=XY(IE1)
      ENDDO
      DENS(INOEUD)=DENS1
C
      RETURN
      END