Method that will be needed for boundary conditions
Method that will be needed for boundary conditions
BCHandler should be done now
Wraper function for Neumann BC
PlanetPhysics uses BCHandler now