subroutine means(wpf,opf,dopf,bpf,dbpf,nopf,ntau,oppsum,bpmean
     &   ,dbpmean)

        implicit double precision (a-h,o-z) 
        include "parameters.dat"

        dimension wpf(nd2), opf(nd1,nd2), dopf(nd1,nd2), bpf(nd1,nd2)
     &   , dbpf(nd1,nd2), bpmean(nd1), dbpmean(nd1), oppsum(nd1)

        do n = 1, ntau
           sumo = zero
           do k = 1, nopf
              sumo = sumo + wpf(k)*opf(n,k)
           end do
           oppsum(n) = sumo
        end do

        do n = 1, ntau
           sumb = zero
           do k = 1, nopf
              sumb = sumb + wpf(k)*opf(n,k)*bpf(n,k)
           end do
           bpmean(n) = sumb
        end do

        do n = 1, ntau
           sumd = zero
           do k = 1, nopf
              sumd = sumd + wpf(k)*(dopf(n,k)*bpf(n,k)
     &         +opf(n,k)*dbpf(n,k))
           end do
           dbpmean(n) = sumd
        end do

        do n = 1, ntau
           bpmean(n) = bpmean(n)/oppsum(n)
        end do

        do n = 1, ntau
           dbpmean(n) = dbpmean(n)/oppsum(n)
        end do

        return
        end