anchor.diff:                                              ((2*mnuc*AtomicNumber_oper(i))/(w**2*(2*J+1)))*p_tot(w,vesc,i)
                   end if
               end do
  +            !!$OMP end do
  +            !!$OMP end parallel
           end if
       end function OMEGA_oper
   end module opermod
  
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:     end do
      !$OMP end do
      !$OMP end parallel
      ! completes integral (2.4) of paper 1501.03729
      capped = 4.d0*pi*Rsun**3*capped
  
      if (capped .gt. 1.d100) then
  