Télécharger pre5.eso

Retour à la liste

Numérotation des lignes :

  1. C PRE5 SOURCE KK2000 14/04/10 21:15:32 8032
  2. SUBROUTINE PRE5()
  3. C************************************************************************
  4. C
  5. C PROJET : CASTEM 2000
  6. C
  7. C NOM : PRE5
  8. C
  9. C DESCRIPTION : Voir PRET
  10. C
  11. C Resolution des Equations de la modele bifluide
  12. C
  13. C Calculation of the left and right states
  14. C at the interface for eact t step
  15. C
  16. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  17. C
  18. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  19. C Modified for two-fluid flow by
  20. C Jose R. Garcia Cascales
  21. C
  22. C************************************************************************
  23. C
  24. C
  25. C APPELES (Outils) : LIRENT, ERREUR
  26. C
  27. C APPELES (Calcul) :
  28. C
  29. C PRE1-------------- PRE12F (1st order)
  30. C |
  31. C |
  32. C --------------- PRE22F (2nd order)
  33. C
  34. C************************************************************************
  35. C
  36. C PHRASE D'APPEL (GIBIANE) :
  37. C
  38. C 1) two-fluid flow
  39. C
  40. C a) 1st order
  41. C
  42. C IALPHF, IUVF, IULF, IPF, ITVF, ITLF, IRVF, IRLF = 'PRET'
  43. C 'MULTIFLU' ORDESP ORDTPS $DOMA
  44. C VCH1 VCH2 VCH3 VCH4 VCH5 VCH6 VCH7 VCH8 ;
  45. C
  46. C CHAM1S CHAM2S CHAM3S CHAM4S CHAM5S CHAM6S CHAM7S CHAM8S = 'PRET'
  47. C 'MULTIFLU' ENTI1 ENTI2 TAB1
  48. C CHPO1 CHPO2 CHPO3 CHPO4 CHPO5 CHPO6 CHPO7 CHPO8 ;
  49. C
  50. C
  51. C ENTREES :
  52. C
  53. C 'MULTIFLU' : key name (two phase flow)
  54. C
  55. C ENTI1 : ordre en espace (=1)
  56. C
  57. C ENTI2 : ordre en temps (=1 si ENTI1 = 1)
  58. C
  59. C TAB1 : TABLE de sous type "DOMAINE"
  60. C
  61. C CHPO1 : CHPOINT "CENTRE" void fraction
  62. C
  63. C CHPO2 : CHPOINT "CENTRE" vapour velocity
  64. C (NDIM components)
  65. C
  66. C CHPO3 : CHPOINT "CENTRE" liquid velocity
  67. C (NDIM components)
  68. C
  69. C CHPO4 : CHPOINT "CENTRE" pressure
  70. C
  71. C CHPO5 : CHPOINT "CENTRE" vapour temperature
  72. C
  73. C CHPO6 : CHPOINT "CENTRE" liquid temperature
  74. C
  75. C CHPO7 : CHPOINT "CENTRE" vapour density
  76. C
  77. C CHPO8 : CHPOINT "CENTRE" liquid density
  78. C
  79. C SORTIES :
  80. C
  81. C CHAM1S : MCHAML "FACEL" void fraction
  82. C (a gauche et a droite de chaque face).
  83. C Only one component ('SCAL')
  84. C
  85. C CHAM2S : MCHAML "FACEL" vapour velocity and the
  86. C director cosines (n,t) in the corresponding face;
  87. C in the 2D case 6 composantes:
  88. C 'UVN' = normal velocity
  89. C 'UVT' = tangent velocity
  90. C 'NX' = n.x
  91. C 'NY' = n.y
  92. C 'TX' = t.x
  93. C 'TY' = t.y
  94. C
  95. C CHAM3S : MCHAML "FACEL" liquid velocity
  96. C in the 2D case 2 composantes:
  97. C 'ULN' = normal velocity
  98. C 'ULT' = tangent velocity
  99. C
  100. C CHAM4S : MCHAML "FACEL" pressure
  101. C Only one component ('SCAL')
  102. C
  103. C CHAM5S : MCHAML "FACEL" vapour temperature
  104. C Only one component ('SCAL')
  105. C
  106. C CHAM6S : MCHAML "FACEL" liquid temperature
  107. C Only one component ('SCAL')
  108. C
  109. C CHAM7S : MCHAML "FACEL" vapour density
  110. C Only one component ('SCAL')
  111. C
  112. C CHAM8S : MCHAML "FACEL" liquid density
  113. C Only one component ('SCAL')
  114. C
  115. C************************************************************************
  116. C
  117. C HISTORIQUE (Anomalies et modifications éventuelles)
  118. C
  119. C HISTORIQUE : Créée le 11.6.98.(Modified February 26th, 2002)
  120. C
  121. C************************************************************************
  122. C
  123. IMPLICIT INTEGER(I-N)
  124. IMPLICIT REAL*8(A-H,O-Z)
  125.  
  126. C
  127. C**** Les variables
  128. C
  129. INTEGER ICOND, IRETOU
  130. & ,ORDESP, ORDTEM
  131. C
  132. C**** Les Includes
  133. C
  134. -INC CCOPTIO
  135. C
  136. C**** Les mots clés
  137. C
  138. C
  139. C**** Initialisation des parametres d'erreur
  140. C
  141. MOTERR(1:40) = ' '
  142. C
  143. C**** Lecture ORDESP
  144. C
  145. ICOND = 1
  146. CALL LIRENT(ORDESP,ICOND,IRETOU)
  147. IF(IERR .NE. 0) GOTO 9999
  148. IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
  149. C
  150. C******* Message d'erreur standard
  151. C -301 0 %m1:40
  152. C
  153. MOTERR(1:40) = 'ENTI1 = ??? '
  154. CALL ERREUR(-301)
  155. C
  156. C******* Message d'erreur standard
  157. C Entier valant: %i1
  158. C -2 0
  159. C
  160. INTERR(1) = ORDESP
  161. CALL ERREUR(-2)
  162. C
  163. C******* Message d'erreur standard
  164. C 21 2
  165. C Données incompatibles
  166. C
  167. CALL ERREUR(21)
  168. GOTO 9999
  169. ENDIF
  170. C
  171. C**** Lecture ORDTEM
  172. C
  173. ICOND = 1
  174. CALL LIRENT(ORDTEM,ICOND,IRETOU)
  175. IF(IERR .NE. 0) GOTO 9999
  176. IF((ORDTEM .LT. 1) .OR. (ORDTEM .GT. ORDESP))THEN
  177. C
  178. C******* Message d'erreur standard
  179. C -301 0 %m1:40
  180. C
  181. MOTERR(1:40) = 'ENTI2 = ??? '
  182. CALL ERREUR(-301)
  183. C
  184. C******* Message d'erreur standard
  185. C Entier valant: %i1
  186. C -2 0
  187. C
  188. INTERR(1) = ORDTEM
  189. CALL ERREUR(-2)
  190. C
  191. C******* Message d'erreur standard
  192. C 21 2
  193. C Données incompatibles
  194. C
  195. CALL ERREUR(21)
  196. GOTO 9999
  197. ENDIF
  198. C
  199. IF((ORDESP .EQ. 1) .AND. (ORDTEM .EQ. 1))THEN
  200. C
  201. C******* ORDESP = 1, ORDTEM = 1
  202. C
  203. CALL PRE12F()
  204. ELSEIF((ORDESP .EQ. 2))THEN
  205. C
  206. C******* ORDESP = 2, ORDTEM = 1 ou 2
  207. C
  208. CALL PRE42F(ORDTEM)
  209. ELSE
  210. C
  211. C******* Message d'erreur standard
  212. C 5 3
  213. C Erreur anormale. Contactez votre support
  214. C
  215. CALL ERREUR(5)
  216. ENDIF
  217. C
  218. 9999 CONTINUE
  219. C
  220. RETURN
  221. END
  222.  
  223.  
  224.  
  225.  
  226.  
  227.  
  228.  

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