Télécharger cli282.eso

Retour à la liste

Numérotation des lignes :

  1. C CLI282 SOURCE PV 16/11/17 21:58:48 9180
  2. SUBROUTINE CLI282(NSP,MELEMF,MELEMC,MELECB,MELEFC,MELRES,INORM,
  3. & ICHPVO,ICHPSU,LRECP,LRECV,
  4. & IROC,IVITC,IPC,IYC,ICHLIM,ILIINC,ILIINP,IJAC,IJACO)
  5. C************************************************************************
  6. C
  7. C PROJET : CASTEM 2000
  8. C
  9. C NOM : CLI282
  10. C
  11. C DESCRIPTION : Subroutine appellée par CLIM22
  12. C OPTION: 'INOU' -- Jacobian
  13. C
  14. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  15. C
  16. C AUTEUR : S. Kudriakov, DEN/DM2S/SFME/LTMF
  17. C
  18. C************************************************************************
  19. C
  20. C APPELES (Calcul) :
  21. C
  22. C************************************************************************
  23. C
  24. C HISTORIQUE (Anomalies et modifications éventuelles)
  25. C
  26. C HISTORIQUE :
  27. C
  28. C************************************************************************
  29. C
  30. C----------------------------------------------------
  31. C**** Variables de COOPTIO
  32. C----------------------------------------------------
  33. c INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  34. c & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  35. c & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  36. c & ,IECHO, IIMPI, IOSPI
  37. c & ,IDIM, IFICLE, IPREFI
  38. c & ,MCOORD
  39. c & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  40. c & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  41. c & ,NORINC,NORVAL,NORIND,NORVAD
  42. c & ,NUCROU, IPSAUV
  43. C
  44. IMPLICIT INTEGER(I-N)
  45. INTEGER MELEMF,MELEMC,MELECB,INORM,ICHPVO,ICHPSU, IROC,IVITC,IPC
  46. & ,IGAMC,ICHLIM,ICEL,NFAC,IFAC,MELRES,IJACO
  47. & ,NGF,NGC,NLF,NLC,NLCB
  48. & ,ILIINC,ILIINP,IJAC,II,JJ,K
  49. & ,MP, NBEL, NBME, NBSOUS, NKID, NKMT, NMATRI, NP, NRIGE
  50. & ,NSP,I, IYC,J, LRECP,LRECV,KV
  51. REAL*8 VOLU,SURF,RC,PC,UXC,UYC,GAMC,CNX,CNY,CTX,CTY
  52. & ,PSRF,COEF
  53. & ,BR1,BOT,TOP,DUNDUY,GAMF,UNC
  54. & ,COEF5,DPSRUN,DUXDUN,DRDUN,DUYDUN,DUNDUX,DPDUN
  55. & ,GM1,USGM1,RF,UNF,UTF,UXF,UYF,HTF,PF,SF,ECIN
  56. & ,TVECT(2),NVECT(2),WVEC_L(4),WVEC_R(4)
  57. CHARACTER*(8) TYPE
  58. -INC CCOPTIO
  59. -INC SMLMOTS
  60. -INC SMELEME
  61. POINTEUR MELEFC.MELEME
  62. -INC SMLENTI
  63. POINTEUR MLEMC.MLENTI, MLEMCB.MLENTI,MLEMF.MLENTI
  64. -INC SMCHPOI
  65. POINTEUR MPNORM.MPOVAL, MPVOL.MPOVAL, MPSURF.MPOVAL, MPRC.MPOVAL,
  66. & MPVC.MPOVAL, MPPC.MPOVAL, MPLIM.MPOVAL, MPYC.MPOVAL
  67. POINTEUR CELL.IZAFM
  68. C-------------------------------------------------------
  69. -INC SMLREEL
  70. POINTEUR MLRECP.MLREEL, MLRECV.MLREEL
  71. C-------------------------------------------------------
  72. C********* Les Jacobians ******************************
  73. C-------------------------------------------------------
  74. SEGMENT JACEL
  75. REAL*8 JAC(3+NSP,3+NSP)
  76. ENDSEGMENT
  77. POINTEUR JLL.JACEL, WL.JACEL, JRR.JACEL, JTT.JACEL, JT2.JACEL,
  78. & JFL.JACEL, JFR.JACEL, JR.JACEL, JPL.JACEL, JTL.JACEL
  79. C-------------------------------------------------------------
  80. C******* Les fractionines massiques **************************
  81. C-------------------------------------------------------------
  82. SEGMENT FRAMAS
  83. REAL*8 YET(NSP)
  84. ENDSEGMENT
  85. POINTEUR YC.FRAMAS, YF.FRAMAS
  86. C-------------------------------------------------------
  87. C********** Les CP's and CV's ***********************
  88. C-------------------------------------------------------
  89. SEGMENT GCONST
  90. REAL*8 GC(NSP)
  91. ENDSEGMENT
  92. POINTEUR CP.GCONST, CV.GCONST
  93. C-------------------------------------------------------------
  94. C********** Segments for the vectors ***********************
  95. C-------------------------------------------------------------
  96. SEGMENT VECEL
  97. REAL*8 VV(NSP)
  98. ENDSEGMENT
  99. POINTEUR DGDYC.VECEL
  100. C--------------------------------------------------
  101. SEGINI JTT
  102. SEGINI JTL
  103. SEGINI JLL
  104. SEGINI JPL
  105. SEGINI JRR
  106. SEGINI JFL
  107. SEGINI JFR
  108. SEGINI WL
  109. SEGINI JR
  110. SEGINI JT2
  111. C----------------------------------------------------
  112. C**** KRIPAD pour la correspondance global/local
  113. C----------------------------------------------------
  114. CALL KRIPAD(MELEMC,MLEMC)
  115. CALL KRIPAD(MELECB,MLEMCB)
  116. CALL KRIPAD(MELEMF,MLEMF)
  117. C----------------------------------------------------
  118. C**** CHPOINTs de la table DOMAINE
  119. C----------------------------------------------------
  120. CALL LICHT(INORM,MPNORM,TYPE,ICEL)
  121. CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL)
  122. CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL)
  123. C----------------------------------------------------
  124. C**** CHPOINTs des variables
  125. C----------------------------------------------------
  126. CALL LICHT(IROC,MPRC,TYPE,ICEL)
  127. CALL LICHT(IVITC,MPVC,TYPE,ICEL)
  128. CALL LICHT(IPC,MPPC,TYPE,ICEL)
  129. CALL LICHT(IYC,MPYC,TYPE,ICEL)
  130. CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL)
  131. C--------------------------------------------------------
  132. C**** Boucle sur le face pour le calcul des invariants de
  133. C Riemann et du flux
  134. C--------------------------------------------------------
  135. SEGACT MELEFC
  136. NFAC=MELEFC.NUM(/2)
  137. C---------------------------------
  138. C**** Objet MATRIK
  139. C---------------------------------
  140. NRIGE = 7
  141. NMATRI = 1
  142. NKID = 9
  143. NKMT = 7
  144. C---------------------------------
  145. SEGINI MATRIK
  146. IJACO = MATRIK
  147. MATRIK.IRIGEL(1,1) = MELRES
  148. MATRIK.IRIGEL(2,1) = MELRES
  149. C---------------------------------
  150. C**** Matrice non symetrique
  151. C---------------------------------
  152. MATRIK.IRIGEL(7,1) = 2
  153. C---------------------------------
  154. NBME = (3+NSP)*(3+NSP)
  155. NBSOUS = 1
  156. SEGINI IMATRI
  157. IF(IJAC.EQ.1)THEN
  158. MLMOTS=ILIINC
  159. ELSEIF(IJAC.EQ.2)THEN
  160. MLMOTS=ILIINP
  161. ENDIF
  162. SEGACT MLMOTS
  163. MATRIK.IRIGEL(4,1) = IMATRI
  164. C-------------------------------------------
  165. DO 1 J=1,(NSP+3)
  166. KV=(J-1)*(3+NSP)
  167. IMATRI.LISPRI(KV+1) = MLMOTS.MOTS(1)
  168. IMATRI.LISPRI(KV+2) = MLMOTS.MOTS(2)
  169. IMATRI.LISPRI(KV+3) = MLMOTS.MOTS(3)
  170. IMATRI.LISPRI(KV+4) = MLMOTS.MOTS(4)
  171. DO 2 I=1,(NSP-1)
  172. IMATRI.LISPRI(KV+4+I) = MLMOTS.MOTS(4+I)
  173. 2 CONTINUE
  174. 1 CONTINUE
  175. C-----------------------------------------------
  176. SEGDES MLMOTS
  177. MLMOTS=ILIINC
  178. SEGACT MLMOTS
  179. C-----------------------------------------------
  180. DO 3 J=1,(NSP+3)
  181. KV=(J-1)*(3+NSP)
  182. IMATRI.LISDUA(KV+1) = MLMOTS.MOTS(j)
  183. IMATRI.LISDUA(KV+2) = MLMOTS.MOTS(j)
  184. IMATRI.LISDUA(KV+3) = MLMOTS.MOTS(j)
  185. IMATRI.LISDUA(KV+4) = MLMOTS.MOTS(j)
  186. DO 4 I=1,(NSP-1)
  187. IMATRI.LISDUA(KV+4+I) = MLMOTS.MOTS(j)
  188. 4 CONTINUE
  189. 3 CONTINUE
  190. C-----------------------------------------------
  191. C-----------------------------------------------
  192. SEGDES MLMOTS
  193. NBEL = NFAC
  194. NBSOUS = 1
  195. NP = 1
  196. MP = 1
  197. C-----------------------------------------------------------
  198. C-----------------------------------------------------------
  199. DO 5 I=1,NBME
  200. SEGINI CELL
  201. IMATRI.LIZAFM(1,I) = CELL
  202. 5 CONTINUE
  203. C---------------------------------
  204. C---------------------------------
  205. C**** Fin definition MATRIK
  206. C---------------------------------
  207. DO IFAC=1,NFAC,1
  208. NGF=MELEFC.NUM(1,IFAC)
  209. NGC=MELEFC.NUM(2,IFAC)
  210. NLF=MLEMF.LECT(NGF)
  211. NLC=MLEMC.LECT(NGC)
  212. NLCB=MLEMCB.LECT(NGF)
  213. VOLU=MPVOL.VPOCHA(NLC,1)
  214. SURF=MPSURF.VPOCHA(NLF,1)
  215. C In CASTEM les normales sont sortantes
  216. CNX=MPNORM.VPOCHA(NLF,1)
  217. CNY=MPNORM.VPOCHA(NLF,2)
  218. CTX=-1.0D0*CNY
  219. CTY=CNX
  220. C----------------------------------------------
  221. SEGINI CP, CV
  222. MLRECP = LRECP
  223. MLRECV = LRECV
  224. SEGACT MLRECP, MLRECV
  225. DO 10 I=1,(NSP-1)
  226. CP.GC(I)=MLRECP.PROG(I)
  227. CV.GC(I)=MLRECV.PROG(I)
  228. 10 CONTINUE
  229. CP.GC(NSP)=MLRECP.PROG(NSP)
  230. CV.GC(NSP)=MLRECV.PROG(NSP)
  231. C---------------------------------
  232. C Variables au centre
  233. C---------------------------------
  234. RC=MPRC.VPOCHA(NLC,1)
  235. PC=MPPC.VPOCHA(NLC,1)
  236. UXC=MPVC.VPOCHA(NLC,1)
  237. UYC=MPVC.VPOCHA(NLC,2)
  238. SEGINI YC
  239. SEGACT MPYC
  240. DO 100 I=1,(NSP-1)
  241. YC.YET(I)=MPYC.VPOCHA(NLC,I)
  242. 100 CONTINUE
  243. C---------------------------------
  244. C Variables à la face
  245. C---------------------------------
  246. SEGACT MPLIM
  247. HTF=MPLIM.VPOCHA(NLCB,1)
  248. SF=MPLIM.VPOCHA(NLCB,2)
  249. SEGINI YF
  250. DO 101 I=1,(NSP-1)
  251. YF.YET(I)=MPLIM.VPOCHA(NLCB,2+I)
  252. 101 CONTINUE
  253. PF=MPLIM.VPOCHA(NLCB,NSP+2)
  254. c-------------------------------------------------------------
  255. c Computing GAMMA at the face
  256. c-------------------------------------------------------------
  257. top=0.0D0
  258. bot=0.0D0
  259. do 102 i=1,(nsp-1)
  260. top=top+yf.yet(i)*(cp.gc(i)-cp.gc(nsp))
  261. bot=bot+yf.yet(i)*(cv.gc(i)-cv.gc(nsp))
  262. 102 continue
  263. top=cp.gc(nsp)+top
  264. bot=cv.gc(nsp)+bot
  265. GAMF=top/bot
  266. GM1=GAMF-1.0D0
  267. USGM1=1.0D0/GM1
  268. c-------------------------------------------------------------
  269. c Computing GAMMA at the cell-centre
  270. c-------------------------------------------------------------
  271. top=0.0D0
  272. bot=0.0D0
  273. do 103 i=1,(nsp-1)
  274. top=top+yc.yet(i)*(cp.gc(i)-cp.gc(nsp))
  275. bot=bot+yc.yet(i)*(cv.gc(i)-cv.gc(nsp))
  276. 103 continue
  277. top=cp.gc(nsp)+top
  278. bot=cv.gc(nsp)+bot
  279. GAMC=top/bot
  280. C-------------------------------------------------------------
  281. SEGINI DGDYC
  282. do 41 i=1,(nsp-1)
  283. dgdyc.vv(i)=(cp.gc(i)-cp.gc(nsp)-
  284. & GAMC*(cv.gc(i)-cv.gc(nsp)))/bot
  285. 41 continue
  286. c-------------------------------------------------------------
  287. c matrix wl(i,j) represents the derivative of the i-component
  288. c of the vector of primitive variables of the left state with
  289. c respect to the j-component of the vector of the conservative
  290. c variables of the left state.
  291. c
  292. c Here: (rho, ux, uy, p, Y_1,...,Y_(nsp-1)) -
  293. c vector of primitive variables;
  294. c (rho, rho ux, rho uy, rho e, rho Y_1,..., rho Y_(nsp-1)) -
  295. c vector of conservative variables.
  296. c-------------------------------------------------------------
  297. wl.jac(1,1)=1.0d0
  298. wl.jac(1,2)=0.0d0
  299. wl.jac(1,3)=0.0d0
  300. wl.jac(1,4)=0.0d0
  301. do 83 i=1,(nsp-1)
  302. wl.jac(1,4+i)=0.0d0
  303. 83 continue
  304. c------------------------------
  305. wl.jac(2,1)=-UXC/RC
  306. wl.jac(2,2)=1.0d0/RC
  307. wl.jac(2,3)=0.0d0
  308. wl.jac(2,4)=0.0d0
  309. do 84 i=1,(nsp-1)
  310. wl.jac(2,4+i)=0.0d0
  311. 84 continue
  312. c------------------------------
  313. wl.jac(3,1)=-UYC/RC
  314. wl.jac(3,2)=0.0d0
  315. wl.jac(3,3)=1.0d0/RC
  316. wl.jac(3,4)=0.0d0
  317. do 85 i=1,(nsp-1)
  318. wl.jac(3,4+i)=0.0d0
  319. 85 continue
  320. c------------------------------
  321. br1=0.0d0
  322. do 86 i=1,(nsp-1)
  323. br1=br1+dgdyc.vv(i)*yc.yet(i)
  324. 86 continue
  325. br1=br1*PC/(RC*(GAMC-1.0D0))
  326. wl.jac(4,1)=(GAMC-1.0D0)*(UXC*UXC+UYC*UYC)/2.0d0-br1
  327. wl.jac(4,2)=-UXC*(GAMC-1.0D0)
  328. wl.jac(4,3)=-UYC*(GAMC-1.0D0)
  329. wl.jac(4,4)=(GAMC-1.0D0)
  330. do 87 i=1,(nsp-1)
  331. wl.jac(4,4+i)=dgdyc.vv(i)*PC/(RC*(GAMC-1.0D0))
  332. 87 continue
  333. c------------------------------
  334. do 88 i=1,(nsp-1)
  335. do 89 j=1,4
  336. wl.jac(4+i,j)=0.0d0
  337. if(j.eq.1) wl.jac(4+i,j)=-yc.yet(i)/RC
  338. 89 continue
  339. c------------
  340. do 890 j=5,(4+nsp-1)
  341. wl.jac(4+i,j)=0.0d0
  342. if(4+i.eq.j) then
  343. wl.jac(4+i,j)=1.0d0/RC
  344. endif
  345. 890 continue
  346. 88 continue
  347. c------------------------------------------------
  348. UNC=(UXC*CNX)+(UYC*CNY)
  349. IF(UNC .LT. 0.0D0) THEN
  350. c write(*,*) 'impl_insiiiiii', UNC, NLCB
  351. C-----------------------------------------
  352. C Option 'INSU' - inlet BC
  353. C-----------------------------------------
  354. CNX=-1.0D0*CNX
  355. CNY=-1.0D0*CNY
  356. CTX=-1.0D0*CNY
  357. CTY=CNX
  358. UNC=(UXC*CNX)+(UYC*CNY)
  359. UNF=UNC
  360. UTF=0.0D0
  361. UXF=UNF*CNX+UTF*CTX
  362. UYF=UNF*CNY+UTF*CTY
  363. C------------------------------------------------
  364. C******* Densite, vitesse, pression sur le bord
  365. C------------------------------------------------
  366. ECIN=0.5D0*((UXF*UXF)+(UYF*UYF))
  367. PSRF=(GM1/GAMF)*(HTF-ECIN)
  368. RF=PSRF/SF
  369. RF=RF**(1.0D0/GM1)
  370. PF=SF*(RF**GAMF)
  371. C------------------------------
  372. C******* Derivatives w.r.t. PC
  373. C------------------------------
  374. DPSRUN=-1*(GM1/GAMF)*UNC
  375. DRDUN=USGM1*RF/PSRF*DPSRUN
  376. DPDUN=GAMF*PSRF*DRDUN
  377. DUXDUN=CNX
  378. DUYDUN=CNY
  379. C-------------------------------
  380. DUNDUX=CNX
  381. DUNDUY=CNY
  382. C--------------------------------------------------------------
  383. C------------------------------
  384. C******* Derivatives
  385. C------------------------------
  386. wvec_l(1)=RF
  387. wvec_l(2)=UXF
  388. wvec_l(3)=UYF
  389. wvec_l(4)=PF
  390. C--------------------------
  391. wvec_r(1)=RC
  392. wvec_r(2)=UXC
  393. wvec_r(3)=UYC
  394. wvec_r(4)=PC
  395. C--------------------------
  396. nvect(1)=CNX
  397. nvect(2)=CNY
  398. tvect(1)=CTX
  399. tvect(2)=CTY
  400. call copms2(nsp,jfl,jfr,wvec_l,wvec_r,nvect,tvect,
  401. & yf,mpyc,lrecp,lrecv,nlc)
  402. C--------------------------------------------------------------
  403. COEF=SURF/VOLU
  404. JLL=JFL
  405. JRR=JFR
  406. SEGACT JLL
  407. SEGACT JRR
  408. C-------------------------------------------------------------
  409. C Taking in to account the dependance of the variables
  410. C at the face on the variables at the centre (UNC)
  411. C-------------------------------------------------------------
  412. DO 105 I=1,(3+NSP)
  413. DO 1050 J=1,(3+NSP)
  414. IF(J .EQ. 2) THEN
  415. COEF5=(JLL.JAC(I,1)*DRDUN)+(JLL.JAC(I,2)*DUXDUN)+
  416. & (JLL.JAC(I,3)*DUYDUN)+(JLL.JAC(I,4)*DPDUN)
  417. JR.JAC(I,J)=(COEF5*DUNDUX+JRR.JAC(I,J))*COEF
  418. ELSEIF(J .EQ. 3) THEN
  419. COEF5=(JLL.JAC(I,1)*DRDUN)+(JLL.JAC(I,2)*DUXDUN)+
  420. & (JLL.JAC(I,3)*DUYDUN)+(JLL.JAC(I,4)*DPDUN)
  421. JR.JAC(I,J)=(COEF5*DUNDUY+JRR.JAC(I,J))*COEF
  422. ELSE
  423. JR.JAC(I,J)=JRR.JAC(I,J)*COEF
  424. ENDIF
  425. 1050 CONTINUE
  426. 105 CONTINUE
  427. C-------------------------------------------------------------
  428. do 114 i=1,(3+nsp)
  429. do 115 j=1,(3+nsp)
  430. jt2.jac(i,j)=0.0d0
  431. do 116 k=1,(3+nsp)
  432. jt2.jac(i,j)=jt2.jac(i,j)+
  433. & jr.jac(i,k)*wl.jac(k,j)/COEF
  434. 116 continue
  435. 115 continue
  436. 114 continue
  437. C----------------------------------------------------------
  438. C******* Jacobian with respect to conservative variables
  439. C----------------------------------------------------------------
  440. IF(IJAC.EQ.1)THEN
  441. DO 9 II = 1,(3+NSP)
  442. DO 15 JJ = 1,(3+NSP)
  443. KV = (II-1)*(3+NSP)
  444. C----------------------------------
  445. CELL = IMATRI.LIZAFM(1,KV+JJ)
  446. CELL.AM(IFAC,1,1) = JT2.JAC(II,JJ)*COEF
  447. 15 CONTINUE
  448. 9 CONTINUE
  449. ELSEIF(IJAC.EQ.2)THEN
  450. DO 20 II = 1,(3+NSP)
  451. DO 25 JJ = 1,(3+NSP)
  452. KV = (II-1)*(3+NSP)
  453. C----------------------------------
  454. CELL = IMATRI.LIZAFM(1,KV+JJ)
  455. CELL.AM(IFAC,1,1) = JR.JAC(II,JJ)
  456. 25 CONTINUE
  457. 20 CONTINUE
  458. ENDIF
  459. c--------------------------------------------------
  460. ELSE
  461. C---------------------------------------------
  462. C Option --- 'OUTP' - outlet BC
  463. C---------------------------------------------
  464. C------------------------------
  465. C******* Derivatives
  466. C------------------------------
  467. wvec_l(1)=RC
  468. wvec_l(2)=UXC
  469. wvec_l(3)=UYC
  470. wvec_l(4)=PC
  471. C--------------------------
  472. wvec_r(1)=RC
  473. wvec_r(2)=UXC
  474. wvec_r(3)=UYC
  475. wvec_r(4)=PF
  476. C--------------------------
  477. nvect(1)=CNX
  478. nvect(2)=CNY
  479. tvect(1)=CTX
  480. tvect(2)=CTY
  481. call copmsp(nsp,jpl,jll,wvec_l,wvec_r,nvect,tvect,
  482. & mpyc,lrecp,lrecv,nlc,nlc)
  483. C----------------------------------------
  484. COEF=-SURF/VOLU
  485. C----------------------------------------
  486. JTT=JLL
  487. JTL=JPL
  488. SEGACT JTT
  489. SEGACT JTL
  490. C----------------------------------------------------------------
  491. C******* Jacobian with respect to conservative variables
  492. C----------------------------------------------------------------
  493. IF(IJAC.EQ.1)THEN
  494. DO 99 II = 1,(3+NSP)
  495. DO 155 JJ = 1,(3+NSP)
  496. KV = (II-1)*(3+NSP)
  497. C----------------------------------
  498. CELL = IMATRI.LIZAFM(1,KV+JJ)
  499. CELL.AM(IFAC,1,1) = JTT.JAC(II,JJ)*COEF
  500. 155 CONTINUE
  501. 99 CONTINUE
  502. ELSEIF(IJAC.EQ.2)THEN
  503. DO 200 II = 1,(3+NSP)
  504. DO 255 JJ = 1,(3+NSP)
  505. KV = (II-1)*(3+NSP)
  506. C----------------------------------
  507. CELL = IMATRI.LIZAFM(1,KV+JJ)
  508. CELL.AM(IFAC,1,1) = JTL.JAC(II,JJ)*COEF
  509. 255 CONTINUE
  510. 200 CONTINUE
  511. ENDIF
  512. C---------------------------------------------------------------
  513. ENDIF
  514. C----------------------------------------------------------------
  515. ENDDO
  516. C
  517. SEGDES MELEFC
  518. C
  519. SEGDES MLEMC
  520. SEGDES MLEMCB
  521. SEGDES MLEMF
  522. C
  523. SEGDES MPNORM
  524. SEGDES MPVOL
  525. SEGDES MPSURF
  526. SEGDES MPRC
  527. SEGDES MPPC
  528. SEGDES MPVC
  529. SEGDES MPYC
  530. SEGDES MPLIM
  531. SEGDES YC
  532. SEGDES YF
  533. SEGDES CP
  534. SEGDES CV
  535. SEGDES JLL
  536. SEGDES JRR
  537. SEGDES JTT
  538. SEGDES JTL
  539. SEGDES JFL
  540. SEGDES JFR
  541. SEGDES JR
  542. SEGDES WL
  543. SEGDES DGDYC
  544. SEGDES MATRIK
  545. DO 80 II=1,NBME
  546. CELL = IMATRI.LIZAFM(1,II)
  547. SEGDES CELL
  548. 80 CONTINUE
  549. SEGDES IMATRI
  550. C---------------------------------------------
  551. 9999 CONTINUE
  552. RETURN
  553. END
  554.  
  555.  
  556.  
  557.  
  558.  
  559.  
  560.  
  561.  
  562.  
  563.  
  564.  
  565.  
  566.  
  567.  

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