Adding higher order Dirichlet boundary conditions

o Renamed "dirichlet_newo2" to "dirichlet_o3" since it extrapolates
  the boundary point to third order

o Added "dirichlet_o4" which goes to 4th-order. This then allows
  second-order accurate calculation of second derivatives near the boundary.