Télécharger tocon0.eso

Retour à la liste

tocon0
C TOCON0    SOURCE    CB215821  20/11/25    13:41:05     10792          
      SUBROUTINE TOCON0
C
C--------------------------------------------------------------------
C Gestion de la condensation en paroi : évaluation du coef. d'échange
C thermique en paroi, du débit de condensation, du titre vapeur
C--------------------------------------------------------------------
C
C---------------------------
C Phrase d'appel (GIBIANE) :
C---------------------------
C
C TSAT1 XV3 ALAT0 ALAT3 ECHAN1 DVAP1 = VARI TOCOND TAB1 MAIL1 NOMMODE;
C ou
C TSAT1 XV3 = VARI TOCOND TAB1 ;
C
C------------------------
C Opérandes et résultat :
C------------------------
C
C  TAB1   : TABLE SET contenant l'ensemble des data 0D
C  MAIL1  : Maillage contenant l'ensemble des points faces associés 
C           au modele
C  NOMMODE: Mot définissant le modele de condensation utilisé
C  ALAT0  : Enthalpie liquide
C  ALAT3  : Débit de condensation / Flux thermique
C  ECHAN1 : Coefficient d'échange global en paroi
C  TSAT1  : Température de saturation
C  XV3    : Titre vapeur
C
C--------------------------------------------------------------------
C
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)     
      CHARACTER*4  NOMTOT(1)
      CHARACTER*8  TYPE,TYPS,MTYPI,MTYPR,MODCON(4)
      CHARACTER*72 CHARR,MOT1
      LOGICAL LOGI,LOGR
      DATA MODCON/'TAGAMI  ','UCHIDA  ','CHILTON ','COPAIN  '/
 
-INC PPARAM
-INC CCOPTIO
-INC SMCHPOI
-INC SMELEME
-INC SMTABLE
      POINTEUR MPOVA7.MPOVAL,MPOVA8.MPOVAL,MPOVA9.MPOVAL,MPOVAA.MPOVAL,
     &         MPOVAB.MPOVAL,MPOVAC.MPOVAL,MPOVAD.MPOVAL,MPOVAE.MPOVAL,
     &         MPOVAF.MPOVAL,MPOVAG.MPOVAL,
     &         MPOMIA.MPOVAL,MPOMIB.MPOVAL,MPOMIC.MPOVAL,MPOMID.MPOVAL,
     &         MPOMIE.MPOVAL,MPOMIF.MPOVAL,MPOMIG.MPOVAL,MPOMIH.MPOVAL,
     &         MPOMII.MPOVAL,MPOMIJ.MPOVAL,MPODV1.MPOVAL,
     &         MPOUT1.MPOVAL,MPOUT2.MPOVAL,MPOUT3.MPOVAL,MPOUT4.MPOVAL
C
      TYPS  = '        '
      TMAX  = 500.D0
      HMAX  = 5000.D0
      BETA  = 0.92D0
      GAMA  = 0.5D0
      TTAG  = 0.D0
      TAGA  = 0.D0
      TIME2 = 0.D0
C
C- Lecture et controles des données d'entrée,
C
C MTAB1 : Table SET
C Lecture de la table SET
      CALL LIRTAB('SET',MTAB1,1,IRETOU)
      IF (IERR.NE.0) RETURN
C Lecture du maillage de points envoyé
      CALL LIROBJ('MAILLAGE',IPT2,0,IRETOU)
      IF (IERR.NE.0) RETURN
C Lecture du modele de condensation
      ICOR = 0
      CALL LIRMOT(MODCON,4,ICOR,0)
      IF (IERR.NE.0) RETURN
C
C MTAB2 : Table INCO
C MTAB3 : Table COMPONENT
C MTAB4 : Table GEOINF
C MTAB5 : Table DATA
C MTAB6 : Table PHYSICAL
C MTAB7 : Table CELLAY
C
      MTYPI = 'MOT     '
      MTYPR = 'TABLE   '
      CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'INCO',LOGI,IRETI,
     &                  MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB2)
      IF (IERR.NE.0) RETURN
      CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'COMPONENT',LOGI,IRETI,
     &                  MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB3)
      IF (IERR.NE.0) RETURN
      CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'GEOINF',LOGI,IRETI,
     &                  MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB4)
      IF (IERR.NE.0) RETURN
      CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'DATA',LOGI,IRETI,
     &                  MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB5)
      IF (IERR.NE.0) RETURN
      CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'PHYSICAL',LOGI,IRETI,
     &                  MTYPR,IVALR,XVALR,CHARR,LOGR,MTAB6)
      IF (IERR.NE.0) RETURN
 
C --------------------------------------------------------------------C
C Données utilisées quelque soit le cas traité
C condensation ou non !
 
C
C IPT1 : Maillage POI1 des CELL
      TYPE = 'MAILLAGE'
      CALL ACMO(MTAB4,'CELL',TYPE,IPT1)
      IF (IERR.NE.0) RETURN
C
C MPOVA1 : Pression
C MPOVA2 : Température
C MPOVA3 : Densité
C MPOVA4 : Densité d'H2O
C MPOVA5 : Titre vapeur
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'PT',TYPE,MCHPO1)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPO1,MPOVA1,TYPS,IGEO1)
      IF (IGEO1.NE.IPT1) THEN
         INDIC  =  1
         NBCOMP = -1
         CALL QUEPOI(MCHPO1,IPT1,INDIC,NBCOMP,NOMTOT)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPO1,MPOVA1,TYPS,IGEO1)
      ENDIF
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'TGAS',TYPE,MCHPO2)
      IF (IERR.NE.0) RETURN
      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
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'RGAS',TYPE,MCHPO3)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
      IF (IGEO1.NE.IPT1) THEN
         INDIC  =  1
         NBCOMP = -1
         CALL QUEPOI(MCHPO3,IPT1,INDIC,NBCOMP,NOMTOT)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPO3,MPOVA3,TYPS,IGEO1)
      ENDIF
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'RH2O',TYPE,MCHPO4)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1)
      IF (IGEO1.NE.IPT1) THEN
         INDIC  =  1
         NBCOMP = -1
         CALL QUEPOI(MCHPO4,IPT1,INDIC,NBCOMP,NOMTOT)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPO4,MPOVA4,TYPS,IGEO1)
      ENDIF
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'XV',TYPE,MCHPO5)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPO5,MPOVA5,TYPS,IGEO1)
      IF (IGEO1.NE.IPT1) THEN
         INDIC  =  1
         NBCOMP = -1
         CALL QUEPOI(MCHPO5,IPT1,INDIC,NBCOMP,NOMTOT)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPO5,MPOVA5,TYPS,IGEO1)
      ENDIF
C
C RELINT : Facteur de relaxation pour les paramètres physiques internes
C
      MTYPI = 'MOT     '
      MTYPR = 'FLOTTANT'
      CALL ACCTAB(MTAB1,MTYPI,IVALI,XVALI,'RELINT',LOGI,IRETI,
     &                  MTYPR,IVALR,RELINT,CHARR,LOGR,IRETR)
      IF (IERR.NE.0) RETURN
 
C
C- Répartition vapeur/eau liquide au sein de chaque compartiment et
C- données thermodynamiques associées à la vapeur
C
C MPOVA6 : Titre vapeur actualisé
C MPOVA7 : Densité vapeur
C MPOVA8 : Pression partielle de vapeur
C MPOVA9 : Température de saturation
C
      CALL TOCON1(TMAX,RELINT,MPOVA2,MPOVA4,MPOVA5,
     &           MPOVA6,MPOVA7,MPOVA8,MPOVA9)
C
C 
C --------------------------------------------------------------------
C  DEBUT CALCUL DE LA CONDENSATION
C --------------------------------------------------------------------
C
C
      IF (IPT2.NE.0.AND.ICOR.NE.0) THEN
C
C IPT3 : Maillage POI1 des WALL
C IPT4 : Maillage SEG2 de CELLAYW (connectivites CELL-LAYERW)
      TYPE = 'MAILLAGE'
      CALL ACMO(MTAB4,'WALL',TYPE,IPT3)
      IF (IERR.NE.0) RETURN
      CALL ACMO(MTAB4,'CELLAYW',TYPE,IPT4)
      IF (IERR.NE.0) RETURN
C
C Données supplémentaires pour TAGAMI
C TIME2  : Temps physique en fin de pas de temps
C TTAG   : Temps physique caractéristique pour TAGAMI-UCHIDA
C TAGA   : Coefficient constant intervenant dans TAGAMI-UCHIDA
      IF (ICOR.EQ.1) THEN
         MTYPI = 'MOT     '
         MTYPR = 'FLOTTANT'
         CALL ACCTAB(MTAB2,MTYPI,IVALI,XVALI,'TIME2',LOGI,IRETI,
     &                     MTYPR,IVALR,TIME2,CHARR,LOGR,IRETR)
         IF (IERR.NE.0) RETURN
         CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'TTAG',LOGI,IRETI,
     &                     MTYPR,IVALR,TTAG,CHARR,LOGR,IRETR)
         IF (IERR.NE.0) RETURN
         CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'TAGA',LOGI,IRETI,
     &                     MTYPR,IVALR,TAGA,CHARR,LOGR,IRETR)
         IF (IERR.NE.0) RETURN
      ENDIF
C
C MPODV1 : Coeff de diffusion de la vapeur dans le mélange
C          (formule de RALOC, de Fuller avec loi de Blanc (FULLERB),
C          avec loi de Blanc et trace de vapeur dans le mélange (FULLERT).
C
      IF (ICOR.EQ.1 .OR. ICOR.EQ.2) THEN
        CALL ECROBJ('TABLE   ',MTAB1)
        CALL TODIFF
c      CALL LIROBJ('CHPOINT ',MCHDV1,1,IRET)
c      IF (IERR.NE.0) RETURN
c      CALL LRCHT(MCHDV1,MPODV1,TYPS,IGEO1)
c      IF (IGEO1.NE.IPT1) THEN
c         INDIC  =  1
c         NBCOMP = -1
c         CALL QUEPOI(MCHDV1,IPT1,INDIC,NBCOMP,NOMTOT)
c         IF (IERR.NE.0) RETURN
c         CALL LRCHT(MCHDV1,MPODV1,TYPS,IGEO1)
c      ENDIF
      ENDIF
      IF (ICOR.EQ.3 .OR. ICOR.EQ.4) THEN    
         MTYPI = 'MOT     '
         MTYPR = 'MOT     '
         CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'DIFVAP',LOGI,IRETI,
     &                     MTYPR,IVALR,XVALR,CHARR,LOGR,IRETR)
         IF (IERR.NE.0) RETURN
C
         IF (CHARR(1:5) .EQ. 'RALOC') THEN
            CALL ECROBJ('TABLE   ',MTAB1)
            CALL TODIFF
C Activation de mpova3 désactivé dans fullerb
            segact mpova3
         ENDIF
C
         IF (CHARR(1:7) .EQ. 'FULLERB') THEN
            CALL ECROBJ('TABLE   ',MTAB1)
            CALL FULLERB
C Activation de mpova3 désactivé dans fullerb
            segact mpova3
         ENDIF
C
         IF (CHARR(1:7) .EQ. 'FULLERT') THEN
            CALL ECROBJ('TABLE   ',MTAB1)
            CALL FULLERT
         ENDIF
C
         IF (CHARR(1:7) .EQ. 'FULLERX') THEN
            CALL ECROBJ('TABLE   ',MTAB1)
            CALL FULLERX
         ENDIF
C
         IF (CHARR(1:7) .EQ. 'FULLERY') THEN
            CALL ECROBJ('TABLE   ',MTAB1)
            CALL FULLERY
C Activation de mpova3 désactivé dans fullerb
            segact mpova3
         ENDIF
      ENDIF
      CALL LIROBJ('CHPOINT ',MCHDV1,1,IRET)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHDV1,MPODV1,TYPS,IGEO1)
      IF (IGEO1.NE.IPT1) THEN
         INDIC  =  1
         NBCOMP = -1
         CALL QUEPOI(MCHDV1,IPT1,INDIC,NBCOMP,NOMTOT)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHDV1,MPODV1,TYPS,IGEO1)
      ENDIF
C      ENDIF
C
C MPOVAA : Coeff d'echange par conduction
C MPOVAB : Température des murs
C MPOVAC : Aire de la surface de contact CELL/WALL
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB5,'ECHAC',TYPE,MCHPOA1)
      IF (IERR.NE.0) RETURN
      CALL REDUIR(MCHPOA1,IPT2,MCHPOA)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPOA,MPOVAA,TYPS,IGEO1)
      IF (IERR.NE.0) RETURN
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'TWAL',TYPE,MCHPOB)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPOB,MPOVAB,TYPS,IGEO1)
      IF (IGEO1.NE.IPT3) THEN
         INDIC  =  1
         NBCOMP = -1
         CALL QUEPOI(MCHPOB,IPT3,INDIC,NBCOMP,NOMTOT)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPOB,MPOVAB,TYPS,IGEO1)
      ENDIF
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB5,'AREAC',TYPE,MCHPOC1)
      IF (IERR.NE.0) RETURN
      CALL REDUIR(MCHPOC1,IPT2,MCHPOC)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPOC,MPOVAC,TYPS,IGEO1)
      IF (IERR.NE.0) RETURN
C
C- Evaluation du coefficient d'échange en paroi
C- et du débit de condensation
C
      IF (ICOR.EQ.1 .OR. ICOR.EQ.2) THEN
C
C MPOVAD : Facteur correctif du coef d'échange pour relaxaion béton
C
         TYPE = 'CHPOINT '
         CALL ACMO(MTAB6,'UCHI',TYPE,MCHPOD1)
         IF (IERR.NE.0) RETURN
         CALL REDUIR(MCHPOD1,IPT2,MCHPOD)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPOD,MPOVAD,TYPS,IGEO1)
         IF (IERR.NE.0) RETURN
C
C MPOUT1 : Coefficient d'échange global en paroi
C MPOUT2 : Débit de condensation / flux thermique
C MPOUT3 : Enthalpie liquide
C MPOUT4 : 0.
C
         CALL UCHIDA(ICOR,IPT1,IPT2,IPT3,IPT4,
     &               MPOVA7,MPOVA3,MPOVA4,MPOVAA,MPOVA1,
     &               MPOVA8,MPOVA9,MPOVA2,MPOVAB,MPOVAC,
     &               MPOVAD,HMAX,BETA,GAMA,TTAG,TIME2,TAGA,
     &               MPOUT1,MPOUT2,MPOUT3,MPOUT4)
      ENDIF
C
      IF (ICOR.EQ.3 .OR. ICOR.EQ.4) THEN
C
C MPOVAD : Ancien flux thermique global
C
         TYPE = 'CHPOINT '
         CALL ACMO(MTAB2,'QW',TYPE,MCHPOD1)
         IF (IERR.NE.0) RETURN
         CALL REDUIR(MCHPOD1,IPT2,MCHPOD)
         IF (IERR.NE.0) RETURN
         CALL LRCHT(MCHPOD,MPOVAD,TYPS,IGEO1)
         IF (IERR.NE.0) RETURN
C
         MTYPI = 'MOT     '
         MTYPR = 'FLOTTANT'
         CALL ACCTAB(MTAB6,MTYPI,IVALI,XVALI,'GRAVITY',LOGI,IRETI,
     &                     MTYPR,IVALR,G1,CHARR,LOGR,IRETR)
         IF (IERR.NE.0) RETURN
C
C MPOMIA : Lambda du mélange
C MPOMIB : Mu du mélange
C MPOMIC : Nu du mélange
C MPOMID : Constante des gaz parfait des incondensables
C MPOMIE : Densité des incondensables
C MPOMIG : Capacité thermique du mélange
C MPOMIH : Constante des gaz parfait du mélange
C MPOMII : Masse molaire des incondensables
C MPOMIJ : Masse molaire de la vapeur d'eau
C
         CALL TOMIXTC(MTAB2,MTAB3,MTAB6,IPT1,MPOVA3,
     &               MPOMIA,MPOMIB,MPOMIC,MPOMID,MPOMIE,MPOMIG,MPOMIH,
     &               MPOMII,MPOMIJ)
C
C MPOUT1 : Coefficient d'échange global en paroi
C MPOUT2 : Débit de condensation / flux thermique
C MPOUT3 : Enthalpie liquide
C MPOUT4 : Coefficient d'échange convectif en paroi
C
        IF (ICOR.EQ.3) THEN
         CALL CHILTO(IPT1,IPT2,IPT3,IPT4,
     &               MPOVA1,MPOVA2,MPOVA3,MPOVA7,MPOVA8,MPOVA9,
     &               MPOVAA,MPOVAB,MPOVAC,MPOVAD,
     &               MPOMIA,MPOMIC,MPOMID,MPODV1,
     &               G1,HMAX,
     &               MPOUT1,MPOUT2,MPOUT3,MPOUT4)
        ENDIF
        IF (ICOR.EQ.4) THEN
         CALL COPAIN(IPT1,IPT2,IPT3,IPT4,
     &               MPOVA1,MPOVA2,MPOVA3,MPOVA7,MPOVA8,MPOVA9,
     &               MPOVAA,MPOVAB,MPOVAC,MPOVAD,
     &               MPOMIA,MPOMIC,MPOMID,MPODV1,
     &               G1,BETA,HMAX,MPOMIG,MPOMIH,MPOMII,MPOMIJ,
     &               MPOUT1,MPOUT2,MPOUT3,MPOUT4)
        ENDIF
      ENDIF
C
C- Si le débit de condensation n'est pas nul,
C- Relaxation du coeff d'échange global et du débit de condensation
C
C MPOVAE : Ancien Coefficient d'échange global en paroi
C MPOVAF : Ancien Débit de condensation / flux thermique
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'ECHA',TYPE,MCHPOE1)
      IF (IERR.NE.0) RETURN
      CALL REDUIR(MCHPOE1,IPT2,MCHPOE)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPOE,MPOVAE,TYPS,IGEO1)
      IF (IERR.NE.0) RETURN
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'ALAT',TYPE,MCHPOF1)
      IF (IERR.NE.0) RETURN
      CALL REDUIR(MCHPOF1,IPT2,MCHPOF)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPOF,MPOVAF,TYPS,IGEO1)
      IF (IERR.NE.0) RETURN
C
      TYPE = 'CHPOINT '
      CALL ACMO(MTAB2,'HT',TYPE,MCHPOG1)
      IF (IERR.NE.0) RETURN
      CALL REDUIR(MCHPOG1,IPT2,MCHPOG)
      IF (IERR.NE.0) RETURN
      CALL LRCHT(MCHPOG,MPOVAG,TYPS,IGEO1)
      IF (IERR.NE.0) RETURN
C
      SEGACT IPT2
      NBLAY = IPT2.NUM(/2)
      SEGDES IPT2
      DO 10 I=1,NBLAY
       IF (MPOUT2.VPOCHA(I,1) .LT. 0.D0) THEN
         MPOUT1.VPOCHA(I,1) = RELINT * MPOUT1.VPOCHA(I,1)
     &                      + (1.D0 - RELINT) * MPOVAE.VPOCHA(I,1)
         MPOUT2.VPOCHA(I,1) = RELINT * MPOUT2.VPOCHA(I,1)
     &                      + (1.D0 - RELINT) * MPOVAF.VPOCHA(I,1)
         MPOUT4.VPOCHA(I,1) = RELINT * MPOUT4.VPOCHA(I,1)
     &                      + (1.D0 - RELINT) * MPOVAG.VPOCHA(I,1)
       ENDIF
 10   CONTINUE
C
C- Ecriture chpo Dvap
C
      CALL ECROBJ('CHPOINT',MCHDV1)
C
C- Ibid chpo contenant le coefficient d'échange global
C
      SEGINI, MCHPOI=MCHPO1
      NC = 1
      SEGINI MSOUPO
      IPCHP(1)  = MSOUPO
      NOCOMP(1) = 'SCAL'
      IGEOC     = IPT2
      IPOVAL    = MPOUT1
      NOHARM(1) = NIFOUR
      SEGDES MCHPOI,MSOUPO,MPOUT1
      CALL ECROBJ('CHPOINT',MCHPOI)
C
C- Ibid chpo contenant le débit de condensation
C
      SEGINI, MCHPOI=MCHPO1
      NC = 1
      SEGINI MSOUPO
      IPCHP(1)  = MSOUPO
      NOCOMP(1) = 'SCAL'
      IGEOC     = IPT2
      IPOVAL    = MPOUT2
      NOHARM(1) = NIFOUR
      SEGDES MCHPOI,MSOUPO,MPOUT2
      CALL ECROBJ('CHPOINT',MCHPOI)      
C
C- Ibid chpo contenant l'enthalpie liquide
C
      SEGINI, MCHPOI=MCHPO1
      NC = 1
      SEGINI MSOUPO
      IPCHP(1)  = MSOUPO
      NOCOMP(1) = 'SCAL'
      IGEOC     = IPT2
      IPOVAL    = MPOUT3
      NOHARM(1) = NIFOUR
      SEGDES MCHPOI,MSOUPO,MPOUT3
      CALL ECROBJ('CHPOINT',MCHPOI)
C
C- Ibid chpo contenant hconv
C
      SEGINI, MCHPOI=MCHPO1
      NC = 1
      SEGINI MSOUPO
      IPCHP(1)  = MSOUPO
      NOCOMP(1) = 'SCAL'
      IGEOC     = IPT2
      IPOVAL    = MPOUT4
      NOHARM(1) = NIFOUR
      SEGDES MCHPOI,MSOUPO,MPOUT4
      CALL ECROBJ('CHPOINT',MCHPOI)
C
C Menage des chpo utilisé pour le calcul de la condensation      
C
      SEGDES MPOVAA,MPOVAB,MPOVAC,MPOVAD,MPOVAE,MPOVAF,MPOVAG
 
      ENDIF
C --------------------------------------------------------------------C
C  FIN CALCUL DE LA CONDENSATION
C --------------------------------------------------------------------C
 
C
C- Création/remplissage du chpo contenant le titre vapeur
C
      SEGINI, MCHPOI=MCHPO1
      NC = 1
      SEGINI MSOUPO
      IPCHP(1)  = MSOUPO
      NOCOMP(1) = 'SCAL'
      IGEOC     = IPT1
      IPOVAL    = MPOVA6
      NOHARM(1) = NIFOUR
      SEGDES MCHPOI,MSOUPO,MPOVA6
      CALL ECROBJ('CHPOINT',MCHPOI)
C
C- Ibid chpo contenant Tsat
C
      SEGINI, MCHPOI=MCHPO1
      NC = 1
      SEGINI MSOUPO
      IPCHP(1)  = MSOUPO
      NOCOMP(1) = 'SCAL'
      IGEOC     = IPT1
      IPOVAL    = MPOVA9
      NOHARM(1) = NIFOUR
      SEGDES MCHPOI,MSOUPO,MPOVA9
      CALL ECROBJ('CHPOINT',MCHPOI)
C
C- Ménage
C
      SEGDES MPOVA1,MPOVA2,MPOVA3,MPOVA4,MPOVA5
      SEGSUP MPOVA7,MPOVA8
C
      RETURN
      END
 
 
 
 

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