anchor.diff:        end if
         allocate(u_int_res(nlines))
   
  +      maxcapped = maxcap(mx_in)
         capped = 0.d0
         !Loop over the different elements
  -    !$OMP parallel do
         do eli = 1, niso
           capiso(eli) = 0.d0
   
  
positive.code:           call dsntdqagse(integrand,vdist_over_u,umin,umax, &
            epsabs,epsrel,limit,u_int_res(ri),abserr,neval,ier,alist,blist,rlist,elist,iord,last)
            u_int_res(ri) = u_int_res(ri) * 2.d0 * sigma_N * NAvo * tab_starrho(ri)*tab_mfr(ri,eli) * (muplus/mx_in)**2
            capped = capped + tab_r(ri)**2*u_int_res(ri)*tab_dr(ri)
            capiso(eli) = capiso(eli) + tab_r(ri)**2*u_int_res(ri)*tab_dr(ri)
  
            if (isnan(capped)) then
              capped = 0.d0
  
negative.code:         end do
  
        end do
  
        capped = 4.d0*pi*Rsun**3*capped
        capiso(:) = 4.d0*pi*Rsun**3*capiso(:)
  