def METHOD_NAME(args):
    """
    Run the test.
    """
    global NR
    # Tolerance to require for agreement with CODE
    TOLERANCE = 1e-2
    success = True
    workdir = pathlib.Path(__file__).parent.absolute()
    T = np.array([1e1,1e2,1e3,1e4])
    nT = T.size
    jKinetic = np.zeros((nT, NR))
    jFluid   = np.zeros((nT, NR))
    for i in range(0, nT):
        print('Checking T = {} eV... '.format(T[i]), end="")
        try:
            jKinetic[i,:], jFluid[i,:], eps = runT(T[i])
        except Exception as e:
            print('\x1B[1;31mFAIL\x1B[0m')
            print(e)
            #continue
            return False
        # Compare conductivity right away
        Delta = np.amax(np.abs(jKinetic[i,:] / jFluid[i,:] - 1.0))
        print("Delta = {:.3f}%".format(Delta*100))
        if Delta > TOLERANCE:
            dreamtests.print_error("DREAM conductivity deviates from CODE at T = {} eV".format(T[i]))
            success = False
    # Save
    if args['save']:
        with h5py.File('{}/DREAM-conductivities-trapping.h5'.format(workdir), 'w') as f:
            f['T'] = T
            f['jFluid']   = jFluid
            f['jKinetic'] = jKinetic
    if args['plot']:
        cmap = GeriMap.get()
        # Compare conductivities
        #plt.figure(figsize=(9,6))
        plt.figure(figsize=(5,3))
        legs = []
        legh = []
        hN = None
        for i in range(0, nT):
            clr = cmap(i/nT)
            h,  = plt.semilogy(eps, jKinetic[:,i], color=clr, linewidth=2)
            hN, = plt.semilogy(eps, jFluid[:,i], 'x', color=clr, markersize=10, markeredgewidth=3)
            if T[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 = .6
            plt.text(0.5, jKinetic[0,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("All conductivities match those of CODE.")
    return success