def METHOD_NAME(T, Z=1, n=5e19):
    """
    Generate appropriate DREAM settings.
    T:    Electron temperature.
    Z:    Effective charge of plasma.
    n:    Electron density.
    """
    a = 0.35
    nr = 20
    ds = DREAMSettings()
    ds.collisions.lnlambda = Collisions.LNLAMBDA_THERMAL
    # Get initial electric field profile
    r0, E0, lmbd = getInitialElectricField()
    ds.eqsys.E_field.setType(Efield.TYPE_SELFCONSISTENT)
    ds.eqsys.E_field.setInitialProfile(efield=E0, radius=r0*a)
    ds.eqsys.E_field.setBoundaryCondition(bctype=Efield.BC_TYPE_SELFCONSISTENT, inverse_wall_time=0, R0=1)
    ds.eqsys.n_i.addIon(name='Ion', Z=Z, n=n/Z, iontype=IonSpecies.IONS_PRESCRIBED_FULLY_IONIZED)   # Imaginary ion with charge Z
    ds.eqsys.T_cold.setPrescribedData(T)
    ds.eqsys.n_re.setDreicer(False)
    ds.eqsys.n_re.setAvalanche(False)
    ds.eqsys.j_ohm.setConductivityMode(JOhm.CONDUCTIVITY_MODE_BRAAMS)
    ds.hottailgrid.setEnabled(False)
    ds.runawaygrid.setEnabled(False)
    ds.radialgrid.setB0(1)
    ds.radialgrid.setMinorRadius(a)
    ds.radialgrid.setWallRadius(a)
    ds.radialgrid.setNr(nr)
    # Expected e-folding time scale is mu0*sigma
    t0  = scipy.constants.mu_0 * evaluateBraamsConductivity(n, T, Z) / (lmbd[0]/a)**2
    dt  = 1e-2 * t0
    ds.timestep.setTmax(10*t0)
    ds.timestep.setDt(dt)
    ds.other.include('fluid/conductivity')
    return ds, lmbd