Esope source files

   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
 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
 124 :             CALL ERREUR(912)
 125 :             GOTO 20
 126 :          ENDIF
 127 :          IF (IRET.NE.0) THEN
 128 : C            460 2 Pas de convergence dans les itérations internes
 128 :              CALL ERREUR(460)
 129 :             GOTO 20
 130 :          ENDIF
 131 :  10   CONTINUE
 132 :  20   CONTINUE
 133 :       SEGDES MPOVA1
 134 :       SEGDES MPOVA2
 135 :       SEGDES MPOVA3
 136 :       SEGSUP MPOVAL
 137 : C
 138 :       RETURN                                                                    
 139 :       END                                                                       
 140 :  
 141 :