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], \
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], \
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)
