anchor.diff:      double precision, intent(in) :: mx_in, sigma_0_SD_in,sigma_0_SI_in
       double precision :: capped_SD,capped_SI
   
  -
       call captn_general(mx_in,sigma_0_SD_in,1,0,0,capped_SD)
       call captn_general(mx_in,sigma_0_SI_in,29,0,0,capped_SI)
   
  
positive.code: !   Sticking with notation of 1504.04378. Cite that paper. Or 1605.06502 it's even better.
  !   Reference q0 is 40 MeV, and v0 is 220 km/s.
  
      module capmod
      implicit none
      double precision, parameter :: pi=3.141592653, NAvo=6.0221409d23,GMoverR = 1.908e15
  
negative.code: 
      subroutine captn_init(solarmodel)
      use capmod
      use iso_c_binding, only: c_ptr
      character (len=300) solarmodel
  !    common solarmodel
  !    external solarmodel
      print*,"Capgen initializing"
      print*,solarmodel
      if  (.not. allocated(tab_r)) then !
          call get_solar_params(solarmodel,nlines)
      else
  