Hello, I noticed that bonds are sometimes stored in top_global->moltype->ilist[F_CONSTR] instead of top_global->moltype->ilist[F_BONDS], What causes this? Regards, Vedran -- Vedran Miletić http://vedranmileti.ch/