Télécharger primit.eso

Retour à la liste

Numérotation des lignes :

  1. C PRIMIT SOURCE BECC 11/05/18 21:15:36 6973
  2. SUBROUTINE PRIMIT()
  3. C************************************************************************
  4. C
  5. C PROJET : CASTEM 2000
  6. C
  7. C NOM : PRIMIT (OPERATEUR GIBIANE)
  8. C
  9. C DESCRIPTION : Calcul des variables primitives (et du "gamma")
  10. C pour les sivantes models du gaz:
  11. C
  12. C 1) gaz parfait mono-espece ("calorically perfect")
  13. C
  14. C 2) gaz parfait multi-especes ("calorically perfect")
  15. C
  16. C 3) gas parfait mono/multi-especes ("thermally perfect")
  17. C
  18. C 5) Modele bifluide ("two-fluid flow, water + air")
  19. C
  20. C 6) Modele bifluide pour la combustion ("RDEM model")
  21. C
  22. C 7) Modele bifluide ("Ghost Fluid method for the poor")
  23. C
  24. C LANGAGE : FORTRAN 77 + ESOPE 2000 (avec extensions CISI)
  25. C
  26. C AUTEUR : A. BECCANTINI, DRN/DMT/SEMT/TTMF
  27. C
  28. C************************************************************************
  29. C
  30. C
  31. C APPELES (E/S) : LIRCHA, OPTLI, ERREUR
  32.  
  33. C
  34. C APPELES (Calcul) :
  35. C
  36. C PRIMIT
  37. C |
  38. C |----- PRIM1E ---- PR1ECA
  39. C |
  40. C |----- PRIMME ---- PRMECA
  41. C |
  42. C |----- PRIMI1 ---- PRIMI2 ---- PRITEM
  43. C |
  44. C |----- PRIJAC
  45. C |
  46. C |----- PRIM2F
  47. C |
  48. C |----- PRIDEM
  49. C |
  50. C |----- PRIGFM
  51. C
  52. C
  53. C
  54. C************************************************************************
  55. C
  56. C PHRASE D'APPEL (GIBIANE) : voir PRIMIE, PRIMME, PRIDET
  57. C
  58. C
  59. C************************************************************************
  60. C
  61. C HISTORIQUE (Anomalies et modifications éventuelles)
  62. C
  63. C HISTORIQUE : Créée le 12.1.98.
  64. C
  65. C Modif: 14.12.98 model gas parfait mono/multi-especes
  66. C ("thermally perfect")
  67. C
  68. C Modif: 20/11/05 model bifluid
  69. C
  70. C Modif: 06/11/09 Discrete Equation Method approach
  71. C
  72. C Modif: 30/04/10 Ghost Fluid Method for the poor
  73. C
  74. C************************************************************************
  75. C
  76. C
  77. IMPLICIT INTEGER(I-N)
  78. INTEGER IP,IRET, ICOND,NBOPT,NBOPT0
  79. PARAMETER(NBOPT=7)
  80. CHARACTER*(8) TOPTIO(NBOPT),MOTOP
  81. C
  82. DATA TOPTIO/'PERFMONO','PERFMULT','PERFTEMP','CONSPRIM',
  83. & 'TWOFLUID','DEM ','GFMP '/
  84. C
  85. C**** Les Includes
  86. C
  87. -INC CCOPTIO
  88. C
  89. C**** Lecture de l'option
  90. C
  91. ICOND = 0
  92. CALL LIRCHA(MOTOP,ICOND,IRET)
  93. IF(IRET.EQ.0) THEN
  94. CALL PRIMI0
  95. RETURN
  96. ENDIF
  97. NBOPT0 = NBOPT
  98. CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT0)
  99. IF(IP .EQ. 0) THEN
  100. C
  101. C******* Message d'erreur standard
  102. C 251 2
  103. C Tentative d'utilisation d'une option non implémentée
  104. C
  105. CALL ERREUR(251)
  106. ELSEIF(IP .EQ. 1)THEN
  107. C
  108. C******* Gaz 'calorically perfect' monoespece
  109. C
  110. CALL PRIM1E()
  111. IF(IERR .NE. 0) GOTO 9999
  112. ELSEIF(IP .EQ. 2)THEN
  113. C
  114. C******* Gaz 'calorically perfect' multiespeces
  115. C
  116. CALL PRIMME()
  117. IF(IERR .NE. 0) GOTO 9999
  118. ELSEIF(IP .EQ. 3)THEN
  119. C
  120. C******* Gaz 'thermally perfect' mono/multiespeces
  121. C
  122. CALL PRIMI1()
  123. IF(IERR .NE. 0) GOTO 9999
  124. ELSEIF(IP .EQ. 4)THEN
  125. C
  126. C******* Gaz 'calorically perfect', la matrice jacobienne des
  127. C variables conservatives par rapport aux variables primitives
  128. C
  129. C
  130. CALL PRIJAC()
  131. IF(IERR .NE. 0) GOTO 9999
  132. ELSEIF(IP .EQ. 5)THEN
  133. C
  134. C******* Modele bifluide
  135. C
  136. CALL PRIM2F()
  137. IF(IERR .NE. 0) GOTO 9999
  138. ELSEIF(IP .EQ. 6)THEN
  139. C
  140. C******* Modele bifluide
  141. C
  142. CALL PRIDEM()
  143. IF(IERR .NE. 0) GOTO 9999
  144. ELSEIF(IP .EQ. 7)THEN
  145. C
  146. C******* Ghost fluid method for the poor
  147. C
  148. CALL PRIGFM()
  149. IF(IERR .NE. 0) GOTO 9999
  150. ELSE
  151. C
  152. C******* Message d'erreur standard
  153. C 5 2
  154. C Erreur anormale.contactez votre support
  155. C
  156. CALL ERREUR(5)
  157. ENDIF
  158. C
  159. 9999 CONTINUE
  160. RETURN
  161. END
  162.  
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  

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