Télécharger la plaquette Cast3M


fr en

Télécharger rvst2.procedur

Retour à la liste

Numérotation des lignes : 

  1. * RVST2     PROCEDUR  BLUMENFE  05/04/04    21:15:02     5052           
  2. *                                                                    ;
  3. * Calcolo rho vapore alla saturazione per:  0 < T < 20
  4. *                                           ------------
  5. *                                                                    ;
  6.  'DEBPROC' rvst2                                                     ;
  7.    argu tn*chpoint                                                   ;
  8.    si (((mini tn) < 0.0) ou ((maxi tn) > 20.0))                    ;
  9.     mess '=========================================================' ;
  10.     mess '===RVSAT: Interpolazione fuori limiti: 0 < T[C] < 20===' ;
  11.     mess '=========================================================' ;
  12.    finsi                                                             ;
  13.    a0   =  4.86256716e-03                                            ;
  14.    a1   =  2.98409029e-04                                            ;
  15.    a2   =  2.10179146e-05                                            ;
  16.    a3   = -7.96684413e-07                                            ;
  17.    a4   =  2.77512067e-08                                            ;
  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.  
  37.  
  38.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales