Télécharger pre1.eso

Retour à la liste

Numérotation des lignes :

  1. C PRE1 SOURCE CHAT 05/01/13 02:23:27 5004
  2. SUBROUTINE PRE1()
  3. C************************************************************************
  4. C
  5. C PROJET : CASTEM 2000
  6. C
  7. C NOM : PRE1
  8. C
  9. C DESCRIPTION : Voir PRET
  10. C
  11. C Resolution des Equations d'Euler, VF,
  12. C "cell-centred approach"
  13. C
  14. C Calcul des etats "gauche et droit" aux
  15. C interfaces pour le modele de gaz ideal
  16. C mono-espece "calorically perfect".
  17. C
  18. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  19. C
  20. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  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-------------- PRE11 (1er ordre en espace)
  30. C |
  31. C |
  32. C --------------- PRE12 (2eme ordre en espace)
  33. C
  34. C************************************************************************
  35. C
  36. C PHRASE D'APPEL (GIBIANE) :
  37. C
  38. C 1) gaz ideal mono-espece
  39. C
  40. C a) 1er ordre en espace (1er ordre en temps)
  41. C
  42. C CHAM1S CHAM2S CHAM3S CHAM4S = 'PRET' 'PERFMONO' ENTI1 ENTI2 MOD1
  43. C CHPO1 CHPO2 CHPO3 CHPO4 ;
  44. C
  45. C
  46. C ENTREES :
  47. C
  48. C 'PERFMONO' : mot clé (gaz ideal)
  49. C
  50. C ENTI1 : ordre en espace (=1)
  51. C
  52. C ENTI2 : ordre en temps (=1 si ENTI1 = 1)
  53. C
  54. C MOD1 : OBJET MODELE de TYPE NAVIER_STOKES
  55. C
  56. C CHPO1 : CHPOINT "CENTRE" contenant la masse volumique.
  57. C
  58. C CHPO2 : CHPOINT "CENTRE" contenant la vitesse
  59. C ( NDIM composantes);
  60. C
  61. C CHPO3 : CHPOINT "CENTRE" contenant la pression
  62. C
  63. C CHPO4 : CHPOINT "CENTRE" contenant les "gamma" du gaz
  64. C
  65. C SORTIES :
  66. C
  67. C CHAM1S : MCHAML "FACEL" contenant la masse
  68. C volumique
  69. C (a gauche et a droite de chaque face).
  70. C Une seule composante ('SCAL')
  71. C
  72. C CHAM2S : MCHAML "FACEL" contenant la vitesse et
  73. C les cosinus directeurs de (n,t) dans le
  74. C repere (x,y);
  75. C dans le cas 2D 6 composantes:
  76. C 'UN' = vitesse normale
  77. C 'UT' = vitesse tangentiel
  78. C 'NX' = n.x
  79. C 'NY' = n.y
  80. C 'TX' = t.x
  81. C 'TY' = t.y
  82. C
  83. C CHAM3S : MCHAML "FACEL" contenant la pression;
  84. C Une seule composante ('SCAL')
  85. C
  86. C CHAM4S : MCHAML "FACEL" contenant les "gamma" du gaz
  87. C Une seule composante ('SCAL')
  88. C
  89. C
  90. C b) 2eme ordre en espace
  91. C 1er ordre en temps ou 2eme ordre en temps
  92. C
  93. C CHAM1S CHAM2S CHAM3S CHAM4S = 'PRET' 'PERFMONO' ENTI1 ENTI2 TAB1
  94. C CHPO1 CHPO2 CHPO3
  95. C CHPO4 CHPO5 CHPO6
  96. C CHPO7 CHPO8 CHPO9
  97. C CHPO10 (FLOT1);
  98. C
  99. C
  100. C ENTREES :
  101. C
  102. C 'PERFMONO' : mot clé (gas ideal)
  103. C
  104. C ENTI1 : ordre en espace (=2)
  105. C
  106. C ENTI2 : ordre en temps (=1 ou 2)
  107. C
  108. C TAB1 : TABLE de sous type "DOMAINE"
  109. C
  110. C CHPO1 : CHPOINT "CENTRE" contenant la masse
  111. C volumique; une seule composante ('SCAL')
  112. C
  113. C CHPO2 : CHPOINT "CENTRE" contenant le gradient de
  114. C la masse volumique (2 composantes en 2D,
  115. C 'P1DX', 'P1DY').
  116. C
  117. C CHPO3 : CHPOINT "CENTRE" contenant le limiteur du
  118. C gradient de la masse volumique
  119. C (une seule composante 'P1 ')
  120. C
  121. C CHPO4 : CHPOINT "CENTRE" contenant la vitesse
  122. C (2 composantes en 2D, 'UX' ey 'UY' );
  123. C
  124. C CHPO5 : CHPOINT "CENTRE" contenant le gradient de la
  125. C vitesse (4 composantes en 2D, 'P1DX','P1DY'
  126. C 'P2DX','P2DY').
  127. C
  128. C CHPO6 : CHPOINT "CENTRE" contenant le limiteur du
  129. C gradient de la vitesse (2 composantes en 2D
  130. C 'P1', 'P2').
  131. C
  132. C CHPO7 : CHPOINT "CENTRE" contenant la pression;
  133. C une seule composante ('SCAL')
  134. C
  135. C CHPO8 : CHPOINT "CENTRE" contenant le gradient de
  136. C la pression (2 composantes en 2D,
  137. C 'P1DX', 'P1DY').
  138. C
  139. C CHPO9 : CHPOINT "CENTRE" contenant le limiteur de
  140. C gradient de la pression
  141. C (1 composantes, 'P1' ).
  142. C
  143. C CHPO10 : CHPOINT "CENTRE" contenant le "gamma" du gaz
  144. C
  145. C FLOT1 : FLOTTANT qui contient l'increment du temps
  146. C pour la prediction; pour avoir le deuxieme
  147. C ordre en temps (ENTI2 = 2)
  148. C FLOT1 = (increment du temps / 2)
  149. C
  150. C
  151. C SORTIES : voir le cas a)
  152. C
  153. C************************************************************************
  154. C
  155. C HISTORIQUE (Anomalies et modifications éventuelles)
  156. C
  157. C HISTORIQUE : Créée le 11.6.98.
  158. C
  159. C************************************************************************
  160. C
  161. C
  162. C**** Variables de COOPTIO
  163. C
  164. C INTEGER IPLLB, IERPER, IERMAX, IERR, INTERR
  165. C & ,IOTER, IOLEC, IOIMP, IOCAR, IOACQ
  166. C & ,IOPER, IOSGB, IOGRA, IOSAU, IORES
  167. C & ,IECHO, IIMPI, IOSPI
  168. C & ,IDIM
  169. C & ,MCOORD
  170. C & ,IFOMOD, NIFOUR, IFOUR, NSDPGE, IONIVE
  171. C & ,NGMAXY, IZROSF, ISOTYP, IOSCR,LTEXLU
  172. C & ,NORINC,NORVAL,NORIND,NORVAD
  173. C & ,NUCROU, IPSAUV
  174. C
  175. C**** Les variables
  176. C
  177. IMPLICIT INTEGER(I-N)
  178. INTEGER ICOND, IRETOU
  179. & ,ORDESP, ORDTEM
  180. C
  181. C**** Les Includes
  182. C
  183. -INC CCOPTIO
  184. C
  185. C**** Les mots clés
  186. C
  187. C
  188. C**** Initialisation des parametres d'erreur
  189. C
  190. MOTERR(1:40) = ' '
  191. C
  192. C**** Lecture ORDESP
  193. C
  194. ICOND = 1
  195. CALL LIRENT(ORDESP,ICOND,IRETOU)
  196. IF(IERR .NE. 0) GOTO 9999
  197. IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
  198. C
  199. C******* Message d'erreur standard
  200. C -301 0 %m1:40
  201. C
  202. MOTERR(1:40) = 'ENTI1 = ??? '
  203. CALL ERREUR(-301)
  204. C
  205. C******* Message d'erreur standard
  206. C Entier valant: %i1
  207. C -2 0
  208. C
  209. INTERR(1) = ORDESP
  210. CALL ERREUR(-2)
  211. C
  212. C******* Message d'erreur standard
  213. C 21 2
  214. C Données incompatibles
  215. C
  216. CALL ERREUR(21)
  217. GOTO 9999
  218. ENDIF
  219. C
  220. C**** Lecture ORDTEM
  221. C
  222. ICOND = 1
  223. CALL LIRENT(ORDTEM,ICOND,IRETOU)
  224. IF(IERR .NE. 0) GOTO 9999
  225. IF((ORDTEM .LT. 1) .OR. (ORDTEM .GT. ORDESP))THEN
  226. C
  227. C******* Message d'erreur standard
  228. C -301 0 %m1:40
  229. C
  230. MOTERR(1:40) = 'ENTI2 = ??? '
  231. CALL ERREUR(-301)
  232. C
  233. C******* Message d'erreur standard
  234. C Entier valant: %i1
  235. C -2 0
  236. C
  237. INTERR(1) = ORDTEM
  238. CALL ERREUR(-2)
  239. C
  240. C******* Message d'erreur standard
  241. C 21 2
  242. C Données incompatibles
  243. C
  244. CALL ERREUR(21)
  245. GOTO 9999
  246. ENDIF
  247. C
  248. IF((ORDESP .EQ. 1) .AND. (ORDTEM .EQ. 1))THEN
  249. C
  250. C******* ORDESP = 1, ORDTEM = 1
  251. C
  252. CALL PRE11()
  253. ELSEIF((ORDESP .EQ. 2))THEN
  254. C
  255. C******* ORDESP = 2, ORDTEM = 1 ou 2
  256. C
  257. CALL PRE12(ORDTEM)
  258. ELSE
  259. C
  260. C******* Message d'erreur standard
  261. C 5 3
  262. C Erreur anormale. Contactez votre support
  263. C
  264. CALL ERREUR(5)
  265. ENDIF
  266. C
  267. 9999 CONTINUE
  268. C
  269. RETURN
  270. END
  271.  
  272.  
  273.  
  274.  
  275.  
  276.  
  277.  
  278.  

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