From a0e7ffff2bf4b3f857f18319ede91819629bbb74 Mon Sep 17 00:00:00 2001
From: Le Roux Erwan <erwan.le-roux@irstea.fr>
Date: Thu, 19 Nov 2020 15:02:45 +0100
Subject: [PATCH] add script to compute longitude and latitude for Ribatet.
---
.../scm_models_data/abstract_study.py | 11 +++++++++++
.../case_studies/ribatet/__init__.py | 0
.../nico2.py => ribatet/day_for_the_maxima.py} | 0
.../scm_models_data/safran/safran.py | 16 ++++++++--------
4 files changed, 19 insertions(+), 8 deletions(-)
create mode 100644 extreme_data/meteo_france_data/scm_models_data/case_studies/ribatet/__init__.py
rename extreme_data/meteo_france_data/scm_models_data/case_studies/{nico/nico2.py => ribatet/day_for_the_maxima.py} (100%)
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 548c81d7..083f3c9e 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 00000000..e69de29b
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 df4fdc01..b10fa0b8 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()
--
GitLab