Télécharger cli163.eso

Retour à la liste

Numérotation des lignes :

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

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