Secondly, I do not know if you realize this or not, but if you interept the answer (even though its functioning on cluster), you should be capable to get the NR residuals for the last 3 (I usually use 3 but you can enter a lot more) iterations. The solution https://a-solution.info/