Next, I do not know if you recognize this or not, but if you interept the solution (even though its running on cluster), you should be capable to get the NR residuals for the last 3 (I ordinarily use 3 but you can enter all the more) iterations. The solution https://a-solution.info/