anchor.diff:              end if
               tally = tally + W_array(Wtype,iso,tau,tauprime,k) * yConverse**(k-1) * G
           end do
  +        !$OMP end do
  +        !$OMP end parallel
           sumW = tally
       end function sumW
       
  
positive.code:     ! completes integral (2.3) in paper 1501.03729 (gives dC/dV as fn of radius)
      !$OMP parallel default(none) shared(nlines, vesc_halo, epsabs, epsrel, limit, u_int_res, capped, tab_r, tab_dr) &
      !$OMP private(ri, result, abserr, neval, ier, alist, blist, rlist, elist, iord, last)
      ! threadnums = omp_get_num_threads()
      ! threadid = omp_get_thread_num()
      ! print *, "There are", threadnums, "lights!"
      !$OMP do
      do ri=1,nlines !loop over the star
          result = 0.d0
  
negative.code:                     p_tot = p_tot + 1./mnuc**2 * RS1D(tau,taup,c,j_chi,coupling_Array) * sumW(w,vesc,i,tau,taup,8,1)
                  end if
              end do
              !$OMP end do
              !$OMP end parallel
          end do
          !$OMP end do
          !$OMP end parallel
          p_tot = p_tot  * hbar**2 * c0**2
      end function p_tot
  
  