Index: wflow-py/wflow/wflow_lib.py
===================================================================
diff -u -rc3ba6c0bd01725069cb0e3f596fd5475847591c1 -r6a7bc8aa1448ca6ea1197143f25524e413ac27b5
--- wflow-py/wflow/wflow_lib.py	(.../wflow_lib.py)	(revision c3ba6c0bd01725069cb0e3f596fd5475847591c1)
+++ wflow-py/wflow/wflow_lib.py	(.../wflow_lib.py)	(revision 6a7bc8aa1448ca6ea1197143f25524e413ac27b5)
@@ -51,6 +51,8 @@
 
 
 
+
+
 Verbose=0
 
 def lddcreate_save(lddname, dem, force, corevolume=1E35, catchmentprecipitation=1E35, corearea=1E35, outflowdepth=1E35):
@@ -586,8 +588,8 @@
     try:
         c = xcor.ndim
     except:
-        xcor = array([xcor])
-        ycor = array([ycor])
+        xcor = numpy.array([xcor])
+        ycor = numpy.array([ycor])
 
     # Loop over points and "burn in" map
     for n in range(0,xcor.size):
@@ -821,8 +823,8 @@
     resY    = geotrans[5]
     cols = ds.RasterXSize
     rows = ds.RasterYSize
-    x = linspace(originX+resX/2,originX+resX/2+resX*(cols-1),cols)
-    y = linspace(originY+resY/2,originY+resY/2+resY*(rows-1),rows)
+    x = numpy.linspace(originX+resX/2,originX+resX/2+resX*(cols-1),cols)
+    y = numpy.linspace(originY+resY/2,originY+resY/2+resY*(rows-1),rows)
     # Retrieve raster
     RasterBand = ds.GetRasterBand(1) # there's only 1 band, starting from 1
     data = RasterBand.ReadAsArray(0,0,cols,rows)