diff --git a/extreme_data/meteo_france_data/scm_models_data/abstract_study.py b/extreme_data/meteo_france_data/scm_models_data/abstract_study.py
index 548c81d730d2ddd042bcaae2dafb0d299c32eec8..083f3c9eb6b7e8f7816b177c9473894d6395156c 100644
--- a/extreme_data/meteo_france_data/scm_models_data/abstract_study.py
+++ b/extreme_data/meteo_france_data/scm_models_data/abstract_study.py
@@ -458,6 +458,17 @@ class AbstractStudy(object):
         data = dict(zip(columns, [longitude[self.flat_mask], latitude[self.flat_mask]]))
         return pd.DataFrame(data=data, index=self.study_massif_names, columns=columns)
 
+    @property
+    def _save_excel_with_longitutde_and_latitude(self):
+        any_ordered_dict = list(self.year_to_dataset_ordered_dict.values())[0]
+        print(any_ordered_dict.variables.keys())
+        longitude = np.array(any_ordered_dict.variables['LON'])[self.flat_mask]
+        latitude = np.array(any_ordered_dict.variables['LAT'])[self.flat_mask]
+        data = [longitude, latitude]
+        df = pd.DataFrame(data=data, index=['Longitude', 'Latitude'], columns=self.study_massif_names).transpose()
+        print(df.head())
+        df.to_csv('S2M_latitude_and_longitude_for_the_centroid_of_each_massif.csv')
+
     @property
     def missing_massif_name(self):
         return set(self.all_massif_names(self.reanalysis_path, self.dbf_filename)) - set(
diff --git a/extreme_data/meteo_france_data/scm_models_data/case_studies/ribatet/__init__.py b/extreme_data/meteo_france_data/scm_models_data/case_studies/ribatet/__init__.py
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/extreme_data/meteo_france_data/scm_models_data/case_studies/nico/nico2.py b/extreme_data/meteo_france_data/scm_models_data/case_studies/ribatet/day_for_the_maxima.py
similarity index 100%
rename from extreme_data/meteo_france_data/scm_models_data/case_studies/nico/nico2.py
rename to extreme_data/meteo_france_data/scm_models_data/case_studies/ribatet/day_for_the_maxima.py
diff --git a/extreme_data/meteo_france_data/scm_models_data/safran/safran.py b/extreme_data/meteo_france_data/scm_models_data/safran/safran.py
index df4fdc012970e100ae2e741154a2a875d321bfba..b10fa0b82d4ffaf46c81747481d00d7e59557f2f 100644
--- a/extreme_data/meteo_france_data/scm_models_data/safran/safran.py
+++ b/extreme_data/meteo_france_data/scm_models_data/safran/safran.py
@@ -207,13 +207,13 @@ class SafranTemperature(Safran):
 
 
 if __name__ == '__main__':
-    altitude = 600
+    # altitude = 1800
+    altitude = 900
     year_min = 1959
-    year_max = 1962
-    study = SafranDateFirstSnowfall(altitude=altitude, year_min=year_min, year_max=year_max)
+    year_max = 2019
+    study = SafranSnowfall(altitude=altitude, year_min=year_min, year_max=year_max)
     print(study.study_massif_names)
-    print(study.massif_name_to_annual_maxima)
-    print(study.year_to_daily_time_serie_array[1959].shape)
-    print(study.massif_name_to_daily_time_series['Vanoise'].shape)
-    # print(study.year_to_annual_maxima[1959])
-    # print(study.ordered_years)
+    # print(study.massif_name_to_annual_maxima)
+    # print(study.year_to_daily_time_serie_array[1959].shape)
+    # print(study.massif_name_to_daily_time_series['Vanoise'].shape)
+    study._save_excel_with_longitutde_and_latitude()