C RTADET    SOURCE    CB215821  17/11/30    21:17:12     9639           
      SUBROUTINE RTADET(XYZPT,IDIMC,ITRNOE,NBNMAX,ITRTRI,NBCMAX,
     >                     NBE,COORD,SPH,
     >                     IADET,NBADET,NADMAX,ZERO,iarr)
C     **********************************************************************
C     OBJET : RECHERCHE DES TRIANGLES A DETRUIRE (IE "NON-DELAUNAY")
C             A L'AJOUT D'UN POINT
C     EN ENTREE   :
C        XYZPT : COORDONNEES DU POINT AJOUTE
C        IDIMC  : DIMENSION DE L'ESPACE
C        ITRNOE,NBNMAX,ITRTRI,NBCMAX,NBE,COORD : LA TRIANGULATION
C        SPH   : LES SPHERES CIRCONSCRITES AUX TRIANGLES
C        NBADETMNAX : TAILLE DU TABLEAU IADET
C        ZERO  : PRECISION DU TEST "POINT DANS SPHERE"
C     EN SORTIE   :
C        IADET    : TABLEAU DES TRIANGLES "NON-DELAUNAY"
C        NBADET   : NOMBRE DE TRIANGLES "  "    "   "
C     **********************************************************************
      IMPLICIT INTEGER(I-N)
      REAL*8       XYZPT(*)
      INTEGER    IDIMC
      INTEGER    ITRNOE(*),NBNMAX,ITRTRI(*),NBCMAX
      INTEGER    NBE,IADET(*),NBADET,NADMAX,iarr
      REAL*8       COORD(*),SPH(*),ZERO
C
C
C     --- POUR LE DEBUG ---
C
C      COMMON /DEBUG/ ITRACE, ITEST, ierrOR, IMESS
C      INTEGER       ITRACE, ITEST, ierrOR
C      CHARACTER*256 IMESS
C     ---------------------------------------------------------------------
C     --- POUR LES STATS ---
C
C      COMMON /STATS/ ICARD(100)
C      INTEGER       ICARD
C     ---------------------------------------------------------------------
C     --- VARIABLES INTERNES ---
      INTEGER   I,NT,IPTDS,IPTDSC,IPTDS2,ITRACE
      REAL*8      SPHC(4)
      INTEGER  SPPOIN
      EXTERNAL SPPOIN
C
      ITRACE = 0
      DO i = 1, 4
        SPHC(i) = 0.D0
      ENDDO
      iarr = 0
      NBADET = 0
      DO 30 I=1,NBE
C     ----- ON PREND LE DERNIER NOEUD ---
        NT =ITRNOE((I-1)*NBNMAX+1+IDIMC)
        IF ( NT.EQ. 0 ) GO TO 30
        IPTDS = SPPOIN(IDIMC,COORD((NT-1)*IDIMC+1),XYZPT,
     >       SPH((I-1)*(IDIMC+1)+1),ZERO)
C
        IF( IPTDS.EQ.1 )THEN
          IADET(1) = I
          NBADET = 1
          CALL RTCONN(XYZPT,IDIMC,ITRNOE,NBNMAX,ITRTRI,NBCMAX,
     >                   COORD,SPH,IADET,NBADET,NADMAX,ZERO,iarr)
          RETURN
        ENDIF
C
   30 CONTINUE

      RETURN
      END