attempt to modularise form construction
adaptive timestepping
some changes to default set up
plot frequency / plot to file
increased eps to stop failures
changes to mapping for P2 functions
now can set input for q, h and phi
working mms (su untested) and taylor test
new functional for deposit gradient