mfix的Negative gas density报错解决

错误很难定位，因为编译正常，而是运行过程中会告知出现Negative gas density，并且不收敛，没有其他错误信息。最后通过一步步定位发现是由于 IC_EP_g(1) < EP_star 导致的，也就是固相的堆积体积分数大于这个临界值导致，只要把 EP_star 调小一些，就好了。

从官方tutorial可以看出，EP_star只是决定了在稀相区和密相去用的不同的应力模型，至于为什么初始条件不能低于这个值，可能得进一步研究代码才能知道了。如果有大神能够指出问题所在，将不甚感激。