Click here to donate to keep PhreeqcUsers open

SOLUTION_MASTER_SPECIES Doc Doc 0 Doc 12SOLUTION_SPECIESDoc = Doc log_k 0 SOLUTION 0 pH ...temp ...Doc 1O(0) 0.2N(5) 0.1Na ......RATESDOC_oxidation 10 k_doc_o2 = 1.57e-920 k_doc_NO3 = 1.67e-11 40 mo2 = MOL("O2")41 mNO3 = TOT("N(5)")50 rate = k_doc_o2 * mo2/(1e-10 + mo2)60 rate = rate + k_doc_NO3 * mNO3/(1.55e-4 + mNO3)70 moles = rate * TIME80 SAVE molesKINETICSDOC_oxidation; -formula Doc -1 CH2O +1 Ntg +1 #here I do not know how to make Ntg production conditional on NO3 consumption

SOLUTION_MASTER_SPECIES Doc Doc 0 Doc 12SOLUTION_SPECIESDoc = Doc log_k 0 SOLUTION 0Doc 1O(0) 0.2N(5) 0.1RATESDOC_oxidation10 k_doc_o2 = 1.57e-920 k_doc_NO3 = 1.67e-11 40 mo2 = MOL("O2")41 mNO3 = TOT("N(5)")50 rate = k_doc_o2 * mo2/(1e-10 + mo2)60 rate = rate + k_doc_NO3 * mNO3/(1.55e-4 + mNO3) * 1e-6/(1e-6 + mO2)70 moles = rate * TIME80 SAVE molesKINETICSDOC_oxidation; -formula Doc -1 CH2O +1 -step 1 1e1 1e2 1e3 1e4 1e5 1e6 1e7INCREMENTAL_REACTIONSUSER_GRAPH 1 -headings time O2 NO3 N2 -axis_scale x_axis auto auto auto auto log -initial_solutions false -connect_simulations true -plot_concentration_vs x -start10 GRAPH_X TOTAL_TIME20 GRAPH_Y MOL("O2"), TOT("N(5)"), MOL("N2") -end -active true

Doc -1.25 N -1 Ntg +0.5 CH2O +1.25