  We investigate the statistics of flat-top solitary wave parameters in the
presence of weak multiplicative dissipative disorder. We consider first
propagation of solitary waves of the cubic-quintic nonlinear Schr\"odinger
equation (CQNLSE) in the presence of disorder in the cubic nonlinear gain. We
show by a perturbative analytic calculation and by Monte Carlo simulations that
the probability density function (PDF) of the amplitude $\eta$ exhibits
loglognormal divergence near the maximum possible amplitude $\eta_{m}$, a
behavior that is similar to the one observed earlier for disorder in the linear
gain [A. Peleg et al., Phys. Rev. E {\bf 72}, 027203 (2005)]. We relate the
loglognormal divergence of the amplitude PDF to the super-exponential approach
of $\eta$ to $\eta_{m}$ in the corresponding deterministic model with
linear/nonlinear gain. Furthermore, for solitary waves of the derivative CQNLSE
with weak disorder in the linear gain both the amplitude and the group velocity
$\beta$ become random. We therefore study analytically and by Monte Carlo
simulations the PDF of the parameter $p$, where
$p=\eta/(1-\varepsilon_s\beta/2)$ and $\varepsilon_s$ is the self-steepening
coefficient. Our analytic calculations and numerical simulations show that the
PDF of $p$ is loglognormally divergent near the maximum $p$-value.
