docs: def lnprior(pars):
      """
      Return probability of parameter values according to prior knowledge.
      Parameter limits should be done here through uniform prior ditributions
      """
      # Limit norm and B to be positive
      logprob = (
          naima.uniform_prior(pars[0], 0.0, np.inf)
          + naima.uniform_prior(pars[1], -1, 5)
          + naima.uniform_prior(pars[3], 0, np.inf)
      )
  
      return logprob