bottleneck detection now obeys direction option
@@ -153,6 +153,19 @@ def map_mab(coords, mask, output, *args, **kwargs):
break
if bottleneck:
+ if direction[n] < 0:
+ if coord == flipdifflist[n]:
+ holder = bottleneck_base + 2 * n
+ special = True
+ break
+
+ if direction[n] > 0:
+ if coord == difflist[n]:
+ holder = bottleneck_base + 2 * n
+ special = True
+ break
+
+ if direction[n] == 0:
if coord == difflist[n]:
holder = bottleneck_base + 2 * n
special = True
