Last data update: 2014.03.03

R: K_HSO4
K_HSO4R Documentation

K_HSO4

Description

PUBLIC function: calculates the dissociation constant of HSO4

Usage

K_HSO4(S, t, p=0, khso4="dickson")

Arguments

S

salinity in practical salinity units (i.e. no unit)

t

temperature in degrees centigrade

p

gauge pressure (total pressure minus atmospheric pressure) in bars

khso4

"dickson" or "khoo": specifies the S, t, dependency to be used. Default is "dickson". (see section below for references)

Value

the dissociation constant of HSO4 in mol/kg-solution on the free proton pH scale

Author(s)

Andreas F. Hofmann (a.hofmann@nioo.knaw.nl)

References

"dickson": Dickson1990, DOE1994, Zeebe2001; "khoo": Khoo1977, Roy1993, Millero1995

Examples

K_HSO4(35, 15)
K_HSO4(35, 15, 10)
plot(K_HSO4(35, 1:25), xlab="temperature / degC")

Results


R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> library(AquaEnv)
Loading required package: minpack.lm
Loading required package: deSolve

Attaching package: 'deSolve'

The following object is masked from 'package:graphics':

    matplot

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/AquaEnv/K_HSO4.Rd_%03d_medium.png", width=480, height=480)
> ### Name: K_HSO4
> ### Title: K_HSO4
> ### Aliases: K_HSO4
> ### Keywords: misc
> 
> ### ** Examples
> 
> K_HSO4(35, 15)
[1] 0.1468661
attr(,"unit")
[1] "mol/kg-soln"
attr(,"pH scale")
[1] "free"
> K_HSO4(35, 15, 10)
[1] 0.147927
attr(,"unit")
[1] "mol/kg-soln"
attr(,"pH scale")
[1] "free"
> plot(K_HSO4(35, 1:25), xlab="temperature / degC")
> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>