The only way that I can think of is to recalculate the pe in USER_PUNCH. The following is an example for iron. You would have to code an equation for each redox couple. David SOLUTION 1 temp 10 pH 5 pe 4 Fe(2) 1 Fe(3) 1 USER_PUNCH -heading pe_Fe(2)/Fe(3) 10 pe_fe = LA("Fe+3") - LA("Fe+2") - LK_SPECIES("Fe+3") 20 PUNCH pe_fe END David Parkhurst (dlpark@xxxxxxxx) U.S. Geological Survey Box 25046, MS 413 Denver Federal Center Denver, CO 80225 Project web page: https://wwwbrr.cr.usgs.gov/projects/GWC_coupled
Please note that some U.S. Geological Survey (USGS) information accessed through this page may be preliminary in nature and presented prior to final review and approval by the Director of the USGS. This information is provided with the understanding that it is not guaranteed to be correct or complete and conclusions drawn from such information are the sole responsibility of the user.
Any use of trade, product, or firm names in this publication is for descriptive purposes only and does not imply endorsement by the U.S. Government.
The URL of this page is:
Last modified: $Date: 2005-09-13 21:04:21 -0600 (Tue, 13 Sep 2005) $
Visitor number 1608 since Jan 22, 1998.