Esope source files

   1 : C TRJTET    SOURCE    CHAT      05/01/13    03:51:19     5004           
   2 :       SUBROUTINE TRJTET(IZNOEU,IZTRAV,IZAPAR,NPAPAR,J,IEL,INOELO,TLI)
   3 : C                                                                               
   4 : C**********************************************************************         
   5 : C                                                                               
   6 : C   OBJET  : CE SOUS-PROGRAMME CALCULE L'APPARTENANCE DE POINTS A UN            
   7 : C   -----    ELEMENT TETRAEDRE D'UN MAILLAGE M1.
   8 : C                                                                               
   9 : C   ARGUMENTS:                                                                  
  10 : C   ---------                                                                   
  11 : C                                                                               
  12 : C   ENTREE  :  IZNOEU = POINTEUR DU SEGMENT CONTENANT LES COORDONNEES           
  13 : C                       ET LES NUMEROS DES NOEUDS DE L ELEMENT COURANT          
  14 : C              IZTRAV = POINTEUR DU SEGMENT CONTENANT LES COORDONNEES           
  15 : C                       DES POINTS DONT ON CHERCHE L'APPARTENANCE               
  16 : C              IZAPAR = POINTEUR DU SEGMENT QUI CONTIENDRA LES                  
  17 : C                       APPARTENANCES                                           
  18 : C              NPAPAR = NOMBRE DE POINTS APPARTENANT A UN ELEMENT               
  19 : C              IEL    = NUMERO GLOBAL DE L'ELEMENT DANS LE MAILLAGE             
  20 : C              J      = NUMERO DU POINT TRAITE                                  
  21 : C                                                                               
  22 : C***********************************************************************        
  23 : C                                                                               
  24 :       IMPLICIT INTEGER(I-N)
  25 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  26 :       DIMENSION TLI(*)
  27 :       DIMENSION X(3)
  28 : C                                                                               
  29 : C                                                                               
  30 :       SEGMENT IZAPAR 
  31 :            INTEGER IAPAR(4,NPT2)
  32 :       ENDSEGMENT
  33 : C                                                                               
  34 :       SEGMENT IZTRAV
  35 :            REAL*8  COOR(NDIM,NPART)
  36 :       ENDSEGMENT
  37 :       SEGMENT  IZNOEU 
  38 :             REAL*8 XELE(IDIM,NOEL)  
  39 :             INTEGER NOEGLO(NOEL)
  40 :       ENDSEGMENT
  41 : C                                                                               
  42 : C                                                                               
  43 : C                                                                               
  44 : C                                                                               
  45 :       IDIM=XELE(/1)
  46 :       NOEL=XELE(/2)
  47 :       INOELO=0
  48 : C                                                                               
  49 : C                                                                               
  50 :       IF(IAPAR(1,J).NE.0) GO TO 9999                                            
  51 :       X(1)=COOR(1,J)                                                            
  52 :       X(2)=COOR(2,J)                                                            
  53 :       X(3)=COOR(3,J)                                                            
  53 :       CALL APATET(X,XELE,IELEM,IFATET,IARTET,INOTET,TLI)
  54 :       IF(IELEM.EQ.0) GO TO 9999
  55 : C                                                                               
  56 : C---  LE POINT APPARTIENT A L'ELEMENT                                           
  57 : C                                                                               
  58 :       NPAPAR=NPAPAR+1                                                           
  59 :       IAPAR(1,J)=IEL                                                            
  60 : C                                                                               
  61 : C---  APPARTENANCE A UNE FACE DU TETRAEDRE
  62 : C                                                                               
  63 :       IF(IFATET.EQ.0) GO TO 9999                                                
  64 :       IAPAR(2,J)=IFATET
  65 : C                                                                               
  66 : C---  APPARTENANCE A UNE ARETE DU TETRAEDRE
  67 : C                                                                               
  68 :       IF(IARTET.EQ.0) GO TO 9999                                                
  69 :       IAPAR(3,J)=IARTET
  70 : C                                                                               
  71 : C---  APPARTENANCE A UN NOEUD DU TETRAEDRE
  72 : C                                                                               
  73 :       IF(INOTET.EQ.0) GO TO 9999                                                
  74 :       INOELO = INOTET
  75 :       IAPAR(4,J)=NOEGLO(INOTET)
  76 :  9999 CONTINUE
  77 :       RETURN
  78 :       END                                                                       
  79 :  
  80 :