# HG changeset patch # User greg # Date 1539872222 14400 # Node ID b7c9837047d9843ec8e557873f48618a69d6e5d2 # Parent e2d12bafe9d74f2a23b8fe5215462fe52bb3affb Uploaded diff -r e2d12bafe9d7 -r b7c9837047d9 insect_phenology_model.R --- a/insect_phenology_model.R Wed Aug 08 14:16:58 2018 -0400 +++ b/insect_phenology_model.R Thu Oct 18 10:17:02 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