Note that in the code get_GDD.R you get an error here:
climdat.wstart <- merge(addstartdat, climdat, by.x="uniqueID", by.y="stn", all.y=TRUE) # currently not losing any data ...
Error in fix.by(by.y, y) : 'by' must specify a uniquely valid column
> climdat.wstart <- climdat.wstart[which(is.na(climdat.wstart$expstartdate)==FALSE),]
So I did not use GDD_plotting.R so much as I could not update the relevant csv files but see gddbyutah_fromGDD_plottingcode.pdf for what the output looked like (values seem pretty similar to my new estimates on QUICK glance)
I push this today, see GDD_plottinglimcues.R.
Note that in the code get_GDD.R you get an error here:
So I did not use GDD_plotting.R so much as I could not update the relevant csv files but see gddbyutah_fromGDD_plottingcode.pdf for what the output looked like (values seem pretty similar to my new estimates on QUICK glance)