Add options for time adaptivity

My first crack at integration with the libMesh TwostepTimeSolver