Télécharger la plaquette Cast3M


fr en

Télécharger flmvl2.eso

Retour à la liste

Numérotation des lignes : 

flmvl2

  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