- wow, it runs now, but doesn't converge. Need to look into interpolation of "rho", which should be nearest Neighbor here.