Making gradients output optional since that may not be needed
and can require a lot of memory for topologically structured
grids. Also adding in divergence option.