C QOND      SOURCE    CHAT      05/01/13    02:40:30     5004           
      SUBROUTINE QOND
C************************************************************************       
C
C OBJET   : To compute the condensed mass flux density of H2O at
C           the wall from a mixture of air and vapour
C
C SYNTAXE : M = QOND CP ALFAB ALFAT H TP PTOT XVAP XH2O <BETA>
C
C M     : CHPOINT SCAL (condensed mass flux density of H2O (kg/m2/s))
C CP    : FLOTTANT     (vapour specific heat (J/kg/K))
C ALFAB : FLOTTANT     (brownian diffusivity (m2/s))
C ALFAT : FLOTTANT     (thermal diffusivity (m2/s))
C H     : CHPOINT SCAL (heat exchange coefficient at the wall (W/m2/K))
C TP    : CHPOINT SCAL (wall temperature (K))
C PTOT  : CHPOINT SCAL (total pressure of the mixture (Pa))
C XVAP  : CHPOINT SCAL (vapour quality ((vapour mass)/(water mass)))
C XH2O  : CHPOINT SCAL                 ((water mass)/(mixture mass))
C BETA  : FLOTTANT     (coefficient (option, default beta=1))
C
C
C M = BETA * H/CP *(ALFAB/ALFAT)**.666 *Ln((PTOT-PSAT(TP))/(PTOT-PVAP))
C M = 0   IF PVAP < PSAT(TP)
C
C IMPORTANT:   All the CHPOINTs MUST! have the same SPG
C ---------    S.I. units MUST! be used
C              M will be defined on the same SPG of H,TP,....
C
C************************************************************************       
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)                                                 
      CHARACTER*8 TYPE
C***                                                                            
-INC CCOPTIO                                                                    
-INC SMCHPOI
      POINTEUR    MH.MCHPOI,  MTP.MCHPOI,MPTOT.MCHPOI,
     1         MXVAP.MCHPOI,MXH2O.MCHPOI,MFVAP.MCHPOI
      POINTEUR    LH.MSOUPO,  LTP.MSOUPO,LPTOT.MSOUPO,
     1         LXVAP.MSOUPO,LXH2O.MSOUPO,LFVAP.MSOUPO
      POINTEUR    NH.MPOVAL,  NTP.MPOVAL,NPTOT.MPOVAL,
     1         NXVAP.MPOVAL,NXH2O.MPOVAL,NFVAP.MPOVAL

C     Reading operands

      CALL LIRREE(CP,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIRREE(ALFAB,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIRREE(ALFAT,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT ',MH,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT ',MTP,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT ',MPTOT,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT ',MXVAP,1,IRET)
      IF(IRET.EQ.0) RETURN
      CALL LIROBJ('CHPOINT ',MXH2O,1,IRET)
      IF(IRET.EQ.0) RETURN
      SEGACT MXH2O
      BETA=1.0
      CALL LIRREE(BETA,0,IRET)

C     Verfying operands

      CALL LICHT(MH,NH,TYPE,IGEO)
      CALL LICHT(MTP,NTP,TYPE,IGEO0)
      IF(IGEO.NE.IGEO0) THEN
C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
         MOTERR(1: 8) = ' TP '
         CALL ERREUR(788)
        RETURN
      ENDIF
      CALL LICHT(MPTOT,NPTOT,TYPE,IGEO0)
      IF(IGEO.NE.IGEO0) THEN
C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
         MOTERR(1: 8) = ' PT '
         CALL ERREUR(788)
        RETURN
      ENDIF
      CALL LICHT(MXVAP,NXVAP,TYPE,IGEO0)
      IF(IGEO.NE.IGEO0) THEN
C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
         MOTERR(1: 8) = ' YVAP'
         CALL ERREUR(788)
        RETURN
      ENDIF
      CALL LICHT(MXH2O,NXH2O,TYPE,IGEO0)
      IF(IGEO.NE.IGEO0) THEN
C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
         MOTERR(1: 8) = ' YH2O '
         CALL ERREUR(788)
        RETURN
      ENDIF

C     Creating the chpoint result

      NC=1
      CALL CRCHPT(TYPE,IGEO,NC,MFVAP)
      CALL LICHT(MFVAP,NFVAP,TYPE,IGEO)
      N=NFVAP.VPOCHA(/1)

C     Computing the condensed mass flux

      CALL XQOND(CP,ALFAB,ALFAT,N,BETA,
     1           NH.VPOCHA,   NTP.VPOCHA,  NPTOT.VPOCHA,
     2           NXVAP.VPOCHA,NXH2O.VPOCHA,NFVAP.VPOCHA,
     3           IRET)


      IF(IRET.NE.0) THEN
C     Tache impossible. Probablement données erronées
         CALL ERREUR(26)
        RETURN
      ENDIF

      CALL ECROBJ('CHPOINT ',MFVAP)

      SEGDES NFVAP,MH,NH,NTP,NPTOT,NXVAP,NXH2O

      RETURN
100   FORMAT(2X,'Les CHPOINTs n''ont pas le meme SPG')
      END