Télécharger pre1.eso

Retour à la liste

Numérotation des lignes :

  1. C PRE1 SOURCE PV 20/03/31 21:15:13 10567
  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.  
  184. -INC PPARAM
  185. -INC CCOPTIO
  186. -INC SMCOORD
  187. segact mcoord
  188. C
  189. C**** Les mots clés
  190. C
  191. C
  192. C**** Initialisation des parametres d'erreur
  193. C
  194. MOTERR(1:40) = ' '
  195. C
  196. C**** Lecture ORDESP
  197. C
  198. ICOND = 1
  199. CALL LIRENT(ORDESP,ICOND,IRETOU)
  200. IF(IERR .NE. 0) GOTO 9999
  201. IF((ORDESP .LT. 1) .OR. (ORDESP .GT. 2))THEN
  202. C
  203. C******* Message d'erreur standard
  204. C -301 0 %m1:40
  205. C
  206. MOTERR(1:40) = 'ENTI1 = ??? '
  207. CALL ERREUR(-301)
  208. C
  209. C******* Message d'erreur standard
  210. C Entier valant: %i1
  211. C -2 0
  212. C
  213. INTERR(1) = ORDESP
  214. CALL ERREUR(-2)
  215. C
  216. C******* Message d'erreur standard
  217. C 21 2
  218. C Données incompatibles
  219. C
  220. CALL ERREUR(21)
  221. GOTO 9999
  222. ENDIF
  223. C
  224. C**** Lecture ORDTEM
  225. C
  226. ICOND = 1
  227. CALL LIRENT(ORDTEM,ICOND,IRETOU)
  228. IF(IERR .NE. 0) GOTO 9999
  229. IF((ORDTEM .LT. 1) .OR. (ORDTEM .GT. ORDESP))THEN
  230. C
  231. C******* Message d'erreur standard
  232. C -301 0 %m1:40
  233. C
  234. MOTERR(1:40) = 'ENTI2 = ??? '
  235. CALL ERREUR(-301)
  236. C
  237. C******* Message d'erreur standard
  238. C Entier valant: %i1
  239. C -2 0
  240. C
  241. INTERR(1) = ORDTEM
  242. CALL ERREUR(-2)
  243. C
  244. C******* Message d'erreur standard
  245. C 21 2
  246. C Données incompatibles
  247. C
  248. CALL ERREUR(21)
  249. GOTO 9999
  250. ENDIF
  251. C
  252. IF((ORDESP .EQ. 1) .AND. (ORDTEM .EQ. 1))THEN
  253. C
  254. C******* ORDESP = 1, ORDTEM = 1
  255. C
  256. CALL PRE11()
  257. ELSEIF((ORDESP .EQ. 2))THEN
  258. C
  259. C******* ORDESP = 2, ORDTEM = 1 ou 2
  260. C
  261. CALL PRE12(ORDTEM)
  262. ELSE
  263. C
  264. C******* Message d'erreur standard
  265. C 5 3
  266. C Erreur anormale. Contactez votre support
  267. C
  268. CALL ERREUR(5)
  269. ENDIF
  270. C
  271. 9999 CONTINUE
  272. C
  273. RETURN
  274. END
  275.  
  276.  
  277.  
  278.  
  279.  
  280.  
  281.  
  282.  
  283.  

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