Télécharger cli123.eso

Retour à la liste

Numérotation des lignes :

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

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