A nonsmooth Newton multigrid method for a hybrid, shallow model of marine ice sheets