Télécharger trjnoe.eso

Retour à la liste

Numérotation des lignes :

trjnoe
  1. C TRJNOE SOURCE CHAT 05/01/13 03:50:16 5004
  2. SUBROUTINE TRJNOE(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL1,IEL,
  3. * INOELO,IZPART,IZUN,IZCOU,ITP,IFORML,XY)
  4. C
  5. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  6. C
  7. C LORSQU UNE PARTICULE EST LACHEE SUR UN NOEUD DU MAILLAGE
  8. C GEOMETRIQUE ON DETERMINE SI ELLE EST OU NON DANS L ELEMENT IEL1
  9. C
  10. C IZUN VITESSE OU FLUX
  11. C XY COORDONNEES DE REFERENCE DU POINT CONSIDERE
  12. C
  13. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  14. C
  15. C
  16. IMPLICIT INTEGER(I-N)
  17. IMPLICIT REAL*8 (A-H,O-Z)
  18. C
  19.  
  20. -INC PPARAM
  21. -INC CCOPTIO
  22. C
  23. C
  24. SEGMENT IZPART
  25. INTEGER NLEPA(NPART),NUMPA(NPART)
  26. REAL*8 COORPA(NDIM,NPART)
  27. ENDSEGMENT
  28. POINTEUR IZREF.IZPART
  29. SEGMENT IZCOU
  30. REAL*8 DTCO(NEL),COU
  31. ENDSEGMENT
  32. C
  33. SEGMENT IZTRAV
  34. REAL*8 COOR(NDIM,NPART)
  35. ENDSEGMENT
  36. SEGMENT IZNOEU
  37. REAL*8 XELE(IDIM,NOEL)
  38. INTEGER NOEGLO(NOEL)
  39. ENDSEGMENT
  40. SEGMENT IZAPAR
  41. INTEGER IAPAR(4,NPART)
  42. ENDSEGMENT
  43. SEGMENT IZUN
  44. REAL*8 UN(I1,I2,I3)
  45. ENDSEGMENT
  46. SEGMENT HYTRAV
  47. REAL*8 XEL(3,NOEL),SHP(6,NOEL),SHY(12,NOEL)
  48. REAL*8 ZJAC(NDIM,NDIM)
  49. ENDSEGMENT
  50. C
  51. C
  52. DIMENSION XY(3),UELEM(3),VECT(3)
  53. DIMENSION TLI(4)
  54. C
  55. C
  56. NDIM=COORPA(/1)
  57. COUR=0.001D0
  58. INOEL1= INOELO
  59. IF(IFORML.EQ.1)THEN
  60. C FORMULATION ELEMENTS FINIS
  61. DO 10 ID=1,NDIM
  62. UELEM(ID)=UN(ID,INOELO,IEL)
  63. 10 CONTINUE
  64. ELSEIF(IFORML.EQ.2)THEN
  65. C FORMULATION HYBRIDE ON VA CALCULER LA VITESSE
  66. C REELLE AU NOEUD DU MAILLAGE
  67.  
  68. CALL VHYNOE(UELEM,IZNOEU,IZUN,NDIM,XY,ITP,IEL)
  69. ENDIF
  70. DO 70 ID=1,NDIM
  71. COOR(ID,IPART)=COORPA(ID,IPART)+
  72. * UELEM(ID)*DTCO(IEL1)*COUR
  73. C write(6,*)' coor trjnoe ',COOR(ID,IPART),COORPA(ID,IPART),
  74. C * UELEM(ID),DTCO(IEL1),COUR
  75. 70 CONTINUE
  76. CALL INITI(IAPAR(1,IPART),4,0)
  77. C*** TRIANGLES
  78. C
  79. IF(ITP.EQ.4.OR.ITP.EQ.6.OR.ITP.EQ.7)THEN
  80. CALL TRJTRI(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL1,
  81. * INOELO,TLI)
  82. C
  83. C*** QUADRANGLES
  84. C
  85. ELSEIF(ITP.EQ.8.OR.ITP.EQ.11)THEN
  86. CALL TRJQUA(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL1,
  87. * INOELO,TLI)
  88. C*** PRISMES
  89. C
  90. ELSEIF(ITP.EQ.16)THEN
  91. CALL TRJPRI(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL,
  92. * INOELO,TLI,ITRI)
  93. C
  94. C*** CUBES
  95. C
  96. ELSEIF(ITP.EQ.14)THEN
  97. CALL TRJCUB(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL,
  98. * INOELO,TLI)
  99. C
  100. C*** TETRAEDRE
  101. C
  102. ELSEIF(ITP.EQ.23)THEN
  103. CALL TRJTET(IZNOEU,IZTRAV,IZAPAR,NPAPAR,IPART,IEL,
  104. * INOELO,TLI)
  105. ELSE
  106. CALL ERREUR(16)
  107. ENDIF
  108. DO 80 ID=1,NDIM
  109. COOR(ID,IPART)=COORPA(ID,IPART)
  110. 80 CONTINUE
  111. INOELO=INOEL1
  112. C
  113. END
  114.  
  115.  
  116.  

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