Télécharger prija2.eso

Retour à la liste

Numérotation des lignes :

  1. C PRIJA2 SOURCE PV 16/11/17 22:01:09 9180
  2. SUBROUTINE PRIJA2(MELEMC,MLMCON,MLMPRI,IRN,IVN,IPN,IGAMN,IJACO)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : PRIJA2
  9. C
  10. C DESCRIPTION : Voir PRIJA1
  11. C Cas deux dimensions, gaz "calorically perfect"
  12. C
  13. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  14. C
  15. C AUTEUR : A. BECCANTINI, DM2S/SFME/LTMF
  16. C
  17. C************************************************************************
  18. C
  19. C ENTREES
  20. C
  21. C MELEMC : SPG des CHPOINTs
  22. C
  23. C MLMPRI : liste des inconnues primitives (P,v,T)
  24. C
  25. C MLMCON : liste des inconnues conservatives (rho, rhov, rhoet)
  26. C
  27. C IRN : CHPOINT CENTRE contenant la densité
  28. C
  29. C IVN : CHPOINT CENTRE contenant la vitesse
  30. C
  31. C IPN : CHPOINT CENTRE contenant la pression
  32. C
  33. C IGAM : CHPOINT CENTRE contenant le gamma
  34. C
  35. C SORTIES
  36. C
  37. C IJACO : pointeur de la MATRIK jacobienne
  38. C
  39. C************************************************************************
  40. C
  41. C HISTORIQUE (Anomalies et modifications éventuelles)
  42. C
  43. C HISTORIQUE : crée le 29.05.02
  44. C
  45. C************************************************************************
  46. C
  47. C
  48. C N.B.: On suppose qu'on a déjà controllé RO, P > 0
  49. C GAMMA \in (1,3)
  50. C Si non il faut le faire!!!
  51. C
  52. C************************************************************************
  53. C
  54. C
  55. C**** Variables de COOPTIO
  56. C
  57. C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  58. C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  59. C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  60. C & ,IECHO, IIMPI, IOSPI
  61. C & ,IDIM
  62. C & ,MCOORD
  63. C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  64. C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  65. C & ,NORINC,NORVAL,NORIND,NORVAD
  66. C & ,NUCROU, IPSAUV, IFICLE, IPREFI
  67. C
  68. IMPLICIT INTEGER(I-N)
  69. INTEGER IRN, IVN, IPN, IGAMN, IJACO, IGEOM, ICEN
  70. & , MP, NP, NBEL, NBME, NBSOUS
  71. & , NRIGE, NMATRI, NKID, NKMT
  72. C
  73. REAL*8 RHO,UX, UY, UZ, GAMMA
  74. CHARACTER*8 TYPE
  75. C
  76. C**** LES INCLUDES
  77. C
  78.  
  79. -INC PPARAM
  80. -INC CCOPTIO
  81. -INC SMCHPOI
  82. -INC SMELEME
  83. -INC SMLMOTS
  84. -INC SMLENTI
  85. POINTEUR MPVN.MPOVAL, MPRN.MPOVAL, MPGAM.MPOVAL
  86. POINTEUR MELEMC.MELEME
  87. POINTEUR RR.IZAFM, RUX.IZAFM, RUY.IZAFM, RUZ.IZAFM, RP.IZAFM,
  88. & UXR.IZAFM, UXUX.IZAFM, UXUY.IZAFM, UXUZ.IZAFM,
  89. & UXP.IZAFM,
  90. & UYR.IZAFM, UYUX.IZAFM, UYUY.IZAFM, UYUZ.IZAFM,
  91. & UYP.IZAFM,
  92. & UZR.IZAFM, UZUX.IZAFM, UZUY.IZAFM, UZUZ.IZAFM,
  93. & UZP.IZAFM,
  94. & RETR.IZAFM, RETUX.IZAFM, RETUY.IZAFM, RETUZ.IZAFM,
  95. & RETP.IZAFM
  96. POINTEUR MLMPRI.MLMOTS, MLMCON.MLMOTS
  97. C
  98. CALL LICHT(IVN,MPVN,TYPE,IGEOM)
  99. CALL LICHT(IRN,MPRN,TYPE,IGEOM)
  100. CALL LICHT(IGAMN,MPGAM,TYPE,IGEOM)
  101. C
  102. C SEGACT MPVN*MOD
  103. C SEGACT MPPN*MOD
  104. C SEGACT MPRN*MOD
  105. C
  106. C**** Maillage des inconnues primales = Maillage des inconnues primales
  107. C = MELEMC
  108. C
  109. C
  110. NRIGE = 7
  111. NMATRI = 1
  112. NKID = 9
  113. NKMT = 7
  114. C
  115. SEGINI MATRIK
  116. IJACO = MATRIK
  117. MATRIK.IRIGEL(1,1) = MELEMC
  118. MATRIK.IRIGEL(2,1) = MELEMC
  119. C
  120. C**** Matrice non symetrique
  121. C
  122. MATRIK.IRIGEL(7,1) = 2
  123. C
  124. NBME = 25
  125. NBSOUS = 1
  126. SEGINI IMATRI
  127. MATRIK.IRIGEL(4,1) = IMATRI
  128. C
  129. SEGACT MLMPRI
  130. C-----------------------------------------------
  131. IMATRI.LISPRI(1) = MLMPRI.MOTS(1)
  132. IMATRI.LISPRI(2) = MLMPRI.MOTS(2)
  133. IMATRI.LISPRI(3) = MLMPRI.MOTS(3)
  134. IMATRI.LISPRI(4) = MLMPRI.MOTS(4)
  135. IMATRI.LISPRI(5) = MLMPRI.MOTS(5)
  136. C-----------------------------------------------
  137. IMATRI.LISPRI(6) = MLMPRI.MOTS(1)
  138. IMATRI.LISPRI(7) = MLMPRI.MOTS(2)
  139. IMATRI.LISPRI(8) = MLMPRI.MOTS(3)
  140. IMATRI.LISPRI(9) = MLMPRI.MOTS(4)
  141. IMATRI.LISPRI(10) = MLMPRI.MOTS(5)
  142. C-----------------------------------------------
  143. IMATRI.LISPRI(11) = MLMPRI.MOTS(1)
  144. IMATRI.LISPRI(12) = MLMPRI.MOTS(2)
  145. IMATRI.LISPRI(13) = MLMPRI.MOTS(3)
  146. IMATRI.LISPRI(14) = MLMPRI.MOTS(4)
  147. IMATRI.LISPRI(15) = MLMPRI.MOTS(5)
  148. C-----------------------------------------------
  149. IMATRI.LISPRI(16) = MLMPRI.MOTS(1)
  150. IMATRI.LISPRI(17) = MLMPRI.MOTS(2)
  151. IMATRI.LISPRI(18) = MLMPRI.MOTS(3)
  152. IMATRI.LISPRI(19) = MLMPRI.MOTS(4)
  153. IMATRI.LISPRI(20) = MLMPRI.MOTS(5)
  154. C-----------------------------------------------
  155. IMATRI.LISPRI(21) = MLMPRI.MOTS(1)
  156. IMATRI.LISPRI(22) = MLMPRI.MOTS(2)
  157. IMATRI.LISPRI(23) = MLMPRI.MOTS(3)
  158. IMATRI.LISPRI(24) = MLMPRI.MOTS(4)
  159. IMATRI.LISPRI(25) = MLMPRI.MOTS(5)
  160. C-----------------------------------------------
  161. SEGDES MLMPRI
  162. C
  163. SEGACT MLMCON
  164. C-----------------------------------------------
  165. IMATRI.LISDUA(1) = MLMCON.MOTS(1)
  166. IMATRI.LISDUA(2) = MLMCON.MOTS(1)
  167. IMATRI.LISDUA(3) = MLMCON.MOTS(1)
  168. IMATRI.LISDUA(4) = MLMCON.MOTS(1)
  169. IMATRI.LISDUA(5) = MLMCON.MOTS(1)
  170. C-----------------------------------------------
  171. IMATRI.LISDUA(6) = MLMCON.MOTS(2)
  172. IMATRI.LISDUA(7) = MLMCON.MOTS(2)
  173. IMATRI.LISDUA(8) = MLMCON.MOTS(2)
  174. IMATRI.LISDUA(9) = MLMCON.MOTS(2)
  175. IMATRI.LISDUA(10) = MLMCON.MOTS(2)
  176. C-----------------------------------------------
  177. IMATRI.LISDUA(11) = MLMCON.MOTS(3)
  178. IMATRI.LISDUA(12) = MLMCON.MOTS(3)
  179. IMATRI.LISDUA(13) = MLMCON.MOTS(3)
  180. IMATRI.LISDUA(14) = MLMCON.MOTS(3)
  181. IMATRI.LISDUA(15) = MLMCON.MOTS(3)
  182. C-----------------------------------------------
  183. IMATRI.LISDUA(16) = MLMCON.MOTS(4)
  184. IMATRI.LISDUA(17) = MLMCON.MOTS(4)
  185. IMATRI.LISDUA(18) = MLMCON.MOTS(4)
  186. IMATRI.LISDUA(19) = MLMCON.MOTS(4)
  187. IMATRI.LISDUA(20) = MLMCON.MOTS(4)
  188. C-----------------------------------------------
  189. IMATRI.LISDUA(21) = MLMCON.MOTS(5)
  190. IMATRI.LISDUA(22) = MLMCON.MOTS(5)
  191. IMATRI.LISDUA(23) = MLMCON.MOTS(5)
  192. IMATRI.LISDUA(24) = MLMCON.MOTS(5)
  193. IMATRI.LISDUA(25) = MLMCON.MOTS(5)
  194. C-----------------------------------------------
  195. SEGDES MLMCON
  196. C
  197. SEGACT MELEMC
  198. NBEL = MELEMC.NUM(/2)
  199. SEGDES MELEMC
  200. NBSOUS = 1
  201. NP = 1
  202. MP = 1
  203. SEGINI RR , RUX , RUY , RUZ, RP ,
  204. & UXR , UXUX , UXUY , UXUZ, UXP ,
  205. & UYR , UYUX , UYUY , UYUZ, UYP ,
  206. & UZR , UZUX , UZUY , UZUZ, UZP ,
  207. & RETR , RETUX , RETUY , RETUZ, RETP
  208. C
  209. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  210. C Primale = IMATRI.LISPRI(4) = 'P'
  211. C -> IMATRI.LIZAFM(1,4) = RP
  212. C
  213. C-----------------------------------------------
  214. IMATRI.LIZAFM(1,1) = RR
  215. IMATRI.LIZAFM(1,2) = RUX
  216. IMATRI.LIZAFM(1,3) = RUY
  217. IMATRI.LIZAFM(1,4) = RUZ
  218. IMATRI.LIZAFM(1,5) = RP
  219. C-----------------------------------------------
  220. IMATRI.LIZAFM(1,6) = UXR
  221. IMATRI.LIZAFM(1,7) = UXUX
  222. IMATRI.LIZAFM(1,8) = UXUY
  223. IMATRI.LIZAFM(1,9) = UXUZ
  224. IMATRI.LIZAFM(1,10) = UXP
  225. C-----------------------------------------------
  226. IMATRI.LIZAFM(1,11) = UYR
  227. IMATRI.LIZAFM(1,12) = UYUX
  228. IMATRI.LIZAFM(1,13) = UYUY
  229. IMATRI.LIZAFM(1,14) = UYUZ
  230. IMATRI.LIZAFM(1,15) = UYP
  231. C------------------------------------------------
  232. IMATRI.LIZAFM(1,16) = UZR
  233. IMATRI.LIZAFM(1,17) = UZUX
  234. IMATRI.LIZAFM(1,18) = UZUY
  235. IMATRI.LIZAFM(1,19) = UZUZ
  236. IMATRI.LIZAFM(1,20) = UZP
  237. C-----------------------------------------------
  238. IMATRI.LIZAFM(1,21) = RETR
  239. IMATRI.LIZAFM(1,22) = RETUX
  240. IMATRI.LIZAFM(1,23) = RETUY
  241. IMATRI.LIZAFM(1,24) = RETUZ
  242. IMATRI.LIZAFM(1,25) = RETP
  243. C-----------------------------------------------
  244. C
  245. DO ICEN = 1, NBEL, 1
  246. RHO=MPRN.VPOCHA(ICEN,1)
  247. UX=MPVN.VPOCHA(ICEN,1)
  248. UY=MPVN.VPOCHA(ICEN,2)
  249. UZ=MPVN.VPOCHA(ICEN,3)
  250. GAMMA=MPGAM.VPOCHA(ICEN,1)
  251. C
  252. C********AB.AM(IFAC,IPRIM,IDUAL)
  253. C A = nom de l'inconnu duale (Ro,rUX,rUY,RET)
  254. C B = nom de l'inconnu primale (Ro,UX,UY,P)
  255. C
  256. C******* Dual RN
  257. C
  258. RR.AM(ICEN,1,1) = 1.0D0
  259. RUX.AM(ICEN,1,1) = 0.0D0
  260. RUY.AM(ICEN,1,1) = 0.0D0
  261. RUZ.AM(ICEN,1,1) = 0.0D0
  262. RP.AM(ICEN,1,1) = 0.0D0
  263. C
  264. C******* Dual RUXN
  265. C
  266. UXR.AM(ICEN,1,1) = UX
  267. UXUX.AM(ICEN,1,1) = RHO
  268. UXUY.AM(ICEN,1,1) = 0.0D0
  269. UXUZ.AM(ICEN,1,1) = 0.0D0
  270. UXP.AM(ICEN,1,1) = 0.0D0
  271. C
  272. C******* Dual RUYN
  273. C
  274. UYR.AM(ICEN,1,1) = UY
  275. UYUX.AM(ICEN,1,1) = 0.0D0
  276. UYUY.AM(ICEN,1,1) = RHO
  277. UYUZ.AM(ICEN,1,1) = 0.0D0
  278. UYP.AM(ICEN,1,1) = 0.0D0
  279. C
  280. C******* Dual RUZN
  281. C
  282. UZR.AM(ICEN,1,1) = UZ
  283. UZUX.AM(ICEN,1,1) = 0.0D0
  284. UZUY.AM(ICEN,1,1) = 0.0D0
  285. UZUZ.AM(ICEN,1,1) = RHO
  286. UZP.AM(ICEN,1,1) = 0.0D0
  287. C
  288. C********Dual RETN
  289. C
  290. RETR.AM(ICEN,1,1) = 0.5D0*((UX*UX)+(UY*UY)+(UZ*UZ))
  291. RETUX.AM(ICEN,1,1) = RHO*UX
  292. RETUY.AM(ICEN,1,1) = RHO*UY
  293. RETUZ.AM(ICEN,1,1) = RHO*UZ
  294. RETP.AM(ICEN,1,1) = 1.0D0/(GAMMA - 1.0D0)
  295. C
  296. ENDDO
  297. C
  298. SEGDES MPVN
  299. SEGDES MPRN
  300. SEGDES MPGAM
  301. C
  302. SEGDES MATRIK
  303. SEGDES IMATRI
  304. C
  305. SEGINI RR , RUX , RUY , RUZ, RP ,
  306. & UXR , UXUX , UXUY , UXUZ, UXP ,
  307. & UYR , UYUX , UYUY , UYUZ, UYP ,
  308. & UZR , UZUX , UZUY , UZUZ, UZP ,
  309. & RETR , RETUX , RETUY , RETUZ, RETP
  310. C
  311. RETURN
  312. END
  313.  
  314.  
  315.  
  316.  
  317.  
  318.  
  319.  
  320.  
  321.  

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