Télécharger la plaquette Cast3M


fr en

Télécharger deto1.eso

Retour à la liste

Numérotation des lignes : 

deto1

  1. C DETO1     SOURCE    CB215821  25/04/23    21:15:14     12247          
  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  IFOUR : 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.  
  43. -INC PPARAM
  44. -INC CCOPTIO
  45. -INC SMCHPOI
  46. -INC SMCOORD
  47.  
  48.       DIMENSION AN(*),AH2(*),AH2O(*),AN2(*),AO2(*)
  49.  
  50. C
  51. C- Construction du chapeau des CHAMPOINs résultats
  52. C
  53.       MPOVAL = IPOVAL
  54.       SEGACT MPOVAL
  55. C-------------------------------------------------> CJ
  56.       NAT    = 1
  57.       NSOUPO = 1
  58.       NC     = 5
  59.       N      = VPOCHA(/1)
  60.       SEGINI MCHPO1
  61.       SEGINI MSOUP1
  62.       SEGINI MPOVA1
  63.       MCHPO1.JATTRI(1) = 2
  64.       MCHPO1.IPCHP(1)  = MSOUP1
  65.       MCHPO1.IFOPOI    = IFOUR
  66.       MSOUP1.NOCOMP(1) = 'RCJ '
  67.       MSOUP1.NOCOMP(2) = 'TCJ '
  68.       MSOUP1.NOCOMP(3) = 'PCJ '
  69.       MSOUP1.NOCOMP(4) = 'TAUX'
  70.       MSOUP1.NOCOMP(5) = 'VCJ '
  71.       MSOUP1.IGEOC     = IGEOC
  72.       MSOUP1.IPOVAL    = MPOVA1
  73.       MSOUP1.NOHARM(1) = NIFOUR
  74.       MSOUP1.NOHARM(2) = NIFOUR
  75.       MSOUP1.NOHARM(3) = NIFOUR
  76.       MSOUP1.NOHARM(4) = NIFOUR
  77.       MSOUP1.NOHARM(5) = NIFOUR
  78.       SEGDES MSOUP1
  79.       SEGDES MCHPO1
  80. C-------------------------------------------------> ZND
  81.       NC = 3
  82.       SEGINI MCHPO2
  83.       SEGINI MSOUP2
  84.       SEGINI MPOVA2
  85.       MCHPO2.JATTRI(1) = 2
  86.       MCHPO2.IPCHP(1)  = MSOUP2
  87.       MCHPO2.IFOPOI    = IFOUR
  88.       MSOUP2.NOCOMP(1) = 'RZND'
  89.       MSOUP2.NOCOMP(2) = 'TZND'
  90.       MSOUP2.NOCOMP(3) = 'PZND'
  91.       MSOUP2.IGEOC     = IGEOC
  92.       MSOUP2.IPOVAL    = MPOVA2
  93.       MSOUP2.NOHARM(1) = NIFOUR
  94.       MSOUP2.NOHARM(2) = NIFOUR
  95.       MSOUP2.NOHARM(3) = NIFOUR
  96.       SEGDES MSOUP2
  97.       SEGDES MCHPO2
  98. C-------------------------------------------------> AICC
  99.       SEGINI MCHPO3
  100.       SEGINI MSOUP3
  101.       SEGINI MPOVA3
  102.       MCHPO3.JATTRI(1) = 2
  103.       MCHPO3.IPCHP(1)  = MSOUP3
  104.       MCHPO3.IFOPOI    = IFOUR
  105.       MSOUP3.NOCOMP(1) = 'RAIC'
  106.       MSOUP3.NOCOMP(2) = 'TAIC'
  107.       MSOUP3.NOCOMP(3) = 'PAIC'
  108.       MSOUP3.IGEOC     = IGEOC
  109.       MSOUP3.IPOVAL    = MPOVA3
  110.       MSOUP3.NOHARM(1) = NIFOUR
  111.       MSOUP3.NOHARM(2) = NIFOUR
  112.       MSOUP3.NOHARM(3) = NIFOUR
  113.       SEGDES MSOUP3
  114.       SEGDES MCHPO3
  115.       SEGDES MSOUPO
  116. C
  117. C- Pour chacun des points du MSOUPO calcul des états et construction
  118. C- des CHAMPOINs résultats, la code sous DETO2 étant en FORTRAN pur
  119. C
  120.       DO 10 I=1,N
  121.          CALL DETO2(VPOCHA(I,1),VPOCHA(I,2),VPOCHA(I,3),
  122.      1      VPOCHA(I,4),VPOCHA(I,5),VPOCHA(I,6),
  123.      2      MPOVA1.VPOCHA(I,1),MPOVA1.VPOCHA(I,2),MPOVA1.VPOCHA(I,3),
  124.      3      MPOVA1.VPOCHA(I,4),MPOVA1.VPOCHA(I,5),
  125.      4      MPOVA2.VPOCHA(I,1),MPOVA2.VPOCHA(I,2),MPOVA2.VPOCHA(I,3),
  126.      5      MPOVA3.VPOCHA(I,1),MPOVA3.VPOCHA(I,2),MPOVA3.VPOCHA(I,3),
  127.      6      AN,AH2,AH2O,AN2,AO2,TINT,IINT,NBC,NBT,IRET)
  128.          IF (MPOVA1.VPOCHA(I,1) . LE . 0.D0) THEN
  129. C           912 2 m1:4 est négatif au %i1 ème point
  130.             MOTERR(1:4) = 'RCJ '
  131.             INTERR(1)   = I
  132.             CALL ERREUR(912)
  133.             GOTO 20
  134.          ENDIF
  135.          IF (IRET.NE.0) THEN
  136. C            460 2 Pas de convergence dans les itérations internes
  137.              CALL ERREUR(460)
  138.             GOTO 20
  139.          ENDIF
  140.  10   CONTINUE
  141.  20   CONTINUE
  142.       SEGDES MPOVA1
  143.       SEGDES MPOVA2
  144.       SEGDES MPOVA3
  145.       SEGSUP MPOVAL
  146. C
  147.       RETURN
  148.       END
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  

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