1 : C PR12F SOURCE BECC 05/12/02 21:15:52 5258 2 : SUBROUTINE PR12f(IM1, ICH1, ICH2, ICH3, 3 : & ICH4, ICH5, ICH6, ICH7, 4 : & ICH8, ICH9, 5 : & Cp, Cvm, 6 : & OCH1, OCH2, OCH3, 7 : & OCH4, OCH5, OCH6, 8 : & OCH7, OCH8, OCH9, 9 : & LOGNEG, LOGBOR, MESERR, 10 : & VALER, VAL1, VAL2) 11 : C************************************************************************ 12 : C 13 : C PROJET : CASTEM 2000 14 : C 15 : C NOM : PR12f 16 : C 17 : C DESCRIPTION : Primitive variables from conserved variables 18 : C 19 : C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI) 20 : C 21 : C AUTEUR : Jose R. Garcia-Cascales, 22 : C Universidad Politecnica de Cartagena, 23 : C jr.garcia@upct.es 24 : C 25 : C************************************************************************ 26 : C 27 : C APPELES (E/S) : 28 : C 29 : C************************************************************************ 30 : C 31 : C ENTREES : 32 : C 33 : C IM1 : MELEME contenant les centres des ELTs 34 : C 35 : C ICH1 : CHPOINT contenant la masse volumique (mass, gas/vapour) 36 : C 37 : C ICH2 : CHPOINT contenant la masse volumique (mass, liquid) 38 : C 39 : C ICH3 : CHPOINT contenant les dèbits (momentum, gas/vapour) 40 : C ( NDIM composantes); 41 : C 42 : C ICH4 : CHPOINT contenant les dèbits (momentum, liquid) 43 : C ( NDIM composantes); 44 : C 45 : C ICH5 : CHPOINT contenat l'énergie totale per 46 : C unité de volume (alpha rho et) (gas/vapour); 47 : C 48 : C ICH6 : CHPOINT contenat l'énergie totale per 49 : C unité de volume (alpha rho et) (liquid); 50 : C 51 : C ICH7 : CHPOINT containing void fraction of the previous 52 : C time step, in order to avoid problems when a 53 : C phase disappears 54 : C 55 : C ICH8 : CHPOINT containing vapour temperature of the previous 56 : C time step, in order to avoid problems when a 57 : C phase disappears 58 : C 59 : C ICH9 : CHPOINT containing liquid temperature of the previous 60 : C time step, in order to avoid problems when the liquid 61 : C phase disappears 62 : C 63 : C Cp : Parameter characterizing the pressure 64 : C correction term 65 : C 66 : C Cvm : Parameter characterizing virtual mass correction, 67 : C it is included in the new CATHARE pressure 68 : C conrrection term 69 : C 70 : C SORTIES : OCH1 : CHPOINT void fraction (alpha); 71 : C 72 : C OCH2 : CHPOINT gas/vapour velocity (uvx, uvy, uvz); 73 : C 74 : C OCH3 : CHPOINT liquid velocity (ulx, uly, ulz); 75 : C 76 : C OCH4 : CHPOINT pressure; 77 : C 78 : C OCH5 : CHPOINT vapour temperature (Tv); 79 : C 80 : C OCH6 : CHPOINT liquid temperature (Tl); 81 : C 82 : C OCH7 : CHPOINT vapour density (rv); 83 : C 84 : C OCH8 : CHPOINT liquid density (rl); 85 : C 86 : C OCH9 : CHPOINT interfacial pressure correction term (pi) 87 : C 88 : C LOGNEG : (LOGICAL): si .TRUE. une pression ou une densité 89 : C negative a été detectée -> le programme s'arrete 90 : C (sa valeur stockée en MESERR(1) et VALER(1)) 91 : C 92 : C LOGBOR : (LOGICAL) 93 : C gamma a été detecté dehor GAMMIN et GAMMAX 94 : C (sa valeur stockée en MESERR(2) et VALER(2),VAL1,VAL2) 95 : C 96 : C MESERR, 97 : C VALER, 98 : C VAL1, 99 : C VAL2 : pour message d'erreur 100 : C 101 : C 102 : C 103 : C************************************************************************ 104 : C 105 : C HISTORIQUE (Anomalies et modifications éventuelles) 106 : C 107 : C HISTORIQUE : Créée le 20.11.05. 108 : C 109 : C************************************************************************ 110 : C 111 : C 112 : C**** Variables de COOPTIO 113 : C 114 : INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR 115 : & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ 116 : & ,IOPER, IOSGB, IOGRA, IOSAU, IORES 117 : & ,IECHO, IIMPI, IOSPI 118 : & ,IDIM, IFICLE, IPREFI 119 : & ,MCOORD 120 : & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE 121 : & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU 122 : & ,NORINC,NORVAL,NORIND,NORVAD 123 : & ,NUCROU, IPSAUV 124 : C 125 : C**** Les variables 126 : C 127 : INTEGER NLCE,N1,IGEOMC 128 : C 129 : INTEGER IM1 130 : & ,ICH1, ICH2, ICH3 131 : & ,ICH4, ICH5, ICH6 132 : & ,ICH7, ICH8, ICH9 133 : & ,OCH1,OCH2,OCH3 134 : & ,OCH4,OCH5,OCH6 135 : & ,OCH7,OCH8,OCH9 136 : C 137 : REAL*8 VALER(2),VAL1,VAL2 138 : & ,alpha, uvx, uvy, uvz, uv2, ulx, uly, ulz, ul2 139 : & ,p, Tv, Tl, rv, rl, ev, el 140 : & ,oldav, oldal, oldTv, oldTl 141 : & ,A, B, delta 142 : & ,gamv, cpv, gaml, cpl, pil 143 : & ,FFV, FFL 144 : & ,rm, rvm, rlm, Cp, cvm, pi 145 : PARAMETER(gamv = 1.4D0, cpv = 1008.D0, 146 : & gaml = 2.8D0, cpl = 4186.D0, pil = 8.5D8) 147 : C 148 : CHARACTER*(8) TYPE 149 : CHARACTER*(40) MESERR(2) 150 : LOGICAL LOGNEG, LOGBOR 151 : C 152 : INTEGER NAT, NSOUPO, N, NC 153 : C 154 : C**** Les includes 155 : C 156 : -INC CCOPTIO 157 : -INC SMCHPOI 158 : POINTEUR MCHPO5.MCHPOI, MCHPO6.MCHPOI, 159 : & MCHPO7.MCHPOI, MCHPO8.MCHPOI, MCHPO9.MCHPOI 160 : C 161 : POINTEUR MSOUP6.MSOUPO, MSOUP7.MSOUPO, MSOUP8.MSOUPO, 162 : & MSOUP9.MSOUPO 163 : C 164 : POINTEUR MPOW1.MPOVAL, MPOW2.MPOVAL, MPOW3.MPOVAL, 165 : & MPOW4.MPOVAL, MPOW5.MPOVAL, MPOW6.MPOVAL, 166 : & MPOW7.MPOVAL, MPOW8.MPOVAL, MPOW9.MPOVAL 167 : C 168 : POINTEUR MPOV1.MPOVAL, MPOV2.MPOVAL, MPOV3.MPOVAL, 169 : & MPOV4.MPOVAL, MPOV5.MPOVAL, MPOV6.MPOVAL, 170 : & MPOV7.MPOVAL, MPOV8.MPOVAL, MPOV9.MPOVAL 171 : -INC SMELEME 172 : C 173 : C 174 : C**** Initialisation des variables pour la gestion des erreurs pas ici, 175 : C mais avant, i.e. 176 : C 177 : LOGNEG = .FALSE. 178 : LOGBOR = .FALSE. 179 : MESERR(1) = ' ' 180 : MESERR(2) = ' ' 181 : VALER(1) = 0.0D0 182 : VALER(2) = 0.0D0 183 : VAL1 = 0.0D0 184 : VAL2 = 0.0D0 185 : 186 : C 187 : C**** Activation du MELEME "CENTRE" 188 : C 189 : IPT1 = IM1 190 : SEGACT IPT1 191 : N1 = IPT1.NUM(/2) 192 : SEGDES IPT1 193 : C 194 : C**** Reading MPOVA of each MCHAMPOIN 195 : C MPOW1 = alpha*rv 196 : C MPOW2 = (1 - alpha)*rl 197 : C MPOW3 = alpha*rv*uvx & alpha*rv*uvy & alpha*rv*uvz 198 : C MPOW4 = (1 - alpha)*rl*ulx & (1 - alpha)*rl*uly & (1 - alpha)*rl*ulz 199 : C MPOW5 = alpha*rv*Ev 200 : C MPOW6 = (1 - alpha)*rl*El 201 : C 202 : C MPOW7: It is the old void fraction, we will use it when 203 : C alpha approaches to zero 204 : C MPOW8: It is the old vapour temperature, we will use it when 205 : C alpha approaches to zero 206 : C MPOW9: It is the old liquid temperature, we will use it when 207 : C alpha approaches to one 208 : C 208 : CALL LICHT(ICH1,MPOW1,TYPE,IGEOMC) 208 : CALL LICHT(ICH2,MPOW2,TYPE,IGEOMC) 208 : CALL LICHT(ICH3,MPOW3,TYPE,IGEOMC) 208 : CALL LICHT(ICH4,MPOW4,TYPE,IGEOMC) 208 : CALL LICHT(ICH5,MPOW5,TYPE,IGEOMC) 208 : CALL LICHT(ICH6,MPOW6,TYPE,IGEOMC) 208 : CALL LICHT(ICH7,MPOW7,TYPE,IGEOMC) 208 : CALL LICHT(ICH8,MPOW8,TYPE,IGEOMC) 208 : CALL LICHT(ICH9,MPOW9,TYPE,IGEOMC) 209 : C 210 : C SEGACT MPOW1 211 : C SEGACT MPOW2 212 : C SEGACT MPOW3 213 : C SEGACT MPOW4 214 : C SEGACT MPOW5 215 : C SEGACT MPOW6 216 : C SEGACT MPOW7 217 : C SEGACT MPOW8 218 : C SEGACT MPOW9 219 : CC 220 : C**** LICHT active les MPOVALs en *MO 221 : C 222 : C 223 : C**** Privitive variables 224 : C 225 : C void fraction: alpha 226 : 227 : NAT = 1 228 : NSOUPO=1 229 : N = N1 230 : NC = 1 231 : SEGINI MCHPO1 232 : MCHPO1.JATTRI(1)=2 233 : MCHPO1.IFOPOI=IFOMOD 234 : SEGINI MSOUP1 235 : MCHPO1.IPCHP(1)=MSOUP1 236 : SEGINI MPOV1 237 : MSOUP1.NOCOMP(1)='SCAL' 238 : MSOUP1.IGEOC = IM1 239 : MSOUP1.IPOVAL = MPOV1 240 : 241 : C vapour velocities: uvx, uvy, uvz 242 : 243 : 244 : NC = IDIM 245 : SEGINI MCHPO2 246 : MCHPO2.JATTRI(1)=2 247 : MCHPO2.IFOPOI=IFOMOD 248 : SEGINI MSOUP2 249 : MCHPO2.IPCHP(1)=MSOUP2 250 : SEGINI MPOV2 251 : MSOUP2.NOCOMP(1)='UVX' 252 : MSOUP2.NOCOMP(2)='UVY' 253 : IF(IDIM .EQ. 3)THEN 254 : MSOUP2.NOCOMP(3)='UVZ' 255 : ENDIF 256 : MSOUP2.IGEOC = IM1 257 : MSOUP2.IPOVAL = MPOV2 258 : 259 : C liquid velocities: ulx, uly, ulz 260 : 261 : 262 : NC = IDIM 263 : SEGINI MCHPO3 264 : MCHPO3.JATTRI(1)=2 265 : MCHPO3.IFOPOI=IFOMOD 266 : SEGINI MSOUP3 267 : MCHPO3.IPCHP(1)=MSOUP3 268 : SEGINI MPOV3 269 : MSOUP3.NOCOMP(1)='ULX' 270 : MSOUP3.NOCOMP(2)='ULY' 271 : IF(IDIM .EQ. 3)THEN 272 : MSOUP3.NOCOMP(3)='ULZ' 273 : ENDIF 274 : MSOUP3.IGEOC = IM1 275 : MSOUP3.IPOVAL = MPOV3 276 : 277 : C pressure: p 278 : 279 : NC = 1 280 : SEGINI MCHPO4 281 : MCHPO4.JATTRI(1)=2 282 : MCHPO4.IFOPOI=IFOMOD 283 : SEGINI MSOUP4 284 : MCHPO4.IPCHP(1)=MSOUP4 285 : SEGINI MPOV4 286 : MSOUP4.NOCOMP(1)='SCAL' 287 : MSOUP4.IGEOC = IM1 288 : MSOUP4.IPOVAL = MPOV4 289 : 290 : C vapour temperature: Tv 291 : 292 : NC = 1 293 : SEGINI MCHPO5 294 : MCHPO5.JATTRI(1)=2 295 : MCHPO5.IFOPOI=IFOMOD 296 : SEGINI MSOUP5 297 : MCHPO5.IPCHP(1)=MSOUP5 298 : SEGINI MPOV5 299 : MSOUP5.NOCOMP(1)='SCAL' 300 : MSOUP5.IGEOC = IM1 301 : MSOUP5.IPOVAL = MPOV5 302 : 303 : C liquid temperature: Tl 304 : 305 : NC = 1 306 : SEGINI MCHPO6 307 : MCHPO6.JATTRI(1)=2 308 : MCHPO6.IFOPOI=IFOMOD 309 : SEGINI MSOUP6 310 : MCHPO6.IPCHP(1)=MSOUP6 311 : SEGINI MPOV6 312 : MSOUP6.NOCOMP(1)='SCAL' 313 : MSOUP6.IGEOC = IM1 314 : MSOUP6.IPOVAL = MPOV6 315 : 316 : C vapour density: rv 317 : 318 : NC = 1 319 : SEGINI MCHPO7 320 : MCHPO7.JATTRI(1)=2 321 : MCHPO7.IFOPOI=IFOMOD 322 : SEGINI MSOUP7 323 : MCHPO7.IPCHP(1)=MSOUP7 324 : SEGINI MPOV7 325 : MSOUP7.NOCOMP(1)='SCAL' 326 : MSOUP7.IGEOC = IM1 327 : MSOUP7.IPOVAL = MPOV7 328 : 329 : C liquid density: rl 330 : 331 : NC = 1 332 : SEGINI MCHPO8 333 : MCHPO8.JATTRI(1)=2 334 : MCHPO8.IFOPOI=IFOMOD 335 : SEGINI MSOUP8 336 : MCHPO8.IPCHP(1)=MSOUP8 337 : SEGINI MPOV8 338 : MSOUP8.NOCOMP(1)='SCAL' 339 : MSOUP8.IGEOC = IM1 340 : MSOUP8.IPOVAL = MPOV8 341 : 342 : C pressure correction term, pi 343 : 344 : NC = 1 345 : SEGINI MCHPO9 346 : MCHPO9.JATTRI(1)=2 347 : MCHPO9.IFOPOI=IFOMOD 348 : SEGINI MSOUP9 349 : MCHPO9.IPCHP(1)=MSOUP9 350 : SEGINI MPOV9 351 : MSOUP9.NOCOMP(1)='SCAL' 352 : MSOUP9.IGEOC = IM1 353 : MSOUP9.IPOVAL = MPOV9 354 : C 355 : C**** "Recapitulando" 356 : C 357 : C We activate MPOV1 - MPOV9 358 : C 359 : C**** "bucle" for the calculation of the primitives at 360 : C each centre 361 : 362 : 363 : DO NLCE = 1, N1 364 : 365 : oldav = MPOW7.VPOCHA(NLCE,1) 366 : oldTv = MPOW8.VPOCHA(NLCE,1) 367 : oldal = 1.d0 - oldav 368 : oldTl = MPOW9.VPOCHA(NLCE,1) 369 : 370 : IF(IDIM .EQ. 3)THEN 371 : if(oldav .LT. 0.01d0)then 372 : A = (gamv - 1.D0)*cpv*oldTv*MPOW1.VPOCHA(NLCE,1)/gamv 373 : else 374 : A = (gamv - 1.D0)*(MPOW5.VPOCHA(NLCE,1) - (MPOW3 375 : $ .VPOCHA(NLCE,1)**2 + MPOW3.VPOCHA(NLCE,2)**2 + 376 : & MPOW3.VPOCHA(NLCE,3)**2)/ 377 : & (2.d0*(MPOW1.VPOCHA(NLCE,1)+ 1.D-20))) 378 : endif 379 : if(oldav .GT. 0.99d0)then 380 : B = (gaml - 1.D0)*MPOW2.VPOCHA(NLCE,1)* 381 : & (cpl*MPOW9.VPOCHA(NLCE,1)/gaml + pil/1000.d0) 382 : else 383 : B = (gaml - 1.D0)*(MPOW6.VPOCHA(NLCE,1) - (MPOW4 384 : $ .VPOCHA(NLCE,1)**2 + MPOW4.VPOCHA(NLCE,2)**2)+ 385 : & MPOW4.VPOCHA(NLCE,3)**2/(2.D0 386 : $ *(MPOW2.VPOCHA(NLCE,1)+ 1.D-20))) 387 : endif 388 : ELSE 389 : if(oldav .LT. 0.01d0)then 390 : A = (gamv - 1.D0)*cpv*oldTv*MPOW1.VPOCHA(NLCE,1)/gamv 391 : else 392 : A = (gamv - 1.D0)*(MPOW5.VPOCHA(NLCE,1) - (MPOW3 393 : $ .VPOCHA(NLCE,1)**2 + MPOW3.VPOCHA(NLCE,2)**2)/(2.d0 394 : $ *(MPOW1.VPOCHA(NLCE,1)+ 1.D-20))) 395 : endif 396 : if(oldav .GT. 0.99d0)then 397 : B = (gaml - 1.D0)*MPOW2.VPOCHA(NLCE,1)* 398 : & (cpl*MPOW9.VPOCHA(NLCE,1)/gaml + pil/1000.d0) 399 : else 400 : 401 : B = (gaml - 1.D0)*(MPOW6.VPOCHA(NLCE,1) - (MPOW4 402 : $ .VPOCHA(NLCE,1)**2 + MPOW4.VPOCHA(NLCE,2)**2)/(2.D0 403 : $ *(MPOW2.VPOCHA(NLCE,1)+ 1.D-20))) 404 : endif 405 : ENDIF 406 : delta = (- A - B + gaml*pil)**2.D0 + 4.D0*gaml*pil*A 407 : 408 : p = 1.D0/2.D0*(A + B - gaml*pil + sqrt(delta)) 409 : 410 : alpha = A/p 410 : call FUNCF(alpha, FFV) 411 : uvx = FFV*MPOW3.VPOCHA(NLCE,1)/(MPOW1.VPOCHA(NLCE,1) + 1.D-20) 412 : uvy = FFV*MPOW3.VPOCHA(NLCE,2)/(MPOW1.VPOCHA(NLCE,1) + 1.D-20) 413 : IF(IDIM .EQ. 3)THEN 414 : uvz = FFV*MPOW3.VPOCHA(NLCE,3)/(MPOW1.VPOCHA(NLCE,1) + 1.D 415 : & -20) 416 : uv2 = uvx*uvx + uvy*uvy + uvz*uvz 417 : ELSE 418 : uv2 = uvx*uvx + uvy*uvy 419 : ENDIF 420 : 420 : call FUNCF((1.d0 - alpha),FFL) 421 : ulx = FFL*MPOW4.VPOCHA(NLCE,1)/(MPOW2.VPOCHA(NLCE,1) + 1.D-20) 422 : uly = FFL*MPOW4.VPOCHA(NLCE,2)/(MPOW2.VPOCHA(NLCE,1) + 1.D-20) 423 : IF(IDIM .EQ. 3)THEN 424 : ulz = FFL*MPOW4.VPOCHA(NLCE,3)/(MPOW2.VPOCHA(NLCE,1) + 1.D 425 : & -20) 426 : ul2 = ulx*ulx + uly*uly + ulz*ulz 427 : ELSE 428 : ul2 = ulx*ulx + uly*uly 429 : ENDIF 430 : 431 : if((alpha .GT. 0.999d0) .OR. (alpha .LT. 0.01d0))then 432 : C 433 : C Water faucet test 434 : C Tv = 323.15D0 435 : C Tl = 323.15D0 436 : C 437 : C Shock tube test 438 : C Tv = 308.15d0 439 : C Tl = 308.15d0 440 : C 441 : C Phase separation test or slosing of a water column 442 : Tv = 298.15D0 443 : Tl = 298.15D0 444 : C For the oscillating manometer should be 323 but to be practical I have left it in 298.15 445 : ev = cpv*Tv/gamv 446 : el = cpl*Tl/gaml*(p + gaml*pil)/(p + pil) 447 : else 448 : ev = MPOW5.VPOCHA(NLCE,1)/(MPOW1.VPOCHA(NLCE,1) + 1.D-20) 449 : & - uv2/2.D0 450 : el = MPOW6.VPOCHA(NLCE,1)/(MPOW2.VPOCHA(NLCE,1) + 1.D-20) 451 : & - ul2/2.D0 452 : Tv = gamv*ev/cpv 453 : Tl = gaml*el/(cpl*(1.D0 + pil*(gaml - 1.D0)/(p + pil))) 454 : endif 455 : rv = gamv*p/((gamv - 1.D0)*cpv*Tv) 456 : rl = gaml*(p + pil)/((gaml - 1.D0)*cpl*Tl) 457 : C 458 : C Bestion´s formula 459 : C 460 : rm = alpha*rv + (1.d0 - alpha)*rl 461 : rvm = rv + Cvm*rm 462 : rlm = rl + Cvm*rm 463 : 464 : pi = p - Cp*(alpha*rvm*(1.D0- alpha)*rlm)/ 465 : & (alpha*rlm + (1.D0- alpha)*rvm)* 466 : & (sqrt(uv2) - sqrt(ul2))**2 467 : 468 : MPOV1.VPOCHA(NLCE,1) = alpha 469 : MPOV2.VPOCHA(NLCE,1) = uvx 470 : MPOV2.VPOCHA(NLCE,2) = uvy 471 : MPOV3.VPOCHA(NLCE,1) = ulx 472 : MPOV3.VPOCHA(NLCE,2) = uly 473 : IF(IDIM .EQ. 3)THEN 474 : MPOV2.VPOCHA(NLCE,3) = uvz 475 : MPOV3.VPOCHA(NLCE,3) = ulz 476 : ENDIF 477 : MPOV4.VPOCHA(NLCE,1) = p 478 : MPOV5.VPOCHA(NLCE,1) = Tv 479 : MPOV6.VPOCHA(NLCE,1) = Tl 480 : MPOV7.VPOCHA(NLCE,1) = rv 481 : MPOV8.VPOCHA(NLCE,1) = rl 482 : MPOV9.VPOCHA(NLCE,1) = pi 483 : 484 : ENDDO 485 : C 486 : OCH1 = MCHPO1 487 : OCH2 = MCHPO2 488 : OCH3 = MCHPO3 489 : OCH4 = MCHPO4 490 : OCH5 = MCHPO5 491 : OCH6 = MCHPO6 492 : OCH7 = MCHPO7 493 : OCH8 = MCHPO8 494 : OCH9 = MCHPO9 495 : C 496 : SEGDES MPOV1 497 : SEGDES MPOV2 498 : SEGDES MPOV3 499 : SEGDES MPOV4 500 : SEGDES MPOV5 501 : SEGDES MPOV6 502 : SEGDES MPOV7 503 : SEGDES MPOV8 504 : SEGDES MPOV9 505 : 506 : SEGDES MPOW1 507 : SEGDES MPOW2 508 : SEGDES MPOW3 509 : SEGDES MPOW4 510 : SEGDES MPOW5 511 : SEGDES MPOW6 512 : SEGDES MPOW7 513 : SEGDES MPOW8 514 : SEGDES MPOW9 515 : 516 : SEGDES MSOUP1 517 : SEGDES MSOUP2 518 : SEGDES MSOUP3 519 : SEGDES MSOUP4 520 : SEGDES MSOUP5 521 : SEGDES MSOUP6 522 : SEGDES MSOUP7 523 : SEGDES MSOUP8 524 : SEGDES MSOUP9 525 : 526 : SEGDES MCHPO1 527 : SEGDES MCHPO2 528 : SEGDES MCHPO3 529 : SEGDES MCHPO4 530 : SEGDES MCHPO5 531 : SEGDES MCHPO6 532 : SEGDES MCHPO7 533 : SEGDES MCHPO8 534 : SEGDES MCHPO9 535 : 536 : RETURN 537 : END 538 : 539 : 540 : 541 : 542 : 543 : 544 : 545 : 546 : 547 : 548 : 549 : 550 : 551 :
© Cast3M 2003 - All rights reserved.
Disclaimer