anchor.diff:  r = tab_r*Rsun ! cm
   phi = -tab_vesc**2/2.d0 ! erg/g
   mxg = mdm*1.782662d-24
  -
  +! open(95,"randphi.dat")
  +! write(95,*) r , phi
  +! close(95)
   ! WIMP number density in isothermal approximation
  -nx_isothermal = exp(-mxg*phi/kB/T_x) 
  +nx_isothermal = exp(-mxg*(phi-phi(1))/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
   
  
positive.code: !	-nx_isothermal: Calculates the WIMP density in the Spergel-Press scheme
  ! 	-Etrans_sp: calculates the WIMP transported energy (eps_x) given the WIMP temperature (Tx)
  !	-Tx_integral: to be used in newtons_meth
  !	-newtons_meth: solves Tx_integral=0 which defines Tx 
  
  ! All units are cgs except tab_r and tab_dr
  ! I apologize for the long function calls.
  
negative.code: module spergelpressmod
  use capmod
  implicit none
  
  double precision, parameter :: kB=1.38064852d-16, mnucg=1.6726219e-24
  
  contains
  