anchor.diff:          y(1) = 0.0
           kn_start = 1.0E-10 
           kn_end = 2.0*kFn*R_a
  -        h = -0.1
  +        h = 0.0
   
           rpar(1) = temperature
           rpar(2) = ionic% Z
  
positive.code:         &   int_default_max_step_size,int_default_max_steps, &
          & rtol,atol,itol, null_solout, iout, work, lwork, iwork, liwork,  &
          &   num_rpar, rpar, num_ipar, ipar, lout, idid)
  
          deallocate(work, iwork)
      
  
negative.code:         !   end function eone
          ! end interface
          real(dp), parameter :: um1 = 2.8, um2 = 12.973, onesixth = 1.0/6.0, fourthird = 4.0*onethird
          real(dp), parameter :: sw_max = log(huge(1.0_dp))
          real(dp) :: electroncompton
          real(dp) :: aB
          real(dp) :: eifac
  