diff --git a/src/seismic_hazard_forecasting.py b/src/seismic_hazard_forecasting.py
index 100dc14..e09aa19 100644
--- a/src/seismic_hazard_forecasting.py
+++ b/src/seismic_hazard_forecasting.py
@@ -837,7 +837,7 @@ verbose: {verbose}")
 
             fig, ax = plt.subplots(figsize=((width + 40) / 100.0, (height + 20) / 100.0),
                                    dpi=100)  # Increase fig size for labels
-            ax.imshow(gradient, aspect='auto', cmap=cmap.reversed(),
+            ax.imshow(gradient, aspect='auto', cmap="YlOrRd_r",
                       extent=[0, 1, vmin, vmax])  # Note: extent order is different for vertical
             ax.set_xticks([])  # Remove x-ticks for vertical colorbar
             num_ticks = 11  # Show more ticks