time start = 0 days, end = 1e6 years
temperature = 200
decouple Acetic\_acid(aq)
decouple Methane(aq)
swap CO2(aq) for HCO3-
1 kg free H2O
total molality CO2(aq) = 1e-10
total molality Methane(aq) = 1e-10
total molality Na+ = .1
total molality Acetic\_acid(aq) = .1
total molality Cl- = .001
pH = 7
kinetic redox-1 \
rxn = "Acetate  + H+  -> Methane(aq)  + CO2(aq)" \
rate\_con = 1.073e-13 \
rate\_law = 'rate\_con * (molality("Acetate") \
+ molality("NaCH3COO(aq)") + molality("Na(CH3COO)2-") \
+ molality("Acetic\_acid(aq)"))'