请各位高手一定帮忙看看这个二维流体程序.初写程序,没经验,不知道问题出那里?

试写程序

我写的这个程序很怪,可以单独算X方向和Y方向的Riemann问题,也可以算圆型爆炸问题,但算2维的Riemann问题就无能为力,实在不知道问题出在来那里?请各位高手一定帮忙看看.先谢谢了.
程序说明,初始化子程序INIT(整数变量),此整数变量取0和1分别是两个2维的Riemann问题事例,取2和3分别是X方向和Y方向的1维Riemann问题事例,取4是一个圆型爆炸问题事例.
子程序HYDRO(整数变量),此整数变量取1是用HLLE计算数值流,取2是用Lax_Freiderich计算数值流.
QQ:554354797
userpar.inc
parameter(nx=100,ny=100,md=4,eqN=4)