diff --git a/projects/exceeding_snow_loads/section_discussion/main_comparison_with_eurocode_global.py b/projects/exceeding_snow_loads/section_discussion/main_comparison_with_eurocode_global.py
index d71b0e53a02c08947332181f0503f86c0f70ca22..8ee2cf457880b7dbca69a5e42cc4c2af618031c7 100644
--- a/projects/exceeding_snow_loads/section_discussion/main_comparison_with_eurocode_global.py
+++ b/projects/exceeding_snow_loads/section_discussion/main_comparison_with_eurocode_global.py
@@ -16,7 +16,7 @@ def max_graph_annual_maxima_comparison():
     we also choose them because they belong to a different climatic area
     :return:
     """
-    save_to_file = False
+    save_to_file = True
     study_classes = [
                         CrocusSnowDensityAtMaxofSwe,
                         CrocusDifferenceSnowLoad,
@@ -30,10 +30,16 @@ def max_graph_annual_maxima_comparison():
     for study_class in study_classes:
         ylim, yticks = study_class_to_ylim_and_yticks[study_class]
 
+        # marker_altitude_massif_name = [
+        #                                   ('magenta', 900, 'Ubaye'),
+        #                                   ('darkmagenta', 1800, 'Vercors'),
+        #                                   ('mediumpurple', 2700, 'Beaufortain'),
+        #                               ][:]
+
         marker_altitude = [
-                              ('deepskyblue', 900),
-                              ('dodgerblue', 1800),
-                              ('blue', 2700),
+                              ('magenta', 900),
+                              ('darkmagenta', 1800),
+                              ('mediumpurple', 2700),
                           ][:]
         ax = plt.gca()
         for color, altitude in marker_altitude:
diff --git a/projects/exceeding_snow_loads/section_results/main_result_trends_and_return_levels.py b/projects/exceeding_snow_loads/section_results/main_result_trends_and_return_levels.py
index c77b83745d16967fe3cff614a0486451458410b8..14ca7aba93ae53fa8182868eb701ffc56105e859 100644
--- a/projects/exceeding_snow_loads/section_results/main_result_trends_and_return_levels.py
+++ b/projects/exceeding_snow_loads/section_results/main_result_trends_and_return_levels.py
@@ -73,9 +73,9 @@ def intermediate_result(altitudes, massif_names=None,
 
         # Plots
         # plot_trend_map(altitude_to_visualizer)
-        # plot_trend_curves(altitude_to_visualizer={a: v for a, v in altitude_to_visualizer.items() if a >= 900})
+        plot_trend_curves(altitude_to_visualizer={a: v for a, v in altitude_to_visualizer.items() if a >= 900})
         # plot_uncertainty_massifs(altitude_to_visualizer)
-        plot_uncertainty_histogram(altitude_to_visualizer)
+        # plot_uncertainty_histogram(altitude_to_visualizer)
         # plot_selection_curves(altitude_to_visualizer)
         # plot_intensity_against_gumbel_quantile_for_3_examples(altitude_to_visualizer)
 
@@ -89,11 +89,11 @@ def major_result():
                            ConfidenceIntervalMethodFromExtremes.ci_mle][1:]
     # massif_names = ['Beaufortain', 'Vercors']
     massif_names = None
-    study_classes = paper_study_classes[:]
+    study_classes = paper_study_classes[:1]
     # study_classes = [CrocusSnowLoad3Days, CrocusSnowLoad5Days, CrocusSnowLoad7Days][::-1]
-    altitudes = [300, 600, 900, 1800, 2700][:2]
-    altitudes = [300, 600, 900, 1200, 1500, 1800]
-    # altitudes = paper_altitudes
+    altitudes = [300, 600, 900, 1800, 2700][:]
+    altitudes = [300, 600, 900, 1200, 1500, 1800][2:]
+    altitudes = paper_altitudes
     # altitudes = [900, 1800, 270{{0][:1]
     for study_class in study_classes:
         print('new stuy class', study_class)
diff --git a/projects/exceeding_snow_loads/section_results/plot_trend_curves.py b/projects/exceeding_snow_loads/section_results/plot_trend_curves.py
index 01ec6441c0118d596b3c5c182bc992cc773dd3dd..bb302a29d58788fc9c098bdeb79148b09d5787b0 100644
--- a/projects/exceeding_snow_loads/section_results/plot_trend_curves.py
+++ b/projects/exceeding_snow_loads/section_results/plot_trend_curves.py
@@ -47,7 +47,7 @@ def plot_trend_curves(altitude_to_visualizer: Dict[int, StudyVisualizerForNonSta
     # parameters
     width = 150
     size = 20
-    legend_fontsize = 20
+    legend_fontsize = 35
     color = 'white'
     labelsize = 15
     linewidth = 3