The problem is with the -chemistry_dimensions specification in GRID. You should include XY instead of Z, -chemistry_dimensions XY Any dimension not included in -chemistry_dimensions (Z in this case), must have exactly two nodes. Flow and transport is always 3D, but the idea is to save some calculation time in the chemistry if there is symmetry in the problem. In your case, phast will calculate chemistry on the lower XY plane of nodes and copy the results to the upper plane of nodes, however, you must have complete symmetry in flow and chemistry to use this option. You could also avoid the error message by running with -chemistry_dimensions XYZ, and then chemistry will be run for all active nodes. In your input file, X and Z were not included in -chemistry_dimensions, so PHAST wanted there to be only 2 nodes in the X and Y directions. Sorry about the obscure error message. That part of the program must have been written early on, and the keywords and options have changed since that error message was coded. David 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 2255 since Jan 22, 1998.