Commit 30a24b44 authored by cmaffeo2's avatar cmaffeo2
Browse files

Fixed processing of push bonds when no crossovers exist between two adjacent helices

parent 39736a40
......@@ -999,7 +999,8 @@ class atomicModel():
if hid2 not in helixIds:
# print("WARNING: writeENM: helix %d not in helixIds" % hid2) # xo[%d] dict" % (hid2,hid))
dist = min([np.abs(z-zidx) for z in xo[hid][hid2]])
dists = [minDist,] + [np.abs(z-zidx) for z in xo[hid][hid2]]
dist = min(dists)
key = (hid2,zidx,isFwd)
if dist >= minDist and key in ntDict:
nt2 = ntDict[key]
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment