Read data from an EQ6 output file

Share:

Description

Extract computational results for aqueous species, solid phases, mineral saturation states, or speciation summaries at each step of reaction progress in an EQ6 output file. The results are written to a comma-separated value file that can be read by other programs.

Usage

1
  eqdata(file, species, property = "log act", outfile = TRUE)

Arguments

file

character, path to EQ6 output file

species

character, name(s) of species or minerals

property

character, property to get

outfile

logical or character, file for saving results

Details

The first argument, file, is the name of the EQ6 (Wolery, 1992; Wolery and Daveler, 1992) output file. species indicates the aqueous species, solid phases, minerals, or basis species for which you want values; multiple names can be provided except for basis species, which can be a single value. property indicates the property to retrieve. Specifying a value other than one listed below will cause an error.

  • Aqueous species: conc, log conc, log g, or log act

  • Solid phases: log moles, moles, grams, or volume, cc

  • Minerals (saturation states): affinity, kcal

  • Basis species (speciation): molal conc or per cent

The result of the function is a data frame (returned invisibly), with columns zi (reaction progress), T (temperature in °C), aH2O (activity of water) and one column for each of the requested species or, for speciation of basis species, one column for each unique species found in all of the speciation summary blocks for that basis species. Values are listed as NA (not available) for species or phases that are not present in the EQ6 output at any of the increments of reaction progress.

If outfile is TRUE, the result is saved in a file named like file.property.csv, in the same directory as file. The name of the outfile can be provided to override this naming scheme, or this argument can be set to FALSE or NULL, to turn off writing the result to a file.

Thanks to Peter Canovas and Everett Shock for helping to test the code and offering ideas for improvements. The function has been tested with output files generated by EQ3/6 version 7.1 running on a Unix platform.

References

Wolery, T. J. (1992) EQ3/6, A Software Package for Geochemical Modeling of Aqueous Systems: Package Overview and Installation Guide (Version 7.0). Lawrence Livermore National Laboratory, UCRL-MA-110662 PT I. http://www.wipp.energy.gov/library/cra/2009_cra/references/Others/Wolery_1992_EQ36_A_Software_Package_for_Geochemical_Modeling_of_Aqueous_Systems_ERMS241375.pdf

Wolery, T. J. and Daveler, S. A. (1992) EQ6, A Computer Program for Reaction Path Modeling of Aqueous Geochemical Systems: Theoretical Manual, User's Guide, and Related Documentation (Version 7.0). Lawrence Livermore National Laboratory, UCRL-MA-110662 PT IV. http://www.wipp.energy.gov/library/cra/2009_cra/references/Others/Wolery_Daveler_1992_EQ36_A_Computer_Program_for_Reaction_Path_Modeling_of_Aqueous_Geochemical_Systems_ERMS241379.pdf

Examples

1
2
3
4
5
6
7
## Not run: 
# if an EQ6 output file named "rainbow2.6o" is in the current 
# working directory, the following command will output values 
# of log act (logarithm of activity) for the selected aqueous 
# species to a file named rainbow2.6o.log act.csv
eqdata("rainbow2.6o",c("h+","sio2,aq","h2,aq"),"log act") 
## End(Not run)

Want to suggest features or report bugs for rdrr.io? Use the GitHub issue tracker.