W 70 "Domain width"
dW 9 "Perfectly matched layer width"
c0 1 "Speed of sound"
rho0 1 "Background density"
p0 1[atm]/(c0^2*rho0) "Background pressure"
alpha log(2)/9 "Gaussian pulse width"
beta 0.04 "Gaussian pulse velocity factor"