$$$$ CALLM * CALLM PROCEDUR MAGN 08/03/13 21:15:01 6072 DEBPROC CALLM ; 'ARGU' T/'CHPOINT' TF/'FLOTTANT' TL/'LISTREEL' ; *-------------------------------------------------* * calcul de la conductivité thermique des gaz * * suivants: H2, He , O2, N2, Co2, Co, Vapeur H2O * * (W/m/K) * * * * Entree T en K * * Sortie Lambda en (W/m/K) * * * * PS : le resultat est du meme type que l'entree * *-------------------------------------------------* * hydrogene CST01 = 1.65D-3; CST11 = -40.13D0; CST21 = 1.800D0; * helium CST02 = 5.464D-3; CST12 = 76.6D0 ; CST22 = 1.616D0 ; * oxygene CST03 = 0.429D-3; CST13 = 41.27D0 ; CST23 = 1.746D0 ; * azote CST04 = 0.0417D-3; CST14 = -133.D0 ; CST24 = 2.042D0 ; * co2 CST05 = 25.438D-3; CST15 = 2432.95D0 ; CST25 = 1.319D0 ; * co CST06 = 1.931D-3; CST16 = 218.17D0 ; CST26 = 1.541D0 ; * vapeur CST07 = -84.442D-3; CST17 = -5827.6D0 ; CST27 = 1.251D0 ; * 'SI' ( 'EXIST' T ) ; RLAH2 = ((T ** CST21) * CST01) / (T + CST11); RLAHE = (((T ** CST22) * CST02) / (T + CST12)); RLAO2 = (((T ** CST23) * CST03) / (T + CST13)); RLAN2 = (((T ** CST24) * CST04) / (T + CST14)); RLACO2 = (((T ** CST25) * CST05) / (T + CST15)); RLACO = (((T ** CST26) * CST06) / (T + CST16)); RLAVAP = (((T ** CST27) * CST07) / (T + CST17)); RLAAIR = (0.8D0 * RLAN2) + (0.2D0 * RLAO2); 'RESPROC' RLAH2 RLAHE RLAO2 RLAN2 RLACO2 RLACO RLAVAP RLAAIR ; 'QUITTER' CALLM ; 'FINSI' ; 'SI' ( 'EXIST' TF) ; RLAH2 = ((TF ** CST21) * CST01) / (TF + CST11); RLAHE = (((TF ** CST22) * CST02) / (TF + CST12)); RLAO2 = (((TF ** CST23) * CST03) / (TF + CST13)); RLAN2 = (((TF ** CST24) * CST04) / (TF + CST14)); RLACO2 = (((TF ** CST25) * CST05) / (TF + CST15)); RLACO = (((TF ** CST26) * CST06) / (TF + CST16)); RLAVAP = (((TF ** CST27) * CST07) / (TF + CST17)); RLAAIR = (0.8D0 * RLAN2) + (0.2D0 * RLAO2); 'RESPROC' RLAH2 RLAHE RLAO2 RLAN2 RLACO2 RLACO RLAVAP RLAAIR ; 'QUITTER' CALLM ; 'FINSI' ; 'SI' ( 'EXIST' TL) ; T00 = TL / TL ; DST01 = CST01 * T00 ; DST11 = CST11 * T00 ; DST21 = CST21 * T00 ; DST02 = CST02 * T00 ; DST12 = CST12 * T00 ; DST22 = CST22 * T00 ; DST03 = CST03 * T00 ; DST13 = CST13 * T00 ; DST23 = CST23 * T00 ; DST04 = CST04 * T00 ; DST14 = CST14 * T00 ; DST24 = CST24 * T00 ; DST05 = CST05 * T00 ; DST15 = CST15 * T00 ; DST25 = CST25 * T00 ; DST06 = CST06 * T00 ; DST16 = CST16 * T00 ; DST26 = CST26 * T00 ; DST07 = CST07 * T00 ; DST17 = CST17 * T00 ; DST27 = CST27 * T00 ; RLAH2 = ((TL ** DST21) * DST01) / (TL + DST11); RLAHE = (((TL ** DST22) * DST02) / (TL + DST12)); RLAO2 = (((TL ** DST23) * DST03) / (TL + DST13)); RLAN2 = (((TL ** DST24) * DST04) / (TL + DST14)); RLACO2 = (((TL ** DST25) * DST05) / (TL + DST15)); RLACO = (((TL ** DST26) * DST06) / (TL + DST16)); RLAVAP = (((TL ** DST27) * DST07) / (TL + DST17)); RLAAIR = (0.8D0 * RLAN2) + (0.2D0 * RLAO2); 'RESPROC' RLAH2 RLAHE RLAO2 RLAN2 RLACO2 RLACO RLAVAP RLAAIR ; 'QUITTER' CALLM ; 'FINSI' ; 'ERRE' 5 ; FINPROC ;