solve_ functions should probably reset initial condition #11
Description
Repeatedly calling solve_* at the moment without re-initializing the SubductionProblem class will increment the solution after each call. These should be reset at the start of each solve_ call.