Add utility function to calculate bin edges needed for pcolormesh
Update plotnll to use getBinEdges