Mercurial > repos > greg > insect_phenology_model

comparison insect_phenology_model.R @ 53:b7c9837047d9 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded
author greg
date Thu, 18 Oct 2018 10:17:02 -0400
parents 3f827fe31756
children 5bb1d76c29ca
comparison
equal deleted inserted replaced
52:e2d12bafe9d7 53:b7c9837047d9
310 end_date_ytd_row = which(temperature_data_frame$DATE==end_date); 310 end_date_ytd_row = which(temperature_data_frame$DATE==end_date);
311 if (length(end_date_ytd_row) > 0) { 311 if (length(end_date_ytd_row) > 0) {
312 end_date_ytd_row = end_date_ytd_row[1]; 312 end_date_ytd_row = end_date_ytd_row[1];
313 # The end date is contained within the input_ytd data. 313 # The end date is contained within the input_ytd data.
314 end_doy_ytd = as.integer(temperature_data_frame$DOY[end_date_ytd_row]); 314 end_doy_ytd = as.integer(temperature_data_frame$DOY[end_date_ytd_row]);
315 if (end_doy_ytd > end_date_ytd_row + 1) {
316 # The input year-to-date dataset is missing 1 or more
317 # days of data.
318 days_missing = end_doy_ytd - end_date_ytd_row;
319 msg = cat("The year-to-date dataset is missing ", days_missing, " days of data.\n");
320 stop_err(msg);
321 }
315 } else { 322 } else {
316 end_date_ytd_row = 0; 323 end_date_ytd_row = 0;
317 } 324 }
318 } else { 325 } else {
319 # We're plotting an entire year. 326 # We're plotting an entire year.
333 # Properly set the end_date to be Dec 31 of the year. 340 # Properly set the end_date to be Dec 31 of the year.
334 end_date = paste(year, "12", "31", sep="-"); 341 end_date = paste(year, "12", "31", sep="-");
335 # Save the first DOY to later check if start_date is Jan 1. 342 # Save the first DOY to later check if start_date is Jan 1.
336 start_doy_ytd = as.integer(temperature_data_frame$DOY[1]); 343 start_doy_ytd = as.integer(temperature_data_frame$DOY[1]);
337 end_doy_ytd = as.integer(temperature_data_frame$DOY[num_ytd_rows]); 344 end_doy_ytd = as.integer(temperature_data_frame$DOY[num_ytd_rows]);
345 if (end_doy_ytd > end_date_ytd_row + 1) {
346 # The input year-to-date dataset is missing 1 or more
347 # days of data.
348 days_missing = end_doy_ytd - end_date_ytd_row;
349 msg = cat("The year-to-date dataset is missing ", days_missing, " days of data.\n");
350 stop_err(msg);
351 }
338 } 352 }
339 } else { 353 } else {
340 # We're processing only the 30 year normals data, so create an empty 354 # We're processing only the 30 year normals data, so create an empty
341 # data frame for containing temperature data after it is converted 355 # data frame for containing temperature data after it is converted
342 # from the 30 year normals format to the year-to-date format. 356 # from the 30 year normals format to the year-to-date format.