Subroutine PLANE

Constrain a particle to the plane $ax+by+cz=d$. data should contain (/a, b, c, d/). Periodic boundary conditions are not applied (the plane only cuts the cell once, if at all).

The function used is $C = \frac{\mathbf{r} \cdot \mathbf{n} - d}{\vert\mathbf{n}\vert}$ where $\mathbf{n} = (a,b,c)$

pos, velo, data -- real(dp), dimension(:), intent(in)

t -- real(dp), intent(in)

C -- real(dp), intent(out)

dC_dr -- real(dp), dimension(size(pos)), intent(out)

dC_dt -- real(dp), intent(out)

BONDLENGTH_DIFF subroutine