Added doxumentation for the Fox-Kemper scheme

- Re-wrote most of the documentation for mixed_layer_restrat(),
  adding equations, symbols for parameters and references.
- This is prelude to updating how we do the buoyancy gradient
  upscaling.