Esope source files

   1 : C DZVDT0    SOURCE    TTMF3     12/07/05    21:15:18     7425           
   2 :       FUNCTION DZVDT0(P,T)
   3 : C--------------------------------------------------------------------
   4 : C Derivee partielle de Z par rapport a T
   5 : C Issue de la fonction 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 DZVDT0 : Derivee de facteur compressibilite Z / Pvap (en J/kg/K)
  10 : C               d/dT (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  = ((2*B0*T+3*B1)*T+4*B2)*P
  24 :       VAL1=R*T*T*T*T*T
  25 :       VAL2=VAL0/VAL1
  26 : 
  27 : 
  28 : C     ECRITURE DE EPSILON(Pv,T)     
  29 : C     *************************************************                 
  30 :       EPSILON= P*((B0*T+B1)*T+B2)
  31 :       EPSILON= (EPSILON/(T*T*T*T))+((R*R)/4)
  32 :       EPSILON = sqrt(EPSILON)
  33 :             
  34 : C     ECRITURE DE VAL3 
  35 : C     *************************************************                 
  36 :       VAL3=(R/2)+EPSILON                           
  37 :       
  38 : C     ECRITURE DE VAL4
  39 : C     *************************************************                     
  40 :       S1 = ((B0*T+B1)*T+B2)*P
  41 :       S2 = 2*T*EPSILON*T*T*T
  42 :       S3 = (S1/S2)/(VAL3*VAL3)
  43 :       S4 = 1/(VAL3)
  44 :       VAL4=S3-S4
  45 :       
  46 : C     ECRITURE DE LA FORMULE BRUTE     
  47 : C     *************************************************                         
  48 :       DZVDT0 =VAL4*VAL2 
  49 :       RETURN
  50 :       END
  51 :  
  52 :  
  53 :  
  54 :