今天是第二次找恶心着急地找bug，上一次是将程序改为MPI程序。关于MPI程序调试的经验可以见上次的文章
今天的目标是
（1）清晰地了解各种边界处理方式，进而清晰地了解Color gradient模型关于三相点的处理方式。（2）统一化个性化自己的程序，提高通用性和可读性。（3）增加f_pc，使边界处理更清晰。（4）迁移采用可判断固体点的方式（事实证明，聚式和散式都可以）
得到的教训
（1）调程序可以采用部分替代法，如果用正确的程序块替代能得到正确结论，极大可能是这一块出错，应该先把这部分内容一点一滴地分析清楚，再考虑大面积分布调试法
（2）另一种笨办法是分步调试法，从第一步循环开始调试，输出宏观量、分布函数（包括f0），中间量f_pc（碰撞后的分布函数）。一块块调试。
