Jump to content
Geochemist's Workbench Support Forum

EmilyP.

Members
  • Posts

    1
  • Joined

  • Last visited

Everything posted by EmilyP.

  1. Greetings, I am a GWB newbie, and I'm trying to model chemical changes in a basaltic rock when I run water through it at elevated temperatures and CO2(g) fugacities. I am trying to use the flush model, and possibly also a kinetic reaction path (I'm trying to plot changes over time, e.g. 1000 yrs, if I control the rate at which water is moving through the rock, e.g. 1 kg of water/year moving through ~ 1 kg of rock). I am repeatedly getting one of 2 error messages: either "Residual's too large, 621th interation" if I have something set for the charge balance (I have tried Cl-, Na+, both with a set value, based on N-MORB concentrations, or with no value put in - the only thing that changes is which iteration I get stuck on). or "N-R didn't converge after 999 its., maximum residual = 1.#R, Xi = 0.0000" if I turn off the charge balance completely. How can I fix my charge balance? Or, even, is my entire approach wrong? I'm attaching the script of the code I'm using for this model. Thank you for any help or advice! Sincerely, Emily P. EmilyP - GWB Script for Exp 1v3.txt
×
×
  • Create New...