Esope source files

   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 :