/[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.1 by mlosch, Wed May 9 09:46:25 2018 UTC revision 1.3 by mlosch, Wed May 9 11:43:12 2018 UTC
# Line 10  from getopt import gnu_getopt as getopt Line 10  from getopt import gnu_getopt as getopt
10  # parse command-line arguments  # parse command-line arguments
11  try:  try:
12      optlist,args = getopt(sys.argv[1:], ':', ['verbose'])      optlist,args = getopt(sys.argv[1:], ':', ['verbose'])
     print(len(args))  
13      assert len(args) == 1      assert len(args) == 1
14  except (AssertionError):  except (AssertionError):
15      sys.exit(__doc__)      sys.exit(__doc__)
# Line 40  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.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22