header should print as intended now

main.py

@@ -235,6 +235,10 @@ def main(n_holes, n_particles, ref=[], d=1.0, g=0.5, pb=0.0, verbose=1, flow_dat
 
     iters = 0
     convergence = 0
+
+
+    print("iter,      \t    s, \t         E, \t         ||eta1B||, \t         ||eta2B||")
+        
     while solver.successful() and solver.t < sfinal:
 
         #print('solver success,', solver.successful())
@@ -244,8 +248,6 @@ def main(n_holes, n_particles, ref=[], d=1.0, g=0.5, pb=0.0, verbose=1, flow_dat
         s_vals.append(solver.t)
         E_vals.append(Es)
 
-        print("iter,\t s, \t E, \t ||eta1B||, \t ||eta2B||")
-
         if iters %10 == 0 and verbose: 
             norm_eta1B = np.linalg.norm(np.ravel(wg.eta1B))
             norm_eta2B = np.linalg.norm(np.ravel(wg.eta2B))