Télécharger la plaquette Cast3M


fr en

Télécharger astiff.eso

Retour à la liste

Numérotation des lignes : 

astiff

  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