Télécharger callm.procedur

Retour à la liste

Numérotation des lignes :

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

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