[guided]We need a weight that rewards smallness on $K$ while penalising the support of the $\bar\partial$-error. Set
\begin{align*}
L := \operatorname{supp}a.
\end{align*}
If $L=\varnothing$, there is no error support to separate. In that case define
\begin{align*}
\psi: \mathbb C &\to \mathbb R \\
z &\mapsto -1,
\end{align*}
and set $\eta:=1$. Then $\psi$ is smooth and subharmonic, $\psi\le -\eta$ on $K$, and the condition $\psi\ge\eta$ on $L$ is vacuous. We therefore assume from now on in this construction that $L\ne\varnothing$.
The set $L$ is compact, and from the previous step it is disjoint from $\widehat K_{0,\mathrm{poly}}$. Since each point of $L$ lies outside the polynomial hull of $K_0$, polynomial convexity first gives only strict separation: for every $\zeta \in L$, there is a polynomial map
\begin{align*}
h_\zeta: \mathbb C \to \mathbb C
\end{align*}
such that
\begin{align*}
T_\zeta:=|h_\zeta(\zeta)|>S_\zeta:=\sup_{z\in K_0}|h_\zeta(z)|.
\end{align*}
To get the stronger numerical constants needed later, we first amplify the strict inequality by taking powers. Choose $m_\zeta\ge1$ so that
\begin{align*}
T_\zeta^{m_\zeta}>4S_\zeta^{m_\zeta};
\end{align*}
if $S_\zeta=0$, this holds with $m_\zeta=1$, while if $S_\zeta>0$ it follows from $(T_\zeta/S_\zeta)^{m_\zeta}\to\infty$. Now choose $\lambda_\zeta>0$ with
\begin{align*}
\frac{2}{T_\zeta^{m_\zeta}}<\lambda_\zeta\le \frac{1}{2S_\zeta^{m_\zeta}},
\end{align*}
interpreting the upper bound as $+\infty$ when $S_\zeta=0$. The interval is non-empty because $T_\zeta^{m_\zeta}>4S_\zeta^{m_\zeta}$. Define
\begin{align*}
q_\zeta: \mathbb C &\to \mathbb C \\
z &\mapsto \lambda_\zeta h_\zeta(z)^{m_\zeta}.
\end{align*}
Then the upper bound on $\lambda_\zeta$ gives
\begin{align*}
\sup_{z \in K_0}|q_\zeta(z)| \le \frac{1}{2},
\end{align*}
and the lower bound gives
\begin{align*}
|q_\zeta(\zeta)| > 2.
\end{align*}
By continuity, the inequality $|q_\zeta|>1$ holds on an open neighbourhood of $\zeta$. Compactness of $L$ gives finitely many points $\zeta_1,\dots,\zeta_N \in L$ such that for every $z \in L$ there is some $j$ with
\begin{align*}
|q_{\zeta_j}(z)|>1.
\end{align*}
The finite number $N$ may be large, so the estimates must be normalised before summing. Choose an integer $M\ge1$ such that
\begin{align*}
N4^{-M}<\frac14,
\end{align*}
and define polynomial maps
\begin{align*}
r_j: \mathbb C &\to \mathbb C \\
z &\mapsto q_{\zeta_j}(z)^M
\end{align*}
for $j\in\{1,\dots,N\}$. On $K_0$, each $|q_{\zeta_j}|<1/2$, hence
\begin{align*}
\sum_{j=1}^N |r_j(z)|^2
=
\sum_{j=1}^N |q_{\zeta_j}(z)|^{2M}
<
N4^{-M}
<
\frac14.
\end{align*}
On $L$, at least one $|q_{\zeta_j}(z)|$ is greater than $1$, so the corresponding $|r_j(z)|$ is also greater than $1$.
Now define
\begin{align*}
\Psi: \mathbb C &\to \mathbb R \\
z &\mapsto \log\left(\sum_{j=1}^N |r_j(z)|^2+\frac{1}{4}\right).
\end{align*}
The function $\Psi$ is smooth because the quantity inside the logarithm is everywhere at least $1/4$. It is subharmonic because logarithms of sums of squared moduli of holomorphic functions, with a positive constant added, are subharmonic.
On $K_0$ we have
\begin{align*}
\sum_{j=1}^N |r_j(z)|^2+\frac14 < \frac12,
\end{align*}
and hence
\begin{align*}
\Psi(z)<\log\left(\frac12\right)
\end{align*}
for $z \in K_0$. On the other hand, if $z \in L$, then one of the polynomials $r_j$ has modulus greater than $1$, so
\begin{align*}
\sum_{j=1}^N |r_j(z)|^2+\frac14 > \frac54,
\end{align*}
and therefore
\begin{align*}
\Psi(z)>\log\left(\frac54\right).
\end{align*}
Finally choose constants $A>0$ and $B \in \mathbb R$ so that
\begin{align*}
\psi(z):=A\Psi(z)-B
\end{align*}
satisfies, for some $\eta>0$,
\begin{align*}
\psi(z)\le -\eta \quad \text{on } K,
\qquad
\psi(z)\ge \eta \quad \text{on } L.
\end{align*}
Affine rescaling by a positive factor preserves subharmonicity, so $\psi$ is still a smooth subharmonic function.[/guided]