interpolation is key. Old 2 dimensional one was not right at the
boundary. New one is. Now things converge!