C BMUG      SOURCE    CHAT      05/01/12    21:41:35     5004           
      FUNCTION BMUG(T,PHI,T0)
      IMPLICIT INTEGER(I-N)
      IMPLICIT REAL*8 (A-H,O-Z)
C
C     operateur FUITE
C>>>  calcule la viscosite du melange air vapeur
C       T  : temperature du fluide en ›K
C      PHI : rapport pression partielle vapeur/pression totale
C       T0 : temperature de reference
C
      BMUG=(1.72+4.4E-3*(T-T0))*(1-PHI)+(.845+4.07E-3*(T-T0-10))*PHI
      BMUG=BMUG*1E-5
      RETURN
      END