bugfix --- actually use the nmax and mmax values

i.e. don't kill my computer trying to compute all the many neighbour
channel flow polynomials