Esope source files

   1 : C DRVDP0    SOURCE    TTMF3     12/07/05    21:15:13     7425           
   2 :       FUNCTION DRVDP0(P,T)
   3 : C--------------------------------------------------------------------
   4 : C derivee partielle de la masse volumique de vapeur par rapport Pvap
   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 DRVDP0 : Derivee partielle de la masse volumique de la vapeur 
  10 : C               par rapport a la pression de vapeur Pvap   (en J/kg/K)
  11 : C               d/dP (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 = 1/(T*VAL0)
  31 :       VAL2 = P*((B0*T+B1)*T+B2)
  32 :       VAL3 = 2*EPSILON*T*T*T*T*VAL0
  33 :       VAL4 = 1-(VAL2/VAL3)
  34 :       VAL5 = VAL4*VAL1
  35 :       
  36 : C     ECRITURE DE LA FORMULE BRUTE     
  37 : C     *************************************************                         
  38 :       DRVDP0 = VAL5
  39 :       RETURN
  40 :       END
  41 :  
  42 :  
  43 :  
  44 :