diff --git a/cgyro/src/cgyro_error_estimate.F90 b/cgyro/src/cgyro_error_estimate.F90
index 732b13602..b8769a579 100644
--- a/cgyro/src/cgyro_error_estimate.F90
+++ b/cgyro/src/cgyro_error_estimate.F90
@@ -13,7 +13,8 @@ subroutine cgyro_error_estimate
   use cgyro_globals
   use cgyro_io
   use timer_lib
-
+  use cgyro_step
+  
   implicit none
 
   integer :: i_f,itor
@@ -175,7 +176,7 @@ subroutine cgyro_error_estimate
   integration_error = integration_error/norm
 
   if (i_proc == 0) then
-     print '(7(1pe12.5))',delta_t_gk,deltah2,delta_t_last_step,delta_t_last,var_error,integration_error
+     print '(7(1pe12.5))',var_error,integration_error,sum(abs(h_x)),sum(abs(rhs(:,:,:,1)))
   endif
 
   ! Trigger code shutdown on large collisionless error