anchor.diff:          print*,"infinite amount of dark matter in the Sun. Best to look into that."
         end if
   
  -      maxcapped = maxcap(mx_in)
  -      if (capped .gt. maxcapped) then
  -        capped = maxcapped
  -      end if
       end subroutine captn_general
   
   
  
positive.code: 
      subroutine captn_general(mx_in,sigma_0,niso,nq_in,nv_in,spin_in,capped,capiso)
        use capmod
        use omp_lib
        implicit none
        integer, intent(in):: nq_in, nv_in, niso, spin_in
        ! integer, intent(in):: spin_in
  
negative.code:         end do
  
        end do
  
        capped = 4.d0*pi*Rsun**3*capped
        capiso(:) = 4.d0*pi*Rsun**3*capiso(:)
  