Fields computed using curl discontinuous

Hi everyone,

when I try to derive the H field from an E field grid function using curl, I get discontinuities in the field. See the image below for the E and H field.

gfu_h = 1/(mu0*w)*curl(gfu_e)

Is this an expected behaviour or am I doing something wrong?

The problem was with the order of the fes. This only occurs for order=1.