1 : C CLI153 SOURCE MAGN 09/07/30 21:15:50 6458 2 : SUBROUTINE CLI153(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 : CLI153 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 : & RUZ.IZAFM, UXUZ.IZAFM, UYUZ.IZAFM, RETUZ.IZAFM, 44 : & UZR.IZAFM, UZUX.IZAFM, UZUY.IZAFM, UZRET.IZAFM, 45 : & UZUZ.IZAFM 46 : C 47 : C**** Variables de COOPTIO 48 : C 49 : C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR 50 : C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ 51 : C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES 52 : C & ,IECHO, IIMPI, IOSPI 53 : C & ,IDIM, IFICLE, IPREFI 54 : C & ,MCOORD 55 : C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE 56 : C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU 57 : C & ,NORINC,NORVAL,NORIND,NORVAD 58 : C & ,NUCROU, IPSAUV 59 : C 60 : INTEGER MELEMF,MELEMC,MELECB,INORM,ICHPVO,ICHPSU, IROC,IVITC,IPC 61 : & ,IGAMC,ICHLIM,ICEL,NFAC,IFAC,MELRES,IJACO 62 : & ,NGF,NGC,NLF,NLC,NLCB 63 : & ,ILIINC,ILIINP,IJAC 64 : & ,MP, NBEL, NBME, NBSOUS, NKID, NKMT, NMATRI, NP, NRIGE 65 : REAL*8 VOLU,SURF,RC,PC,UXC,UYC,GAMC,CNX,CNY,CTX,CTY,PF 66 : & ,CNZ,CTZ,CT2X,CT2Y,CT2Z,UZC 67 : & ,VECG(5),VECD(5),DFDUG(5,5),DFDUD(5,5) 68 : & ,NVEC(3),TVEC(3),TVEC2(3) 69 : & ,DRORO,DROUX,DROUY,DROP 70 : & ,DUXRO,DUXUX,DUXUY,DUXP 71 : & ,DUYRO,DUYUX,DUYUY,DUYP 72 : & ,DPRO,DPUX,DPUY,DPP 73 : & ,DROUZ,DUXUZ,DUYUZ,DPUZ 74 : & ,DUZRO,DUZUX,DUZUY,DUZUZ,DUZP 75 : & ,COEF1,COEF2,COEF3,GM1,COEF,COEF4 76 : CHARACTER*(8) TYPE 77 : C 78 : C 79 : C**** KRIPAD pour la correspondance global/local 80 : C 80 : CALL KRIPAD(MELEMC,MLEMC) 81 : C SEGINI MLEMC 81 : CALL KRIPAD(MELECB,MLEMCB) 82 : C SEGINI MLEMCB 82 : CALL KRIPAD(MELEMF,MLEMF) 83 : C SEGINI MLEMF 84 : C 85 : C**** CHPOINTs de la table DOMAINE 86 : C 86 : CALL LICHT(INORM,MPNORM,TYPE,ICEL) 86 : CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL) 86 : CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL) 87 : C 88 : C**** LICHT active les MPOVALs en *MOD 89 : C 90 : C SEGACT MPNORM*MOD 91 : C SEGACT MPOVSU*MOD 92 : C SEGACT MPOVOL*MOD 93 : C 94 : C 95 : C**** CHPOINTs des variables 96 : C 96 : CALL LICHT(IROC,MPRC,TYPE,ICEL) 96 : CALL LICHT(IVITC,MPVC,TYPE,ICEL) 96 : CALL LICHT(IPC,MPPC,TYPE,ICEL) 96 : CALL LICHT(IGAMC,MPGAMC,TYPE,ICEL) 96 : CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL) 97 : C 98 : C SEGACT *MOD 99 : C SEGACT *MOD 100 : C SEGACT *MOD 101 : C SEGACT *MOD 102 : C SEGACT *MOD 103 : C 104 : C 105 : C**** Boucle sur le face pour le calcul des invariants de 106 : C Riemann et du flux 107 : C 108 : SEGACT MELEFC 109 : NFAC=MELEFC.NUM(/2) 110 : C 111 : C**** Objet MATRIK 112 : C 113 : NRIGE = 7 114 : NMATRI = 1 115 : NKID = 9 116 : NKMT = 7 117 : C 118 : SEGINI MATRIK 119 : IJACO = MATRIK 120 : MATRIK.IRIGEL(1,1) = MELRES 121 : MATRIK.IRIGEL(2,1) = MELRES 122 : C 123 : C**** Matrice non symetrique 124 : C 125 : MATRIK.IRIGEL(7,1) = 2 126 : C 127 : NBME = 25 128 : NBSOUS = 1 129 : SEGINI IMATRI 130 : IF(IJAC.EQ.1)THEN 131 : MLMOTS=ILIINC 132 : ELSEIF(IJAC.EQ.2)THEN 133 : MLMOTS=ILIINP 134 : ENDIF 135 : SEGACT MLMOTS 136 : MATRIK.IRIGEL(4,1) = IMATRI 137 : C 138 : IMATRI.LISPRI(1) = MLMOTS.MOTS(1) 139 : IMATRI.LISPRI(2) = MLMOTS.MOTS(2) 140 : IMATRI.LISPRI(3) = MLMOTS.MOTS(3) 141 : IMATRI.LISPRI(4) = MLMOTS.MOTS(4) 142 : IMATRI.LISPRI(5) = MLMOTS.MOTS(5) 143 : C 144 : IMATRI.LISPRI(6) = MLMOTS.MOTS(1) 145 : IMATRI.LISPRI(7) = MLMOTS.MOTS(2) 146 : IMATRI.LISPRI(8) = MLMOTS.MOTS(3) 147 : IMATRI.LISPRI(9) = MLMOTS.MOTS(4) 148 : IMATRI.LISPRI(10) = MLMOTS.MOTS(5) 149 : C 150 : IMATRI.LISPRI(11) = MLMOTS.MOTS(1) 151 : IMATRI.LISPRI(12) = MLMOTS.MOTS(2) 152 : IMATRI.LISPRI(13) = MLMOTS.MOTS(3) 153 : IMATRI.LISPRI(14) = MLMOTS.MOTS(4) 154 : IMATRI.LISPRI(15) = MLMOTS.MOTS(5) 155 : C 156 : IMATRI.LISPRI(16) = MLMOTS.MOTS(1) 157 : IMATRI.LISPRI(17) = MLMOTS.MOTS(2) 158 : IMATRI.LISPRI(18) = MLMOTS.MOTS(3) 159 : IMATRI.LISPRI(19) = MLMOTS.MOTS(4) 160 : IMATRI.LISPRI(20) = MLMOTS.MOTS(5) 161 : C 162 : IMATRI.LISPRI(21) = MLMOTS.MOTS(1) 163 : IMATRI.LISPRI(22) = MLMOTS.MOTS(2) 164 : IMATRI.LISPRI(23) = MLMOTS.MOTS(3) 165 : IMATRI.LISPRI(24) = MLMOTS.MOTS(4) 166 : IMATRI.LISPRI(25) = MLMOTS.MOTS(5) 167 : C 168 : SEGDES MLMOTS 169 : MLMOTS=ILIINC 170 : SEGACT MLMOTS 171 : C 172 : C 173 : IMATRI.LISDUA(1) = MLMOTS.MOTS(1) 174 : IMATRI.LISDUA(2) = MLMOTS.MOTS(1) 175 : IMATRI.LISDUA(3) = MLMOTS.MOTS(1) 176 : IMATRI.LISDUA(4) = MLMOTS.MOTS(1) 177 : IMATRI.LISDUA(5) = MLMOTS.MOTS(1) 178 : C 179 : IMATRI.LISDUA(6) = MLMOTS.MOTS(2) 180 : IMATRI.LISDUA(7) = MLMOTS.MOTS(2) 181 : IMATRI.LISDUA(8) = MLMOTS.MOTS(2) 182 : IMATRI.LISDUA(9) = MLMOTS.MOTS(2) 183 : IMATRI.LISDUA(10) = MLMOTS.MOTS(2) 184 : C 185 : IMATRI.LISDUA(11) = MLMOTS.MOTS(3) 186 : IMATRI.LISDUA(12) = MLMOTS.MOTS(3) 187 : IMATRI.LISDUA(13) = MLMOTS.MOTS(3) 188 : IMATRI.LISDUA(14) = MLMOTS.MOTS(3) 189 : IMATRI.LISDUA(15) = MLMOTS.MOTS(3) 190 : C 191 : IMATRI.LISDUA(16) = MLMOTS.MOTS(4) 192 : IMATRI.LISDUA(17) = MLMOTS.MOTS(4) 193 : IMATRI.LISDUA(18) = MLMOTS.MOTS(4) 194 : IMATRI.LISDUA(19) = MLMOTS.MOTS(4) 195 : IMATRI.LISDUA(20) = MLMOTS.MOTS(4) 196 : C 197 : IMATRI.LISDUA(21) = MLMOTS.MOTS(5) 198 : IMATRI.LISDUA(22) = MLMOTS.MOTS(5) 199 : IMATRI.LISDUA(23) = MLMOTS.MOTS(5) 200 : IMATRI.LISDUA(24) = MLMOTS.MOTS(5) 201 : IMATRI.LISDUA(25) = MLMOTS.MOTS(5) 202 : C 203 : SEGDES MLMOTS 204 : NBEL = NFAC 205 : NBSOUS = 1 206 : NP = 1 207 : MP = 1 208 : SEGINI RR , RUX , RUY , RUZ ,RRET , 209 : & UXR , UXUX , UXUY , UXUZ ,UXRET , 210 : & UYR , UYUX , UYUY , UYUZ , UYRET , 211 : & UZR , UZUX , UZUY , UZUZ , UZRET , 212 : & RETR , RETUX , RETUY, RETUZ , RETRET 213 : C 214 : C**** Duale = IMATRI.LISDUA(1) = 'RN' 215 : C Primale = IMATRI.LISPRI(1) = 'RN' 216 : C -> IMATRI.LIZAFM(1,1) = RR 217 : C 218 : IMATRI.LIZAFM(1,1) = RR 219 : IMATRI.LIZAFM(1,2) = RUX 220 : IMATRI.LIZAFM(1,3) = RUY 221 : IMATRI.LIZAFM(1,4) = RUZ 222 : IMATRI.LIZAFM(1,5) = RRET 223 : C 224 : IMATRI.LIZAFM(1,6) = UXR 225 : IMATRI.LIZAFM(1,7) = UXUX 226 : IMATRI.LIZAFM(1,8) = UXUY 227 : IMATRI.LIZAFM(1,9) = UXUZ 228 : IMATRI.LIZAFM(1,10) = UXRET 229 : C 230 : IMATRI.LIZAFM(1,11) = UYR 231 : IMATRI.LIZAFM(1,12) = UYUX 232 : IMATRI.LIZAFM(1,13) = UYUY 233 : IMATRI.LIZAFM(1,14) = UYUZ 234 : IMATRI.LIZAFM(1,15) = UYRET 235 : C 236 : IMATRI.LIZAFM(1,16) = UZR 237 : IMATRI.LIZAFM(1,17) = UZUX 238 : IMATRI.LIZAFM(1,18) = UZUY 239 : IMATRI.LIZAFM(1,19) = UZUZ 240 : IMATRI.LIZAFM(1,20) = UZRET 241 : C 242 : IMATRI.LIZAFM(1,21) = RETR 243 : IMATRI.LIZAFM(1,22) = RETUX 244 : IMATRI.LIZAFM(1,23) = RETUY 245 : IMATRI.LIZAFM(1,24) = RETUZ 246 : IMATRI.LIZAFM(1,25) = RETRET 247 : C 248 : SEGDES MATRIK 249 : SEGDES IMATRI 250 : C 251 : C**** Fin definition MATRIK 252 : C 253 : DO IFAC=1,NFAC,1 254 : NGF=MELEFC.NUM(1,IFAC) 255 : NGC=MELEFC.NUM(2,IFAC) 256 : NLF=MLEMF.LECT(NGF) 257 : NLC=MLEMC.LECT(NGC) 258 : NLCB=MLEMCB.LECT(NGF) 259 : VOLU=MPVOL.VPOCHA(NLC,1) 260 : SURF=MPSURF.VPOCHA(NLF,1) 261 : C In CASTEM les normales sont sortantes 262 : CNX=MPNORM.VPOCHA(NLF,1) 263 : CNY=MPNORM.VPOCHA(NLF,2) 264 : CNZ=MPNORM.VPOCHA(NLF,3) 265 : CTX=MPNORM.VPOCHA(NLF,4) 266 : CTY=MPNORM.VPOCHA(NLF,5) 267 : CTZ=MPNORM.VPOCHA(NLF,6) 268 : CT2X=MPNORM.VPOCHA(NLF,7) 269 : CT2Y=MPNORM.VPOCHA(NLF,8) 270 : CT2Z=MPNORM.VPOCHA(NLF,9) 271 : C Variables au centre 272 : RC=MPRC.VPOCHA(NLC,1) 273 : PC=MPPC.VPOCHA(NLC,1) 274 : UXC=MPVC.VPOCHA(NLC,1) 275 : UYC=MPVC.VPOCHA(NLC,2) 276 : UZC=MPVC.VPOCHA(NLC,3) 277 : GAMC=MPGAMC.VPOCHA(NLC,1) 278 : C Variables à la face 279 : PF=MPLIM.VPOCHA(NLCB,1) 280 : C 281 : C******* Jacobian with respect to primitive variables 282 : C 283 : VECG(1)=RC 284 : VECG(2)=UXC 285 : VECG(3)=UYC 286 : VECG(4)=UZC 287 : VECG(5)=PC 288 : C 289 : VECD(1)=RC 290 : VECD(2)=UXC 291 : VECD(3)=UYC 292 : VECD(4)=UZC 293 : VECD(5)=PF 294 : C 295 : NVEC(1)=CNX 296 : NVEC(2)=CNY 297 : NVEC(3)=CNZ 298 : TVEC(1)=CTX 299 : TVEC(2)=CTY 300 : TVEC(3)=CTZ 301 : TVEC2(1)=CT2X 302 : TVEC2(2)=CT2Y 303 : TVEC2(3)=CT2Z 304 : C 304 : CALL CONJP4(DFDUG,DFDUD,VECG,VECD,NVEC,TVEC,TVEC2,GAMC) 305 : C 306 : COEF=-1*SURF/VOLU 307 : DRORO=(DFDUG(1,1)+DFDUD(1,1))*COEF 308 : DROUX=(DFDUG(1,2)+DFDUD(1,2))*COEF 309 : DROUY=(DFDUG(1,3)+DFDUD(1,3))*COEF 310 : DROUZ=(DFDUG(1,4)+DFDUD(1,4))*COEF 311 : DROP=DFDUG(1,5)*COEF 312 : C 313 : DPRO=(DFDUG(5,1)+DFDUD(5,1))*COEF 314 : DPUX=(DFDUG(5,2)+DFDUD(5,2))*COEF 315 : DPUY=(DFDUG(5,3)+DFDUD(5,3))*COEF 316 : DPUZ=(DFDUG(5,4)+DFDUD(5,4))*COEF 317 : DPP=DFDUG(5,5)*COEF 318 : C 319 : DUXRO=(DFDUG(2,1)+DFDUD(2,1))*COEF 320 : DUXUX=(DFDUG(2,2)+DFDUD(2,2))*COEF 321 : DUXUY=(DFDUG(2,3)+DFDUD(2,3))*COEF 322 : DUXUZ=(DFDUG(2,4)+DFDUD(2,4))*COEF 323 : DUXP =DFDUG(2,5)*COEF 324 : C 325 : DUYRO=(DFDUG(3,1)+DFDUD(3,1))*COEF 326 : DUYUX=(DFDUG(3,2)+DFDUD(3,2))*COEF 327 : DUYUY=(DFDUG(3,3)+DFDUD(3,3))*COEF 328 : DUYUZ=(DFDUG(3,4)+DFDUD(3,4))*COEF 329 : DUYP =DFDUG(3,5)*COEF 330 : C 331 : DUZRO=(DFDUG(4,1)+DFDUD(4,1))*COEF 332 : DUZUX=(DFDUG(4,2)+DFDUD(4,2))*COEF 333 : DUZUY=(DFDUG(4,3)+DFDUD(4,3))*COEF 334 : DUZUZ=(DFDUG(4,4)+DFDUD(4,4))*COEF 335 : DUZP =DFDUG(4,5)*COEF 336 : C 337 : C******* Jacobian with respect to conservative variables 338 : C 339 : IF(IJAC.EQ.1)THEN 340 : GM1=GAMC-1.0D0 341 : C 342 : COEF1=-1.0D0*UXC/RC 343 : COEF2=-1.0D0*UYC/RC 344 : COEF3=-1.0D0*UZC/RC 345 : COEF4=0.5D0*GM1*((UXC*UXC)+(UYC*UYC)+(UZC*UZC)) 346 : C 347 : RR.AM(IFAC,1,1)=DRORO+(DROUX*COEF1)+(DROUY*COEF2)+(DROUZ 348 : $ *COEF3)+(DROP*COEF4) 349 : RUX.AM(IFAC,1,1)=(DROUX/RC)-((UXC*GM1)*DROP) 350 : RUY.AM(IFAC,1,1)=(DROUY/RC)-((UYC*GM1)*DROP) 351 : RUZ.AM(IFAC,1,1)=(DROUZ/RC)-((UZC*GM1)*DROP) 352 : RRET.AM(IFAC,1,1)=GM1*DROP 353 : C 354 : UXR.AM(IFAC,1,1)=DUXRO+(DUXUX*COEF1)+(DUXUY*COEF2)+(DUXUZ 355 : $ *COEF3)+(DUXP*COEF4) 356 : UXUX.AM(IFAC,1,1)=(DUXUX/RC)-((UXC*GM1)*DUXP) 357 : UXUY.AM(IFAC,1,1)=(DUXUY/RC)-((UYC*GM1)*DUXP) 358 : UXUZ.AM(IFAC,1,1)=(DUXUZ/RC)-((UZC*GM1)*DUXP) 359 : UXRET.AM(IFAC,1,1)=GM1*DUXP 360 : C 361 : UYR.AM(IFAC,1,1)=DUYRO+(DUYUX*COEF1)+(DUYUY*COEF2)+(DUYUZ 362 : $ *COEF3)+(DUYP*COEF4) 363 : UYUX.AM(IFAC,1,1)=(DUYUX/RC)-((UXC*GM1)*DUYP) 364 : UYUY.AM(IFAC,1,1)=(DUYUY/RC)-((UYC*GM1)*DUYP) 365 : UYUZ.AM(IFAC,1,1)=(DUYUZ/RC)-((UZC*GM1)*DUYP) 366 : UYRET.AM(IFAC,1,1)=GM1*DUYP 367 : C 368 : UZR.AM(IFAC,1,1)=DUZRO+(DUZUX*COEF1)+(DUZUY*COEF2)+(DUZUZ 369 : $ *COEF3)+(DUZP*COEF4) 370 : UZUX.AM(IFAC,1,1)=(DUZUX/RC)-((UXC*GM1)*DUZP) 371 : UZUY.AM(IFAC,1,1)=(DUZUY/RC)-((UYC*GM1)*DUZP) 372 : UZUZ.AM(IFAC,1,1)=(DUZUZ/RC)-((UZC*GM1)*DUZP) 373 : UZRET.AM(IFAC,1,1)=GM1*DUZP 374 : C 375 : RETR.AM(IFAC,1,1)=DPRO+(DPUX*COEF1)+(DPUY*COEF2)+(DPUZ 376 : $ *COEF3)+(DPP*COEF4) 377 : RETUX.AM(IFAC,1,1)=(DPUX/RC)-((UXC*GM1)*DPP) 378 : RETUY.AM(IFAC,1,1)=(DPUY/RC)-((UYC*GM1)*DPP) 379 : RETUZ.AM(IFAC,1,1)=(DPUZ/RC)-((UZC*GM1)*DPP) 380 : RETRET.AM(IFAC,1,1)=GM1*DPP 381 : C 382 : ELSEIF(IJAC.EQ.2)THEN 383 : RR.AM(IFAC,1,1)=DRORO 384 : RUX.AM(IFAC,1,1)=DROUX 385 : RUY.AM(IFAC,1,1)=DROUY 386 : RUZ.AM(IFAC,1,1)=DROUZ 387 : RRET.AM(IFAC,1,1)=DROP 388 : C 389 : UXR.AM(IFAC,1,1)=DUXRO 390 : UXUX.AM(IFAC,1,1)=DUXUX 391 : UXUY.AM(IFAC,1,1)=DUXUY 392 : UXUZ.AM(IFAC,1,1)=DUXUZ 393 : UXRET.AM(IFAC,1,1)=DUXP 394 : C 395 : UYR.AM(IFAC,1,1)=DUYRO 396 : UYUX.AM(IFAC,1,1)=DUYUX 397 : UYUY.AM(IFAC,1,1)=DUYUY 398 : UYUZ.AM(IFAC,1,1)=DUYUZ 399 : UYRET.AM(IFAC,1,1)=DUYP 400 : C 401 : UZR.AM(IFAC,1,1)=DUZRO 402 : UZUX.AM(IFAC,1,1)=DUZUX 403 : UZUY.AM(IFAC,1,1)=DUZUY 404 : UZUZ.AM(IFAC,1,1)=DUZUZ 405 : UZRET.AM(IFAC,1,1)=DUZP 406 : C 407 : RETR.AM(IFAC,1,1)=DPRO 408 : RETUX.AM(IFAC,1,1)=DPUX 409 : RETUY.AM(IFAC,1,1)=DPUY 410 : RETUZ.AM(IFAC,1,1)=DPUZ 411 : RETRET.AM(IFAC,1,1)=DPP 412 : ENDIF 413 : ENDDO 414 : C 415 : SEGDES MELEFC 416 : C 417 : SEGSUP MLEMC 418 : SEGSUP MLEMCB 419 : SEGSUP MLEMF 420 : C 421 : SEGDES MPNORM 422 : SEGDES MPVOL 423 : SEGDES MPSURF 424 : SEGDES MPRC 425 : SEGDES MPPC 426 : SEGDES MPVC 427 : SEGDES MPGAMC 428 : SEGDES MPLIM 429 : C 430 : SEGDES RR , RUX , RUY , RUZ ,RRET , 431 : & UXR , UXUX , UXUY , UXUZ ,UXRET , 432 : & UYR , UYUX , UYUY , UYUZ , UYRET , 433 : & UZR , UZUX , UZUY , UZUZ , UZRET , 434 : & RETR , RETUX , RETUY, RETUZ , RETRET 435 : C 436 : 9999 CONTINUE 437 : RETURN 438 : END 439 : 440 : 441 : 442 : 443 : 444 : 445 : 446 : 447 :
© Cast3M 2003 - All rights reserved.
Disclaimer