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 = "CHEMICAL"
    wn.add_pattern(
        "NewPattern", [1.0]
    )  # start_time=0, end_time=wn.options.time.duration)
    wn.add_source("Source1", "River", "CONCEN", 100, "NewPattern")
    # WQ = wntr.scenario.Waterquality('CHEM', ['River'], 'CONCEN', 100, 0, -1)
    sim = wntr.sim.EpanetSimulator(wn)
    results = sim.run_sim(version=2.0)
    expected = 91661.72 * (1e-6 / 0.001)  # Node '159' at hour 6
    error = abs(
        (results.node["quality"].loc[6 * 3600, "159"] - expected) / expected
    )
    self.assertLess(error, 0.0001)  # 0.01% error