test-levelPlot.R

In taywater/pwdgsi: Philadelphia Water Department Green Stormwater Infrastructure

library(pool)

conn_old <- dbPool(
  drv = RPostgres::Postgres(),
  host = "PWDMARSDBS1",
  port = 5434,
  dbname = "mars_data",
  user= Sys.getenv("shiny_uid"),
  password = Sys.getenv("shiny_pwd"),
  timezone = NULL)

rain <- old_rain(conn_old, "1267-2-1", "gage", "2024-03-27 15:15:00", "2024-03-28 04:30:00")
waterlevel_ft <- seq(.4, 3, by = 2.6/27) |> sort(decreasing = TRUE)

old_mars <- old_levelPlot(event = rain$gage_event_uid[[1]], 
                               structure_name = "CS2", 
                               storage_depth_ft = 6, 
                               obs_datetime = rain$dtime_est, 
                               obs_level_ft = waterlevel_ft,
                               #level_names = c("lvl1", "lvl2", "lvl3", "lvl4"),
                               orifice_show = FALSE,
                               orifice_height_ft = 1,
                               datetime_2 = rain$dtime_est[8:14],
                               level_ft_2 = waterlevel_ft[8:14],
                               datetime_3 = rain$dtime_est[15:21],
                               level_ft_3 = waterlevel_ft[15:21],
                               datetime_4 = rain$dtime_est[22:28],
                               level_ft_4 = waterlevel_ft[22:28],
                               metrics_show = FALSE)

poolClose(conn_old)

conn_sand <- dbPool(
  drv = RPostgres::Postgres(),
  host = "PWDMARSDBS1",
  port = 5434,
  dbname = "sandbox_dtime",
  user= Sys.getenv("shiny_uid"),
  password = Sys.getenv("shiny_pwd"),
  timezone = NULL)
 
rain_new <- marsFetchRainfallData(conn_sand, "1267-2-1", "gage", "2024-03-27 15:15:00", "2024-03-28 04:30:00") |>
  # Filter out the 00:00:00 value missing in old rain data for matching
  dplyr::filter(dtime != lubridate::ymd_hms("2024-03-28 00:00:00", tz = "America/New_York"))

new_mars <- marsWaterLevelPlot(event = rain_new$gage_event_uid[[1]], 
                              structure_name = "CS2", 
                              storage_depth_ft = 6, 
                              obs_datetime = rain_new$dtime, 
                              obs_level_ft = waterlevel_ft,
                              #level_names = c("lvl1", "lvl2", "lvl3", "lvl4"),
                              orifice_show = TRUE,
                              orifice_height_ft = 1,
                              datetime_2 = rain_new$dtime[8:14],
                              level_ft_2 = waterlevel_ft[8:14],
                              datetime_3 = rain_new$dtime[15:21],
                              level_ft_3 = waterlevel_ft[15:21],
                              datetime_4 = rain_new$dtime[22:28],
                              level_ft_4 = waterlevel_ft[22:28],
                              metrics_show = FALSE)

poolClose(conn_sand)