Esope source files

   1 : C QOND      SOURCE    CHAT      05/01/13    02:40:30     5004           
   2 :       SUBROUTINE QOND
   3 : C************************************************************************       
   4 : C
   5 : C OBJET   : To compute the condensed mass flux density of H2O at
   6 : C           the wall from a mixture of air and vapour
   7 : C
   8 : C SYNTAXE : M = QOND CP ALFAB ALFAT H TP PTOT XVAP XH2O 
   9 : C
  10 : C M     : CHPOINT SCAL (condensed mass flux density of H2O (kg/m2/s))
  11 : C CP    : FLOTTANT     (vapour specific heat (J/kg/K))
  12 : C ALFAB : FLOTTANT     (brownian diffusivity (m2/s))
  13 : C ALFAT : FLOTTANT     (thermal diffusivity (m2/s))
  14 : C H     : CHPOINT SCAL (heat exchange coefficient at the wall (W/m2/K))
  15 : C TP    : CHPOINT SCAL (wall temperature (K))
  16 : C PTOT  : CHPOINT SCAL (total pressure of the mixture (Pa))
  17 : C XVAP  : CHPOINT SCAL (vapour quality ((vapour mass)/(water mass)))
  18 : C XH2O  : CHPOINT SCAL                 ((water mass)/(mixture mass))
  19 : C BETA  : FLOTTANT     (coefficient (option, default beta=1))
  20 : C
  21 : C
  22 : C M = BETA * H/CP *(ALFAB/ALFAT)**.666 *Ln((PTOT-PSAT(TP))/(PTOT-PVAP))
  23 : C M = 0   IF PVAP < PSAT(TP)
  24 : C
  25 : C IMPORTANT:   All the CHPOINTs MUST! have the same SPG
  26 : C ---------    S.I. units MUST! be used
  27 : C              M will be defined on the same SPG of H,TP,....
  28 : C
  29 : C************************************************************************       
  30 :       IMPLICIT INTEGER(I-N)
  31 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  32 :       CHARACTER*8 TYPE
  33 : C***                                                                            
  34 : -INC CCOPTIO                                                                    
  35 : -INC SMCHPOI
  36 :       POINTEUR    MH.MCHPOI,  MTP.MCHPOI,MPTOT.MCHPOI,
  37 :      1         MXVAP.MCHPOI,MXH2O.MCHPOI,MFVAP.MCHPOI
  38 :       POINTEUR    LH.MSOUPO,  LTP.MSOUPO,LPTOT.MSOUPO,
  39 :      1         LXVAP.MSOUPO,LXH2O.MSOUPO,LFVAP.MSOUPO
  40 :       POINTEUR    NH.MPOVAL,  NTP.MPOVAL,NPTOT.MPOVAL,
  41 :      1         NXVAP.MPOVAL,NXH2O.MPOVAL,NFVAP.MPOVAL
  42 : 
  43 : C     Reading operands
  44 : 
  44 :       CALL LIRREE(CP,1,IRET)
  45 :       IF(IRET.EQ.0) RETURN
  45 :       CALL LIRREE(ALFAB,1,IRET)
  46 :       IF(IRET.EQ.0) RETURN
  46 :       CALL LIRREE(ALFAT,1,IRET)
  47 :       IF(IRET.EQ.0) RETURN
  47 :       CALL LIROBJ('CHPOINT ',MH,1,IRET)
  48 :       IF(IRET.EQ.0) RETURN
  48 :       CALL LIROBJ('CHPOINT ',MTP,1,IRET)
  49 :       IF(IRET.EQ.0) RETURN
  49 :       CALL LIROBJ('CHPOINT ',MPTOT,1,IRET)
  50 :       IF(IRET.EQ.0) RETURN
  50 :       CALL LIROBJ('CHPOINT ',MXVAP,1,IRET)
  51 :       IF(IRET.EQ.0) RETURN
  51 :       CALL LIROBJ('CHPOINT ',MXH2O,1,IRET)
  52 :       IF(IRET.EQ.0) RETURN
  53 :       SEGACT MXH2O
  54 :       BETA=1.0
  54 :       CALL LIRREE(BETA,0,IRET)
  55 : 
  56 : C     Verfying operands
  57 : 
  57 :       CALL LICHT(MH,NH,TYPE,IGEO)
  57 :       CALL LICHT(MTP,NTP,TYPE,IGEO0)
  58 :       IF(IGEO.NE.IGEO0) THEN
  59 : C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  60 :          MOTERR(1: 8) = ' TP '
  60 :          CALL ERREUR(788)
  61 :         RETURN
  62 :       ENDIF
  62 :       CALL LICHT(MPTOT,NPTOT,TYPE,IGEO0)
  63 :       IF(IGEO.NE.IGEO0) THEN
  64 : C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  65 :          MOTERR(1: 8) = ' PT '
  65 :          CALL ERREUR(788)
  66 :         RETURN
  67 :       ENDIF
  67 :       CALL LICHT(MXVAP,NXVAP,TYPE,IGEO0)
  68 :       IF(IGEO.NE.IGEO0) THEN
  69 : C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  70 :          MOTERR(1: 8) = ' YVAP'
  70 :          CALL ERREUR(788)
  71 :         RETURN
  72 :       ENDIF
  72 :       CALL LICHT(MXH2O,NXH2O,TYPE,IGEO0)
  73 :       IF(IGEO.NE.IGEO0) THEN
  74 : C     Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  75 :          MOTERR(1: 8) = ' YH2O '
  75 :          CALL ERREUR(788)
  76 :         RETURN
  77 :       ENDIF
  78 : 
  79 : C     Creating the chpoint result
  80 : 
  81 :       NC=1
  81 :       CALL CRCHPT(TYPE,IGEO,NC,MFVAP)
  81 :       CALL LICHT(MFVAP,NFVAP,TYPE,IGEO)
  82 :       N=NFVAP.VPOCHA(/1)
  83 : 
  84 : C     Computing the condensed mass flux
  85 : 
  85 :       CALL XQOND(CP,ALFAB,ALFAT,N,BETA,
  86 :      1           NH.VPOCHA,   NTP.VPOCHA,  NPTOT.VPOCHA,
  87 :      2           NXVAP.VPOCHA,NXH2O.VPOCHA,NFVAP.VPOCHA,
  88 :      3           IRET)
  89 : 
  90 : 
  91 :       IF(IRET.NE.0) THEN
  92 : C     Tache impossible. Probablement données erronées
  92 :          CALL ERREUR(26)
  93 :         RETURN
  94 :       ENDIF
  95 : 
  95 :       CALL ECROBJ('CHPOINT ',MFVAP)
  96 : 
  97 :       SEGDES NFVAP,MH,NH,NTP,NPTOT,NXVAP,NXH2O
  98 : 
  99 :       RETURN
 100 : 100   FORMAT(2X,'Les CHPOINTs n''ont pas le meme SPG')
 101 :       END
 102 :  
 103 :  
 104 :  
 105 :  
 106 :