In [1]:
import numpy as np
from matplotlib import pyplot as plt
import matplotlib as mpl
import datetime as dt
mpl.rc('xtick', labelsize=14)
mpl.rc('ytick', labelsize=14)
mpl.rc('legend', fontsize=16)
mpl.rc('axes', titlesize=16)
mpl.rc('figure', titlesize=16)
mpl.rc('axes', labelsize=16)
mpl.rc('font', size=16)
mpl.rc('legend', numpoints= 1)
mpl.rc('lines', markersize= 8)
%matplotlib inline
/home/eolson/anaconda3/lib/python3.5/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')
/home/eolson/anaconda3/lib/python3.5/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')
N¶
In [2]:
zz_rate_kapa_diat = 1.0
zz_rate_kapa_myri = 0.5
zz_rate_kapa_nano = 0.300
zz_rate_k_diat = 2.0
zz_rate_k_myri = 0.5
zz_rate_k_nano = 0.1
In [3]:
def OUP(zz_NO,zz_NH,zz_rate_kapa,zz_rate_k):
Oup_cell = zz_NO * zz_rate_kapa / (zz_rate_k + zz_NO * zz_rate_kapa + zz_NH)
return Oup_cell
In [4]:
def HUP(zz_NO,zz_NH,zz_rate_kapa,zz_rate_k):
Hup_cell = zz_NH / (zz_rate_k + zz_NO * zz_rate_kapa + zz_NH)
return Hup_cell
In [5]:
NH=np.arange(0,20,.1)
NO=5
plt.plot(NH,HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat),'b-')
plt.plot(NH,HUP(NO,NH,.9,zz_rate_k_diat),'-',color='purple')
plt.plot(NH,HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri),'m-')
plt.plot(NH,HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano),'r-')
plt.plot(NH,OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat),'b--')
plt.plot(NH,OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri),'m--')
plt.plot(NH,OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano),'r--')
Out[5]:
[<matplotlib.lines.Line2D at 0x7f7a5c554f98>]
In [6]:
NH=np.arange(0,20,.1)
NO=1
plt.plot(NH,HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat),'b-')
plt.plot(NH,HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri),'m-')
plt.plot(NH,HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano),'r-')
plt.plot(NH,OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat),'b--')
plt.plot(NH,OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri),'m--')
plt.plot(NH,OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano),'r--')
Out[6]:
[<matplotlib.lines.Line2D at 0x7f7a5c1f6c18>]
In [7]:
NH=np.arange(0,5,.1)
NO=.5
plt.plot(NH,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r-')
Out[7]:
[<matplotlib.lines.Line2D at 0x7f7a5745cdd8>]
In [8]:
NO=np.arange(0,10,.1)
NH=.1
plt.plot(NO,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b-')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m-')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r-')
plt.plot(NO,0.7E-04*(HUP(NO,NH,zz_rate_kapa_diat,5*zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,5*zz_rate_k_diat)),'b--')
plt.plot(NO,0.2E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m--')
plt.plot(NO,0.2E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r--')
Out[8]:
[<matplotlib.lines.Line2D at 0x7f7a57432198>]
In [9]:
NH=np.arange(0,10,.1)
NO=0
plt.plot(NH,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r-')
NO=np.arange(0,10,.1)
NH=0
plt.plot(NO,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b--')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m--')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r--')
Out[9]:
[<matplotlib.lines.Line2D at 0x7f7a573ffb38>]
In [10]:
NH=np.arange(0,10,.1)
NO=0
plt.plot(NH,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,2*zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,2*zz_rate_kapa_nano,zz_rate_k_nano)),'r-')
NO=np.arange(0,10,.1)
NH=0
plt.plot(NO,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b--')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m--')
plt.plot(NO,0.19E-04*(HUP(NO,NH,2*zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,2*zz_rate_kapa_nano,zz_rate_k_nano)),'r--')
Out[10]:
[<matplotlib.lines.Line2D at 0x7f7a573cd400>]
In [11]:
NO=np.arange(0,30,.1)
NH=.2
plt.plot(NO,0.565E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b-')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m-')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r-')
plt.plot(NO,0.62E-04*(HUP(NO,NH,2,9)+OUP(NO,NH,2,9)),'b--')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m--')
plt.plot(NO,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r--')
Out[11]:
[<matplotlib.lines.Line2D at 0x7f7a573197b8>]
In [12]:
NH=np.arange(0,10,.1)
NO=1
plt.plot(NH,0.564E-04*(HUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)+OUP(NO,NH,zz_rate_kapa_diat,zz_rate_k_diat)),'b-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m-')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)+OUP(NO,NH,zz_rate_kapa_nano,zz_rate_k_nano)),'r-')
plt.plot(NH,0.564E-04*(HUP(NO,NH,1.5*zz_rate_kapa_diat,1.5*zz_rate_k_diat)+OUP(NO,NH,1.5*zz_rate_kapa_diat,1.5*zz_rate_k_diat)),'b--')
plt.plot(NH,0.19E-04*(HUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)+OUP(NO,NH,zz_rate_kapa_myri,zz_rate_k_myri)),'m--')
plt.plot(NH,0.19E-04*(HUP(NO,NH,1.1*zz_rate_kapa_nano,.5*zz_rate_k_nano)+OUP(NO,NH,1.1*zz_rate_kapa_nano,.5*zz_rate_k_nano)),'r--')
plt.plot(NH,0.19E-04*(HUP(NO,NH,.5*zz_rate_kapa_nano,.5*zz_rate_k_nano)+OUP(NO,NH,.5*zz_rate_kapa_nano,.5*zz_rate_k_nano)),'k-.')
Out[12]:
[<matplotlib.lines.Line2D at 0x7f7a572f1048>]
In [ ]: