diff --git a/shapepipe/modules/MCCD_package/mccd_plot_utilities.py b/shapepipe/modules/MCCD_package/mccd_plot_utilities.py
index 05b72546c..a50920b90 100644
--- a/shapepipe/modules/MCCD_package/mccd_plot_utilities.py
+++ b/shapepipe/modules/MCCD_package/mccd_plot_utilities.py
@@ -120,7 +120,7 @@ def plot_meanshapes(starcat_path, output_path, nb_pixel, w_log,
                     plot_histograms=True):
     r"""Plot meanshapes and histograms."""
     # READ FULL STARCAT
-    starcat = fits.open(starcat_path)
+    starcat = fits.open(starcat_path, memmap=False)
 
     auto_colorbar = False
     colorbar_ampl = 1.
@@ -794,7 +794,7 @@ def rho_stats(starcat_path, output_path, rho_def='HSC',
     depending on the desired definition to use for tho statistics.
     """
     # Read starcat
-    starcat = fits.open(starcat_path)
+    starcat = fits.open(starcat_path, memmap=False)
 
     rho_stats_fun = None
 
diff --git a/shapepipe/modules/MCCD_package/shapepipe_auxiliary_mccd.py b/shapepipe/modules/MCCD_package/shapepipe_auxiliary_mccd.py
index 32e95825c..1dfc343f9 100644
--- a/shapepipe/modules/MCCD_package/shapepipe_auxiliary_mccd.py
+++ b/shapepipe/modules/MCCD_package/shapepipe_auxiliary_mccd.py
@@ -202,7 +202,7 @@ def mccd_fit_pipeline(trainstar_path, file_number_string, mccd_parser,
     use_SNR_weight = mccd_parser.get_extra_kw('use_SNR_weight')
 
     # Open fits file
-    starcat = fits.open(trainstar_path)
+    starcat = fits.open(trainstar_path, memmap=False)
 
     mccd.auxiliary_fun.mccd_fit(
         starcat=starcat[1],
@@ -227,7 +227,7 @@ def mccd_validation_pipeline(teststar_path, mccd_model_path, mccd_parser,
     # Get MCCD parameters
     save_extension = '.fits'
     mccd_val_kw = mccd_parser.get_val_kw()
-    testcat = fits.open(teststar_path)
+    testcat = fits.open(teststar_path, memmap=False)
 
     # Check if there is the fitted model
     if os.path.isfile(mccd_model_path):
@@ -259,7 +259,7 @@ def mccd_interpolation_pipeline(mccd_model_path, galcat_path, pos_params,
     # Import MCCD model
     mccd_model = mccd.mccd_quickload(mccd_model_path)
     # Open galaxy catalog
-    galcat = fits.open(galcat_path)
+    galcat = fits.open(galcat_path, memmap=False)
 
     # Extract positions
     x_pos = galcat[2].data[pos_params[0]]
diff --git a/shapepipe/modules/mccd_merge_starcat_runner.py b/shapepipe/modules/mccd_merge_starcat_runner.py
index 9e8b5b911..36318e3ff 100644
--- a/shapepipe/modules/mccd_merge_starcat_runner.py
+++ b/shapepipe/modules/mccd_merge_starcat_runner.py
@@ -60,7 +60,7 @@ def mccd_merge_starcat_runner(input_file_list, run_dirs, file_number_string,
                 my_mask[i, j] = True
 
     for name in input_file_list:
-        starcat_j = fits.open(name[0])
+        starcat_j = fits.open(name[0], memmap=False)
 
         stars = np.copy(starcat_j[hdu_table].data['VIGNET_LIST'])
         stars[stars < -1e6] = 0