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

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