R/alkIO.R

In pmcelhany/aws: Instant Ocean Chemistry Caclulator

#' Instant Ocean alkaliniy from salinity with only DI water addition at 400 uatm
#' pCO2
#' @description  this function is just a placeholder assuming a linear
#'   relationship with alk at 30 psu = 3100 and alk at 35 psu = 3200; need to
#'   develop an empircal curve in the range 25-60 psu
#' @param salinity Instant Ocean water salinity (psu)
#' @param slope empirically derived slope of Instant Ocean alklaliniy  vs
#'   salinity relationship. Default is 3200/35.
#' @param intercept empirically derived slope of Instant Ocean alklaliniy  vs
#'   salinity relationship. Default is 0.
#' @return alkalinity of Instant Ocean (umol/kg)
#' @export
alkIO <- function(salinity, slope = 3200/35, intercept = 0){
  alk <- salinity * slope + intercept
  return(alk)
}