Télécharger la plaquette Cast3M


fr en

Télécharger rvsat.procedur

Retour à la liste

Numérotation des lignes : 

  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.  
  37.  

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