{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import netCDF4 as nc\n",
    "import datetime as dt\n",
    "import subprocess\n",
    "import requests\n",
    "import matplotlib.pyplot as plt\n",
    "import cmocean\n",
    "import numpy as np\n",
    "import os\n",
    "import glob\n",
    "import dateutil as dutil\n",
    "from salishsea_tools import viz_tools\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "with nc.Dataset('/ocean/eolson/MEOPAR/NEMO-forcing/grid/mesh_mask201702_noLPE.nc') as fm:\n",
    "    tmask=np.copy(fm.variables['tmask'])\n",
    "    umask=np.copy(fm.variables['umask'])\n",
    "    vmask=np.copy(fm.variables['vmask'])\n",
    "    navlon=np.copy(fm.variables['nav_lon'])\n",
    "    navlat=np.copy(fm.variables['nav_lat'])\n",
    "    e3t_0=np.copy(fm.variables['e3t_0'])\n",
    "    e3u_0=np.copy(fm.variables['e3u_0'])\n",
    "    e3v_0=np.copy(fm.variables['e3v_0'])\n",
    "    e1t=np.copy(fm.variables['e1t'])\n",
    "    e2t=np.copy(fm.variables['e2t'])\n",
    "    e1v=np.copy(fm.variables['e1v'])\n",
    "    e2u=np.copy(fm.variables['e2u'])\n",
    "    A=fm.variables['e1t'][0,:,:]*fm.variables['e2t'][0,:,:]*tmask[0,0,:,:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "t0=dt.datetime(2015,2,6) # 1st start date of run\n",
    "#te=dt.datetime(2016,12,1)# last start date of runfnum=18\n",
    "stm=np.shape(tmask)\n",
    "SiN=2.0\n",
    "nlen=28\n",
    "#nlen=37\n",
    "dlist=[t0+dt.timedelta(days=ii*10) for ii in range(0,nlen)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "sdir0='/results/SalishSea/nowcast-green/'\n",
    "sdir1='/results/SalishSea/hindcast_v1/'\n",
    "sdir3='/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_Z3/'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "tmaskC=np.copy(tmask)\n",
    "tmaskC[:,:,370:490,:12]=0\n",
    "tmaskC[:,:,887:,30:70]=0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/results/SalishSea/nowcast-green/06feb15/SalishSea_1h_20150206_20150206_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/16feb15/SalishSea_1h_20150216_20150216_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/26feb15/SalishSea_1h_20150226_20150226_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/08mar15/SalishSea_1h_20150308_20150308_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/18mar15/SalishSea_1h_20150318_20150318_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/28mar15/SalishSea_1h_20150328_20150328_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/07apr15/SalishSea_1h_20150407_20150407_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/17apr15/SalishSea_1h_20150417_20150417_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/27apr15/SalishSea_1h_20150427_20150427_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/07may15/SalishSea_1h_20150507_20150507_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/17may15/SalishSea_1h_20150517_20150517_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/27may15/SalishSea_1h_20150527_20150527_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/06jun15/SalishSea_1h_20150606_20150606_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/16jun15/SalishSea_1h_20150616_20150616_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/26jun15/SalishSea_1h_20150626_20150626_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/06jul15/SalishSea_1h_20150706_20150706_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/16jul15/SalishSea_1h_20150716_20150716_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/26jul15/SalishSea_1h_20150726_20150726_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/05aug15/SalishSea_1h_20150805_20150805_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/15aug15/SalishSea_1h_20150815_20150815_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/25aug15/SalishSea_1h_20150825_20150825_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/04sep15/SalishSea_1h_20150904_20150904_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/14sep15/SalishSea_1h_20150914_20150914_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/24sep15/SalishSea_1h_20150924_20150924_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/04oct15/SalishSea_1h_20151004_20151004_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/14oct15/SalishSea_1h_20151014_20151014_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/24oct15/SalishSea_1h_20151024_20151024_ptrc_T.nc\n",
      "/results/SalishSea/nowcast-green/03nov15/SalishSea_1h_20151103_20151103_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/06feb15/SalishSea_1h_20150206_20150206_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/16feb15/SalishSea_1h_20150216_20150216_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/26feb15/SalishSea_1h_20150226_20150226_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/08mar15/SalishSea_1h_20150308_20150308_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/18mar15/SalishSea_1h_20150318_20150318_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/28mar15/SalishSea_1h_20150328_20150328_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/07apr15/SalishSea_1h_20150407_20150407_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/17apr15/SalishSea_1h_20150417_20150417_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/27apr15/SalishSea_1h_20150427_20150427_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/07may15/SalishSea_1h_20150507_20150507_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/17may15/SalishSea_1h_20150517_20150517_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/27may15/SalishSea_1h_20150527_20150527_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/06jun15/SalishSea_1h_20150606_20150606_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/16jun15/SalishSea_1h_20150616_20150616_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/26jun15/SalishSea_1h_20150626_20150626_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/06jul15/SalishSea_1h_20150706_20150706_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/16jul15/SalishSea_1h_20150716_20150716_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/26jul15/SalishSea_1h_20150726_20150726_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/05aug15/SalishSea_1h_20150805_20150805_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/15aug15/SalishSea_1h_20150815_20150815_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/25aug15/SalishSea_1h_20150825_20150825_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/04sep15/SalishSea_1h_20150904_20150904_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/14sep15/SalishSea_1h_20150914_20150914_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/24sep15/SalishSea_1h_20150924_20150924_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/04oct15/SalishSea_1h_20151004_20151004_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/14oct15/SalishSea_1h_20151014_20151014_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/24oct15/SalishSea_1h_20151024_20151024_ptrc_T.nc\n",
      "/results/SalishSea/hindcast_v1/03nov15/SalishSea_1h_20151103_20151103_ptrc_T.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150206-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150216-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150226-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150308-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150318-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150328-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150407-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150417-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150427-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150507-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150517-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150527-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150606-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150616-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150626-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150706-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150716-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150726-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150805-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150815-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150825-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150904-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150914-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20150924-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151004-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151014-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151024-*.nc\n",
      "/data/eolson/MEOPAR/SS36runs/CedarRuns/spring2015_HCMZ/SalishSea_1h_*_ptrc_T_20151103-*.nc\n"
     ]
    }
   ],
   "source": [
    "tlist=dlist\n",
    "SiGlobalTot=dict()\n",
    "SiTot=dict()\n",
    "BSiTot=dict()\n",
    "DiatTot=dict()\n",
    "changeSiGlobalTot=dict()\n",
    "for idir in (sdir0,sdir1,sdir3):\n",
    "    fformat1='%d%b%y/'\n",
    "    if idir.startswith('/data/eolson/MEOPAR/SS36runs/CedarRuns/'):\n",
    "        fformatT='SalishSea_1h_*_ptrc_T_%Y%m%d-*.nc'\n",
    "        fformatP='SalishSea_1h_*_ptrc_T_%Y%m%d-*.nc'\n",
    "    elif idir==sdir0:\n",
    "        fformatT='SalishSea_1h_%Y%m%d_%Y%m%d_ptrc_T.nc'\n",
    "        fformatP='SalishSea_1h_%Y%m%d_%Y%m%d_grid_T.nc'\n",
    "    elif idir==sdir1:\n",
    "        fformatT='SalishSea_1h_%Y%m%d_%Y%m%d_ptrc_T.nc'\n",
    "        fformatP='SalishSea_1h_%Y%m%d_%Y%m%d_aux_T_.nc'\n",
    "    sumSi=np.zeros((len(tlist),stm[2],stm[3]))\n",
    "    sumBSi=np.zeros((len(tlist),stm[2],stm[3]))\n",
    "    sumDiat=np.zeros((len(tlist),stm[2],stm[3]))\n",
    "    ind=-1\n",
    "    for idt0 in tlist:\n",
    "        ind=ind+1\n",
    "        cdir=idt0.strftime(fformat1).lower()\n",
    "        iffT=idt0.strftime(fformatT)\n",
    "        iffP=idt0.strftime(fformatP)\n",
    "        if idir.startswith('/data/eolson/MEOPAR/SS36runs/CedarRuns/'):\n",
    "            sffT=idir+iffT\n",
    "            sffP=idir+iffP\n",
    "        elif idir.startswith('/results/'):\n",
    "            sffT=idir+cdir+iffT\n",
    "            sffP=idir+cdir+iffP\n",
    "        f=nc.Dataset(glob.glob(sffT)[0])\n",
    "        print(sffT)\n",
    "        fP=nc.Dataset(glob.glob(sffP)[0])\n",
    "        if idir==sdir0:\n",
    "            e3t=np.expand_dims((1+fP.variables['sossheig'][0,:,:]/np.sum(e3t_0*tmask,1)),0)*e3t_0\n",
    "        elif idir==sdir1:\n",
    "            e3t=fP.variables['e3t'][:,:,:,:]\n",
    "        Vol=A*e3t\n",
    "        sumSi[ind,:,:]=1e-3*np.sum(tmaskC[0,:,:,:]*Vol[0,:,:,:]*f.variables['silicon'][0,:,:,:],0) #mmol/m3*m3*10^-3=mol\n",
    "        sumBSi[ind,:,:]=1e-3*np.sum(tmaskC[0,:,:,:]*Vol[0,:,:,:]*f.variables['biogenic_silicon'][0,:,:,:],0) #mmol/m3*m3*10^-3=mol\n",
    "        sumDiat[ind,:,:]=SiN*1e-3*np.sum(tmaskC[0,:,:,:]*Vol[0,:,:,:]*f.variables['diatoms'][0,:,:,:],0) #mmol/m3*m3*10^-3=mol\n",
    "    SiGlobalTot[idir]=np.sum(np.sum(sumSi+sumBSi+sumDiat,2),1)\n",
    "    SiTot[idir]=np.sum(np.sum(sumSi,2),1)\n",
    "    BSiTot[idir]=np.sum(np.sum(sumBSi,2),1)\n",
    "    DiatTot[idir]=np.sum(np.sum(sumDiat,2),1)\n",
    "    changeSiGlobalTot[idir]=[SiGlobalTot[idir][ii+1]-SiGlobalTot[idir][ii] for ii in range(0,len(tlist)-1)]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "collapsed": false
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7f8bd4591b70>]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEDCAYAAAA849PJAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJzt3XmcTfUbwPHPdxa77GtMopBUYmwRRbL+7FsbilApKoVS\nSMkWbRIiWpU1kew7Ycgu+27sYxnbLPf5/fEdWZoZd2buzLkz93m/Xvc1dzn3e54zd+Y893xXIyIo\npZTyPX5OB6CUUsoZmgCUUspHaQJQSikfpQlAKaV8lCYApZTyUZoAlFLKR3l9AjDGjDPGnDDGbHFj\n22rGmPXGmChjTPNbXmtrjNkVc2ubfBErpVTq4PUJABgP1HFz24NAO+CnG580xuQE+gAVgQpAH2NM\nDs+FqJRSqY/XJwARWQqcufE5Y0wxY8yfxph1xphlxpiSMdvuF5FNgOuWYmoD80TkjIiEAfNwP6ko\npVSaFOB0AIk0GugsIruMMRWBr4Aa8Wx/J3DohseHY55TSimfleoSgDEmC/AIMMkYc+3p9M5FpJRS\nqVOqSwDYaquzIlImAe85Ajx2w+NCwGIPxqSUUqmO17cB3EpEzgP7jDEtAIz10G3eNgd40hiTI6bx\n98mY55RSymd5fQIwxvwMrAJKGGMOG2PaA88A7Y0xG4GtQKOYbcsbYw4DLYBRxpitACJyBugPrI25\nfRDznFJK+Syj00ErpZRv8vorAKWUUsnDqxuBc+fOLUWKFHE6DKWUSjXWrVt3SkTyuLOtVyeAIkWK\nEBIS4nQYSimVahhjDri7rVYBKaWUj9IEoJRSPkoTgFJK+ShNAEop5aM0ASillI/SBKCUUj5KE4BS\nSvkorx4HoJRSiRXtimbGjhnsCdtDsRzFuCfnPRTNUZTM6TI7HZrX0ASglEpTLkVeYvyG8QxbNYw9\nYXv+83rBrAX/TQg33orlKEa2DNkciNg5mgCUUmnCyYsn+XLNl4xYO4LTl09T8c6KDHpiEI/f/Tj7\nwvax+8xu9oTtYfeZ3ew+s5s/d/9JaHjoTWXkyZSH0nlLc3+e+ymdt7S9n/d+smfI7tBRJS9NAEqp\nVG3X6V0MWzWM8RvHcyXqCg1LNKR75e5UDarKtVUDc2bMSbmC5f7z3osRF9kbtvffpLDj9A62ntzK\n+I3jCY8I/3e7QncUsgkhT+l/E8N9ee4jU2CmFDvO5KAJQCmVKq06tIohK4cw/Z/pBPoH0ubBNrz5\nyJuUzF3S7TIyp8vMA/ke4IF8D9z0vIhw8NxBtpzYYm8n7c9F+xZxNfoqAH7Gj2p3VaNFqRY0va8p\n+bPk9+jxpQSvXg8gODhYdDI4pRTYk/L+s/tZdXgVI9aOYOWhleTIkIOXy79MlwpdUuQEHOWKYs+Z\nPWw5sYV1oeuY9s80/jn1DwZzUzIokLVAsscSF2PMOhEJdmtbTyQAY0wd4DPAH/hGRAbe8vozQA/A\nABeAl0Rk4+3K1QSglG8SEfad3ce6o+tYF2pv60PXc+ayXcivSPYivF7pdV54+AWypMviaJxbT25l\n8rbJTNo2iW0nt2EwVA2qSotSLWhWqhkFsxZM0ZhSNAEYY/yBnUAt4DB2ycWnRGTbDds8AmwXkTBj\nTF2gr4hUvF3ZmgCU8g0Hzh5g9ZHV/57w14euJ+xKGACBfoE8kO8ByhUoZ28Fy1EmfxkC/LyvBnvb\nyW1M2jqJSdsmsfXkVgyGKkFVaH5fc5598FlyZcqV7DGkdAKojD2h14553AtARD6OY/scwBYRufN2\nZWsCUCrtG7pyKG/NewuAdP7peCDvA/+e6MsVKEfpvKVJH5De4SgTbvvJ7f9eGWw+sZki2Yvw5zN/\nUiJ3iWTdb0ongOZAHRHpEPP4OaCiiHSJY/vuQMlr28fyekegI0BQUFC5AwfcXttAKZWKiAi9F/Zm\nwPIBNC/VnF5Ve1E6b2nS+adzOjSPW3VoFY0mNsIlLmY+PZNKhSol274SkgBSdCoIY8zjQHtse0Cs\nRGS0iASLSHCePG6taqaUSmVc4uKVP15hwPIBvFj2RSY2m0jZAmXT5MkfoHLhyqxqv4rsGbJTY0IN\nft/xu9MhAZ5JAEeAwjc8LhTz3E2MMQ8C3wCNROS0B/arlEqFIqMjeW7ac4wMGcnbj7zNqAaj8Pfz\ndzqsZFcsZzFWtl/J/Xnvp/EvjRmzbozTIXkkAawF7jXG3G2MSQe0BmbcuIExJgiYCjwnIjs9sE+l\nVCp0OfIyTX5pwk+bf+Ljmh8zqNagfwdr+YK8mfOyqO0iaherTceZHem7uC9OdsVPcjO6iEQZY7oA\nc7DdQMeJyFZjTOeY178G3gdyAV/FfNhR7tZRKaXShnNXztFwYkOWHVjG1/W/plNwJ6dDckSWdFn4\nrfVvdJrZiX5L+nHk/BFGNhjpSK8mj+xRRP4A/rjlua9vuN8BiLXRVymV9p28eJI6P9Zh0/FN/NTs\nJ1qXbu10SI4K9A9kbMOxFMxakI+WfcSxi8eY2Gxiis9UqusBKKWS1aFzh6g2vhrbTm7jt9a/+fzJ\n/xpjDB/W+JCv6n3FH7v+oOZ3NTl16VSKxqAJQCmVbHae3knVb6ty9MJR5j47l3r31nM6JK/zUvmX\nmNJyChuPb+SRsY+wL2xfiu1bE4BSKllsOLaBR799lMuRl1nUdhGP3vWo0yF5rcYlGzP/ufmcunSK\nymMr83fo3ymyX00ASimPm7tnLo+Nf4z0/ulZ9vwyyhYo63RIXq9KUBVWvLCC9AHpefKHJ2+ajjq5\neN9kGkqpVOv81fO8NfctRq8fTak8pZj9zGyCsgU5HVaqcV+e+1jVfhUbj21MkUnuNAEopTxi3p55\ndPi9A4fPH6Z75e588PgHZAzM6HRYqU7BrAVTbAZRTQBKqSS58Vt/iVwlWP78cioXrux0WMoNmgCU\nUok2f+982s9oz6Fzh/RbfyqkCUAplWAXrl6g+9zujF4/muK5irPihRX6rT8V0gSglEoQ/dafdmgC\nUEq55cLVC7w17y1GrRul3/rTCE0ASqnbWn14NU9NeYr9Z/fzZuU36f94f/3WnwboQDClVJxc4mLQ\n8kFU/bYqLnGx7PllDH1yaLwnfxH45Rfo0QOiolIwWJVgegWglIrVsfBjPDftOebvnU/zUs0Z878x\nZM+QPd73hIbCSy/Bb7/ZxxkyQL9+KRCsShRNAEqp//hz95+0mdaG8IhwRjcYTYeyHeJduEUEvvsO\nunWDK1dg6FDYtAn694fHHoPHH0+52JX7NAEopf4VER3BuwveZeiqoZTOW5rFzRdTKk+peN9z6BB0\n6gSzZ0PVqjB2LBQvDuHhsHo1PPMMbNgAefOm0EEot2kbgFIKgN1ndlNlXBWGrhrKy8Evs6bDmnhP\n/iIwZgzcfz8sWQKff25/Fi9uX8+SBX79Fc6cgbZtweXyTJxbtsD5854py9dpAlBK8eOmH3l41MPs\nObOHqS2nMqL+iHgbevftg1q1oGNHKF8eNm+GV18Fv1vOKA8+CMOHw59/wrBhSY9zwgR46CGoUQMu\nXkx6eb5OE4BSPiw8Ipy209vy7LRnKZO/DBs6b6DJfU3i3N7lgi+/hAcegDVr4OuvYf58KFo07n10\n7gzNmkGvXrZKKLFGj4Z27aBMGfj7b3j2Wc9dVfgqjyQAY0wdY8wOY8xuY0zPWF43xpjPY17fZIzR\nycGVctjZK2cpP6Y8P2z6gT7V+7Co7aJ4p27+5x/bmPvqq7auf8sWW/cfT9swYF//5hsoVAhat4az\nZxMe6+ef233VqwcrVtiriunTbVdTlQQikqQb4A/sAYoC6YCNQKlbtqkHzAYMUAlY7U7Z5cqVk8TY\nt3+PXDgWnqj3KuUrnpv6nPj385e5u+fGuU1UlMhvv4nUri0CItmyiYwbJ+JyJXx/f/0lEhAg0rx5\nwt4/cKDdd5MmIlev2udcLpFXXrHPjxqV8FjSMiBE3D1/u7thnAVAZWDODY97Ab1u2WYU8NQNj3cA\nBW5XdmISwOHdB8Xv9YJS4pVH5MDZAwl+v1K+YPLWyUJfpM+iPrG+fuyYyIcfigQF2bNEwYIi/frZ\n55Ni8GBb3siRt9/W5RLp29du37q1SETEza9HRorUrSvi7y8yN+4c5nNSOgE0B7654fFzwJe3bDMT\nqHrD4wVAcBzldQRCgJCgoKAEH3xEVITcX+c14Z1MkvnDTDJ0xVCJjI5McDlKpVWhF0Il16BcEjw6\nWCKirp9VXS6RpUvtyTYw0J4datYUmTLlvyffxIqOtift9OlFNm6MezuXS6RnTxtDu3b2SiQ2586J\nPPCAyB13iGzd6pkYU7uEJACvawQWkdEiEiwiwXny5Enw+wP9A/mmUy8YsY27Dhen+7zuBI8OZvXh\nJLQ+KZVGiAgv/v4iFyMv8l3j7wj0D+T8efjqK9tjp1o125//lVdsnf/8+dC0KQQGemb/fn62J0/O\nnNCqVew9eUTg9ddh4EDbgDx2LPj7x17eHXfAzJmQKRPUrw/Hj3smTl/hiQRwBCh8w+NCMc8ldBuP\nqdQ4P9WzuTj39W/8UvpjTl06ReWxlXll1iucu3IuuXarlNcb9/c4Zu6cycCaAwnKdB/dusGdd9oT\nfmCgbaw9csQ2spYokTwx5MkDP/4IO3ZAly43v+Zy2akkPvsMuna1ienWrqW3CgqCGTPsyb9xY7h8\nOXniTpPcvVSI64YdTbwXuJvrjcD337JNfW5uBF7jTtmJbQQWEZk99ZKAyLdFP5Dzl89J19ldxa+f\nn+Qfml8mbp4orsS0YimViu05s0eyDMgiNSbUkJ27ouWBB0SMEXnuOdtAm9L/Eu+/b6t4vv/ePo6K\nstU9YKt/EhrPlCn2eFq2tFVNvoqUbAOw+6MesBPbG+jdmOc6A51j7htgRMzrm4mj/v/WW1ISgMsl\n8lDhU1KSbRL9w08iIhJyJETKjior9EXq/FBH9pzZk+jylUpNoqKj5NFxj8odH98h46cdkOzZRXLk\nEJkzx7mYIiNFqlUTyZzZ1t+3bm3PSP36JT4ZXWtkfvddz8aamqR4AkiuW1ISgIjIzz9GC4hMy9Ve\n5NIlERGJjI6UT1d9KlkGZJEMH2aQj5Z+JFejriZpP0p5u8HLBwt9kab9JogxIg89JLJ3r9NRiRw+\nLJIrl20UBpFBg5JWnssl0qGDLevbbz0SYqqTkARg7PbeKTg4WEJCQhL9/qgoKHHXZfIc3ciqD+Zj\n3uv972uHzx+m659dmbp9Knkz5+XRoEepGlSVqkFVKZO/DAF+Ok+eShs2H99M8Jhgcp1uQOink3nm\nGcPo0bbh1Bv88Qe0bAkDBsBrryW9vMhIqFsXli6FuXPtbKS+xBizTkSC3do2LScAsEPVX3oJFmeo\nQ/U946BgwZten71rNj9v+ZnlB5ez7+w+ADIHZqZy4cpULWwTQsVCFcmSLkuS4lDKCRHRETz4ZQV2\nHwtFRmxhWP88vPba7UfvprSoKAjw4Heus2ehcmXbMLxqVfI1aHsjTQA3uHIFihSO4uHTC5jd7hcY\nNy7ObQ+fP8yKgytYfnA5yw4uY9PxTQiCv/Hn4QIP82jQo9QqWos699SJd250pbxFs6/eYerJj8k2\nawa/Dfkf1as7HVHK2bsXKla0XU0bNLDTUNStCxnT+EqWmgBuMXCgnYjqbx6mzLqxUNa9qYjOXTnH\nqsOr/k0Ia46s4UrUFTo83IER9UeQzj9dkmNTKjlER8Pz76/k+4BHyXPkedb3/YZChZyOKuX984+d\nR2jyZDh50k5R3aiRTQZPPgnp0uC/sCaAW5w9C0FBQv2Iafxc8TNYvDhR18BXo67Sf2l/Plr2EVWD\nqjKl5RTyZtZVLpR3OXMGWj4bzoJ7y5Dljmj2vrWRPHfc4XRYjoqKgkWL7FrFU6bYc0L27HaQW6tW\ndnppT1ZBOSkhCcDrRgInh+zZ4aWXDL9GNGbP0sMwbVqiykkfkJ4Pa3zIz81+JuRoCOXHlGfDsQ0e\njlapxNu6FYKDYVHgW5ice5nVYYLPn/zBntxr1bID3Y4ft6OH//c/mDQJateGAgVsW2FSpqtOjXwi\nAYBdqzQg0DA05wB46y24ejXRZbUu3Zrlzy8n2hVNlXFVmLxtsgcjVSpxFi6ERx6Bs3lm4yr7NW9W\nfpNqd1VzOiyvky6dnTbiu+/gxAmYOhVq1rSPK1e26xz4Cp9JAAUKQLt2hm8vNOfY3ou2YjAJyhUs\nR0jHEB7M9yAtJrWgz6I+uERXp1DO+P57qFMH8hc/RGCz9tyf53761+jvdFheL0MGaNIEJk6Ew4ft\nNBVvvGHnI/IFPpMAwH7xj4z257NiX8CHH9r0nwT5s+RncdvFtCvTjg+WfkCLSS0IjwgHbCPcpElw\n4YInIlcqdiL2T7lN22iKtPyS0Kb3Ex51ju+bfE+GgAxOh5eq5Mhhf5crVth2Ap/g7ogxJ25JHQkc\nm5YtRe7IEiVn/XOKdO7skTJdLpcMXzVc/Pr5yYMjH5R9Yfvkq6/saMQOHTyyC6X+IyJCpH17EfJt\nlNw9Kwp9kSe/f1KnOEmCqCg7vfTdd4tcueJ0NIlDap4OOrn16AHnw/35uuK3dpHRzZuTXKYxhm6V\nujH7mdkcPHeQcqPK89aXS8mQwU5lu0HbiZWHnT8P9RpdZuz+Xvh1LofJuZcfmvzAn8/8SdEc8SzQ\nq+Ll7w+ffGIXvf/iC6ejSX4+0Q30VrVrw8YNLvZH3EmG8g/AnDkeGxq58/ROyg9ryHn/PfR8cARj\nOnekdGnbBU3HjilPOHoUHm07n72lOkPOPTxf5nmG1BpCrky5nA4tzahfH5Yvh927bbtAaqLdQG+j\nZ084fsKPCbV+gHnz7GQkHhK6pTjnP1lNMVOLgVs7UbFHX5YsSXTPU6VusmzdKe59uy17q9bizoJ+\nLGizgHGNxunJ38OGDLEjiPv1czqS5OWTVwAiUKkSnDol7AgoTcCVcDthyC3zBCVURASUKWMXpNi0\nOZquC17k2w3fkm/zQDKt78H27ZA+vYcOQvkUEaHnz98zZNMbSPpzdCjZky9avKsNvcnolVdg1Chb\nS3zffU5H4z69ArgNY2xbwN69hinPz7RDJ+vUscMDk2D4cNi+3dYdZs3iz5j/jeGp0k9x/IGe7Mv7\nOZ995qEDUD5l/9n9lB5ci8G72pLxYgnmNdvAmKf668k/mfXtC5kzQ/fuTkeSjNxtLXbilhy9gK6J\njhYpUUKkTBkR17z5IunSiVSt+u+6AQm1f79IpkwijRvf/HxEVIQ0mdhE6Iukf2SUHDvmgeCVz9hz\nZo9k61dI6HmHlHhmpJwJ8+GlrhxwbYEZJxfOSSi0F9Dt+fnB22/bHjrzXDXhhx9sB+DWre3EIQnU\ntav9+e+3/I0boW5dAlevZWLziVQvWI+rtTrTasD3njsIlaYdOHuA8l/W4NzFS9Q+uoxN4zqTI7vP\n/ss64rXX4O674c037dieNMfdTOHELTmvAERsP98777SrI506JSIjRth0/8ILCVqTbsYMuXk1o4kT\nRTJmtE8WKCBy7JhcjrwsQb1rCu/7yaCZvyTL8ai04/C5w5Lnw6JCj+xSrdU6iYhwOiLfNWmS/Vce\nPdrpSNxDSi0JCeQE5gG7Yn7miGWbwsAiYBuwFejqbvnJnQBE7Mk7fXqR4sVF9uyR6ytV9+rl1vsv\nXhQpUkSkVCmRq5eiRHr0sO+vUsVeN2bIIFKzpkhUlBw+ES4BHauKeT9Apm2fnrwHplKt0AuhUnBA\ncaFXVglutFouX3Y6It/mctl/57x5Rc6fdzqa20vJBDAY6BlzvycwKJZtCgBlY+5nxS4eX8qd8lMi\nAYiILFsmkjOn/YDXrHaJdOxofzWffnrb977zjt10yczzIrVr2wedO4tcjVln+Jtv7HN9+4qIyLAR\n54QOFSSgXzr5c9efyXlYKhU6EX5CggaVEt7JLPfXXS4XLjgdkRIRWb3a/hu/847TkdxeSiaAHUAB\nuX6i3+HGe34DarlTfkolABGRf/6xw78zZRKZMS1KpGlT++v58cc437Ntm0hgoEibhmEixYrZB6NG\n3byRyyXSpo2IMSLz5klkpMh9D5+RdK+WkQwfZpCFexcm85Gp1OLUxVNSbOiDwrsZpdgTiyQszOmI\n1I2eecbWFuzf73Qk8UvJBHD2hvvmxsdxbF8EOAjcEc82HYEQICQoKCj5fkuxOHZMJDhYxM9P5KvP\nIkQee0wkICDWLgAul8jjj4tkz3xVjmcqIpI/v8iKFbEXHB5u64jy5BE5ckTmzRMh0wnJ26+UZP4o\nsyw/sDyZj0x5u7DLYVJyeFmhd3opVH2uHD/udETqVgcO2Brdp592OpL4eTQBAPOBLbHcGt16wgfC\n4iknC7AOaOpucCl5BXBNeLhIgwb2N9Oj22WJfrCMSObM9hrwBj98Fy0gMpJOIhUqiBw+HH/BW7fa\ny4tq1UQiI6VhQ5HM+UKl6PB75Y6P75C1R9Ym41Epb3buyjl58POKwnuBkrfKLDl0yOmIVFzefdee\nG/76y+lI4uZ1VUBAIDAHeCMh5TuRAEREIiNtNT6IPNXkkly5u6RIrly2nkhEwg6ck3zpz0h5VktU\n2xfE7Va677+Xaw3MO3faGqPWnQ5KkU+LSI6BOeTv0L+T8aiUN7pw9YKUG1FFeC9AclSaLrt3Ox2R\nis/58yL58ok88kiCOgqmqJRMAENuaQQeHMs2BvgO+DSh5TuVAETsh/vxx/Y3VL3CJTmT+16RoCCR\nhQulS/bvxY8oCXn7l4T/Fbz4oi101ix54w3bNPD7sr1SaFghyfxRZpm0dVLyHJDyOhcjLkrlUY8J\n7/tJloq/ypYtTkek3DFmjP0X/vVXpyOJXUomgFzAgphuoPOBnDHPFwT+iLlfFRBgE7Ah5lbPnfKd\nTADX/Pij/aZequgl2Z+5lIRQVvyIki5NblPlE5dLl+zAg5w5JWzzIcmdW6R6dZHD545IpW8qCX2R\nnvN6SlR0lEePQ3mXy5GXpfo3tYQ+RjKU/1FCQpyOSLnr2poBRYq4f/GfklIsAST3zRsSgIjIwoUi\n2bKJFMh1RUpnOyj58kTJ2bNJKHDHDpGsWUUqVZKRX0QKiEyZInIl8op0+r3Tvwt7nL502mPHoLyH\ny+WS+t83FvoigRXGydKlTkekEmrePHv2LFHCDhTzpuqghCQAHVfuhscft7NEBGRKz5ZzhRn2qT/Z\nsiWhwOLF4Ztv4K+/6LCnF6VL2wmnJCo9Xzf4mtENRrN4/2KCRwez8dhGjx2H8g6j141h1p7p+M37\nhN8/eJ5HH3U6IpVQTzwBM2bYBWRatIDy5WHu3FS4lrC7mcKJm7dcAVwTGmrr/TyW7bt0EQGZ32+Z\ngMiQIddfWnVolRT8pKBk/DCj/Lz5Zw/tUDltX9g+yfBBFqFNDRnxlU7sltpFRYlMmCBy1132iuDx\nx53vIUQCrgB8cj0Ar3H1KlStCrt28ei9oVyIzHjT8pHHwo/RYlILlh9czpuV32TgEwMJ8AtwLl6V\nJC5xUWN8TZbtWUfROZvZvuouAvTjTBOuXrUrzPbvDydPQuPGdoH5++9P+Vh0PYDUIn16+PVXABqF\njmLjRjh06PrL+bPkZ0GbBXQp34VPVn1C7R9qc+rSKYeCVUk1Ys0IlhxcjGv2MIa9ryf/tCR9enj1\nVdizBz74ABYuhAcfhHbtYP9+p6OLmyYAp919N0yYQP0jo4D/rk6Zzj8dX9T7gm8bfcuKgysIHh3M\n+tD1DgSqkmLX6V30mN+DdAfqUjVzexo0cDoilRyyZoX33rOJ4PXXYeJE2+TXtatdd8rbaALwBo0a\nUbJRSYr672fm77FXybUr047lLyzHJS6qjKvC9xt1XYHUItoVTbvf2iFR6YmYPIYhgw3GOB2VSk65\nc8PQobBrF7RtCyNGQOfOTkf1X5oAvIR55mnqR89gwXwXly/Hvk1wwWBCOoZQqVAl2kxvQ8ffO3I5\nMo6NldcY/tdwVh5aiWvW5zR78k4qVXI6IpVSCheGMWPs+sIzZsD5805HdDNNAN6iQQMaZFjA5av+\nLFoU92Z5M+dl3nPz6FW1F2PWj6Hy2MrsPL0z5eJUCbLt5DZ6L+zN3VcbEbX+WQYMcDoi5YRWrWxD\n8YwZTkdyM00A3iJjRqo3yUlmwpn5W/xrzwX4BTCg5gD+ePoPDp8/TLnR5fhlyy8pFKhyV5QrinbT\n25ExIAsHvxpF506G4sWdjko5oVIlezUwcaLTkdxME4AXSf9sC2oxj1nTrro1oKTuvXX5u9PfPJjv\nQVpPac3Ls17mStSV5A9UuWXQ8kGsPbqWEru+ImN0Pt5/3+mIlFP8/KBlSztYLCzM6Wiu0wTgTWrV\nokHmxRw8mYktW9x7S+FshVncdjHdK3dnZMhIqoyrwp4ze5I3TnVbG49tpN+SftTM35LV41rSvTvk\ny+d0VMpJrVpBZCRMm+Z0JNdpAvAmgYHUa5IegFnTItx/m38gQ54cwozWM9gXto+yo8syZduU5IpS\n3UZEdARtp7clR8YcXJo0gnz54M03nY5KOS04GIoWhV+8qLZWE4CXKdC+HuUIYeZP5xL83v+V+B/r\nO62nZO6SNJ/UnK6zu3I16moyRKni8+HSD9l4fCMv5h/NqgW56dMHsmRxOirlNGPsVcCCBXa0sDfQ\nBOBtHn2U+lmWsmpHTk6fTvjbi2QvwrLnl9GtYjc+X/M5j377KPvC9nk+ThWrkKMhDFg2gGceeI5p\nHzeieHHo0MHpqJS3aNUKoqNh6lSnI7E0AXgbf38a/M/gwp8/p1xMVBHp/NMxvM5wprScws7TOyk7\nuiwrDq7wcKDqVleirtB2elvyZclHhdOfsW0bDBgAgYFOR6a8xYMPQokS3tMbSBOAFyrXtSr5OMbM\nsceTVE7T+5qyruM68mbOy5M/PMnCfQs9FKGKTZ9Ffdh2chsjao9lUN8cVKoETZs6HZXyJteqgZYs\ngdBQp6PRBOCV/CoEUy/rMv5cn4eoqKSVVSxnMZa0W8Ld2e+m/k/1mb1rtmeCVICdTn3pgaU0/Lkh\ng1cOpsPDHdg+ow5Hj8LgweiUD+o/WrWy6wZMnux0JJoAvJMx1H8yirNRWVn5eyIaAm6RP0t+Frdb\nzH2576OLu1MHAAAgAElEQVTRxEZM2+5F/dBSqShXFL9s+YWK31Sk+vjqrDq8ivervc/7FT5n4EBo\n2BBd6EXFqlQpKF3aO3oDJSkBGGNyGmPmGWN2xfzMEc+2/saYv40xM5OyT19R6+2HCSSCWSP2e6S8\n3Jlys7DtQsoVLEeLSS34efPPHinX14RHhPPZX59xz+f30HpKa85eOcvI+iM50O0A/R7vx7BBGQkP\nh48/djpS5c1atbKrDN44/bsTknoF0BNYICL3YheH7xnPtl2B7Uncn8+4o0JJqmVZz8yVcebUBMue\nITtzn51L1aCqPDP1Gcb9Pc5jZad1Ry8cpdf8XhQeXphuc7pROFthprWaxvZXttM5uDOZAjOxb5+d\n9fGFF+y3PKXi0qqV/TlpkrNxJDUBNAImxNyfADSObSNjTCGgPvBNEvfnUxpUD2fb5aLsW3bYY2Vm\nTZ+VP575g1rFatF+RntGrBnhsbLToi0ntvD8b89T5NMiDF45mCeKPsGq9qtY9vwyGpdsjL+fPwCX\nL0OvXhAQAP36ORy08nr33gtlyzrfGyipCSCfiFxryz4GxDXY/VPgbcB1uwKNMR2NMSHGmJCT3jJa\nwiH13ygBwKyhnr1wyhSYiRmtZ9CwREO6zO7CJys/8Wj5aUFEdASvzX6NB0Y+wK9bf+Xp4p34qdJO\nmkZNYulPlXj1VWjUyP4T58kDmTLZOt3XX4eCBZ2OXqUGrVrB2rWwd69zMdx2TWBjzHwgfywvvQtM\nEJHsN2wbJiI31VkYYxoA9UTkZWPMY0B3EXFrPaQ0vyawG0pkPMDdgYf583wVj5cdGR3Js9Oe5det\nv/LBYx/Qu1pvjHZb4eiFo7SY1IKVh1aSddtrhM/sg1zKedM22bLZ2R1vvN1zDzRvji71qNyyf79d\nEPDjj6FnfJXnCZSQNYFv+6cqIk/Es6PjxpgCIhJqjCkAnIhlsypAQ2NMPSADcIcx5gcRedadAH1d\ng0qn+XJxMOEh/5AluKRHyw70D+Snpj+RMSAj7y9+n0uRlxhQc4BPJ4GlB5bSclJLwiPCKb5xIuGr\nW9Htzesn+aAg+zNrVqcjValdkSJQsaK9cvRkAkiIpFYBzQDaxtxvC/x26wYi0ktEColIEaA1sFBP\n/u6r/0oRIkjPgiHJsw6wv58/4xqNo3O5zgxcMZBuf3bjdleFaZGIMHzVcGpMqEG2DNn4Ong1O6e1\nomdPu8j3iy9CnTq2cVdP/spTWrWCDRtgp0NrOiU1AQwEahljdgFPxDzGGFPQGPNHvO9UbqnaMCd3\n+F9k1p9+uLVIQCL4GT++qv/Vv/MH1f+pPgfOHkiWfXmj8IhwnpryFG/MfYP/lfgfazqsYerX95Mj\nh+3Ro1RyadnSDhZ0akxAkhKAiJwWkZoicq+IPCEiZ2KePyoi9WLZfrG79f/KSpcOnixzgpnnqyEh\n65JtP8YYhtUexhd1v2DpgaXc/9X9fL76c6Jd8a9OllzCw1NmPztP76TSN5WYtG0SA2sOZGrLqRw/\nmI3p0+HllyFz5pSJQ/mmO++EqlWd6w2kI4FTgQYv5CWUgvz92dJk3Y8xhi4VurD15a1Uu6saXf/s\nSpVxVdhyws3VaTxABPr2tY2syb1+6vR/phM8OpjjF48z59k59KjaA2MMw4fbxPvqq8m7f6XAVgNt\n24bbi0B5kiaAVKBu88wYXMyaEQ2u2/akTbK7st/FrKdn8WPTH9kTtoeHRz3MewvfS/blJiMjoX17\n248+XTq7iEqE++viuC3aFU2v+b1o8ksTSuYuybqO63iiqO3rcOIEjB8PbdroCl4qZTRvbpeMdKIa\nSBNAKpA3L1S45wwzL1SD5ctTZJ/GGJ5+4Gm2v7Kdp0o/xYfLPqTM12VYdmBZsuwvPNzOn/Ptt/D+\n+3airN27YdQoz+7n5MWT1PmxDgNXDKRj2Y4sfX4pQdmC/n19xAi4cgXeeMOz+1UqLvnywWOP2QSQ\n4v0vRMRrb+XKlRNlfdD7qhii5VibtxzZ/5zdc6TIp0WEvkin3zvJ2ctnPVZ2aKhI2bIi/v4iY8bY\n51wukZo1RXLlEgkL88x+lh9YLoWGFZL0/dPL2PVj//P6xYt2fw0bemZ/Srlr1CgREFm/PullASHi\n5jlWrwBSiQZN0yH4MXvKJVtXksKeLPYkW17awhuV3mDM+jGU+qoU0/+ZnuRy//kHKle2P2fMuL56\nljEwdCicOWMXVUkKl7gYsmII1cdXJ51/Ola2X8kLD/+3e8/48XD6NHTvnrT9KZVQTZuCv3/KVwPd\ndiSwk3Qk8HUiUCj3ZSqfmcXkPzJD3bqOxbL2yFo6/N6BTcc3cVe2u8gQkIEAv4CbboH+gf95Lk+m\nPDQs0ZDaxWqTMTAjK1bYap+AAJg1yy6afat27eDnn2HHDjtwJqHOXD5D2+ltmblzJs3ua8bYhmPJ\nliHbf7aLjrYrNeXODatW6Tz+KuXVrWu/CO3dm7S/v4SMBNYEkIp0bB/NxHEXOfV0V9L9+O1tt4+K\ngnHj4LffbN12Yk6gcYmMjuTLNV+yLnQdUa4oolxRRLoi/73/73PR15/bf3Y/YVfCyByYmYcy1mfN\nhGbcFVGPuTOzULRo7Ps5fBiKF4fGjeGnnxIW41+H/6LV5FaEXgjlkyc/oUuFLnGOcp4yxTbGTZ4M\nzZol8JehlAeMHw/PPw9//WVHCCdWQhKA4/X88d20DeBm06fbesL5GeqLXLoU53Yul8isWSKlStnt\n/fxESpYUOXMmBYONRURUhMzbM0+qDOgsdM8r9EUy9M8gjSc2lh82/hBnu0Lv3vY4Vq92bz8ul0uG\nrRwmAR8ESJFPi8iaw2tus71IxYoixYqJREUl9KiU8oywMJF06URefz1p5aBtAGlTzZqQPjCaWVdq\nwB+xD7T++2944gmoX992oZwyBRYutJeVjRvD1aspHPQN/E0gf371BCveGUmj3UeZ+9QSXiz3ImuP\nrOXZac+Sd2heGvzUgG///pYzl8/8+76337Y9obp3v30vibDLYTT5pQlvzH2D+vfWZ33H9ZS/s3y8\n71mxAlavtj1//P09caRKJVz27FC7Nvz6a4r09rbczRRO3PQK4L9qP+mSe/13i9Svf9PX1YMHRdq0\nETHG9mT5/HORq1evv++nn+y36NatRaKjUz7uK1dEWrWyMbzyys3ftKNd0bLy4Ep5c86bctfwu4S+\niF8/P3ngqwfk+enPy4g1I6TH56sF/ysybVrc+1hzeI0U+bSIBHwQIMNWDhOXy+VWbA0b2t/ZxYtJ\nPEilkuiHH+z/yLJliS+DBFwBOH6Sj++mCeC/vvjCfmo7uFfkoYfk3NT50quXSIYMIunTi/ToIXI2\njh6aAwfa9/bsmbIxi4i0b2/3PWiQrXKJi8vlkpAjIdJnUR+p+0NdyT04t9AXe3svUNK/WlY6TO8o\nY9aNkb9D/5aIqAhxuVzy+V+fS+AHgRI0PEhWHVrldlzbt9u43n/fAwepVBKdP2//l7t0SXwZCUkA\n2gicyuzbB0WLwqCnN5JlzhT6nu7CSfLybP0wPhyRg7vuivu9InZ+m6+/hpEjoXPnlIl5zhw7k2aP\nHjBwYMLeKyIcPHeQkKMhTFwWwuSVa8lYLITLcg6A9P7pCcoWxK4zu2hQvAETGk8gZ8actyn1uo4d\n4fvv4cABW82klNOaNbPVkkeOJK5KUnsBpXH332/nDgF4rNghhp5oQ7nwJbbPZP/+doapOERF2baA\n2bNt76AGyTw13/nzNt6sWWH9esiQIfFlidj2jQ0bhfnr9vDP+bWEHA1h04lN1LunHl0rdcXPuN+s\ndfw43HWX/bV9/XXi41LKk9asgbNnbZtfcicAx6t54rtpFVDsvvxSpHx5kRkzYqpTzpwR6d7ddiHI\nmFHk3XdFzp2L8/0XLoiUKyeSKZPI2rXJG2vHjrYX0ir3a2XitX69befo0SPpZfXubcvauTPpZSnl\nLdA2AB+1d6/IU0/ZjzVPHpERI0QiImLdNDRU5K67RPLlE9m3L3nCmTfPhtK9u2fLbdPGtnfs35/4\nMsLDRXLmFGnSxHNxKeUNEpIAtBtoWnL33Xa01Jo1dumqV16B0qVjnVc5f35bDXT1KtSrB2Fhng3l\nwgU7s2fx4nZFLU/68EM7UvKddxJfxrhxdpoJnfZB+TJNAGlR+fKwaJE98fv5QaNGMG/efza77z6Y\nPh327IEmTTw7RqBHDzh0yM7umTGj58oFuybvm2/aXLd2bcLfHxUFw4fDI4/Ym1K+ShNAWmUM/O9/\ndsHRoCDo3TvWUVTVq9uT9JIldhi6JwagLFxoexl165Z8J9gePdwfHHarqVNtbyr99q98nSaAtC59\nejvB/po18PvvsW7y9NPw8cd20rXevZO2u/BwW/Vzzz22qia5ZM1qF45ZujRhK4eJwJAhcO+9diI6\npXxZkhKAMSanMWaeMWZXzM8ccWyX3Rgz2RjzjzFmuzGmclL2qxKoTRt7Rn7vvTi/4vfoAZ062UTQ\no0fiV+Lq1cv2qR83DjJlSkLMbujQAUqWtFNFxDVDtssFly/DuXN2ta/ff4eQEFuFpNM+KF+X1CuA\nnsACEbkXWBDzODafAX+KSEngIWB7EverEiIw0C60u2mTne4yFsbAl1/Ciy/C4MFQoYLdPCGWLLFl\nvPoqPPpo0sO+nYAA+21+5047OC4oyDZu58hhF3MPCLAn+UyZ7Dwr+fLZ5pA8eWxOVMrXJWkgmDFm\nB/CYiIQaYwoAi0WkxC3bZAM2AEUlgTvTgWAeFB0NDz1kW0C3bLFnxzj8/rtNBGFhtgdP9+63/7Z8\n8aItXsQmjsyZPRx/HESgTx/bkJ0una3xiu9nunS2jbxMmZSJT6mUlmIjgY0xZ0Uke8x9A4Rde3zD\nNmWA0cA27Lf/dUBXEbkYR5kdgY4AQUFB5Q4cOJDo+NQtpk6148zHj4e2bePd9NQpO1XElClQpQpM\nmADFisW9fbdu8NlntvPRY495NGqlVAIkJAHctgrIGDPfGLMlllujG7eL+XYfWzYJAMoCI0XkYeAi\ncVcVISKjRSRYRILz5MnjzjEodzVpAmXL2tbT21Ty584NkybZeXK2bLHf7keNir3HzfLl8PnndtiB\nnvyVSj1umwBE5AkRKR3L7TfgeEzVDzE/T8RSxGHgsIisjnk8GZsQVEozxnbN2bfPttK6sfmzz8Lm\nzXbd3s6d7aCxo0evb3PpErzwgp1TJ6ETvSmlnJXURuAZwLW6hLbAb7duICLHgEPGmGttAzWx1UHK\nCXXq2M75/fvb7jFuKFzYzuj55Ze2obd0aZg40b723nuwaxeMHQtZsiRj3Eopj0tqAhgI1DLG7AKe\niHmMMaagMebGJateBX40xmwCygADkrhflVjGwEcf2a/xo0a5/TY/P1vFs2GDnd7hqafsItbDh9sr\ngxo1kjFmpVSy0OmgfdUTT9juOnv3Jvire1QUDBpke5YWLGjbCLJmTZ4wlVIJ49FGYJVGffghnDwJ\nX3yR4LcGBMC778L27bYBWE/+SqVOmgB8VaVKdjWYwYPt6hOJcM89tn1AKZU6aQLwZf3725P/sGFO\nR6KUcoAmAF9Wpgy0aGFbck+dcjoapVQK0wTg6/r1s535Bw92OhKlVArTBODr7rsPnnnGdvIPDXU6\nGqVUCtIEoOxsapGRMECHZyjlSzQBKDvL2wsv2IFhOvmeUj5DE4Cy3nvPDvft39/pSJRSKUQTgLIK\nFbJzOowfbyf3UUqleZoA1HW9etmVX0aOdDoSpVQK0ASgrsuXD2rXtstGevEcUUopz9AEoG7WvDkc\nOgRr1jgdiVIqmWkCUDdr2NAuIh/H4vFKqbRDE4C6WfbsUKuWXQ9Sq4GUStM0Aaj/atHCjgdYt87p\nSJRSyUgTgPqvhg3tpP+TJjkdiVJJI6JXsvEIcDoA5YVy5rQrhk2ebFd6N8bpiJSyzp61f5cHD8KF\nCxAefvtb8eLw1VdQvbrT0XudJCUAY0xO4BegCLAfaCkiYbFs9zrQARBgM/C8iFxJyr5VMmveHDp0\ngL//hrJlnY5G+TIRWLkSxoyBX3+Fy5ftl5IsWW6+Zc0KefNC0aLXn8uUCX75BR57zP49Dx4MOXI4\nfUTeQ0QSfQMGAz1j7vcEBsWyzZ3APiBjzONfgXbulF+uXDlRDjl1SsTfX6RnT6cjUb7q1CmR4cNF\nSpWyFTlZsoh06iQSEiLicrlfzsWLIm+/bf+e8+YV+fnnhL0/lQFCxM1zeFLbABoBE2LuTwAax7Fd\nAJDRGBMAZAKOJnG/KrnlygU1auigMJWyRGDJEjtF+Z13wuuv22/y33xjpyv/+msoVy5h1ZKZMsGg\nQRASAkFB8NRTUL8+7N+fbIeRWiQ1AeQTkWuTyB8D8t26gYgcAYYCB4FQ4JyIzI2rQGNMR2NMiDEm\n5OTJk0kMTyVJixawezds3Oh0JCqtO3kShg6FkiVtdc2sWbbKZsMGWL0a2re3iSApypSBv/6Czz6D\npUvh/vvhk08gKsr9MkRg506bkLp1s/GlZre7RADmA1tiuTUCzt6ybVgs788BLATyAIHAdOBZdy5P\ntArIYSdO2Mvmd991OhKVVu3bZ6t10qWz1TxVqoiMH2+rbZLTgQMiDRrYfT78sK1Wik10tMjmzSJf\nfinSsqVI/vzX+hWJ+PmJZMggMmFC8saaQCSgCiipbQA7gAIx9wsAO2LZpgUw9obHbYCv3ClfE4AX\nqFlTpHjxNF1nqhywa5fICy+IBASIBAbaJLBlS8rG4HKJTJokUqCAPZm//rrI2bMi69aJDBsm0rix\nSK5c10/4hQqJPPOMyKhRIv/8I3LsmMhjj9nXXn5Z5OrVlI0/DimZAIZwcyPw4Fi2qQhsxdb9G2xb\nwavulK8JwAuMHGn/TDZtcjoSlRZs3y7y7LP2hJs+vcirr4ocOuRsTGFhIp07279zY66f8IsWFWnX\nTuTbb0X27o39S1BkpMhbb9ntK1USOXw4xcO/VUomgFzAAmBXTFVRzpjnCwJ/3LBdP+CfmKqj74H0\n7pSvCcALHDtm/1nfe8/pSFRqtnmzSKtW9gSbMaPIG2+IHD3qdFQ3W75cpHt3kR9/THhSmjTJ9lLK\nm1dk0aJkCc9dCUkAxm7vnYKDgyUkJMTpMNTjj8Px47Btm9ORqNRmwwa7ytzUqZA5M3TpAm+8Yfvr\npzXbt0OTJrbjxKBB9jgdGERpjFknIsHubKsjgdXttWgBr7wCW7fanhNK3eryZdtN89gxewsNhTlz\n4Pff4Y47oHdv22smVy6nI00+991np1F/4QXo3t32Xho71g5Q81KaANTtNW1qv7lNnqwJwJetXw/z\n5v33RH/sGJw//9/tc+SAfv3gtdfsLLO+4I477BxaQ4dCz56wZYu9+ilZ0unIYqVVQMo91avD6dP2\nD1r5nvnzoUEDuHrV9sfPnx8KFLj556338+SxS4z6qoULoXVruHLFrrXdtGmK7FargJTnNW9uv8lt\n324vdZXvWLzYzhBbvDjMnWtP7ur2atSwU6o3bw7Nmtmroffe86rJFXU6aOWeZs3sT10pzLcsW2an\nTbj7bnsVoCf/hClc2I46btsW+vSB99/3qqlVNAEo9xQsCFWqaALwJatWQb16UKgQLFiQNnvupIT0\n6WHcODu1xYcf2gZxL0kCmgCU+1q0gE2b7FwoKm0LCYE6dSBfPluXrd/8k8bPD0aNgo4dYcAAePdd\nr0gCmgCU+7QayDf8/bddFzpXLli0yM7KqZLOzw9GjoROneDjj6FXL8eTgCYA5b5ChaByZV0qMi3b\nvNme/LNmtd/8Cxd2OqK0xc/Prk720kt2sFiPHo4mAe0FpBKmeXN480072vGee5yORnnStm1Qs6at\ns160CIoUcTqitMnPD0aMsL2BhgyxCWDwYEd6B+kVgEqY5s3tT60GSlt27rQnf39/+82/WDGnI0rb\njIEvv7Qj7IcOtSOHHbgS0ASgEiYoCCpU0GqgtGTPHttnPTra9vYpUcLpiHyDMfDFF/DqqzBsmL2y\nTuEkoFVAKuFatIC33oK9e+0C3Cr12r/fnvyvXLHVPqVKOR2RbzHGrlDm5wfDh4PLZX+mUHWQXgGo\nhNPeQGnD33/b5RcvXLCDvB54wOmIfJMx9qTfrZtNBt26pdiVgCYAlXB33w3BwZoAUisR2x2xUiW7\nHu68eXa9XOUcY2w10Ouvw+efQ9euKZIENAGoxGneHNautVUIKvU4dw5atYKXX7aNvhs2QLlyTkel\nwCaBTz6xbQErV8LFi8m/S50NVCXKnj22G+grr9i5Yvz8bA8SP7+b79/4XL582q/cSevXQ8uWNml/\n9JFtx/HT74BeRwQuXbIL6CSCzgaqkl+xYrY30IgR9uYOY2wD8nvvQenSyRtfWnLhAhw9amfjTEzj\noIgdfHRtJa4lS+y8Tso7GZPok39CJSkBGGNaAH2B+4AKIhLr13VjTB3gM8Af+EZEBiZlv8pL/Pmn\nHRDmctkuhC5X/PdXrbLd3n791VYhvf++NjzezpUr8MQTdqWpkiVt9U3Llu731jl3zk5CNnmyndht\nwgTInTt5Y1aph7uLB8d2w574SwCLgeA4tvEH9gBFgXTARqCUO+XrovBp0OnTIr17i2TNKgIiTZuK\nbNjgdFTeyeUSee45+3t6+22Rxx6zi6qDSOnSIv37i+zcGff7Q0JEihYV8fcXGTxYJDo65WJXjiEB\ni8InqQJQRLaLyI7bbFYB2C0ie0UkApgINErKflUqljOnXSR8/35bFTR/vu2B0qSJ7ZaorvvkE/j+\ne/v7GjTI9tM/csT2EsmWzf7+iheHsmVh4EDYt8++T8SOMn3kEYiIsPPRa32/io27mSK+G/FfATTH\nVvtce/wc8GU8ZXUEQoCQoKCg5EmRynucOSPSp49Itmz2m23Dhvabq6/74w8RPz+RFi3slUBsDh0S\nGTZMpGJF+7sDkfLlRWrXtvfr1xc5dSpl41aOw5NXAMaY+caYLbHckuVbvIiMFpFgEQnOkydPcuxC\neZMcOaBvX3tF0K+f/bYaHAyNGsHZs05H54x//rFryT74IHz7bdwNv4UK2X7jf/1lv/0PGmTbWpYs\nsZOMzZhhp3RWKg4e6QZqjFkMdJdYGoGNMZWBviJSO+ZxLwAR+fh25Wo3UB907pyt4ujTB955x66g\n5EvOnoWKFSEszC7KEhSU8DJcLq3u8WEJ6QaaEn8la4F7jTF3G2PSAa2BGSmwX5UaXavbbtLE1mOf\nP+90RCknOtp+89+3D6ZOTdzJH/Tkr9yWpL8UY0wTY8xhoDIwyxgzJ+b5gsaYPwBEJAroAswBtgO/\nisjWpIWt0rxevezVwMiRTkeScnr0gDlzbJ/9qlWdjkb5AB0JrLzXk0/aNYj37YOMGZ2OJnl99x20\nbQtdutixEkolkrdVASmVOO+8A8eP24bQtOyvv+DFF+20zMOGOR2N8iGaAJT3ql7dzlg5ZAhERjod\nTfI4csS2dxQqZEdIBwY6HZHyIZoAlPcyxrYF7N8PEyc6HY3nXb5sT/7h4fDbb9plU6U4TQDKuzVo\nYCeOGzjQdm9MSbNn2/lzDh/2fNkittpn7Vr44QedHE85QhOA8m5+ftCzJ2zbBr//nnL7XbfOTlg3\nezbUqgUnTni2/MGD4ccf7TiHRjozinKGJgDl/Vq1squQDRiQMkvlHTxorzxy57b18gcOQO3anhuZ\nPHCgTWqtWtmGbqUcoglAeb+AAHj7bTsl8qJFybuv8+ftAjeXLsEff9j1C6ZNg61boW5dW1+fWCL2\nxN+rFzz9tJ3oLYUW/1YqNpoAVOrQrh3kz2+vApJLZKQ94f/zD0yZAvffb5+vXds2Qq9da6trrlxJ\neNnR0fDSS3a+npdesid/7fGjHKYJQKUOGTLYFa0WLLBXAp4mYgdhzZ1rRx8/8cTNrzdtascjLFxo\nk0RCuqVGRMAzz8CoUfbb/4gROl2D8gr6V6hSj86dIXt2+Pi28wgm3NChMHq0raLp0CH2bZ57zk7T\nMHOmvR8dfftyL12Cxo3hl19sw++AAVrto7yGJgCVemTNCq++CtOn215BnjJ5sm1jaNnSLpYen5de\nsifyX36BTp3i75p67hzUqWOXzhw92i7KopQX0QSgUpfXXoNMmWxPGk/46y/7bb5yZRg/3r2qmbfe\ngt69YexYWy0VW8+kkyfh8cftOsg//2z7/CvlZTQBqNQld27o2BF++smOEE6KffugYUMoWNCOxE3I\nhHMffABdu8Jnn9nF7W906BA8+qhtTJ4xw3b3VMoLaQJQqc+bb9pv6kOGJL6MsDA7yjcqynb3TOjq\nc8bA8OHQvr0dzDV4sH1+1y47lXNoqJ3auW7dxMeoVDILcDoApRKsUCFo0wbGjbPfvvPlS9j7IyKg\nWTPYswfmzYMSJRIXhzG2Z094uJ3L/8wZ21PI5bLjFcqWTVy5SqUQvQJQqVOPHvZE/umnCXufiG28\nXbTI1uFXr560OPz9bZ/+Bg1sH/906WDZMj35q1RBE4BKne69187VM2KEe1M0HDliJ11r0cI29vbp\nYxt/PSEwECZNsl1JV6yAkiU9U65SyUxXBFOp14YN8PDDtuvmrXPqnDoFixfbgVsLF8KOHfb5HDls\nvf3gwdofX6VJCVkRTNsAVOpVpoxtZB0+HJ5/HkJCrp/wN22y22TJAtWqXV9x66GHdBSuUjGSlACM\nMS2AvsB9QAUR+c/XdWNMYeA7IB8gwGgR+Swp+1XqX++8Y7tcFixoH2fIAFWq2J45NWpAcLDOuaNU\nHJJ6BbAFaAqMimebKOBNEVlvjMkKrDPGzBMRDw7lVD6ralXbJz8y0p7wK1WySUApdVtJSgAish3A\nxFOXKiKhQGjM/QvGmO3AnYAmAOUZ773ndARKpUopWhlqjCkCPAysjmebjsaYEGNMyMmTJ1MqNKWU\n8jm3vQIwxswH8sfy0rsi8pu7OzLGZAGmAN1E5Hxc24nIaGA02F5A7pavlFIqYW6bAETkidttczvG\nmEB0tmMAAAOJSURBVEDsyf9HEZma1PKUUkolXbJXARnbQDAW2C4iw5J7f0oppdyTpARgjGlijDkM\nVAZmGWPmxDxf0BjzR8xmVYDngBrGmA0xt3pJiloppVSSJbUX0DRgWizPHwXqxdxfDuiQS6WU8jI6\nJFIppXyUJgCllPJRXj0ZnDHmJHAgkW/PDZzyYDjeJq0fH6T9Y9TjS/288RjvEhG3Vjjy6gSQFMaY\nEHdnxEuN0vrxQdo/Rj2+1C+1H6NWASmllI/SBKCUUj4qLSeA0U4HkMzS+vFB2j9GPb7UL1UfY5pt\nA1BKKRW/tHwFoJRSKh6aAJRSykeluQRgjKljjNlhjNltjOnpdDzJwRiz3xizOWZepf8sw5naGGPG\nGWNOGGO23PBcTmPMPGPMrpifOZyMManiOMa+xpgjaWGOLGNMYWPMImPMNmPMVmNM15jn08TnGM/x\nperPME21ARhj/IGdQC3gMLAWeCqtLT9pjNkPBIuItw1ASRRjTDUgHPhORErHPDcYOCMiA2MSeQ4R\n6eFknEkRxzH2BcJFZKiTsXmCMaYAUODGpV+BxkA70sDnGM/xtSQVf4Zp7QqgArBbRPaKSAQwEWjk\ncEzqNkRkKXDmlqcbARNi7k/A/rOlWnEcY5ohIqEisj7m/gXg2tKvaeJzjOf4UrW0lgDuBA7d8Pgw\naeBDioUA840x64wxHZ0OJpnki1lPGuAYkM/JYJLRq8aYTTFVRKmyeuRWtyz9muY+x1iWtk21n2Fa\nSwC+oqqIlAHqAq/EVC+kWWLrKdNOXeV1I4GiQBkgFPjE2XCSLr6lX9PC5xjL8aXqzzCtJYAjQOEb\nHheKeS5NEZEjMT9PYNdjqOBsRMnieEy967X61xMOx+NxInJcRKJFxAWMIZV/jnEs/ZpmPsfYji+1\nf4ZpLQGsBe41xtxtjEkHtAZmOByTRxljMsc0QmGMyQw8CWyJ/12p0gygbcz9tsBvDsaSLK6dGGM0\nIRV/jvEs/ZomPse4ji+1f4ZpqhcQQEw3rE8Bf2CciHzkcEgeZYwpyvVV2AKAn1L7MRpjfgYew06t\nexzoA0wHfgWCsFOCtxSRVNuIGscxPoatOhBgP9DphvryVMUYUxVYBmwGXDFPv4OtJ0/1n2M8x/cU\nqfgzTHMJQCmllHvSWhWQUkopN2kCUEopH6UJQCmlfJQmAKWU8lGaAJRSykdpAlBKKR+lCUAppXzU\n/wFI2pcHFaH4zgAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7f8bd470ca20>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(SiGlobalTot[sdir0]-SiGlobalTot[sdir0][0],'r-')\n",
    "plt.plot(SiGlobalTot[sdir1]-SiGlobalTot[sdir1][0],'b-')\n",
    "plt.plot(SiGlobalTot[sdir3]-SiGlobalTot[sdir3][0],'g-')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "anaconda-cloud": {},
  "kernelspec": {
   "display_name": "Python [conda env:python36]",
   "language": "python",
   "name": "conda-env-python36-py"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}