at "$z_{(2)}(i)$" and "$z_{(2)}(1+i)$" : (Pari/GP)
\\ input to Pari/GP \\ output from Pari/GP
\\ -------------------------\\-----------------------------
t=exp(LambertW(log(I))) \\ %666 = 1.36062 + 1.11944*I
t^t \\ %667 = 1.00000*I
t=exp(LambertW(log(1+I))) \\ %668 = 1.39402 + 0.577732*I
t^t \\ %669 = 1.00000 + 1.00000*I
For $t = z_{(3)}(x)$ I used my series-representation. The value $x=\exp(-\pi/2)$ seems too much out of radius of convergence, so I could not yet obtain a result. But for $x=\exp(\pi/2) \approx 4.81048$ I could obtain $t \approx 1.77053$ and $t^{t^t} \approx 4.81048 $ with error of about $-0.0000000045...$. Here I had to use a procedure for summing divergent series (because the found powerseries at this argument diverges strongly), which is an adaption of the Noerlund-summation with a manually adapted "order" for the summation, but likely the better known Borel-summation should have been possible too.