anchor.diff:  module spergelpressmod
   use capmod
   implicit none
  -   
  +
   double precision, parameter :: kB=1.38064852d-16, mnucg=1.6726219e-24
   
   contains
  
positive.code: r = tab_r*Rsun ! cm
  phi = -tab_vesc**2/2.d0 ! erg/g
  mxg = mdm*1.782662d-24
  
  ! WIMP number density in isothermal approximation
  nx_isothermal = exp(-mxg*phi/kB/T_x) 
  n_0 = Nwimps/trapz(r, 4.d0*pi*r**2.d0*nx_isothermal, nlines) ! Normalize so that integral(nx) = Nwimps
  nx_isothermal = n_0*nx_isothermal
  
  
negative.code: 
  function Tx_integral(T_x, sigma_N, Nwimps, niso)
  implicit none
  ! Calculates the Tx defining integral
  
  integer, intent(in) :: niso
  double precision, intent(in) :: T_x, Nwimps
  