next up previous
Next: taxonomy Up: CHNOSZ examples Previous: util.seq

util.units

## Don't show: data(thermo)
thermo$obigt: 1809 aqueous, 3368 total species
## End(Don't show) ## examples using convert # temperature (Kelvin) to degrees C convert(273.15, "C")
[1] 0
# temperature (degrees C) to Kelvin convert(100, "K")
[1] 373.15
# Gibbs energy (cal mol-1) to/from logK convert(1000, "logK")
[1] -0.7330067
convert(1000, "logK", T=373.15)
[1] -0.5856785
convert(1, "G")
[1] -1364.244
# Eh (volt) to pe convert(-1, "pe")
[1] -16.90445
convert(-1, "pe", T=373.15)
[1] -13.5068
# logfO2 to E0 (volt) convert(-80, "E0")
info.character: found H2O(liq), also available in gas subcrt: 4 species at 298.15 K and 1 bar (wet) [1] -0.3681732
convert(-80, "E0", pH=5)
info.character: found H2O(liq), also available in gas subcrt: 4 species at 298.15 K and 1 bar (wet) [1] -0.2498612
convert(-80, "E0", pH=5, logaH2O=-5)
info.character: found H2O(liq), also available in gas subcrt: 4 species at 298.15 K and 1 bar (wet) [1] -0.1019711
# calorie to/from joule convert(10, "J")
[1] 41.84
convert(10, "cal")
[1] 2.390057
# cm3bar to calories convert(10, "calories")
[1] 0.2390057
## examples showing unit settings T.units("C")
[1] "C"
T1in <- envert(25, "C") # no convertion T1out <- outvert(313.15, "K") # K to C T.units("K")
[1] "K"
T2in <- envert(298.15, "C") # K to C T2out <- outvert(-233.15, "C") # C to K # these are the same temperature (25 deg C) stopifnot(all.equal(T1in, T2in)) # these are numerically equivalent (40 deg C / 40 K) stopifnot(all.equal(T1out, T2out)) T.units("C") # return to default
[1] "C"


next up previous
Next: taxonomy Up: CHNOSZ examples Previous: util.seq