Download callm.procedur

Back to the list

   1 : * CALLM     PROCEDUR  MAGN      08/03/13    21:15:01     6072           
   2 : DEBPROC CALLM ;
   3 : 'ARGU' T/'CHPOINT'
   4 :        TF/'FLOTTANT'
   5 :        TL/'LISTREEL' ;
   6 : 
   7 : *-------------------------------------------------*
   8 : * calcul de la conductivité thermique des gaz     *
   9 : * suivants: H2, He , O2, N2, Co2, Co, Vapeur H2O  *
  10 : * (W/m/K)                                         *
  11 : *                                                 *
  12 : * Entree T en K                                   *
  13 : * Sortie Lambda en  (W/m/K)                       *
  14 : *                                                 *
  15 : * PS : le resultat est du meme type que l'entree  *
  16 : *-------------------------------------------------*
  17 : * hydrogene
  18 :       CST01 = 1.65D-3;
  19 :       CST11 = -40.13D0;
  20 :       CST21 = 1.800D0;
  21 : * helium
  22 :       CST02 = 5.464D-3;
  23 :       CST12 = 76.6D0 ;
  24 :       CST22 = 1.616D0 ;
  25 : * oxygene
  26 :       CST03 = 0.429D-3;
  27 :       CST13 = 41.27D0 ;
  28 :       CST23 = 1.746D0  ;
  29 : * azote
  30 :       CST04 = 0.0417D-3;
  31 :       CST14 = -133.D0 ;
  32 :       CST24 = 2.042D0  ; 
  33 : * co2
  34 :       CST05 = 25.438D-3;
  35 :       CST15 = 2432.95D0 ; 
  36 :       CST25 = 1.319D0   ;
  37 : * co
  38 :       CST06 = 1.931D-3;
  39 :       CST16 = 218.17D0 ;  
  40 :       CST26 = 1.541D0   ;
  41 : * vapeur
  42 :       CST07 =  -84.442D-3;
  43 :       CST17 =  -5827.6D0 ;
  44 :       CST27 =  1.251D0 ;
  45 : *
  46 : 'SI' ( 'EXIST' T ) ;
  47 :       RLAH2 = ((T ** CST21) * CST01) / (T + CST11);
  48 :       RLAHE = (((T ** CST22) * CST02) / (T + CST12));
  49 :       RLAO2 = (((T ** CST23) * CST03) / (T + CST13));
  50 :       RLAN2 = (((T ** CST24) * CST04) / (T + CST14));
  51 :       RLACO2 = (((T ** CST25) * CST05) / (T + CST15));
  52 :       RLACO = (((T ** CST26) * CST06) / (T + CST16));
  53 :       RLAVAP = (((T ** CST27) * CST07) / (T + CST17));
  54 :       RLAAIR = (0.8D0 * RLAN2) + (0.2D0 * RLAO2);
  55 :     'RESPROC' RLAH2 RLAHE RLAO2 RLAN2 RLACO2 RLACO RLAVAP RLAAIR ;
  56 :     'QUITTER' CALLM ;
  57 : 'FINSI' ;
  58 : 'SI' ( 'EXIST' TF) ;
  59 :       RLAH2 = ((TF ** CST21) * CST01) / (TF + CST11);
  60 :       RLAHE = (((TF ** CST22) * CST02) / (TF + CST12));
  61 :       RLAO2 = (((TF ** CST23) * CST03) / (TF + CST13));
  62 :       RLAN2 = (((TF ** CST24) * CST04) / (TF + CST14));
  63 :       RLACO2 = (((TF ** CST25) * CST05) / (TF + CST15));
  64 :       RLACO = (((TF ** CST26) * CST06) / (TF + CST16));
  65 :       RLAVAP = (((TF ** CST27) * CST07) / (TF + CST17));
  66 :       RLAAIR = (0.8D0 * RLAN2) + (0.2D0 * RLAO2);
  67 :     'RESPROC' RLAH2 RLAHE RLAO2 RLAN2 RLACO2 RLACO RLAVAP RLAAIR ;
  68 :     'QUITTER' CALLM ;
  69 : 'FINSI' ;
  70 : 'SI' ( 'EXIST' TL) ;
  71 :       T00 = TL / TL ;
  72 :       DST01 = CST01 * T00 ;
  73 :       DST11 = CST11 * T00 ;
  74 :       DST21 = CST21 * T00 ;
  75 :       DST02 = CST02 * T00 ;
  76 :       DST12 = CST12 * T00 ;
  77 :       DST22 = CST22 * T00 ;
  78 :       DST03 = CST03 * T00 ;
  79 :       DST13 = CST13 * T00 ;
  80 :       DST23 = CST23 * T00 ;
  81 :       DST04 = CST04 * T00 ;
  82 :       DST14 = CST14 * T00 ;
  83 :       DST24 = CST24 * T00 ;
  84 :       DST05 = CST05 * T00 ;
  85 :       DST15 = CST15 * T00 ;
  86 :       DST25 = CST25 * T00 ;
  87 :       DST06 = CST06 * T00 ;
  88 :       DST16 = CST16 * T00 ;
  89 :       DST26 = CST26 * T00 ;
  90 :       DST07 = CST07 * T00 ;
  91 :       DST17 = CST17 * T00 ;
  92 :       DST27 = CST27 * T00 ;
  93 :       RLAH2 = ((TL ** DST21) * DST01) / (TL + DST11);
  94 :       RLAHE = (((TL ** DST22) * DST02) / (TL + DST12));
  95 :       RLAO2 = (((TL ** DST23) * DST03) / (TL + DST13));
  96 :       RLAN2 = (((TL ** DST24) * DST04) / (TL + DST14));
  97 :       RLACO2 = (((TL ** DST25) * DST05) / (TL + DST15));
  98 :       RLACO = (((TL ** DST26) * DST06) / (TL + DST16));
  99 :       RLAVAP = (((TL ** DST27) * DST07) / (TL + DST17));
 100 :       RLAAIR = (0.8D0 * RLAN2) + (0.2D0 * RLAO2);
 101 :     'RESPROC' RLAH2 RLAHE RLAO2 RLAN2 RLACO2 RLACO RLAVAP RLAAIR ;
 102 :     'QUITTER' CALLM ;
 103 : 'FINSI' ;
 104 : 'ERRE' 5 ;
 105 : FINPROC                ;
 106 :  

© Cast3M 2003 - All rights reserved.
Disclaimer