Hi, Just a quick question - in do_coul() in mkinl_interactions.c, should the line nflop = 1; actually read nflop += 1; or not ? I don't think it makes much difference either way, but ... Thanks, Vishal