is_cyclic: determine which cells contain issues

[AquaPore]

Hi,

Thanks for this amazing software.

I am using the Wflow model and from the derived Local Drain Direction it is checking if there are issues in cyclic, by using is_cyclic . I am writing to determine if it may be possible to determine which cells are causing issues?

Many thanks for any help you may provide,
Joseph

    if is_cyclic(graph)
        error("""One or more cycles detected in flow graph.
            The provided local drainage direction map may be unsound.
            Verify that each active flow cell flows towards a pit.
            """)
    end

[henrik-wolf]

You should be able to run Graphs.simplecycles (or probably Graphs.simplecycles_iter if your graph is on the larger side) to get some idea about the cycles?