--- a/globalneighbors/grid.py	Sun Jun 25 11:37:52 2017 -0700
+++ b/globalneighbors/grid.py	Sun Jun 25 12:28:36 2017 -0700
@@ -59,4 +59,12 @@
 
     def neighbor_points(self, i, j):
         """return all points in a lat-lon box"""
-        raise NotImplementedError('TODO') # -> record TODO items
+        geoids = set()
+        geoids.update(self[(i,j)])
+        for ii, jj in self.neighbors(i,j):
+            geoids.update(self[ii,jj])
+        return geoids
+
+class GriddedLocations(object):
+
+    def __init__(self,