Bug fix in grid computation.

Author: evanmason

setup.py

@@ -1,6 +1,7 @@
 # -*- coding: utf-8 -*-
 from setuptools import setup, find_packages, Extension
 from Cython.Distutils import build_ext as cython_build_ext
+import numpy
 
 setup(
     name="pyeddytracker",
@@ -23,9 +24,10 @@
         'build_ext': cython_build_ext,
     },
     ext_modules=[Extension("py_eddy_tracker.tools",
-                           ["src/py_eddy_tracker/tools.pyx"])],
+                           ["src/py_eddy_tracker/tools.pyx"],
+			   include_dirs=[numpy.get_include()])],
     setup_requires=[
-        'numpy>=1.9'],
+        'numpy>=1.8'],
     install_requires=[
         'numpy>=1.9',
         'matplotlib>=1.2.1',

src/py_eddy_tracker/global_tracking.py

@@ -26,12 +26,14 @@ def create_variable(self, handler_nc, kwargs_variable,
             var.setncattr(attr, attr_value)
 
         var[:] = data
+            # Variable names are "track_0000", "track_0001", etc.
 
         var.setncattr('min', var[:].min())
         var.setncattr('max', var[:].max())
 
     def write_netcdf(self):
         """Write a netcdf with eddy
+        Write eddy property data to tracks.
         """
         eddy_size = None
         for key in VAR_DESCR:

src/py_eddy_tracker/py_eddy_tracker.py

@@ -224,6 +224,7 @@ def half_interp(self, h_one, h_two):
         """
         Speed up frequent operations of type 0.5 * (arr[:-1] + arr[1:])
         """
+        h_one = np.copy(h_one)
         h_one += h_two
         h_one *= 0.5
         return h_one

src/scripts/make_eddy_track_CLS.py

@@ -24,7 +24,7 @@
 
 make_eddy_track_CLS.py
 
-Version 1.4.2
+Version 2.0.3
 
 
 ===============================================================================