fmegaro = fmeg fmegarm = fmeg fmegato = fmeg fmegatm = fmeg fmegh = fmeg totc = [0.045, 0.105, 0.104, 0.105, 0.074, 0.056] #mol total organic carbon totw = [0.013588, 0.005854, 0.011788, 0.006379, 0.010690, 0.006620] #kg water tots = [1.41, 9.15, 3.21, 8.62, 4.31 , 8.38] #g soil dry mass ph = [5.02, 4.84, 5.21, 4.54, 5.23, 4.95] cace = [6.370, 2.799, 0.057, 2.668, 1.026, 1.839] #mM acetate cfe2 = [0.789, 22.234, 0.102, 22.973, 15.665, 7.177] #mM Fe hdsp = [0.042528, 0.042528, 0.044005, 0.044005, 0.043575, 0.043575] #head space volume L index = 0 writephrq('co', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fbuf, fgluco, fsom1co, fsom2co, fsom3, fferbco, fmegaco, fmegh, glu_stoich, o2scalar) writephrq('co', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fbuf, fgluco, fsom1co, fsom2co, fsom3, fferbco, fmegaco, fmegh, glu_stoich, o2scalar) writephrq('co', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fbuf, fgluco, fsom1co, fsom2co, fsom3, fferbco, fmegaco, fmegh, glu_stoich, o2scalar) index = 1 writephrq('cm', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3cm, fbuf, fglucm, fsom1cm, fsom2cm, fsom3, fferbcm, fmegacm, fmegh, glu_stoich, o2scalar) writephrq('cm', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3cm, fbuf, fglucm, fsom1cm, fsom2cm, fsom3, fferbcm, fmegacm, fmegh, glu_stoich, o2scalar) writephrq('cm', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \
fmegaro = fmeg fmegarm = fmeg fmegato = fmeg fmegatm = fmeg fmegh = fmeg totc = [0.045, 0.105, 0.104, 0.105, 0.074, 0.056] #mol total organic carbon totw = [0.013588, 0.005854, 0.011788, 0.006379, 0.010690, 0.006620] #kg water tots = [1.41, 9.15, 3.21, 8.62, 4.31 , 8.38] #g soil dry mass ph = [5.02, 4.84, 5.21, 4.54, 5.23, 4.95] cace = [6.370, 2.799, 0.057, 2.668, 1.026, 1.839] #mM acetate cfe2 = [0.789, 22.234, 0.102, 22.973, 15.665, 7.177] #mM Fe hdsp = [0.042528, 0.042528, 0.044005, 0.044005, 0.043575, 0.043575] #head space volume L index = 0 writephrq('co', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) index = 1 writephrq('cm', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3cm, fldoccm, fsom1cm, fsom2cm, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbcm, fmegacm, fmegh, sldoc) writephrq('cm', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3cm, fldoccm, fsom1cm, fsom2cm, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbcm, fmegacm, fmegh, sldoc) writephrq('cm', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \
fldoco = fldoc ffe3o = ffe3 ffe3co = ffe3o * 2.0 fldocco = fldoco fsom1co = fsom1o fsom2co = fsom2o fferbco = fferbo fmegaco = fmeg fmegh = fmeg totc = [0.045, 0.105, 0.104, 0.105, 0.074, 0.056] #mol total organic carbon totw = [0.013588, 0.005854, 0.011788, 0.006379, 0.010690, 0.006620] #kg water tots = [1.41, 9.15, 3.21, 8.62, 4.31 , 8.38] #g soil dry mass ph = [5.02, 4.84, 5.21, 4.54, 5.23, 4.95] cace = [6.370, 2.799, 0.057, 2.668, 1.026, 1.839] #mM acetate cfe2 = [0.789, 22.234, 0.102, 22.973, 15.665, 7.177] #mM Fe hdsp = [0.042528, 0.042528, 0.044005, 0.044005, 0.043575, 0.043575] #head space volume L index = 0 writephrq('co', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc)
fmegh = fmega fferbp = fmicp * 5 fmegap = fmicp * 2 fmeghp = fmegap totc = [0.1611, 0.0208, 0.0779] #mol total organic carbon totw = [0.0031, 0.0071, 0.008] #kg water tots = [6.94, 2.87, 2.02] #g soil dry mass ph = [4.24, 4.86, 4.95] cace = [5.75, 2.15, 10.76] #mM acetate T and P measured value cfe2 = [22.09, 20.24, 17.45] #mM Fe hdsp = [0.059, 0.059, 0.059] #head space volume L index = 0 writephrq('fa', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3, fldoc, fsom1, fsom2, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferb, fmega, fmegh, sldoc) writephrq('fa', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3, fldoc, fsom1, fsom2, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferb, fmega, fmegh, sldoc) writephrq('fa', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3, fldoc, fsom1, fsom2, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferb, fmega, fmegh, sldoc) index = 1 writephrq('ft', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3, fldoc, fsom1, fsom2, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferb, fmega, fmegh, sldoc) writephrq('ft', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3, fldoc, fsom1, fsom2, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferb, fmega, fmegh, sldoc) writephrq('ft', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \
0.006620] #kg water tots = [1.41, 9.15, 3.21, 8.62, 4.31, 8.38] #g soil dry mass ph = [5.02, 4.84, 5.21, 4.54, 5.23, 4.95] cace = [6.370, 2.799, 0.057, 2.668, 1.026, 1.839] #mM acetate cfe2 = [0.789, 22.234, 0.102, 22.973, 15.665, 7.177] #mM Fe hdsp = [0.042528, 0.042528, 0.044005, 0.044005, 0.043575, 0.043575] #head space volume L #consider organic acids, pH, and Fe(II) measurement microcosms have 20 g rather than 15 g wet soil. totc = np.asarray(totc) * 20.0 / 15.0 totw = np.asarray(totw) * 20.0 / 15.0 tots = np.asarray(tots) * 20.0 / 15.0 hdsp = np.asarray(hdsp) - 5.0e-3 # assume bulk density = 1 index = 0 writephrq('co', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) index = 1 writephrq('cm', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3cm, fldoccm, fsom1cm, fsom2cm, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbcm, fmegacm, fmegh, sldoc) writephrq('cm', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3cm, fldoccm, fsom1cm, fsom2cm, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbcm, fmegacm, fmegh, sldoc) writephrq('cm', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \
fldoco = fldoc ffe3o = ffe3 ffe3co = ffe3o * 2.0 fldocco = fldoco fsom1co = fsom1o fsom2co = fsom2o fferbco = fferbo fmegaco = fmeg fmegh = fmeg totc = [0.045, 0.105, 0.104, 0.105, 0.074, 0.056] #mol total organic carbon totw = [0.013588, 0.005854, 0.011788, 0.006379, 0.010690, 0.006620] #kg water tots = [1.41, 9.15, 3.21, 8.62, 4.31, 8.38] #g soil dry mass ph = [5.02, 4.84, 5.21, 4.54, 5.23, 4.95] cace = [6.370, 2.799, 0.057, 2.668, 1.026, 1.839] #mM acetate cfe2 = [0.789, 22.234, 0.102, 22.973, 15.665, 7.177] #mM Fe hdsp = [0.042528, 0.042528, 0.044005, 0.044005, 0.043575, 0.043575] #head space volume L index = 0 writephrq('co', -2, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 4, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc) writephrq('co', 8, totc[index], totw[index], tots[index], hdsp[index], ph[index], cace[index], cfe2[index], \ ffe3co, fldocco, fsom1co, fsom2co, fsom3, kfera, kferh, kmega, kmegh, kdecay, fferbco, fmegaco, fmegh, sldoc)