gvmu is initialised with a Maxwellian weighting for flux tube simulations, with the Maxwellian evaluated at ia = 1 we are undoing that weighting here, so also need to use ia = 1