Removed ideal gas specific SineWavePotential scripts and incorporated them in the isothermal script. Ideal gas now seems to work.

Added 2D and 3D SineWavePotential tests. Need to run them at higher resolution and for longer, but they seem to produce reasonable results.

Added a maximal smoothing length that can be used to prevent particles from becoming too big in vacuum cases near the box edges. Will break the accuracy of the hydrodynamic scheme but is necessary in some situations.

Added some extra isothermal eos flags in Gizmo/hydro.h. Made sure gpart mass is updated when conserved mass changes. Fixed a number of possible issues with the isothermal Riemann solver.

Particle velocities are now correct. Period. Noh still does not work properly, but increasing resolution helps (so it is probably just a resolution issue).

Added some extra vacuum checks. Made sure isothermal Riemann solver crashes on vacuum (since isothermal vacuum cannot exist). Changed force.v_full to v instead of v_full and tried to explain why. Still have a problem with dt in acceleration computation which needs to be solved; will do this tomorrow.

Improved initial density guess in isothermal Riemann solver, which leads to (faster) convergence. Added some extra isothermal eos flags in Gizmo hydro.h. Temporarily disabled primitive variable drift and changed the way the actual particle movement is calculated.

Added some extra checks to correctly handle vacuum + some debug checks. 2D Noh test now seems to work with moving particles.

Cleaned up GIZMO gravity a bit. Still doesn't work. Added a gradient interpolation limiter that should solve the Noh problem.

Compute the total mass in the system and add a variable to the gpart to record how much mass they interacted with via self-gravity.

Repaired GIZMO_SPH after the time integration changes that had been made.

The most important change is the addition of a new method, `hydro_timestep_extra`, that is called after the new time step for a particle has been computed, and that passes the physical time step for the next step to the particle.
In `GIZMO_SPH`, this method is used to
 - store the physical time step in `p->force.dt`. This value is used during the force loop to exchange fluxes between particles
 - reset a new `p->force.active` flag to an inactive value (`0`). This flag is set to an active value (`1`) in `hydro_init_part`, and is used during the force loop to decide whether or not to exchange flux with particle `j` in the asymmetric version of the interaction method.

For all other hydro schemes, this method is empty.

Other changes:
 - made sure `benchmarkInteractions` compiles with `GIZMO_SPH`
 - added a compilation flag, `GIZMO_FIX_PARTICLES`, that disables particle movement. If possible, this should be replaced by a corresponding parameter
 - added some sanity checks on density and pressure values, and a proper treatment of vacuum in the gradient prediction step

This version successfully runs most of the hydro examples (including the Sod tests, Sedov blasts, Kelvin-Helmholtz test, square tests, Gresho vortex), but fails to run the 2D and 3D Noh problem (1D works).

See merge request !315

Streamlined GIZMO_FIX_PARTICLES usage (tried to add parameter, but this does not work, as space_init has no access to the hydro_properties). Moved conserved variable initialization to hydro_first_init_part, and removed it from hydro_convert_quantities. Added some sanity checks to maintain positive density and pressure.

Cleaned up Gizmo fixes. Added a hydro_timestep_extra method and added an empty version for all hydro schemes.