Télécharger prija1.eso

Retour à la liste

Numérotation des lignes :

  1. C PRIJA1 SOURCE PV 16/11/17 22:01:08 9180
  2. SUBROUTINE PRIJA1(MELEMC,MLMCON,MLMPRI,IRN,IVN,IPN,IGAMN,IJACO)
  3. C
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : PRIJA1
  9. C
  10. C DESCRIPTION : Voir PRIJAC
  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, GAMMA
  74. CHARACTER*8 TYPE
  75. C
  76. C**** LES INCLUDES
  77. C
  78. -INC CCOPTIO
  79. -INC SMCHPOI
  80. -INC SMELEME
  81. -INC SMLMOTS
  82. -INC SMLENTI
  83. POINTEUR MPVN.MPOVAL, MPRN.MPOVAL, MPGAM.MPOVAL
  84. POINTEUR MELEMC.MELEME
  85. POINTEUR RR.IZAFM, RUX.IZAFM, RUY.IZAFM, RP.IZAFM,
  86. & RUXR.IZAFM, RUXUX.IZAFM, RUXUY.IZAFM, RUXP.IZAFM,
  87. & RUYR.IZAFM, RUYUX.IZAFM, RUYUY.IZAFM, RUYP.IZAFM,
  88. & RETR.IZAFM, RETUX.IZAFM, RETUY.IZAFM, RETP.IZAFM
  89. POINTEUR MLMPRI.MLMOTS, MLMCON.MLMOTS
  90. C
  91. CALL LICHT(IVN,MPVN,TYPE,IGEOM)
  92. CALL LICHT(IRN,MPRN,TYPE,IGEOM)
  93. CALL LICHT(IGAMN,MPGAM,TYPE,IGEOM)
  94. C
  95. C SEGACT MPVN*MOD
  96. C SEGACT MPPN*MOD
  97. C SEGACT MPRN*MOD
  98. C
  99. C**** Maillage des inconnues primales = Maillage des inconnues primales
  100. C = MELEMC
  101. C
  102. C
  103. NRIGE = 7
  104. NMATRI = 1
  105. NKID = 9
  106. NKMT = 7
  107. C
  108. SEGINI MATRIK
  109. IJACO = MATRIK
  110. MATRIK.IRIGEL(1,1) = MELEMC
  111. MATRIK.IRIGEL(2,1) = MELEMC
  112. C
  113. C**** Matrice non symetrique
  114. C
  115. MATRIK.IRIGEL(7,1) = 2
  116. C
  117. NBME = 16
  118. NBSOUS = 1
  119. SEGINI IMATRI
  120. MATRIK.IRIGEL(4,1) = IMATRI
  121. C
  122. SEGACT MLMPRI
  123. IMATRI.LISPRI(1) = MLMPRI.MOTS(1)
  124. IMATRI.LISPRI(2) = MLMPRI.MOTS(2)
  125. IMATRI.LISPRI(3) = MLMPRI.MOTS(3)
  126. IMATRI.LISPRI(4) = MLMPRI.MOTS(4)
  127. IMATRI.LISPRI(5) = MLMPRI.MOTS(1)
  128. IMATRI.LISPRI(6) = MLMPRI.MOTS(2)
  129. IMATRI.LISPRI(7) = MLMPRI.MOTS(3)
  130. IMATRI.LISPRI(8) = MLMPRI.MOTS(4)
  131. IMATRI.LISPRI(9) = MLMPRI.MOTS(1)
  132. IMATRI.LISPRI(10) = MLMPRI.MOTS(2)
  133. IMATRI.LISPRI(11) = MLMPRI.MOTS(3)
  134. IMATRI.LISPRI(12) = MLMPRI.MOTS(4)
  135. IMATRI.LISPRI(13) = MLMPRI.MOTS(1)
  136. IMATRI.LISPRI(14) = MLMPRI.MOTS(2)
  137. IMATRI.LISPRI(15) = MLMPRI.MOTS(3)
  138. IMATRI.LISPRI(16) = MLMPRI.MOTS(4)
  139. SEGDES MLMPRI
  140. C
  141. SEGACT MLMCON
  142. IMATRI.LISDUA(1) = MLMCON.MOTS(1)
  143. IMATRI.LISDUA(2) = MLMCON.MOTS(1)
  144. IMATRI.LISDUA(3) = MLMCON.MOTS(1)
  145. IMATRI.LISDUA(4) = MLMCON.MOTS(1)
  146. IMATRI.LISDUA(5) = MLMCON.MOTS(2)
  147. IMATRI.LISDUA(6) = MLMCON.MOTS(2)
  148. IMATRI.LISDUA(7) = MLMCON.MOTS(2)
  149. IMATRI.LISDUA(8) = MLMCON.MOTS(2)
  150. IMATRI.LISDUA(9) = MLMCON.MOTS(3)
  151. IMATRI.LISDUA(10) = MLMCON.MOTS(3)
  152. IMATRI.LISDUA(11) = MLMCON.MOTS(3)
  153. IMATRI.LISDUA(12) = MLMCON.MOTS(3)
  154. IMATRI.LISDUA(13) = MLMCON.MOTS(4)
  155. IMATRI.LISDUA(14) = MLMCON.MOTS(4)
  156. IMATRI.LISDUA(15) = MLMCON.MOTS(4)
  157. IMATRI.LISDUA(16) = MLMCON.MOTS(4)
  158. SEGDES MLMCON
  159. C
  160. SEGACT MELEMC
  161. NBEL = MELEMC.NUM(/2)
  162. SEGDES MELEMC
  163. NBSOUS = 1
  164. NP = 1
  165. MP = 1
  166. SEGINI RR , RUX , RUY , RP ,
  167. & RUXR , RUXUX , RUXUY , RUXP ,
  168. & RUYR , RUYUX , RUYUY , RUYP ,
  169. & RETR , RETUX , RETUY , RETP
  170. C
  171. C**** Duale = IMATRI.LISDUA(1) = 'RN'
  172. C Primale = IMATRI.LISPRI(4) = 'P'
  173. C -> IMATRI.LIZAFM(1,4) = RP
  174. C
  175. IMATRI.LIZAFM(1,1) = RR
  176. IMATRI.LIZAFM(1,2) = RUX
  177. IMATRI.LIZAFM(1,3) = RUY
  178. IMATRI.LIZAFM(1,4) = RP
  179. IMATRI.LIZAFM(1,5) = RUXR
  180. IMATRI.LIZAFM(1,6) = RUXUX
  181. IMATRI.LIZAFM(1,7) = RUXUY
  182. IMATRI.LIZAFM(1,8) = RUXP
  183. IMATRI.LIZAFM(1,9) = RUYR
  184. IMATRI.LIZAFM(1,10) = RUYUX
  185. IMATRI.LIZAFM(1,11) = RUYUY
  186. IMATRI.LIZAFM(1,12) = RUYP
  187. IMATRI.LIZAFM(1,13) = RETR
  188. IMATRI.LIZAFM(1,14) = RETUX
  189. IMATRI.LIZAFM(1,15) = RETUY
  190. IMATRI.LIZAFM(1,16) = RETP
  191. C
  192. DO ICEN = 1, NBEL, 1
  193. RHO=MPRN.VPOCHA(ICEN,1)
  194. UX=MPVN.VPOCHA(ICEN,1)
  195. UY=MPVN.VPOCHA(ICEN,2)
  196. GAMMA=MPGAM.VPOCHA(ICEN,1)
  197. C
  198. C********AB.AM(IFAC,IPRIM,IDUAL)
  199. C A = nom de l'inconnu duale (Ro,rUX,rUY,RET)
  200. C B = nom de l'inconnu primale (Ro,UX,UY,P)
  201. C
  202. C******* Dual RN
  203. C
  204. RR.AM(ICEN,1,1) = 1.0D0
  205. RUX.AM(ICEN,1,1) = 0.0D0
  206. RUY.AM(ICEN,1,1) = 0.0D0
  207. RP.AM(ICEN,1,1) = 0.0D0
  208. C
  209. C******* Dual RUXN
  210. C
  211. RUXR.AM(ICEN,1,1) = UX
  212. RUXUX.AM(ICEN,1,1) = RHO
  213. RUXUY.AM(ICEN,1,1) = 0.0D0
  214. RUXP.AM(ICEN,1,1) = 0.0D0
  215. C
  216. C******* Dual RUYN
  217. C
  218. RUYR.AM(ICEN,1,1) = UY
  219. RUYUX.AM(ICEN,1,1) = 0.0D0
  220. RUYUY.AM(ICEN,1,1) = RHO
  221. RUYP.AM(ICEN,1,1) = 0.0D0
  222. C
  223. C********Dual RETN
  224. C
  225. RETR.AM(ICEN,1,1) = 0.5D0*((UX*UX)+(UY*UY))
  226. RETUX.AM(ICEN,1,1) = RHO*UX
  227. RETUY.AM(ICEN,1,1) = RHO*UY
  228. RETP.AM(ICEN,1,1) = 1.0D0/(GAMMA - 1.0D0)
  229. C
  230. ENDDO
  231. C
  232. SEGDES MPVN
  233. SEGDES MPRN
  234. SEGDES MPGAM
  235. C
  236. SEGDES MATRIK
  237. SEGDES IMATRI
  238. C
  239. SEGDES RR , RUX , RUY , RP ,
  240. & RUXR , RUXUX , RUXUY , RUXP ,
  241. & RUYR , RUYUX , RUYUY , RUYP ,
  242. & RETR , RETUX , RETUY , RETP
  243. C
  244. RETURN
  245. END
  246.  
  247.  
  248.  
  249.  
  250.  
  251.  
  252.  
  253.  
  254.  

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