diff --git a/gusto/time_discretisation.py b/gusto/time_discretisation.py
index a1221c6b9..03a42fe81 100644
--- a/gusto/time_discretisation.py
+++ b/gusto/time_discretisation.py
@@ -959,9 +959,6 @@ def solve_stage(self, x0, stage):
         if self.increment_form:
             self.x1.assign(x0)
 
-            # Use x0 as a first guess (otherwise may not converge)
-            self.x_out.assign(x0)
-
             for i in range(stage):
                 self.x1.assign(self.x1 + self.dt*self.butcher_matrix[stage-1, i]*self.k[i])
             for evaluate in self.evaluate_source: