Télécharger cli191.eso

Retour à la liste

Numérotation des lignes :

  1. C CLI191 SOURCE PV 11/03/07 21:15:31 6885
  2. SUBROUTINE CLI191(MELEMF,MELEMC,MELECB,MELEFC,INORM,ICHPVO,ICHPSU,
  3. & IROC,IVITC,IPC,IGAMC,ICHLIM,ICHRES,ICHRLI)
  4. C************************************************************************
  5. C
  6. C PROJET : CASTEM 2000
  7. C
  8. C NOM : CLI191
  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 : S. KUDRIAKOV, 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. IMPLICIT REAL*8(A-H,O-Z)
  30. -INC CCOPTIO
  31. -INC SMLMOTS
  32. -INC SMELEME
  33. POINTEUR MELEFC.MELEME
  34. -INC SMLENTI
  35. POINTEUR MLEMC.MLENTI, MLEMCB.MLENTI,MLEMF.MLENTI
  36. -INC SMCHPOI
  37. POINTEUR MPNORM.MPOVAL, MPVOL.MPOVAL, MPSURF.MPOVAL, MPRC.MPOVAL,
  38. & MPVC.MPOVAL, MPPC.MPOVAL, MPGAMC.MPOVAL, MPLIM.MPOVAL,
  39. & MPRES.MPOVAL, MPRLI.MPOVAL
  40. C
  41. C**** Variables de COOPTIO
  42. C
  43. C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  44. C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  45. C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  46. C & ,IECHO, IIMPI, IOSPI
  47. C & ,IDIM, IFICLE, IPREFI
  48. C & ,MCOORD
  49. C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  50. C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  51. C & ,NORINC,NORVAL,NORIND,NORVAD
  52. C & ,NUCROU, IPSAUV
  53. C
  54. INTEGER MELEMF,MELEMC,MELECB,INORM,ICHPVO,ICHPSU, IROC,IVITC,IPC
  55. & ,IGAMC,ICHLIM,ICHRES,ICHRLI,ICEL,NFAC,IFAC
  56. & ,NGF,NGC,NLF,NLC,NLCB
  57. REAL*8 VOLU,SURF,GAMC,CNX,CNY,CNZ,CTX,CTY,CTZ,GM1,USGM1
  58. & ,CT2X,CT2Y,CT2Z,UTF,UT2F
  59. & ,PC,PF,TF,ECIN,P,RHO,UN,UT,UT2,UX,UY,UZ
  60. & ,RHOC,UXC,UYC,UZC,UNC,RAIR,SF
  61.  
  62. CHARACTER*(8) TYPE
  63. C
  64. C
  65. C**** KRIPAD pour la correspondance global/local
  66. C
  67. CALL KRIPAD(MELEMC,MLEMC)
  68. C SEGINI MLEMC
  69. CALL KRIPAD(MELECB,MLEMCB)
  70. C SEGINI MLEMCB
  71. CALL KRIPAD(MELEMF,MLEMF)
  72. C SEGINI MLEMF
  73. C
  74. C**** CHPOINTs de la table DOMAINE
  75. C
  76. CALL LICHT(INORM,MPNORM,TYPE,ICEL)
  77. CALL LICHT(ICHPVO,MPVOL,TYPE,ICEL)
  78. CALL LICHT(ICHPSU,MPSURF,TYPE,ICEL)
  79. C
  80. C**** LICHT active les MPOVALs en *MOD
  81. C
  82. C SEGACT MPNORM*MOD
  83. C SEGACT MPOVSU*MOD
  84. C SEGACT MPOVOL*MOD
  85. C
  86. C
  87. C**** CHPOINTs des variables
  88. C
  89. CALL LICHT(IROC,MPRC,TYPE,ICEL)
  90. CALL LICHT(IVITC,MPVC,TYPE,ICEL)
  91. CALL LICHT(IPC,MPPC,TYPE,ICEL)
  92. CALL LICHT(IGAMC,MPGAMC,TYPE,ICEL)
  93. CALL LICHT(ICHLIM,MPLIM,TYPE,ICEL)
  94. CALL LICHT(ICHRES,MPRES,TYPE,ICEL)
  95. CALL LICHT(ICHRLI,MPRLI,TYPE,ICEL)
  96. C
  97. C SEGACT *MOD
  98. C SEGACT *MOD
  99. C SEGACT *MOD
  100. C SEGACT *MOD
  101. C SEGACT *MOD
  102. C SEGACT *MOD
  103. C SEGACT *MOD
  104. C
  105. C
  106. C**** Boucle sur le face pour le calcul des invariants de
  107. C Riemann et du flux
  108. C
  109. SEGACT MELEFC
  110. NFAC=MELEFC.NUM(/2)
  111. CNZ=0.0D0
  112. CTZ=0.0D0
  113. CT2X=0.0D0
  114. CT2Y=0.0D0
  115. CT2Z=0.0D0
  116. UZC=0.0D0
  117. DO IFAC=1,NFAC,1
  118. NGF=MELEFC.NUM(1,IFAC)
  119. NGC=MELEFC.NUM(2,IFAC)
  120. NLF=MLEMF.LECT(NGF)
  121. NLC=MLEMC.LECT(NGC)
  122. NLCB=MLEMCB.LECT(NGF)
  123. VOLU=MPVOL.VPOCHA(NLC,1)
  124. SURF=MPSURF.VPOCHA(NLF,1)
  125. C In CASTEM les normales sont sortantes
  126. CNX=-1*MPNORM.VPOCHA(NLF,1)
  127. CNY=-1*MPNORM.VPOCHA(NLF,2)
  128. IF(IDIM.EQ.2)THEN
  129. CTX=-1.0D0*CNY
  130. CTY=CNX
  131. ELSE
  132. CNZ=-1*MPNORM.VPOCHA(NLF,3)
  133. CTX=-1*MPNORM.VPOCHA(NLF,4)
  134. CTY=-1*MPNORM.VPOCHA(NLF,5)
  135. CTZ=-1*MPNORM.VPOCHA(NLF,6)
  136. CT2X=-1*MPNORM.VPOCHA(NLF,7)
  137. CT2Y=-1*MPNORM.VPOCHA(NLF,8)
  138. CT2Z=-1*MPNORM.VPOCHA(NLF,9)
  139. ENDIF
  140. GAMC=MPGAMC.VPOCHA(NLC,1)
  141. GM1=GAMC-1.0D0
  142. USGM1=1.0D0/GM1
  143. C Variables au centre
  144. PC=MPPC.VPOCHA(NLC,1)
  145. RHOC=MPRC.VPOCHA(NLC,1)
  146. UXC=MPVC.VPOCHA(NLC,1)
  147. UYC=MPVC.VPOCHA(NLC,2)
  148. IF(IDIM.EQ.3) UZC=MPVC.VPOCHA(NLC,3)
  149. UNC=(UXC*CNX)+(UYC*CNY)+(UZC*CNZ)
  150. C Variables in the tank (stagnation values)
  151. PF=MPLIM.VPOCHA(NLCB,1)
  152. RF=MPLIM.VPOCHA(NLCB,2)
  153. UTF=0.0D0
  154. UT2F=0.0D0
  155. C
  156. C******* Variables à l'interface
  157. C
  158. c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  159. C******* entropy in the tank
  160. SF = PF/(RF**GAMC)
  161. c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
  162. RHO=((2.0/(GAMC + 1.0))**USGM1)*RF
  163. P = SF * (RHO**GAMC)
  164. UN=(GAMC*P/RHO)**0.5
  165. UT=UTF
  166. UT2=UT2F
  167. C
  168. C******* On calcule U
  169. C
  170. UX=UN*CNX+UT*CTX+UT2*CT2X
  171. UY=UN*CNY+UT*CTY+UT2*CT2Y
  172. UZ=UN*CNZ+UT*CTZ+UT2*CT2Z
  173. C
  174. ECIN=0.5D0*((UX*UX)+(UY*UY)+(UZ*UZ))
  175. C
  176. C******* Densite, vitesse, pression sur le bord
  177. C
  178. MPRLI.VPOCHA(NLCB,1)=RHO
  179. MPRLI.VPOCHA(NLCB,2)=UX
  180. MPRLI.VPOCHA(NLCB,3)=UY
  181. IF(IDIM.EQ.3) MPRLI.VPOCHA(NLCB,4)=UZ
  182. MPRLI.VPOCHA(NLCB,IDIM+2)=P
  183. C
  184. C
  185. C******* Residuum (son SPG a le meme ordre que MELEFC)
  186. C
  187. MPRES.VPOCHA(IFAC,1)=RHO*UN*SURF/VOLU
  188. MPRES.VPOCHA(IFAC,2)=((RHO*UN*UX)+(P*CNX))*SURF/VOLU
  189. MPRES.VPOCHA(IFAC,3)=((RHO*UN*UY)+(P*CNY))*SURF/VOLU
  190. IF(IDIM.EQ.3)MPRES.VPOCHA(IFAC,4)=((RHO*UN*UZ)+
  191. & (P*CNZ))*SURF/VOLU
  192. MPRES.VPOCHA(IFAC,IDIM+2)=(RHO*UN*((GAMC*USGM1*P/RHO)+ECIN))
  193. & *SURF/VOLU
  194. ENDDO
  195. C
  196. SEGDES MELEFC
  197. C
  198. SEGSUP MLEMC
  199. SEGSUP MLEMCB
  200. SEGSUP MLEMF
  201. C
  202. SEGDES MPNORM
  203. SEGDES MPVOL
  204. SEGDES MPSURF
  205. SEGDES MPRC
  206. SEGDES MPPC
  207. SEGDES MPVC
  208. SEGDES MPGAMC
  209. SEGDES MPLIM
  210. SEGDES MPRES
  211. SEGDES MPRLI
  212. C
  213. 9999 CONTINUE
  214. RETURN
  215. END
  216.  
  217.  
  218.  
  219.  
  220.  
  221.  
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  

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