set up the vpa grid points and integration weights set up the mu grid points and integration weights
maxwell_fac = 1 unless radially global