Télécharger centj1.eso

Retour à la liste

Numérotation des lignes :

  1. C CENTJ1 SOURCE CHAT 05/01/12 21:52:02 5004
  2. SUBROUTINE CENTJ1(RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY,
  3. & DFRO,DFRUN,DFRUT,DFRET)
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : CENTJ1
  9. C
  10. C DESCRIPTION : Calcul du Jacobien du flux numerique centré
  11. C Contribution de l'onde de gauche à droite
  12. C Voir CENTJ0
  13. C
  14. C LANGAGE : FORTRAN 77
  15. C
  16. C AUTEUR : A. BECCANTINI DRN/DMT/SEMT/LTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES : MPFUNT
  21. C
  22. C************************************************************************
  23. C
  24. C**** Entrées:
  25. C
  26. C RHO = densité
  27. C
  28. C UX, UY = vitesse
  29. C
  30. C P = pression
  31. C
  32. C GAM = le "gamma"
  33. C
  34. C CNX, CNY = normale à la face
  35. C
  36. C CTX, CTY = tangente à la face
  37. C
  38. C
  39. C**** Sorties:
  40. C
  41. C DFRO(4) = derivés partielles de F_RHO+ par raport a RHO, UX,
  42. C UY, P
  43. C
  44. C DFRUN(4) = derivés partielles de F_RHOUN+ par raport a RHO, UX,
  45. C UY, P
  46. C
  47. C DFRUT(4) = derivés partielles de F_RHOUT+ par raport a RHO, UX,
  48. C UY, P
  49. C
  50. C DFRET(4) = derivés partielles de F_RHOET+ par raport a RHO,
  51. C UX, UY, P
  52. C
  53. C************************************************************************
  54. C
  55. C HISTORIQUE (Anomalies et modifications éventuelles)
  56. C
  57. C HISTORIQUE : Créé le 04.09.00
  58. C
  59. C************************************************************************
  60. C
  61. C N.B.: Toutes les variables sont declarées explicitement
  62. C
  63. C
  64. IMPLICIT INTEGER(I-N)
  65. REAL*8 RHO,UX,UY,P,GAM,CNX,CNY,CTX,CTY
  66. & ,DFRO(4),DFRUN(4),DFRUT(4),DFRET(4)
  67. & ,UN,UT,CELL,GSGM1P,ECIN
  68. C
  69. UN = UX * CNX + UY * CNY
  70. UT = UX * CTX + UY * CTY
  71. CELL = RHO * UN
  72. C
  73. DFRO(1) = 0.5D0 * UN
  74. DFRO(2) = 0.5D0 * RHO * CNX
  75. DFRO(3) = 0.5D0 * RHO * CNY
  76. DFRO(4) = 0.0D0
  77. C
  78. DFRUN(1) = DFRO(1) * UN
  79. DFRUN(2) = CELL * CNX
  80. DFRUN(3) = CELL * CNY
  81. DFRUN(4) = 0.5D0
  82. C
  83. DFRUT(1) = DFRO(1) * UT
  84. DFRUT(2) = (DFRO(2) * UT) + (0.5D0 * CELL * CTX)
  85. DFRUT(3) = (DFRO(3) * UT) + (0.5D0 * CELL * CTY)
  86. DFRUT(4) = 0.0D0
  87. C
  88. GSGM1P=(GAM/(GAM-1.0D0))*P
  89. ECIN=0.5D0*((UX*UX)+(UY*UY))
  90. DFRET(1) = 0.5D0*ECIN*UN
  91. DFRET(2) = 0.5D0*((CNX*(GSGM1P+(RHO*ECIN)))+
  92. & (CELL*UX))
  93. DFRET(3) = 0.5D0*((CNY*(GSGM1P+(RHO*ECIN)))+
  94. & (CELL*UY))
  95. DFRET(4) = 0.5D0*UN*(GAM/(GAM-1.D0))
  96. C
  97. RETURN
  98. END
  99.  
  100.  
  101.  
  102.  
  103.  

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