def METHOD_NAME(self):
    inp_file = join(datadir, "Net3.inp")
    wn = wntr.network.WaterNetworkModel(inp_file)
    wn.options.time.hydraulic_timestep = 15*60
    wn.options.time.quality_timestep = 15*60
    wn.options.time.report_timestep = 15*60
    wn.options.quality.parameter = "TRACE"
    wn.options.quality.trace_node = "121"
    # WQ = wntr.scenario.Waterquality('TRACE', ['121'])
    sim = wntr.sim.EpanetSimulator(wn)
    results = sim.run_sim(version=2.0)
    # print(results.node.keys())
    expected = 91.66  # Node '159' at hour 6
    error = abs(
        float(results.node["quality"].loc[6 * 3600, "159"] - expected)
        / float(expected)
    )
    self.assertLess(error, 0.0001)  # 0.01% error