Télécharger cli113.eso

Retour à la liste

Numérotation des lignes :

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

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