[constants]
gamma = 1.4
S = 13.5
M = 0.4
R = 1.5

[solver]
system = euler

[solver-time-integrator]
scheme = rk4
tstart = 0
tend = 50
dt = 0.005
controller = none

[solver-interfaces]
riemann-solver = rusanov

[solver-interfaces-line]
flux-pts = gauss-legendre

[soln-plugin-nancheck]
nsteps = 50

[soln-plugin-writer]
dt_out = 10
basedir = .
basename = euler_vortex_2d-{t:.1f}

[soln-ics]
rho = pow(1 - S*S*M*M*(gamma - 1)*exp(2*%(f)s)/(8*pi*pi), 1/(gamma - 1))
u = S*y*exp(%(f)s)/(2*pi*R)
v = 1 - S*x*exp(%(f)s)/(2*pi*R)
p = 1/(gamma*M*M)*pow(1 - S*S*M*M*(gamma - 1)*exp(2*%(f)s)/(8*pi*pi), gamma/(gamma - 1))
f = ((1 - x*x - y*y)/(2*R*R))