Fixed a bug
Only first order is 100% working, 2nd order does not work during subsonic flow, Cl is non zero for 2nd order subsonic, have to recheck it
Added additional test
Changed the definition of res_rho in jst