Télécharger hvs0.eso

Retour à la liste

Numérotation des lignes :

  1. C HVS0 SOURCE SERRE 16/11/25 21:15:00 9165
  2. FUNCTION HVS0(P,T)
  3. C
  4. C--------------------------------------------------------------------
  5. * Enthalpie vapeur spécifique fonction de P et de T (issue de RALOC).
  6. C--------------------------------------------------------------------
  7. C
  8. C E/ P : Pression partielle de vapeur (en Pa)
  9. C E/ T : Température (en K)
  10. C /S HVS0 : Enthalpie spécifique de la vapeur (en J/kg)
  11. C
  12. C--------------------------------------------------------------------
  13. IMPLICIT INTEGER(I-N)
  14. IMPLICIT REAL*8 (A-H,O-Z)
  15. C
  16. A0 = 2009.9595D3
  17. A1 = 1.748383D3
  18. A2 = 0.277768D0/2.D0
  19. A3 = 0.275D-3/3.D0
  20. B0 = -4512.70756D0
  21. B1 = 4245585.9D0
  22. B2 = -1584624342.D0
  23. R = 461.523D0
  24. C
  25. RAC1 = ((B0*T + B1) * T + B2) * P / (T**4.D0) + (R*R/4.D0)
  26. IF (RAC1 .GT. 0.D0) THEN
  27. RAC1 = RAC1 ** 0.5D0
  28. ELSE
  29. RAC1 = 0.D0
  30. ENDIF
  31. ZVAP1 = RAC1/R + 0.5D0
  32. C
  33. DNUM1 = (2.D0*B0*T + (3.D0*B1)) * T + 4.D0*B2
  34. DENOM1 = (B0*T + B1) * T + B2
  35. VPMRT1 = R * T * (ZVAP1 - 1.D0)
  36. C
  37. HVS0 = ((A3*T+A2)*T+A1)*T+A0 + DNUM1/DENOM1*VPMRT1
  38. C
  39. END
  40.  
  41.  
  42.  
  43.  
  44.  

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