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

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