Télécharger cli153.eso

Retour à la liste

Numérotation des lignes :

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

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