Télécharger cli182.eso

Retour à la liste

Numérotation des lignes :

  1. C CLI182 SOURCE PV 16/11/17 21:58:41 9180
  2. SUBROUTINE CLI182(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 : CLI182
  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,GAMC,CNX,CNY,CTX,CTY,GM1
  62. & ,RC,UXC,UYC,UNC,CC,ACEL,BCEL,CCEL
  63. & ,PC,PSRF,RHOUF,P,PSTAR,RHO,UN,UT,UX,UY,ECIN
  64. & ,DADR,DADP,DBDR,DBDP,DBDUX,DBDUY,DCDP
  65. & ,DPSDA,DPSDB,DPSDC,DPDR,DPDP,DPDUX,DPDUY
  66. & ,DUXDP,DUYDP,DECINP,DF2DP,DF3DP,DF4DP
  67. & ,DRORO,DROUX,DROUY,DROP
  68. & ,DUXRO,DUXUX,DUXUY,DUXP
  69. & ,DUYRO,DUYUX,DUYUY,DUYP
  70. & ,DPRO,DPUX,DPUY,DPP
  71. & ,COEF1,COEF2,COEF3,COEF
  72. C & ,YCAC,YCAC2,YCAC3,XCAC,F1,F2,F3,F4,USGM1
  73. CHARACTER*(8) TYPE
  74. C
  75. C
  76. C**** KRIPAD pour la correspondance global/local
  77. C
  78. CALL KRIPAD(MELEMC,MLEMC)
  79. C SEGINI MLEMC
  80. CALL KRIPAD(MELECB,MLEMCB)
  81. C SEGINI MLEMCB
  82. CALL KRIPAD(MELEMF,MLEMF)
  83. C SEGINI MLEMF
  84. C
  85. C**** CHPOINTs de la table DOMAINE
  86. C
  87. CALL LICHT(INORM,MPNORM,TYPE,ICEL)
  88. 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
  100. CALL LICHT(IROC,MPRC,TYPE,ICEL)
  101. CALL LICHT(IVITC,MPVC,TYPE,ICEL)
  102. CALL LICHT(IPC,MPPC,TYPE,ICEL)
  103. CALL LICHT(IGAMC,MPGAMC,TYPE,ICEL)
  104. CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL)
  105. C
  106. C SEGACT *MOD
  107. C SEGACT *MOD
  108. C SEGACT *MOD
  109. C SEGACT *MOD
  110. C SEGACT *MOD
  111. C
  112. C
  113. C**** Boucle sur le face pour le calcul des invariants de
  114. C Riemann et du flux
  115. C
  116. SEGACT MELEFC
  117. NFAC=MELEFC.NUM(/2)
  118. C
  119. C**** Objet MATRIK
  120. C
  121. NRIGE = 7
  122. NMATRI = 1
  123. NKID = 9
  124. NKMT = 7
  125. C
  126. SEGINI MATRIK
  127. IJACO = MATRIK
  128. MATRIK.IRIGEL(1,1) = MELRES
  129. MATRIK.IRIGEL(2,1) = MELRES
  130. C
  131. C**** Matrice non symetrique
  132. C
  133. MATRIK.IRIGEL(7,1) = 2
  134. C
  135. NBME = 16
  136. NBSOUS = 1
  137. SEGINI IMATRI
  138. IF(IJAC.EQ.1)THEN
  139. MLMOTS=ILIINC
  140. ELSEIF(IJAC.EQ.2)THEN
  141. MLMOTS=ILIINP
  142. ENDIF
  143. SEGACT MLMOTS
  144. MATRIK.IRIGEL(4,1) = IMATRI
  145. C
  146. IMATRI.LISPRI(1) = MLMOTS.MOTS(1)
  147. IMATRI.LISPRI(2) = MLMOTS.MOTS(2)
  148. IMATRI.LISPRI(3) = MLMOTS.MOTS(3)
  149. IMATRI.LISPRI(4) = MLMOTS.MOTS(4)
  150. IMATRI.LISPRI(5) = MLMOTS.MOTS(1)
  151. IMATRI.LISPRI(6) = MLMOTS.MOTS(2)
  152. IMATRI.LISPRI(7) = MLMOTS.MOTS(3)
  153. IMATRI.LISPRI(8) = MLMOTS.MOTS(4)
  154. IMATRI.LISPRI(9) = MLMOTS.MOTS(1)
  155. IMATRI.LISPRI(10) = MLMOTS.MOTS(2)
  156. IMATRI.LISPRI(11) = MLMOTS.MOTS(3)
  157. IMATRI.LISPRI(12) = MLMOTS.MOTS(4)
  158. IMATRI.LISPRI(13) = MLMOTS.MOTS(1)
  159. IMATRI.LISPRI(14) = MLMOTS.MOTS(2)
  160. IMATRI.LISPRI(15) = MLMOTS.MOTS(3)
  161. IMATRI.LISPRI(16) = MLMOTS.MOTS(4)
  162. C
  163. SEGDES MLMOTS
  164. MLMOTS=ILIINC
  165. SEGACT MLMOTS
  166. C
  167. IMATRI.LISDUA(1) = MLMOTS.MOTS(1)
  168. IMATRI.LISDUA(2) = MLMOTS.MOTS(1)
  169. IMATRI.LISDUA(3) = MLMOTS.MOTS(1)
  170. IMATRI.LISDUA(4) = MLMOTS.MOTS(1)
  171. IMATRI.LISDUA(5) = MLMOTS.MOTS(2)
  172. IMATRI.LISDUA(6) = MLMOTS.MOTS(2)
  173. IMATRI.LISDUA(7) = MLMOTS.MOTS(2)
  174. IMATRI.LISDUA(8) = MLMOTS.MOTS(2)
  175. IMATRI.LISDUA(9) = MLMOTS.MOTS(3)
  176. IMATRI.LISDUA(10) = MLMOTS.MOTS(3)
  177. IMATRI.LISDUA(11) = MLMOTS.MOTS(3)
  178. IMATRI.LISDUA(12) = MLMOTS.MOTS(3)
  179. IMATRI.LISDUA(13) = MLMOTS.MOTS(4)
  180. IMATRI.LISDUA(14) = MLMOTS.MOTS(4)
  181. IMATRI.LISDUA(15) = MLMOTS.MOTS(4)
  182. IMATRI.LISDUA(16) = MLMOTS.MOTS(4)
  183. C
  184. SEGDES MLMOTS
  185. NBEL = NFAC
  186. NBSOUS = 1
  187. NP = 1
  188. MP = 1
  189. SEGINI RR , RUX , RUY , RRET ,
  190. & UXR , UXUX , UXUY , UXRET ,
  191. & UYR , UYUX , UYUY , UYRET ,
  192. & RETR , RETUX , RETUY , RETRET
  193. C
  194. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  195. C Primale = IMATRI.LISPRI(1) = 'RN'
  196. C -> IMATRI.LIZAFM(1,1) = RR
  197. C
  198. IMATRI.LIZAFM(1,1) = RR
  199. IMATRI.LIZAFM(1,2) = RUX
  200. IMATRI.LIZAFM(1,3) = RUY
  201. IMATRI.LIZAFM(1,4) = RRET
  202. IMATRI.LIZAFM(1,5) = UXR
  203. IMATRI.LIZAFM(1,6) = UXUX
  204. IMATRI.LIZAFM(1,7) = UXUY
  205. IMATRI.LIZAFM(1,8) = UXRET
  206. IMATRI.LIZAFM(1,9) = UYR
  207. IMATRI.LIZAFM(1,10) = UYUX
  208. IMATRI.LIZAFM(1,11) = UYUY
  209. IMATRI.LIZAFM(1,12) = UYRET
  210. IMATRI.LIZAFM(1,13) = RETR
  211. IMATRI.LIZAFM(1,14) = RETUX
  212. IMATRI.LIZAFM(1,15) = RETUY
  213. IMATRI.LIZAFM(1,16) = RETRET
  214. C
  215. SEGDES MATRIK
  216. SEGDES IMATRI
  217. C
  218. C**** Fin definition MATRIK
  219. C
  220. DO IFAC=1,NFAC,1
  221. NGF=MELEFC.NUM(1,IFAC)
  222. NGC=MELEFC.NUM(2,IFAC)
  223. NLF=MLEMF.LECT(NGF)
  224. NLC=MLEMC.LECT(NGC)
  225. NLCB=MLEMCB.LECT(NGF)
  226. VOLU=MPVOL.VPOCHA(NLC,1)
  227. SURF=MPSURF.VPOCHA(NLF,1)
  228. C In CASTEM les normales sont sortantes
  229. CNX=-1*MPNORM.VPOCHA(NLF,1)
  230. CNY=-1*MPNORM.VPOCHA(NLF,2)
  231. CTX=-1*CNY
  232. CTY=CNX
  233. C Variables au centre
  234. GAMC=MPGAMC.VPOCHA(NLC,1)
  235. GM1=GAMC-1.0D0
  236. PC=MPPC.VPOCHA(NLC,1)
  237. RC=MPRC.VPOCHA(NLC,1)
  238. UXC=MPVC.VPOCHA(NLC,1)
  239. UYC=MPVC.VPOCHA(NLC,2)
  240. UNC=(UXC*CNX)+(UYC*CNY)
  241. CC=GAMC*PC/RC
  242. CC=CC**0.5D0
  243. C Variables à la face
  244. RHOUF=MPLIM.VPOCHA(NLCB,1)
  245. PSRF=MPLIM.VPOCHA(NLCB,2)
  246. C
  247. C******* Variables à l'interface
  248. C
  249. ACEL=CC/GAMC
  250. BCEL=ACEL - UNC
  251. CCEL=(RHOUF / PC) * PSRF
  252. PSTAR=BCEL+(((BCEL*BCEL) + (4*ACEL*CCEL))**0.5D0)
  253. PSTAR=PSTAR/(2*ACEL)
  254. P=PSTAR*PC
  255. RHO=P/PSRF
  256. UN=RHOUF/RHO
  257. UT=0.0D0
  258. UX=UN*CNX+UT*CTX
  259. UY=UN*CNY+UT*CTY
  260. ECIN=0.5D0*((UX*UX)+(UY*UY))
  261. C
  262. C******* Derivatives of ACEL, BCEL, CCEL with respect to the
  263. C variables at centers
  264. C
  265. DADR = -0.5D0*CC/(RC*GAMC)
  266. DADP = 0.5D0*CC/(PC*GAMC)
  267. DBDR = DADR
  268. DBDP = DADP
  269. DBDUX = -1*CNX
  270. DBDUY = -1*CNY
  271. DCDP = -1*CCEL/PC
  272. C
  273. C******* Derivatives of PSTAR with respect to ACEL,BCEL,CCEL
  274. C
  275. DPSDA = -1*(PSTAR*PSTAR)/((2*ACEL*PSTAR)-BCEL)
  276. DPSDB = PSTAR/((2*ACEL*PSTAR)-BCEL)
  277. DPSDC = 1/((2*ACEL*PSTAR)-BCEL)
  278. C
  279. C******* Derivatives of PC*PSTAR with respect to RC,PC,UXC,UYC
  280. C
  281. DPDR=(DPSDA*DADR)+(DPSDB*DBDR)
  282. DPDR=DPDR*PC
  283. DPDP=(DPSDA*DADP)+(DPSDB*DBDP)+(DPSDC*DCDP)
  284. DPDP=(DPDP*PC)+PSTAR
  285. DPDUX=(DPSDB*DBDUX)*PC
  286. DPDUY=(DPSDB*DBDUY)*PC
  287. CC
  288. CC******* Test 1
  289. CC
  290. CC We check dpdp
  291. CC
  292. C YCAC=P
  293. C XCAC=PC
  294. C PC=PC*(1+1.0D-4)
  295. C CC=GAMC*PC/RC
  296. C CC=CC**0.5D0
  297. CC
  298. C ACEL=CC/GAMC
  299. C BCEL=ACEL - UNC
  300. C CCEL=(RHOUF / PC) * PSRF
  301. C PSTAR=BCEL+(((BCEL*BCEL) + (4*ACEL*CCEL))**0.5D0)
  302. C PSTAR=PSTAR/(2*ACEL)
  303. C P=PSTAR*PC
  304. C write(*,*) ((P - YCAC)/(PC - XCAC)), DPDP
  305. CC
  306. CC We check dpdr
  307. CC
  308. C YCAC=P
  309. C XCAC=RC
  310. C RC=RC*(1+1.0D-4)
  311. C CC=GAMC*PC/RC
  312. C CC=CC**0.5D0
  313. CC
  314. C ACEL=CC/GAMC
  315. C BCEL=ACEL - UNC
  316. C CCEL=(RHOUF / PC) * PSRF
  317. C PSTAR=BCEL+(((BCEL*BCEL) + (4*ACEL*CCEL))**0.5D0)
  318. C PSTAR=PSTAR/(2*ACEL)
  319. C P=PSTAR*PC
  320. C write(*,*) ((P - YCAC)/(RC - XCAC)), DPDR
  321. CC
  322. CC We check dpdux
  323. CC
  324. C YCAC=P
  325. C XCAC=UXC
  326. C UXC=UXC*(1+1.0D-4)
  327. C UNC=(UXC*CNX)+(UYC*CNY)
  328. C CC=GAMC*PC/RC
  329. C CC=CC**0.5D0
  330. CC
  331. C ACEL=CC/GAMC
  332. C BCEL=ACEL - UNC
  333. C CCEL=(RHOUF / PC) * PSRF
  334. C PSTAR=BCEL+(((BCEL*BCEL) + (4*ACEL*CCEL))**0.5D0)
  335. C PSTAR=PSTAR/(2*ACEL)
  336. C P=PSTAR*PC
  337. C write(*,*) ((P - YCAC)/(UXC - XCAC)), DPDUX
  338. CC
  339. CC We check dpduy
  340. CC
  341. C YCAC=P
  342. C XCAC=UYC
  343. C UYC=UYC*(1+1.0D-4)
  344. C UNC=(UXC*CNX)+(UYC*CNY)
  345. C CC=GAMC*PC/RC
  346. C CC=CC**0.5D0
  347. CC
  348. C ACEL=CC/GAMC
  349. C BCEL=ACEL - UNC
  350. C CCEL=(RHOUF / PC) * PSRF
  351. C PSTAR=BCEL+(((BCEL*BCEL) + (4*ACEL*CCEL))**0.5D0)
  352. C PSTAR=PSTAR/(2*ACEL)
  353. C P=PSTAR*PC
  354. C write(*,*) ((P - YCAC)/(UYC - XCAC)), DPDUY
  355. C
  356. CC
  357. CC******* Fin test 1
  358. C
  359. C
  360. C******* The interfacial state is thus given by
  361. C RHOUF,PSRF,P and UT (UT=0)
  362. C Let us now compute the derivative of the residuum with
  363. C respect to P.
  364. C The residuum is given by
  365. C
  366. C F1=RHOUF*SURF/VOLU
  367. C F2=(RHOUF*UX+P*CNX)*SURF/VOLU
  368. C F3=(RHOUF*UY+P*CNY)*SURF/VOLU
  369. C F4=(RHOUF*((GAMC*USGM1*PSRF)+ECIN))*SURF/VOLU
  370. C
  371. C In these expression: UX, UY, P and ECIN depends on P
  372. C
  373. COEF=SURF/VOLU
  374. DUXDP=-1*UX/P
  375. DUYDP=-1*UY/P
  376. DECINP=-2*ECIN/P
  377. DF2DP=((RHOUF*DUXDP)+CNX)*COEF
  378. DF3DP=((RHOUF*DUYDP)+CNY)*COEF
  379. DF4DP=(RHOUF*DECINP)*COEF
  380.  
  381. CC
  382. CC******* Test 2
  383. CC
  384. CC We check DF2P
  385. CC
  386. C USGM1=1/GM1
  387. C F2=(RHOUF*UX+P*CNX)*COEF
  388. C F3=(RHOUF*UY+P*CNY)*COEF
  389. C F4=(RHOUF*((GAMC*USGM1*PSRF)+ECIN))*COEF
  390. CC
  391. C YCAC=F2
  392. C YCAC2=F3
  393. C YCAC3=F4
  394. C XCAC=P
  395. C P=P*(1+1.0D-4)
  396. C RHO=P/PSRF
  397. C UN=RHOUF/RHO
  398. C UT=0.0D0
  399. C UX=UN*CNX+UT*CTX
  400. C UY=UN*CNY+UT*CTY
  401. C ECIN=0.5D0*((UX*UX)+(UY*UY))
  402. C
  403. C F2=(RHOUF*UX+P*CNX)*COEF
  404. C F3=(RHOUF*UY+P*CNY)*COEF
  405. C F4=(RHOUF*((GAMC*USGM1*PSRF)+ECIN))*COEF
  406. C write(*,*) ((F2 -YCAC)/(P-XCAC)),DF2DP
  407. C write(*,*) ((F3 -YCAC2)/(P-XCAC)),DF3DP
  408. C write(*,*) ((F4 -YCAC3)/(P-XCAC)),DF4DP
  409. CC
  410. CC******* Fin test 2
  411. CC
  412. C
  413. C******* Jacobian with respect to primitive variables
  414. C
  415. C DF1D...
  416. C
  417. DRORO=0.0D0
  418. DROUX=0.0D0
  419. DROUY=0.0D0
  420. DROP=0.0D0
  421. C
  422. C DF4D...
  423. C
  424. DPRO=DF4DP*DPDR
  425. DPUX=DF4DP*DPDUX
  426. DPUY=DF4DP*DPDUY
  427. DPP=DF4DP*DPDP
  428. C
  429. C DF2D...
  430. C
  431. DUXRO=DF2DP*DPDR
  432. DUXUX=DF2DP*DPDUX
  433. DUXUY=DF2DP*DPDUY
  434. DUXP =DF2DP*DPDP
  435. C
  436. C DF3D
  437. C
  438. DUYRO=DF3DP*DPDR
  439. DUYUX=DF3DP*DPDUX
  440. DUYUY=DF3DP*DPDUY
  441. DUYP =DF3DP*DPDP
  442. C
  443. C******* Jacobian with respect to conservative variables
  444. C
  445. IF(IJAC.EQ.1)THEN
  446.  
  447. COEF1=-1.0D0*UXC/RC
  448. COEF2=-1.0D0*UYC/RC
  449. COEF3=0.5D0*GM1*((UXC*UXC)+(UYC*UYC))
  450. C
  451. RR.AM(IFAC,1,1)=DRORO+(DROUX*COEF1)+(DROUY*COEF2)+(DROP
  452. $ *COEF3)
  453. RUX.AM(IFAC,1,1)=(DROUX/RC)-((UXC*GM1)*DROP)
  454. RUY.AM(IFAC,1,1)=(DROUY/RC)-((UYC*GM1)*DROP)
  455. RRET.AM(IFAC,1,1)=GM1*DROP
  456. C
  457. UXR.AM(IFAC,1,1)=DUXRO+(DUXUX*COEF1)+(DUXUY*COEF2)+(DUXP
  458. $ *COEF3)
  459. UXUX.AM(IFAC,1,1)=(DUXUX/RC)-((UXC*GM1)*DUXP)
  460. UXUY.AM(IFAC,1,1)=(DUXUY/RC)-((UYC*GM1)*DUXP)
  461. UXRET.AM(IFAC,1,1)=GM1*DUXP
  462. C
  463. UYR.AM(IFAC,1,1)=DUYRO+(DUYUX*COEF1)+(DUYUY*COEF2)+(DUYP
  464. $ *COEF3)
  465. UYUX.AM(IFAC,1,1)=(DUYUX/RC)-((UXC*GM1)*DUYP)
  466. UYUY.AM(IFAC,1,1)=(DUYUY/RC)-((UYC*GM1)*DUYP)
  467. UYRET.AM(IFAC,1,1)=GM1*DUYP
  468. C
  469. RETR.AM(IFAC,1,1)=DPRO+(DPUX*COEF1)+(DPUY*COEF2)+(DPP
  470. $ *COEF3)
  471. RETUX.AM(IFAC,1,1)=(DPUX/RC)-((UXC*GM1)*DPP)
  472. RETUY.AM(IFAC,1,1)=(DPUY/RC)-((UYC*GM1)*DPP)
  473. RETRET.AM(IFAC,1,1)=GM1*DPP
  474. C
  475. ELSEIF(IJAC.EQ.2)THEN
  476. RR.AM(IFAC,1,1)=DRORO
  477. RUX.AM(IFAC,1,1)=DROUX
  478. RUY.AM(IFAC,1,1)=DROUY
  479. RRET.AM(IFAC,1,1)=DROP
  480. C
  481. UXR.AM(IFAC,1,1)=DUXRO
  482. UXUX.AM(IFAC,1,1)=DUXUX
  483. UXUY.AM(IFAC,1,1)=DUXUY
  484. UXRET.AM(IFAC,1,1)=DUXP
  485. C
  486. UYR.AM(IFAC,1,1)=DUYRO
  487. UYUX.AM(IFAC,1,1)=DUYUX
  488. UYUY.AM(IFAC,1,1)=DUYUY
  489. UYRET.AM(IFAC,1,1)=DUYP
  490. C
  491. RETR.AM(IFAC,1,1)=DPRO
  492. RETUX.AM(IFAC,1,1)=DPUX
  493. RETUY.AM(IFAC,1,1)=DPUY
  494. RETRET.AM(IFAC,1,1)=DPP
  495. ENDIF
  496. ENDDO
  497. C
  498. SEGDES MELEFC
  499. C
  500. SEGSUP MLEMC
  501. SEGSUP MLEMCB
  502. SEGSUP MLEMF
  503. C
  504. SEGDES MPNORM
  505. SEGDES MPVOL
  506. SEGDES MPSURF
  507. SEGDES MPRC
  508. SEGDES MPPC
  509. SEGDES MPVC
  510. SEGDES MPGAMC
  511. SEGDES MPLIM
  512. C
  513. SEGDES RR , RUX , RUY , RRET ,
  514. & UXR , UXUX , UXUY , UXRET ,
  515. & UYR , UYUX , UYUY , UYRET ,
  516. & RETR , RETUX , RETUY , RETRET
  517. C
  518. 9999 CONTINUE
  519. RETURN
  520. END
  521.  
  522.  
  523.  
  524.  
  525.  
  526.  
  527.  
  528.  
  529.  
  530.  
  531.  
  532.  
  533.  

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