diff --git a/src/invert4geom/inversion.py b/src/invert4geom/inversion.py
index 7d5f9fa..771074e 100644
--- a/src/invert4geom/inversion.py
+++ b/src/invert4geom/inversion.py
@@ -1101,7 +1101,7 @@ def run_inversion(
     return results
 
 
-def run_inversion_workflow(  # equivalent to monte_carlo_full_workflow
+def run_inversion_workflow(
     grav_df: pd.DataFrame,
     create_starting_topography: bool = False,
     create_starting_prisms: bool = False,
diff --git a/src/invert4geom/optimization.py b/src/invert4geom/optimization.py
index eaae032..521fc96 100644
--- a/src/invert4geom/optimization.py
+++ b/src/invert4geom/optimization.py
@@ -1649,7 +1649,6 @@ def optimize_inversion_zref_density_contrast(
             calculate_regional_misfit=True,
             zref=best_zref,
             density_contrast=best_density_contrast,
-            plot_convergence=True,
             fname=fname,
             **new_kwargs,
         )
@@ -1787,10 +1786,6 @@ def optimize_inversion_zref_density_contrast_kfolds(
     test_dfs, train_dfs = cross_validation.kfold_df_to_lists(testing_training_df)
     log.info("Constraints split into %s folds", len(test_dfs))
 
-    # for i in range(len(test_dfs)):
-    # log.info("points in fold %s train set:", len(train_dfs[i]))
-    # log.info("point in fold %s test set:", len(test_dfs[i]))
-
     regional_grav_kwargs = kwargs.pop("regional_grav_kwargs", None)
 
     starting_topography_kwargs = kwargs.pop("starting_topography_kwargs", None)