Thanks for your reply.A couple things to try. If this is right at initialization, you should look at your initial state information. What are you using to force Icepack? You can try reducing the timestep (dt). You can also try increasing the Picard iterations (nit_max). This requires a change to the code. However, changing the number of iterations is unlikely to help. The other thing you might try is to try a roundoff perturbation to an icepack parameter. For example, you could add 1.0e-8 to the r_snw value.
Dave