Secondly, I don't know if you realize this or not, but for those who interept the answer (even though its running on cluster), you should be in a position to get the NR residuals for the last 3 (I typically use 3 however you can enter even more) iterations. The https://a-solution.info/