I am trying to simulate contact angle of water on a quartz slab using LAMMPS. I have created and equilibriated the quartz slab and the water cluster separately. Now I want to create a composite geometry that has the water cluster placed on the quartz. For that I wrote a Python script that modifies the output trajectory files by bringing the COMs of both the geometries to (0,0,0) and then puts the water cluster on the slab by adjusting the z coordinates. By that does not work as LAMMPS gives an "Out of range atoms error". Is there any other way/software with which this composite geometry can be created ?