Index: wflow-py/wflow/wflow_lib.py
===================================================================
diff -u -r3e39e84af48f1bcb5ec0d243748147be223674f2 -r9ce029ed9592f4d55fd733ffef41ac5a6ace1583
--- wflow-py/wflow/wflow_lib.py	(.../wflow_lib.py)	(revision 3e39e84af48f1bcb5ec0d243748147be223674f2)
+++ wflow-py/wflow/wflow_lib.py	(.../wflow_lib.py)	(revision 9ce029ed9592f4d55fd733ffef41ac5a6ace1583)
@@ -42,6 +42,7 @@
 import configparser
 import osgeo.gdal as gdal
 from osgeo.gdalconst import *
+import pcraster
 from pcraster import *
 from pcraster.framework import *
 import scipy
@@ -986,7 +987,7 @@
 
     x = pcr2numpy(xcoordinate(defined(in_map)), numpy.nan)
     y = pcr2numpy(ycoordinate(defined(in_map)), numpy.nan)
-    XX = pcr2numpy(celllength(), 0.0)
+    cell_length = float(celllength())
 
     # simple check to use both floats and numpy arrays
     try:
@@ -1001,8 +1002,8 @@
             print(n)
         diffx = x - xcor[n]
         diffy = y - ycor[n]
-        col_ = numpy.absolute(diffx) <= (XX[0, 0] * tolerance)  # cellsize
-        row_ = numpy.absolute(diffy) <= (XX[0, 0] * tolerance)  # cellsize
+        col_ = numpy.absolute(diffx) <= (cell_length  * tolerance)  # cellsize
+        row_ = numpy.absolute(diffy) <= (cell_length  * tolerance)  # cellsize
         point = point + numpy2pcr(Scalar, ((col_ * row_) * (n + 1)), numpy.nan)
 
     return ordinal(point)