Work toward improving BC's: put in a boundary_conditions_flux method
Implement boundary_conditions_emf methods