Télécharger drvdt3.eso

Retour à la liste

Numérotation des lignes :

drvdt3
  1. C DRVDT3 SOURCE CHAT 05/01/12 22:56:05 5004
  2. SUBROUTINE DRVDT3(MLREE1,MLREE2,MLREE3,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/ 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 la densité de la vapeur d'eau
  16. C par rapport à la température (en kg/m3/K)
  17. C
  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 = DRVDT0(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 = DRVDT0(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