anchor.diff:          !   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
  
positive.code:         y(1) = 0.0
          kn_start = 1.0E-10 
          kn_end = 2.0*kFn*R_a
          h = -0.1
  
          rpar(1) = temperature
          rpar(2) = ionic% Z
  
negative.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)
  
          if (idid < 0) then
              write(error_unit,*) 'n_imp: dop853 returned with error',idid
              ierr = idid
          end if
  
          deallocate(work, iwork)
      
  