Télécharger vlhj0.eso

Retour à la liste

Numérotation des lignes :

vlhj0
  1. C VLHJ0 SOURCE CHAT 05/01/13 04:07:44 5004
  2. SUBROUTINE VLHJ0(RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY,
  3. & DFROC,DFRUNC,DFRUTC,DFRETC)
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : VLHJ0
  9. C
  10. C DESCRIPTION : Calcul du Jacobien du flux numerique VLH
  11. C Contribution de l'onde de gauche à droite
  12. C
  13. C LANGAGE : FORTRAN 77
  14. C
  15. C AUTEUR : A. BECCANTINI DRN/DMT/SEMT/LTMF
  16. C
  17. C************************************************************************
  18. C
  19. C APPELES : VLHJ1
  20. C
  21. C************************************************************************
  22. C
  23. C**** Entrées:
  24. C
  25. C RHO = densité
  26. C
  27. C UX, UY = vitesse
  28. C
  29. C P = pression
  30. C
  31. C RHOET = rho * (ein + ecin)
  32. C
  33. C GAM = le "gamma"
  34. C
  35. C CNX, CNY = normale à la face
  36. C
  37. C CTX, CTY = tangente à la face
  38. C
  39. C
  40. C**** Sorties:
  41. C
  42. C DFROC(4) = derivés partielles de F_RHO+ par raport a RHO, RHOUX,
  43. C RHOUY, RHOET
  44. C
  45. C DFRUNC(4)= derivés partielles de F_RHOUN+ par raport a RHO, RHOUX,
  46. C RHOUY, RHOET
  47. C
  48. C DFRUNT(4)= derivés partielles de F_RHOUT+ par raport a RHO, RHOUX,
  49. C RHOUY, RHOET
  50. C
  51. C DFRETC(4)= derivés partielles de F_RHOET+ par raport a RHO,
  52. C RHOUX, RHOUY, RHOET
  53. C
  54. C************************************************************************
  55. C
  56. C HISTORIQUE (Anomalies et modifications éventuelles)
  57. C
  58. C HISTORIQUE : Créé le 04.09.00
  59. C
  60. C************************************************************************
  61. C
  62. C N.B.: Toutes les variables sont declarées explicitement
  63. C
  64. IMPLICIT INTEGER(I-N)
  65. REAL*8 RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY
  66. & ,DFROC(4),DFRUNC(4),DFRUTC(4),DFRETC(4)
  67. & ,DFROP(4),DFRUNP(4),DFRUTP(4),DFRETP(4)
  68. & ,USRO,UXSRO,UYSRO,GM1,GM1EC,GM1UX,GM1UY
  69. C
  70. C**** Jacobian with respect to the primitive variables
  71. C
  72. CALL VLHJ1(RHO,UX,UY,P,RHOET,GAM,CNX,CNY,CTX,CTY,
  73. & DFROP,DFRUNP,DFRUTP,DFRETP)
  74. C
  75. C**** Jacobian with respect to the conservative variables
  76. C
  77. USRO = 1.0D0 / RHO
  78. UXSRO = UX * USRO
  79. UYSRO = UY * USRO
  80. GM1 = GAM - 1.0D0
  81. GM1EC = GM1 * 0.5D0 * (UX * UX + UY * UY)
  82. GM1UX = GM1 * UX
  83. GM1UY = GM1 * UY
  84. C
  85. DFROC(1) = DFROP(1) - UXSRO * DFROP(2) - UYSRO * DFROP(3)
  86. & + GM1EC * DFROP(4)
  87. DFROC(2) = USRO * DFROP(2) - GM1UX * DFROP(4)
  88. DFROC(3) = USRO * DFROP(3) - GM1UY * DFROP(4)
  89. DFROC(4) = GM1 * DFROP(4)
  90. C
  91. DFRUNC(1) = DFRUNP(1) - UXSRO * DFRUNP(2) - UYSRO * DFRUNP(3)
  92. & + GM1EC * DFRUNP(4)
  93. DFRUNC(2) = USRO * DFRUNP(2) - GM1UX * DFRUNP(4)
  94. DFRUNC(3) = USRO * DFRUNP(3) - GM1UY * DFRUNP(4)
  95. DFRUNC(4) = GM1 * DFRUNP(4)
  96. C
  97. DFRUTC(1) = DFRUTP(1) - UXSRO * DFRUTP(2) - UYSRO * DFRUTP(3)
  98. & + GM1EC * DFRUTP(4)
  99. DFRUTC(2) = USRO * DFRUTP(2) - GM1UX * DFRUTP(4)
  100. DFRUTC(3) = USRO * DFRUTP(3) - GM1UY * DFRUTP(4)
  101. DFRUTC(4) = GM1 * DFRUTP(4)
  102. C
  103. DFRETC(1) = DFRETP(1) - UXSRO * DFRETP(2) - UYSRO * DFRETP(3)
  104. & + GM1EC * DFRETP(4)
  105. DFRETC(2) = USRO * DFRETP(2) - GM1UX * DFRETP(4)
  106. DFRETC(3) = USRO * DFRETP(3) - GM1UY * DFRETP(4)
  107. DFRETC(4) = GM1 * DFRETP(4)
  108. C
  109. RETURN
  110. END
  111.  
  112.  
  113.  
  114.  

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