Esope source files

   1 : C DHVDT0    SOURCE    SERRE     12/07/10    21:15:06     7431           
   2 :       FUNCTION DHVDT0(P,T)
   3 : 
   4 : C--------------------------------------------------------------------
   5 : C Derivee partielle de enthalpie specifique vapeur par rapport a la
   6 : C Temperature T
   7 : C Issue de la fonction enthalpie specifique vapeur HVS0(P,T)
   8 : C--------------------------------------------------------------------
   9 : C  E/  P    : Pression partielle de vapeur (en Pa)
  10 : C  E/  T    : Température (en K)
  11 : C   /S DHVDT0 : Derivee de enthalpie specifique (en J/kg/K)
  12 : C             d/dT (Hvap(Pvap,T)) 
  13 : C--------------------------------------------------------------------
  14 :       IMPLICIT INTEGER(I-N)
  15 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  16 : C
  17 : C     COEFFICIENTS                
  18 :       
  19 :       A1=1.748383D3
  20 :       A2=0.277768D0
  21 :       A3=0.275D-3
  22 :       
  23 :       
  24 :       B0 = -4512.70756D0
  25 :       B1 = 4245585.9D0
  26 :       B2 = -1584624342.D0 
  27 :       R  = 461.523D0
  28 : 
  29 : C     premier terme de la formule
  30 :       VAL0  = A1+(A3*T+A2)*T
  31 : 
  32 : C     ECRITURE DE EPSILON(Pv,T)     
  33 : C     *************************************************                       
  34 :       EPSILON= P*((B0*T+B1)*T+B2)
  35 :       EPSILON= (EPSILON/(T*T*T*T))+((R*R)/4)
  36 :       EPSILON = SQRT(EPSILON)
  37 : 
  38 : C     ECRITURE DE VAL1
  39 : C     *************************************************                 
  40 :       VAL1=(-R/2)+EPSILON                           
  41 :       
  42 : C     ECRITURE DE VAL2
  43 : C     *************************************************                                     
  44 :       S1=2*B0*T*T*T
  45 :       S2=3*B1*T*T
  46 :       S3=4*B2*T
  47 :       S4=S1+S2+S3
  48 :       S5= (B0*T+B1)*T+B2
  49 :       VAL2=S4/S5
  50 :       
  51 : C     ECRITURE DE BETA1    
  52 : C     *************************************************                 
  53 :       X1=(6*B0*T*T)+(6*B1*T)+4*B2
  54 :       X2= (B0*T*T)+(B1*T)+B2
  55 :       X3=(2*B0*T*T*T)+(3*B1*T*T)+(4*B2*T)
  56 :       X4=(2*B0*T)+B1
  57 :       BETA1=((X1*X2)-(X3*X4))/(X2*X2)
  58 :     
  59 : C     ECRITURE DE BETA2   
  60 : C     *************************************************                 
  61 :       Y1=(((2*B0*T+3*B1)*T)+4*B2)*(-P)
  62 :       Y2=2*EPSILON*(T*T*T*T*T)
  63 :       BETA2=Y1/Y2
  64 :      
  65 : C     ECRITURE DE LA FORMULE BRUTE     
  66 : C     *************************************************                         
  67 :       DHVDT0 =VAL0 + VAL1*BETA1+VAL2*BETA2 
  68 :       RETURN 
  69 :       END
  70 :  
  71 :  
  72 :  
  73 :  
  74 :