Added loop
1/2 new model

Added a while loop to find what x max goes to the sigma max, added two
user functions for eta_0 and a function for BCs, a check that our BC
function is far enough away from the shore, and a loop to find when we
need to swap methods.