Télécharger flmvl2.eso

Retour à la liste

Numérotation des lignes :

  1. C FLMVL2 SOURCE CHAT 05/01/13 00:03:52 5004
  2. C FLMVLH SOURCE BECC 97/12/23 22:13:47 3017
  3. SUBROUTINE FLMVL2(NESP,GAMD,ROD,MD,AD,UTD,UVD,HTD,Y,FLUX)
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : FLMVLH
  9. C
  10. C DESCRIPTION : Voir FLUVLH
  11. C
  12. C Onde de Droite à Gauche
  13. C
  14. C LANGUAGE : FORTRAN 77
  15. C
  16. C AUTEUR : A. BECCANTINI DRN/DMT/SEMT/TTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES : AUCUN
  21. C
  22. C************************************************************************
  23. C
  24. C**** Entrées:
  25. C
  26. C NESP = nombre d'especes in FRAMAS
  27. C
  28. C GAMD = le "gamma" (à droite)
  29. C
  30. C ROD = la densité
  31. C
  32. C MD = nombre de Mach (normale)
  33. C
  34. C AD = vitesse du son
  35. C
  36. C UTD = vitesse tangentielle
  37. C
  38. C Y = table des fractiones massiques
  39. C
  40. C**** Sorties:
  41. C
  42. C FLUX = table du flux "droite -> gauche" à
  43. C l'interface
  44. C
  45. C************************************************************************
  46. C
  47. C HISTORIQUE (Anomalies et modifications éventuelles)
  48. C
  49. C HISTORIQUE : Créé le 6.1.98
  50. C
  51. C************************************************************************
  52. C
  53. C N.B.: Toutes les variables sont declarées explicitement
  54. C
  55. CC
  56. IMPLICIT INTEGER(I-N)
  57. INTEGER I1, NESP
  58. REAL*8 GAMD,ROD,MD,AD,UTD,UVD,HTD
  59. & ,GM1D,MMINUS, CELL, CELLM
  60. & ,Y(*),FLUX(*)
  61. C
  62. C**** N.B. :
  63. C
  64. C FLUX(2) = MMINUS * AD * ROD * UND + PMINUS * PD
  65. C
  66. C PD = ROD * AD * AD / GAMD
  67. C
  68. C ABS(MD) < 1 -> PMINUS = - MMINUS * (2 + MD)
  69. C
  70. C ABS(MD) > 1 -> PMINUS = - MMINUS / MD
  71. C
  72. GM1D = GAMD - 1.0D0
  73. IF(ABS(MD) .LE. 1.0D0)THEN
  74. MMINUS = 0.5D0 * (MD - 1.0D0)
  75. MMINUS = -1.0D0 * MMINUS * MMINUS
  76. CELL = (GM1D * MD - 2.0D0)*AD/GAMD
  77. ELSE
  78. MMINUS = 0.5D0 * (MD - ABS(MD))
  79. CELL = (MD + 1.0D0/GAMD/MD)*AD
  80. ENDIF
  81. C
  82. CELLM = MMINUS * AD * ROD
  83. FLUX(1) = CELLM
  84. FLUX(2) = CELLM * CELL
  85. FLUX(3) = CELLM * UTD
  86. FLUX(4) = CELLM * UVD
  87. FLUX(5) = CELLM * HTD
  88. C
  89. C**** N.B. NESP = 0 -> no boucle
  90. C
  91. DO I1 = 1, NESP
  92. FLUX(5+I1) = CELLM * Y(I1)
  93. ENDDO
  94. C
  95. RETURN
  96. END
  97.  
  98.  
  99.  
  100.  

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