added 2nd order accuracy in space (defect correction) for implicit scheme. Changed defaults.
started to fix crank-nicholson and bdf2 time integration methods