I'm wanting to find the area of intersection or overlap between two, 2-D fields with scattered positions. I can get this result in a roundabout way by using two instances of Map. Map will give as a bi-product a set of invalid positions which can be used to mask out areas of non-intersection on each respective field. The problem is that the process is too slow to be practical (> 24 hr run time !) because while the positions of overlap are being calculated, the data is also being interpolated, which is something I don't need.
I'm wondering if there is a way to get Map to not do the interpolation (I have tried sending in fields which contain no data, but got no improvement). Or, does anyone out there have a better method ? I would think that there would have to be others that need to do this same kind of thing. Thanks in advance for any comments or suggesions. -Sharon
