anchor.diff:          w = 4.0*um2*kF2/qD2*w1
       
           sw = s*w
  -        if (w < 1.0e4) then
  +        if (sw < sw_max) then
               fac = exp(sw)*(eone(sw)-eone(sw+w))
               L1 = 0.5*(log(1.0+1.0/s) + (1.0-exp(-w))*s/(s+1.0) - fac*(1.0+sw))
               L2 = 0.5*(1.0-(1.0-exp(-w))/w + s*(s/(1.0+s)*(1.0-exp(-w))  &
  
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
  