I have searched all around for an answer to this, what I would imagine to be a rather natural question, but to no avail. How does one extend the Lindelöf maximum principle to sets of measure 0? Here's a precise formulation, which has a more restrictive condition on the domain of the harmonic function:
Let $f: \frac{1}{2}\mathbb{D} \to \mathbb{R}$ be harmonic and bounded. If $\limsup_{r\nearrow 1/2} f(re^{i\theta}) \leq 0$ for almost all $\theta$, then show $f \leq 0$ on $\frac{1}{2}\mathbb{D}$.
My attempt is as follows; I have tried to follow the usual technique used to prove the theorem when the exceptional set is a single point: let $E$ be the exceptional subset of the boundary, and let $S_{1, \epsilon}, S_{2, \epsilon}, \ldots, S_{N_\epsilon, \epsilon}$ be a collection of connected subsets of $\partial \frac{1}{2}\mathbb{D}$ whose lengths sum to $\epsilon$ and such that $E\subseteq \bigcup_{n=1}^{N_\epsilon} S_{n, \epsilon}$ (from the definition of Lebesgue measure). These naturally give rise to points $z_{i, \epsilon}$ and radii $r_{i, \epsilon}$ such that $S_{i, \epsilon} \subseteq \{|z-z_{i, \epsilon}| < r_{i, \epsilon}\}$ and such that $\sum_i r_{i, \epsilon} \leq \epsilon/2$. Then the function $$f_\epsilon :=f(z) + \sum_{i=1}^{N_\epsilon} \frac{1}{-\log(r_{i, \epsilon})} \log|z-z_{i, \epsilon}|$$ is a harmonic function which is less than 0 as one approaches any boundary point, and so is 0 throughout the disk. The problem is though that our modified function does not in general converge pointwise to $f$ as $\epsilon \searrow 0$ (since there is a priori no bound on $N_\epsilon$), which is what makes the usual proof of the Lindelöf maximum principle for a single point work.
Any pointers on how to adjust this proof or go about it differently would be greatly appreciated.
The other answer is probably better than this one; potential theory, polar sets, etc get to the heart of the matter. But for the disk specifically one can give a very elementary argument. Here "elementary" means "just using things in for example Rudin Real and Complex Analysis".
If $u$ is a bounded harmonic function in $\mathbb D$ then $f(e^{it})=\lim_{r\to1^-}u(re^{it})$ exists for almost every $t$. It's clear that $f$ is measurable, so $f\in L^\infty(\Bbb T)$. And in fact $u=P[f]$, the Poisson integral of $f$. Since the Poisson kernel is positive it's clear that $f\le0$ almost everywhere implies $u\le 0$.
(That's a proof for bounded harmonic functions in $\Bbb D$; I'm totally stumped on how to transfer it to a proof for bounded harmonic functions in $\frac12\Bbb D$, sorry.)