bugfix plus some work in progress
Gave velocity argument of get_hydro_state_at_point function a default of 0 for Gadget2 and Fi.
thermal model and tests scripts