************************************************************************
*
*     F3total.f:
*
*     This function returns the value of the inclusive structure function
*     F3.
*
************************************************************************
      function F3total(x)
*
      implicit none
**
*     Input Variables
*
      double precision x
**
*     Internal Variables
*
      double precision F3light,F3charm,F3bottom,F3top
**
*     Output Variables
*
      double precision F3total
*
      F3total = F3light(x)
     1        + F3charm(x)
     2        + F3bottom(x)
     3        + F3top(x)
*
      return
      end
c$$$*
c$$$************************************************************************
c$$$      function F3total(x)
c$$$*
c$$$      implicit none
c$$$*
c$$$      include "../commons/grid.h"
c$$$      include "../commons/StructureFunctions.h"
c$$$      include "../commons/TMC.h"
c$$$      include "../commons/TimeLike.h"
c$$$**
c$$$*     Input Variables
c$$$*
c$$$      double precision x
c$$$**
c$$$*     Internal Variables
c$$$*
c$$$      integer n
c$$$      integer alpha
c$$$      double precision w_int_gen
c$$$      double precision tau,xi
c$$$      double precision c1,c2
c$$$      double precision tol
c$$$      parameter(tol=1d-10)
c$$$**
c$$$*     Output Variables
c$$$*
c$$$      double precision F3total
c$$$*
c$$$*
c$$$      if(TMC)then
c$$$         tau = 1d0 + 4d0 * rhop * x**2
c$$$         xi  = 2d0 * x / ( 1d0 + dsqrt(tau) )
c$$$*
c$$$         c1 = x**2 / xi**2 / tau
c$$$         c2 = 4d0 * rhop * x**3 / tau**1.5d0
c$$$*
c$$$         if(xi.lt.xmin(1)-tol.or.xi.gt.xmax+tol)then
c$$$            write(6,*) "In F3total.f:"
c$$$            write(6,*) "Invalid value of x =",xi
c$$$            call exit(-10)
c$$$         endif
c$$$         if (xi.lt.xmin(1)) xi = xmin(1)
c$$$         if (xi.gt.xmax) xi = 1d0
c$$$*
c$$$*     Interpolation
c$$$*
c$$$         F3total = 0d0
c$$$         n = inter_degree(0)
c$$$         do alpha=0,nin(0)
c$$$            F3total = F3total + w_int_gen(n,alpha,xi) 
c$$$     1           * ( c1 * F3(7,0,alpha) + c2 * I3(7,0,alpha) )
c$$$         enddo
c$$$         if(dabs(F3total).le.1d-14) F3total = 0d0
c$$$      else
c$$$         if(x.lt.xmin(1)-tol.or.x.gt.xmax+tol)then
c$$$            write(6,*) "In F3total.f:"
c$$$            write(6,*) "Invalid value of x =",x
c$$$            call exit(-10)
c$$$         endif
c$$$         if (x.lt.xmin(1)) x = xmin(1)
c$$$         if (x.gt.xmax) x = 1d0
c$$$*
c$$$*     Interpolation
c$$$*
c$$$         F3total = 0d0
c$$$         n = inter_degree(0)
c$$$         do alpha=0,nin(0)
c$$$            F3total = F3total + w_int_gen(n,alpha,x) * F3(7,0,alpha)
c$$$         enddo
c$$$         if(dabs(F3total).le.1d-14) F3total = 0d0
c$$$      endif
c$$$*
c$$$      if(Timelike) F3total = F3total / x
c$$$*
c$$$      return
c$$$      end