@@ -0,0 +1,871 @@
+import pytest%0A%0Afrom landlab import RasterModelGrid%0Afrom landlab.components import FlowAccumulator%0Afrom landlab.plot import analyze_channel_network_and_plot%0Afrom landlab.plot.channel_profile import channel_nodes%0A%0Adef test_route_to_multiple_error_raised():%0A    mg = RasterModelGrid((10, 10))%0A    z = mg.add_zeros('node', 'topographic__elevation')%0A    z += mg.x_of_node + mg.y_of_node%0A    fa = FlowAccumulator(mg, flow_director='MFD')%0A    fa.run_one_step()%0A%0A    with pytest.raises(NotImplementedError):%0A        analyze_channel_network_and_plot(mg)%0A%0A    with pytest.raises(NotImplementedError):%0A        channel_nodes(mg,%0A                      mg.at_node%5B'topographic__steepest_slope'%5D,%0A                      mg.at_node%5B'drainage_area'%5D,%0A                      mg.at_node%5B'flow__receiver_node'%5D,%0A                      number_of_channels=1,%0A                      threshold=1.0)%0A

test of channel profile