Harmonic approximation by finite sums of moduli

Let $h(B_d)$ denote the space of real-valued harmonic functions on the unit ball $B_d$ of ${\mathbb{R}}^d$, $d\ge 2$. Given a radial weight $w$ on $B_d$, consider the following problem: construct a finite family ${f_1, f_2, \dots, f_J}$in$h(B_d)$suchthatthesum$|f_1| + |f_2|+\dots + |f_J|$ is equivalent to $w$. We solve the problem for weights $w$ with a doubling property. Moreover, if $d$ is even, then we characterize those $w$ for which the problem has a solution.