Procedures | Cast3M

   1 : * PSATT     PROCEDUR  MAGN      99/12/07    21:15:11     3724           
   2 : 'DEBPROC' PSATT ;
   3 : 'ARGU' T/'CHPOINT'
   4 :        TF/'FLOTTANT'
   5 :        TL/'LISTREEL' ;
   6 : *
   7 : *--------------------------------------------------------------------
   8 : * Température de saturation de la vapeur pour une pression partielle
   9 : * de vapeur donnée.
  10 : *--------------------------------------------------------------------
  11 : *  E/   T    : Température de la vapeur (en K)
  12 : *   /S  PS1  : PSAT(T) (en Pa)
  13 : *--------------------------------------------------------------------
  14 : * La donnée d'entrée est un CHPOINT, un FLOTTANT ou un LISTREEL.
  15 : * Le résultat est du meme type que la donnée d'entrée.
  16 : *--------------------------------------------------------------------
  17 : *
  18 : A0 = 4.3820662D-8 ;
  19 : A1 = 1.9770073D-6 ;
  20 : A2 = 6.7955758D-5 ;
  21 : A3 =-9.3874734D-7 ;
  22 : A4 = 1.2653874D-8 ;
  23 : 'SI' ('EXIST' T) ;
  24 :      T1  = T - 273.15 ;
  25 :      PS1 = (((A4*T1+A3)*T1+A2)*T1+A1)*T1+A0 ;
  26 :      PS1 = PS1 * 1.E+5 ;
  27 :     'RESPROC' PS1 ;
  28 :     'QUITTER' PSATT ;
  29 : 'FINSI' ;
  30 : 'SI' ('EXIST' TF) ;
  31 :      T1  = TF - 273.15 ;
  32 :      PS1 = (((A4*T1+A3)*T1+A2)*T1+A1)*T1+A0 ;
  33 :      PS1 = PS1 * 1.E+5 ;
  34 :     'RESPROC' PS1 ;
  35 :     'QUITTER' PSATT ;
  36 : 'FINSI' ;
  37 : 'SI' ('EXIST' TL) ;
  38 :      T00 = TL / TL  ;
  39 :      A00 = A0 * T00 ;
  40 :      A10 = A1 * T00 ;
  41 :      A20 = A2 * T00 ;
  42 :      A30 = A3 * T00 ;
  43 :      A40 = A4 * T00 ;
  44 :      T1  = TL - (273.15 * T00) ;
  45 :      PS1 = (((A40*T1+A30)*T1+A20)*T1+A10)*T1+A00 ;
  46 :      PS1 = PS1 * 1.E+5 ;
  47 :     'RESPROC' PS1 ;
  48 :     'QUITTER' PSATT ;
  49 : 'FINSI' ;
  50 : 'FINPROC' ;
  51 :