changed degY to Y in dfdk and H for degenerate modes. This makes Jacobian spherically symmetric, at the expense of calculating an extra term in the Jacobian.