diff --git a/docs/source/quickstart.rst b/docs/source/quickstart.rst
index 57c3116..1fbaaeb 100644
--- a/docs/source/quickstart.rst
+++ b/docs/source/quickstart.rst
@@ -12,7 +12,7 @@ Users can run a PyFlowline simulation in the following steps:
     * `cython` (optional, for performance)
     * `cartopy` (optional, for visualization)
     * `matplotlib` (optional, for visualization)
-    * `utm` (optional, for visualization)
+    
 2. Clone the latest PyFlowline repository from https://github.com/changliao1025/pyflowline. Or Install the PyFlowline through Conda for a released version.
 3. Download the additional large files (DEM and MPAS mesh) and move them under the `data/susquehanna/input` folder.
 4. Change the `sFilename_mesh_netcdf`, `sFilename_basins`, and `sFilename_flowline_filter` to the actual paths,
diff --git a/examples/susquehanna/run_simulation_hexagon.py b/examples/susquehanna/run_simulation_hexagon.py
index 9e0a3d3..11dc525 100644
--- a/examples/susquehanna/run_simulation_hexagon.py
+++ b/examples/susquehanna/run_simulation_hexagon.py
@@ -2,13 +2,16 @@
 from pathlib import Path
 from os.path import realpath
 
-from pyflowline.pyflowline_read_model_configuration_file import pyflowline_read_model_configuration_file
 
 
 #===================================
 #set up workspace path
 #===================================
 sPath_parent = str(Path(__file__).parents[2]) # data is located two dir's up
+import sys
+sys.path.append(sPath_parent)
+from pyflowline.pyflowline_read_model_configuration_file import pyflowline_read_model_configuration_file
+
 sPath_data = realpath( sPath_parent +  '/data/susquehanna' )
 sWorkspace_input =  str(Path(sPath_data)  /  'input')
 sWorkspace_output=  str(Path(sPath_data)  /  'output')
@@ -42,8 +45,6 @@
 oPyflowline.flowline_simplification()
 aCell = oPyflowline.mesh_generation()
 oPyflowline.reconstruct_topological_relationship(aCell)
-oPyflowline.analyze()
-oPyflowline.evaluate()
 oPyflowline.export()
 iCase_index= iCase_index+1
            
diff --git a/examples/susquehanna/run_simulation_latlon.py b/examples/susquehanna/run_simulation_latlon.py
index 94b8b8a..37917d6 100644
--- a/examples/susquehanna/run_simulation_latlon.py
+++ b/examples/susquehanna/run_simulation_latlon.py
@@ -2,13 +2,14 @@
 from pathlib import Path
 from os.path import realpath
 
-from pyflowline.pyflowline_read_model_configuration_file import pyflowline_read_model_configuration_file
-
-
 #===================================
 #set up workspace path
 #===================================
 sPath_parent = str(Path(__file__).parents[2]) # data is located two dir's up
+sPath_parent = str(Path(__file__).parents[2]) # data is located two dir's up
+import sys
+sys.path.append(sPath_parent)
+from pyflowline.pyflowline_read_model_configuration_file import pyflowline_read_model_configuration_file
 sPath_data = realpath( sPath_parent +  '/data/susquehanna' )
 sWorkspace_input =  str(Path(sPath_data)  /  'input')
 sWorkspace_output=  str(Path(sPath_data)  /  'output')
@@ -41,8 +42,6 @@
 oPyflowline.flowline_simplification()
 aCell = oPyflowline.mesh_generation()
 oPyflowline.reconstruct_topological_relationship(aCell)
-oPyflowline.analyze()
-oPyflowline.evaluate()
 oPyflowline.export()
 
 
diff --git a/examples/susquehanna/run_simulation_mpas.py b/examples/susquehanna/run_simulation_mpas.py
index 9980468..a5f5393 100644
--- a/examples/susquehanna/run_simulation_mpas.py
+++ b/examples/susquehanna/run_simulation_mpas.py
@@ -39,8 +39,6 @@
 oPyflowline.flowline_simplification()
 aCell = oPyflowline.mesh_generation()
 oPyflowline.reconstruct_topological_relationship(aCell)
-oPyflowline.analyze()
-oPyflowline.evaluate()
 oPyflowline.export()
 
 print('Finished')
diff --git a/examples/susquehanna/run_simulation_square.py b/examples/susquehanna/run_simulation_square.py
index b8cc1a4..ca55aae 100644
--- a/examples/susquehanna/run_simulation_square.py
+++ b/examples/susquehanna/run_simulation_square.py
@@ -2,13 +2,15 @@
 from pathlib import Path
 from os.path import realpath
 
-from pyflowline.pyflowline_read_model_configuration_file import pyflowline_read_model_configuration_file
-
 
 #===================================
 #set up workspace path
 #===================================
 sPath_parent = str(Path(__file__).parents[2]) # data is located two dir's up
+sPath_parent = str(Path(__file__).parents[2]) # data is located two dir's up
+import sys
+sys.path.append(sPath_parent)
+from pyflowline.pyflowline_read_model_configuration_file import pyflowline_read_model_configuration_file
 sPath_data = realpath( sPath_parent +  '/data/susquehanna' )
 sWorkspace_input =  str(Path(sPath_data)  /  'input')
 sWorkspace_output=  str(Path(sPath_data)  /  'output')
@@ -41,8 +43,6 @@
 oPyflowline.flowline_simplification()
 aCell = oPyflowline.mesh_generation()
 oPyflowline.reconstruct_topological_relationship(aCell)
-oPyflowline.analyze()
-oPyflowline.evaluate()
 oPyflowline.export()