diff --git a/CHANGELOG.rst b/CHANGELOG.rst
index fa71e7e..e84a3eb 100644
--- a/CHANGELOG.rst
+++ b/CHANGELOG.rst
@@ -5,7 +5,7 @@ Change Log
 in development
 **************
 
-0.6.0  (19-12-2024)
+0.5.9  (19-12-2024)
 ********************
 * Fixed numpy 2.0 compatibility issues
 * Fixed bug in strahler stream order method when more than two streams meet.
diff --git a/pyflwdir/__init__.py b/pyflwdir/__init__.py
index 9a199b6..324059c 100644
--- a/pyflwdir/__init__.py
+++ b/pyflwdir/__init__.py
@@ -1,11 +1,26 @@
 """Fast methods to work with hydro- and topography data in pure Python."""
 
 # version number without 'v' at start
-__version__ = "0.5.8"
+__version__ = "0.5.9"
 
 # submodules
-from .flwdir import *
-from .pyflwdir import *
-from .dem import *
-from .core_nextxy import *
 from . import gis_utils, regions
+
+# public functions
+from .core_nextxy import read_nextxy
+from .dem import fill_depressions, slope
+from .flwdir import Flwdir, from_dataframe
+from .pyflwdir import FlwdirRaster, from_array, from_dem
+
+__all__ = [
+    "Flwdir",
+    "FlwdirRaster",
+    "from_array",
+    "from_dataframe",
+    "from_dem",
+    "read_nextxy",
+    "gis_utils",
+    "regions",
+    "slope",
+    "fill_depressions",
+]
diff --git a/pyflwdir/dem.py b/pyflwdir/dem.py
index 53a7448..3de8455 100644
--- a/pyflwdir/dem.py
+++ b/pyflwdir/dem.py
@@ -230,8 +230,8 @@ def slope(elevtn, nodata=-9999.0, latlon=False, transform=gis_utils.IDENTITY):
     """Returns the local gradient
 
     The slope is calculated on the basis of the dem in a 3 x 3 cell window, using 2nd order partial derivatives.
-    The slope [m/m] is given as the increase in height per distance in horizontal direction.    
-    
+    The slope [m/m] is given as the increase in height per distance in horizontal direction.
+
     Parameters
     ----------
     elevnt : 1D array of float