Purpose

---

  This procedure computes the curl of a 3d vector field sampled on a
  uniform grid.

Inputs

---

  x: The x component of the vector field
  y: The y component of the vector field
  z: The z component of the vector field

Keyword Parameters

---

  order: Sets the precision of the lagrange interpolation when
  computing partial derivatives. see pdiv.pro.

Outputs

---

  cx: The x component of the curl
  cy: The y component of the curl
  cz: The z component of the curl

Modification History

---

  2010-08-03: Written by Chris Beaumont