diff --git a/imagine/likelihoods/ensemble_likelihood/ensemble_likelihood.py b/imagine/likelihoods/ensemble_likelihood/ensemble_likelihood.py
index 90cf4ee5e90b80754aaa2263ec7eb1adc90678a9..1835a48e960bd27fbd76751b945841dc5bd842d9 100644
--- a/imagine/likelihoods/ensemble_likelihood/ensemble_likelihood.py
+++ b/imagine/likelihoods/ensemble_likelihood/ensemble_likelihood.py
@@ -46,10 +46,11 @@ class EnsembleLikelihood(Likelihood):
         obs_mean = observable.ensemble_mean().val.get_full_data()
 
         u_val = obs_val - obs_mean
-
+        self.logger.debug(('obs_val', obs_val))
+        self.logger.debug(('obs_mean', obs_mean))
         # compute quantities for OAS estimator
         mu = np.vdot(u_val, u_val)/k/n
-        self.logger.debug("mu: %f" % mu)
+        self.logger.debug(("mu", mu))
 
         alpha = (np.einsum(u_val, [0, 1], u_val, [2, 1])**2).sum()
         alpha /= k**2
diff --git a/imagine/observers/hammurapy/hammurapy.py b/imagine/observers/hammurapy/hammurapy.py
index 137939e7548d2c67552f8f5686b4eddb0bb0c03e..eb94573559215159ed3cdbf52b574ad96b4dd004 100644
--- a/imagine/observers/hammurapy/hammurapy.py
+++ b/imagine/observers/hammurapy/hammurapy.py
@@ -129,6 +129,8 @@ class Hammurapy(Observer):
         # know that the distribution strategy is the same for the
         # randam samples and the magnetic field itself
         random_seed = magnetic_field.random_seed.data[local_ensemble_index]
+        self.logger.debug("Local ensemble index %i uses random seed %i" %
+                          (local_ensemble_index, random_seed))
         custom_parameters += [['./MagneticField/Random', 'seed',
                                random_seed]]