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

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


# ### SST -- data are from CentralSoG_SST_Data.ipynb

# In[2]:


NPGO_C_SST=[ 6.41475627,  6.48755541,  7.4785702 ,  9.25347408, 12.63790333,
       15.98057834, 18.54948443, 19.06939349, 15.93989528, 12.02364599,
        9.0698088 ,  6.8426298 ]


# In[3]:


NPGO_C_SST_SEM=[0.21721438, 0.33888804, 0.33143871, 0.15566903, 0.06565611,
       0.28891357, 0.22434856, 0.30571173, 0.24001654, 0.23185627,
       0.14982848, 0.25240781]


# In[4]:


NPGO_W_SST=[ 6.54155677,  6.32416194,  7.75323583, 10.39172893, 14.8838251 ,
       17.22590977, 19.5343776 , 19.48116382, 16.07562534, 12.52031955,
        9.00751125,  7.22206359]


# In[5]:


NPGO_W_SST_SEM=[0.10085818, 0.4890846 , 0.37887922, 0.15664109, 0.40595597,
       0.20241134, 0.40119105, 0.67801848, 0.50271186, 0.33871933,
       0.18173575, 0.33064517]


# In[6]:


Spring_C_SST=[ 9.78365537, 10.22426138,  9.4735278 ,  9.67848561]


# In[7]:


Spring_W_SST=[11.72465165, 11.13348375, 10.47916834, 10.70108273]


# In[8]:


Summer_C_SST=[17.50388524, 17.87908419, 18.36135449, 17.72161776]


# In[9]:


Summer_W_SST=[17.94220457, 19.22289418, 19.45471819, 18.36878465]


# ### Salinity -- data are from CentralSoG_SSS_Data.ipynb

# In[10]:


NPGO_C_Sal= [25.50095144, 24.93931884, 25.66983048, 24.71296806, 18.50090254,
       13.74989501, 15.83386163, 18.77032177, 22.11775614, 23.80291518,
       23.92629601, 24.63696603]


# In[11]:


NPGO_C_Sal_SEM=[0.25090128, 0.34827499, 0.22707056, 0.53063788, 1.29902834,
       1.59029185, 2.24347327, 0.97587449, 0.3152079 , 0.52615727,
       0.41454268, 0.26008036]


# In[12]:


NPGO_W_Sal=[24.16407734, 24.02108627, 23.46290609, 21.96345317, 15.89424179,
       14.72397722, 16.48597299, 19.74857647, 22.7618917 , 23.66667069,
       22.93491948, 23.81730112]


# In[13]:


NPGO_W_Sal_SEM=[0.8826463 , 0.99281972, 0.94991001, 0.90580449, 1.61168061,
       1.24163564, 1.76138473, 1.28528857, 0.81950455, 0.49974175,
       0.69538148, 1.10459364]


# In[14]:


Spring_C_Sal=[23.34870235, 23.97076454, 21.99257067, 22.53289721]


# In[15]:


Spring_W_Sal=[18.47972066, 20.29101052, 22.52859016, 20.46148007]


# In[16]:


Summer_C_Sal=[18.59444785, 18.55219596, 14.78432114, 12.54113959]


# In[17]:


Summer_W_Sal=[20.58331684, 16.97778168, 16.03245907, 14.35114464]


# ### Halocline Strength -- data are from CentralSoG_Halocline_Data.ipynb

# In[18]:


NPGO_C_Halo= [0.36154052, 0.49809309, 0.35279773, 0.59444239, 2.17599837,
       2.83692228, 2.06920697, 1.88577881, 1.31863951, 0.84843678,
       0.58477922, 0.50166281]


# In[19]:


NPGO_C_Halo_SEM=[0.0425931 , 0.06372294, 0.0536515 , 0.10031701, 0.37381109,
       0.46136078, 0.39131849, 0.18637298, 0.09991184, 0.06102145,
       0.08902362, 0.0766228]


# In[20]:


NPGO_W_Halo=[0.60806095, 0.64845463, 0.75646377, 1.24462616, 2.71714942,
       2.62485502, 2.29233551, 1.82686781, 1.266163  , 0.91992592,
       1.10688201, 0.64399621]


# In[21]:


NPGO_W_Halo_SEM=[0.10139684, 0.21724346, 0.11031644, 0.11941914, 0.33992397,
       0.31827148, 0.46300932, 0.33271079, 0.15843871, 0.11012063,
       0.17339882, 0.12728148]


# In[22]:


Spring_C_Halo=[1.02274943, 0.78320406, 1.33834988, 1.02001461]


# In[23]:


Spring_W_Halo=[1.787881  , 1.68205798, 1.24675077, 1.57429606]


# In[24]:


Summer_C_Halo=[1.76022024, 1.88017547, 2.33848005, 3.07700164]


# In[25]:


Summer_W_Halo=[1.3562591 , 2.24361753, 2.58961049, 2.80259067]


# ### PAR -- data are from CentralSoG_SurfacePAR_Data.ipynb

# In[26]:


NPGO_C_PAR= [12.04400661, 23.23528714, 36.51879487, 54.5113959 , 76.02100952,
       87.81722026, 99.97068933, 86.76086059, 61.21366767, 31.11335206,
       15.38539501, 10.86040996]


# In[27]:


NPGO_C_PAR_SEM=[0.98927229, 0.80121518, 2.21281148, 2.65062751, 4.36595316,
       6.5682177 , 3.54313251, 2.28383506, 2.1619708 , 1.18555839,
       0.89343873, 1.10119006]


# In[28]:


NPGO_W_PAR=[12.32470128, 26.11127048, 45.67482302, 63.15013291, 91.98044567,
       93.73015807, 98.88877244, 82.78904517, 55.05779156, 33.1923411 ,
       18.17157652,  9.76882253]


# In[29]:


NPGO_W_PAR_SEM=[1.47123463, 0.57532029, 3.05257932, 2.50705727, 4.39118173,
       5.75917672, 2.30639185, 2.69061531, 2.1721261 , 0.98056129,
       1.51775549, 0.60773262]


# In[30]:


Spring_C_PAR=[61.90657111, 52.45042499, 52.51039674, 55.86754087]


# In[31]:


Spring_W_PAR=[66.25929667, 68.03363515, 68.91267142, 64.53493222]


# In[32]:


Summer_C_PAR=[91.57734924, 97.41059098, 94.13498441, 82.94210228]


# In[33]:


Summer_W_PAR=[94.68697684, 93.73494871, 92.82362585, 85.96508285]


# ### Wind -- data are from HRDPSWinds_Data.ipynb

# In[34]:


NPGO_C_Wind= [5.255556106567383, 4.943120002746582 ,5.280871868133545, 4.541804313659668, 4.1618571281433105, 4.005492687225342, 4.127147674560547 ,4.005497932434082 ,4.025290012359619, 4.468635082244873, 5.065437316894531, 5.1902337074279785]


# In[35]:


NPGO_C_Wind_SEM=[0.22892044484615326, 0.23120790719985962, 0.19691279530525208 ,0.45478391647338867, 0.13002191483974457, 0.12630566954612732, 0.38189953565597534, 0.2072964608669281 ,0.23421116173267365, 0.2180320769548416, 0.26443806290626526, 0.5637392401695251]


# In[36]:


NPGO_W_Wind=[4.4879350662231445 ,4.2724103927612305 ,3.360086679458618 ,3.922926187515259, 3.1792898178100586, 3.588247537612915, 3.384995937347412, 3.394467830657959, 3.5997471809387207, 3.6979005336761475, 4.087904453277588, 4.875214576721191]


# In[37]:


NPGO_W_Wind_SEM=[0.595543622970581, 0.614099383354187, 0.2084141969680786 ,0.1963309645652771, 0.1585114598274231, 0.19811184704303741 ,0.1936599761247635, 0.2502512037754059, 0.17527444660663605, 0.1652664840221405, 0.3087571859359741, 0.3759503662586212]


# In[38]:


Spring_C_Wind=[5.280871868133545, 4.541804313659668, 4.1618571281433105]


# In[39]:


Spring_W_Wind=[3.360086679458618, 3.922926187515259, 3.1792898178100586]


# In[40]:


Summer_C_Wind=[4.005492687225342, 4.127147674560547, 4.005497932434082]


# In[41]:


Summer_W_Wind=[3.588247537612915, 3.384995937347412, 3.394467830657959]


# ### Fraser River Flow -- data are from FraserRiverFlow_Data.ipynb

# In[42]:


NPGO_C_FR= [991.5 ,  983.25,  979.  , 1024.  , 1077.  , 1012.5 ,  979.5 ,
        977.  ,  985.75, 1018.5 , 1066.25, 1114.  , 1071.  , 1102.5 ,
       1125.  , 1145.  , 1126.25, 1029.25,  952.  ,  929.5 ,  914.  ,
        883.  ,  883.  ,  899.75,  910.75,  913.5 ,  920.75,  928.5 ,
        941.  ,  922.  ,  901.75,  895.5 ,  910.5 ,  914.  ,  945.  ,
        966.  ,  960.75,  956.25,  934.75,  920.5 ,  910.5 ,  900.25,
        907.5 ,  925.25,  937.75,  935.25,  912.  ,  905.  ,  894.75,
        877.75,  857.25,  847.25,  864.75,  848.75,  850.  ,  840.75,
        814.5 ,  810.5 ,  824.  ,  816.5 ,  836.75,  835.5 ,  832.5 ,
        851.75,  874.25,  865.5 ,  859.5 ,  859.25,  865.  ,  878.25,
        896.25,  900.75,  894.75,  888.5 ,  905.75,  924.25,  911.5 ,
        919.25,  921.  ,  916.75,  920.  ,  916.  ,  917.25,  925.  ,
        928.75,  944.75,  974.5 , 1024.5 , 1056.25, 1132.5 , 1159.5 ,
       1154.75, 1111.  , 1119.75, 1106.25, 1104.25, 1081.  , 1099.  ,
       1100.25, 1109.  , 1139.75, 1162.25, 1187.5 , 1231.5 , 1275.75,
       1397.5 , 1497.5 , 1610.  , 1672.5 , 1777.5 , 1895.  , 2002.5 ,
       2102.5 , 2212.5 , 2307.5 , 2387.5 , 2452.5 , 2557.5 , 2697.5 ,
       2952.5 , 3150.  , 3280.  , 3315.  , 3355.  , 3460.  , 3520.  ,
       3560.  , 3607.5 , 3632.5 , 3687.5 , 3752.5 , 3850.  , 3960.  ,
       4250.  , 4480.  , 4645.  , 4842.5 , 5117.5 , 5562.5 , 6075.  ,
       6472.5 , 6757.5 , 6907.5 , 7037.5 , 7170.  , 7245.  , 7260.  ,
       7285.  , 7352.5 , 7442.5 , 7487.5 , 7552.5 , 7667.5 , 7755.  ,
       7580.  , 7530.  , 7627.5 , 7705.  , 7670.  , 7670.  , 7795.  ,
       8002.5 , 7987.5 , 7855.  , 7740.  , 7587.5 , 7560.  , 7695.  ,
       7895.  , 7742.5 , 7637.5 , 7717.5 , 7985.  , 8127.5 , 8080.  ,
       7945.  , 7850.  , 7952.5 , 7982.5 , 7932.5 , 7785.  , 7562.5 ,
       7507.5 , 7575.  , 7562.5 , 7422.5 , 7297.5 , 7267.5 , 7200.  ,
       7102.5 , 6985.  , 7025.  , 7115.  , 7112.5 , 6925.  , 6660.  ,
       6432.5 , 6280.  , 6225.  , 6100.  , 5907.5 , 5745.  , 5642.5 ,
       5542.5 , 5542.5 , 5515.  , 5402.5 , 5242.5 , 5092.5 , 4947.5 ,
       4827.5 , 4712.5 , 4635.  , 4597.5 , 4552.5 , 4500.  , 4417.5 ,
       4245.  , 4062.5 , 3952.5 , 3872.5 , 3825.  , 3777.5 , 3690.  ,
       3622.5 , 3547.5 , 3537.5 , 3492.5 , 3360.  , 3265.  , 3227.5 ,
       3177.5 , 3152.5 , 3130.  , 3100.  , 3105.  , 3190.  , 3277.5 ,
       3210.  , 3050.  , 2947.5 , 2882.5 , 2770.  , 2730.  , 2712.5 ,
       2715.  , 2657.5 , 2577.5 , 2522.5 , 2457.5 , 2382.5 , 2360.  ,
       2320.  , 2235.  , 2137.5 , 2085.  , 2032.5 , 1995.  , 1982.5 ,
       1995.  , 1957.5 , 1952.5 , 1957.5 , 1917.5 , 1850.  , 1787.5 ,
       1780.  , 1770.  , 1917.5 , 2197.5 , 2210.  , 2202.5 , 2175.  ,
       2175.  , 2170.  , 2127.5 , 2045.  , 2065.  , 1970.  , 1967.5 ,
       1890.  , 1822.5 , 1817.5 , 1802.5 , 1792.5 , 1717.5 , 1660.  ,
       1655.  , 1737.5 , 1730.  , 1750.  , 1750.  , 1825.  , 1777.5 ,
       1730.  , 1782.5 , 1815.  , 1795.  , 1750.  , 1707.5 , 1680.  ,
       1657.5 , 1677.5 , 1695.  , 1682.5 , 1662.5 , 1587.5 , 1567.5 ,
       1597.5 , 1640.  , 1612.5 , 1625.  , 1630.  , 1630.  , 1627.5 ,
       1582.5 , 1647.5 , 1720.  , 1670.  , 1655.  , 1652.5 , 1642.5 ,
       1637.5 , 1580.  , 1502.5 , 1415.  , 1350.  , 1285.  , 1245.  ,
       1232.5 , 1382.5 , 1364.5 , 1335.  , 1315.  , 1320.  , 1357.5 ,
       1337.5 , 1310.  , 1272.5 , 1247.5 , 1205.  , 1181.  , 1165.  ,
       1160.  , 1170.  , 1222.5 , 1242.5 , 1192.5 , 1177.5 , 1122.5 ,
       1047.5 , 1004.75,  976.5 ,  943.5 ,  922.5 ,  920.25,  939.75,
        994.5 ,  969.75,  955.75,  954.  ,  995.75, 1014.75,  989.25,
        982.75]


# In[43]:


NPGO_C_FR_Index=[ 1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,
        14,  15,  16,  17,  18,  19,  20,  21,  22,  23,  24,  25,  26,
        27,  28,  29,  30,  31,  32,  33,  34,  35,  36,  37,  38,  39,
        40,  41,  42,  43,  44,  45,  46,  47,  48,  49,  50,  51,  52,
        53,  54,  55,  56,  57,  58,  59,  60,  61,  62,  63,  64,  65,
        66,  67,  68,  69,  70,  71,  72,  73,  74,  75,  76,  77,  78,
        79,  80,  81,  82,  83,  84,  85,  86,  87,  88,  89,  90,  91,
        92,  93,  94,  95,  96,  97,  98,  99, 100, 101, 102, 103, 104,
       105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117,
       118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130,
       131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143,
       144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156,
       157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169,
       170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182,
       183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195,
       196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208,
       209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221,
       222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234,
       235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247,
       248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260,
       261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 273,
       274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286,
       287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299,
       300, 301, 302, 303, 304, 305, 306, 307, 308, 309, 310, 311, 312,
       313, 314, 315, 316, 317, 318, 319, 320, 321, 322, 323, 324, 325,
       326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 336, 337, 338,
       339, 340, 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, 351,
       352, 353, 354, 355, 356, 357, 358, 359, 360, 361, 362, 363, 364,
       365]


# In[44]:


#NPGO_C_FR_SEM=[]


# In[45]:


NPGO_W_FR=[1060.  , 1038.  , 1075.75, 1198.  , 1151.25, 1149.  , 1139.75,
       1166.75, 1140.  , 1137.5 , 1135.  , 1116.5 , 1089.25, 1015.5 ,
        988.25, 1001.25,  969.25, 1027.75, 1049.25, 1021.25, 1051.75,
       1068.  , 1085.25, 1197.25, 1218.25, 1183.75, 1166.25, 1184.25,
       1221.5 , 1288.25, 1318.75, 1714.  , 1476.25, 1362.25, 1263.25,
       1195.75, 1314.5 , 1371.75, 1309.5 , 1213.5 , 1171.75, 1136.5 ,
       1105.  , 1114.25, 1154.75, 1171.75, 1194.25, 1209.5 , 1225.25,
       1215.  , 1191.  , 1161.5 , 1141.5 , 1127.75, 1119.  , 1105.25,
       1074.75, 1055.5 , 1041.25, 1030.75, 1018.25,  995.5 ,  976.  ,
        947.5 ,  942.  ,  940.  ,  941.  ,  944.  ,  936.25,  929.75,
        956.  ,  984.  , 1073.  , 1151.5 , 1215.25, 1284.  , 1332.75,
       1340.75, 1315.75, 1330.25, 1339.  , 1393.25, 1433.  , 1431.75,
       1448.75, 1462.75, 1562.75, 1683.5 , 1964.  , 2157.  , 2218.  ,
       2250.  , 2281.  , 2223.75, 2142.  , 2092.5 , 2069.75, 2058.5 ,
       2046.5 , 2052.75, 2099.25, 2095.  , 2092.5 , 2112.5 , 2092.5 ,
       2052.5 , 2085.  , 2095.  , 2247.5 , 2387.5 , 2582.5 , 2770.  ,
       2957.5 , 3205.  , 3547.5 , 3775.  , 3880.  , 3992.5 , 4037.5 ,
       4090.  , 4262.5 , 4352.5 , 4392.5 , 4495.  , 4595.  , 4745.  ,
       4915.  , 4990.  , 5100.  , 5227.5 , 5442.5 , 5742.5 , 6022.5 ,
       6295.  , 6555.  , 6697.5 , 6840.  , 6992.5 , 7100.  , 7245.  ,
       7360.  , 7452.5 , 7615.  , 7787.5 , 7897.5 , 7980.  , 8030.  ,
       8067.5 , 8100.  , 8080.  , 7992.5 , 7927.5 , 7730.  , 7685.  ,
       7707.5 , 7682.5 , 7605.  , 7400.  , 7112.5 , 6835.  , 6597.5 ,
       6425.  , 6325.  , 6175.  , 6025.  , 5967.5 , 6017.5 , 6070.  ,
       6067.5 , 6145.  , 6130.  , 6072.5 , 6077.5 , 6117.5 , 6165.  ,
       6200.  , 6182.5 , 6165.  , 6227.5 , 6297.5 , 6262.5 , 6102.5 ,
       5997.5 , 5932.5 , 5842.5 , 5897.5 , 6040.  , 6060.  , 5982.5 ,
       5790.  , 5587.5 , 5485.  , 5392.5 , 5285.  , 5235.  , 5207.5 ,
       5140.  , 5032.5 , 4880.  , 4777.5 , 4852.5 , 4965.  , 4977.5 ,
       4875.  , 4762.5 , 4702.5 , 4570.  , 4460.  , 4410.  , 4350.  ,
       4217.5 , 4102.5 , 4002.5 , 3925.  , 3857.5 , 3805.  , 3780.  ,
       3777.5 , 3732.5 , 3632.5 , 3497.5 , 3390.  , 3330.  , 3255.  ,
       3187.5 , 3080.  , 3000.  , 2937.5 , 2915.  , 2850.  , 2805.  ,
       2802.5 , 2787.5 , 2785.  , 2732.5 , 2635.  , 2550.  , 2505.  ,
       2462.5 , 2432.5 , 2382.5 , 2360.  , 2312.5 , 2282.5 , 2235.  ,
       2167.5 , 2160.  , 2170.  , 2187.5 , 2240.  , 2230.  , 2147.5 ,
       2070.  , 2077.5 , 2105.  , 2127.5 , 2130.  , 2125.  , 2130.  ,
       2140.  , 2197.5 , 2182.5 , 2195.  , 2150.  , 2112.5 , 2177.5 ,
       2172.5 , 2287.5 , 2352.5 , 2350.  , 2362.5 , 2430.  , 2432.5 ,
       2370.  , 2292.5 , 2242.5 , 2150.  , 2062.5 , 2005.  , 1952.5 ,
       1925.  , 1882.5 , 1940.  , 2095.  , 2080.  , 2072.5 , 2257.5 ,
       2417.5 , 2377.5 , 2270.  , 2182.5 , 2137.5 , 2125.  , 2082.5 ,
       2027.5 , 1987.5 , 1922.5 , 1912.5 , 1910.  , 1880.  , 1977.5 ,
       2005.  , 1930.  , 1947.5 , 2000.  , 2237.5 , 2212.5 , 2227.5 ,
       2250.  , 2475.  , 2710.  , 2765.  , 2652.5 , 2520.  , 2517.5 ,
       2622.5 , 2602.5 , 2410.  , 2250.  , 2140.  , 2192.5 , 2180.  ,
       2225.  , 2255.  , 2220.  , 2100.  , 1992.5 , 1917.5 , 1835.  ,
       1835.  , 1985.  , 1890.  , 1749.75, 1631.75, 1552.  , 1483.5 ,
       1471.75, 1502.5 , 1497.5 , 1465.  , 1485.  , 1542.5 , 1612.5 ,
       1550.  , 1515.  , 1532.5 , 1537.5 , 1525.  , 1485.  , 1442.5 ,
       1425.  , 1402.5 , 1382.5 , 1422.5 , 1435.  , 1400.  , 1357.5 ,
       1312.5 , 1272.5 , 1220.  , 1180.  , 1163.  , 1184.5 , 1197.5 ,
       1160.25]


# In[46]:


NPGO_W_FR_Index=[ 1,   2,   3,   4,   5,   6,   7,   8,   9,  10,  11,  12,  13,
        14,  15,  16,  17,  18,  19,  20,  21,  22,  23,  24,  25,  26,
        27,  28,  29,  30,  31,  32,  33,  34,  35,  36,  37,  38,  39,
        40,  41,  42,  43,  44,  45,  46,  47,  48,  49,  50,  51,  52,
        53,  54,  55,  56,  57,  58,  59,  60,  61,  62,  63,  64,  65,
        66,  67,  68,  69,  70,  71,  72,  73,  74,  75,  76,  77,  78,
        79,  80,  81,  82,  83,  84,  85,  86,  87,  88,  89,  90,  91,
        92,  93,  94,  95,  96,  97,  98,  99, 100, 101, 102, 103, 104,
       105, 106, 107, 108, 109, 110, 111, 112, 113, 114, 115, 116, 117,
       118, 119, 120, 121, 122, 123, 124, 125, 126, 127, 128, 129, 130,
       131, 132, 133, 134, 135, 136, 137, 138, 139, 140, 141, 142, 143,
       144, 145, 146, 147, 148, 149, 150, 151, 152, 153, 154, 155, 156,
       157, 158, 159, 160, 161, 162, 163, 164, 165, 166, 167, 168, 169,
       170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180, 181, 182,
       183, 184, 185, 186, 187, 188, 189, 190, 191, 192, 193, 194, 195,
       196, 197, 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 208,
       209, 210, 211, 212, 213, 214, 215, 216, 217, 218, 219, 220, 221,
       222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234,
       235, 236, 237, 238, 239, 240, 241, 242, 243, 244, 245, 246, 247,
       248, 249, 250, 251, 252, 253, 254, 255, 256, 257, 258, 259, 260,
       261, 262, 263, 264, 265, 266, 267, 268, 269, 270, 271, 272, 273,
       274, 275, 276, 277, 278, 279, 280, 281, 282, 283, 284, 285, 286,
       287, 288, 289, 290, 291, 292, 293, 294, 295, 296, 297, 298, 299,
       300, 301, 302, 303, 304, 305, 306, 307, 308, 309, 310, 311, 312,
       313, 314, 315, 316, 317, 318, 319, 320, 321, 322, 323, 324, 325,
       326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 336, 337, 338,
       339, 340, 341, 342, 343, 344, 345, 346, 347, 348, 349, 350, 351,
       352, 353, 354, 355, 356, 357, 358, 359, 360, 361, 362, 363, 364,
       365]


# In[47]:


Spring_C_FR=[836.75,  835.5 ,  832.5 ,  851.75,  874.25,  865.5 ,  859.5 ,
        859.25,  865.  ,  878.25,  896.25,  900.75,  894.75,  888.5 ,
        905.75,  924.25,  911.5 ,  919.25,  921.  ,  916.75,  920.  ,
        916.  ,  917.25,  925.  ,  928.75,  944.75,  974.5 , 1024.5 ,
       1056.25, 1132.5 , 1159.5 , 1154.75, 1111.  , 1119.75, 1106.25,
       1104.25, 1081.  , 1099.  , 1100.25, 1109.  , 1139.75, 1162.25,
       1187.5 , 1231.5 , 1275.75, 1397.5 , 1497.5 , 1610.  , 1672.5 ,
       1777.5 , 1895.  , 2002.5 , 2102.5 , 2212.5 , 2307.5 , 2387.5 ,
       2452.5 , 2557.5 , 2697.5 , 2952.5 , 3150.  , 3280.  , 3315.  ,
       3355.  , 3460.  , 3520.  , 3560.  , 3607.5 , 3632.5 , 3687.5 ,
       3752.5 , 3850.  , 3960.  , 4250.  , 4480.  , 4645.  , 4842.5 ,
       5117.5 , 5562.5 , 6075.  , 6472.5 , 6757.5 , 6907.5 , 7037.5 ,
       7170.  , 7245.  , 7260.  , 7285.  , 7352.5 , 7442.5 , 7487.5 ,
       7552.5]


# In[48]:


Spring_W_FR=[1018.25,  995.5 ,  976.  ,  947.5 ,  942.  ,  940.  ,  941.  ,
        944.  ,  936.25,  929.75,  956.  ,  984.  , 1073.  , 1151.5 ,
       1215.25, 1284.  , 1332.75, 1340.75, 1315.75, 1330.25, 1339.  ,
       1393.25, 1433.  , 1431.75, 1448.75, 1462.75, 1562.75, 1683.5 ,
       1964.  , 2157.  , 2218.  , 2250.  , 2281.  , 2223.75, 2142.  ,
       2092.5 , 2069.75, 2058.5 , 2046.5 , 2052.75, 2099.25, 2095.  ,
       2092.5 , 2112.5 , 2092.5 , 2052.5 , 2085.  , 2095.  , 2247.5 ,
       2387.5 , 2582.5 , 2770.  , 2957.5 , 3205.  , 3547.5 , 3775.  ,
       3880.  , 3992.5 , 4037.5 , 4090.  , 4262.5 , 4352.5 , 4392.5 ,
       4495.  , 4595.  , 4745.  , 4915.  , 4990.  , 5100.  , 5227.5 ,
       5442.5 , 5742.5 , 6022.5 , 6295.  , 6555.  , 6697.5 , 6840.  ,
       6992.5 , 7100.  , 7245.  , 7360.  , 7452.5 , 7615.  , 7787.5 ,
       7897.5 , 7980.  , 8030.  , 8067.5 , 8100.  , 8080.  , 7992.5 ,
       7927.5]


# In[49]:


Summer_C_FR=[7667.5, 7755. , 7580. , 7530. , 7627.5, 7705. , 7670. , 7670. ,
       7795. , 8002.5, 7987.5, 7855. , 7740. , 7587.5, 7560. , 7695. ,
       7895. , 7742.5, 7637.5, 7717.5, 7985. , 8127.5, 8080. , 7945. ,
       7850. , 7952.5, 7982.5, 7932.5, 7785. , 7562.5, 7507.5, 7575. ,
       7562.5, 7422.5, 7297.5, 7267.5, 7200. , 7102.5, 6985. , 7025. ,
       7115. , 7112.5, 6925. , 6660. , 6432.5, 6280. , 6225. , 6100. ,
       5907.5, 5745. , 5642.5, 5542.5, 5542.5, 5515. , 5402.5, 5242.5,
       5092.5, 4947.5, 4827.5, 4712.5, 4635.]


# In[50]:


Summer_W_FR=[7730. , 7685. , 7707.5, 7682.5, 7605. , 7400. , 7112.5, 6835. ,
       6597.5, 6425. , 6325. , 6175. , 6025. , 5967.5, 6017.5, 6070. ,
       6067.5, 6145. , 6130. , 6072.5, 6077.5, 6117.5, 6165. , 6200. ,
       6182.5, 6165. , 6227.5, 6297.5, 6262.5, 6102.5, 5997.5, 5932.5,
       5842.5, 5897.5, 6040. , 6060. , 5982.5, 5790. , 5587.5, 5485. ,
       5392.5, 5285. , 5235. , 5207.5, 5140. , 5032.5, 4880. , 4777.5,
       4852.5, 4965. , 4977.5, 4875. , 4762.5, 4702.5, 4570. , 4460. ,
       4410. , 4350. , 4217.5, 4102.5, 4002.5, 3925. , 3857.5, 3805. ,
       3780. , 3777.5, 3732.5, 3632.5, 3497.5, 3390. , 3330. , 3255. ,
       3187.5, 3080. , 3000. , 2937.5, 2915. , 2850. , 2805. , 2802.5,
       2787.5, 2785. , 2732.5, 2635. , 2550. , 2505. , 2462.5, 2432.5,
       2382.5, 2360. , 2312.5, 2282.5]


# ### Nitrate -- data are from CentralSoG_0-10mNitrate_Data.ipynb

# In[51]:


NPGO_C_N= [24.3980628 , 23.67905457, 22.59146629, 14.48463374,  5.38950038,
        2.51176774,  1.97124377,  3.54658716,  8.74831538, 17.18580439,
       21.6594439 , 23.57027296]


# In[52]:


NPGO_C_N_SEM=[0.11667523, 0.11190603, 0.20137107, 1.23414193, 0.51169805,
       0.58081648, 0.51456234, 0.5983553 , 0.45166945, 0.24010748,
       0.22360953, 0.1472202]


# In[53]:


NPGO_W_N=[22.82788288, 22.19948811, 16.63627192,  6.90341105,  3.57964221,
        1.85326166,  1.53406205,  2.710954  ,  9.13413621, 16.421376  ,
       20.70562188, 22.49668357]


# In[54]:


NPGO_W_N_SEM=[0.36677255, 0.25092139, 0.68853362, 0.50627303, 0.44464386,
       0.24919049, 0.63251747, 1.03646634, 0.85953584, 0.45853951,
       0.25025754, 0.12691501]


# In[55]:


Spring_C_N=[12.46289449, 13.83035024, 16.02460427, 14.30295155]


# In[56]:


Spring_W_N=[7.57758729, 8.8820203 , 9.86462511, 9.83486753]


# In[57]:


Summer_C_N=[3.48006453, 3.20434749, 1.67761007, 2.34410947]


# In[58]:


Summer_W_N=[3.69164181, 1.13548262, 1.43987888, 1.86403363]


# ### Silicon -- data are from CentralSoG_0-10mSilicon_Data.ipynb

# In[59]:


NPGO_C_Si= [49.32323968, 50.37285438, 49.70968201, 34.4677553 , 12.64501118,
       12.50659466, 17.0418997 , 19.41833774, 26.0765079 , 35.45568785,
       43.45061737, 47.87295499]


# In[60]:


NPGO_C_Si_SEM=[0.41073373, 0.50735594, 0.36454512, 2.70717051, 2.41409781,
       1.4875741 , 2.96242003, 4.76044281, 3.73065234, 2.24611888,
       1.19261423, 0.90892549]


# In[61]:


NPGO_W_Si=[48.15701592, 49.07171776, 41.72099818, 17.42549544,  8.36526166,
       15.26847673, 21.98395522, 27.37391357, 33.72818989, 39.00032263,
       43.47364324, 46.34268646]


# In[62]:


NPGO_W_Si_SEM=[0.2275976 , 0.33726364, 1.34646478, 1.47386724, 0.67683366,
       0.9930824 , 2.5463565 , 2.47506929, 2.05869284, 1.55110053,
       0.78331193, 0.5134743]


# In[63]:


Spring_C_Si=[28.15483335, 29.61777566, 37.88429131, 33.43969765]


# In[64]:


Spring_W_Si=[19.10670492, 22.81162635, 21.68596449, 26.41137794]


# In[65]:


Summer_C_Si=[8.64361296, 10.58562337, 19.52603833, 26.5338348]


# In[66]:


Summer_W_Si=[15.81506055, 21.94194099, 20.0092991 , 28.40216005]


# In[67]:


fig, ax = plt.subplots(8,2,figsize=(10, 15),gridspec_kw={'width_ratios': [4, 1]})
bbox = {'boxstyle': 'round', 'facecolor': 'w', 'alpha': 0.9}
cmap = plt.get_cmap('tab10')
palette = [cmap(0), cmap(0.2), 'k', cmap(0.1), cmap(0.3)]
xticks=['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov',"Dec"]
border = 0.25

fig.tight_layout(pad=2.5)
c='k'

### SST

ax[0,0].errorbar(xticks, NPGO_C_SST, yerr= NPGO_C_SST_SEM, capsize=3,label='NPGO+ cold',linewidth=2)
ax[0,0].errorbar(xticks, NPGO_W_SST,yerr= NPGO_W_SST_SEM, capsize=3,linestyle='--',label='NPGO- warm',color='r',linewidth=2)
ax[0,0].set_title('Sea Surface Temperature',fontsize=14)
ax[0,0].legend(frameon=False,fontsize=12,loc=2)
ax[0,0].set_ylim(5,25)
ax[0,0].set_ylabel('Degrees C',fontsize=12)
ax[0,0].xaxis.set_tick_params(labelsize=12)
ax[0,0].yaxis.set_tick_params(labelsize=12)
ax[0,0].set_xticklabels([])
ax[0,0].text(-.4, 26.5, '(a)', fontsize=12,  color='k')

#ax[0,1].bp1 = color_boxplot(Spring_C_SST, 'royalblue', [1])
ax[0,1].boxplot(Spring_C_SST,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[0,1].boxplot(Spring_W_SST,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[0,1].boxplot(Summer_C_SST,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[0,1].boxplot(Summer_W_SST,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[0,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[0,1].set_ylim(5,25)
ax[0,1].set_ylabel('', fontsize=12)
ax[0,1].xaxis.set_tick_params(labelsize=12)
ax[0,1].yaxis.set_tick_params(labelsize=12)
#ax.legend([bp1["boxes"], bp2["boxes"], ['A', 'B'], loc='upper right')


### SSS

ax[1,0].errorbar(xticks, NPGO_C_Sal, yerr= NPGO_C_Sal_SEM, capsize=3,label='NPGO+ coldest',linewidth=2)
ax[1,0].errorbar(xticks, NPGO_W_Sal,yerr= NPGO_W_Sal_SEM, capsize=3,linestyle='--',label='NPGO- warmest',color='r',linewidth=2)
ax[1,0].set_title('Sea Surface Salinity',fontsize=14)
ax[1,0].legend((),frameon=False)
ax[1,0].set_ylim(5,30)
ax[1,0].set_ylabel('g/kg',fontsize=12)
ax[1,0].xaxis.set_tick_params(labelsize=12)
ax[1,0].yaxis.set_tick_params(labelsize=12)
ax[1,0].set_xticklabels([])
ax[1,0].text(-.4, 32, '(b)', fontsize=12,  color='k')


ax[1,1].boxplot(Spring_C_Sal,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[1,1].boxplot(Spring_W_Sal,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[1,1].boxplot(Summer_C_Sal,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[1,1].boxplot(Summer_W_Sal,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[1,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[1,1].set_ylim(5,30)
ax[1,1].set_ylabel('',fontsize=12)
ax[1,1].xaxis.set_tick_params(labelsize=12)
ax[1,1].yaxis.set_tick_params(labelsize=12)


### Halocline Strength
ax[2,0].errorbar(xticks, NPGO_C_Halo, yerr= NPGO_C_Halo_SEM, capsize=3,label='NPGO+ coldest',linewidth=2)
ax[2,0].errorbar(xticks, NPGO_W_Halo,yerr= NPGO_W_Halo_SEM, capsize=3,linestyle='--',label='NPGO- warmest',color='r',linewidth=2)
ax[2,0].set_title('Halocline Strength',fontsize=14)
ax[2,0].legend((),frameon=False,loc=1)
ax[2,0].set_ylim(0,4)
ax[2,0].set_ylabel('g/kg m$^{-1}$',fontsize=12)
#ax.invert_yaxis()
ax[2,0].xaxis.set_tick_params(labelsize=12)
ax[2,0].yaxis.set_tick_params(labelsize=12)
ax[2,0].set_xticklabels([])
ax[2,0].text(-.4, 4.3, '(c)', fontsize=12,  color='k')

ax[2,1].boxplot(Spring_C_Halo,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[2,1].boxplot(Spring_W_Halo,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[2,1].boxplot(Summer_C_Halo,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[2,1].boxplot(Summer_W_Halo,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[2,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[2,1].set_ylim(0,4)
ax[2,1].set_ylabel('',fontsize=12)
ax[2,1].xaxis.set_tick_params(labelsize=12)
ax[2,1].yaxis.set_tick_params(labelsize=12)


### Surface PAR

ax[3,0].errorbar(xticks, NPGO_C_PAR, yerr= NPGO_C_PAR_SEM, capsize=3,label='NPGO+ coldest',linewidth=2)
ax[3,0].errorbar(xticks, NPGO_W_PAR,yerr= NPGO_W_PAR_SEM, capsize=3,linestyle='--',label='NPGO- warmest',color='r',linewidth=2)
ax[3,0].set_title('Surface PAR',fontsize=14)
ax[3,0].legend((),frameon=False)
ax[3,0].set_ylim(0,125)
ax[3,0].set_ylabel('uE/m$^{2}$/s',fontsize=12)
ax[3,0].xaxis.set_tick_params(labelsize=12)
ax[3,0].yaxis.set_tick_params(labelsize=12)
ax[3,0].set_xticklabels([])
ax[3,0].text(-.4, 133, '(d)', fontsize=12,  color='k')


ax[3,1].boxplot(Spring_C_PAR,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[3,1].boxplot(Spring_W_PAR,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[3,1].boxplot(Summer_C_PAR,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[3,1].boxplot(Summer_W_PAR,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[3,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[3,1].set_ylim(0,125)
ax[3,1].set_ylabel('',fontsize=12)
ax[3,1].xaxis.set_tick_params(labelsize=12)
ax[3,1].yaxis.set_tick_params(labelsize=12)


### Wind Speed


ax[4,0].errorbar(xticks, NPGO_C_Wind, yerr= NPGO_C_Wind_SEM, capsize=3,label='NPGO+ coldest',linewidth=2)
ax[4,0].errorbar(xticks, NPGO_W_Wind,yerr= NPGO_W_Wind_SEM, capsize=3,linestyle='--',label='NPGO- warmest',color='r',linewidth=2)
ax[4,0].set_title('Wind Speed',fontsize=14)  #(with 2015)
ax[4,0].legend((),frameon=False)
ax[4,0].set_ylim(2,6)
ax[4,0].set_ylabel('m s$^{-1}$',fontsize=12)
ax[4,0].xaxis.set_tick_params(labelsize=12)
ax[4,0].yaxis.set_tick_params(labelsize=12)
ax[4,0].set_xticklabels([])
ax[4,0].text(-.4, 6.3, '(e)', fontsize=12,  color='k')


ax[4,1].boxplot(Spring_C_Wind,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[4,1].boxplot(Spring_W_Wind,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[4,1].boxplot(Summer_C_Wind,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[4,1].boxplot(Summer_W_Wind,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[4,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[4,1].set_ylim(2,6)
ax[4,1].set_ylabel('')
ax[4,1].xaxis.set_tick_params(labelsize=12)
ax[4,1].yaxis.set_tick_params(labelsize=12)


ax[5,0].plot(NPGO_C_FR_Index,NPGO_C_FR, linestyle='-',label='Cold Years',linewidth=2);
ax[5,0].plot(NPGO_W_FR_Index,NPGO_W_FR, color='r',linestyle='--',label='Warm Years',linewidth=2);
ax[5,0].set_title('Fraser River Flow',fontsize=14)
ax[5,0].legend((),frameon=False,bbox_to_anchor=(1, 1))
ax[5,0].set_ylim(0,10000)
ax[5,0].set_ylabel('m$^{3}$ s$^{-1}$',fontsize=12)
#ax.set_xlabel('Day of Year',fontsize=14)
ax[5,0].set_xlim(0,365)
ax[5,0].xaxis.set_tick_params(labelsize=12)
ax[5,0].yaxis.set_tick_params(labelsize=12)
ax[5,0].set_xticks([1,32,60,91,121,152,182,213,244,274,305,335])
#ax.set_xticklabels(['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'],fontsize=12)
ax[5,0].set_xticklabels([])
ax[5,0].text(4, 10600, '(f)', fontsize=12,  color='k')


ax[5,1].boxplot(Spring_C_FR,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[5,1].boxplot(Spring_W_FR,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[5,1].boxplot(Summer_C_FR,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[5,1].boxplot(Summer_W_FR,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[5,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[5,1].set_ylim(0,10000)
ax[5,1].set_ylabel('',fontsize=12)
ax[5,1].xaxis.set_tick_params(labelsize=12)
ax[5,1].yaxis.set_tick_params(labelsize=12)


ax[6,0].errorbar(xticks, NPGO_C_N, yerr= NPGO_C_N_SEM, capsize=3,label='NPGO+ coldest',linewidth=2)
ax[6,0].errorbar(xticks, NPGO_W_N,yerr= NPGO_W_N_SEM, capsize=3,linestyle='--',label='NPGO- warmest',color='r',linewidth=2)
ax[6,0].set_title('Depth Averaged Nitrate (0-10 m)',fontsize=14)
ax[6,0].legend((),frameon=False)
ax[6,0].set_ylim(0,30)
ax[6,0].set_ylabel('\u03bcmol N m$^{-3}$',fontsize=12)
ax[6,0].xaxis.set_tick_params(labelsize=12)
ax[6,0].yaxis.set_tick_params(labelsize=12)
ax[6,0].set_xticklabels([])
ax[6,0].text(-.4, 32, '(g)', fontsize=12,  color='k')


ax[6,1].boxplot(Spring_C_N,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[6,1].boxplot(Spring_W_N,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[6,1].boxplot(Summer_C_N,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[6,1].boxplot(Summer_W_N,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[6,1].set(xticks=[1.25,2.75], xticklabels=['',''])
ax[6,1].set_ylim(0,30)
ax[6,1].set_ylabel('',fontsize=12)
ax[6,1].xaxis.set_tick_params(labelsize=12)
ax[6,1].yaxis.set_tick_params(labelsize=12)


ax[7,0].errorbar(xticks, NPGO_C_Si, yerr= NPGO_C_Si_SEM, capsize=3,label='NPGO+ coldest',linewidth=2)
ax[7,0].errorbar(xticks, NPGO_W_Si,yerr= NPGO_W_Si_SEM, capsize=3,linestyle='--',label='NPGO- warmest',color='r',linewidth=2)
ax[7,0].set_title('Depth Averaged Silicon (0-10 m)',fontsize=14)
ax[7,0].legend((),frameon=False)
ax[7,0].set_ylim(0,60)
ax[7,0].set_ylabel('\u03bcmol N m$^{-3}$',fontsize=12)
ax[7,0].xaxis.set_tick_params(labelsize=12)
ax[7,0].yaxis.set_tick_params(labelsize=12)
#ax.set_xticklabels([])
ax[7,0].text(-.4, 63, '(h)', fontsize=12,  color='k')

ax[7,1].boxplot(Spring_C_Si,positions=[1], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[7,1].boxplot(Spring_W_Si,positions=[1.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[7,1].boxplot(Summer_C_Si,positions=[2.5], patch_artist=True,showmeans=False,boxprops=dict(facecolor='royalblue',color='royalblue'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[7,1].boxplot(Summer_W_Si,positions=[3], patch_artist=True,showmeans=False,boxprops=dict(facecolor='r',color='r'),capprops=dict(color=c),whiskerprops=dict(color=c),flierprops=dict(color=c, markeredgecolor=c),\
                medianprops=dict(color=c),widths=0.4)
ax[7,1].set(xticks=[1.25,2.75], xticklabels=['Spring','Summer'])
ax[7,1].set_ylim(0,60)
ax[7,1].set_ylabel('',fontsize=12)
ax[7,1].xaxis.set_tick_params(labelsize=12)
ax[7,1].yaxis.set_tick_params(labelsize=12)

#plt.savefig('Fig3a_SSTboxplots.png', bbox_inches='tight',dpi=1000,transparent=True)

#plt.show()


# In[ ]: