今天对程序大修了一下，顺便把所有算例测试了一遍，突然发现二维浅水方程有些算例出现了明显的错误。
这次突然出现的错误让我有点措手不及，因为一直没有修改过浅水方程求解器，所以这些算例很久没有测试过了。硬着头皮从斜率限制器，数值通量，计算项等等所有计算函数里面寻找 bug，但是找了一下午仍然一无所获。看哪里都好像不太对，但是不管怎么改错误仍然会出现。
好在上个月我写算例时候当时是能够正确计算的，于是利用 git reset 命令开始返回过去提交的代码版本，开始一个个测试，最后终于在6月2号版本程序里得到了正确的计算结果。这也就是说，计算的 bug 在2号后的版本里就开始出现了。下面就是比较后面提交的程序做了哪些修改。感觉我当时还挺勤奋的，两天时间改了一千多行程序，看的我是头晕眼花。又经过整整一晚的搜索，终于在判断单元干湿的函数中找到了可能存在的逻辑错误。把修改的代码用正确代码替换，终于得到正确计算结果。
真的是感谢 git 这款优秀的版本控制系统，节省了我不知道多少时间。