--- MITgcm_contrib/mlosch/optim_m1qn3/plotfc.py	2018/05/09 11:36:59	1.2
+++ MITgcm_contrib/mlosch/optim_m1qn3/plotfc.py	2018/05/09 11:43:12	1.3
@@ -39,9 +39,11 @@
 iters, simuls, fc = get_output(fname, "f=")
 # sort out restarts
 iters0 = np.asarray(iters)
-for k,it in enumerate(iters[0:]):
-    if iters0[k]<iters0[k-1]:
-        iters0[k:] = iters0[k:]+(iters0[k-1]-iters0[k]+1)
+for k,it in enumerate(iters[1:]):
+    print(k,it)
+    kp1=k+1
+    if iters0[kp1]<iters0[kp1-1]:
+        iters0[kp1:] = iters0[kp1:]+(iters0[k]-iters0[kp1]+1)
         
 fig, ax1 = plt.subplots()
 ax1.semilogy(iters0,(np.asarray(fc)-fc[-1]),'bx-')