Télécharger deto1.eso

Retour à la liste

Numérotation des lignes :

  1. C DETO1 SOURCE CHAT 05/01/12 22:44:11 5004
  2. SUBROUTINE DETO1(MSOUPO,MCHPO1,MCHPO2,MCHPO3,
  3. 1 AN,AH2,AH2O,AN2,AO2,TINT,IINT,NBC,NBT)
  4. C-----------------------------------------------------------------------
  5. C Calcul des points CJ, ZND et AICC pour un mélange
  6. C-----------------------------------------------------------------------
  7. C
  8. C---------------------------
  9. C Parametres Entree/Sortie :
  10. C---------------------------
  11. C
  12. C E/ MSOUPO : MSOUPO du CHPO décrivant le mélange (segment actif)
  13. C /S MCHPO1 : Etat CJ
  14. C /S MCHPO2 : Etat ZND
  15. C /S MCHPO3 : Etat AICC
  16. C E/ AN(NBT) : Masse molaire des especes (en gramme)
  17. C E/ AH2(NBT) : Coeff du devt en T des propriétés de h2
  18. C E/ AH2O(NBT) : Coeff du devt en T des propriétés de h2o
  19. C E/ AN2(NBT) : Coeff du devt en T des propriétés de n2
  20. C E/ AO2(NBT) : Coeff du devt en T des propriétés de o2
  21. C E/ TINT : Température seuil (K)
  22. C E/ IINT : Nombre de coeff nécessaire à T donnée
  23. C E/ NBC : Nombre de constituants du mélange
  24. C E/ NBT : Nombre de coeff total
  25. C
  26. C----------------------
  27. C Variables en COMMON :
  28. C----------------------
  29. C
  30. C IFOMOD : cf CCOPTIO.INC
  31. C
  32. C-----------------------------------------------------------------------
  33. C
  34. C Langage : ESOPE + FORTRAN77
  35. C
  36. C Auteurs : F.DABBENE et E.STUDER 02/99
  37. C
  38. C-----------------------------------------------------------------------
  39. IMPLICIT INTEGER(I-N)
  40. IMPLICIT REAL*8 (A-H,O-Z)
  41. C
  42. -INC CCOPTIO
  43. -INC SMCHPOI
  44. C
  45. C- Construction du chapeau des CHAMPOINs résultats
  46. C
  47. MPOVAL = IPOVAL
  48. SEGACT MPOVAL
  49. C-------------------------------------------------> CJ
  50. NAT = 1
  51. NSOUPO = 1
  52. NC = 5
  53. N = VPOCHA(/1)
  54. SEGINI MCHPO1
  55. SEGINI MSOUP1
  56. SEGINI MPOVA1
  57. MCHPO1.JATTRI(1) = 2
  58. MCHPO1.IPCHP(1) = MSOUP1
  59. MCHPO1.IFOPOI = IFOMOD
  60. MSOUP1.NOCOMP(1) = 'RCJ '
  61. MSOUP1.NOCOMP(2) = 'TCJ '
  62. MSOUP1.NOCOMP(3) = 'PCJ '
  63. MSOUP1.NOCOMP(4) = 'TAUX'
  64. MSOUP1.NOCOMP(5) = 'VCJ '
  65. MSOUP1.IGEOC = IGEOC
  66. MSOUP1.IPOVAL = MPOVA1
  67. MSOUP1.NOHARM(1) = NIFOUR
  68. MSOUP1.NOHARM(2) = NIFOUR
  69. MSOUP1.NOHARM(3) = NIFOUR
  70. MSOUP1.NOHARM(4) = NIFOUR
  71. MSOUP1.NOHARM(5) = NIFOUR
  72. SEGDES MSOUP1
  73. SEGDES MCHPO1
  74. C-------------------------------------------------> ZND
  75. NC = 3
  76. SEGINI MCHPO2
  77. SEGINI MSOUP2
  78. SEGINI MPOVA2
  79. MCHPO2.JATTRI(1) = 2
  80. MCHPO2.IPCHP(1) = MSOUP2
  81. MCHPO2.IFOPOI = IFOMOD
  82. MSOUP2.NOCOMP(1) = 'RZND'
  83. MSOUP2.NOCOMP(2) = 'TZND'
  84. MSOUP2.NOCOMP(3) = 'PZND'
  85. MSOUP2.IGEOC = IGEOC
  86. MSOUP2.IPOVAL = MPOVA2
  87. MSOUP2.NOHARM(1) = NIFOUR
  88. MSOUP2.NOHARM(2) = NIFOUR
  89. MSOUP2.NOHARM(3) = NIFOUR
  90. SEGDES MSOUP2
  91. SEGDES MCHPO2
  92. C-------------------------------------------------> AICC
  93. SEGINI MCHPO3
  94. SEGINI MSOUP3
  95. SEGINI MPOVA3
  96. MCHPO3.JATTRI(1) = 2
  97. MCHPO3.IPCHP(1) = MSOUP3
  98. MCHPO3.IFOPOI = IFOMOD
  99. MSOUP3.NOCOMP(1) = 'RAIC'
  100. MSOUP3.NOCOMP(2) = 'TAIC'
  101. MSOUP3.NOCOMP(3) = 'PAIC'
  102. MSOUP3.IGEOC = IGEOC
  103. MSOUP3.IPOVAL = MPOVA3
  104. MSOUP3.NOHARM(1) = NIFOUR
  105. MSOUP3.NOHARM(2) = NIFOUR
  106. MSOUP3.NOHARM(3) = NIFOUR
  107. SEGDES MSOUP3
  108. SEGDES MCHPO3
  109. SEGDES MSOUPO
  110. C
  111. C- Pour chacun des points du MSOUPO calcul des états et construction
  112. C- des CHAMPOINs résultats, la code sous DETO2 étant en FORTRAN pur
  113. C
  114. DO 10 I=1,N
  115. CALL DETO2(VPOCHA(I,1),VPOCHA(I,2),VPOCHA(I,3),
  116. 1 VPOCHA(I,4),VPOCHA(I,5),VPOCHA(I,6),
  117. 2 MPOVA1.VPOCHA(I,1),MPOVA1.VPOCHA(I,2),MPOVA1.VPOCHA(I,3),
  118. 3 MPOVA1.VPOCHA(I,4),MPOVA1.VPOCHA(I,5),
  119. 4 MPOVA2.VPOCHA(I,1),MPOVA2.VPOCHA(I,2),MPOVA2.VPOCHA(I,3),
  120. 5 MPOVA3.VPOCHA(I,1),MPOVA3.VPOCHA(I,2),MPOVA3.VPOCHA(I,3),
  121. 6 AN,AH2,AH2O,AN2,AO2,TINT,IINT,NBC,NBT,IRET)
  122. IF (MPOVA1.VPOCHA(I,1) . LE . 0.D0) THEN
  123. C 912 2 m1:4 est négatif au %i1 ème point
  124. MOTERR(1:4) = 'RCJ '
  125. INTERR(1) = I
  126. CALL ERREUR(912)
  127. GOTO 20
  128. ENDIF
  129. IF (IRET.NE.0) THEN
  130. C 460 2 Pas de convergence dans les itérations internes
  131. CALL ERREUR(460)
  132. GOTO 20
  133. ENDIF
  134. 10 CONTINUE
  135. 20 CONTINUE
  136. SEGDES MPOVA1
  137. SEGDES MPOVA2
  138. SEGDES MPOVA3
  139. SEGSUP MPOVAL
  140. C
  141. RETURN
  142. END
  143.  
  144.  
  145.  

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