R/shinypbsol.R

In bentrueman/shinypbcusol: Calculate Lead and Copper Solubility Using Shiny

#' Calculate equilibrium lead solubility using `pbcusol` and `shiny`
#'
#' @description This function is used in deploying the `shinypbcusol` `shiny` app.
#'
#' @param ... Argument passed to function. (Currently ignored.)
#'
#' @importFrom dplyr %>%
#' @importFrom rlang .data
#' @importFrom vroom vroom_write
#'
#' @export
shinypbsol <- function(...) {
  sliderwidth <- "700px"
  ui <- shiny::fluidPage(
    theme = bslib::bs_theme(version = 4, bootswatch = "flatly"),
    shiny::titlePanel("Calculate equilibrium lead solubility"),
    shiny::fluidRow(
      shiny::column(5,
        shiny::sliderInput("pH", label = "pH", min = 4, max = 11, value = 7.3, step = .1, width = sliderwidth),
        shiny::sliderInput(
          "alkalinity", label = shiny::HTML("Alkalinity (mg CaCO<sub>3</sub> L<sup>-1</sup>)"),
          min = 0, max = 150, value = 20, step = .1, width = sliderwidth
        ),
        shiny::sliderInput(
          "phosphate", label = shiny::HTML("Orthophosphate (mg P L<sup>-1</sup>)"),
          min = 0, max = 3, value = .16, step = .1, width = sliderwidth
        ),
        shiny::sliderInput(
          "chloride", label = shiny::HTML("Chloride (mg L<sup>-1</sup>)"),
          min = 0, max = 500, value = 10, step = 1, width = sliderwidth
        ),
        shiny::sliderInput(
          "sulfate", label = shiny::HTML("Sulfate (mg L<sup>-1</sup>)"),
          min = 0, max = 500, value = 10, step = 1, width = sliderwidth
        ),
        shiny::sliderInput(
          "calcium", label = shiny::HTML("Calcium (mg L<sup>-1</sup>)"),
          min = 0, max = 200, value = 10, step = 1, width = sliderwidth
        ),
        shiny::sliderInput(
          "magnesium", label = shiny::HTML("Magnesium (mg L<sup>-1</sup>)"),
          min = 0, max = 200, value = 10, step = 1, width = sliderwidth
        ),
        shiny::sliderInput(
          "humic",
          label = shiny::HTML(
            "Humic/fulvic substances<sup>1</sup> (mg DOC L<sup>-1</sup>)"
          ),
          min = 0, max = 5, value = 0, step = .1, width = sliderwidth
        ),
        shiny::sliderInput(
          "aluminum",
          label = shiny::HTML(
            "Aluminum<sup>2</sup> (mg L<sup>-1</sup>)"
          ),
          min = 0, max = 1, value = 0, step = .002, width = sliderwidth
        )
      ),
      shiny::column(3, align = "left",
        shiny::HTML("&ensp;&nbsp;&nbsp;Equilibrium lead solubility at 25&deg;C:"),
        shiny::tableOutput("pbsol"),
        shiny::plotOutput("plot", width = 450, height = 765),
        shiny::downloadButton("download")
      )
    ),
    shiny::tags$footer(shiny::HTML(
      "<sup>1</sup>Experimental: see doi.org/10.1039/D1EW00903F<br>
      <sup>2</sup>Experimental: see doi.org/10.1021/acsestwater.1c00320"
    ))
  )

  server <- function(input, output, session) {

    out_leadsol <- shiny::reactive(eq_sol_tbl(input))
    out_minteq.v4 <- shiny::reactive(eq_sol_tbl(input, "minteq.v4", dbase = phreeqc::minteq.v4.dat))

    table1 <- shiny::reactive(
      tibble::tibble(
        `Calculated DIC` = glue::glue(
          "{round(pbcusol::calculate_dic(input$pH, input$alkalinity), 1)} ppm C"
        ),
        LEADSOL = glue::glue("{round(out_leadsol())} ppb Pb"),
        minteq.v4 = glue::glue("{round(out_minteq.v4())} ppb Pb")
      )
    )

    output$plot <- shiny::renderPlot({
      list(
        "LEADSOL" = out_leadsol(),
        "minteq.v4" = out_minteq.v4()
      ) %>%
        plot_kinetics()
    })

    output$pbsol <- shiny::renderTable({table1()})

    output$download <- shiny::downloadHandler(
      filename = function() {
        "shinypbsol-output.tsv"
      },
      content = function(file) {
        out_list <- list(
          "LEADSOL" = out_leadsol(),
          "minteq.v4" = out_minteq.v4()
        )
        out_table <- plot_kinetics(out_list, return = "table") %>%
          dplyr::mutate(
            pH = input$pH,
            alkalinity_mg_caco3_l = input$alkalinity,
            phosphate_mg_p_l = input$phosphate,
            chloride_mg_l = input$chloride,
            sulfate_mg_so4_l = input$sulfate,
            calcium_mg_l = input$calcium,
            magnesium_mg_l = input$magnesium,
            humic_substances_mg_c_l = input$humic,
            aluminum_mg_l = input$aluminum
          )
        vroom::vroom_write(out_table, file)
        # readr::write_csv(out_table, file)
      }
    )

  }

  shiny::shinyApp(ui, server)
}