So.. I've discovered that I can do something like<div><br></div><div>z = array([x,y])-m.coors</div><div>zz = (z*z).sum(axis=1)<br></div><div>np.where(zz&lt;0.1)<br></div><div><br></div><div><div>to get the index of the point closest to the field point (with sorting or something if I get more than one.)</div><div><br></div><div>Next I can figure out from m.conns which nodes are connected to that one, and go on from there?</div><div><br></div><div>is this a reasonable path?</div><div><br></div><div>thanks,</div><div>-steve</div><div><br></div><br><blockquote class="gmail_quote" style="margin: 0;margin-left: 0.8ex;border-left: 1px #ccc solid;padding-left: 1ex;"><div>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?</div><div><br></div></blockquote></div>