Temperature-dependent stable phase domains of Zr–C–N (Zr–ZrC–ZrC0.96–Zr(C1−xNx)−ZrN, x = 0.1–0.9) systems from ab initio calculations and experimental results