Télécharger dhvdp3.eso

Retour à la liste

Numérotation des lignes :

dhvdp3
  1. C DHVDP3 SOURCE CHAT 05/01/12 22:49:08 5004
  2. SUBROUTINE DHVDP3(MLREE1,MLREE2,MLREE3,ISAT)
  3. C--------------------------------------------------------------------
  4. C Enthalpie vapeur spécifique en fonction de P et de T:
  5. C dérivée partielle par rapport à la pression partielle de vapeur
  6. C
  7. C Dans le cas ou seule la pression est fournie, on considère que la
  8. C pression P transmise est Psat et on calcule Tsat afin d'évaluer
  9. C la fonction derivée
  10. C--------------------------------------------------------------------
  11. C
  12. C E/ MLREE1 : Pression partielle de vapeur (en Pa)
  13. C E/ MLREE2 : Température (en K)
  14. C E/ ISAT : 0 si T est donné, 1 sinon (T=TSAT(P))
  15. C /S MLREE3 : Dérivée partielle de l'enthalpie spécifique de la
  16. C vapeur par rapport à la pression partielle de
  17. C vapeur (en J/kg/Pa)
  18. C-----------------------------------------------------------------------
  19. C
  20. C Langage : ESOPE + FORTRAN77
  21. C 10/2002
  22. C
  23. C-----------------------------------------------------------------------
  24. IMPLICIT INTEGER(I-N)
  25. IMPLICIT REAL*8 (A-H,O-Z)
  26. C
  27. -INC SMLREEL
  28. C
  29. JG1 = MLREE1.PROG(/1)
  30. IF (ISAT.EQ.0) THEN
  31. DO 10 I=1,JG1
  32. X1 = MLREE1.PROG(I)
  33. X2 = MLREE2.PROG(I)
  34. X3 = DHVDP0(X1,X2)
  35. MLREE3.PROG(I) = X3
  36. 10 CONTINUE
  37. ELSE
  38. DO 20 I=1,JG1
  39. X1 = MLREE1.PROG(I)
  40. X2 = TSATP0(X1)
  41. X3 = DHVDP0(X1,X2)
  42. MLREE3.PROG(I) = X3
  43. 20 CONTINUE
  44. ENDIF
  45. C
  46. RETURN
  47. END
  48.  
  49.  
  50.  
  51.  

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