Télécharger la plaquette Cast3M


fr en

Télécharger trjtet.eso

Retour à la liste

Numérotation des lignes : 

trjtet

  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)
  54.       CALL APATET(X,XELE,IELEM,IFATET,IARTET,INOTET,TLI)
  55.       IF(IELEM.EQ.0) GO TO 9999
  56. C
  57. C---  LE POINT APPARTIENT A L'ELEMENT
  58. C
  59.       NPAPAR=NPAPAR+1
  60.       IAPAR(1,J)=IEL
  61. C
  62. C---  APPARTENANCE A UNE FACE DU TETRAEDRE
  63. C
  64.       IF(IFATET.EQ.0) GO TO 9999
  65.       IAPAR(2,J)=IFATET
  66. C
  67. C---  APPARTENANCE A UNE ARETE DU TETRAEDRE
  68. C
  69.       IF(IARTET.EQ.0) GO TO 9999
  70.       IAPAR(3,J)=IARTET
  71. C
  72. C---  APPARTENANCE A UN NOEUD DU TETRAEDRE
  73. C
  74.       IF(INOTET.EQ.0) GO TO 9999
  75.       INOELO = INOTET
  76.       IAPAR(4,J)=NOEGLO(INOTET)
  77.  9999 CONTINUE
  78.       RETURN
  79.       END
  80.  
  81.  
  82.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales