text: # comment out N00,N01,Nt assignments in run_rom.py
  # replace with N00 += a, N01 += b, Nt += c
  
  combs = [(-1,-1,0),(-1,-1,-1),(-1,-1,-2),(-1,-1,-3),\
           (-2,-1,0),(-2,-1,-1),(-2,-1,-2),(-2,-1,-3),\
           (-2,-2,0),(-2,-2,-1),(-2,-2,-2),(-2,-2,-3),\
           (-2,-3,0),(-2,-3,-1),(-2,-3,-2),(-2,-3,-3),\
           (-3,-3,0),(-3,-3,-1),(-3,-3,-2),(-3,-3,-3),\
           (-3,-4,0),(-3,-4,-1),(-3,-4,-2),(-3,-4,-3),\
           (-4,-4,0),(-4,-4,-1),(-4,-4,-2),(-4,-4,-3),\
           (-4,-5,0),(-4,-5,-1),(-4,-5,-2),(-4,-5,-3),\
           (-5,-5,0),(-5,-5,-1),(-5,-5,-2),(-5,-5,-3),\
                     (-5,-6,-1),(-5,-6,-2),(-5,-6,-3),\
                                (-6,-6,-2),(-6,-6,-3)]
  
  for (a,b,c) in combs:
      exec(open('./run_rom.py').read())
      print('Nb = {:d}, Nt = {:d}'.format(Nb,Nt))
      exec(open('./recon.py').read())
  