Change `domain_stats` to use `spatial_average`

[brianpm]

New feature type

New plot and/or plot enhancement

What is this new feature?

During documentation of plotting_functions.py, I noted that the function domain_stats is making some assumptions that could be avoided by using spatial_average instead.

Assistance required?

No, I will make a PR when the feature is ready

Extra info

No response