diff --git a/python/testtools/dorie/testtools/dorie_run/ode.py b/python/testtools/dorie/testtools/dorie_run/ode.py
index bacd2e69641b171e3bd76aed68d1b31e0ff92843..c5bc7db4662eec85cdde1e2e64fa96c001bb2527 100644
--- a/python/testtools/dorie/testtools/dorie_run/ode.py
+++ b/python/testtools/dorie/testtools/dorie_run/ode.py
@@ -108,6 +108,11 @@ def evaluate(iniinfo,runtime):
     plt.savefig("{}/flux_residual.png".format(iniinfo["output.outputPath"]))
     plt.close()
 
+    # print data into file
+    filename = "{}/residuals.txt".format(iniinfo["output.outputPath"])
+    f1,f2,f3 = res_flux.T
+    np.savetxt(filename,np.vstack((y,res_head,f1,f2,f3)).T,header="y,head,flux1,flux2,flux3")
+
     head_tol = 1E-5 if not "_ode.head_abstol" in iniinfo else float(iniinfo["_ode.head_abstol"])
     flux_tol = abs(influx) * 1e-5 if not "_ode.flux_abstol" in iniinfo else float(iniinfo["_ode.flux_abstol"])
     print("tolerances: head L2 error: {:.2e}, flux L2 error: {:.2e}".format(head_tol,flux_tol))