anchor.diff:          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)
  +        call transgen(sigma_0,nwimpsin,num_isotopes,nq(j),nv(j),spin_dependency,Etrans,Etranstot)
           ! print*, "Etranstot: ", Etranstot !FIXME units?
   
           ! print*,"Calling fastevap."
  
positive.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
  
      ! Initialise capture calculations
      call captn_init(modfile,0.4d0,235.d0,235.d0,550.d0)
  
negative.code:     double precision :: maxcap, nwimpsin, evapRate
      double precision, allocatable :: Etrans(:)
      double precision :: EtransTot
      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]
  