Télécharger astiff.eso

Retour à la liste

Numérotation des lignes :

  1. C ASTIFF SOURCE BECC 11/05/26 21:15:01 6981
  2. SUBROUTINE ASTIFF(
  3. & RHO, P, GAMMA, PC, A,
  4. & LOGDEB, LOGAN)
  5. *
  6. *************************************************************************
  7. *
  8. * project : CAST3M; Europlexus...
  9. *
  10. * name : a_stiff
  11. *
  12. * description : euler equations for a mixture of stiffened gases
  13. * flux in the non-reactive case.
  14. *
  15. * computation of the sound speed.
  16. *
  17. * language : fortran 77
  18. *
  19. * author : a. beccantini den/dm2s/sfme/ltmf
  20. *
  21. *************************************************************************
  22. *
  23. * called by :
  24. *
  25. *
  26. *************************************************************************
  27. *
  28. ***** input
  29. *
  30. * pc, gam = properties of the stiffened gas
  31. *
  32. * rho, p, = density, pressure
  33. *
  34. * logdeb = debugging ?
  35. *
  36. ***** output
  37. *
  38. * a = sound speed
  39. *
  40. * logan = anomaly ?
  41. *
  42. *
  43. *************************************************************************
  44. *
  45. * 26/11/2009 created
  46. * 25/05/2011 evolution in CAST3M
  47. *
  48. *************************************************************************
  49. *
  50. * n.b.: all variables are declared
  51. *
  52. C IMPLICIT NONE
  53. IMPLICIT INTEGER(I-N)
  54. REAL*8 RHO, P, PC, GAMMA, A
  55. *
  56. LOGICAL LOGDEB, LOGAN
  57. *
  58. IF (LOGDEB ) THEN
  59. IF (GAMMA .LE. 1.0D0) THEN
  60. WRITE(*,*) 'GAMMA =', GAMMA
  61. LOGAN = .TRUE.
  62. GOTO 9999
  63. ENDIF
  64. IF ((P .LE. 0.0D0) .OR. (PC .LE. 0.0D0)) THEN
  65. WRITE(*,*) 'P = ', P, ' PC = ', PC
  66. LOGAN = .TRUE.
  67. GOTO 9999
  68. ENDIF
  69. IF ((RHO .LE. 0.0D0)) THEN
  70. WRITE(*,*) 'RHO = ', RHO
  71. LOGAN = .TRUE.
  72. GOTO 9999
  73. ENDIF
  74. ENDIF
  75. A = GAMMA * (P + PC) / RHO
  76. A = A ** 0.5D0
  77. *
  78. 9999 RETURN
  79. END
  80.  
  81.  

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