Adding compute_method to junction boundaries
Experimenting with getting the beast to run. Problem
is that parts are going dry. Just recently forced a stage compute_method
for junction boundaries. The way these are imposed is not good
for very shallow / steep channels.