Made all model arrays (i.e. xx, yy, zz, znum, and znum_r) into arrays of size natoms*2 so that we never have to reallocate them. This will save time.