also group each sorbate pair automatically, and add a handy switch to print xyz files to visulize model structure