diff --git a/pybamm/discretisations/discretisation.py b/pybamm/discretisations/discretisation.py
index 847f5b662f..bcfe9949da 100644
--- a/pybamm/discretisations/discretisation.py
+++ b/pybamm/discretisations/discretisation.py
@@ -1030,6 +1030,10 @@ def check_initial_conditions_rhs(self, model):
         y0 = model.concatenated_initial_conditions
         # Individual
         for var in model.rhs.keys():
+            print('rhs')
+            print(model.rhs[var])
+            print('init')
+            print(model.initial_conditions[var])
             assert (
                 model.rhs[var].shape == model.initial_conditions[var].shape
             ), pybamm.ModelError(
diff --git a/pybamm/expression_tree/symbol.py b/pybamm/expression_tree/symbol.py
index 912b102b6d..06b9693557 100644
--- a/pybamm/expression_tree/symbol.py
+++ b/pybamm/expression_tree/symbol.py
@@ -713,6 +713,7 @@ def shape(self):
         try:
             y = np.linspace(0.1, 0.9, int(1e4))
             evaluated_self = self.evaluate(0, y, y, u="shape test")
+            print('evaluated self is ',evaluated_self)
         # If that fails, fall back to calculating how big y should really be
         except ValueError:
             state_vectors_in_node = [
@@ -726,7 +727,7 @@ def shape(self):
                 )
                 # Pick a y that won't cause RuntimeWarnings
                 y = np.linspace(0.1, 0.9, min_y_size)
-            evaluated_self = self.evaluate(0, y)
+            evaluated_self = self.evaluate(0, y, y)
 
         # Return shape of evaluated object
         if isinstance(evaluated_self, numbers.Number):
diff --git a/pybamm/solvers/base_solver.py b/pybamm/solvers/base_solver.py
index 33f94d65d0..ce8c115ea6 100644
--- a/pybamm/solvers/base_solver.py
+++ b/pybamm/solvers/base_solver.py
@@ -194,7 +194,6 @@ def report(string):
             else:
                 # Process with CasADi
                 report(f"Converting {name} to CasADi")
-                print(u_casadi)
                 func = func.to_casadi(t_casadi, y_casadi, u=u_casadi)
                 if use_jacobian:
                     report(f"Calculating jacobian for {name} using CasADi")
diff --git a/tests/unit/test_parameters/test_parameters_cli.py b/tests/unit/test_parameters/test_parameters_cli.py
index 469a100bd0..85584b7c84 100644
--- a/tests/unit/test_parameters/test_parameters_cli.py
+++ b/tests/unit/test_parameters/test_parameters_cli.py
@@ -119,3 +119,11 @@ def test_list_params(self):
         # ./input/parameters/lithium-ion/cathodes/tmp_dir
         # but must not intefere with existing input dir if it exists
         # in the current dir...
+
+if __name__ == "__main__":
+    print("Add -v for more debug output")
+    import sys
+
+    if "-v" in sys.argv:
+        debug = True
+    unittest.main()
diff --git a/tests/unit/test_solvers/test_casadi_solver.py b/tests/unit/test_solvers/test_casadi_solver.py
index e63e88ce59..71632874f2 100644
--- a/tests/unit/test_solvers/test_casadi_solver.py
+++ b/tests/unit/test_solvers/test_casadi_solver.py
@@ -337,5 +337,4 @@ def test_model_solver_with_dvdt(self):
     if "-v" in sys.argv:
         debug = True
     pybamm.settings.debug_mode = True
-    pybamm.set_logging_level("DEBUG")
     unittest.main()