def METHOD_NAME(self):
    inp_file = join(netdir, "Net1.inp")
    wn = wntr.network.WaterNetworkModel(inp_file)
    earthquake = wntr.scenario.Earthquake((40, 55), 5, 10000.0)
    """
    Junction 12 is at (50,70)
    """
    R = earthquake.distance_to_epicenter(wn, wntr.network.Junction)
    expected = np.sqrt(np.power(50 - 40, 2) + np.power(70 - 55, 2))
    self.assertLess(np.abs(R["12"] - expected), 1e-5)
    wn = wntr.morph.scale_node_coordinates(wn, 100)
    R = earthquake.distance_to_epicenter(wn, wntr.network.Junction)
    expected = np.sqrt(np.power(5000 - 40, 2) + np.power(7000 - 55, 2))
    self.assertLess(np.abs(R["12"] - expected), 1e-5)