Esope source files

   1 : C BCALQ0    SOURCE    DURA      06/07/11    21:16:14     5481           
   2 :       SUBROUTINE BCALQ0(PE,PVE,TE,PS,XL,E,B,RUG,Q0,RECU,XKUL,XKUT1,
   3 :      * XKUT2,XKUT3,XKUT4)
   4 :       IMPLICIT INTEGER(I-N)
   5 :       IMPLICIT REAL*8 (A-H,O-Z)
   6 : C
   7 : C     operateur FUITE
   8 : C>>>  calcule un debit d'initialisation en supposant le fluide
   9 : C     isotherme a la temperature d'injection TE
  10 : C     RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4 : coef lois de frot utilisateur
  11 : C
  11 :       CALL BPHYS(T0,P0,RA,RV,CA,CV,CL,XLAT,ROL,XKL,XKT,REL)
  12 :       Q=10.
  13 :       RO = PE/RA/TE
  14 :       PHI=0.001
  15 :       XMU= BMUG(TE,PHI,T0)
  16 :       DO 1  I=1,20
  17 :       RE = Q*2*E/XMU
  18 :       Q =BKFRO(RE,REL,XKL,XKT,2*E,RUG,RECU,XKUL,XKUT1,XKUT2,XKUT3,XKUT4)
  19 :       Q =4*E*RO*(PE-PS)/Q/XL
  20 :       Q = SQRT(Q)
  21 :   1   CONTINUE
  22 :       Q0=Q
  23 :       RETURN
  24 :       END
  25 :  
  26 :  
  27 :  
  28 :