Merge pull request #96 from CIAT-DAPA/develop

fix season check to do the addition of 29 day in case of leap year

src/aclimate_resampling/resampling.py

@@ -292,11 +292,12 @@ def forecast_station(self, station, prob, daily_data_root, output_root, year_for
 
       # Adjust the year if the forecast period is 'tri' if necessary
       for idx2, sublist2 in enumerate(seasons_1):
-         if season == sublist2: 
-            year_forecast = year_forecast+1
+        if any(s in sublist2 for s in season): 
+            year_forecast += 1
+            break
 
       # Check if year of forecast is a leap year for February
-      leap_forecast = ((year_forecast %400 == 0) and (year_forecast % 100 == 0)) or ((year_forecast%4==0) and (year_forecast%100!=0))
+      leap_forecast =  (year_forecast % 4 == 0) and (year_forecast % 100 != 0 or year_forecast % 400 == 0)
 
       # Filter the February data for leap years
       clim_feb = clim.loc[clim['month'] == 2]