20170921
8:20 11:40 2:50 + 0:30 = 3:20
Case 1
重新测试均匀流算例，发散，不知道哪里改错了。
查看在另一台服务器上的程序发现，昨天测试算例有误。
cell.fT[i][j] = aTP*cell.fBP[i][j] - bTP*cell.fT[i][j] + cell.DtSlashVolume*cell.fFlux[i][j];
cell.gT[i][j] = aTP*cell.gBP[i][j] - bTP*cell.gT[i][j] + cell.DtSlashVolume*cell.gFlux[i][j];
cell.fT[i][j] = aTP*cell.fBP[i][j] - bTP*cell.fT[i][j] + 4.0*dt*cell.fFlux[i][j];
cell.gT[i][j] = aTP*cell.gBP[i][j] - bTP*cell.gT[i][j] + 4.0*dt*cell.gFlux[i][j];
昨天的算例固定cell volume为0.25，所以dt/volume = 4.0*dt，但是1、2行的内容应该注释掉，昨天测试不发散的算例都没有注释1、2行，这就意味着，fT被更新了两次，结果错误。
Case 2
注释上述1、2行代码，更新网格文件，保证格点坐标无误差，重新测试均匀流，初始化为激波结构上游平衡态，164步之后监测点发散，与之前结果不同的是，fT, fh在46步到145值恒定为9.184982020200943e-68 ，gT, gh在52步到145值恒定为5.487069988109647e-68，之后出现振荡，进而发散。
14:20 18:00 3:10
Case 3
Case 31 ubuntu001
- 出现振荡的分布函数基本都位于速度空间的边缘，绝对值较大，在更新界面通量时，并没有乘以时间，而是在更新fT的时候将通量乘以dt/V，推测是fh乘以大数引入了机器误差，所以将dt移到xi_n_dS中，限制与fh相乘的值。
- 跳过宏观量更新，将xi_n_dS的值乘以小数dt，结果发散。
- 跳过宏观量更新，修改界面法矢量求解方法，不归一化，220步发散。
- 修改监测点。
Case 32 ubuntu002
- 跳过宏观量更新，界面分布函数和格心分布函数设为定值，系数保持恒定，fT在160步出现振荡，gT在219步未发生振荡。
- 添加通量输出函数。
- 添加监测点周围单元fT,gT的输出。
- 重新测试，监测点位置不佳，更换监测点（12.75,0.25,199,3）。
- 将xi_n_dS的绝对值设为1，符号保持不变，
if(xi·n >=0) 1 else -1
Case 33 ubuntu003
- 跳过通量更新，保持宏观量更新，系数保持恒定，2500步分布函数保持恒定。
20:20 22:30
