eqdata: Read data from an EQ6 output file

Description Usage Arguments Details References Examples

View source: R/eqdata.R


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.


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



character, path to EQ6 output file


character, name(s) of species or minerals


character, property to get


logical or character, file for saving results


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.

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.


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


## 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)

CHNOSZ documentation built on July 27, 2017, 3:01 a.m.