Télécharger fstral.eso

Retour à la liste

Numérotation des lignes :

  1. C FSTRAL SOURCE BECC 11/05/26 21:15:29 6981
  2. SUBROUTINE FSTRAL(
  3. & PC, GAM,
  4. & P, RHO_L, P_L, U_L,
  5. & U, RHO, D,
  6. & LOGDEB, LOGAN)
  7. *
  8. *************************************************************************
  9. *
  10. * project : CAST3M, EUROPLEXUS...
  11. *
  12. * name : fstral
  13. *
  14. * description : euler equations for a mixture of stiffened gases
  15. * flux in the non-reactive case.
  16. *
  17. * left rarefaction curve as function of p
  18. *
  19. * language : fortran 77
  20. *
  21. * author : a. beccantini den/dm2s/sfme/ltmf
  22. *
  23. *************************************************************************
  24. *
  25. * called by :
  26. *
  27. *
  28. *************************************************************************
  29. *
  30. ***** input
  31. *
  32. * pc, gam = properties of the gas
  33. *
  34. * p = pressure in the right side of the left
  35. * rarefaction
  36. *
  37. * rho_l, p_l, u_l
  38. * = density, pressure, velocity on the left
  39. *
  40. * logdeb = debugging ?
  41. *
  42. ***** output
  43. *
  44. * u = velocity in the right side of the left
  45. * rarefaction
  46. *
  47. * d = rarefaction speed
  48. *
  49. * logan = anomaly ?
  50. *
  51. *************************************************************************
  52. *
  53. * 26/11/2009 created
  54. * 25/05/2011 evolution in CAST3M
  55. *
  56. *************************************************************************
  57. *
  58. * n.b.: all variables are declared
  59. *
  60. C IMPLICIT NONE
  61. IMPLICIT INTEGER(I-N)
  62. REAL*8
  63. & PC, GAM,
  64. & P, RHO_L, P_L, U_L,
  65. & U, RHO, D
  66. *
  67. REAL*8 VAR
  68. *
  69. * debugging ?
  70. *
  71. LOGICAL LOGDEB, LOGAN
  72. *
  73. * a left rarefaction in the x-frame is a right rarefaction in the
  74. * -x-frame !
  75. *
  76. VAR = -1.0D0 * U_L
  77. CALL FSTRAR(PC, GAM, P, RHO_L, P_L, VAR,
  78. & U, RHO, D, LOGDEB, LOGAN)
  79. U = -1.0D0 * U
  80. D = -1.0D0 * D
  81. *
  82. IF (LOGAN) THEN
  83. WRITE(*,*) 'SUBROUTINE FSTRAL.F'
  84. WRITE(*,*) 'ANOMALY DETECTED.'
  85. GOTO 9999
  86. ENDIF
  87. *
  88. 9999 RETURN
  89. END
  90.  
  91.  

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