a

mrdna/readers/segmentmodel_from_oxdna_pinyi.py

@@ -127,7 +127,7 @@ def mrdna_model_from_oxdna(coordinate_file, topology_file,virt2nuc=None,get_nt_p
         
         non_stack_ind,=np.where(nt_prop["stack"]==-1)
         for i in non_stack_ind:
-            zid=int(nt_prop.loc[i][zid])+int(nt_prop.loc[i]["fwd"])*2-1
+            zid=int(nt_prop.loc[i]["zid"])+int(nt_prop.loc[i]["fwd"])*2-1
             try:
                 nt_prop["stack"][i]=bp_map[(nt_prop.loc[i]["vh"],str(zid),nt_prop.loc[i]["fwd"])]
             except:

mrdna/readers/segmentmodel_from_scadnano.py

@@ -122,7 +122,7 @@ def gen_prop_table(fname):
     nt_prop["bp"]=bp
     non_stack_ind,=np.where(nt_prop["stack"]==-1)
     for i in non_stack_ind:
-        zid=int(nt_prop.loc[i][zid])+int(nt_prop.loc[i]["fwd"])*2-1
+        zid=int(nt_prop.loc[i]["zid"])+int(nt_prop.loc[i]["fwd"])*2-1
         try:
             nt_prop["stack"][i]=bp_map[(nt_prop.loc[i]["vh"],str(zid),nt_prop.loc[i]["fwd"])]
         except: