Procedures | Cast3M

   1 : * CALCP     PROCEDUR  MAGN      99/12/09    21:15:07     3728           
   2 : DEBPROC CALCP ;
   3 : 'ARGU' T/'CHPOINT'
   4 :        TF/'FLOTTANT'
   5 :        TL/'LISTREEL' ;
   6 : *-------------------------------------------------*
   7 : * calcul de la capacite calorifique a pression    *
   8 : * constante pour les incondensables en fonction   *
   9 : * de la temperature  (J/kg/K)                     *
  10 : *                                                 *
  11 : * Entree T en Celsius                             *
  12 : * Sortie Cp en J/kg/K                             *
  13 : *                                                 *
  14 : * PS : le resultat est du meme type que l'entree  *
  15 : *-------------------------------------------------*
  16 : A0 = 14.346 ;
  17 : A1 = 8.16D-4 ;
  18 : A2 = 4.35D-7 ;
  19 : B0 = 5.193 ;
  20 : B1 = 0.0 ;
  21 : B2 = 0.0 ;
  22 : C0 = 0.917 ;
  23 : C1 = 1.404D-4 ;
  24 : C2 = 3.0D-8 ;
  25 : D0 = 1.0309 ;
  26 : D1 = 2.044D-4 ;
  27 : D2 = -3.3D-8 ;
  28 : E0 = 0.836  ;
  29 : E1 = 7.432D-4 ;
  30 : E2 = -2.511D-7 ;
  31 : F0 = 1.0302 ;
  32 : F1 = 2.396D-4 ;
  33 : F2 = -4.950D-8 ;
  34 : 
  35 : 'SI' ( 'EXIST' T ) ;
  36 :     Cph2  = ( ((A2 * T) * T ) + ( T * A1 ) + A0 ) * 1.D3        ;
  37 :     Cphe  = ( ((B2 * T) * T ) + ( T * B1 ) + B0 ) * 1.D3        ;
  38 :     Cpo2  = ( ((C2 * T) * T ) + ( T * C1 ) + C0 ) * 1.D3        ;
  39 :     Cpn2  = ( ((D2 * T) * T ) + ( T * D1 ) + D0 ) * 1.D3        ;
  40 :     Cpco2 = ( ((E2 * T) * T ) + ( T * E1 ) + E0 ) * 1.D3        ;
  41 :     Cpco  = ( ((F2 * T) * T ) + ( T * F1 ) + F0 ) * 1.D3        ;
  42 :     Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  43 :     'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  44 :     'QUITTER' CALCP ;
  45 : 'FINSI' ;
  46 : 'SI' ( 'EXIST' TF) ;
  47 :     Cph2  = ( ((A2 * TF) * TF ) + ( TF * A1 ) + A0 ) * 1.D3        ;
  48 :     Cphe  = ( ((B2 * TF) * TF ) + ( TF * B1 ) + B0 ) * 1.D3        ;
  49 :     Cpo2  = ( ((C2 * TF) * TF ) + ( TF * C1 ) + C0 ) * 1.D3        ;
  50 :     Cpn2  = ( ((D2 * TF) * TF ) + ( TF * D1 ) + D0 ) * 1.D3        ;
  51 :     Cpco2 = ( ((E2 * TF) * TF ) + ( TF * E1 ) + E0 ) * 1.D3        ;
  52 :     Cpco  = ( ((F2 * TF) * TF ) + ( TF * F1 ) + F0 ) * 1.D3        ;
  53 :     Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  54 :     'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  55 :     'QUITTER' CALCP ;
  56 : 'FINSI' ;
  57 : 'SI' ( 'EXIST' TL) ;
  58 :     T00 = TL / TL ;
  59 :     A00 = A0 * T00 ;
  60 :     A10 = A1 * T00 ;
  61 :     A20 = A2 * T00 ;
  62 :     B00 = B0 * T00 ;
  63 :     B10 = B1 * T00 ;
  64 :     B20 = B2 * T00 ;
  65 :     C00 = C0 * T00 ;
  66 :     C10 = C1 * T00 ;
  67 :     C20 = C2 * T00 ;
  68 :     D00 = D0 * T00 ;
  69 :     D10 = D1 * T00 ;
  70 :     D20 = D2 * T00 ;
  71 :     E00 = E0 * T00 ;
  72 :     E10 = E1 * T00 ;
  73 :     E20 = E2 * T00 ;
  74 :     F00 = F0 * T00 ;
  75 :     F10 = F1 * T00 ;
  76 :     F20 = F2 * T00 ;
  77 : 
  78 :     Cph2  = ( ((A20* TL) * TL ) + ( TL * A10) + A00) * 1.D3        ;
  79 :     Cphe  = ( ((B20* TL) * TL ) + ( TL * B10) + B00) * 1.D3        ;
  80 :     Cpo2  = ( ((C20* TL) * TL ) + ( TL * C10) + C00) * 1.D3        ;
  81 :     Cpn2  = ( ((D20* TL) * TL ) + ( TL * D10) + D00) * 1.D3        ;
  82 :     Cpco2 = ( ((E20* TL) * TL ) + ( TL * E10) + E00) * 1.D3        ;
  83 :     Cpco  = ( ((F20* TL) * TL ) + ( TL * F10) + F00) * 1.D3        ;
  84 :     Cpair = (0.8*Cpn2) + (0.2*Cpo2) ;
  85 :     'RESPROC' Cph2 Cphe Cpo2 Cpn2 Cpco2 Cpco Cpair ;
  86 :     'QUITTER' CALCP ;
  87 : 
  88 : 'FINSI' ;
  89 : 'ERRE' 5 ;
  90 : FINPROC                ;
  91 : 
  92 : 
  93 :