Procedures | Cast3M

   1 : * RVST2     PROCEDUR  BLUMENFE  05/04/04    21:15:02     5052           
   2 : * RVST2     PROCEDUR        04/07/15              
   3 : *                                                                    ;
   4 : * Calcolo rho vapore alla saturazione per:  0 < T < 20
   5 : *                                           ------------
   6 : *                                                                    ;
   7 :  'DEBPROC' rvst2                                                     ;
   8 :    argu tn*chpoint                                                   ;
   9 :    si (((mini tn) < 0.0) ou ((maxi tn) > 20.0))                    ;
  10 :     mess '=========================================================' ;
  11 :     mess '===RVSAT: Interpolazione fuori limiti: 0 < T[C] < 20===' ;
  12 :     mess '=========================================================' ;
  13 :    finsi                                                             ;
  14 :    a0   =  4.86256716e-03                                            ;
  15 :    a1   =  2.98409029e-04                                            ;
  16 :    a2   =  2.10179146e-05                                            ;
  17 :    a3   = -7.96684413e-07                                            ;
  18 :    a4   =  2.77512067e-08                                            ;
  19 :    t1   =  kops       tn           '*' a1                            ;
  20 :    t2   =  kops (kops tn '**' 2.0) '*' a2                            ;
  21 :    t3   =  kops (kops tn '**' 3.0) '*' a3                            ;
  22 :    t4   =  kops (kops tn '**' 4.0) '*' a4                            ;
  23 :    rvs  =  kops a0 '+' (kops t1 '+' (kops t2 '+' (kops t3 '+' t4)))  ;
  24 :    rvs  =  kops rvs '|<' 0.0                                         ;
  25 : *
  26 :    s0   =  a1                                                        ;
  27 :    s1   =  kops       tn           '*' (a2 * 2.0)                    ;
  28 :    s2   =  kops (kops tn '**' 2.0) '*' (a3 * 3.0)                    ;
  29 :    s3   =  kops (kops tn '**' 3.0) '*' (a4 * 4.0)                    ;
  30 :    rvp  =  kops s0 '+' (kops s1 '+' (kops s2 '+' s3))                ;
  31 :    rvp  =  kops rvp '|<' 0.0                                         ;
  32 : *
  33 : *finproc rvs                                                         ;
  34 :  'FINPROC' rvs rvp                                                   ;
  35 : 
  36 : 
  37 :  
  38 :