Esope source files

   1 : C DRVDT0    SOURCE    TTMF3     12/07/05    21:15:14     7425           
   2 :       FUNCTION DRVDT0(P,T)
   3 : C--------------------------------------------------------------------
   4 : C derivee partielle de la masse volumique de vapeur par rapport a T
   5 : C Issue de la fonction ROVAP0(P,T)
   6 : C--------------------------------------------------------------------
   7 : C  E/  P    : Pression partielle de vapeur (en Pa)
   8 : C  E/  T    : Température (en K)
   9 : C   /S DRVDT0 : Derivee partielle de la masse volumique de la vapeur
  10 : C               par  rapport a la temperature T   (en J/kg/K)
  11 : C               d/dT (RHOvap(Pvap,T)) 
  12 : C--------------------------------------------------------------------
  13 :       IMPLICIT INTEGER(I-N)
  14 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  15 : C
  16 : C      COEFFICIENTS          
  17 :          
  18 :       B0 = -4512.70756D0
  19 :       B1 = 4245585.9D0
  20 :       B2 = -1584624342.D0 
  21 :       R  = 461.523D0
  22 : 
  23 : C     ECRITURE DE EPSILON(Pv,T)     
  24 : C     *************************************************                   
  25 :       EPSILON= P*((B0*T+B1)*T + B2)
  26 :       EPSILON= EPSILON/(T*T*T*T)+((R*R)/4)
  27 :       EPSILON = sqrt(EPSILON)
  28 :       
  29 :       VAL0 = (R/2+EPSILON)
  30 :       VAL1=P*((2*B0*T+3*B1)*T+4*B2)
  31 :       VAL2=2*T*T*T*T*T*EPSILON*VAL0*VAL0
  32 :       VAL3=VAL1/VAL2
  33 :       VAL4=1/(T*VAL0)
  34 :       VAL5=VAL3-VAL4
  35 :       
  36 : C     ECRITURE DE LA FORMULE BRUTE     
  37 : C     *************************************************                         
  38 :       DRVDT0 = (VAL5*P)/T
  39 :       RETURN
  40 :       END
  41 :  
  42 :  
  43 :  
  44 :