First working version of code. Single precision or double precision
selectable at compile time. Sheppard filter still needs work. Still need
to add variable time stepping.