C TOMIXT    SOURCE    PV        05/04/25    21:15:48     5077           
      SUBROUTINE TOMIXT(MTAB2,MTAB3,MTAB5,IPT1,MPOVA1,
     &                  MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE)
C--------------------------------------------------------------------
C Evaluation des propriétés du mélange pour une composition donnée
C--------------------------------------------------------------------
C
C E/  MTAB2  : Table INCO
C E/  MTAB3  : Table COMPONENT
C E/  MTAB5  : Table PHYSICAL
C E/  IPT1   : Maillage POI1 des CELLs
C E/  MPOVA1 : Densité du mélange
C  /S MPOVAA : Lambda du mélange
C  /S MPOVAB : Mu du mélange
C  /S MPOVAC : Nu du mélange
C  /S MPOVAD : Constante des gaz parfait des incondensables
C  /S MPOVAE : Densité des incondensables
C
C---------------------------------------------------------------------
C
C Langage : ESOPE + FORTRAN77
C Auteurs : F.DABBENE 01/00
C
C---------------------------------------------------------------------
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)                 
      CHARACTER*4  NOMTOT(1)
      CHARACTER*8  TYPE,TYPS,MOT1,MOT2,MOT3,MOT4,MOT5
      POINTEUR MPOVAA.MPOVAL,MPOVAB.MPOVAL,MPOVAC.MPOVAL,
     &         MPOVAD.MPOVAL,MPOVAE.MPOVAL
C
-INC CCOPTIO
-INC SMCHPOI
-INC SMELEME
-INC SMTABLE
C
      SEGACT IPT1
      NBCELL = IPT1.NUM(/2)
      SEGACT MTAB3
      IDIM3 = MTAB3.MLOTAB
      N    = NBCELL
      NC   = 1
      SEGINI MPOVAA
      SEGINI MPOVAB
      SEGINI MPOVAC
      SEGINI MPOVAD
      SEGINI MPOVAE
C MCHPO2 : Densité du constituant i
C MCHPO3 : Lambda du constituant i
C MCHPO4 : Mu du constituant i
      DO 40 I=1,IDIM3
         CALL ACEM(MTAB3,I,MOT1)
         MOT2 = 'R'   // MOT1(1:7)
         MOT3 = 'LBD' // MOT1(1:5)
         MOT4 = 'MU'  // MOT1(1:6)
         MOT5 = 'R_'  // MOT1(1:6)
C
         TYPE = 'CHPOINT '
         CALL ACMO(MTAB2,MOT2,TYPE,MCHPO2)
         IF (IERR.NE.0) RETURN
         CALL ACMF(MTAB5,MOT3,FLO3)
CC       CALL ACMO(MTAB5,MOT3,TYPE,MCHPO3)
         IF (IERR.NE.0) RETURN
         CALL ACMF(MTAB5,MOT4,FLO4)
CC       CALL ACMO(MTAB5,MOT4,TYPE,MCHPO4)
         IF (IERR.NE.0) RETURN
         CALL ACMF(MTAB5,MOT5,RI)
         IF (IERR.NE.0) RETURN
C
         CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1)
         IF (IGEO1.NE.IPT1) THEN
            INDIC  =  1
            NBCOMP = -1
            CALL QUEPOI(MCHPO2,IPT1,INDIC,NBCOMP,NOMTOT)
            IF (IERR.NE.0) RETURN
            CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1)
         ENDIF
CC       CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
C        IF (IGEO1.NE.IPT1) THEN
C           INDIC  =  1
C           NBCOMP = -1
C           CALL QUEPOI(MCHPO3,IPT1,INDIC,NBCOMP,NOMTOT)
C           IF (IERR.NE.0) RETURN
C           CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
C        ENDIF
C        CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1)
C        IF (IGEO1.NE.IPT1) THEN
C           INDIC  =  1
C           NBCOMP = -1
C           CALL QUEPOI(MCHPO4,IPT1,INDIC,NBCOMP,NOMTOT)
C           IF (IERR.NE.0) RETURN
C           CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1)
CC       ENDIF
C
         DO 30 J=1,NBCELL
            MPOVAA . VPOCHA(J,1) =  MPOVAA . VPOCHA(J,1) +
     &      MPOVA2 . VPOCHA(J,1) *  FLO3
CC   &      MPOVA2 . VPOCHA(J,1) *  MPOVA3 . VPOCHA(J,1)
            MPOVAB . VPOCHA(J,1) =  MPOVAB . VPOCHA(J,1) +
     &      MPOVA2 . VPOCHA(J,1) *  FLO4
CC   &      MPOVA2 . VPOCHA(J,1) *  MPOVA4 . VPOCHA(J,1)
            IF (MOT1(1:3) . NE . 'H2O') THEN
               MPOVAD . VPOCHA(J,1) =  MPOVAD . VPOCHA(J,1) +
     &         RI * MPOVA2 . VPOCHA(J,1)
               MPOVAE . VPOCHA(J,1) =  MPOVAE . VPOCHA(J,1) +
     &         MPOVA2 . VPOCHA(J,1)
            ENDIF
 30      CONTINUE
         SEGDES MPOVA2
CC       SEGDES MPOVA2,MPOVA3,MPOVA4
 40   CONTINUE
      DO 50 J=1,NBCELL
         MPOVAA . VPOCHA(J,1) = MPOVAA . VPOCHA(J,1)
     &                        / MPOVA1 . VPOCHA(J,1)
         MPOVAB . VPOCHA(J,1) = MPOVAB . VPOCHA(J,1)
     &                        / MPOVA1 . VPOCHA(J,1)
         MPOVAC . VPOCHA(J,1) = MPOVAB . VPOCHA(J,1)
     &                        / MPOVA1 . VPOCHA(J,1)
         MPOVAD . VPOCHA(J,1) = MPOVAD . VPOCHA(J,1)
     &                        / MPOVAE . VPOCHA(J,1)
          X1 = MPOVAA . VPOCHA(J,1)
          X2 = MPOVAB . VPOCHA(J,1)
          X3 = MPOVAC . VPOCHA(J,1)
          X4 = MPOVAD . VPOCHA(J,1)
          X5 = MPOVAE . VPOCHA(J,1)
          X6 = MPOVAE . VPOCHA(J,1)
 50   CONTINUE
C
      RETURN
      END