1 : C CLI122 SOURCE MAGN 09/07/30 21:15:38 6458 2 : SUBROUTINE CLI122(MELEMF,MELEMC,MELECB,MELEFC,MELRES,INORM,ICHPVO, 3 : & ICHPSU,IROC,IVITC,IPC,IGAMC,ICHLIM,ILIINC,ILIINP,IJAC,IJACO) 4 : C************************************************************************ 5 : C 6 : C PROJET : CASTEM 2000 7 : C 8 : C NOM : CLI122 9 : C 10 : C DESCRIPTION : Subroutine appellée par CLIM11 11 : C 12 : C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec estensions CISI) 13 : C 14 : C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/LTMF 15 : C 16 : C************************************************************************ 17 : C 18 : C APPELES (Calcul) : 19 : C 20 : C************************************************************************ 21 : C 22 : C HISTORIQUE (Anomalies et modifications éventuelles) 23 : C 24 : C HISTORIQUE : 25 : C 26 : C************************************************************************ 27 : C 28 : IMPLICIT INTEGER(I-N) 29 : -INC CCOPTIO 30 : -INC SMLMOTS 31 : -INC SMELEME 32 : POINTEUR MELEFC.MELEME 33 : -INC SMLENTI 34 : POINTEUR MLEMC.MLENTI, MLEMCB.MLENTI,MLEMF.MLENTI 35 : -INC SMCHPOI 36 : POINTEUR MPNORM.MPOVAL, MPVOL.MPOVAL, MPSURF.MPOVAL, MPRC.MPOVAL, 37 : & MPVC.MPOVAL, MPPC.MPOVAL, MPGAMC.MPOVAL, MPLIM.MPOVAL 38 : -INC SMMATRIK 39 : POINTEUR RR.IZAFM, RUX.IZAFM, RUY.IZAFM, RRET.IZAFM, 40 : & UXR.IZAFM, UXUX.IZAFM, UXUY.IZAFM, UXRET.IZAFM, 41 : & UYR.IZAFM, UYUX.IZAFM, UYUY.IZAFM, UYRET.IZAFM, 42 : & RETR.IZAFM, RETUX.IZAFM, RETUY.IZAFM, RETRET.IZAFM 43 : C 44 : C**** Variables de COOPTIO 45 : C 46 : C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR 47 : C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ 48 : C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES 49 : C & ,IECHO, IIMPI, IOSPI 50 : C & ,IDIM, IFICLE, IPREFI 51 : C & ,MCOORD 52 : C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE 53 : C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU 54 : C & ,NORINC,NORVAL,NORIND,NORVAD 55 : C & ,NUCROU, IPSAUV 56 : C 57 : INTEGER MELEMF,MELEMC,MELECB,INORM,ICHPVO,ICHPSU, IROC,IVITC,IPC 58 : & ,IGAMC,ICHLIM,ICEL,NFAC,IFAC,MELRES,IJACO 59 : & ,NGF,NGC,NLF,NLC,NLCB 60 : & ,ILIINC,ILIINP,IJAC 61 : & ,MP, NBEL, NBME, NBSOUS, NKID, NKMT, NMATRI, NP, NRIGE 62 : REAL*8 VOLU,SURF,RC,PC,UXC,UYC,GAMC,CNX,CNY,CTX,CTY 63 : & ,RF,PF,UXF,UYF 64 : & ,UNC,UNF,ASONC,ASONF,ASON 65 : * & ,UTF,SF 66 : & ,UTC,SC 67 : & ,GM1,USGM1,DSGM1,G1,G3,ASON2,S,UT,UN,RHO,P,UX,UY 68 : & ,DUNDG3,DASDG3,DRHDG3,DPDG3,DRHDS,DPDS 69 : * & ,CELL,EPS,CACCA 70 : & ,DUXDG3,DUYDG3,DUXDUT,DUYDUT 71 : & ,DFRDG3,DFRDS,DFMXG3,DFMXS,DFMXUT,DFMYG3,DFMYS,DFMYUT 72 : & ,DFEDG3,DFEDS,DFEDUT 73 : & ,DG3DR,DG3DP,DG3DUX,DG3DUY,DSDR,DSDP,DUTDUX,DUTDUY,COEF 74 : & ,DRORO,DROUX,DROUY,DROP 75 : & ,DUXRO,DUXUX,DUXUY,DUXP 76 : & ,DUYRO,DUYUX,DUYUY,DUYP 77 : & ,DPRO,DPUX,DPUY,DPP 78 : & ,COEF1,COEF2,COEF3 79 : CHARACTER*(8) TYPE 80 : C 81 : C 82 : C**** KRIPAD pour la correspondance global/local 83 : C 83 : CALL KRIPAD(MELEMC,MLEMC) 84 : C SEGINI MLEMC 84 : CALL KRIPAD(MELECB,MLEMCB) 85 : C SEGINI MLEMCB 85 : CALL KRIPAD(MELEMF,MLEMF) 86 : C SEGINI MLEMF 87 : C 88 : C**** CHPOINTs de la table DOMAINE 89 : C 89 : CALL LICHT(INORM,MPNORM,TYPE,ICEL) 89 : CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL) 89 : CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL) 90 : C 91 : C**** LICHT active les MPOVALs en *MOD 92 : C 93 : C SEGACT MPNORM*MOD 94 : C SEGACT MPOVSU*MOD 95 : C SEGACT MPOVOL*MOD 96 : C 97 : C 98 : C**** CHPOINTs des variables 99 : C 99 : CALL LICHT(IROC,MPRC,TYPE,ICEL) 99 : CALL LICHT(IVITC,MPVC,TYPE,ICEL) 99 : CALL LICHT(IPC,MPPC,TYPE,ICEL) 99 : CALL LICHT(IGAMC,MPGAMC,TYPE,ICEL) 99 : CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL) 100 : C 101 : C SEGACT *MOD 102 : C SEGACT *MOD 103 : C SEGACT *MOD 104 : C SEGACT *MOD 105 : C SEGACT *MOD 106 : C 107 : C 108 : C**** Boucle sur le face pour le calcul des invariants de 109 : C Riemann et du flux 110 : C 111 : SEGACT MELEFC 112 : NFAC=MELEFC.NUM(/2) 113 : C 114 : C**** Objet MATRIK 115 : C 116 : NRIGE = 7 117 : NMATRI = 1 118 : NKID = 9 119 : NKMT = 7 120 : C 121 : SEGINI MATRIK 122 : IJACO = MATRIK 123 : MATRIK.IRIGEL(1,1) = MELRES 124 : MATRIK.IRIGEL(2,1) = MELRES 125 : C 126 : C**** Matrice non symetrique 127 : C 128 : MATRIK.IRIGEL(7,1) = 2 129 : C 130 : NBME = 16 131 : NBSOUS = 1 132 : SEGINI IMATRI 133 : IF(IJAC.EQ.1)THEN 134 : MLMOTS=ILIINC 135 : ELSEIF(IJAC.EQ.2)THEN 136 : MLMOTS=ILIINP 137 : ENDIF 138 : SEGACT MLMOTS 139 : MATRIK.IRIGEL(4,1) = IMATRI 140 : C 141 : IMATRI.LISPRI(1) = MLMOTS.MOTS(1) 142 : IMATRI.LISPRI(2) = MLMOTS.MOTS(2) 143 : IMATRI.LISPRI(3) = MLMOTS.MOTS(3) 144 : IMATRI.LISPRI(4) = MLMOTS.MOTS(4) 145 : IMATRI.LISPRI(5) = MLMOTS.MOTS(1) 146 : IMATRI.LISPRI(6) = MLMOTS.MOTS(2) 147 : IMATRI.LISPRI(7) = MLMOTS.MOTS(3) 148 : IMATRI.LISPRI(8) = MLMOTS.MOTS(4) 149 : IMATRI.LISPRI(9) = MLMOTS.MOTS(1) 150 : IMATRI.LISPRI(10) = MLMOTS.MOTS(2) 151 : IMATRI.LISPRI(11) = MLMOTS.MOTS(3) 152 : IMATRI.LISPRI(12) = MLMOTS.MOTS(4) 153 : IMATRI.LISPRI(13) = MLMOTS.MOTS(1) 154 : IMATRI.LISPRI(14) = MLMOTS.MOTS(2) 155 : IMATRI.LISPRI(15) = MLMOTS.MOTS(3) 156 : IMATRI.LISPRI(16) = MLMOTS.MOTS(4) 157 : C 158 : SEGDES MLMOTS 159 : MLMOTS=ILIINC 160 : SEGACT MLMOTS 161 : C 162 : IMATRI.LISDUA(1) = MLMOTS.MOTS(1) 163 : IMATRI.LISDUA(2) = MLMOTS.MOTS(1) 164 : IMATRI.LISDUA(3) = MLMOTS.MOTS(1) 165 : IMATRI.LISDUA(4) = MLMOTS.MOTS(1) 166 : IMATRI.LISDUA(5) = MLMOTS.MOTS(2) 167 : IMATRI.LISDUA(6) = MLMOTS.MOTS(2) 168 : IMATRI.LISDUA(7) = MLMOTS.MOTS(2) 169 : IMATRI.LISDUA(8) = MLMOTS.MOTS(2) 170 : IMATRI.LISDUA(9) = MLMOTS.MOTS(3) 171 : IMATRI.LISDUA(10) = MLMOTS.MOTS(3) 172 : IMATRI.LISDUA(11) = MLMOTS.MOTS(3) 173 : IMATRI.LISDUA(12) = MLMOTS.MOTS(3) 174 : IMATRI.LISDUA(13) = MLMOTS.MOTS(4) 175 : IMATRI.LISDUA(14) = MLMOTS.MOTS(4) 176 : IMATRI.LISDUA(15) = MLMOTS.MOTS(4) 177 : IMATRI.LISDUA(16) = MLMOTS.MOTS(4) 178 : C 179 : SEGDES MLMOTS 180 : NBEL = NFAC 181 : NBSOUS = 1 182 : NP = 1 183 : MP = 1 184 : SEGINI RR , RUX , RUY , RRET , 185 : & UXR , UXUX , UXUY , UXRET , 186 : & UYR , UYUX , UYUY , UYRET , 187 : & RETR , RETUX , RETUY , RETRET 188 : C 189 : C**** Duale = IMATRI.LISDUA(1) = 'RN' 190 : C Primale = IMATRI.LISPRI(1) = 'RN' 191 : C -> IMATRI.LIZAFM(1,1) = RR 192 : C 193 : IMATRI.LIZAFM(1,1) = RR 194 : IMATRI.LIZAFM(1,2) = RUX 195 : IMATRI.LIZAFM(1,3) = RUY 196 : IMATRI.LIZAFM(1,4) = RRET 197 : IMATRI.LIZAFM(1,5) = UXR 198 : IMATRI.LIZAFM(1,6) = UXUX 199 : IMATRI.LIZAFM(1,7) = UXUY 200 : IMATRI.LIZAFM(1,8) = UXRET 201 : IMATRI.LIZAFM(1,9) = UYR 202 : IMATRI.LIZAFM(1,10) = UYUX 203 : IMATRI.LIZAFM(1,11) = UYUY 204 : IMATRI.LIZAFM(1,12) = UYRET 205 : IMATRI.LIZAFM(1,13) = RETR 206 : IMATRI.LIZAFM(1,14) = RETUX 207 : IMATRI.LIZAFM(1,15) = RETUY 208 : IMATRI.LIZAFM(1,16) = RETRET 209 : C 210 : SEGDES MATRIK 211 : SEGDES IMATRI 212 : C 213 : C**** Fin definition MATRIK 214 : C 215 : DO IFAC=1,NFAC,1 216 : NGF=MELEFC.NUM(1,IFAC) 217 : NGC=MELEFC.NUM(2,IFAC) 218 : NLF=MLEMF.LECT(NGF) 219 : NLC=MLEMC.LECT(NGC) 220 : NLCB=MLEMCB.LECT(NGF) 221 : VOLU=MPVOL.VPOCHA(NLC,1) 222 : SURF=MPSURF.VPOCHA(NLF,1) 223 : C In CASTEM les normales sont sortantes 224 : CNX=MPNORM.VPOCHA(NLF,1) 225 : CNY=MPNORM.VPOCHA(NLF,2) 226 : CTX=-1.0D0*CNY 227 : CTY=CNX 228 : C Variables au centre 229 : RC=MPRC.VPOCHA(NLC,1) 230 : PC=MPPC.VPOCHA(NLC,1) 231 : UXC=MPVC.VPOCHA(NLC,1) 232 : UYC=MPVC.VPOCHA(NLC,2) 233 : GAMC=MPGAMC.VPOCHA(NLC,1) 234 : C Variables à la face 235 : RF=MPLIM.VPOCHA(NLCB,1) 236 : UXF=MPLIM.VPOCHA(NLCB,2) 237 : UYF=MPLIM.VPOCHA(NLCB,3) 238 : PF=MPLIM.VPOCHA(NLCB,IDIM+2) 239 : C 240 : C******* On calcule UN, UT ASON, S 241 : C 242 : UNC=(UXC*CNX)+(UYC*CNY) 243 : UTC=(UXC*CTX)+(UYC*CTY) 244 : UNF=(UXF*CNX)+(UYF*CNY) 245 : C UTF=(UXF*CTX)+(UYF*CTY) 246 : C 247 : ASONC=(GAMC*PC/RC)**0.5D0 248 : ASONF=(GAMC*PF/RF)**0.5D0 249 : C 250 : SC=PC/(RC**GAMC) 251 : C SF=PF/(RF**GAMC) 252 : C 253 : C******* Densite, vitesse, pression sur le bord 254 : C 255 : GM1=(GAMC-1.0D0) 256 : USGM1=1.0D0/GM1 257 : DSGM1=2.0D0*USGM1 258 : G1=UNF-(DSGM1*ASONF) 259 : G3=UNC+(DSGM1*ASONC) 260 : UN=0.5D0*(G1+G3) 261 : ASON=(0.5D0*(G3-G1)) 262 : ASON=ASON/DSGM1 263 : ASON2=ASON*ASON 264 : S=SC 265 : UT=UTC 266 : RHO=ASON2/(GAMC*S) 267 : RHO=RHO**USGM1 268 : P=RHO*ASON2/GAMC 269 : UX=(UN*CNX)+(UT*CTX) 270 : UY=(UN*CNY)+(UT*CTY) 271 : C 272 : C******* Derivatives 273 : C 274 : DUNDG3=0.5D0 275 : DASDG3=0.5D0/DSGM1 276 : DRHDG3=GAMC*S 277 : DRHDG3=1.0D0/DRHDG3 278 : DRHDG3=DRHDG3**USGM1 279 : DRHDG3=DRHDG3*DSGM1 280 : DRHDG3=DRHDG3*(ASON**((3.0D0-GAMC)/(GAMC-1.0D0))) 281 : DRHDG3=DRHDG3*DASDG3 282 : DRHDS=-1.0D0*USGM1*RHO/S 283 : DPDG3=((ASON2/GAMC)*DRHDG3)+(((2*ASON*RHO)/GAMC)*DASDG3) 284 : DPDS=(ASON2/GAMC)*DRHDS 285 : CC 286 : CC******* Test 287 : CC 288 : CC ...DG3 289 : C EPS=1.0D-6 290 : C G3=G3*(1+EPS) 291 : C CELL=UN 292 : C UN=0.5D0*(G1+G3) 293 : C write(*,*) DUNDG3, ((UN - CELL)/(EPS*G3)) 294 : C CELL=ASON 295 : C ASON=(0.5D0*(G3-G1)) 296 : C ASON=ASON/DSGM1 297 : C write(*,*) DASDG3, ((ASON - CELL)/(EPS*G3)) 298 : C ASON2=ASON*ASON 299 : C S=SC 300 : C UT=UTC 301 : C CELL=RHO 302 : C RHO=ASON2/(GAMC*S) 303 : C RHO=RHO**USGM1 304 : C write(*,*) DRHDG3, ((RHO - CELL)/(EPS*G3)) 305 : C CELL=P 306 : C P=RHO*ASON2/GAMC 307 : C write(*,*) DPDG3, ((P - CELL)/(EPS*G3)) 308 : CC 309 : CC******* Test 310 : CC 311 : CC ...DS 312 : C EPS=1.0D-6 313 : C S=S*(1+EPS) 314 : C CELL=RHO 315 : C RHO=ASON2/(GAMC*S) 316 : C RHO=RHO**USGM1 317 : C write(*,*) DRHDS, ((RHO - CELL)/(EPS*S)) 318 : C CELL=P 319 : C P=RHO*ASON2/GAMC 320 : C write(*,*) DPDS, ((P - CELL)/(EPS*S)) 321 : CC 322 : CC************************************************************ 323 : CC 324 : DUXDG3=DUNDG3*CNX 325 : DUYDG3=DUNDG3*CNY 326 : DUXDUT=CTX 327 : DUYDUT=CTY 328 : C 329 : DFRDG3=(DRHDG3*UN)+(RHO*DUNDG3) 330 : DFRDS=(DRHDS*UN) 331 : DFMXG3=(DRHDG3*UN*UX)+(RHO*DUNDG3*UX)+ 332 : & (RHO*UN*DUXDG3)+(DPDG3*CNX) 333 : DFMXS=(DRHDS*UN*UX)+(DPDS*CNX) 334 : DFMXUT=(RHO*UN*DUXDUT) 335 : DFMYG3=(DRHDG3*UN*UY)+(RHO*DUNDG3*UY)+ 336 : & (RHO*UN*DUYDG3)+(DPDG3*CNY) 337 : DFMYS=(DRHDS*UN*UY)+(DPDS*CNY) 338 : DFMYUT=(RHO*UN*DUYDUT) 339 : DFEDG3=(DUNDG3*GAMC*USGM1*P) + (UN*GAMC*USGM1*DPDG3) + 340 : & (0.5D0*DRHDG3*UN*((UN*UN)+(UT*UT))) + 341 : & (0.5D0*RHO*DUNDG3*((UN*UN)+(UT*UT))) + 342 : & (RHO*UN*UN*DUNDG3) 343 : DFEDS=(UN*GAMC*USGM1*DPDS) + 344 : & (0.5D0*DRHDS*UN*((UN*UN)+(UT*UT))) 345 : DFEDUT=RHO*UN*UT 346 : C 347 : C******* Jacobian with respect to primitive variables 348 : C 349 : DG3DR=-1*USGM1*ASONC/RC 350 : DG3DP=USGM1*ASONC/PC 351 : DG3DUX=CNX 352 : DG3DUY=CNY 353 : DSDR=-1*GAMC*S/RC 354 : DSDP=S/PC 355 : DUTDUX=CTX 356 : DUTDUY=CTY 357 : C 358 : COEF=-1*SURF/VOLU 359 : C 360 : DRORO=((DFRDG3*DG3DR)+(DFRDS*DSDR))*COEF 361 : DROUX=(DFRDG3*DG3DUX)*COEF 362 : DROUY=(DFRDG3*DG3DUY)*COEF 363 : DROP=((DFRDG3*DG3DP)+(DFRDS*DSDP))*COEF 364 : C 365 : DUXRO=((DFMXG3*DG3DR)+(DFMXS*DSDR))*COEF 366 : DUXUX=((DFMXG3*DG3DUX)+(DFMXUT*DUTDUX))*COEF 367 : DUXUY=((DFMXG3*DG3DUY)+(DFMXUT*DUTDUY))*COEF 368 : DUXP=((DFMXG3*DG3DP)+(DFMXS*DSDP))*COEF 369 : C 370 : DUYRO=((DFMYG3*DG3DR)+(DFMYS*DSDR))*COEF 371 : DUYUX=((DFMYG3*DG3DUX)+(DFMYUT*DUTDUX))*COEF 372 : DUYUY=((DFMYG3*DG3DUY)+(DFMYUT*DUTDUY))*COEF 373 : DUYP=((DFMYG3*DG3DP)+(DFMYS*DSDP))*COEF 374 : CC 375 : DPRO=((DFEDG3*DG3DR)+(DFEDS*DSDR))*COEF 376 : DPUX=((DFEDG3*DG3DUX)+(DFEDUT*DUTDUX))*COEF 377 : DPUY=((DFEDG3*DG3DUY)+(DFEDUT*DUTDUY))*COEF 378 : DPP=((DFEDG3*DG3DP)+(DFEDS*DSDP))*COEF 379 : CC 380 : CC******* Jacobian with respect to conservative variables 381 : CC 382 : IF(IJAC.EQ.1)THEN 383 : 384 : COEF1=-1.0D0*UXC/RC 385 : COEF2=-1.0D0*UYC/RC 386 : COEF3=0.5D0*GM1*((UXC*UXC)+(UYC*UYC)) 387 : C 388 : RR.AM(IFAC,1,1)=DRORO+(DROUX*COEF1)+(DROUY*COEF2)+(DROP 389 : $ *COEF3) 390 : RUX.AM(IFAC,1,1)=(DROUX/RC)-((UXC*GM1)*DROP) 391 : RUY.AM(IFAC,1,1)=(DROUY/RC)-((UYC*GM1)*DROP) 392 : RRET.AM(IFAC,1,1)=GM1*DROP 393 : C 394 : UXR.AM(IFAC,1,1)=DUXRO+(DUXUX*COEF1)+(DUXUY*COEF2)+(DUXP 395 : $ *COEF3) 396 : UXUX.AM(IFAC,1,1)=(DUXUX/RC)-((UXC*GM1)*DUXP) 397 : UXUY.AM(IFAC,1,1)=(DUXUY/RC)-((UYC*GM1)*DUXP) 398 : UXRET.AM(IFAC,1,1)=GM1*DUXP 399 : C 400 : UYR.AM(IFAC,1,1)=DUYRO+(DUYUX*COEF1)+(DUYUY*COEF2)+(DUYP 401 : $ *COEF3) 402 : UYUX.AM(IFAC,1,1)=(DUYUX/RC)-((UXC*GM1)*DUYP) 403 : UYUY.AM(IFAC,1,1)=(DUYUY/RC)-((UYC*GM1)*DUYP) 404 : UYRET.AM(IFAC,1,1)=GM1*DUYP 405 : C 406 : RETR.AM(IFAC,1,1)=DPRO+(DPUX*COEF1)+(DPUY*COEF2)+(DPP 407 : $ *COEF3) 408 : RETUX.AM(IFAC,1,1)=(DPUX/RC)-((UXC*GM1)*DPP) 409 : RETUY.AM(IFAC,1,1)=(DPUY/RC)-((UYC*GM1)*DPP) 410 : RETRET.AM(IFAC,1,1)=GM1*DPP 411 : C 412 : ELSEIF(IJAC.EQ.2)THEN 413 : RR.AM(IFAC,1,1)=DRORO 414 : RUX.AM(IFAC,1,1)=DROUX 415 : RUY.AM(IFAC,1,1)=DROUY 416 : RRET.AM(IFAC,1,1)=DROP 417 : C 418 : UXR.AM(IFAC,1,1)=DUXRO 419 : UXUX.AM(IFAC,1,1)=DUXUX 420 : UXUY.AM(IFAC,1,1)=DUXUY 421 : UXRET.AM(IFAC,1,1)=DUXP 422 : C 423 : UYR.AM(IFAC,1,1)=DUYRO 424 : UYUX.AM(IFAC,1,1)=DUYUX 425 : UYUY.AM(IFAC,1,1)=DUYUY 426 : UYRET.AM(IFAC,1,1)=DUYP 427 : C 428 : RETR.AM(IFAC,1,1)=DPRO 429 : RETUX.AM(IFAC,1,1)=DPUX 430 : RETUY.AM(IFAC,1,1)=DPUY 431 : RETRET.AM(IFAC,1,1)=DPP 432 : ENDIF 433 : ENDDO 434 : C 435 : SEGDES MELEFC 436 : C 437 : SEGSUP MLEMC 438 : SEGSUP MLEMCB 439 : SEGSUP MLEMF 440 : C 441 : SEGDES MPNORM 442 : SEGDES MPVOL 443 : SEGDES MPSURF 444 : SEGDES MPRC 445 : SEGDES MPPC 446 : SEGDES MPVC 447 : SEGDES MPGAMC 448 : SEGDES MPLIM 449 : C 450 : SEGDES RR , RUX , RUY , RRET , 451 : & UXR , UXUX , UXUY , UXRET , 452 : & UYR , UYUX , UYUY , UYRET , 453 : & RETR , RETUX , RETUY , RETRET 454 : C 455 : 9999 CONTINUE 456 : RETURN 457 : END 458 : 459 : 460 : 461 : 462 : 463 : 464 : 465 : 466 :
© Cast3M 2003 - All rights reserved.
Disclaimer