Acoustics with radial symmetric initial conditions. The solution should remain radially symmetric. First run the code in the 1drad subdirectory to compute the "true solution" and then setplot.py contains code to produce a scatter plot of the computed 2d pressure vs. distance from origin and compare with the 1d solution.