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?

thanks,
-steve


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?