+ call int_bounds(nres_molec(2)-2,ibond_nucl_start,ibond_nucl_end)
+ ibond_nucl_start=ibond_nucl_start+nnt_molec(2)-1
+ ibond_nucl_end=ibond_nucl_end+nnt_molec(2)-1
+ print *,"NUCLibond",ibond_nucl_start,ibond_nucl_end
+ print *, "before devision",nnt_molec(2),nct_molec(2)-nnt_molec(2)
+ call int_bounds(nct_molec(2)-nnt_molec(2),ibondp_nucl_start,ibondp_nucl_end)
+ ibondp_nucl_start=ibondp_nucl_start+nnt_molec(2)
+ ibondp_nucl_end=ibondp_nucl_end+nnt_molec(2)
+ print *,"NUCLibond2",ibondp_nucl_start,ibondp_nucl_end
+
+
+ call int_bounds1(nres_molec(1)-1,ivec_start,ivec_end)