Télécharger cli123.eso

Retour à la liste

Numérotation des lignes :

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

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