u0 uj*exp(-log(2)*(y/b)^2) Background x-velocity v0 0[m/s] Background y-velocity rho0 1/(-0.5*(gamma-1)/(gamma*p0)*(u0-uj)*u0+1/rhoj*u0/uj+1/rho_inf*(uj-u0)/uj) Background density T0 p0/(rho0*Rs) Background temperature S A*exp(-Bx*x^2-By*y^2) Gaussian source