# This example equilibrates pure water with calcite and writes the
# dump results to file.
phrLoadDatabaseString(phreeqc.dat)
phrSetDumpFileOn(TRUE)
phrSetDumpFileName(file.path(tempdir(), "phreeqc.dump"))
input <- c(
'SOLUTION 1 Pure water ',
'EQUILIBRIUM_PHASES 1 ',
' Calcite 0 10 ',
'SAVE solution 1 ',
'SAVE equilibrium_phases 1 ',
'DUMP ',
' -solution 1 ',
' -equilibrium_phases 1 '
)
if (is.null(phrRunString(input))) {
cat(paste("see ", phrGetDumpFileName(), ".\n", sep = ""))
}
Run the code above in your browser using DataLab