Commit 9eff7dc3 authored by Jorn Bruggeman's avatar Jorn Bruggeman

more carbon fluxes in plot_trends.py

parent e3b9c6fe
......@@ -79,12 +79,18 @@ for r in rs:
pylab.plot(div/pCO2s,delta13C_cals,'-')
pylab.figure(fig_fluxes.number)
pylab.subplot(1,3,1)
pylab.subplot(2,3,1)
pylab.plot(pCO2s,FB_ius/(FC_ius+FB_ius),'-')
pylab.subplot(1,3,2)
pylab.subplot(2,3,2)
pylab.plot(pCO2s,FB_xus/(FC_xus+FB_xus),'-')
pylab.subplot(1,3,3)
pylab.subplot(2,3,3)
pylab.plot(pCO2s,FB_vus/(FC_vus+FB_vus),'-')
pylab.subplot(2,3,4)
pylab.plot(pCO2s,FB_ius/(F_FIX+F_CAL),'-')
pylab.subplot(2,3,5)
pylab.plot(pCO2s,FB_xus/F_FIX,'-')
pylab.subplot(2,3,6)
pylab.plot(pCO2s,FB_vus/F_CAL,'-')
labels = [u'r = %i \u00b5m' % r for r in rs]
......@@ -118,22 +124,40 @@ pylab.savefig('delta13C_trends.png',dpi=150)
pylab.figure(fig_fluxes.number)
pylab.subplot(1,3,1)
pylab.subplot(2,3,1)
pylab.legend(labels)
pylab.grid(True)
pylab.xlabel('pCO2 (ppm)')
pylab.ylabel('FB_iu/(FC_iu+FB_iu) (-)')
pylab.subplot(1,3,2)
pylab.subplot(2,3,2)
pylab.legend(labels)
pylab.grid(True)
pylab.xlabel('pCO2 (ppm)')
pylab.ylabel('FB_xu/(FC_xu+FB_xu) (-)')
pylab.subplot(1,3,3)
pylab.subplot(2,3,3)
pylab.legend(labels)
pylab.grid(True)
pylab.xlabel('pCO2 (ppm)')
pylab.ylabel('FB_vu/(FC_vu+FB_vu) (-)')
pylab.subplot(2,3,4)
pylab.legend(labels)
pylab.grid(True)
pylab.xlabel('pCO2 (ppm)')
pylab.ylabel('FB_iu/(F_FIX+F_CAL) (-)')
pylab.subplot(2,3,5)
pylab.legend(labels)
pylab.grid(True)
pylab.xlabel('pCO2 (ppm)')
pylab.ylabel('FB_xu/F_FIX (-)')
pylab.subplot(2,3,6)
pylab.legend(labels)
pylab.grid(True)
pylab.xlabel('pCO2 (ppm)')
pylab.ylabel('FB_vu/F_CAL (-)')
pylab.show()
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment