Add operators that support gradient calculations for computing the strain on collocated grids

When computing the strain on collocated grids, we need a few extra gradient operators. These essentially compute the gradient of quantities that live at the "corners" of cell-centered control volumes.