Télécharger qond.eso

Retour à la liste

Numérotation des lignes :

  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 <BETA>
  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.  
  45. CALL LIRREE(CP,1,IRET)
  46. IF(IRET.EQ.0) RETURN
  47. CALL LIRREE(ALFAB,1,IRET)
  48. IF(IRET.EQ.0) RETURN
  49. CALL LIRREE(ALFAT,1,IRET)
  50. IF(IRET.EQ.0) RETURN
  51. CALL LIROBJ('CHPOINT ',MH,1,IRET)
  52. IF(IRET.EQ.0) RETURN
  53. CALL LIROBJ('CHPOINT ',MTP,1,IRET)
  54. IF(IRET.EQ.0) RETURN
  55. CALL LIROBJ('CHPOINT ',MPTOT,1,IRET)
  56. IF(IRET.EQ.0) RETURN
  57. CALL LIROBJ('CHPOINT ',MXVAP,1,IRET)
  58. IF(IRET.EQ.0) RETURN
  59. CALL LIROBJ('CHPOINT ',MXH2O,1,IRET)
  60. IF(IRET.EQ.0) RETURN
  61. SEGACT MXH2O
  62. BETA=1.0
  63. CALL LIRREE(BETA,0,IRET)
  64.  
  65. C Verfying operands
  66.  
  67. CALL LICHT(MH,NH,TYPE,IGEO)
  68. CALL LICHT(MTP,NTP,TYPE,IGEO0)
  69. IF(IGEO.NE.IGEO0) THEN
  70. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  71. MOTERR(1: 8) = ' TP '
  72. CALL ERREUR(788)
  73. RETURN
  74. ENDIF
  75. CALL LICHT(MPTOT,NPTOT,TYPE,IGEO0)
  76. IF(IGEO.NE.IGEO0) THEN
  77. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  78. MOTERR(1: 8) = ' PT '
  79. CALL ERREUR(788)
  80. RETURN
  81. ENDIF
  82. CALL LICHT(MXVAP,NXVAP,TYPE,IGEO0)
  83. IF(IGEO.NE.IGEO0) THEN
  84. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  85. MOTERR(1: 8) = ' YVAP'
  86. CALL ERREUR(788)
  87. RETURN
  88. ENDIF
  89. CALL LICHT(MXH2O,NXH2O,TYPE,IGEO0)
  90. IF(IGEO.NE.IGEO0) THEN
  91. C Indice %m1:8 : L'objet %m9:16 n'a pas le bon support géométrique
  92. MOTERR(1: 8) = ' YH2O '
  93. CALL ERREUR(788)
  94. RETURN
  95. ENDIF
  96.  
  97. C Creating the chpoint result
  98.  
  99. NC=1
  100. CALL CRCHPT(TYPE,IGEO,NC,MFVAP)
  101. CALL LICHT(MFVAP,NFVAP,TYPE,IGEO)
  102. N=NFVAP.VPOCHA(/1)
  103.  
  104. C Computing the condensed mass flux
  105.  
  106. CALL XQOND(CP,ALFAB,ALFAT,N,BETA,
  107. 1 NH.VPOCHA, NTP.VPOCHA, NPTOT.VPOCHA,
  108. 2 NXVAP.VPOCHA,NXH2O.VPOCHA,NFVAP.VPOCHA,
  109. 3 IRET)
  110.  
  111.  
  112. IF(IRET.NE.0) THEN
  113. C Tache impossible. Probablement données erronées
  114. CALL ERREUR(26)
  115. RETURN
  116. ENDIF
  117.  
  118. CALL ECROBJ('CHPOINT ',MFVAP)
  119.  
  120. SEGDES NFVAP,MH,NH,NTP,NPTOT,NXVAP,NXH2O
  121.  
  122. RETURN
  123. 100 FORMAT(2X,'Les CHPOINTs n''ont pas le meme SPG')
  124. END
  125.  
  126.  
  127.  
  128.  
  129.  
  130.  

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