playing around with things, solver is not super stable, but I can get it to work. Needs a newton solver
added in convergence on tube diameter
cleaning up prints