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 PPARAM
  88. -INC CCOPTIO
  89. C
  90. C**** Lecture de l'option
  91. C
  92. ICOND = 0
  93. CALL LIRCHA(MOTOP,ICOND,IRET)
  94. IF(IRET.EQ.0) THEN
  95. CALL PRIMI0
  96. RETURN
  97. ENDIF
  98. NBOPT0 = NBOPT
  99. CALL OPTLI(IP,TOPTIO,MOTOP,NBOPT0)
  100. IF(IP .EQ. 0) THEN
  101. C
  102. C******* Message d'erreur standard
  103. C 251 2
  104. C Tentative d'utilisation d'une option non implémentée
  105. C
  106. CALL ERREUR(251)
  107. ELSEIF(IP .EQ. 1)THEN
  108. C
  109. C******* Gaz 'calorically perfect' monoespece
  110. C
  111. CALL PRIM1E()
  112. IF(IERR .NE. 0) GOTO 9999
  113. ELSEIF(IP .EQ. 2)THEN
  114. C
  115. C******* Gaz 'calorically perfect' multiespeces
  116. C
  117. CALL PRIMME()
  118. IF(IERR .NE. 0) GOTO 9999
  119. ELSEIF(IP .EQ. 3)THEN
  120. C
  121. C******* Gaz 'thermally perfect' mono/multiespeces
  122. C
  123. CALL PRIMI1()
  124. IF(IERR .NE. 0) GOTO 9999
  125. ELSEIF(IP .EQ. 4)THEN
  126. C
  127. C******* Gaz 'calorically perfect', la matrice jacobienne des
  128. C variables conservatives par rapport aux variables primitives
  129. C
  130. C
  131. CALL PRIJAC()
  132. IF(IERR .NE. 0) GOTO 9999
  133. ELSEIF(IP .EQ. 5)THEN
  134. C
  135. C******* Modele bifluide
  136. C
  137. CALL PRIM2F()
  138. IF(IERR .NE. 0) GOTO 9999
  139. ELSEIF(IP .EQ. 6)THEN
  140. C
  141. C******* Modele bifluide
  142. C
  143. CALL PRIDEM()
  144. IF(IERR .NE. 0) GOTO 9999
  145. ELSEIF(IP .EQ. 7)THEN
  146. C
  147. C******* Ghost fluid method for the poor
  148. C
  149. CALL PRIGFM()
  150. IF(IERR .NE. 0) GOTO 9999
  151. ELSE
  152. C
  153. C******* Message d'erreur standard
  154. C 5 2
  155. C Erreur anormale.contactez votre support
  156. C
  157. CALL ERREUR(5)
  158. ENDIF
  159. C
  160. 9999 CONTINUE
  161. RETURN
  162. END
  163.  
  164.  
  165.  
  166.  
  167.  
  168.  
  169.  
  170.  
  171.  
  172.  
  173.  
  174.  
  175.  
  176.  
  177.  
  178.  

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