code: def log10norm(x, mu, sigma=1.0):
      
      return stats.lognorm(sigma * np.log(10), scale=mu).pdf(x)
docs: Scale scipy lognorm from natural log to base 10
      x     : input parameter
      mu    : mean of the underlying log10 gaussian
      sigma : variance of underlying log10 gaussian