diff --git a/examples/example_16/ulens_model_fit.py b/examples/example_16/ulens_model_fit.py
index 62fb991f..ffcb4348 100644
--- a/examples/example_16/ulens_model_fit.py
+++ b/examples/example_16/ulens_model_fit.py
@@ -2356,6 +2356,10 @@ def _parse_results_EMCEE(self):
         if self._yaml_results:
             self._print_yaml_best_model()
 
+        if self._shift_t_0 and self._yaml_results:
+            print("Plots shift_t_0 : {:}".format(self._shift_t_0_value),
+                  **self._yaml_kwargs)
+
     def _extract_posterior_samples_EMCEE(self):
         """
         set self._samples_flat and self._samples for EMCEE
@@ -2500,6 +2504,7 @@ def _shift_t_0_in_samples(self):
                 index = self._fit_parameters.index(name)
                 values = self._samples_flat[:, index]
                 mean = np.mean(values)
+                self._shift_t_0_value = int(mean)
                 try:
                     self._samples_flat[:, index] -= int(mean)
                     if 'trace' in self._plots: