Troubleshooting 2D riemann solver
 

Hi all, I am writing a 2D FV code in fortran 90 to solve the Compressible Euler Equations to solve a standard Riemann problem on a Cartesian grid. I think I have most of it down, but am not sure what to change to fix my flux function call when evaluating the residual. Currently I am doing first order upwind reconstruction with forward Euler time marching, but I wrote it in this way to allow further time schemes and flux functions to be added in the future. Please try to compile the attached code with gfortran and see where I am going wrong with my Cartesian solver. Thanks.


If there are other glaring issues, please point them out, but I mainly wish to have it compile and run successfully to start.