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