Télécharger cli122.eso

Retour à la liste

Numérotation des lignes :

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

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