anchor.diff:      double precision :: maxcap, nwimpsin, evapRate
       double precision, allocatable :: Etrans(:)
       double precision :: EtransTot
  -    integer :: nq(7), nv(7), i, j, nlines
  +    integer :: nq(7), nv(7), i, j, nlines, num_isotopes, spin_dependency
   
       ! Choose velocity and momentum transfer powers in differential cross-section
       nq = [0,-1,1,2,0,0,0]
  
positive.code:         nwimpsin = 1.2d42
          ! nwimpsin = capped_sd*3.d7*4.57d9
          ! print*,"Calling transgen, with nwimpsin = ", nwimpsin
          call transgen(sigma_0,nwimpsin,1,nq(j),nv(j),1,Etrans,Etranstot)
          ! print*, "Etranstot: ", Etranstot !FIXME units?
  
          ! print*,"Calling fastevap."
  
negative.code:     !modfile = "solarmodels/model_gs98_nohead.dat"
      !modfile = "solarmodels/struct_b16_agss09_nohead.dat"
      modfile = "solarmodels/struct_b16_agss09_nohead_Tsmoothed.dat" !temperature smoothed to not nonsense
  
      ! number of isotopes capt'n will loop over in the calculation: up to 29 isotopes
      num_isotopes = 29
  
      ! zero for spin_independent, one for spin_dependent
      spin_dependency = 0
  
      ! Initialise capture calculations
      call captn_init(modfile,0.4d0,235.d0,235.d0,550.d0)
  