fix 366 handling

ecohealthalliance · Nov 6, 2023 · 8a33183 · 8a33183
1 parent 80aa37d
commit 8a33183
Show file tree

Hide file tree

Showing 2 changed files with 146 additions and 89 deletions.
diff --git a/R/calculate_forecasts_anomalies.R b/R/calculate_forecasts_anomalies.R
@@ -73,8 +73,13 @@ calculate_forecasts_anomalies <- function(ecmwf_forecasts_transformed,
     lead_dates <- seq(lead_start_date, lead_end_date, by = "day")
     lead_doys <- yday(lead_dates)
     if(366 %in% lead_doys) {
-      lead_doys <- lead_doys[lead_doys!=366]
-      lead_doys <- c(lead_doys, tail(lead_doys, 1) + 1)
+      if(tail(lead_doys, 1) == 366){
+        lead_doys <- lead_doys[lead_doys!=366]
+        lead_doys <- c(lead_doys, 1)
+      }else{
+        lead_doys <- lead_doys[lead_doys!=366]
+        lead_doys <- c(lead_doys, tail(lead_doys, 1) + 1)
+      }
     }
 
     doy_start <- head(lead_doys, 1)