Radial ring initial condition with 1.6 outside and 1.2 inside the ring in qinit. More commented versions of surface_tension_utils and setrun