[contrasting] add slight modification for the second round of review for NHESS

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:

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)

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