def METHOD_NAME(args):
    """
    Run the test.
    """
    # Tolerance to require for agreement with CODE
    TOLERANCE = 1e-2
    success = True
    workdir = pathlib.Path(__file__).parent.absolute()
    T, Z = np.meshgrid([100, 1000, 10000], [1, 5, 10])
    nZ, nT = T.shape
    sigma_calc = np.zeros((nZ, nT))
    sigma_braams = np.zeros((nZ, nT))
    for i in range(0, nZ):
        for j in range(0, nT):
            print('Checking T = {} eV, Z = {:.1f}... '.format(T[i,j], Z[i,j]), end="")
            try:
                r, _, s_calc, s_braams = runTZ(T[i,j], Z[i,j])
                sigma_calc[i,j] = s_calc[-1,0]
                sigma_braams[i,j] = s_braams[-1,0]
            except Exception as e:
                print('\x1B[1;31mFAIL\x1B[0m')
                print(e)
                sigma_calc[i,j] = 0
                sigma_braams[i,j] = 0
                return False
            # Compare conductivity right away
            Delta = np.abs(sigma_calc[i,j] / sigma_braams[i,j] - 1.0)
            print("Delta = {:.3f}%".format(Delta*100))
            if Delta > TOLERANCE:
                dreamtests.print_error("DREAM conductivity deviates from Braams-Karney at T = {} eV, Z = {}".format(T[i,j], Z[i,j]))
                success = False
    # Save
    if args['save']:
        with h5py.File('{}/DREAM-conductivities.h5'.format(workdir), 'w') as f:
            f['T'] = T
            f['Z'] = Z
            f['sigma_braams'] = sigma_braams
            f['sigma_calc'] = sigma_calc
    if args['plot']:
        cmap = GeriMap.get()
        # Compare conductivities
        plt.figure(figsize=(5,3))
        legs = []
        legh = []
        hN = None
        for i in range(0, nT):
            clr = cmap(i/nT)
            h,  = plt.semilogy(Z[:,i], sigma_braams[:,i], color=clr, linewidth=2)
            hN, = plt.semilogy(Z[:,i], sigma_calc[:,i], 'x', color=clr, markersize=10, markeredgewidth=3)
            if T[0,i] >= 1e3:
                s = r'$T = {:.0f}\,\mathrm{{keV}}$'.format(T[0,i]/1e3)
            else:
                s = r'$T = {:.0f}\,\mathrm{{eV}}$'.format(T[0,i])
            yfac = 1.2
            plt.text(Z[-2,i]+4, sigma_calc[-2,i]*yfac, s, color=clr)
            legs.append(s)
            legh.append(h)
        legs.append('$\mathrm{DREAM}$')
        legh.append(hN)
        plt.xlabel(r'$Z$')
        plt.ylabel(r'$\sigma\ \mathrm{(S/m)}$')
        #plt.legend(legh, legs)
        plt.tight_layout()
        plt.show()
    if success:
        dreamtests.print_ok("The Ampère-Faraday equation yields the expected solution.")
    return success