Next, I don't know if you already know this or not, but for those who interept the answer (even though its running on cluster), you have to be ready to get the NR residuals for the final three (I ordinarily use 3 but you can enter far more) iterations. The https://a-solution.info/