C DZVDT0    SOURCE    TTMF3     12/07/05    21:15:18     7425
      FUNCTION DZVDT0(P,T)
C--------------------------------------------------------------------
C Derivee partielle de Z par rapport a T
C Issue de la fonction ZVAP0(P,T)
C--------------------------------------------------------------------
C  E/  P    : Pression partielle de vapeur (en Pa)
C  E/  T    : Température (en K)
C   /S DZVDT0 : Derivee de facteur compressibilite Z / Pvap (en J/kg/K)
C               d/dT (Zvap(Pvap,T))
C--------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C      COEFFICIENTS

      B0 = -4512.70756D0
      B1 = 4245585.9D0
      B2 = -1584624342.D0
      R  = 461.523D0

C     premier terme de la formule
      VAL0  = ((2*B0*T+3*B1)*T+4*B2)*P
      VAL1=R*T*T*T*T*T
      VAL2=VAL0/VAL1


C     ECRITURE DE EPSILON(Pv,T)
C     *************************************************
      EPSILON= P*((B0*T+B1)*T+B2)
      EPSILON= (EPSILON/(T*T*T*T))+((R*R)/4)
      EPSILON = sqrt(EPSILON)

C     ECRITURE DE VAL3
C     *************************************************
      VAL3=(R/2)+EPSILON

C     ECRITURE DE VAL4
C     *************************************************
      S1 = ((B0*T+B1)*T+B2)*P
      S2 = 2*T*EPSILON*T*T*T
      S3 = (S1/S2)/(VAL3*VAL3)
      S4 = 1/(VAL3)
      VAL4=S3-S4

C     ECRITURE DE LA FORMULE BRUTE
C     *************************************************
      DZVDT0 =VAL4*VAL2
      RETURN
      END