diff --git a/pydda/retrieval/wind_retrieve.py b/pydda/retrieval/wind_retrieve.py
index 1efac38e..8d3973f6 100644
--- a/pydda/retrieval/wind_retrieve.py
+++ b/pydda/retrieval/wind_retrieve.py
@@ -643,8 +643,6 @@ def loss_and_gradient(x):
             )
             winds = {"winds": winds}
             winds, state = solver.run(winds, bounds=bounds)
-            print("Cost function: %4.3f" % winds.theta)
-            print("L^2 norm of gradient: %4.3f" % jnp.linalg.norm(state.grad["winds"]))
             winds = [np.asanyarray(winds["winds"])]
 
         winds = np.reshape(