#!/usr/bin/env python
# coding: utf-8

# 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)
get_ipython().run_line_magic('matplotlib', 'inline')


# ### 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--')


# 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--')


# 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-')


# 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--')


# 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--')


# 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--')


# 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--')


# 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-.')


# In[ ]: