Télécharger drvdt1.eso

Retour à la liste

Numérotation des lignes :

  1. C DRVDT1 SOURCE CHAT 05/01/12 22:55:49 5004
  2. SUBROUTINE DRVDT1(MPOVA1,MPOVA2,MPOVA3,ISAT)
  3. C--------------------------------------------------------------------
  4. C Densité de la vapeur d'eau en fonction de P et de T:
  5. C dérivée partielle par rapport à la température
  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 l'enthalpie vapeur spécifique à saturation : les data sont des CHPOS
  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 la densité de la vapeur d'eau
  16. C par rapport à la température
  17. C (en kg/m3/K)
  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 = DRVDT0(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 = DRVDT0(X1,X2)
  42. MPOVA3.VPOCHA(I,1) = 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