Esope source files

   1 : C DZVDP0    SOURCE    TTMF3     12/07/05    21:15:17     7425           
   2 :       FUNCTION DZVDP0(P,T)
   3 : C--------------------------------------------------------------------
   4 : C derivee partielle de Z par rapport a P
   5 : C Issue de la fonction compressibilite ZVAP0(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 DZVDP0 :Derivee de facteur compressibilite Z / Pvap (en J/kg/Pa)
  10 : C              d/dP (Zvap(Pvap,T)) 
  11 : C--------------------------------------------------------------------
  12 :       IMPLICIT INTEGER(I-N)
  13 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  14 : C
  15 : C      COEFFICIENTS                 
  16 :          
  17 :       B0 = -4512.70756D0
  18 :       B1 = 4245585.9D0
  19 :       B2 = -1584624342.D0 
  20 :       R  = 461.523D0
  21 : 
  22 : C     premier terme de la formule
  23 :       VAL0  =(B0*T+B1)*T+B2
  24 :       VAL1=R*T*T*T*T
  25 :       VAL2=VAL0/VAL1
  26 : 
  27 : C     ECRITURE DE EPSILON(Pv,T)     
  28 : C     *************************************************                 
  29 :       EPSILON= P*((B0*T+B1)*T+B2)
  30 :       EPSILON= (EPSILON/(T*T*T*T)+((R*R)/4))
  31 :       EPSILON = SQRT(EPSILON)  
  32 :      
  33 : C     ECRITURE DE VAL3 
  34 : C     *************************************************                 
  35 :       VAL3=(R/2)+EPSILON                           
  36 :       
  37 : C     ECRITURE DE VAL4
  38 : C     *************************************************                              
  39 :       S1=(VAL0*P)/(EPSILON*T*T*T*T*2)
  40 :       S2=VAL3-S1
  41 :       S3=S2/(VAL3*VAL3)
  42 :       
  43 : C     ECRITURE DE LA FORMULE BRUTE     
  44 : C     *************************************************                         
  45 :       DZVDP0 =S3*VAL2 
  46 :       RETURN
  47 :       END
  48 :  
  49 :  
  50 :  
  51 :