Drawing trials from asymmetric error distribution (Numpy)

I'm running a quick Monte Carlo error analysis. One of the parameters I'm looking at has a non-Gaussian error distribution; let's say Log10[param] = 9.33 +0.21 -0.01. In NumPy, how would I draw trials from a distribution centered at the parameter representing these asymmetric errors?