diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index ad58232..95297ba 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -69,6 +69,7 @@ def main(catalog_file, mc_file, pdf_file, m_file, m_select, mag_label, mc, m_max
     from matplotlib.contour import ContourSet
     import xml.etree.ElementTree as ET
     import json
+    import multiprocessing as mp
 
     logger = getDefaultLogger('igfash')
 
@@ -448,9 +449,10 @@ verbose: {verbose}")
 
         start = timer()
 
-        use_pp = False
+        use_pp = True
 
         if use_pp:         # use dask parallel computing
+            mp.set_start_method("fork", force=True)
             pbar = ProgressBar()
             pbar.register()
             iter = indices