/[MITgcm]/MITgcm_contrib/mlosch/optim_m1qn3/plotfc.py
ViewVC logotype

Diff of /MITgcm_contrib/mlosch/optim_m1qn3/plotfc.py

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.2 by mlosch, Wed May 9 11:36:59 2018 UTC revision 1.3 by mlosch, Wed May 9 11:43:12 2018 UTC
# Line 39  def get_output (fname, mystring): Line 39  def get_output (fname, mystring):
39  iters, simuls, fc = get_output(fname, "f=")  iters, simuls, fc = get_output(fname, "f=")
40  # sort out restarts  # sort out restarts
41  iters0 = np.asarray(iters)  iters0 = np.asarray(iters)
42  for k,it in enumerate(iters[0:]):  for k,it in enumerate(iters[1:]):
43      if iters0[k]<iters0[k-1]:      print(k,it)
44          iters0[k:] = iters0[k:]+(iters0[k-1]-iters0[k]+1)      kp1=k+1
45        if iters0[kp1]<iters0[kp1-1]:
46            iters0[kp1:] = iters0[kp1:]+(iters0[k]-iters0[kp1]+1)
47                    
48  fig, ax1 = plt.subplots()  fig, ax1 = plt.subplots()
49  ax1.semilogy(iters0,(np.asarray(fc)-fc[-1]),'bx-')  ax1.semilogy(iters0,(np.asarray(fc)-fc[-1]),'bx-')
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22