Télécharger derr.eso

Retour à la liste

Numérotation des lignes :

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

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