Add a program called hydro, which should become a 1d Euler solver

Right now it just identifies and adds ghost nodes.
Have 1d Euler solver working: 1st order in space and time, uniform mesh
Show that Euler solver can work over jumps in refinement