New functions `usgs_nearest_gage()` and `usgs_calib_rarray()` to auto populate gage calibration rendering process

[rburghol]

• Get closest DA gage match (requires dh_feature database with usgs_full_drainage features): usgs_nearest_gage(riverseg_feature, riverseg_json)
  • riverseg_feature - RomFeature
  • riverseg_json - te json model
• Create a properly formatted param list for gage_vs_model.Rmd: usgs_calib_rarray(riverseg_json, gage_info, run_id)

Example:

# get river model feature
riverseg ="JA5_7480_0001"
hydrocode=paste0('vahydrosw_wshed_',riverseg)
riverseg_feature <- RomFeature$new(ds,list(hydrocode=hydrocode),TRUE)
riverseg_json <- ds$get_json_prop(riverseg_model$pid)

# Call usgs_nearest_gage() to get best gage
gage = usgs_nearest_gage(riverseg_feature, riverseg_json)
gage_info = dataRetrieval::readNWISsite(gage$gageid)

# populate render params
render_params <- usgs_calib_rarray(riverseg_json, gage_info, 11)
rmarkdown::render(
  paste(github_location,'/hydro-tools/USGS/gage_vs_model.Rmd',sep="/"),
  output_file = paste0(
    '/WorkSpace/modeling/projects/james_river/appomattox/', 
    gage_info$site_no, '_',
    riverseg_json$riverseg$value
  ),
  params = render_params
)

@COBrogan