diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index b1efe81..3daf871 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -446,7 +446,7 @@ verbose: {verbose}")
 
             imls = [dask.delayed(compute_IMT_exceedance)(rx_lat[i], rx_lon[i], distances[i].flatten(), fr, p, lambdas,
                                                          forecast_len, lambdas_perc, m_range, m_pdf, m_cdf, model,
-                                                         log_level=logging.DEBUG, imt=imt, IMT_min=0.0, IMT_max=2.0,
+                                                         log_level=logging.DEBUG, imt=imt, IMT_min=0.0, IMT_max=2.0, rtol=0.1,
                                                          rx_label=i) for i in iter]
             iml = dask.compute(*imls)
             iml_grid_raw.append(list(iml))