The code has now a more stable way of calculating the modified Bessel functions that also work for large kperp. We have absorbed the exp(-z) term directly in the functions BESSI0 and BESSI1.
The code has now a more stable way of calculating the modified Bessel functions that also work for large kperp. We have absorbed the exp(-z) term directly in the functions BESSI0 and BESSI1.