Added count_bonded_adj_atoms function which will create two numpy arrays that will be referenced by saturate function later.