My problem is that I have a bunch of regions, in this case they are tectonic plates. Each region has a complex boundary, described by polygons, in a text file with co-ordinates. I have a long list of points (>100,000) distributed globally and I need to know in which region each point lies. Basically I need an algorithm, that can be easily coded up in fortran. Any helpers?