平面上的超声速流（有粘性边界层和激波）

代码说明：

应用背景程序求解完全的二维N-S方程，考虑了粘性和导热。当来流Mach数为4时，计算结果与书上的结果一致，然而，我的程序不能计算Mach数为25的情况（即书中提供的另一种工况的比较结果），经过仔细检查，我认为自己的程序没毛病，可能是书中求解步长的公式给错了（eq.(10.16) in p.457），还望大家帮助检查一下。最后还要说明的是，该章在指导你写程序时，书上有些关键地方有印刷错误，会误导你，一定要小心。参考文献：Chapter 10, Computational Fluid Dynamics, The Basics with Applications, John D. Anderson, JR. McGraw-Hill, 2002, 4 《Anderson计算流体力学入门》 关键技术         bc.f---------------------边界条件          conver.f-----------------迭代收敛判断          dynvis.f-----------------计算粘性项          mac.f--------------------MacCormack求解器          mdot.f-------------------质量守恒判断          output.f-----------------结果输出          qcx.f,qcy.f--------------计算导热项          tauxx.f,tauyy.f,tauxy.f--计