# Define U

# EL = (Ex + 1i*Ey)/sqrt(2);
# ER = (Ex - 1i*Ey)/sqrt(2);
# E_circular = U * E_cartesian, therefore:

U = [ 1, 1i, 0;
      1, -1i, 0;
      0, 0, sqrt(2) ] / sqrt(2);

select("kerr_attribute");
set("U",U);