using nanmin/nanmax instead of minmax + test

Author: delandmeterp

parcels/grid.py

@@ -52,7 +52,7 @@ def __init__(self, lon, lat, time, time_origin, mesh):
         self.meridional_halo = 0
         self.lat_flipped = False
         self.defer_load = False
-        self.lonlat_minmax = np.array([np.min(lon), np.max(lon), np.min(lat), np.max(lat)], dtype=np.float32)
+        self.lonlat_minmax = np.array([np.nanmin(lon), np.nanmax(lon), np.nanmin(lat), np.nanmax(lat)], dtype=np.float32)
 
     @property
     def ctypes_struct(self):
@@ -203,7 +203,7 @@ def add_periodic_halo(self, zonal, meridional, halosize=5):
                                       self.lat, self.lat[0:halosize] + latshift))
             self.ydim = self.lat.size
             self.meridional_halo = halosize
-        self.lonlat_minmax = np.array([np.min(self.lon), np.max(self.lon), np.min(self.lat), np.max(self.lat)], dtype=np.float32)
+        self.lonlat_minmax = np.array([np.nanmin(self.lon), np.nanmax(self.lon), np.nanmin(self.lat), np.nanmax(self.lat)], dtype=np.float32)
 
 
 class RectilinearZGrid(RectilinearGrid):

tests/test_grids.py

@@ -327,7 +327,7 @@ def sampleSpeed(particle, fieldset, time, dt):
     class MyParticle(ptype[mode]):
         speed = Variable('speed', dtype=np.float32, initial=0.)
 
-    pset = ParticleSet.from_list(field_set, MyParticle, lon=[400], lat=[600])
+    pset = ParticleSet.from_list(field_set, MyParticle, lon=[400, -200], lat=[600, 600])
     pset.execute(pset.Kernel(sampleSpeed), runtime=0, dt=0)
     assert(np.allclose(pset[0].speed, 1000))