Changed BoundaryDirichlet_2ndOrder to calculate the guard cell from a 3rd order accurate interpolation, so that (first) derivatives should be second order accurate. Added a preprocessor flag which, if set, calculates by extrapolation a value for the second guard cells rather than setting them to 1.e60 (#ifdef CHECK) or not setting them at all (#else).