Issues when running the model

[geraec]

Hi,

When trying to run the animation (animate_spme) after running the spme I am having the following error message:
Error using matlab.graphics.axis.Axes/set
While setting the 'YLim' property of 'Axes':
This is not a valid LimitsWithInfs value. Complex inputs are not supported

Error in ylim (line 43)
set(ax,'ylim',val);

Error in animate_spme (line 64)
ylim([min(V), max(V)]);

After a bit of digging I found that the values of solid concentration c_ss_n and c_ss_p are negative therefore I get complex number when computing the voltage.

I managed to run a simulation by reducing the current with a factor of 0.0000009.

Thanks in advance