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

# ## Figure 7 - NPGO Paper

# In[1]:


import numpy as np
import netCDF4 as nc
import matplotlib.pyplot as plt
from salishsea_tools import bio_tools as bt, places
import xarray as xr
import os
import glob
get_ipython().run_line_magic('matplotlib', 'inline')


# ## Data are from the following notebooks:
# Experiments1and2_Nutrients_Data.ipynb\
# Experiments3and4_TempThreshold_Data.ipynb\
# Experiments5and6_Winds_Data.ipynb\
# Experiments7and8_ThermalForcing_Data.ipynb\
# Experiments9and10_Rivers_Data.ipynb

# #### SST Summary

# In[2]:


CY_SST_Orig=[ 6.12548543,  5.80986749,  7.48572357,  9.13112748, 12.73411507,
       16.15264921, 18.19417126, 18.16483525, 15.4019242 , 12.00080142,
        9.30544262,  6.33640542]


# In[3]:


CY_SST_Nuts=[ 6.13830953,  5.82549722,  7.48766262,  9.10256158, 12.68796849,
       16.20069323, 18.13362829, 18.18589498, 15.38127229, 11.96611296,
        9.38018667,  6.34800222]


# In[4]:


CY_SST_Thresh=[ 6.15109604,  5.8556036 ,  7.47975049,  9.10203945, 12.70656031,
       16.08675342, 18.27222176, 18.2109748 , 15.41443395, 11.9993017 ,
        9.32771591,  6.42471353]


# In[5]:


CY_SST_Winds=[ 5.59662968,  5.90259706,  7.60633249,  9.04497093, 13.09915475,
       16.51518926, 20.73906041, 19.60519318, 16.00876067, 11.95839157,
        9.37915939,  6.13304258]


# In[6]:


CY_SST_Thermal=[ 6.55805313,  4.6761362 ,  6.77186194, 10.27311911, 13.59493183,
       17.42428301, 18.35690238, 18.73194292, 16.01386188, 11.39800505,
        8.80717597,  6.78115408]


# In[7]:


CY_SST_Rivers=[ 6.11334258,  5.87921074,  7.48800428,  9.09396314, 12.75033495,
       16.07578677, 17.97372228, 18.18481324, 15.3735647 , 11.97804754,
        9.27182311,  6.43142438]


# #### Surface PAR Summary

# In[8]:


CY_PAR_Orig=[ 14.23656201,  25.43115062,  41.18548663,  59.62425692,
        84.90996976,  94.11623964, 100.57448297,  80.0413251 ,
        60.47416396,  29.24353915,  13.14300199,  10.14804782]


# In[9]:


CY_PAR_Nuts=[ 14.23679829,  25.43165744,  41.19107135,  59.75958975,
        85.18338548,  94.02408369, 100.44978733,  80.01338611,
        60.86782944,  29.30627043,  13.14189992,  10.14616278]


# In[10]:


CY_PAR_Thresh=[ 14.23694213,  25.42996031,  41.18085163,  59.76610693,
        84.90226449,  94.06975087, 100.70847597,  79.87635294,
        60.07499117,  29.22549197,  13.13997812,  10.14688224]


# In[11]:


CY_PAR_Winds=[ 14.22929443,  25.46042814,  40.97263571,  59.96603881,
        84.80457427,  94.07934195, 101.42659579,  80.66420632,
        60.70089652,  29.308822  ,  13.13657237,  10.14639563]


# In[12]:


CY_PAR_Thermal=[ 15.32649944,  26.11869233,  53.62521119,  62.5976493 ,
        90.68618976, 103.30374366,  93.93976964,  79.87383634,
        51.12972378,  31.54723225,  20.31534137,   8.03696216]


# In[13]:


CY_PAR_Rivers=[ 14.23691735,  25.43099464,  41.17917068,  59.83653825,
        85.16592265,  93.9387496 , 100.45387598,  80.0917271 ,
        60.47075832,  29.24338572,  13.14281766,  10.14914155]


# #### Halocline Summary

# In[14]:


CY_Halo_Orig=[0.24051982, 0.67116015, 0.49243159, 0.66998316, 1.90583354,
       2.39576078, 1.47431471, 1.41058524, 1.23433735, 0.80587692,
       0.68233926, 0.53663027]


# In[15]:


CY_Halo_Nuts=[0.23752941, 0.63396048, 0.48768099, 0.61599326, 1.56442715,
       2.35102508, 1.44770793, 1.34585816, 1.23384271, 0.83535316,
       0.66548312, 0.57292803]


# In[16]:


CY_Halo_Thresh=[0.2444939 , 0.61571177, 0.44665777, 0.6627658 , 1.81397565,
       2.34781621, 1.46515038, 1.38864157, 1.20166323, 0.83278059,
       0.68750955, 0.50213517]


# In[17]:


CY_Halo_Winds=[0.43660533, 0.37044598, 0.83569275, 0.73468569, 2.00271344,
       2.85756909, 3.35009795, 1.92678914, 1.44402677, 0.50235841,
       1.0366904 , 0.76952968]


# In[18]:


CY_Halo_Thermal=[0.27519593, 0.58094923, 0.44808606, 0.61038868, 1.6451783 ,
       2.3691709 , 1.48569156, 1.47417957, 1.26240626, 0.86845728,
       0.67721426, 0.53281211]


# In[19]:


CY_Halo_Rivers=[0.25069044, 0.62134257, 0.38932114, 0.95954486, 1.4824841 ,
       1.91582917, 1.28698258, 1.44149623, 1.14937059, 1.10101823,
       0.87432615, 0.51783257]


# ### Nitrate Summary

# In[20]:


CY_Orig_N=[24.75940603, 23.24451257, 21.67931533,  9.49131555,  6.21805258,
        2.60290254,  2.49482242,  5.34246862, 11.184194  , 15.54722773,
       21.87987576, 23.24806845]


# In[21]:


CY_Nuts_N=[23.9644255 , 22.29828934, 20.45139818,  8.97187906,  5.34151849,
        1.96037127,  1.81394895,  4.30505968, 11.52311063, 16.41310313,
       21.46072197, 22.33267238]


# In[22]:


CY_Thresh_N=[24.77422039, 23.29999955, 21.7438187 ,  9.87203474,  6.32937779,
        2.56494218,  2.649082  ,  4.58757776,  7.08316269, 15.02847707,
       21.88930665, 23.35885963]


# In[23]:


CY_Winds_N=[24.13891096, 23.43445999, 20.16081944,  9.53873272,  5.7397893 ,
        2.06754755,  0.95348431,  2.25127683,  8.106326  , 16.05059715,
       20.61628598, 23.32419151]


# In[24]:


CY_Thermal_N=[24.61595001, 23.26098825, 20.61280733,  7.68941071,  5.75685816,
        1.97729292,  2.72072923,  4.003872  , 11.9685279 , 19.67454109,
       23.13616984, 23.96606766]


# In[ ]:





# In[25]:


CY_Rivers_N=[24.74045154, 23.38186726, 21.69020601, 10.37701929,  5.74879707,
        2.70849594,  2.85380337,  5.48641413, 11.13281657, 15.64983911,
       21.55740055, 23.13323183]


# #### Silicon Summary

# In[26]:


CY_Orig_Si=[50.11491397, 51.4414015 , 50.00709959, 24.5001276 ,  9.95727285,
        8.81748288,  7.44068909,  9.67266691, 22.68525289, 29.63323648,
       41.74223596, 46.96035035]


# In[27]:


CY_Nuts_Si=[47.7632301 , 49.11254302, 47.48290684, 23.36581641,  8.2090725 ,
       11.07120201, 11.99591242, 12.30988367, 25.86081355, 32.83781192,
       41.9107794 , 46.15827035]


# In[28]:


CY_Thresh_Si=[50.13888235, 51.47213852, 50.01017314, 25.05295036,  9.56076795,
        6.89618133,  3.9125294 ,  4.51631545,  7.42875662, 19.71116328,
       38.33231021, 45.14306696]


# In[29]:


CY_Winds_Si=[50.06105751, 51.39125156, 49.65665279, 23.79865366,  8.11240424,
       10.9729251 , 13.07850414, 11.59949357, 21.74708727, 35.04296807,
       42.83359454, 47.61766334]


# In[30]:


CY_Thermal_Si=[50.04129223, 51.42741271, 47.97131586, 17.22230011,  7.05583654,
       10.13642248, 13.54628137, 15.74654144, 32.40864067, 40.63494107,
       46.43542205, 50.00465601]


# In[31]:


CY_Rivers_Si=[50.13663507, 51.513422  , 50.16886923, 27.00665494,  8.25802888,
        7.09006067,  6.40147506,  9.4362604 , 22.08336301, 30.14843401,
       41.93311685, 47.46025761]


# #### Diatom Summary

# In[32]:


CY_Diat_Orig=[ 0.23982064,  0.17371493,  3.13439335, 38.20802266, 35.26040455,
       20.78222047, 19.78868161, 10.83708083,  9.15037476,  1.56574763,
        0.69712752,  0.49628067]


# In[33]:


CY_Diat_Nuts=[ 0.23862619,  0.17259154,  3.03130916, 35.88205354, 31.55947094,
       15.65210445, 17.19965861, 10.24639004,  4.86930885,  1.22387145,
        0.67081294,  0.45610902]


# In[34]:


CY_Diat_Thresh=[ 0.23841182,  0.17512783,  3.16494256, 36.63234263, 34.80147671,
       22.14162946, 20.22215462, 14.54490841, 19.63002983,  3.79048587,
        1.16383968,  0.85648795]


# In[35]:


CY_Diat_Winds=[ 0.23412573,  0.1572702 ,  4.01715842, 35.63039506, 38.12295916,
       16.12961134, 13.7174575 ,  7.94355529,  4.27612511,  0.95877994,
        0.58034538,  0.41772271]


# In[36]:


CY_Diat_Thermal=[ 0.24531845,  0.1722018 , 11.15381515, 35.62884681, 32.77864728,
       17.98166513, 15.88818571,  8.06120314,  0.97470466,  0.65459443,
        1.07020377,  0.641704 ]


# In[ ]:





# In[37]:


CY_Diat_Rivers=[ 0.23932773,  0.17097974,  3.09714916, 35.42740652, 32.1912089 ,
       19.12590297, 20.34302593, 11.20385005,  9.36346982,  1.66558406,
        0.72738703,  0.4823474 ]


# In[39]:


fig,ax=plt.subplots(6,1,figsize=(12,18))

fig.tight_layout(pad=2.5)

months=np.arange(1,13)
xticks=['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov',"Dec"]


ax[0].plot(xticks,CY_SST_Nuts,'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='Exp 2 with WY Nutrients')
ax[0].plot(xticks,CY_SST_Thresh,'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='Exp 4 Increased Temp Threshold')
ax[0].plot(xticks,CY_SST_Winds,'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='Exp 6 with WY Winds')
ax[0].plot(xticks,CY_SST_Thermal,'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='Exp 8 with WY Thermal')
ax[0].plot(xticks,CY_SST_Rivers,'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='Exp 10 with WY Rivers')
ax[0].plot(xticks,CY_SST_Orig,'-',color='b',markersize=9,linewidth=4,label='CY Original')

ax[0].legend(frameon=False,loc=2,fontsize=12)
ax[0].set_ylabel('Degrees C',fontsize=14)
ax[0].set_xlabel('')
ax[0].set_ylim(0,25)
ax[0].set_title('SST',fontsize=18)
ax[0].xaxis.set_tick_params(labelsize=14)
ax[0].yaxis.set_tick_params(labelsize=14)
ax[0].set_xticklabels([])
ax[0].text(-.5, 26, '(a)', fontsize=15,  color='k')


ax[1].plot(xticks,CY_PAR_Nuts,'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[1].plot(xticks,CY_PAR_Thresh,'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[1].plot(xticks,CY_PAR_Winds,'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[1].plot(xticks,CY_PAR_Thermal,'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[1].plot(xticks,CY_PAR_Rivers,'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
ax[1].plot(xticks,CY_PAR_Orig,'-',color='b',markersize=9,linewidth=4,label='CY Original')

ax[1].legend((),frameon=False)
ax[1].set_ylabel('uE/m$^{2}$/s',fontsize=14)
ax[1].set_xlabel('')
ax[1].set_ylim(0,150)
ax[1].set_title('Surface PAR',fontsize=18)
ax[1].xaxis.set_tick_params(labelsize=14)
ax[1].yaxis.set_tick_params(labelsize=14)
ax[1].set_xticklabels([])
ax[1].text(-.5, 153, '(b)', fontsize=15,  color='k')


ax[2].plot(xticks,CY_Halo_Nuts,'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[2].plot(xticks,CY_Halo_Thresh,'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[2].plot(xticks,CY_Halo_Winds,'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[2].plot(xticks,CY_Halo_Thermal,'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[2].plot(xticks,CY_Halo_Rivers,'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
ax[2].plot(xticks,CY_Halo_Orig,'-',color='b',markersize=9,linewidth=4,label='CY Original')

ax[2].legend((),frameon=False)
ax[2].set_ylabel('g/kg m$^{-1}$',fontsize=14)
ax[2].set_xlabel('')
ax[2].set_ylim(0,5)
ax[2].set_title('Halocline Strength',fontsize=18)
ax[2].xaxis.set_tick_params(labelsize=14)
ax[2].yaxis.set_tick_params(labelsize=14)
ax[2].set_xticklabels([])
ax[2].text(-.5, 5.1, '(c)', fontsize=15,  color='k')


ax[3].plot(xticks,CY_Nuts_N,'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[3].plot(xticks,CY_Thresh_N,'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[3].plot(xticks,CY_Winds_N,'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[3].plot(xticks,CY_Thermal_N,'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[3].plot(xticks,CY_Rivers_N,'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
ax[3].plot(xticks,CY_Orig_N,'-',color='b',markersize=9,linewidth=4,label='CY Original')

ax[3].legend((),frameon=False)
ax[3].set_ylabel('mmol N m$^{-2}$',fontsize=14)
ax[3].set_xlabel('')
ax[3].set_ylim(0,30)
ax[3].set_title('Nitrate (0-10 m)',fontsize=18)
ax[3].xaxis.set_tick_params(labelsize=14)
ax[3].yaxis.set_tick_params(labelsize=14)
ax[3].set_xticklabels([])
ax[3].text(-.5, 31, '(d)', fontsize=15,  color='k')


ax[4].plot(xticks,CY_Nuts_Si,'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[4].plot(xticks,CY_Thresh_Si,'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[4].plot(xticks,CY_Winds_Si,'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[4].plot(xticks,CY_Thermal_Si,'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[4].plot(xticks,CY_Rivers_Si,'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
ax[4].plot(xticks,CY_Orig_Si,'-',color='b',markersize=9,linewidth=4,label='CY Original')

ax[4].legend((),frameon=False)
ax[4].set_ylabel('mmol N m$^{-2}$',fontsize=14)
ax[4].set_xlabel('')
ax[4].set_ylim(0,60)
ax[4].set_title('Silicon (0-10 m)',fontsize=18)
ax[4].xaxis.set_tick_params(labelsize=14)
ax[4].yaxis.set_tick_params(labelsize=14)
ax[4].set_xticklabels([])
ax[4].text(-.5, 61, '(e)', fontsize=15,  color='k')


ax[5].plot(xticks,CY_Diat_Nuts,'--',marker='^',markersize=9,markerfacecolor='mediumpurple', linewidth=1,color='k',label='CY with WY Nutrients')
ax[5].plot(xticks,CY_Diat_Thresh,'--',marker='s',markersize=9,markerfacecolor='violet', linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[5].plot(xticks,CY_Diat_Winds,'--',marker='*',markersize=11,markerfacecolor='lightskyblue', linewidth=1,color='k',label='CY with WY Winds')
ax[5].plot(xticks,CY_Diat_Thermal,'--',marker='d',markersize=9,markerfacecolor='lightgreen', linewidth=1,color='k',label='CY with WY Thermal')
ax[5].plot(xticks,CY_Diat_Rivers,'--',marker='o',markersize=9,markerfacecolor='burlywood', linewidth=1,color='k',label='CY with WY Rivers')
ax[5].plot(xticks,CY_Diat_Orig,'-',color='b',markersize=9,linewidth=3,label='CY Original')

ax[5].legend((),frameon=False)
ax[5].set_ylabel('mmol N m$^{-2}$',fontsize=14)
ax[5].set_xlabel('')
ax[5].set_ylim(0,40)
ax[5].set_title('Diatoms (0-100 m)',fontsize=18)
ax[5].xaxis.set_tick_params(labelsize=14)
ax[5].yaxis.set_tick_params(labelsize=14)
ax[5].text(-.5, 41, '(f)', fontsize=15,  color='k')

#ax[5].set_xticklabels([])

#fig.suptitle('Warm Year Experiments')

#plt.savefig('Figure7_CYExperiments.png', bbox_inches='tight',dpi=1000,transparent=False)


# In[41]:


fig,ax=plt.subplots(6,1,figsize=(10,19))
months=np.arange(1,13)
xticks=['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov',"Dec"]


fig.tight_layout(pad=2.5)

ax[0].plot(xticks,np.subtract(CY_SST_Nuts,CY_SST_Orig),'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='Exp 2 with WY Nutrients')
ax[0].plot(xticks,np.subtract(CY_SST_Thresh,CY_SST_Orig),'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='Exp 4 Increased Temp Threshold')
ax[0].plot(xticks,np.subtract(CY_SST_Winds,CY_SST_Orig),'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='Exp 6 with WY Winds')
ax[0].plot(xticks,np.subtract(CY_SST_Thermal,CY_SST_Orig),'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='Exp 8 with WY Thermal')
ax[0].plot(xticks,np.subtract(CY_SST_Rivers,CY_SST_Orig),'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='Exp 10 with WY Rivers')
#ax.plot(xticks,CY_Orig,'-',color='r',linewidth=3,label='CY Original')

ax[0].legend(frameon=True,bbox_to_anchor=(1, 1))
ax[0].set_ylabel('\u0394 Degrees C',fontsize=14)
ax[0].set_xlabel('')
ax[0].set_ylim(-6,6)
ax[0].axhline(y=0,linewidth=4,color='b')
ax[0].set_title('Cold Year Experiments - SST',fontsize=18)
ax[0].xaxis.set_tick_params(labelsize=14)
ax[0].yaxis.set_tick_params(labelsize=14)
ax[0].set_xticklabels([])
ax[0].text(-.5, 5, '(a)', fontsize=15,  color='k')


ax[1].plot(xticks,np.subtract(CY_PAR_Nuts,CY_PAR_Orig),'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[1].plot(xticks,np.subtract(CY_PAR_Thresh,CY_PAR_Orig),'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[1].plot(xticks,np.subtract(CY_PAR_Winds,CY_PAR_Orig),'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[1].plot(xticks,np.subtract(CY_PAR_Thermal,CY_PAR_Orig),'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[1].plot(xticks,np.subtract(CY_PAR_Rivers,CY_PAR_Orig),'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
#ax.plot(xticks,CY_Orig,'-',color='r',linewidth=3,label='CY Original')

ax[1].legend((),frameon=False)
ax[1].set_ylabel('\u0394 uE/m$^{2}$/s',fontsize=14)
ax[1].set_xlabel('')
ax[1].set_ylim(-15,15)
ax[1].axhline(y=0,linewidth=4,color='b')
ax[1].set_title('Surface PAR',fontsize=18)
ax[1].xaxis.set_tick_params(labelsize=14)
ax[1].yaxis.set_tick_params(labelsize=14)
ax[1].set_xticklabels([])
ax[1].text(-.5, 12.5, '(b)', fontsize=15,  color='k')




ax[2].plot(xticks,np.subtract(CY_Halo_Nuts,CY_Halo_Orig),'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[2].plot(xticks,np.subtract(CY_Halo_Thresh,CY_Halo_Orig),'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[2].plot(xticks,np.subtract(CY_Halo_Winds,CY_Halo_Orig),'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[2].plot(xticks,np.subtract(CY_Halo_Thermal,CY_Halo_Orig),'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[2].plot(xticks,np.subtract(CY_Halo_Rivers,CY_Halo_Orig),'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
#ax.plot(xticks,CY_Orig,'-',color='r',linewidth=3,label='CY Original')

ax[2].legend((),frameon=False)
ax[2].set_ylabel('\u0394 g/kg m$^{-1}$',fontsize=14)
ax[2].set_xlabel('')
ax[2].set_ylim(-2,2)
ax[2].axhline(y=0,linewidth=4,color='b')
ax[2].set_title('Halocline Strength',fontsize=18)
ax[2].xaxis.set_tick_params(labelsize=14)
ax[2].yaxis.set_tick_params(labelsize=14)
ax[2].set_xticklabels([])
ax[2].text(-.5, 1.7, '(c)', fontsize=15,  color='k')



ax[3].plot(xticks,np.subtract(CY_Nuts_N,CY_Orig_N),'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[3].plot(xticks,np.subtract(CY_Thresh_N,CY_Orig_N),'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[3].plot(xticks,np.subtract(CY_Winds_N,CY_Orig_N),'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[3].plot(xticks,np.subtract(CY_Thermal_N,CY_Orig_N),'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[3].plot(xticks,np.subtract(CY_Rivers_N,CY_Orig_N),'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
#ax.plot(xticks,CY_Orig_N,'-',color='r',linewidth=3,label='CY Original')

ax[3].legend((),frameon=False)
ax[3].set_ylabel('\u0394  \u03BCmol N m$^{-3}$',fontsize=14)
ax[3].set_xlabel('')
ax[3].set_ylim(-5,5)
ax[3].axhline(y=0,linewidth=4,color='b')
ax[3].set_title('Nitrate (0-10 m)',fontsize=18)
ax[3].xaxis.set_tick_params(labelsize=14)
ax[3].yaxis.set_tick_params(labelsize=14)
ax[3].set_xticklabels([])
ax[3].text(-.5, 4.2, '(d)', fontsize=15,  color='k')

ax[4].plot(xticks,np.subtract(CY_Nuts_Si,CY_Orig_Si),'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[4].plot(xticks,np.subtract(CY_Thresh_Si,CY_Orig_Si),'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[4].plot(xticks,np.subtract(CY_Winds_Si,CY_Orig_Si),'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[4].plot(xticks,np.subtract(CY_Thermal_Si,CY_Orig_Si),'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[4].plot(xticks,np.subtract(CY_Rivers_Si,CY_Orig_Si),'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
#ax.plot(xticks,CY_Orig_Si,'-',color='r',linewidth=3,label='CY Original')

ax[4].legend((),frameon=False)
ax[4].set_ylabel('\u0394 \u03BCmol N m$^{-3}$',fontsize=14)
ax[4].set_xlabel('')
ax[4].set_ylim(-30,30)
ax[4].axhline(y=0,linewidth=4,color='b')
ax[4].set_title('Silicon (0-10 m)',fontsize=18)
ax[4].xaxis.set_tick_params(labelsize=14)
ax[4].yaxis.set_tick_params(labelsize=14)
ax[4].set_xticklabels([])
ax[4].text(-.5, 25, '(e)', fontsize=15,  color='k')

ax[5].plot(xticks,np.subtract(CY_Diat_Nuts,CY_Diat_Orig),'--',marker='^',markersize=9,markerfacecolor='mediumpurple',linewidth=1,color='k',label='CY with WY Nutrients')
ax[5].plot(xticks,np.subtract(CY_Diat_Thresh,CY_Diat_Orig),'--',marker='s',markersize=9,markerfacecolor='violet',linewidth=1,color='k',label='CY with Increased Temp Threshold')
ax[5].plot(xticks,np.subtract(CY_Diat_Winds,CY_Diat_Orig),'--',marker='*',markersize=11,markerfacecolor='lightskyblue',linewidth=1,color='k',label='CY with WY Winds')
ax[5].plot(xticks,np.subtract(CY_Diat_Thermal,CY_Diat_Orig),'--',marker='d',markersize=9,markerfacecolor='lightgreen',linewidth=1,color='k',label='CY with WY Thermal')
ax[5].plot(xticks,np.subtract(CY_Diat_Rivers,CY_Diat_Orig),'--',marker='o',markersize=9,markerfacecolor='burlywood',linewidth=1,color='k',label='CY with WY Rivers')
#ax.plot(xticks,CY_Orig,'-',color='r',linewidth=3,label='CY Original')

ax[5].legend((),frameon=False)
ax[5].set_ylabel('\u0394 mmol N m$^{-2}$',fontsize=14)
ax[5].set_xlabel('',fontsize=14)
ax[5].set_ylim(-20,20)
ax[5].axhline(y=0,linewidth=4,color='b')
ax[5].set_title('Diatoms (0-100 m)',fontsize=18)
ax[5].xaxis.set_tick_params(labelsize=14)
ax[5].yaxis.set_tick_params(labelsize=14)
#ax[5].set_xticklabels([])
ax[5].text(-.5, 17, '(f)', fontsize=15,  color='k')

#plt.savefig('FigureS12_CYExperiments.png', bbox_inches='tight',dpi=1000,transparent=False)


# In[ ]: