Procedures | Cast3M

   1 : * RVSAT     PROCEDUR  MAGN      99/12/09    21:15:19     3728           
   2 : *                                                                    ;
   3 : * Calcolo rho vapore alla saturazione per:  20 < T < 200
   4 : *                                           ------------
   5 : *                                                                    ;
   6 :  'DEBPROC' rvsat                                                     ;
   7 :    argu tn*chpoint                                                   ;
   8 :    si (((mini tn) < 20.0) ou ((maxi tn) > 200.0))                    ;
   9 :     mess '=========================================================' ;
  10 :     mess '===RVSAT: Interpolazione fuori limiti: 20 < T[C] < 200===' ;
  11 :     mess '=========================================================' ;
  12 :    finsi                                                             ;
  13 :    a0   =  4.18173362e-02                                            ;
  14 :    a1   = -2.58186762e-03                                            ;
  15 :    a2   =  8.65818031e-05                                            ;
  16 :    a3   = -7.08615945e-07                                            ;
  17 :    a4   =  6.58702608e-09                                            ;
  18 :    t1   =  kops       tn           '*' a1                            ;
  19 :    t2   =  kops (kops tn '**' 2.0) '*' a2                            ;
  20 :    t3   =  kops (kops tn '**' 3.0) '*' a3                            ;
  21 :    t4   =  kops (kops tn '**' 4.0) '*' a4                            ;
  22 :    rvs  =  kops a0 '+' (kops t1 '+' (kops t2 '+' (kops t3 '+' t4)))  ;
  23 :    rvs  =  kops rvs '|<' 0.0                                         ;
  24 : *
  25 :    s0   =  a1                                                        ;
  26 :    s1   =  kops       tn           '*' (a2 * 2.0)                    ;
  27 :    s2   =  kops (kops tn '**' 2.0) '*' (a3 * 3.0)                    ;
  28 :    s3   =  kops (kops tn '**' 3.0) '*' (a4 * 4.0)                    ;
  29 :    rvp  =  kops s0 '+' (kops s1 '+' (kops s2 '+' s3))                ;
  30 :    rvp  =  kops rvp '|<' 0.0                                         ;
  31 : *
  32 : *finproc rvs                                                         ;
  33 :  'FINPROC' rvs rvp                                                   ;
  34 : 
  35 : 
  36 :