compute_bfrac.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	compute_bfrac (time, xl, tburn, uparam, v1, rho1, rho10, bfrac)

Function/Subroutine Documentation

compute_bfrac()

subroutine compute_bfrac	(		time,
			xl,
			tburn,
			uparam,
			v1,
			rho1,
			rho10,
			bfrac )

Definition at line 28 of file compute_bfrac.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   I N P U T   O U T P U T   A r g u m e n t s 
C-----------------------------------------------
      my_real time,xl,tburn,v1,rho1,rho10,uparam(*),bfrac
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ibug, IBFRAC
      my_real vdet,bhe
C-----------------------------------------------
      vdet  = uparam(42)
      bhe   = uparam(44)
      ibfrac= uparam(68)
 
      !--------------------------------!
      ! Calculation of BFRAC in [0,1]  !
      !--------------------------------!
      IF(bfrac < one) THEN
        bfrac = zero
        IF(ibfrac/=1 .AND. time > -tburn)THEN
          xl = v1**third
          bfrac = vdet*(time+tburn)*two_third/xl 
        ENDIF
        IF(ibfrac/=2) bfrac  = max( bfrac , bhe * (one - rho10/rho1) )
        IF(bfrac < em03) THEN
          bfrac = zero
        ELSEIF(bfrac > one) THEN
          bfrac = one
        ENDIF
      ENDIF
      RETURN