In the following cookbook example:

the variance along the profile is extracted:

= prof.variance['gas', 'velocity_magnitude'].value

but is then labelled as the standard deviation when the plot is created:

.loglog(radius, variance, label='Standard Deviation')

So, does "prof.variance" return the standard deviation and not the variance?  Or, is the plot mislabelled?  Should the plot instead be:

plt.loglog(radius, np.sqrt(variance), label='Standard Deviation')