Solution (source code)

= Solution

At the source, $\zeta=z_0=\mathcal B_0/\beta$. The initial width required to lie exactly on the similarity solution is therefore
$$
\boxed{b_0=\frac{3E}{4}z_0
=\frac{3E\mathcal B_0}{4\beta}
=\frac{3E\rho_0c_pT_0\,\mathcal B_0}{4gq_0}}.
$$
The corresponding $Q_0=A z_0^{4/3}$ and $M_0=C z_0^{5/3}$ are the compatible initial fluxes stipulated in the question.