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 40 : CALL LIRCHA(MOTOP,1,IRET) 40 : CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT) 41 : C 42 : C====================================================== 43 : C Branchements sur les différents modèles de combustion 44 : C====================================================== 45 : C 46 : C================================= 47 : C Cas de la déflagration laminaire 48 : C================================= 49 : C 50 : IF (IP .EQ. 1) THEN 51 : C 52 : C Lecture des paramètres 53 : C 53 : CALL LIROBJ('CHPOINT ',MRHO,1,IRET) 54 : IF (IRET.EQ.0) RETURN 54 : CALL LIROBJ('CHPOINT ',MCV,1,IRET) 55 : IF (IRET.EQ.0) RETURN 56 : C 57 : C Lecture de l'option sur la forme des Cv 58 : C 58 : CALL LIRCHA(MOTOP,1,IRET) 58 : CALL OPTLI(IP,COPTIO,MOTOP,NBCV) 59 : IF (IP .EQ. 0) THEN 59 : CALL ERREUR(251) 60 : RETURN 61 : ENDIF 62 : IF ( MOTOP .EQ. 'LINEAIRE') THEN 63 : NCVOPT = 1 64 : ELSEIF( MOTOP .EQ. 'QUADRATI') THEN 65 : NCVOPT = 2 66 : ENDIF 67 : C 67 : CALL LIROBJ('CHPOINT ',MR,1,IRET) 68 : IF (IRET.EQ.0) RETURN 68 : CALL LIROBJ('CHPOINT ',MYH2,1,IRET) 69 : IF (IRET.EQ.0) RETURN 69 : CALL LIROBJ('CHPOINT ',MYO2,1,IRET) 70 : IF (IRET.EQ.0) RETURN 70 : CALL LIROBJ('CHPOINT ',MYN2,1,IRET) 71 : IF (IRET.EQ.0) RETURN 71 : CALL LIROBJ('CHPOINT ',MYH2O,1,IRET) 72 : IF (IRET.EQ.0) RETURN 72 : CALL LIROBJ('CHPOINT ',MTEMP,1,IRET) 73 : IF (IRET.EQ.0) RETURN 73 : CALL LIRREE(Dt,1,IRET) 74 : IF (IRET.EQ.0) RETURN 75 : C 76 : C Activation champoints et vérification des pointeurs 77 : C 77 : CALL LICHT(MRHO,NRHO,TYPE,IGEO) 77 : CALL LICHT(MCV,NCV,TYPE,IGEO0) 78 : IF (IGEO.NE.IGEO0) THEN 79 : WRITE(IOIMP,100) 80 : RETURN 81 : ENDIF 81 : CALL LICHT(MR,NR,TYPE,IGEO0) 82 : IF (IGEO.NE.IGEO0) THEN 83 : WRITE(IOIMP,100) 84 : RETURN 85 : ENDIF 85 : CALL LICHT(MYH2,NYH2,TYPE,IGEO0) 86 : IF (IGEO.NE.IGEO0) THEN 87 : WRITE(IOIMP,100) 88 : RETURN 89 : ENDIF 89 : CALL LICHT(MYO2,NYO2,TYPE,IGEO0) 90 : IF (IGEO.NE.IGEO0) THEN 91 : WRITE(IOIMP,100) 92 : RETURN 93 : ENDIF 93 : CALL LICHT(MYN2,NYN2,TYPE,IGEO0) 94 : IF (IGEO.NE.IGEO0) THEN 95 : WRITE(IOIMP,100) 96 : RETURN 97 : ENDIF 97 : CALL LICHT(MYH2O,NYH2O,TYPE,IGEO0) 98 : IF (IGEO.NE.IGEO0) THEN 99 : WRITE(IOIMP,100) 100 : RETURN 101 : ENDIF 101 : CALL LICHT(MTEMP,NTEMP,TYPE,IGEO0) 102 : IF (IGEO.NE.IGEO0) THEN 103 : WRITE(IOIMP,100) 104 : RETURN 105 : ENDIF 106 : C 107 : C Création et activation des chpoints résultats 108 : C 109 : NC=1 109 : CALL CRCHPT(TYPE,IGEO,NC,MZH2) 109 : CALL LICHT(MZH2,NZH2,TYPE,IGEO) 109 : CALL CRCHPT(TYPE,IGEO,NC,MZO2) 109 : CALL LICHT(MZO2,NZO2,TYPE,IGEO) 109 : CALL CRCHPT(TYPE,IGEO,NC,MZN2) 109 : CALL LICHT(MZN2,NZN2,TYPE,IGEO) 109 : CALL CRCHPT(TYPE,IGEO,NC,MZH2O) 109 : CALL LICHT(MZH2O,NZH2O,TYPE,IGEO) 109 : CALL CRCHPT(TYPE,IGEO,NC,MQ) 109 : CALL LICHT(MQ,NQ,TYPE,IGEO) 110 : C 111 : C Nombre de Points du Maillage 112 : C 113 : NBPTS = NRHO.VPOCHA(/1) 114 : C 115 : C Calcul des nouvelles concentrations et de l'énergie libérée 116 : C 116 : CALL FLAM1(NBPTS,NRHO.VPOCHA,NCV.VPOCHA,NR.VPOCHA,NYH2.VPOCHA, 117 : 1 NYO2.VPOCHA,NYN2.VPOCHA,NYH2O.VPOCHA,NTEMP.VPOCHA, 118 : 2 NZH2.VPOCHA,NZO2.VPOCHA,NZN2.VPOCHA,NZH2O.VPOCHA, 119 : 3 NQ.VPOCHA,Dt,NCVOPT) 120 : C 121 : C Ecriture des chpoints résultats 122 : C 122 : CALL ECROBJ('CHPOINT ',MQ) 122 : CALL ECROBJ('CHPOINT ',MZH2O) 122 : CALL ECROBJ('CHPOINT ',MZN2) 122 : CALL ECROBJ('CHPOINT ',MZO2) 122 : CALL ECROBJ('CHPOINT ',MZH2) 123 : C 124 : C Désactivation des segments 125 : C 126 : SEGDES NZH2,NZO2,NZN2,NZH2O,NQ,NRHO,NCV,NR,NYH2 127 : SEGDES NYO2,NYN2,NYH2O,NTEMP 128 : C 129 : C 130 : C================================== 131 : C Cas de la déflagration turbulente 132 : C================================== 133 : C 134 : ELSEIF( IP .EQ. 2) THEN 135 : C 136 : C Lecture des paramètres 137 : C 137 : CALL LIROBJ('CHPOINT ',MRHO,1,IRET) 138 : IF (IRET.EQ.0) RETURN 138 : CALL LIROBJ('CHPOINT ',MCV,1,IRET) 139 : IF (IRET.EQ.0) RETURN 140 : C 141 : C Lecture de l'option sur la forme des Cv 142 : C 142 : CALL LIRCHA(MOTOP,1,IRET) 142 : CALL OPTLI(IP,COPTIO,MOTOP,NBCV) 143 : IF(IP .EQ. 0) THEN 143 : CALL ERREUR(251) 144 : RETURN 145 : ENDIF 146 : IF ( MOTOP .EQ. 'LINEAIRE') THEN 147 : NCVOPT = 1 148 : ELSEIF( MOTOP .EQ. 'QUADRATI') THEN 149 : NCVOPT = 2 150 : ENDIF 151 : C 151 : CALL LIROBJ('CHPOINT ',MR,1,IRET) 152 : IF (IRET.EQ.0) RETURN 152 : CALL LIROBJ('CHPOINT ',MYH2,1,IRET) 153 : IF (IRET.EQ.0) RETURN 153 : CALL LIROBJ('CHPOINT ',MYO2,1,IRET) 154 : IF (IRET.EQ.0) RETURN 154 : CALL LIROBJ('CHPOINT ',MYN2,1,IRET) 155 : IF (IRET.EQ.0) RETURN 155 : CALL LIROBJ('CHPOINT ',MYH2O,1,IRET) 156 : IF (IRET.EQ.0) RETURN 156 : CALL LIROBJ('CHPOINT ',MTEMP,1,IRET) 157 : IF (IRET.EQ.0) RETURN 157 : CALL LIRREE(Dt,1,IRET) 158 : IF (IRET.EQ.0) RETURN 159 : 160 : C paramètres spécifiques à l'option choisie 161 : 161 : CALL LIROBJ('CHPOINT ',MZK,1,IRET) 162 : IF (IRET.EQ.0) RETURN 162 : CALL LIROBJ('CHPOINT ',MZE,1,IRET) 163 : IF (IRET.EQ.0) RETURN 163 : CALL LIROBJ('CHPOINT ',MYH2u,1,IRET) 164 : IF (IRET.EQ.0) RETURN 164 : CALL LIRREE(C_EBU,1,IRET) 165 : IF (IRET.EQ.0) RETURN 165 : CALL LIRREE(CMIN,0,IRET) 166 : IF (IRET.EQ.0) CMIN=1.D-4 167 : C 168 : C Activation champoints et vérification des pointeurs 169 : C 169 : CALL LICHT(MRHO,NRHO,TYPE,IGEO) 169 : CALL LICHT(MCV,NCV,TYPE,IGEO0) 170 : IF (IGEO.NE.IGEO0) THEN 171 : WRITE(IOIMP,100) 172 : RETURN 173 : ENDIF 173 : CALL LICHT(MR,NR,TYPE,IGEO0) 174 : IF (IGEO.NE.IGEO0) THEN 175 : WRITE(IOIMP,100) 176 : RETURN 177 : ENDIF 177 : CALL LICHT(MYH2,NYH2,TYPE,IGEO0) 178 : IF (IGEO.NE.IGEO0) THEN 179 : WRITE(IOIMP,100) 180 : RETURN 181 : ENDIF 181 : CALL LICHT(MYO2,NYO2,TYPE,IGEO0) 182 : IF (IGEO.NE.IGEO0) THEN 183 : WRITE(IOIMP,100) 184 : RETURN 185 : ENDIF 185 : CALL LICHT(MYN2,NYN2,TYPE,IGEO0) 186 : IF (IGEO.NE.IGEO0) THEN 187 : WRITE(IOIMP,100) 188 : RETURN 189 : ENDIF 189 : CALL LICHT(MYH2O,NYH2O,TYPE,IGEO0) 190 : IF (IGEO.NE.IGEO0) THEN 191 : WRITE(IOIMP,100) 192 : RETURN 193 : ENDIF 193 : CALL LICHT(MTEMP,NTEMP,TYPE,IGEO0) 194 : IF (IGEO.NE.IGEO0) THEN 195 : WRITE(IOIMP,100) 196 : RETURN 197 : ENDIF 197 : CALL LICHT(MZK,NZK,TYPE,IGEO0) 198 : IF (IGEO.NE.IGEO0) THEN 199 : WRITE(IOIMP,100) 200 : RETURN 201 : ENDIF 201 : CALL LICHT(MZE,NZE,TYPE,IGEO0) 202 : IF (IGEO.NE.IGEO0) THEN 203 : WRITE(IOIMP,100) 204 : RETURN 205 : ENDIF 205 : CALL LICHT(MYH2u,NYH2u,TYPE,IGEO0) 206 : IF (IGEO.NE.IGEO0) THEN 207 : WRITE(IOIMP,100) 208 : RETURN 209 : ENDIF 210 : C 211 : C Création et activation des chpoints résultats 212 : C 213 : NC=1 213 : CALL CRCHPT(TYPE,IGEO,NC,MZH2) 213 : CALL LICHT(MZH2,NZH2,TYPE,IGEO) 213 : CALL CRCHPT(TYPE,IGEO,NC,MZO2) 213 : CALL LICHT(MZO2,NZO2,TYPE,IGEO) 213 : CALL CRCHPT(TYPE,IGEO,NC,MZN2) 213 : CALL LICHT(MZN2,NZN2,TYPE,IGEO) 213 : CALL CRCHPT(TYPE,IGEO,NC,MZH2O) 213 : CALL LICHT(MZH2O,NZH2O,TYPE,IGEO) 213 : CALL CRCHPT(TYPE,IGEO,NC,MQ) 213 : CALL LICHT(MQ,NQ,TYPE,IGEO) 214 : C 215 : C Nombre de Points 216 : C 217 : NBPTS = NRHO.VPOCHA(/1) 218 : C 218 : CALL FLAM2(NBPTS,NRHO.VPOCHA,NCV.VPOCHA,NR.VPOCHA,NYH2.VPOCHA, 219 : 1 NYO2.VPOCHA,NYN2.VPOCHA,NYH2O.VPOCHA,NTEMP.VPOCHA, 220 : 2 NZK.VPOCHA,NZE.VPOCHA,NYH2u.VPOCHA, 221 : 3 NZH2.VPOCHA,NZO2.VPOCHA,NZN2.VPOCHA,NZH2O.VPOCHA, 222 : 4 NQ.VPOCHA,C_EBU,Dt,CMIN,NCVOPT) 223 : C 224 : C Ecriture des chpoints résultats 225 : C 225 : CALL ECROBJ('CHPOINT ',MQ) 225 : CALL ECROBJ('CHPOINT ',MZH2O) 225 : CALL ECROBJ('CHPOINT ',MZN2) 225 : CALL ECROBJ('CHPOINT ',MZO2) 225 : CALL ECROBJ('CHPOINT ',MZH2) 226 : C 227 : C Désactivation des segments 228 : C 229 : SEGDES NZH2,NZO2,NZN2,NZH2O,NQ,NRHO,NCV,NR,NYH2 230 : SEGDES NYO2,NYN2,NYH2O,NTEMP,NZK,NZE,NYH2u 231 : C 232 : C=================================== 233 : C Cas de la loi d'Arrhénius générale 234 : C=================================== 235 : C 236 : ELSEIF( IP .EQ. 3) THEN 236 : CALL FLAMB() 237 : C 238 : C=================================== 239 : C Maquette CREBCOM 240 : C=================================== 241 : C 242 : ELSEIF( IP .EQ. 7) THEN 242 : CALL FLAMCR() 243 : IF(IERR.NE.0)RETURN 244 : C 245 : C====================================== 246 : C Maquette CREBCOM, modele non-homogene 247 : C====================================== 248 : C 249 : ELSEIF( IP .EQ. 8) THEN 249 : CALL FLACR2() 250 : IF(IERR.NE.0)RETURN 251 : C 252 : C============================================ 253 : C Combustion catalytique : modèle recombineur 254 : C============================================ 255 : C 256 : ELSEIF((IP .EQ. 4) .OR. (IP .EQ. 5)) THEN 257 : C 258 : C Lecture des paramètres 259 : C YH2MOY : flottant : Fraction massique moyenne de H2 à l'entrée du PAR 260 : C YO2MOY : flottant : Fraction massique moyenne de O2 à l'entrée du PAR 261 : C YN2MOY : flottant : Fraction massique moyenne de N2 à l'entrée du PAR 262 : C YH2OMOY : flottant : Fraction massique moyenne de H2O à l'entrée du PAR 263 : C YHEMOY : flottant : Fraction massique moyenne de HE à l'entrée du PAR 264 : C YCO2MOY : flottant : Fraction massique moyenne de CO2 à l'entrée du PAR 265 : C YCOMOY : flottant : Fraction massique moyenne de CO à l'entrée du PAR 266 : C TMOY : flottant : Température moyenne du gaz à l'entrée du PAR (K) 267 : C TPLA : flottant : Température moyenne des plaques du PAR (K) 268 : C PRESS : flottant : Pression moyenne autour du PAR (Pa) 269 : C DT : flottant : Pas de temps (s) 270 : C KALP : entier : Flag pour correction du coefficient d'échange 271 : C comme dans RALOC (voir fla009.eso) 271 : CALL LIRREE(YH2MOY,1,IRET) 272 : IF (IRET.EQ.0) RETURN 272 : CALL LIRREE(YO2MOY,1,IRET) 273 : IF (IRET.EQ.0) RETURN 273 : CALL LIRREE(YN2MOY,1,IRET) 274 : IF (IRET.EQ.0) RETURN 274 : CALL LIRREE(YH2OMOY,1,IRET) 275 : IF (IRET.EQ.0) RETURN 275 : CALL LIRREE(YHEMOY,1,IRET) 276 : IF (IRET.EQ.0) RETURN 276 : CALL LIRREE(YCO2MOY,1,IRET) 277 : IF (IRET.EQ.0) RETURN 277 : CALL LIRREE(YCOMOY,1,IRET) 278 : IF (IRET.EQ.0) RETURN 278 : CALL LIRREE(TMOY,1,IRET) 279 : IF (IRET.EQ.0) RETURN 279 : CALL LIRREE(TPLA,1,IRET) 280 : IF (IRET.EQ.0) RETURN 280 : CALL LIRREE(PRESS,1,IRET) 281 : IF (IRET.EQ.0) RETURN 281 : CALL LIRREE(Dt,1,IRET) 282 : IF (IRET.EQ.0) RETURN 282 : CALL LIRENT(KALP,1,IRET) 283 : IF (IRET.EQ.0) RETURN 284 : C Lecture du modèle de PAR à choisir parmi 'SIEMENS ' ou 'HEATER ' 285 : C MODEL : CHAR*8 : Nom du modèle de PAR 285 : CALL LIRCHA(MODEL,1,IRET) 286 : IF (IRET.EQ.0) RETURN 286 : CALL OPTLI(IQ,MODREC,MODEL,NBMOD) 287 : IF (IQ.EQ.1) THEN 288 : C Lecture des paramètres A et B du modele recombineur 'SIEMENS ' 289 : C (loi lineaire en pression) 289 : CALL LIRREE(A1,1,IRET) 290 : IF (IRET.EQ.0) RETURN 290 : CALL LIRREE(B1,1,IRET) 291 : IF (IRET.EQ.0) RETURN 291 : CALL LIRREE(XH2ON,1,IRET) 292 : IF (IRET.EQ.0) RETURN 292 : CALL LIRREE(XH2OFF,1,IRET) 293 : IF (IRET.EQ.0) RETURN 293 : CALL LIRREE(BETA,1,IRET) 294 : IF (IRET.EQ.0) RETURN 295 : JG = 5 296 : SEGINI MLREEL 297 : PROG(1) = A1 298 : PROG(2) = B1 299 : PROG(3) = XH2ON 300 : PROG(4) = XH2OFF 301 : PROG(5) = BETA 302 : C 303 : C Calcul des conditions de sortie, de la temperature des plaques et du 304 : C débit dans le PAR (arret si débit ou température des plaques négatif) 304 : CALL FLAMR(YH2MOY,YO2MOY,YN2MOY,YH2OMOY,YHEMOY, 305 : & YCO2MOY,YCOMOY,PRESS,TMOY, 306 : & TS,TPLA,DEB,YH2SOR,YO2SOR,YN2SOR,YH2OSOR,YHESOR, 307 : & YCO2SOR,YCOSOR,XH2ONF,EFF,Dt,KALP,IQ,MLREEL,IHI) 308 : ELSEIF (IQ.EQ.2) THEN 309 : C Lecture de la puissance du recombineur thermique (modèle 'HEATER ') 309 : CALL LIRREE(QHEAT,1,IRET) 310 : IF (IRET.EQ.0) RETURN 310 : CALL LIRREE(FUITE,1,IRET) 311 : IF (IRET.EQ.0) RETURN 311 : CALL LIRREE(HRAPP,1,IRET) 312 : IF (IRET.EQ.0) RETURN 312 : CALL LIRREE(CRAPP,1,IRET) 313 : IF (IRET.EQ.0) RETURN 314 : JG = 4 315 : SEGINI MLREEL 316 : PROG(1) = QHEAT 317 : PROG(2) = 1.D0 - fuite 318 : PROG(3) = hrapp 319 : PROG(4) = crapp 320 : C 321 : C Calcul des conditions de sortie, de la temperature des plaques et du 322 : C débit dans le HEATER (arret si débit ou température des plaques négatif) 322 : CALL FLAMH(YH2MOY,YO2MOY,YN2MOY,YH2OMOY,YHEMOY, 323 : & YCO2MOY,YCOMOY,PRESS,TMOY, 324 : & TS,TPLA,DEB,YH2SOR,YO2SOR,YN2SOR,YH2OSOR,YHESOR, 325 : & YCO2SOR,YCOSOR,Dt,KALP,IQ,MLREEL,IHI) 326 : ENDIF 327 : IF (IHI .LT. 0) THEN 327 : CALL ERREUR(944) 328 : RETURN 329 : ENDIF 330 : C 331 : C Ecriture des paramètres 332 : C XH2ONF : flottant : Flag pour l'allumage (négatif si le PAR est ON) 333 : C EFF : flottant : Efficacité du PAR (0.6 ou 1 suivant deficit en O2) 334 : C DEB : flottant : Débit au travers du PAR (en kg/s) 335 : C TPLA : flottant : Température moyenne des plaques du PAR (en K) 336 : C TS : flottant : Température moyenne du gaz en sortie du PAR (en K) 337 : C YCOSOR : flottant : Fraction massique moyenne de CO en sortie du PAR 338 : C YCO2SOR : flottant : Fraction massique moyenne de CO2 en sortie du PAR 339 : C YHESOR : flottant : Fraction massique moyenne de HE en sortie du PAR 340 : C YH2OSOR : flottant : Fraction massique moyenne de H2O en sortie du PAR 341 : C YN2SOR : flottant : Fraction massique moyenne de N2 en sortie du PAR 342 : C YO2SOR : flottant : Fraction massique moyenne de O2 en sortie du PAR 343 : C YH2SOR : flottant : Fraction massique moyenne de H2 en sortie du PAR 344 : IF (IQ .EQ. 1) THEN 344 : CALL ECRREE(XH2ONF) 344 : CALL ECRREE(EFF) 344 : CALL ECRREE(DEB) 344 : CALL ECRREE(TPLA) 344 : CALL ECRREE(TS) 344 : CALL ECRREE(YCOSOR) 344 : CALL ECRREE(YCO2SOR) 344 : CALL ECRREE(YHESOR) 344 : CALL ECRREE(YH2OSOR) 344 : CALL ECRREE(YN2SOR) 344 : CALL ECRREE(YO2SOR) 344 : CALL ECRREE(YH2SOR) 345 : ELSE 346 : C VERIF : Pour le modèle thermique, les compositions entrée/sortie sont égales 346 : CALL ECRREE(DEB) 346 : CALL ECRREE(TPLA) 346 : CALL ECRREE(TS) 346 : CALL ECRREE(YCOSOR) 346 : CALL ECRREE(YCO2SOR) 346 : CALL ECRREE(YHESOR) 346 : CALL ECRREE(YH2OSOR) 346 : CALL ECRREE(YN2SOR) 346 : CALL ECRREE(YO2SOR) 346 : CALL ECRREE(YH2SOR) 347 : ENDIF 348 : C 349 : C=================================== 350 : C Cinetique de type HEAVYSIDE 351 : C=================================== 352 : C 353 : ELSEIF( IP .EQ. 6) THEN 353 : CALL FLAMA() 354 : C 355 : ELSE 356 : C 357 : C******* Message d'erreur standard 358 : C Tentative d'utilisation d'une option non implémentée 359 : C 251 2 360 : C 360 : CALL ERREUR(251) 361 : RETURN 362 : ENDIF 363 : C 364 : C 365 : RETURN 366 : 100 FORMAT(2X,'Les CHPOINTs n''ont pas le meme SPG') 367 : END 368 : 369 : 370 : 371 : 372 : 373 : 374 : 375 : 376 : 377 : 378 : 379 : 380 : 381 :
© Cast3M 2003 - All rights reserved.
Disclaimer