def METHOD_NAME(args):
    """
    Run the test.
    """
    global nTimeSteps
    # Tolerance to require for agreement with CODE
    TOLERANCE = 2e-2    # 2%
    success = True
    workdir = pathlib.Path(__file__).parent.absolute()
    B, CODEbump = loadCODE('{}/CODE-bumps.mat'.format(workdir))
    B = B[:,0]
    CODEbump = CODEbump[:,0]
    nt     = nTimeSteps
    nB     = B.size
    bumpP  = np.zeros((nB,))
    for i in range(1, nB):
    #for i in [2]:
        print('Checking B = {} T... '.format(B[i]), end="")
        try:
            bumpP[i] = runB(B[i],CODEbump[i])
        except Exception as e:
            print(e)
            bumpP[i] = 0
            return False
        # Compare runaway rates
        Delta = np.abs(bumpP[i] / CODEbump[i] - 1.0)
        print("Delta = {:.3f}%".format(Delta*100))
        if Delta > TOLERANCE:
            dreamtests.print_error("DREAM synchrotron bump location deviates from CODE at B = {} T".format(B[i]))
            dreamtests.print_error("DREAM bump at p = {} mc, CODE bump at p = {} mc".format(bumpP[i], CODEbump[i]))
            success = False
    # Save
    if args['save']:
        with h5py.File('{}/DREAM-bumps.h5'.format(workdir), 'w') as f:
            f['B'] = B
            f['bumpP'] = bumpP
    if args['plot']:
        plt.plot(B[1:], bumpP[1:], '*', markersize=10)
        plt.plot(B[1:], CODEbump[1:], 's', markersize=10)
        plt.legend(['DREAM', 'CODE'])
        plt.xlabel('Magnetic field strength (T)')
        plt.ylabel('Bump momentum ($mc$)')
        plt.tight_layout()
        plt.show()
    if success:
        dreamtests.print_ok("All runaway rates match those of CODE.")
    return success