Compute the gradient of a field defined at the GLL points.
This function computes the derivatives of the N + 1 Lagrange polynomials (N
being the polynomial order) to compute the derivative of a field that is
defined at the GLL points.
RETURNS
Return type: dict[str, npt.NDArray]
Dict containing the derivatives of the fields with respect to each dimension. Each derivative has the same shape as the input field.