This function takes a parameter file (pf) along with the nr, ntheta and nphi points to generate out to rmax, and it grids fields onto those points and returns a dict. center if supplied will be the center, otherwise the most dense point will be chosen. smoothed governs whether regular covering grids or smoothed covering grids will be used.