Exact Gauss-Legendre weights are used for integration over the interval , where ENERGY_MAX.
Beyond , the distribution is assumed to be constant (since GYRO solves for the deviation from a Maxwellian, the full perturbed distribution is therefore assumed to be Maxwellian beyond ).
Very accurate results are obtained with ENERGY_GRID=8 and ENERGY_MAX=5.0.