Télécharger dhvdp1.eso

Retour à la liste

Numérotation des lignes :

dhvdp1
  1. C DHVDP1 SOURCE CB215821 20/11/25 13:25:29 10792
  2. SUBROUTINE DHVDP1(MPOVA1,MPOVA2,MPOVA3,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/ MPOVA1 : Pression partielle de vapeur (en Pa)
  13. C E/ MPOVA2 : Température (en K)
  14. C E/ ISAT : 0 si T est donné, 1 sinon (T=TSAT(P))
  15. C /S MPOVA3 : 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 SMCHPOI
  28. C
  29. N1 = MPOVA1.VPOCHA(/1)
  30. IF (ISAT.EQ.0) THEN
  31. DO 10 I=1,N1
  32. X1 = MPOVA1.VPOCHA(I,1)
  33. X2 = MPOVA2.VPOCHA(I,1)
  34. X3 = DHVDP0(X1,X2)
  35. MPOVA3.VPOCHA(I,1) = X3
  36. 10 CONTINUE
  37. ELSE
  38. DO 20 I=1,N1
  39. X1 = MPOVA1.VPOCHA(I,1)
  40. X2 = TSATP0(X1)
  41. X3 = DHVDP0(X1,X2)
  42. MPOVA3.VPOCHA(I,1) = X3
  43. 20 CONTINUE
  44. ENDIF
  45. C
  46. RETURN
  47. END
  48.  
  49.  
  50.  
  51.  
  52.  

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