So.. I've discovered that I can do something like
z = array([x,y])-m.coors zz = (z*z).sum(axis=1) np.where(zz<0.1)
to get the index of the point closest to the field point (with sorting or something if I get more than one.)
Next I can figure out from m.conns which nodes are connected to that one, and go on from there?
is this a reasonable path?
Is it best to re-map the data onto a rectangular grid, or is there some way
to quickly find the nodes nearest some particular spatial coordinates and interpolate directly from the irregular grid data?