fluxlvl.f

Go to the documentation of this file.

      subroutine fluxlvl(buft,sumg,iflag)
c
c     include the common statements.....................................
      implicit real*8 (a-h,o-z)
      include 'common_all.cmn'
      real*8 buft(4,50,46)
c
c***********************************************************************
      if(iflag.eq.1)then
         m=jjj
         n=nmum1
      else
         m=1
         n=jjjj
      endif
c
      anophi=dfloat(nophi)                                                      
      sumg = 0.0d0
      do 5001 j = m,n                                                   
         psumg = 0.0d0
         do 4003 k = 1,2                                                   
            psumg=buft(k,j,1)+buft(k,j,jpart)+psumg                               
            do 4003 i=2,noph2                                                 
            psumg=psumg+2.0d0*buft(k,j,i)
 4003    continue
         sumg=sumg+psumg*dmus2(j)*6.2831854d0/anophi                          
 5001 continue
c
      sumg=dabs(sumg)                                                      
      return                                                            
      end                                                               
c***********************************************************************