Esope source files

   1 : C TOMIXT    SOURCE    PV        05/04/25    21:15:48     5077           
   2 :       SUBROUTINE TOMIXT(MTAB2,MTAB3,MTAB5,IPT1,MPOVA1,
   3 :      &                  MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE)
   4 : C--------------------------------------------------------------------
   5 : C Evaluation des propriétés du mélange pour une composition donnée
   6 : C--------------------------------------------------------------------
   7 : C
   8 : C E/  MTAB2  : Table INCO
   9 : C E/  MTAB3  : Table COMPONENT
  10 : C E/  MTAB5  : Table PHYSICAL
  11 : C E/  IPT1   : Maillage POI1 des CELLs
  12 : C E/  MPOVA1 : Densité du mélange
  13 : C  /S MPOVAA : Lambda du mélange
  14 : C  /S MPOVAB : Mu du mélange
  15 : C  /S MPOVAC : Nu du mélange
  16 : C  /S MPOVAD : Constante des gaz parfait des incondensables
  17 : C  /S MPOVAE : Densité des incondensables
  18 : C
  19 : C---------------------------------------------------------------------
  20 : C
  21 : C Langage : ESOPE + FORTRAN77
  22 : C Auteurs : F.DABBENE 01/00
  23 : C
  24 : C---------------------------------------------------------------------
  25 :       IMPLICIT INTEGER(I-N)
  26 :       IMPLICIT REAL*8 (A-H,O-Z)                 
  27 :       CHARACTER*4  NOMTOT(1)
  28 :       CHARACTER*8  TYPE,TYPS,MOT1,MOT2,MOT3,MOT4,MOT5
  29 :       POINTEUR MPOVAA.MPOVAL,MPOVAB.MPOVAL,MPOVAC.MPOVAL,
  30 :      &         MPOVAD.MPOVAL,MPOVAE.MPOVAL
  31 : C
  32 : -INC CCOPTIO
  33 : -INC SMCHPOI
  34 : -INC SMELEME
  35 : -INC SMTABLE
  36 : C
  37 :       SEGACT IPT1
  38 :       NBCELL = IPT1.NUM(/2)
  39 :       SEGACT MTAB3
  40 :       IDIM3 = MTAB3.MLOTAB
  41 :       N    = NBCELL
  42 :       NC   = 1
  43 :       SEGINI MPOVAA
  44 :       SEGINI MPOVAB
  45 :       SEGINI MPOVAC
  46 :       SEGINI MPOVAD
  47 :       SEGINI MPOVAE
  48 : C MCHPO2 : Densité du constituant i
  49 : C MCHPO3 : Lambda du constituant i
  50 : C MCHPO4 : Mu du constituant i
  51 :       DO 40 I=1,IDIM3
  51 :          CALL ACEM(MTAB3,I,MOT1)
  52 :          MOT2 = 'R'   // MOT1(1:7)
  53 :          MOT3 = 'LBD' // MOT1(1:5)
  54 :          MOT4 = 'MU'  // MOT1(1:6)
  55 :          MOT5 = 'R_'  // MOT1(1:6)
  56 : C
  57 :          TYPE = 'CHPOINT '
  57 :          CALL ACMO(MTAB2,MOT2,TYPE,MCHPO2)
  58 :          IF (IERR.NE.0) RETURN
  58 :          CALL ACMF(MTAB5,MOT3,FLO3)
  58 : CC       CALL ACMO(MTAB5,MOT3,TYPE,MCHPO3)
  59 :          IF (IERR.NE.0) RETURN
  59 :          CALL ACMF(MTAB5,MOT4,FLO4)
  59 : CC       CALL ACMO(MTAB5,MOT4,TYPE,MCHPO4)
  60 :          IF (IERR.NE.0) RETURN
  60 :          CALL ACMF(MTAB5,MOT5,RI)
  61 :          IF (IERR.NE.0) RETURN
  62 : C
  62 :          CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1)
  63 :          IF (IGEO1.NE.IPT1) THEN
  64 :             INDIC  =  1
  65 :             NBCOMP = -1
  65 :             CALL QUEPOI(MCHPO2,IPT1,INDIC,NBCOMP,NOMTOT)
  66 :             IF (IERR.NE.0) RETURN
  66 :             CALL LRCHT(MCHPO2,MPOVA2,TYPS,IGEO1)
  67 :          ENDIF
  67 : CC       CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
  68 : C        IF (IGEO1.NE.IPT1) THEN
  69 : C           INDIC  =  1
  70 : C           NBCOMP = -1
  70 : C           CALL QUEPOI(MCHPO3,IPT1,INDIC,NBCOMP,NOMTOT)
  71 : C           IF (IERR.NE.0) RETURN
  71 : C           CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
  72 : C        ENDIF
  72 : C        CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1)
  73 : C        IF (IGEO1.NE.IPT1) THEN
  74 : C           INDIC  =  1
  75 : C           NBCOMP = -1
  75 : C           CALL QUEPOI(MCHPO4,IPT1,INDIC,NBCOMP,NOMTOT)
  76 : C           IF (IERR.NE.0) RETURN
  76 : C           CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1)
  77 : CC       ENDIF
  78 : C
  79 :          DO 30 J=1,NBCELL
  80 :             MPOVAA . VPOCHA(J,1) =  MPOVAA . VPOCHA(J,1) +
  81 :      &      MPOVA2 . VPOCHA(J,1) *  FLO3
  82 : CC   &      MPOVA2 . VPOCHA(J,1) *  MPOVA3 . VPOCHA(J,1)
  83 :             MPOVAB . VPOCHA(J,1) =  MPOVAB . VPOCHA(J,1) +
  84 :      &      MPOVA2 . VPOCHA(J,1) *  FLO4
  85 : CC   &      MPOVA2 . VPOCHA(J,1) *  MPOVA4 . VPOCHA(J,1)
  86 :             IF (MOT1(1:3) . NE . 'H2O') THEN
  87 :                MPOVAD . VPOCHA(J,1) =  MPOVAD . VPOCHA(J,1) +
  88 :      &         RI * MPOVA2 . VPOCHA(J,1)
  89 :                MPOVAE . VPOCHA(J,1) =  MPOVAE . VPOCHA(J,1) +
  90 :      &         MPOVA2 . VPOCHA(J,1)
  91 :             ENDIF
  92 :  30      CONTINUE
  93 :          SEGDES MPOVA2
  94 : CC       SEGDES MPOVA2,MPOVA3,MPOVA4
  95 :  40   CONTINUE
  96 :       DO 50 J=1,NBCELL
  97 :          MPOVAA . VPOCHA(J,1) = MPOVAA . VPOCHA(J,1)
  98 :      &                        / MPOVA1 . VPOCHA(J,1)
  99 :          MPOVAB . VPOCHA(J,1) = MPOVAB . VPOCHA(J,1)
 100 :      &                        / MPOVA1 . VPOCHA(J,1)
 101 :          MPOVAC . VPOCHA(J,1) = MPOVAB . VPOCHA(J,1)
 102 :      &                        / MPOVA1 . VPOCHA(J,1)
 103 :          MPOVAD . VPOCHA(J,1) = MPOVAD . VPOCHA(J,1)
 104 :      &                        / MPOVAE . VPOCHA(J,1)
 105 :           X1 = MPOVAA . VPOCHA(J,1)
 106 :           X2 = MPOVAB . VPOCHA(J,1)
 107 :           X3 = MPOVAC . VPOCHA(J,1)
 108 :           X4 = MPOVAD . VPOCHA(J,1)
 109 :           X5 = MPOVAE . VPOCHA(J,1)
 110 :           X6 = MPOVAE . VPOCHA(J,1)
 111 :  50   CONTINUE
 112 : C
 113 :       RETURN
 114 :       END
 115 :  
 116 :  
 117 :