[CHNOSZ-commits] r553 - in pkg/CHNOSZ: . inst
noreply at r-forge.r-project.org
noreply at r-forge.r-project.org
Wed Jul 8 06:34:08 CEST 2020
Author: jedick
Date: 2020-07-08 06:34:08 +0200 (Wed, 08 Jul 2020)
New Revision: 553
Added:
pkg/CHNOSZ/inst/NEWS.Rd
Removed:
pkg/CHNOSZ/inst/NEWS
Modified:
pkg/CHNOSZ/DESCRIPTION
Log:
Convert NEWS to NEWS.Rd
Modified: pkg/CHNOSZ/DESCRIPTION
===================================================================
--- pkg/CHNOSZ/DESCRIPTION 2020-07-07 08:53:15 UTC (rev 552)
+++ pkg/CHNOSZ/DESCRIPTION 2020-07-08 04:34:08 UTC (rev 553)
@@ -1,6 +1,6 @@
-Date: 2020-07-07
+Date: 2020-07-08
Package: CHNOSZ
-Version: 1.3.6-26
+Version: 1.3.6-27
Title: Thermodynamic Calculations and Diagrams for Geochemistry
Authors at R: c(
person("Jeffrey", "Dick", , "j3ffdick at gmail.com", role = c("aut", "cre"),
Deleted: pkg/CHNOSZ/inst/NEWS
===================================================================
--- pkg/CHNOSZ/inst/NEWS 2020-07-07 08:53:15 UTC (rev 552)
+++ pkg/CHNOSZ/inst/NEWS 2020-07-08 04:34:08 UTC (rev 553)
@@ -1,1928 +0,0 @@
-CHANGES IN CHNOSZ 1.3.6-26 (2020-07-07)
----------------------------------------
-
-MAJOR CHANGES
-
-- The naming convention for methane has changed. The name "methane"
- now applies exclusively to the gas; "CH4", as before, represents the
- aqueous species. This behavior is consistent with inorganic gases
- but differs from most organic substances, where the English names
- refer to the aqueous species. This change was made because in the past
- there was no way to use info() or subcrt() with a single character
- argument to identify gaseous methane, which is common in geochemistry.
- (Note that info("CH4", "gas") can still be used.) A hack to revert to
- the old behavior is mod.OBIGT(info("CH4"), name = "methane").
-
-- The all-uppercase abbreviation OBIGT is now used for all file,
- function, and object names. In particular, the active database now
- sits at thermo()$OBIGT, and functions add.OBIGT() and mod.OBIGT()
- replace the previous add.obigt() and mod.obigt().
-
-CHANGES TO OBIGT DATABASE
-
-- Update species parameters for the 2019 Deep Earth Water (DEW) model.
-
-- OBIGT/SUPCRT92.csv: Change lowercase to uppercase "A" (Angstrom) in
- names of minerals with suffixes indicating interlayer spacing
- (amesite, clinochlore, cronstedtite, daphnite).
-
-- OBIGT/inorganic_aq.csv: Update arsenopyrite with data from Naumov et
- al., 1974 (via Unitherm).
-
-- OBIGT/inorganic_cr.csv: Correct syntax for number of H2O in formulas
- of some As-bearing minerals.
-
-- Install bibtex file (docs/OBIGT.bib) for use in logKcalc package.
-
-- Use uncompressed data files (i.e. .csv, not .csv.xz) and rename
- DEW_aq.csv to DEW.csv.
-
-- Remove old data file suffixes in references (e.g. [S92] was used to
- indicate that the data first appeared in SUPCRT92).
-
-- Add GHS for almandine, dickite, fluorphlogopite, glaucophane,
- grunerite, larnite, pyrope (SUPCRT92.csv) and bromellite
- (inorganic_cr.csv) from Robie and Hemingway, 1995.
-
-DEMOS AND VIGNETTES
-
-- Add demo/comproportionation.R: Gibbs energy of sulfur
- comproportionation, after Amend et al., 2020.
-
-- Remove vignette hotspring.Rnw.
-
-- Revise OBIGT.Rmd to reduce the size of the HTML file and make
- deep linking to individual sections work.
-
-- Work in progress: convert equilibrium.Rnw to equilibrium.Rmd.
-
-CHANGES TO FUNCTIONS
-
-- subcrt(): change 'action.unbalanced' argument to 'autobalance',
- which provides the ability to prevent autobalancing.
-
-- Changing the water model with water() now updates the literature
- references in thermo()$OBIGT.
-
-- thermo.refs() now shows CHNOSZ version and date.
-
-- subcrt() uses degree symbol in messages (°C).
-
-- Change thermo$... to thermo()$... in messages and comments.
-
-BUG FIXES
-
-- Fix palply.Rd for new warning "Non-file package-anchored link(s) in
- documentation object" in R CMD check.
-
-TODO
-
-- TODO: fix error in expr.species("AsO4---")
-
-- TODO: make a site for community uploaded data
-
-- TODO: refresh extdata/supcrt/compare.R
-
-- TODO: get arsenopyrite from Perfetti et al., 2008
-
-- TODO: make separate directories for default and optional data files
-
-- TODO: add missing G and H for dickite and other minerals in SUPCRT92.csv
-
-- TODO: demo/saturation.R: put Al+3 first in basis definition (it's
- confusing to have H2O first when the reacitons are balanced on Al)
-
-- TODO: diagram(): make error or warning for NA values of affinity
-
-- TODO: use markup in NEWS
-
-CHANGES IN CHNOSZ 1.3.6 (2020-03-16)
-------------------------------------
-
-- Fix for R-devel (4.0.0): test-add.protein.R had a test that failed
- with new stringsAsFactors = FALSE default for read.table().
-
-CHANGES IN CHNOSZ 1.3.5 (2020-02-16)
-------------------------------------
-
-- Fix for R-devel (4.0.0): remove documentation for unused argument
- '...' in add.protein.Rd.
-
-CHANGES IN CHNOSZ 1.3.4 (2019-12-15)
-------------------------------------
-
-THERMODYNAMIC DATA
-
-- Change data for dawsonite, scheelite, and ferberite from calories to
- Joules (as originally reported).
-
-- Use correct Cp coefficients for ferberite (2nd term in equation of
- Wood and Samson, 2000 is off by a factor of 10). Thanks to Xiangchong
- Liu and David Polya for the bug report and assistance.
-
-- Move more carbonates from SUPCRT92 back into the default database:
- artinite, azurite.
-
-DIAGRAMS
-
-- Improve the handling of solids in equilibrate(). They are now excluded
- from the equilibrium calculation, but their stability fields are
- calculated using the maximum affinity method (with diagram()). Where
- solids are stable, logarithms of activities of aqueous species are set
- to -999, and vice versa. Thanks to Feng Lai for the request and test
- case.
-
-PACKAGE CLEANUP
-
-- Remove files Sce.csv.xz and yeastgfp.csv.xz and functions yeastgfp()
- and yeast.aa(). These have been moved to the JMDplots package
- (https://github.com/jedick/JMDplots).
-
-- Remove protein_refseq.csv.xz (based on RefSeq release 61). An updated
- version, based on RefSeq release 95, is available in JMDplots.
-
-BUG FIXES
-
-- describe.reaction(): revert the Unicode double arrow. Unlike that, a
- simple equals sign is displayed correctly on all systems.
-
-- Fix bug in unexported function obigt2eos(): lambda (exponent in heat
- capacity equation) was incorrectly going through a units conversion.
-
-- subcrt(): properties of minerals are now output at the listed upper
- T limit (or transition temperature, for e.g. cr1 -> cr2). Previously,
- properties were set to NA at, and not only above, the T limit. Thanks
- to Evgeniy Bastrakov for the bug report.
-
-- In mosaic(), set the activities of only aqueous basis species to the
- total activity, which is taken from the incoming basis() definition.
- This fixes a bug where activities of minerals (particularly, native
- sulfur) were unexpectedly changed; their logarithms are now set to 0.
- Thanks to Evgeniy Bastrakov for the bug report.
-
-- Fix a bug where berman() returned NA values for K-feldspar below
- 298.15 K, due to improperly initialized values for the disorder
- properties. Thanks to Kaustubh Hakim for the bug report.
-
-- Use inherits() instead of class() for checking try() errors (don't
- use conditions of length > 1 in if() statements -- fix for R-devel).
-
-- Fix a bug in affinity() with calculations in IS-T (but not T-IS)
- space. This involves a workaround for the limitation in subcrt() that
- only IS (if present) can be the 'grid' variable.
-
-CHANGES IN CHNOSZ 1.3.3 (2019-08-02)
-------------------------------------
-
-NEW FEATURE: ENERGY UNITS
-
-- The csv data files in OBIGT now have a column named "E_units", which
- is used to specify the energy units for each species (cal or J).
-
-- The messages produced by info() (and checkGHS() and checkEOS()) use
- cal or J depending on the units in the database.
-
-- All properties returned by info() reflect the units in the database
- for each species (no conversion is done).
-
-- All properties returned by subcrt() continue to follow the user's
- units setting (from E.units(), default cal). The messages from
- subcrt() now include the units used for the output.
-
-- For backward compatibility, add.obigt() allows reading data from a
- user-supplied file with no "E_units" column. In this case, the units
- are set to cal.
-
-NEW FEATURE: SOLUBILITY UNITS
-
-- convert() now handles the output from solubility(), for converting
- calculated concentrations to ppm or ppb.
-
-- diagram() has a new argument 'levels' for specifying the levels at
- which to draw solubility contours.
-
-- demo/contour.R has been updated to show the 1, 10, 100, 1000 ppb
- contours for Au solubility.
-
-- Add demo/sphalerite.R for calculations of solubility of sphalerite,
- converted to log ppm (after Akinfiev and Tagirov, 2014).
-
-NEW FEATURES: DIAGRAMS
-
-- diagram(): add 'lty.cr' argument to specify line type for boundaries
- between minerals (two crystalline phases). Setting lty.cr = 0
- suppresses the lines, which is useful for plotting only the activities
- of aqueous species (can be used to make solubility diagrams for some
- systems).
-
-- Also add 'lty.aq' to specify the line type between aqueous species.
-
-- The 'names' argument of diagram() recognizes more values, making it
- easier to use in scripts. A value of NULL or TRUE plots the names, and
- NA or FALSE hides the names.
-
-- Add functions lNaCl(), lS(), lT(), lP(), lTP(), and lex() for making
- legend text to describe system conditions.
-
-BUG FIXES
-
-- In mosaic(), constant activities were used for each changing basis
- species, leading to an artificial inflation of the affinities of
- species around the transitions where the basis species become equally
- abundant. Now, affinities of formation of species are calculated using
- activities of basis species that are adjusted (i.e. lowered) according
- to their relative abundances.
-
-- The 'mixing' argument of mosaic() has been removed; it is no longer
- needed with the improved handling of activities of basis species.
-
-- Using the "Alberty" method, nonideal() calculated incorrect values of
- adjusted H, S, and Cp at given ionic strength. The correct equations
- are now used (taken from Alberty, 2001), and checks of the self-
- consistency of the calculated values and comparisons with the
- Debye-Hückel limiting slopes for G, H, and Cp have been added to
- test-nonideal.R.
-
-- For improved calculation of the Debye-Hückel limiting slopes, the
- equation for the Debye-Hückel constant was changed to use the equation
- given by Clarke and Glew, 1980 (valid to 150 °C).
-
-IMPROVEMENTS TO MOSAIC
-
-- In mosaic(), groups of changing basis species are equilibrated for a
- total activity given by the activity of the starting basis species in
- the incoming basis() definition. This provides a means to specify the
- the amount of a conserved element in the basis species.
-
-- equilibrate(): make it possible to process the output of mosaic().
- This feature combines the equilibrium activities of the formed species
- with those of the changing basis species to make an object that can be
- plotted with diagram(). This can be useful for conserving two elements
- in a system: one in the basis species, the other in the formed
- species. Thanks to Kirt Robinson for the feature request and test
- system for acetamide, which has been added to ?equilibrate.
-
-- Add tests to test-mosaic.R to check that activitives produced by
- mosaic() - equilibrate() and mosaic() - solubility() really are
- equilibrium activities, i.e. that the affinities of reactions between
- species are zero everywhere, in particular near the transitions of
- basis species.
-
-- In mosaic(), include the relative abundances in each group of basis
- species in the output ('E.bases').
-
-THERMODYNAMIC DATA
-
-- Correct the value of a1 for acetamide (Shock, 1993). Thanks to Kirt
- Robinson.
-
-- Add properties of aqueous species (dimethylamine, trimethylamine,
- resorcinol, phloroglucinol, cyclohexane carboxylate, and cyclohexane
- carboxylic acid) from LaRowe and Amend, 2019 and HKF parameters for
- dimethylammonium chloride from Shvedov and Tremaine, 1997.
-
-- Add data for zoisite, clinozoisite, and epidote from Gottschalk, 2004.
- Thanks to Evgeniy Bastrakov for contributing the data file.
-
-- Move some carbonates from SUPCRT92 back into default database:
- huntite, hydromagnesite, malachite, nesquehonite.
-
-- Add trisulfur radical ion (S3-) from Pokrovski and Dubessy, 2015.
-
-- Correct entries for aqueous phenanthrene and methylphenanthrene
- isomers to use Joules (Dick et al., 2013).
-
-DOCUMENTATION
-
-- demo/aluminum.R: add calculations using Si(OH)4 from the Akinfiev-
- Diamond model (SiO2 in these reactions is replaced by Si(OH)4 - 2 H2O).
-
-- Modify demo/saturation.R to use the argument recall feature of
- affinity().
-
-- In anintro.Rmd, include an example of using mod.obigt() to add a
- solid (magnesiochromite, with data from Klemme et al., 2000). Thanks
- to Alan Boudreau for the suggestion.
-
-- In anintro.Rmd, add description of add.obigt("AkDi") and change
- mod.obigt() example for aqueous species to use CoCl4-2 (Liu et al.,
- 2011).
-
-OTHER CHANGES
-
-- Add function moles(), to calculate the total number of moles of
- elements in the output from equilibrate().
-
-- subcrt(): stop with an error if length(logact) is not the same as the
- number of species.
-
-- Add 'transect' argument to affinity(), to allow forcing transect mode
- for any number of values of the variables.
-
-- solubility() has a new argument 'codeanal' that when set to TRUE
- indicates the printing of additional diagnostic messages.
-
-- check.obigt(): add 'E_units' column and fix a bug, also involving
- makeup() and entropy(), where many species with differences between
- database and computed values of Gibbs energy were not listed.
-
-CHANGES IN CHNOSZ 1.3.2 (2019-04-20)
-------------------------------------
-
-NEW FEATURES
-
-- Add thermo/stoich.csv.xz (loaded as thermo()$stoich), containing a
- precalculated stochiometric matrix for the default database, to speed
- up retrieve().
-
-- retrieve() now updates the stoichiometric matrix when the database
- changes, instead of performing a full recalculation.
-
-- Add 'ligands' argument to retrieve(), for getting metal-bearing
- species with a range of possible elements in the ligands.
-
-- diagram(): Starting with R version 3.6.0, 'fill' can be the name of
- any available HCL-based color palette (see hcl.pals() for a list).
-
-- Add demo/contour.R for gold solubility contours on logfO2-pH diagram.
-
-- In the output of affinity(), the list of 'vals' (the values of the
- input variables) now has names corresponding to the variables.
-
-- In diagram(), a background grid corresponding to major and/or minor
- tick marks can be added using the 'grid' argument. Or, use
- thermo.axis() after diagram() to overlay a grid.
-
-- Add demo/arsenic.R (Eh-pH diagram after Lu and Zhu, 2011).
-
-- Rename demo/go-IU.R to aluminum.R, and add plots from Tutolo et al.,
- 2014.
-
-- Add argument 'plot.box' (default TRUE) to thermo.plot.new() and
- 'plot.line' (default FALSE) to thermo.axis() to prevent unnecessary
- plotting of axis lines. thermo.axis() also avoids plotting tick
- marks at the plot limits.
-
-THERMODYNAMIC DATA
-
-- OBIGT: use dawsonite Cp coefficients from Tutolo et al., 2014.
-
-- OBIGT: Move H4SiO4 (Stefánsson, 2001) to AS04.csv (optional data).
-
-- Berman: rename .csv files to include whole year. berman()
- now uses the years in the filenames to assemble the data in reverse
- chronological order, and keeps only the most recent data for any
- duplicated minerals.
-
-- Berman: Add data for sudoite, daphnite, Mg-amesite, and Fe-amesite
- (Vidal et al., 1992, 2001, 2005).
-
-- Berman: make "albite" identical to "albite,low", not "albite,high".
-
-- Move transuranic actinide complexes to optional data (SLOP98), as
- entropies of the elements at 298.15 K are not available to check the
- self-consistency of the GHS values for the complexes.
-
-- obigt.Rmd: Use Zenodo DOI for citations of slop98.dat, slop07.dat, and
- slop16.dat.
-
-CHANGES IN CHNOSZ 1.3.1 (2019-03-02)
-------------------------------------
-
-- Fix unexported function obigt2eos(), which was affected by "Assignment
- to empty subset of data.frame fails check of length of value"
- (https://bugs.r-project.org/bugzilla/show_bug.cgi?id=17483). The bug
- caused errors in R-release (3.5.2) but not R-devel, when using info()
- for species with the AkDi model.
-
-- Fixes to maintain compatibility with R version 3.1.0: don't use
- lengths(), and provide explicit colClasses to read.csv() to avoid
- automatic assignment of character for long numeric values.
-
-CHANGES IN CHNOSZ 1.3.0 (2019-02-26)
-------------------------------------
-
-CRAN COMPLIANCE
-
-- Because data(thermo) and data(OBIGT) did not actually create said
- objects in the user's workspace (the expected behavior for data()),
- the package's data are now loaded by two new functions, reset() and
- obigt(). reset() is automatically run upon loading the package, so
- that existing scripts beginning with data(thermo) still work (this
- command now has no effect other than producing a warning).
-
-NEW FEATURES
-
-- Add thermo() as a convenience function to access or modify the
- package's data, especially various computational options.
-
-- Add retrieve() to retrieve all the species having given elements, or
- all species in a given chemical system. Thanks to Evgeniy Bastrakov
- for the suggestion.
-
-- Add AkDi() to calculate thermodynamic properties of aqueous
- nonelectrolytes using the Akinfiev-Diamond model. Thanks to Evgeniy
- Bastrakov for guidance.
-
-THERMODYNAMIC DATA
-
-- Revert to using SiO2(aq) from SUPCRT92 (i.e. Shock et al., 1989) in
- the default database.
-
-- Move SiO2(aq) from Apps and Spycher, 2004 and recalculated HSiO3- to
- new optional data file, OBIGT/AS04.csv.
-
-- Move H2AsO3- from OBIGT/SLOP98.csv to OBIGT/SUPCRT92.csv.
-
-- Add OBIGT/AkDi.csv as optional data file for aqueous species from
- Akinfiev and Diamond, 2003 and Akinfiev and Plyasunov, 2014.
-
-- Add gaseous HF and HCl to as requirements for the Akinfiev-Diamond
- model for dissolved gases.
-
-- Move Au+3, AuCl3-2, and AuCl4-3 back into default database.
-
-- Add/replace data for Ag+ and Cu+ and their complexes from Akinfiev
- and Zotov, 2001 and 2010.
-
-- Update HCl(aq) with data from Tagirov et al., 1997.
-
-- Update Pd+2 and complexes with data from Tagirov et al., 2013.
-
-- Update Pt+2 and complexes with data from Tagirov et al., 2015.
-
-- Update Zn+2 and complexes with data from Akinfiev and Tagirov, 2015.
-
-DOCUMENTATION
-
-- In demo/NaCl.R, indicate region not considered by Shock et al., 1992
- in developing the "g function" applicable to electrolytes in the
- revised HKF model, and note presence of resulting discontinuities (see
- man/examples.Rd).
-
-- Add demo/AkDi.R showing calculations of Henry's constant of dissolved
- gases, after Figure 1 of Akinfiev and Diamond, 2003.
-
-OTHER CHANGES
-
-- In describe.reaction(), change equals sign to reaction double arrow.
-
-- water() and water.SUPCRT92(): add 'P1' argument to choose whether to
- output 1 bar for Psat at temperatures less than 100 degrees C
- (default is TRUE).
-
-- diagram(): 'xlim' and 'ylim' arguments now apply to 2-variable
- diagrams. Thanks to Evgeniy Bastrakov for the suggestion.
-
-CHANGES IN CHNOSZ 1.2.0 (2019-02-09)
-------------------------------------
-
-BUG FIXES
-
-- Fix bug in nonideal() where "Zn" in formula was identified as charge.
- Thanks to Feng Lai for the bug report.
-
-- Fix a bug where subcrt()$reaction$coeffs was incorrect for reactions
- involving minerals with phase transitions. Also ensure that the output
- reaction stoichiometry is correct for duplicated species in reactions.
- Thanks to Grayson Boyer for the bug report.
-
-- For species in the revised HKF model, subcrt() now sets properties to
- NA where the density of H2O is less than 0.35 g/cm3, avoiding the
- output of bogus values in this region. Thanks to Evgeniy Bastrakov.
-
-- For systems where no basis species is present in all formation
- reactions, and the user hasn't provided balance coefficients, stop
- with an error instead of setting the balance cofficients to 1.
- Thanks to Shuang Kong for an example calculation and Tucker Ely for
- a previous suggestion to produce an error here. The affected code is
- in balance(), an unexported function used in equilibrate() and
- diagram() (and now also solubility()).
-
-NEW FEATURE: SOLUBILITY CALCULATIONS
-
-- Add solubility(). Run this after affinity() to calculate the
- solubility of a solid or gas defined as the conserved basis species,
- which is involved in the formation of one or more dissolved species.
-
-- Features include automatic detection of dissociation reactions and
- finding the final ionic strength for dissolution of a mineral into
- pure water (find.IS argument).
-
-- find.IS depends on the new argument recall feature of affinity() (or
- mosaic() if that is used instead). This allows a calculation to be
- re-run with the same settings except for particular additions or
- modifications, in this case modified values of ionic strength.
-
-- Revise demo/solubility.R to show solubility calculations for CO2(gas)
- and calcite as a function of T and pH.
-
-NEW FEATURE: EXPANDED ACTIVITY COEFFICIENT CALCULATIONS
-
-- Reorganize and expand options for activity coefficient calculations
- (set in thermo$opt$nonideal: Bdot, Bdot0, bgamma, bgamma0, or Alberty).
- The previous default, which corresponds to 'bgamma' (T- and
- P-dependent extended term parameter with single ion-size parameter),
- has been replaced by 'Bdot' (T-dependent extended term parameter and
- species-dependent ion-size parameter; see below).
-
-- nonideal() with the 'Bdot' or 'Bdot0' equation uses specific
- ion-size parameters for different ions, in accord with the HCh package
- (Shvarov and Bastrakov, 1999). Parameters are from Table 2.7 of
- Garrels and Christ, 1965.
-
-- nonideal() now calulates activity coefficients of neutral species,
- using the Setchénow equation. Whether the extended-term parameter in
- this equation is taken to be zero or is taken from the value for
- charged species (see above) is controlled by setting
- 'thermo$opt$Setchenow' to bgamma0 (default) or bgamma.
-
-- Add NaCl(), implementing a first-order calculation of the speciation
- of NaCl in water, taking account of activity coefficients and the
- reaction Na+ + Cl- = NaCl(aq).
-
-OTHER NEW FEATURES
-
-- Add dumpdata() for returning/writing all packaged thermodynamic data
- (including default database and optional data files). The file is
- also available on the website (chnosz.net/download/alldata.csv).
-
-- mosaic() has been rewritten to handle more than two changing groups
- of basis species.
-
-- mosaic() gets a new argument 'mixing' (default TRUE), indicating
- whether the Gibbs energy of ideal mixing should be included in the
- calculations with blend = TRUE. As a test of this change, an example
- using mosaic() to calculate the pH-dependent solubility of calcite
- has been added to solubility.Rd.
-
-DOCUMENTATION
-
-- Add demo/gold.R for calculations of Au solubility in hydrothermal
- chloride and sulfide solutions (based on diagrams from Akinfiev and
- Zotov, 2001, Stefánsson and Seward, 2004, and Williams-Jones et al.,
- 2009). This depends on the revised nonideal() and new NaCl() functions
- described above.
-
-- anintro.Rmd: add section on solubility calculations. Add cuprite to
- mosaic diagram example, and note about implications of changing
- balance coefficients to 1. Add 'mklinks.sh' script for adding links
- to CHNOSZ and R documentation.
-
-- Add demo/bison.R (average oxidation state of carbon of metagenome-
- derived proteins in different microbial phyla at Bison Pool)
-
-- demo/saturation.R replaces the previous activity_ratios.R, and
- includes not only activity ratios on the axes but also saturation
- limits for calcite, dolomite, magnesite, and brucite in the
- H2O-CO2-CaO-MgO-SiO2 system.
-
-- Add demo/glycinate.R showing logK of complexation of glycinate with
- divalent and monovalent metals.
-
-- Remove vignette wjd.Rnw (Winding journey down in Gibbs energy).
-
-THERMODYNAMIC DATA
-
-- The Berman data (Berman, 1988 and later additions) have replaced the
- SUPCRT92 data (based on Helgeson et al., 1978) for most minerals in
- the default database (i.e. the one loaded by data(thermo)). Only
- native elements, sulfides, halides, sulfates, and selected carbonates
- and oxides from SUPCRT92 that do not duplicate any in the Berman
- dataset are loaded by default.
-
-- The superseded SUPCRT92 minerals (primarily the silicates and
- Al-bearing minerals, as well as calcite, dolomite, hematite, and
- magnetite) are in the file extdata/OBIGT/SUPCRT92.csv, which can be
- loaded with add.obigt("SUPCRT92"). This is used in some examples for
- comparing the datasets (see ?berman and demo("go-IU")), but should not
- be needed by most users.
-
-- In thermo$obigt, the minerals using the Berman equations are
- identified as those having NA values for all thermodynamic properties
- (columns 8-20). These minerals have the state of "cr" (not "cr_Berman"
- as before).
-
-- The heat capacity equation in berman() has been expanded to
- Cp = k0 + k1T-1/2 + k2T-2 + k3T-3 + k4T-1 + k5T + k6T2 (as used in
- winTWQ).
-
-- Multipliers on values in the Berman data files (extdata/Berman/*.csv)
- have been removed, except for 10^5, 10^5, 10^5, and 10^8 on v1 to v4,
- to be consistent with the winTWQ data format.
-
-- berman() now reads data from a user-supplied file specified in
- thermo$opt$Berman.
-
-- Add mineral data from Berman and Aranovich, 1996 to
- extdata/Berman/testing and an example of loading these data via
- add.obigt() and thermo$opt$Berman (see ?berman).
-
-- Correct swapped v1 <-> v3 and v2 <-> v4 in extdata/Berman/DS10.csv.
-
-- To reduce confusion that "d5" might be a coefficient in the same
- equation as d0 .. d4, rename "d5" (constant that scales disordering
- enthalpy to volume) to "Vad" (notation used in Theriak-Domino manual).
-
-- Change abbreviation of grossular to Grs.
-
-- Make names for dicarboxylic acid anions from Shock, 1995 (in
- organic_aq.csv) consistent with those used in later papers, e.g.
- H-succinate and succinate-2 for the -1 and -2 charged ions. Remove
- duplicated entries for H-succinate and succinate-2 from biotic_aq.csv.
-
-- Update G, H, and S of celestite using properties of the dissolution
- reaction given by Reardon and Armstrong, 1987.
-
-- Add scheelite (CaWO4) and ferberite (FeWO4) using data adopted by
- Wood and Samson, 2000.
-
-- Update aqueous Au species with data from Akinfiev and Zotov, 2001
- and 2010, and Pokrovski et al., 2014.
-
-- Move SUPCRTBL updates (As and Al minerals and aqueous species and
- SiO2(aq)) into default database.
-
-- Values for SiO2(aq) are taken from Apps and Spycher, 2004. Thanks
- to John Apps for providing the document.
-
-- Recalculate GHS for HSiO3- by adding difference from SiO2(aq)
- (Sverjensky et al., 1997) to updated values for SiO2(aq).
-
-- Add test-recalculate.R to check that recalculated values are correctly
- entered in OBIGT.
-
-- Move superseded data for aqueous Au, As, and Al species, SiO2 and
- HSiO3- to OBIGT/SLOP98.csv.
-
-- Add or update glycine, diglycine, and triglycine (zwitterions and
- ions), and diketopiperazine, [Gly] and [UPBB] groups from Kitadai,
- 2014.
-
-- Add data for metal-glycinate complexes from Azadi et al., 2019,
- superseding Shock and Koretsky, 1995.
-
-- For reproducing previous calculations, superseded data for [Gly] and
- [UPBB], as well as [Met] (earlier superseded by LaRowe and Dick,
- 2012), and metal-glycinate complexes have been moved to
- OBIGT/OldAA.csv.
-
-DIAGRAMS
-
-- Lines in 1-D diagram()s can optionally be drawn as splines using the
- method for splinefun() given in the 'spline.method' argument (the
- default of NULL means no splines).
-
-- Add 'srt' argument to diagram() (rotation of line labels).
-
-- Export thermo.axis(), as it is useful for adding major and minor tick
- marks after (above) other plot elements such as legends.
-
-- diagram(): rename 'what' argument to 'type'.
-
-- diagram(): add new type of diagram, 'saturation', which is used to
- plot saturation lines for minerals (where their affinity equals
- zero).
-
-- thermo.plot.new() and, by extension, diagram(), now "listen" to the
- setting of par("mar"). However, if a new plot is being made, or the
- device has the default setting of mar from par(), mar is (re)set to
- the default value in thermo.plot.new().
-
-- diagram(): numeric 'names' now indicates subset of species names to
- include as labels on the plot.
-
-- diagram(): new argument 'family' to select the font family for species
- names (labels on the plot). New arguments 'bold' and 'italic' apply to
- formatted names (e.g. chemical formulas), and 'font' applies to those
- without formatting.
-
-- 'lines' is added to the return list of diagram(), giving the
- coordinates of lines (field boundaries) on 2-D diagrams (these are
- taken from the output of contourLines()).
-
-COMPUTATIONAL OPTIONS
-
-- Add 'exceed.rhomin' argument to subcrt() and affinity() to enable
- output of properties for species in the revised HKF model below 0.35
- g/cm3.
-
-- In equilibrate(), accept a length > 1 'normalize' argument to
- normalize the chemical formulas of only the selected species.
-
-- Add thermo$opt$maxcores (default 2) to specify maximum number of
- cores for parallel calculations with palply().
-
-- Keywords in basis(): Change 'CHNOPS+' to use O2 instead of e-, and add
- 'CHNOPSe' and 'MgCHNOPSe' for sets of basis species that have e-.
-
-- Add 'keep.duplicates' argument to thermo.refs(). Set it to TRUE to
- output a single primary reference for each species, keeping any
- duplicated references (but not including any secondary references in
- thermo$obigt$ref2). Thanks to Evgeniy Bastrakov for the suggestion.
-
-USABILITY ENHANCEMENTS
-
-- To provide better diagnostics for other apps using CHNOSZ, warning
- messages produced by subcrt() are now available in the output of
- affinity(), under 'sout$warnings'.
-
-- Change internal variable names in subcrt() for better readability
- (sinfo -> ispecies, inpho -> iphases, sinph -> phasespecies).
-
-- info(): display mineral names instead of just "cr" in messages (e.g.
- for info("SiO2")).
-
-- info(): species with no matches in default database are also
- searched in optional data files.
-
-- Add C implementation of counting occurrences of all letters in a
- string (src/count_letters.c) to speed up operation of count.aa().
-
-- read.fasta(): add support for file connections created using
- archive::archive_read (https://github.com/jimhester/archive).
-
-- The arguments in expr.species() have been reorganized for more
- flexible and concise usage.
-
-- Reduce message clutter for diagrams showing water stability limits.
-
-- In subcrt(), improve error message for invalid property names.
-
-- In affinity(), make invalid variable names an error.
-
-- Add 'return.sout' argument to affinity(), to return just the values
- calculated with subcrt().
-
-CHANGES IN CHNOSZ 1.1.3 (2017-11-13)
-------------------------------------
-
-MAJOR CHANGES:
-
-- Add berman() function and extdata/Berman/*.csv files for calculating
- thermodynamic properties of minerals using equations of Berman, 1988.
-
-- Calculations related to Berman's (1988) Figs. 1 and 2 for the lambda
- transition of quartz are available in the new demo lambda.R.
-
-- Add functions implementing the Deep Earth Water (DEW) model
- (Sverjensky et al., 2014): water.DEW() and its supporting functions
- calculateDensity(), calculateGibbsofWater(), calculateEpsilon(),
- calculateQ().
-
-- The computational setting for water (thermo$opt$water) can now be set
- using water("DEW"), water("IAPWS"), etc.
-
-- Usage of the DEW model is shown in the new demo DEW.R. This demo also
- depends on the Berman equations (above) and, for the last diagram in
- the demo, the following *four* NEWS items:
-
-- In equilibrate(), it is now possible to combine affinity calculations
- with variable activity of the balancing basis species (loga.balance).
- For example, in the last plot of the DEW demo, the calculation
- involves simultaneously varying temperature and logfO2 (given as
- arguments to affinity()) as well as total concentration of carbon
- (given by the loga.balance argument in equilibrate()).
-
-COMPUTATIONAL IMPROVEMENTS:
-
-- The 'alpha' argument of diagram() can be set to 'balance' to scale
- the values by the balancing component. This is useful for making
- "percent carbon" plots for systems where the species have different
- carbon numbers.
-
-- nonideal() now has three methods for calculating activity
- coefficients: 'Helgeson' and 'Helgeson0', which utilize the
- Debye-Huckel equation with parameters as described by Helgeson, 1969,
- and 'Alberty' (the only method previously available). The 'Helgeson'
- method depends on the following *two* NEWS items (but the 'Helgeson0'
- method omits the B-dot term):
-
-- All three water options (SUPCRT92, IAPWS95, DEW) can be used to
- calculate 'A_DH' and 'B_DH', i.e. A and B coefficients in the extended
- Debye-Huckel equation as described by Helgeson, 1969.
-
-- Add Bdot() to calculate the "B-dot" (or b_gamma) extended term
[TRUNCATED]
To get the complete diff run:
svnlook diff /svnroot/chnosz -r 553
More information about the CHNOSZ-commits
mailing list