Currently I use Gromacs to perform molecular dynamics of tobermorite 11A, but it turns out that this structure has a monoclinic box, I need to solvate my structure only above in the z axis, but I have problems, because gromacs does not allow to solvate with angles, so some molecules of solvent are located outside the box .... how could I do this?