Télécharger la plaquette Cast3M


fr en

Télécharger fstrar.eso

Retour à la liste

Numérotation des lignes : 

fstrar

  1. C FSTRAR    SOURCE    BECC      11/05/26    21:15:30     6981
  2.       SUBROUTINE FSTRAR(
  3.      &     PC, GAM,
  4.      &     P, RHO_R, P_R, U_R,
  5.      &     U, RHO, D,
  6.      &     LOGDEB, LOGAN)
  7. *
  8. *************************************************************************
  9. *
  10. * project           :  CAST3M, EUROPLEXUS...
  11. *
  12. * name              :  fstrar
  13. *
  14. * description       :  euler equations for a mixture of stiffened gases
  15. *                      flux in the non-reactive case.
  16. *
  17. *                      right 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 left side of the right
  35. *                       rarefaction
  36. *
  37. *     rho_r, p_r, u_r
  38. *                     = density, pressure, velocity on the right
  39. *
  40. *     logdeb          = debugging ?
  41. *
  42. ***** output
  43. *
  44. *     u, rho          = velocity and density in the left side of the
  45. *                       right 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_R, P_R, U_R,
  65.      &     U, D
  66. *
  67.       REAL*8 A_R, RHO, A, COEF
  68. *
  69. *     debugging ?
  70. *
  71.       LOGICAL LOGDEB, LOGAN
  72. *
  73.       CALL ASTIFF (RHO_R, P_R, GAM, PC, A_R, LOGDEB, LOGAN)
  74. *      write(*,*) 'ciao', logan, a_r
  75.       IF (LOGAN) GOTO 9999
  76. *
  77. *     density (isentropic transformation)
  78. *
  79. *     a negative p here is catched by astiff after...
  80. *
  81.       RHO = (P + PC) / (P_R + PC)
  82.       RHO = RHO ** (1.0D0 / GAM)
  83.       RHO = RHO * RHO_R
  84. *
  85.       CALL ASTIFF (RHO, P, GAM, PC, A, LOGDEB, LOGAN)
  86.       IF (LOGAN) GOTO 9999
  87. *
  88.       COEF =  2.0D0 / (GAM - 1.0D0)
  89.       U = U_R + (COEF * (A - A_R))
  90. *
  91.       D = U + A
  92. *
  93.  9999 RETURN
  94.       END
  95.  
  96.  
  97.  

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