add

mrdna/readers/segmentmodel_from_oxdna_pinyi.py

@@ -114,6 +114,9 @@ def mrdna_model_from_oxdna(coordinate_file, topology_file,virt2nuc=None,get_nt_p
                 pass
             counter+=1
         nt_prop["bp"]=bp
+        vh_bool=1-(nt_prop["vh"]%2)*2
+        is_scaf_bool=nt_prop["is_scaf"]*2-1
+        nt_prop["fwd"]=np.array((is_scaf_bool.T*vh_bool+1)/2,dtype=bool)
         return nt_prop
     try:
         nt_prop=get_virt2nuc(virt2nuc,top_data)