Esope source files

   1 : C FLA001    SOURCE    SERRE     12/07/10    21:15:09     7431           
   2 :       FUNCTION FLA001(P,T,X,TRAV)
   3 : C------------------------------------------------------------------------
   4 : C Calcul de la densité du mélange en fonction de la pression, de la
   5 : C température et de la composition du mélange (on suppose que la loi des
   6 : C gaz parfaits s'applique aussi à la vapeur d'eau)
   7 : C------------------------------------------------------------------------
   8 : C
   9 : C---------------------------
  10 : C Parametres Entree/Sortie :
  11 : C---------------------------
  12 : C
  13 : C E/  P      : flottant : Pression de référence (Pa)
  14 : C E/  T      : flottant : Température de référence (K)
  15 : C E/  X(7)   : flottant : Fraction molaire de N2, O2, H2 et H2O
  16 : C  /S FLA001 : flottant : Densité du mélange (kg/m3)
  17 : C
  18 : C------------------------------
  19 : C Variables de TRAV utilisées :
  20 : C------------------------------
  21 : C
  22 : C E/  R      : Constante des gaz (= 8.314 J/mol/K)
  23 : C E/  M(7)   : Masse molaire des constituants du mélange (kg/mol) 
  24 : C
  25 : C------------------------------------------------------------------------
  26 : C
  27 : C Langage : ESOPE + FORTRAN 77
  28 : C
  29 : C Mise en oeuvre : H. Paillère (1997, TTMF)
  30 : C
  31 : C------------------------------------------------------------------------
  32 :       IMPLICIT INTEGER(I-N)
  33 :       IMPLICIT REAL*8 (A-H,O-Z)
  34 :       REAL*8 FLA001
  35 :       REAL*8 X(7)
  36 :       segment trav
  37 :            integer iKALP,iMODEL
  38 :            real*8 e,L,Lch,Dh,S,sp,Ck
  39 :            real*8 mc,Cpc
  40 :            real*8 g,R,deltah
  41 :            real*8 M(nbesp),cstmod(ncst)
  42 :            real*8 Cpi(nbesp),al
  43 :            real*8 eps_mh2,eps_dt,eps_con,u
  44 :            real*8 XH2MOY,XO2MOY,XN2MOY,XH2OMOY,PRESSION,TEMPENT
  45 :            real*8 XHEMOY,XCO2MOY,XCOMOY
  46 :       endsegment
  47 : C
  48 :       ibou = cpi(/1)
  49 :       sum  = 0.D0
  50 :       do 10 i=1,ibou
  51 :          sum = sum + x(i)*m(i)
  52 :  10   continue
  53 :       FLA001=P*sum/(R*T)
  54 : C
  55 :       RETURN
  56 :       END
  57 : 
  58 : 
  59 :  
  60 :  
  61 :  
  62 :  
  63 :  
  64 :  
  65 :