I'll expand on Nelson's proof a bit.
You know the Mean Value Property that the value of a harmonic function $u$ at $z$ equals the average of its values on any circle centered at $z$. We can turn this into a fact about averages over balls (which in this case are disks) just by using polar coordinates:
$$\int_{B(z,R)} u(w)\,dw = \int_0^R \int_0^{2\pi} u(z + r e^{i \theta}) r \,d\theta\,dr = \int_0^R 2 \pi u(z) r\,dr = u(z) \pi R^2.$$
Dividing both sides by $\pi R^2$, we see that $u(z)$ equals the average value of $u$ over the ball $B(z,R)$.
Now fix $z,w \in \mathbb{C}$, and consider disks $B(z,R)$, $B(w,R)$ centered at each of them. Let $L = B(z,R) \cap B(w,R)$ be the lens-shaped region lying in both disks. Let $L_z = B(z,r) \backslash L$ be the region in the $z$ disk but not the $w$ disk. (Draw a picture.) One can work out the area $A(L)$ of $L$ by a bit of multivariable calculus (or look it up on Mathworld). A bit more calculus shows that the ratio $\frac{A(L)}{\pi R^2} \to 1$ as $R \to \infty$, so that the lens occupies "most" of the disk. Since $A(L) + A(L_z) = \pi R^2$, we must therefore have $\frac{A(L_z)}{\pi R^2} \to 0$.
Now by the mean value property, for any $R$ we have
$$u(z) = \frac{1}{\pi R^2} \int_{B(z,R)} u(\zeta)\,d\zeta = \frac{1}{\pi R^2} \int_L u(\zeta)\,d\zeta + \frac{1}{\pi R^2} \int_{L_z} u(\zeta)\,d\zeta.$$
For the second term, if $|u|$ is bounded by some constant $M$, we have as $R \to \infty$,
$$\left|\frac{1}{\pi R^2} \int_{L_z} u(\zeta)\,d\zeta \right| \le \frac{M A(L_z)}{\pi R^2} \to 0.$$
Thus
$$u(z) = \lim_{R \to \infty} \frac{1}{\pi R^2} \int_L u(\zeta)\,d\zeta.$$
But if we replace $z$ with $w$ we can make exactly the same argument (replacing $L_z$ by $L_w$ whose area is the same, by symmetry). So we also have
$$u(w) = \lim_{R \to \infty} \frac{1}{\pi R^2} \int_L u(\zeta)\,d\zeta.$$
Thus $u(z) = u(w)$.
a. $e^{iz}$ satisfies $|e^{iz}|=e^{-\text{Im}(y)}<1$
b. $iz$ as Tim Raczkowski told you.
c. $e^{z}$ is bounded on the imaginary axis. If $iz\in\mathbb{R}$ then $|e^{z}|=1$.
Best Answer
$$ u(z) = \arg(z) \quad (0 < \arg(z) < \pi) $$ (the principal value of the argument) is harmonic in the upper half-plane because it is the imaginary part of the (principal value of the) logarithm.