Poisson These routines are meant to solve the Poisson equation in 2D (three dimensions to come) for use in calculating the self-gravity of the fluid in DISCO. -> incomplete