Dear all, I am trying to apply a spherical flat-bottomed potential to maintain some molecules within a certain region. By default, GROMACS uses the initial coordinates as reference. Is there any way to provide different references? Many thanks Salvatore