allocate and initialise kperp2 and dkperp2dr allocate and initialise vperp2 init_bessel sets up arrays needed for gyro-averaging; for a flux tube simulation, this is j0 and j1; for a flux annulus simulation, gyro-averaging is non-local in ky and so more effort is required