anchor.diff:  !   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.
   
  +!    character*300 solarmodel
  +!    common solarmodel
  +
       module capmod
       implicit none
       double precision, parameter :: pi=3.141592653, NAvo=6.0221409d23,GMoverR = 1.908e15
  
positive.code:     !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
      !read in solar parameters from Aldo Serenelli-style files, with header removed
      subroutine get_solar_params(filename,nlines)
      character*50 :: filename
      double precision :: Temp, Pres, Lumi !these aren't used, but dummies are required
      double precision, allocatable :: phi(:) !this is used briefly
      integer :: i,j, nlines,iostatus
  
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
  