Square-root-radius diffusion transform (source code)

= Square-root-radius diffusion transform
{title2=$x=\sqrt r,\quad g=\sqrt r\bar\nu_0\Sigma_1,\quad g_t=3q\bar\nu_0g_{xx}/(4r)$}

For a steady zero-mass-flux background and constant <viscous transport response exponent>, this change removes the first derivative in the radial diffusion operator. If $\Sigma_0\propto r^{-p}$, the <kinematic viscosity> scales as $r^{p-1/2}$, so the transformed diffusivity scales as $x^{2p-3}$. It is constant at $p=3/2$.