Télécharger la plaquette Cast3M


fr en

Télécharger primit.eso

Retour à la liste

Numérotation des lignes : 

primit

  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