# React script, saved Tue May 22 2012 by bfarrell data = thermo.dat verify time start = 0 day, end = 30 day temperature = 25 C decouple CH3COO- H2O = 1 free kg balance on Na+ Na+ = 50 mmolal Cl- = 50 mmolal CH3COO- = 1 umolal pH = 7 swap O2(g) for O2(aq) O2(g) = .2 fugacity swap CO2(g) for HCO3- CO2(g) = -3.5 log fugacity reactants times 10 fix pH fix fugacity of O2(g) fix fugacity of CO2(g) kinetic microbe-1 rxn = "CH3COO- + 2*O2(aq) -> 2*HCO3- + H+" biomass = .001 rate_con = 2e-6 KA = 0 KD = 2e-5 mpower(CH3COO-) = 1 mpowerD(CH3COO-) = 1 ATP_energy = -48 ATP_number = 3 growth_yield = 900 decay_con = 8.33e-6 order1 = .25 kinetic microbe-2 rxn = "CH3COO- + 2*O2(aq) -> 2*HCO3- + H+" biomass = .001 rate_con = 6e-7 KA = 0 KD = 1e-5 mpower(CH3COO-) = 1 mpowerD(CH3COO-) = 1 ATP_energy = -48 ATP_number = 3 growth_yield = 450 decay_con = 1.67e-6 order1 = .25 react 1 kg of H2O react 20 umol of CH3COO- react 20 umol of Na+ flush