diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 0163270..e79289b 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -132,7 +132,7 @@ verbose: {verbose}")
     logger.debug(f"Rbeast version {version('rbeast')}")
     logger.debug(f"Dask version {version('dask')}")
 
-    dask.config.set(scheduler='threads')
+    dask.config.set(scheduler='threads', pool=ThreadPoolExecutor(16))
 
     # run magnitude distribution modeling if selected by user and no magnitude pdf file provided
     if m_select and m_range[0] == None and m_pdf[0] == None: