divergence_field function

Numerical divergence of a vector field