Télécharger flam.eso

Retour à la liste

Numérotation des lignes :

  1. C FLAM SOURCE TTMF3 12/07/05 21:15:44 7425
  2. SUBROUTINE FLAM
  3. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  4. CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
  5.  
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8 (A-H,O-Z)
  8. C REAL*8 MH2,MO2,MN2,MH2O
  9. C DATA Rg,MH2,MO2,MN2,MH2O/8.313D0,2.0D-3,32.0D-3,28.0D-3,18.0D-3/
  10. C DATA RRH2,RRO2,RRN2,RRH2O/4156.5D0,259.8D0,296.9D0,461.8D0/
  11. -INC CCOPTIO
  12. -INC CCREEL
  13. C-
  14. -INC SMCHPOI
  15. -INC SMLREEL
  16. C
  17. POINTEUR MYH2.MCHPOI,MYO2.MCHPOI,MYN2.MCHPOI,MYH2O.MCHPOI,
  18. 1 MTEMP.MCHPOI,MQ.MCHPOI,MRHO.MCHPOI,MCV.MCHPOI,MR.MCHPOI,
  19. 2 MZK.MCHPOI,MZE.MCHPOI,MYH2u.MCHPOI,MTSEUI.MCHPOI,
  20. 3 MZH2.MCHPOI,MZO2.MCHPOI,MZN2.MCHPO,MZH2O.MCHPOI
  21. POINTEUR NYH2.MPOVAL,NYO2.MPOVAL,NYN2.MPOVAL,NYH2O.MPOVAL,
  22. 1 NTEMP.MPOVAL,NQ.MPOVAL,NRHO.MPOVAL,NCV.MPOVAL,NR.MPOVAL,
  23. 2 NZK.MPOVAL,NZE.MPOVAL,NYH2u.MPOVAL,NTSEUI.MPOVAL,
  24. 3 NZH2.MPOVAL,NZO2.MPOVAL,NZN2.MPOVAL,NZH2O.MPOVAL
  25. C
  26. PARAMETER (NBOPT=8)
  27. PARAMETER (NBCV=2)
  28. PARAMETER (NBINT=2)
  29. PARAMETER (NBMOD=2)
  30. CHARACTER*8 MOTOP,TOPTIO(NBOPT),COPTIO(NBCV),TYPE,MODEL
  31. CHARACTER*8 IOPTIO(NBINT),MODREC(NBMOD)
  32. DATA TOPTIO/'LAMINAIR','EBU ','ARRHENIU','RECOMBIN',
  33. 1 'RECOMB0D','HEAVYSID','CREBCOM','CREBCOM2'/
  34. DATA COPTIO/'LINEAIRE','QUADRATI'/
  35. DATA MODREC/'SIEMENS ','HEATER '/
  36. C
  37. C==================================
  38. C Lecture de l'option de combustion
  39. C==================================
  40. C
  41. CALL LIRCHA(MOTOP,1,IRET)
  42. CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT)
  43. C
  44. C======================================================
  45. C Branchements sur les différents modèles de combustion
  46. C======================================================
  47. C
  48. C=================================
  49. C Cas de la déflagration laminaire
  50. C=================================
  51. C
  52. IF (IP .EQ. 1) THEN
  53. C
  54. C Lecture des paramètres
  55. C
  56. CALL LIROBJ('CHPOINT ',MRHO,1,IRET)
  57. IF (IRET.EQ.0) RETURN
  58. CALL LIROBJ('CHPOINT ',MCV,1,IRET)
  59. IF (IRET.EQ.0) RETURN
  60. C
  61. C Lecture de l'option sur la forme des Cv
  62. C
  63. CALL LIRCHA(MOTOP,1,IRET)
  64. CALL OPTLI(IP,COPTIO,MOTOP,NBCV)
  65. IF (IP .EQ. 0) THEN
  66. CALL ERREUR(251)
  67. RETURN
  68. ENDIF
  69. IF ( MOTOP .EQ. 'LINEAIRE') THEN
  70. NCVOPT = 1
  71. ELSEIF( MOTOP .EQ. 'QUADRATI') THEN
  72. NCVOPT = 2
  73. ENDIF
  74. C
  75. CALL LIROBJ('CHPOINT ',MR,1,IRET)
  76. IF (IRET.EQ.0) RETURN
  77. CALL LIROBJ('CHPOINT ',MYH2,1,IRET)
  78. IF (IRET.EQ.0) RETURN
  79. CALL LIROBJ('CHPOINT ',MYO2,1,IRET)
  80. IF (IRET.EQ.0) RETURN
  81. CALL LIROBJ('CHPOINT ',MYN2,1,IRET)
  82. IF (IRET.EQ.0) RETURN
  83. CALL LIROBJ('CHPOINT ',MYH2O,1,IRET)
  84. IF (IRET.EQ.0) RETURN
  85. CALL LIROBJ('CHPOINT ',MTEMP,1,IRET)
  86. IF (IRET.EQ.0) RETURN
  87. CALL LIRREE(Dt,1,IRET)
  88. IF (IRET.EQ.0) RETURN
  89. C
  90. C Activation champoints et vérification des pointeurs
  91. C
  92. CALL LICHT(MRHO,NRHO,TYPE,IGEO)
  93. CALL LICHT(MCV,NCV,TYPE,IGEO0)
  94. IF (IGEO.NE.IGEO0) THEN
  95. WRITE(IOIMP,100)
  96. RETURN
  97. ENDIF
  98. CALL LICHT(MR,NR,TYPE,IGEO0)
  99. IF (IGEO.NE.IGEO0) THEN
  100. WRITE(IOIMP,100)
  101. RETURN
  102. ENDIF
  103. CALL LICHT(MYH2,NYH2,TYPE,IGEO0)
  104. IF (IGEO.NE.IGEO0) THEN
  105. WRITE(IOIMP,100)
  106. RETURN
  107. ENDIF
  108. CALL LICHT(MYO2,NYO2,TYPE,IGEO0)
  109. IF (IGEO.NE.IGEO0) THEN
  110. WRITE(IOIMP,100)
  111. RETURN
  112. ENDIF
  113. CALL LICHT(MYN2,NYN2,TYPE,IGEO0)
  114. IF (IGEO.NE.IGEO0) THEN
  115. WRITE(IOIMP,100)
  116. RETURN
  117. ENDIF
  118. CALL LICHT(MYH2O,NYH2O,TYPE,IGEO0)
  119. IF (IGEO.NE.IGEO0) THEN
  120. WRITE(IOIMP,100)
  121. RETURN
  122. ENDIF
  123. CALL LICHT(MTEMP,NTEMP,TYPE,IGEO0)
  124. IF (IGEO.NE.IGEO0) THEN
  125. WRITE(IOIMP,100)
  126. RETURN
  127. ENDIF
  128. C
  129. C Création et activation des chpoints résultats
  130. C
  131. NC=1
  132. CALL CRCHPT(TYPE,IGEO,NC,MZH2)
  133. CALL LICHT(MZH2,NZH2,TYPE,IGEO)
  134. CALL CRCHPT(TYPE,IGEO,NC,MZO2)
  135. CALL LICHT(MZO2,NZO2,TYPE,IGEO)
  136. CALL CRCHPT(TYPE,IGEO,NC,MZN2)
  137. CALL LICHT(MZN2,NZN2,TYPE,IGEO)
  138. CALL CRCHPT(TYPE,IGEO,NC,MZH2O)
  139. CALL LICHT(MZH2O,NZH2O,TYPE,IGEO)
  140. CALL CRCHPT(TYPE,IGEO,NC,MQ)
  141. CALL LICHT(MQ,NQ,TYPE,IGEO)
  142. C
  143. C Nombre de Points du Maillage
  144. C
  145. NBPTS = NRHO.VPOCHA(/1)
  146. C
  147. C Calcul des nouvelles concentrations et de l'énergie libérée
  148. C
  149. CALL FLAM1(NBPTS,NRHO.VPOCHA,NCV.VPOCHA,NR.VPOCHA,NYH2.VPOCHA,
  150. 1 NYO2.VPOCHA,NYN2.VPOCHA,NYH2O.VPOCHA,NTEMP.VPOCHA,
  151. 2 NZH2.VPOCHA,NZO2.VPOCHA,NZN2.VPOCHA,NZH2O.VPOCHA,
  152. 3 NQ.VPOCHA,Dt,NCVOPT)
  153. C
  154. C Ecriture des chpoints résultats
  155. C
  156. CALL ECROBJ('CHPOINT ',MQ)
  157. CALL ECROBJ('CHPOINT ',MZH2O)
  158. CALL ECROBJ('CHPOINT ',MZN2)
  159. CALL ECROBJ('CHPOINT ',MZO2)
  160. CALL ECROBJ('CHPOINT ',MZH2)
  161. C
  162. C Désactivation des segments
  163. C
  164. SEGDES NZH2,NZO2,NZN2,NZH2O,NQ,NRHO,NCV,NR,NYH2
  165. SEGDES NYO2,NYN2,NYH2O,NTEMP
  166. C
  167. C
  168. C==================================
  169. C Cas de la déflagration turbulente
  170. C==================================
  171. C
  172. ELSEIF( IP .EQ. 2) THEN
  173. C
  174. C Lecture des paramètres
  175. C
  176. CALL LIROBJ('CHPOINT ',MRHO,1,IRET)
  177. IF (IRET.EQ.0) RETURN
  178. CALL LIROBJ('CHPOINT ',MCV,1,IRET)
  179. IF (IRET.EQ.0) RETURN
  180. C
  181. C Lecture de l'option sur la forme des Cv
  182. C
  183. CALL LIRCHA(MOTOP,1,IRET)
  184. CALL OPTLI(IP,COPTIO,MOTOP,NBCV)
  185. IF(IP .EQ. 0) THEN
  186. CALL ERREUR(251)
  187. RETURN
  188. ENDIF
  189. IF ( MOTOP .EQ. 'LINEAIRE') THEN
  190. NCVOPT = 1
  191. ELSEIF( MOTOP .EQ. 'QUADRATI') THEN
  192. NCVOPT = 2
  193. ENDIF
  194. C
  195. CALL LIROBJ('CHPOINT ',MR,1,IRET)
  196. IF (IRET.EQ.0) RETURN
  197. CALL LIROBJ('CHPOINT ',MYH2,1,IRET)
  198. IF (IRET.EQ.0) RETURN
  199. CALL LIROBJ('CHPOINT ',MYO2,1,IRET)
  200. IF (IRET.EQ.0) RETURN
  201. CALL LIROBJ('CHPOINT ',MYN2,1,IRET)
  202. IF (IRET.EQ.0) RETURN
  203. CALL LIROBJ('CHPOINT ',MYH2O,1,IRET)
  204. IF (IRET.EQ.0) RETURN
  205. CALL LIROBJ('CHPOINT ',MTEMP,1,IRET)
  206. IF (IRET.EQ.0) RETURN
  207. CALL LIRREE(Dt,1,IRET)
  208. IF (IRET.EQ.0) RETURN
  209.  
  210. C paramètres spécifiques à l'option choisie
  211.  
  212. CALL LIROBJ('CHPOINT ',MZK,1,IRET)
  213. IF (IRET.EQ.0) RETURN
  214. CALL LIROBJ('CHPOINT ',MZE,1,IRET)
  215. IF (IRET.EQ.0) RETURN
  216. CALL LIROBJ('CHPOINT ',MYH2u,1,IRET)
  217. IF (IRET.EQ.0) RETURN
  218. CALL LIRREE(C_EBU,1,IRET)
  219. IF (IRET.EQ.0) RETURN
  220. CALL LIRREE(CMIN,0,IRET)
  221. IF (IRET.EQ.0) CMIN=1.D-4
  222. C
  223. C Activation champoints et vérification des pointeurs
  224. C
  225. CALL LICHT(MRHO,NRHO,TYPE,IGEO)
  226. CALL LICHT(MCV,NCV,TYPE,IGEO0)
  227. IF (IGEO.NE.IGEO0) THEN
  228. WRITE(IOIMP,100)
  229. RETURN
  230. ENDIF
  231. CALL LICHT(MR,NR,TYPE,IGEO0)
  232. IF (IGEO.NE.IGEO0) THEN
  233. WRITE(IOIMP,100)
  234. RETURN
  235. ENDIF
  236. CALL LICHT(MYH2,NYH2,TYPE,IGEO0)
  237. IF (IGEO.NE.IGEO0) THEN
  238. WRITE(IOIMP,100)
  239. RETURN
  240. ENDIF
  241. CALL LICHT(MYO2,NYO2,TYPE,IGEO0)
  242. IF (IGEO.NE.IGEO0) THEN
  243. WRITE(IOIMP,100)
  244. RETURN
  245. ENDIF
  246. CALL LICHT(MYN2,NYN2,TYPE,IGEO0)
  247. IF (IGEO.NE.IGEO0) THEN
  248. WRITE(IOIMP,100)
  249. RETURN
  250. ENDIF
  251. CALL LICHT(MYH2O,NYH2O,TYPE,IGEO0)
  252. IF (IGEO.NE.IGEO0) THEN
  253. WRITE(IOIMP,100)
  254. RETURN
  255. ENDIF
  256. CALL LICHT(MTEMP,NTEMP,TYPE,IGEO0)
  257. IF (IGEO.NE.IGEO0) THEN
  258. WRITE(IOIMP,100)
  259. RETURN
  260. ENDIF
  261. CALL LICHT(MZK,NZK,TYPE,IGEO0)
  262. IF (IGEO.NE.IGEO0) THEN
  263. WRITE(IOIMP,100)
  264. RETURN
  265. ENDIF
  266. CALL LICHT(MZE,NZE,TYPE,IGEO0)
  267. IF (IGEO.NE.IGEO0) THEN
  268. WRITE(IOIMP,100)
  269. RETURN
  270. ENDIF
  271. CALL LICHT(MYH2u,NYH2u,TYPE,IGEO0)
  272. IF (IGEO.NE.IGEO0) THEN
  273. WRITE(IOIMP,100)
  274. RETURN
  275. ENDIF
  276. C
  277. C Création et activation des chpoints résultats
  278. C
  279. NC=1
  280. CALL CRCHPT(TYPE,IGEO,NC,MZH2)
  281. CALL LICHT(MZH2,NZH2,TYPE,IGEO)
  282. CALL CRCHPT(TYPE,IGEO,NC,MZO2)
  283. CALL LICHT(MZO2,NZO2,TYPE,IGEO)
  284. CALL CRCHPT(TYPE,IGEO,NC,MZN2)
  285. CALL LICHT(MZN2,NZN2,TYPE,IGEO)
  286. CALL CRCHPT(TYPE,IGEO,NC,MZH2O)
  287. CALL LICHT(MZH2O,NZH2O,TYPE,IGEO)
  288. CALL CRCHPT(TYPE,IGEO,NC,MQ)
  289. CALL LICHT(MQ,NQ,TYPE,IGEO)
  290. C
  291. C Nombre de Points
  292. C
  293. NBPTS = NRHO.VPOCHA(/1)
  294. C
  295. CALL FLAM2(NBPTS,NRHO.VPOCHA,NCV.VPOCHA,NR.VPOCHA,NYH2.VPOCHA,
  296. 1 NYO2.VPOCHA,NYN2.VPOCHA,NYH2O.VPOCHA,NTEMP.VPOCHA,
  297. 2 NZK.VPOCHA,NZE.VPOCHA,NYH2u.VPOCHA,
  298. 3 NZH2.VPOCHA,NZO2.VPOCHA,NZN2.VPOCHA,NZH2O.VPOCHA,
  299. 4 NQ.VPOCHA,C_EBU,Dt,CMIN,NCVOPT)
  300. C
  301. C Ecriture des chpoints résultats
  302. C
  303. CALL ECROBJ('CHPOINT ',MQ)
  304. CALL ECROBJ('CHPOINT ',MZH2O)
  305. CALL ECROBJ('CHPOINT ',MZN2)
  306. CALL ECROBJ('CHPOINT ',MZO2)
  307. CALL ECROBJ('CHPOINT ',MZH2)
  308. C
  309. C Désactivation des segments
  310. C
  311. SEGDES NZH2,NZO2,NZN2,NZH2O,NQ,NRHO,NCV,NR,NYH2
  312. SEGDES NYO2,NYN2,NYH2O,NTEMP,NZK,NZE,NYH2u
  313. C
  314. C===================================
  315. C Cas de la loi d'Arrhénius générale
  316. C===================================
  317. C
  318. ELSEIF( IP .EQ. 3) THEN
  319. CALL FLAMB()
  320. C
  321. C===================================
  322. C Maquette CREBCOM
  323. C===================================
  324. C
  325. ELSEIF( IP .EQ. 7) THEN
  326. CALL FLAMCR()
  327. IF(IERR.NE.0)RETURN
  328. C
  329. C======================================
  330. C Maquette CREBCOM, modele non-homogene
  331. C======================================
  332. C
  333. ELSEIF( IP .EQ. 8) THEN
  334. CALL FLACR2()
  335. IF(IERR.NE.0)RETURN
  336. C
  337. C============================================
  338. C Combustion catalytique : modèle recombineur
  339. C============================================
  340. C
  341. ELSEIF((IP .EQ. 4) .OR. (IP .EQ. 5)) THEN
  342. C
  343. C Lecture des paramètres
  344. C YH2MOY : flottant : Fraction massique moyenne de H2 à l'entrée du PAR
  345. C YO2MOY : flottant : Fraction massique moyenne de O2 à l'entrée du PAR
  346. C YN2MOY : flottant : Fraction massique moyenne de N2 à l'entrée du PAR
  347. C YH2OMOY : flottant : Fraction massique moyenne de H2O à l'entrée du PAR
  348. C YHEMOY : flottant : Fraction massique moyenne de HE à l'entrée du PAR
  349. C YCO2MOY : flottant : Fraction massique moyenne de CO2 à l'entrée du PAR
  350. C YCOMOY : flottant : Fraction massique moyenne de CO à l'entrée du PAR
  351. C TMOY : flottant : Température moyenne du gaz à l'entrée du PAR (K)
  352. C TPLA : flottant : Température moyenne des plaques du PAR (K)
  353. C PRESS : flottant : Pression moyenne autour du PAR (Pa)
  354. C DT : flottant : Pas de temps (s)
  355. C KALP : entier : Flag pour correction du coefficient d'échange
  356. C comme dans RALOC (voir fla009.eso)
  357. CALL LIRREE(YH2MOY,1,IRET)
  358. IF (IRET.EQ.0) RETURN
  359. CALL LIRREE(YO2MOY,1,IRET)
  360. IF (IRET.EQ.0) RETURN
  361. CALL LIRREE(YN2MOY,1,IRET)
  362. IF (IRET.EQ.0) RETURN
  363. CALL LIRREE(YH2OMOY,1,IRET)
  364. IF (IRET.EQ.0) RETURN
  365. CALL LIRREE(YHEMOY,1,IRET)
  366. IF (IRET.EQ.0) RETURN
  367. CALL LIRREE(YCO2MOY,1,IRET)
  368. IF (IRET.EQ.0) RETURN
  369. CALL LIRREE(YCOMOY,1,IRET)
  370. IF (IRET.EQ.0) RETURN
  371. CALL LIRREE(TMOY,1,IRET)
  372. IF (IRET.EQ.0) RETURN
  373. CALL LIRREE(TPLA,1,IRET)
  374. IF (IRET.EQ.0) RETURN
  375. CALL LIRREE(PRESS,1,IRET)
  376. IF (IRET.EQ.0) RETURN
  377. CALL LIRREE(Dt,1,IRET)
  378. IF (IRET.EQ.0) RETURN
  379. CALL LIRENT(KALP,1,IRET)
  380. IF (IRET.EQ.0) RETURN
  381. C Lecture du modèle de PAR à choisir parmi 'SIEMENS ' ou 'HEATER '
  382. C MODEL : CHAR*8 : Nom du modèle de PAR
  383. CALL LIRCHA(MODEL,1,IRET)
  384. IF (IRET.EQ.0) RETURN
  385. CALL OPTLI(IQ,MODREC,MODEL,NBMOD)
  386. IF (IQ.EQ.1) THEN
  387. C Lecture des paramètres A et B du modele recombineur 'SIEMENS '
  388. C (loi lineaire en pression)
  389. CALL LIRREE(A1,1,IRET)
  390. IF (IRET.EQ.0) RETURN
  391. CALL LIRREE(B1,1,IRET)
  392. IF (IRET.EQ.0) RETURN
  393. CALL LIRREE(XH2ON,1,IRET)
  394. IF (IRET.EQ.0) RETURN
  395. CALL LIRREE(XH2OFF,1,IRET)
  396. IF (IRET.EQ.0) RETURN
  397. CALL LIRREE(BETA,1,IRET)
  398. IF (IRET.EQ.0) RETURN
  399. JG = 5
  400. SEGINI MLREEL
  401. PROG(1) = A1
  402. PROG(2) = B1
  403. PROG(3) = XH2ON
  404. PROG(4) = XH2OFF
  405. PROG(5) = BETA
  406. C
  407. C Calcul des conditions de sortie, de la temperature des plaques et du
  408. C débit dans le PAR (arret si débit ou température des plaques négatif)
  409. CALL FLAMR(YH2MOY,YO2MOY,YN2MOY,YH2OMOY,YHEMOY,
  410. & YCO2MOY,YCOMOY,PRESS,TMOY,
  411. & TS,TPLA,DEB,YH2SOR,YO2SOR,YN2SOR,YH2OSOR,YHESOR,
  412. & YCO2SOR,YCOSOR,XH2ONF,EFF,Dt,KALP,IQ,MLREEL,IHI)
  413. ELSEIF (IQ.EQ.2) THEN
  414. C Lecture de la puissance du recombineur thermique (modèle 'HEATER ')
  415. CALL LIRREE(QHEAT,1,IRET)
  416. IF (IRET.EQ.0) RETURN
  417. CALL LIRREE(FUITE,1,IRET)
  418. IF (IRET.EQ.0) RETURN
  419. CALL LIRREE(HRAPP,1,IRET)
  420. IF (IRET.EQ.0) RETURN
  421. CALL LIRREE(CRAPP,1,IRET)
  422. IF (IRET.EQ.0) RETURN
  423. JG = 4
  424. SEGINI MLREEL
  425. PROG(1) = QHEAT
  426. PROG(2) = 1.D0 - fuite
  427. PROG(3) = hrapp
  428. PROG(4) = crapp
  429. C
  430. C Calcul des conditions de sortie, de la temperature des plaques et du
  431. C débit dans le HEATER (arret si débit ou température des plaques négatif)
  432. CALL FLAMH(YH2MOY,YO2MOY,YN2MOY,YH2OMOY,YHEMOY,
  433. & YCO2MOY,YCOMOY,PRESS,TMOY,
  434. & TS,TPLA,DEB,YH2SOR,YO2SOR,YN2SOR,YH2OSOR,YHESOR,
  435. & YCO2SOR,YCOSOR,Dt,KALP,IQ,MLREEL,IHI)
  436. ENDIF
  437. IF (IHI .LT. 0) THEN
  438. CALL ERREUR(944)
  439. RETURN
  440. ENDIF
  441. C
  442. C Ecriture des paramètres
  443. C XH2ONF : flottant : Flag pour l'allumage (négatif si le PAR est ON)
  444. C EFF : flottant : Efficacité du PAR (0.6 ou 1 suivant deficit en O2)
  445. C DEB : flottant : Débit au travers du PAR (en kg/s)
  446. C TPLA : flottant : Température moyenne des plaques du PAR (en K)
  447. C TS : flottant : Température moyenne du gaz en sortie du PAR (en K)
  448. C YCOSOR : flottant : Fraction massique moyenne de CO en sortie du PAR
  449. C YCO2SOR : flottant : Fraction massique moyenne de CO2 en sortie du PAR
  450. C YHESOR : flottant : Fraction massique moyenne de HE en sortie du PAR
  451. C YH2OSOR : flottant : Fraction massique moyenne de H2O en sortie du PAR
  452. C YN2SOR : flottant : Fraction massique moyenne de N2 en sortie du PAR
  453. C YO2SOR : flottant : Fraction massique moyenne de O2 en sortie du PAR
  454. C YH2SOR : flottant : Fraction massique moyenne de H2 en sortie du PAR
  455. IF (IQ .EQ. 1) THEN
  456. CALL ECRREE(XH2ONF)
  457. CALL ECRREE(EFF)
  458. CALL ECRREE(DEB)
  459. CALL ECRREE(TPLA)
  460. CALL ECRREE(TS)
  461. CALL ECRREE(YCOSOR)
  462. CALL ECRREE(YCO2SOR)
  463. CALL ECRREE(YHESOR)
  464. CALL ECRREE(YH2OSOR)
  465. CALL ECRREE(YN2SOR)
  466. CALL ECRREE(YO2SOR)
  467. CALL ECRREE(YH2SOR)
  468. ELSE
  469. C VERIF : Pour le modèle thermique, les compositions entrée/sortie sont égales
  470. CALL ECRREE(DEB)
  471. CALL ECRREE(TPLA)
  472. CALL ECRREE(TS)
  473. CALL ECRREE(YCOSOR)
  474. CALL ECRREE(YCO2SOR)
  475. CALL ECRREE(YHESOR)
  476. CALL ECRREE(YH2OSOR)
  477. CALL ECRREE(YN2SOR)
  478. CALL ECRREE(YO2SOR)
  479. CALL ECRREE(YH2SOR)
  480. ENDIF
  481. C
  482. C===================================
  483. C Cinetique de type HEAVYSIDE
  484. C===================================
  485. C
  486. ELSEIF( IP .EQ. 6) THEN
  487. CALL FLAMA()
  488. C
  489. ELSE
  490. C
  491. C******* Message d'erreur standard
  492. C Tentative d'utilisation d'une option non implémentée
  493. C 251 2
  494. C
  495. CALL ERREUR(251)
  496. RETURN
  497. ENDIF
  498. C
  499. C
  500. RETURN
  501. 100 FORMAT(2X,'Les CHPOINTs n''ont pas le meme SPG')
  502. END
  503.  
  504.  
  505.  
  506.  
  507.  
  508.  
  509.  
  510.  
  511.  
  512.  
  513.  
  514.  
  515.  
  516.  
  517.  

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