Télécharger cubksi.eso

Retour à la liste

Numérotation des lignes :

  1. C CUBKSI SOURCE CHAT 07/10/22 21:15:34 5921
  2. SUBROUTINE CUBKSI(NPINT,XXG,BKSI)
  3. implicit real*8(A-H,O-Z)
  4. implicit integer (I-N)
  5.  
  6. * SUBROUTINE CUB_KSI(NPINT,XXG,BKSI)
  7. *
  8. * ------------------------------------------------------------------
  9. *
  10. * DERIVEE P/R VARIABLES CANONIQUE AUX POINTS D INTEGRATION
  11. * POUR LES ELEMENTS CUBES A 8 NOEUDS
  12. * H. BUNG 03-93
  13. * ------------------------------------------------------------------
  14. *
  15. *
  16. * ENTREES
  17. * NPINT : NBRE DE PTS D INTEGRATION
  18. * XG : COOR. CANONIQUES DES PTS D INTEGRATION
  19. *
  20. * SORTIES :
  21. * BKSI(3,8,NPINT) : LES DERIVEES
  22. *
  23. * IMPLICIT NONE
  24. *
  25. *--- VARIABLES GLOBALES
  26. *
  27. * INTEGER NPINT
  28. * REAL *8 XG(3,5),XXG(5),GS(8,4)
  29. * REAL *8 BKSI(3,8,5)
  30. dimension XG(3,5),XXG(*),GS(8,4),BKSI(3,8,*)
  31.  
  32. *
  33. *--- VARIABLES LOCALES
  34. *
  35. * REAL *8 AUX
  36. * REAL *8 UN,UNS8
  37. * INTEGER IP,I
  38. dimension XP(3),XM(3)
  39.  
  40. UN=1.D0
  41. UNS8=0.125D0
  42. *
  43. C MODIF ANTOINE:
  44. C ON REMPLIT LE TABLEAU XG
  45. C
  46. DO IP=1,NPINT
  47. XG(1,IP)=0.d0
  48. XG(2,IP)=0.d0
  49. XG(3,IP)=XXG(IP)
  50. END DO
  51. *
  52. *
  53. DO IP=1,NPINT
  54. *-
  55. DO I=1,3
  56. XP(I) = UN + XG(I,IP)
  57. XM(I) = UN - XG(I,IP)
  58. END DO
  59. *-
  60. AUX= -UNS8*XM(3)
  61. BKSI(1,1,IP)= AUX*XM(2)
  62. BKSI(1,2,IP)=-BKSI(1,1,IP)
  63. BKSI(1,3,IP)=-AUX*XP(2)
  64. BKSI(1,4,IP)=-BKSI(1,3,IP)
  65. BKSI(2,1,IP)= AUX*XM(1)
  66. BKSI(2,2,IP)= AUX*XP(1)
  67. BKSI(2,3,IP)=-BKSI(2,2,IP)
  68. BKSI(2,4,IP)=-BKSI(2,1,IP)
  69. AUX= -UNS8*XP(3)
  70. BKSI(1,5,IP)= AUX*XM(2)
  71. BKSI(1,6,IP)=-BKSI(1,5,IP)
  72. BKSI(1,7,IP)=-AUX*XP(2)
  73. BKSI(1,8,IP)=-BKSI(1,7,IP)
  74. BKSI(2,5,IP)= AUX*XM(1)
  75. BKSI(2,6,IP)= AUX*XP(1)
  76. BKSI(2,7,IP)=-BKSI(2,6,IP)
  77. BKSI(2,8,IP)=-BKSI(2,5,IP)
  78. AUX= -UNS8*XM(2)
  79. BKSI(3,1,IP)= AUX*XM(1)
  80. BKSI(3,2,IP)= AUX*XP(1)
  81. BKSI(3,5,IP)=-BKSI(3,1,IP)
  82. BKSI(3,6,IP)=-BKSI(3,2,IP)
  83. AUX= -UNS8*XP(2)
  84. BKSI(3,3,IP)= AUX*XP(1)
  85. BKSI(3,4,IP)= AUX*XM(1)
  86. BKSI(3,7,IP)=-BKSI(3,3,IP)
  87. BKSI(3,8,IP)=-BKSI(3,4,IP)
  88. *
  89. END DO
  90. *
  91. END
  92.  
  93.  
  94.  
  95.  

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