# HG changeset patch # User greg # Date 1539872191 14400 # Node ID fbe9aa2dc181ab0630774916f61d6beef385c8ba # Parent 35bf606a11dfa45e49aa038c61026c994092af84 Uploaded diff -r 35bf606a11df -r fbe9aa2dc181 insect_phenology_model.R --- a/insect_phenology_model.R Wed Aug 08 14:01:41 2018 -0400 +++ b/insect_phenology_model.R Thu Oct 18 10:16:31 2018 -0400 @@ -312,6 +312,13 @@ end_date_ytd_row = end_date_ytd_row[1]; # The end date is contained within the input_ytd data. end_doy_ytd = as.integer(temperature_data_frame$DOY[end_date_ytd_row]); + if (end_doy_ytd > end_date_ytd_row + 1) { + # The input year-to-date dataset is missing 1 or more + # days of data. + days_missing = end_doy_ytd - end_date_ytd_row; + msg = cat("The year-to-date dataset is missing ", days_missing, " days of data.\n"); + stop_err(msg); + } } else { end_date_ytd_row = 0; } @@ -335,6 +342,13 @@ # Save the first DOY to later check if start_date is Jan 1. start_doy_ytd = as.integer(temperature_data_frame$DOY[1]); end_doy_ytd = as.integer(temperature_data_frame$DOY[num_ytd_rows]); + if (end_doy_ytd > end_date_ytd_row + 1) { + # The input year-to-date dataset is missing 1 or more + # days of data. + days_missing = end_doy_ytd - end_date_ytd_row; + msg = cat("The year-to-date dataset is missing ", days_missing, " days of data.\n"); + stop_err(msg); + } } } else { # We're processing only the 30 year normals data, so create an empty