Sources Esope | Cast3M

qond
C QOND      SOURCE    CB215821  20/11/25    13:38:06     10792          
      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 &lt;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 &lt; 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 PPARAM
-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