Skip to content

Instantly share code, notes, and snippets.

@adcroft
Created January 22, 2016 19:37
Show Gist options
  • Save adcroft/47ee2517303105150442 to your computer and use it in GitHub Desktop.
Save adcroft/47ee2517303105150442 to your computer and use it in GitHub Desktop.
Demo Gibraltar sections for salt.ipynb
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 35,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"import netCDF4\n",
"import numpy\n",
"import matplotlib.pyplot as plt\n",
"%matplotlib inline"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Note: The 2013 World Ocean Atlas has climatologies for several decades. The CM4_2000 tests are using year 2000 radiative forcing but are initialized from the 2005 version of WOA."
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"climatology_bounds,crs,depth,depth_bnds,lat,lat_bnds,lon,lon_bnds,ptemp_an,ptemp_dd,ptemp_gp,ptemp_mn,ptemp_oa,ptemp_sd,ptemp_se,time,\n",
"climatology_bounds,crs,depth,depth_bnds,lat,lat_bnds,lon,lon_bnds,s_an,s_dd,s_gp,s_mn,s_oa,s_sd,s_se,time,"
]
}
],
"source": [
"woa_t = netCDF4.Dataset('/archive/mdt/obs/NOAA-NODC/WOA13/v1.00a/woa13_1995-2004_ptemp_annual_01.nc')\n",
"for v in woa_t.variables:\n",
" print(v, end=',')\n",
"print()\n",
"woa_s = netCDF4.Dataset('/archive/mdt/obs/NOAA-NODC/WOA13/v1.00a/woa13_1995-2004_s_annual_01.nc')\n",
"for v in woa_s.variables:\n",
" print(v, end=',')\n",
"x_woa = woa_s.variables['lon']\n",
"y_woa = woa_s.variables['lat']\n",
"z_woa = woa_s.variables['depth']"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": [
"def find_index(coordinate,value):\n",
" \"\"\"Helper function that returns the index of a particular coordinate value.\"\"\"\n",
" return numpy.abs(coordinate[:]-value).argmin()"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Climatology"
]
},
{
"cell_type": "code",
"execution_count": 78,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"j,i0,i1,shape(q) 127 149 174 (102, 25)\n"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEfCAYAAAC04jrjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXmYXUW57//5Zk4YA8GQABKEcBQUES7gxYFGRSMo4AgO\ngIoeNUfkOl2C3qON/kT0XjioXDkKKOCAcBA4KEFBpBG8R8JsBBRRgiSQQJCEQObk/f1RtZOVnb12\nV+3uld179/t5nn567VpVtarWUG/V+1a9JTPDcRzHcRoxot0FcBzHcYYuLiQcx3GcUlxIOI7jOKW4\nkHAcx3FKcSHhOI7jlOJCwnEcxynFhYTjOI5TigsJx3EcpxQXEo7TD5LmSXp9i+leF4//KOm1g1+6\nDdf6mqRTq8q/k5B0u6R92l2ObmHYCAlJP5L0hKSlkv4s6eTCuWmSZkv6R4zzbUkj47lRks6Q9D8l\nzSykmSdpkaQJhbAPS7o5sTyfkHSnpJWSftDg/Esk/UbSEkl/kXRsyrl4vk/SCknL4t+DTcoxRtJF\nsT7PSrpH0oy6ODtIulrSczHee1LO1eUxPdb1hyn3Z4hh8a+VdOHA7KVm9tvBK9JGJO0EnAD8u6TT\nJc2uO/+XkrB3F35/QNJcSc/Hb+A7krYruV5f/FbGNClT0/eqv/cm4fso/Z6B/wN8uaxsTiZmNiz+\ngH2BcfH4n4AngFfE37OBHwBjgMnAH4BT4rm3AP8tHn8W2DYezwOeAk4vXOPDwM2J5XkbcAzwHeAH\ndedGAQ8B/wMQcDjwHLBXk3PTC+lvBj6UWI4JwJeAF8bfRwHPArsX4lwW/yYArwKWAPv0d67uOjcA\nvwUubfe70MK78wjwui2VroXrfA74bjw+ND4Dxd9TYjkeB0YUwtYDO8ffnwEWAm8ERgK7A9cBc4DR\nddeaBqwDFgPvbOG9qv1u+t40+z7i+Ubf8wHx9zjgaWByu9+dbvgbNiMJM7vfzFYWg4A94/E04HIz\nW21mi4BfArXh6mJgl3i8HbCikP7/AJ8t63H1U56rzew/CS9zPS8GppjZuRa4GfgdcGKTcyfU5aHE\nciw3szPM7O/x93WERuUAAElbAW8H/jXG/R3wn8AJzc5tUhDpeOAZ4KZm5ZI0S9LDsed5f4MR0jxJ\nn5F0XxxF/VTS2P7OxfPrJb2o8PtiSV9JvXaTMu8m6SpJT0paLOnbJfGKqqd5kj4by/qcpAslTZZ0\nfbz+jZK2T7l+ZAZwSzy+ExgN7B9/v4bQaXioLuyvZrZQ0rZAL/AJM7vBzNaZ2aPAuwnfxfvrrnUi\n8F/AJcBJZQVq8l4dmPLe9PN9lH3PL4rnVgJ3AW8qK5+TzrAREgBxCP088CChZ1Ubgp8LHC9pvKRd\ngDcTBAVm9ntgoqRPA782szWFLO8E+ggjjJaLlRhvBKH31EjtMQJ4aV3Y1yQ9Jek2SYclF0aaDOwN\n3B+D9gbWmtnDhWj3xbJMb3Kult+2wBnAp+i/rg8DrzazWpofSdq5cN6AdxE+/j2A/YAPFM43O1dP\nvQqp0bUnNyusgkryF4TGb3dCZ+KnTa5X5O3AGwj3962Ed3EWsBPheX6y2bXreBnwZwAzWw3cDtSe\n+WuBW4Hb4nEtrCZUDiX0vK/apLBmz8cyHVF3rROBHwE/Bt4k6QUpBax7r5q9U5slbZJn2fdMDHt5\nStmc5gwrIWFmM4GtCT2pq4HV8dSthEb2WeAx4I7Yi6mlu9jMzjGzW+qzBL4InCJpUqvFahD2Z+BJ\nSZ+TNFrSGwkf9oQm58YX0p9GaCinAt8Dfl7sRZchaTTh47/YzB6KwVsT7kuRZcA2/Zyr8RXgQjN7\nvKSuGzCzK81sYTy+AvgLcHBdtG+Z2UIzewb4ORt7x9bkXBkbGqDEa9dzMEF18zkzW2Fmq2KvuD8M\n+LaZPRXvy63A783sPjNbRXg3X5GQT43tCfe9xi1sFAivJqj5bi2EvYaNQmISsNjM1jfId2E8D4Ck\nVwMvBK4ws7uBvwLv7a9wDd6rlPemRuk70+R7ruWXMxpzShhWQgIgqmh+B+wKfFySCKOGnxEa4UnA\nDpK+npjf/YTe5CxaM25u1lOKo5VjCXrcJwi98CuA+Wa2tuxcIf0cM3vezNaY2aUEddSRAJLep40G\n7es2FEIaAfwQWAl8olCc54Bt64q4HeEjbHYOSfsDryeM1BrWdZMbIZ0YDZzPSHqGILh3rIu2sHC8\nAtiqybmtm10v4dr9Cf7dgEdLGtj+WFQ4XlH3eyUZZSeo8ooN7G+BV0uaCOxkZn8lqIgOjWH7xjgQ\n1KmT4vOvZwrB7lbjJOAGM/tH/H0ZTVROUPpeNX1v6rNoln/991w4tS3hvjgDZNgJiQKjCTrMHQgf\n+3mxUf0HcDGxUU3kS8BH2Gi7yKGhYDGzuWbWY2aTzOzNBPvJnP7O9Xsxsx+b2Tbx7yiAKCgvIqg6\n3mFm6wpJHgJGSdqrEPZy4I/9nAPoIei1/y7pCYKB9B2S7qwvl6TdCaOefwF2MLOJMZ9UddxmVa37\nvZzQCagxpRZnANd+DHhhVDsNlFbrCWGixT8Vfv+e0Oh+hNBBwMyeJahk/hl4PNodIAiPVcA7NimM\ntDXB1nFT/D2eYKc4LM4qeoIweeLlkvZrWKHy96q/96ZIaser9j3XeAlBheUMkGEhJCTtJOl4SVtL\nGinpTcDxwE1m9jRBp/zxeG57Qu8o+QWLPbXLgU3mqUfj6GbT9+K5kZLGEWYrjZQ0ttjYSHqZpHGS\nJkj6LGHW1cUJ57aT9KZ4fpSk9xGG479sUoXzCQbxo6O6o1i35wn66i/H672aoEP/oZktLzsXk3+X\n8OG+nKD6+XfCrJlGBsWtCA3CYmCEpA+yuZ2lEWWNa334vcD74n2fwUbVy0CufTthNHdWrP84SYcm\npBtsZrPRBoGZrSDYyz7NxhEDBLvEp9moasLMlhJsMN+O781oSdMIo9PH2PgsjwXWEhrfl8e/lxDU\nWCeWlKvhe9XsnarFafZ9NPue4/lxhIkXNybcO6cfhoWQIDQAHyO89P8AvgGcama/iOffTjBWP0XQ\nRa8iqHFy+DKhp1rs+exK+DAb8a+E3u1phBkkK4AvFM6fQOj5LSJMcz2iYDRvdm40wQ7wZKzPvwDH\n1BkJNxB70f9M+OgXFlRRxXnrMwk2jycJRsuPmdmD/Z2Levon498igpphRRTMm2BmDwBnE3q2CwmN\ndNm925CM8p5m/blTCQ3RMwQ9+tUDvDZRzfRWwtTkvxPer3c3TdQku7KyK6zhmdUk7aXAkbFxrHEL\noQdfrMetBBXaJus1zOx/A58nzNZbShiJPAq8vvBenQh838zm1z3T84D31qurEt6rZu8UNP8++vue\n30qYil5UPzotUptL3RHEHuC5hLncF5pZkt2gHSgsNLoH2K9OfeM4g46krwJPmtk3212WdiPp94R1\nQg+0uyzdQMcIiTjU/DNh2uAC4A7gPXW9D8dxHGcQ6SR108HAw2Y2Lw6Bf0pYkek4juNURCcJiV0I\nOsga82ltNpHjOI6TSCcJic7QizmO43QRo9pdgAwWENYz1NiNwgIyAEkuSBzHScLMBrI2Jbu9Gej1\n2kUnCYk7gelxDvfjwHHAZm6prbgqIWfN6g4ZcadmxN2q/ygb2G3zd663t5fe3t5NA3+X8a49nxhv\nVf9RNpDz1pQsM+v9IfTWuyQc1zBqcp4NySnr2sR4JXPVen8AvR9sMc/cuKkMwry63h9Bb72bv5UN\nozbmbUOv7xbW+Q2c/y8x3v8alKu1h44REma2VtIngF8RmomLGs5sKvpjzWmg650/NCPD5+uqjLhj\n+48S6NcLU4GlifGqaKCg/A3bjuAFKCVuPck3qiLK7tW2bDrWhbxGugqBMhjX3wqod+NX1fvSYYxu\ndwG2AB0jJADM7Hrg+qaRio19jgPvjJHEsy8o3WtlM5aNbOSzrDGpVvhnpozvP1Jkm61W9B8JGJnx\n0a/LeGvK4q7dDlbVCYlVY9Pu67qR6UOJdRnDjrWJcUeVtLzLJ65i8R6bSrCRGa30yHUZcdemPbCx\nq9LdSo0qG002EujPJWfb1XRUA9oi3VfHYo8nQ0jkNPxPj0wfdizfxF1QcxoJiZ6engHluW7bwXAr\ntCmpjSnA6pJu/8uOWM2iCWPq4qY9g1UZQ4kcIZETtxHTe57nsbrha5aQGJked+zINP3gyLHpeY7Z\nxInqRvY9ag0Ldtq0zzx2p3T9ZKvukTuB9O5a59Ixi+lSkGT25Mbfz+6Q3vAvGZnuVXhJhgfinAb9\nldyTFO9+9uo/UmRsorEhp+Ffl9G3WJ7xGZUJlM3jpT/XqgRKKllCIkOHUzaaqaes4W8cN73hH5uR\n735hq4shhaRBMVxfkBj3I7jheshQFAxVNfxPZxgwchqpKq6f00ilUlUjnSpQUoUJ5DX8VTyrHNot\nJFI7FLn5NnQR2yV0XQPagK6rY1EVtKzhHiaNyRESz2TEzWnQUskpa6qQqKoxzREoqaOudo8kqhhx\nQO6oIy1uTsM/njT7VW6+3YwbrjuQxYVedk5juoSJGXFzhER6g5ZKjpBKJUeFlNNIVqFuWpVxT3OE\ndOqzylHN5ZA6OoB0IVHVSKKKEWon0moDGj323kKYqzcKuNLMeuO5UwhectcB15nZaQ3Sbw9cyMYt\njT9kYavlQafrhMRDhb1XckYSizNUOHlCYvBHEg9tsr9MczpJ3dQphuuqRhI5pD7XnOefY5PIEWjd\nTKsjCTNbKelwM1suaRRwm6TrCdsNHE3wHr1G0k4lWXwTmG1m74zpcyb8Z9F1QqLYgOcIiecy4ubk\nW4WQyLl+FeT05HNGKOlCoprrD3chMTbjvvpIIjAQdVPctAtgTMyqtk/G12r7eJjZU/XpJG0HvMbM\nTopx1pK+IiqbrhMSRVVMTsOfo8LJGUnkNFJVXL8Km0RV9ovUfKsSEul5dpKQSDeG59TLhURgIFNg\n40ZNdxO2Hz7PzOZI2ht4raQzCevaP2tm9dv97gE8FXe9fDlwF2HTpeVUQNcJiacLs7JzetzLMnx4\n5NgvqmhQqhASObRbSFSxQC6HKgQPVDO7aWRFAi1nCmw3U/Ym3Bv/mhF3Ntw/jgyulrRvzHKimb1S\n0kGEbWTrfSyMImzP+gkzu0PSucAs4Iut1qMZXSckWlU35RhYc/KtQkjkXL+KkUQOg7HwbiAMhV5/\nOoOv7smxM+R1KFxIQLm66aD4V+OSJnmY2VJJNwMzCE5Lr4rhd0haL2nHui1/5wPzzeyO+PtKgpCo\nhK4TEsU1BDm92KpGHVX0OjtJSFTR66+qJ18FVYwOcqhu1OPqJhjQ7KZJwFozWyJpPHAEcBbB4cnr\ngFui6mlM/Z7wZrZQ0mOS9jazhwi7dd7fei2a0zlfWyKPFTysVbGQC2BFjluM9RmNb+LuHouYnJ6n\nk0xOg55KTsNf1XTVVHLqvw3LBv36ncgADNdTgEvitswjgMvNbLak0cD3Jc0lDNdOBJA0FbjAzI6K\n6U8BfixpDPBXoN7/8KDRdUKiaICuamXu8uUZQmJthpDYNi3aivXDwWPM4DByRM7ag7RnldOLbvdU\n0apGMjlqrG6m1QbUzOYS7Ar14WuAeif6mNnjwFGF3/exqUarMrpOSDiO42wpfMV1B/LYoo3qpvU5\nPq1zevw5cXNIHEksW5Juk1i7JtH99ej265hHjkqcsZMYL8RN70mPHZOm7qnCfQZU40JjQoarjRwV\n0mSe7D/SMGA4jOm7Tkg4juNsKXwk0YGsn1dYnT4UnmAFd3j1c+k2kdRRz+qM3nkOI8ZWYGDNKGvq\n6ABgPGlrkXJ63Ntk7M6Tk+/2PJMYb0lynjmjg5cxNzluN9N1DWgDuq+OxcliObVL3V8Z8rbPbPcd\nXpvqwj5DhVaRQBkzLq1BnzAhfWFpXoOeFnfHTV6y5uQ00juyODnuxMR8X5DR8O/GY8lx/xv1i4CH\nJ6NTv+8O3u613U2Y4zhOxzLKhUQHUuwMVuUHL0eNlTNCSWTM1uk96dWj0oY9IzI2uR47Pl2FMyGj\nrONHpBlZc0YHOb3zSYkjhJyRxGQWJcfNyTe1XjkqpJyRxB6PPJEclz3So3YaoztpQX+LdJ+QcBzH\n2UIkjyQ6mGFQxURy7kS6Vw4YlzPOTCvE+IzeeaqRd8y4dAPzhDHp0yonJBqDIV1/X5WeP7XXPZXH\nk/PMsQlMyihr6qgjZyQzZVV6vfhLetSuHkm0d8fbLYILCcdxnFYZBi1o91WxaC/IqV2ObjFjdJA3\nBTStwDk9+dQpoKnTPyF3qmZ6rz91xk4VuntIH0m8IKN3nmMTqGJ2U8692urJ9clxeTY9alfTfS3o\nZnRfFYs1yjEaZxi5t56U3vCNzVDjpO5AmNNIp/qkyslz64y4qY0ZVKNuqiJuFdNqIW91dPKK6+Xp\nefJ8etROnq0zqHRfC7oZw6CKjuM4FeGzmwaGpO8TPBc+aWYvi2E7AJcDuwPzgHeb2ZJ47nTgQ8A6\n4JNmdkMMPxC4mDA2mG1mp5Ze1Hs4m1DNXsiDv4oZ0nvdOb3zKkYyOXnmjSTS71Vq3DErk7PMw7+z\nwDDoZifuYNAyPyDstlRkFnCjme0N3BR/I2kf4Dhgn5jmO5Jqy4XPB042s+nAdEn1eTqO42x5xib+\ndTCVykEzu1XStLrgo4HD4vElQB9BUBwDXBb9qc+T9DBwiKRHgW3MbE5McylwLPDLhhddUzjO6e1k\nxM3ZI2LVyow3JNEl0/KMTY9Se5zVeTZNH3WkbrqT0+OuYiOfnDxz9l3IGXWNWZWWr3JMYjk2iZy4\n3UyLLaikccAtBBEyCrjSzHrjuVOAmQSNynVmdlpJHiOBOwlbmb61tZL0TzsGS5PNrDY9ZBFs2GZt\nKvD7Qrz5wC6EZn9+IXxBDG9Msf3KGWqn+2FjRYar7iw/R4ltfxU7qOWQc/084ZOWb06eOY10auNf\nlYO/icuXJscdmzppKT1LMiZCwT8y4nYzLbagZrZS0uFmtlzSKOA2SdcTWoGjgf3MbI2knZpkcyrw\nANX5lgDarFEzM5Nkg5rpz3o3Hh/UA/v1DGr2juN0Hn19ffT19Q1+xgMwXJtZbeg4hjB534CPAV+L\nGhXM7KlGaSXtChwJfBX4dOul6J92CIlFknaOm3lPgQ0TyRdAYYNq2JUwglgQj4vhC0pzP6Z34/HW\npI8mckYdi3NuW0bcqWnRnluf0XEYkdbrzdm+NYcqVFM5o4OcaaXJxuCs66erkMZWoe7JGUnkXL8q\ng3hF9PT00NPTs+H3GWecMTgZD6AFlTQCuBvYEzjPzOZI2ht4raQzCXf5s2bWyOXuvwGfI3mrstap\n2nDdiGuBk+LxScA1hfDjJY2RtAcwHZhjZguBZyUdEg3ZJxTSOI7jtI9RiX8NMLP1ZrY/oeN7iKR9\nY+yJZvZKghC4oj6dpLcQZozeA6TuBdAyVU+BvYxgpJ4k6THgi8BZwBWSTiZOgQUwswckXUHQsa0F\nZppZTRU1kzAFdjxhCmxjozWwyezEHNX9mv6jtJRvBeaDZ57aPj1yM41mgW0SRxxVkjrqyDGG59hP\nUm0SOdfPWsyWYRdL7vXnrIzO2ZHUdy8NlLSgff8IfymY2VJJNxNmdc4Hrorhd0haL2lHMytajA4F\njpZ0JGFZwLaSLjWzE1uuRxOqnt30npJTbyiJfyZwZoPwu4CXDWLRHMdxBk7J5MWeKeGvxhl/3fS8\npEnAWjNbImk8cAShA/0c8Drglqh6GlMnIDCzzwOfj/kcRlBJVSIgoBuXgqSvc9qUHBce8/uPsoGc\n8hyZFm39ojT3HQDLRqf1zrffIX1a7eqMid9V2TraSdaMrZyRZBXTVXPsDDnX7zCbRGW03oJOAS6J\n01hHAJeb2WxJo4HvS5oLrAZOBJA0FbjAzI5qkNfgTv6po/uERNFHWs6LnPMx56gFWhVazXgwPerq\nFWl2rcdfkmg1B8ZvO/grg3Pi5uSZs/ZgReLDWs2Y5DzXZXxhbd+bIOcbqGYH286jxT6Qmc0FDmgQ\nvoZgd60Pf5zgvaI+/BbCeovKaPdr6TiO07kMgxa0+6pY3FUxZ4lJzqgjx8abs0Aplb/2H2UDiT2+\nFaMnJmc578Xpu8iMmtA5Xc5RiTcrZ8X1+AnphuupO6S/LKpC3ZTTGviK60D3taCbMQyq6DiOUxHd\nZ3LbjO4TEsXOWE5vJ8d2kDOSqMIm8aeMuBVMwX2OSclx/7xX+le0fNs043mO76oq4q7NaBlyDPfr\ndkiPO3VU2hzUrA88Z9Tb4U7rBo3ua0E3o/uquLBwnDNjKadHUJXwSWVeBXlWtPZjxcp0NdZD09Ia\n6WVT0/WIS0hfU7IsUT+5nPHJea7IEFKp1wdYtu28pHjTXvxIcp5brczYme6F6VG7mpw2pkPpPiHh\nOI6zpXB1UwdSdIc1ujTW5uTciapUU6k8nBE3tddf1SryrOnCaTqMJxanG86X7ZXeO18+IVXdlT6S\nyBkd5Ix6UuMuGZue5z6HPJAcd+LjGSvJu5nua0E3YxhU0XEcpyKGQQvafVUs9vJzRhJV+W7K8EJK\nag+1ofPgElJ1pjlvQs69yplanBo3YyT33HPpRvYHp6X1+vNsIuk2maczJgQs2rANy+DEg7yRTM+R\nNyfHTfcP0IG4uslxHMcpZRi0oN1XxbUlx4NKTle6gpFEzggp1SayuP8oLZGzli61rDk2oQpsIvOX\n7JWc5bK9tk6O+/SYHdPjkhb3SV4w6HkCLB6bHvek/qN0Lt3Xgm5GF1axVYNaTsOf4385x3K9Q1q0\nFRllXZkjURLpILfqWVRw/aVrd06Ou2xyhhprcppqaHFOw1+Bugu6XEgMg/UiXSgkHMdxthDDoAXt\nwirm9PKL5IxAcq5RxVTBRelRLbHHtzhjxFGVO6bUEUpV7rdTH1WOMT5DNbZ+cbqJ98lEd/FPTkrv\n8T+xc7on4HkjpiXH7Wq6sAWtZxhU0XEcpyJ8dlMn0mrPPSddVXFTybFzpI4QEu0hAM9kjDqq6PXn\njGSq2FOkqs15cgzyqaaGielK86cX7pIeN2OE0tUuPLqwBa1nGFQxVTWU05jnNNJVCIkcdVfqI36W\n9JYnQ6C0e+vsKozhC0m/VTmPP2Pr8mSBkj65Ku/622c0HS4kNkPSOMJmQWNjLleaWW88dwowk9Ad\nus7MTqtLuxtwKfACwq503zOzb7VWkv7pciHRqn1iOJI+C2bY47fKqdH6znQrJR1uZssljQJuk3Q9\nMAE4GtjPzNZI2qlB8jXAp8zsXklbA3dJutHMMvasTKcLhUSx65bTjRxNuuU0p3uYvn1mOlWomyqY\nKgt01KgjlXbvWw3pI4mcjbeqGnV0MwPwAmtmtcZhDOEDNOBjwNfiNqaY2Wb+FcxsIdHftZk9J+lB\nYCpZGxun04VConXMDkmKJ/284pI4jtMRDKAFlTQCuBvYEzjPzOZI2ht4raQzCRatz5rZnU3ymAa8\nAri99ZI0pwuFRFHFVNXqrJzRQbtXiFUxOsp5bXLibpsWraoRRyd5zE3t9ecYw3NGHTll7WYGMLvJ\nzNYD+0vaDrha0r6ED2aimb1S0kHAFcCLGqWPqqYrgVPNrLIn0oVCovsxe0NyXKmyDobjOCUtaN/c\n8JeCmS2VdDMwA5gPXBXD75C0XtKOZrbJvoGSRgM/A35kZte0XP4EulBItGqTyCEn3xwfFlWQOkKo\nyiZRBYkjDsgbdVTxuuQ8/vRtKtLrlZNnjk2kU+xHVVPSgva8IvzVOOOnm56XNAlYa2ZLJI0HjgDO\nIozRXgfcElVPYxoICAEXAQ+Y2bmDVJNSKhUSZVO1JO0AXA7sTtiM891mtiSmOR34EGH61yfN7IYY\nfiBwMcFUNNvMTm181eKXXtXGKJ204Upqy1eVuikn3yoEVYZAqeKx5giJHCN36vKHKlacg6ubarSu\nbpoCXCJpJDACuNzMZscRwvclzQVWAycCSJoKXGBmRwGvAt4P/EHSPTG/083sl61XpJyqRxINp2oB\nHwRuNLNvSDoNmAXMkrQPcBywD7AL8GtJ083MgPOBk6NxZ7akGVXdFMdxnCRanN1kZnOBAxqErwFO\naBD+OHBUPL6NIFi2CJUKiZKpWrsQ5gEfFqNdAvQRBMUxwGXxRs2T9DBwiKRHgW3MbE5McylwLNCP\nkKhK1dNuFVIOqd3DnFehKstt6n3NGXFU4Ko9J8t2jyRyVEg5Dd4wcEeRxDC4D1vMJlE3VWuymdW8\n1C2CDX6HpwK/LySbTxAqa+JxjQUxfDOKRl3pqoEX3HEcp4wutOrWs0WqGFVNPyNM1VoW7C4BMzNJ\nNljX6u3tLfxaQ9BcDTY51sBdK7h+DrsnxsupU05PPucVq8ImUYWhIeNeVTUFNrVaOY81lw7al7Sv\nr4++vr7Bz3gYCAkFdX+FFwiGmF8A19cs8ZL+BPSY2UJJU4CbzezFkmYBmNlZMd4vgS8Bj8Y4L4nh\n7wEOM7OP1V3LivWRfpxVVrP3Jdbp4oxc04VEztTWVKQ/JMbsJCFR1Uys1HtQUcur/qNsIFXNkVPU\nnBXXgD2eF38oIQkzy7njjfIwuzcx7v4M+HrtolLjR5OpWteyccOqk4BrCuHHSxojaQ9gOjAn2jae\nlXRIzPOEQhrHcZz2MDLxr4OperDUcKoWYT7wFZJOJk6BBTCzByRdATxAGKjPLAwNZhKmwI4nTIHt\nkJlN7Z4uW8Uy4qpo93TdKvxsZYx6LMcFe2K9lmXkuaxK3VSXMgzUTVXPbmo2VauhbsXMzgTObBB+\nF/CywSud4zjOAPE9rh3HcZxShkELOgyq6DiOUxHDoAUdBlWsgm703VSVnr+qmVBVXL8KqvrEUutV\n5RTkjC1Mu5Vh0IJ2dRVTp7Q6juO0gnX4zKUUulpIOI7jVMm6YdCCDoMqtptOUTflUJWfpyquX0X9\n263CgqGhmnN1kwsJx3Ecp5RVY8ckxlxdaTmqxIWE4zhOi6wb2f1GCRcSjuM4LbKu031uJOBConLa\n7e6i3TaJbrQJDIXPpt33wAFY26KQkDQOuIWwZnsUcKWZ9cZzpxDcEK0DrjOz0xqknwGcS/AMdaGZ\nfb2lgiSpAfAeAAAfFElEQVQwFN52x3GcjmRdi02oma2UdLiZLZc0CrhN0vXABMKmbPuZ2RpJO9Wn\njVuenkdwbbQAuEPStWb2YMsVaYILCcdxnBYZiLrJzGoeJccQhoYGfAz4WtydEzN7qkHSg4GHzWwe\ngKSfEnb1dCHRmbR7CmybPZtWkm+7VXhD4bNp9z1wYGBCQtII4G5gT+A8M5sjaW/gtZLOBFYCnzWz\nO+uS7gI8Vvg9Hzik5YL0w1B42x3HcTqSVTSeAjunbyVz+lY2TWtm64H9JW0HXC1pX0KbPNHMXinp\nIOAK4EX1SQdc8AyShYSkrYCVZrauwvI4juN0DGU2iQN7tubAno1b/f3fM5aW5mFmSyXdDMwgjAqu\niuF3SFovaUcze7qQZAGwW+H3bjFdJZQKiWgcOQ54H3AQYTXIWEmLCduRftfMHq6qYN1Du9UC7b5+\nDqmquXbP7Gm3ChE667l2L62qmyRNAtaa2RJJ44EjCJuxPQe8Drglqp7G1AkIgDuB6ZKmAY8T2un3\ntFSQBJqNJH4D3ETYSe7+2ghC0o7A4cDXJV1jZj+sqnBDFbMPJ8eVLqywJI7jtJMB2CSmAJfEzvgI\n4HIzmy1pNPB9SXMJHfMTASRNBS4ws6PMbK2kTwC/IkyBvaiqmU0A2rg7aN0JaYyZNV1LLml0zQo/\nFJBkZfVpFzlCIkf4pF//x4OeZ/tNWe0eSQwF2n8PzN7a7iK0jCTMTAPMw+bYS5PiHqw/Dvh67aL0\nay8KCEkTgRdS2NLbzO4eSgLCcRxnS9PqOolOot8aSvoK8AHgb8D6wqnDKyqTM6h0oxzP0cd360fc\njc+183C3HIHjgD37Uz05jalCheQ4ztBgdckU2G4iRUj8EZgILKq4LI7jOB1Fq76bOokUIXEmcLek\n+4FVMczM7OjqiuUMHqmqmRy1TLtVHTlG23ZPFc25V+02Rrf7uXYebpMIXAp8nTCiqNkkhtYUIsdx\nnDYw1G0ScQ3GW4DXAFMJbpnnErzL3p+SR4qQeN7MvtVyKR3HcbqUoSwkJJ0BvBXoA34PPAWMA/YG\nzooC5NNm9odm+aQIiVslfQ24lo3qJszs7n4K2NBfuqQdgMuB3YF5wLvNbElMczrwIYIf9U+a2Q0x\n/EDg4ljB2WZ2akK5HcdxKmWI2yTmmNmXSs6dLWkym7r3aEiKkDiAoF56ZV140ymwTfylvwO40cy+\nIek0YBYwS9I+hJlU+xC8HP5a0vS4Ou584OToJXG2pBlm9suEsjvJVKW7r0Jn20m685z7WsXU3nbb\nZLqb1YxtdxFKMbPrACS9y8z+o3iuENbvhKR+3zQz6xlAIRv5Sz8aOCyGX0IYCs0i+EO/LC7Qmyfp\nYeAQSY8C25jZnJjmUuBYwIWE4zhtZSirmwp8HviPhLCGNHPwdwLw4+jOttH5vYApZnZrkzwa+Uuf\nbGY16bUImByPpxL0ZjXmE0YUa9jUw+GCGO44jtNWhrK6SdKbgSOBXSR9C6i5BdmGjOF4s5HEjsA9\nku4meB1cTLAJ7AW8FniaMAIopYG/9JfWnTdJPlOqUtrtWdXVHdVQxX3tJDXe0GCIT4F9HLiLoKW5\nqxC+DPhUaibNfDedK+k8gtvaVwH7EaZPPQicYGZ/T71IwV/6m4BFknY2s4WSpgBPxmj1PtJ3JYwg\nFsTjYviCsmv19vZuOO7p6aGnpye1mI7jdCl9fX309fUNer5DWd1kZvcB9yl4+RwBvJig8v9TjgeN\nUi+wA6WBv/RfEfyl9wBPm9nXJc0CtjezmuH6J4T9W3cBfg3sFUcbtwOfBOYA1wHfamS4HopeYNuN\ndH5izHYv5HLaT95IwuzjFZWjegbLC+z59oGkuB/XxW3zAivpKODfCf73IOx091Ezm52SvsqxUpm/\n9N8DV0g6mTgFFsDMHpB0BfAAYSw9s9DizyRMgR1PmALrRmvHcdrOUB5JFDgHOLy2SZykPYHZ8a9f\nKhMSZjaXMH22PvwfwBtK0pxJcANSH34X8LLBLqNTxPXRjtuPclk1hKfAFni2bhfRvwHPpiYeMfjl\ncRzHGR6sY2TSXz2Sxkm6XdK9kv4oqTeG90qaL+me+Dej0XUlfSqmmyvpJ5KaSau74vqyD0j6AGH7\n6TslvV3S2/urY8p+EuMIC+CmFeKbmX25v7SO4zjdTKvqpiaLjQ04x8zOKUsraRfgFOAlZrZK0uXA\n8YR1Z40YR5ggVFufVnPPUdta8KpmZU1RN/0nsIQwhWplQnxnSFGFCmFIT/sbJrhqaCgwkHUSJYuN\nYeN6hmaMAiZIWgdMoMmMT7NE63qTC/XHLmb2poFcxHEcpxsZyDqJksXGbwZOkXQiYX3aZ2q+7WqY\n2QJJZwN/JyxL+JWZ/brJdfYEzgX+O0EQ/T/gU2b2t7I0RVJq+P8k7defp0DHcZzhRpm6aV7fozza\n92jTtA0WG+9L8FNXU+V/BTgbOLmYTtJEgnujacBS4D8kvc/MflxyqZ8A5wE1+8NxwGXAIU0LGGnm\nlmNuPBwJfFDSI2y66dB+KRdwuhFXdXQOPmutSsqExG49L2K3nhdt+P3bM24rzaOw2HiGmZ1dC5d0\nIfDzBkneADxiZk/HeFcBhwJlQmK8mf2w8PtHkj5XWqA6mo0kakYNY3Mdma9Ycxxn2LOqxT2uGyw2\nPoKwx8POZrYwRnsbYYOgeh4FXhnTrSQIjTkN4tW4Pm7DcFn8fVwM2wE2LEsopZlbjnmxMj80sxPq\nKvhD4IRG6RzHcYYLA7BJlC02vlTS/oSO+CPARwEkTQUuMLOjou3iSoI9Y238/70m1zou5vfPJeEv\n2ixFgZQabuKUL07XOjAhneM4TlczgCmwZYuNTyyJ/zhwVOF3L9CbeK1prZSxRjObxOeB04HxkpYV\nTq2hudRyHGcDbhPoZoayWw5Jr+lnK4dtgd2jwCqlmbrpTOBMSWeZWVOX4I7jOMORobyfBPAOSd8A\nriesc3uK4P9uL4Kj1WnAZ/rLJEXddLqkdwCvBtYDt5nZ1a2V2XEcp3sYyvtJmNn/kLQjwWPGuwh2\nkNp2D981s/IpVwVSavgdwmKPywiznD4m6Qgzm9lSyZ0tjKs7uhOfhjwUGMrqJoA4TfZ7DMBEkCIk\nDgf2qW1jKuligjtvx3GcYc3qFqfAdhIpQuJh4IWEvR+Ixw+XxnYcxxkmDHGbxKCQIiS2BR6UNIcw\np/Zg4A5JPyesvD66ygI6juMMVYayTWKwSKnhF5uc85XXXUOOjrv7P4zBo5NsB26/ymUo2yQkHQw8\nZmZPxN8nEYzY84De/lZa1+j3azezPknTCPtN/1rSBGCUmSXvbOS0D7NPJ8ULM+Ucx8lhKAsJ4LvA\n6wEkvRY4C/gE8AqCIfudKZmkbDr0z8BHgB0Is5x2JXgqfH0rpXYcx+kWhrhNYkRhtHAcYdrrz4Cf\nSbovNZMUvcG/EOwQvwcws4ckvSC3tE430UkqlG7FVUNDgSFukxgpabSZrSE4ASz6bkoueErEVXGL\nPGCD7ya3RTiOM+wZ4lNgLwNukbQYWA7cCiBpOmG30SRShMQtkr5A2CrvCGAmjX2cO47jDCuGsrrJ\nzL4q6TfAzsANtbVuhEXRp6TmkyIkZhF2RppLcFs7G7gwr7iO4zjdxxBXN2Fm/9Ug7KGcPFJmN62T\ndA1wjZk9mZO54zg5tNvO4LamXIb47KZBYUTZCQV6oz7rz8CfJS2W9CXVDBSO4zjDmHWMTPrrZEqF\nBPAp4FXAQWY20cwmEmY5vSqecxzHGda0KiQkjZN0u6R7Jf1RUm8M75U0X9I98W9Go+tK2l7SlZIe\nlPSApFdWVcdm6qYTgSPM7KlagJn9TdL7gBuBc1IuELfnuxOYb2ZvjfuqXg7sTlj5924zWxLjng58\nCFgHfNLMbojhBwIXA+OA2WZ2ak4lnRTarepwqsFVSFWyirEtpTOzlZION7PlccbobZKuJ8wcPcfM\n+mtfv0loC98Z02/VUkESaDaSGFUUEDViWI615lSC19jatNlZwI1mtjdwU/yNpH0ICz72AWYA3ymo\ntc4HTjaz6cD0MunqOI6zJRmIusnMlsfDMcBoNraRTdX5krYDXmNm34/5rDWzpYNTo81pJiSadS2T\nup2SdgWOJMyGqlX8aOCSeHwJcGw8Pga4zMzWmNk8gqfZQyRNAbYxszkx3qWFNI7jOG1jIEJC0ghJ\n9wKLCFNUa23cKZLuk3SRpO0bJN0DeErSDyTdLemC6C6pEpqNCPar29u6yPjE/P8N+BzBk2yNyWa2\nKB4vAibH46nEVd2R+cAuBIE0vxC+IIY7zjDG1UhDgbJ1Eiv75rCyb07DczXiuoX948jgakn7ErQm\nX45RvgKcTViCUGQUcADwCTO7Q9K5BI1MM2esLdNsj+sBmeQlvQV40szukdRTcg2TNKirt3t7ezcc\n9/T00NPT8NKO4wwj+vr66OvrG/R8y9ZJjO45lNE9h274vfSM/1uah5ktlXQzMMPMzq6FS7qQxguX\n5xNsvHfE31cS1fZVUOVKkEOBoyUdSTA4byvph8AiSTub2cKoSqqtvVgA7FZIvyvhZiyIx8XwBWUX\nLQoJx3Ec2LzDeMYZZwxKvq1Ob5U0CVhrZkskjQeOAM6qtY0x2tsIi5g3Ibadj0naOy6MewNwf2s1\n6J9mNokBYWafN7PdzGwP4HjgN2Z2AnAtcFKMdhJwTTy+Fjhe0hhJewDTgTnxhj0r6ZBoyD6hkMZx\nHKdtDMAmMQX4TfTGOodgk5gNfEPSH2L4YcTlBpKmSrqukP4U4Mcx3n7AmVXVcUuuKa+plc4CrpB0\nMnEKLICZPSDpCsJMqLXATDOrpZlJmAI7njDt65dbsNyOs4VwO0OnsWp1aw7+zGwuwa5QH35iSfzH\ngaMKv+8DDmrp4ploYzvc+UiybqrPlkT6aruL4HSYkDD7UruL0DKSMLMBeY6QZFs/v9kqgYY8t9VO\nA75euxja3qkcx3GGMOvWdrbLjRRcSDhDiE7pSQ+Fzya1DJ1yTzsTFxKO4zhOKWvXuJBwHMdxSli/\nrvub0O6voeM4TlW4uslxGjHc9dxDof5DoQwOK7u/Ce3+GjqO41TFMJDVLiQcx3FaxYWEM3wYBm+7\n4ww2w+CzcSHhVE7qylxpcJyuOc4WYxhs6OhCwnEcp1XWtbsA1eNCwnEcp1Vc3eQ4juOUsrLdBage\nFxKOM4Rw+02HMQxGEpVtOuQ4jtP1rE38q0PSOEm3S7pX0h8l9cbwXknzJd0T/2aUXVrSyBin0Ran\ng4aPJBzHcVqlxZGEma2UdLiZLZc0CrhN0vWEzdnOMbNzErI5lbBJ2zatlSINFxJOS3TyhjOOM2gM\nYAqsmS2Ph2OA0WzcvbPfzYkk7QocCXwV+HTrpegfVzc5juO0yrrEvwZIGiHpXmARYY/rOfHUKZLu\nk3SRpO1LrvxvwOeA9YNUk1JcSDiO47RKizYJADNbb2b7A7sCh0jaFzgf2APYH3gCOLs+naS3AE+a\n2T0kjDoGiqubHKC71Uc+Y8ipjLIpsH/ug4f6krIws6WSbgZmmNkGoSDpQqCRUfpQ4GhJRwLjgG0l\nXWpmJ+YUPRUXEo7jOK1SZrjesyf81fjFph0QSZOAtWa2RNJ44AjgLEk7m9nCGO1twNz6rM3s88Dn\nYz6HAZ+tSkCACwnHcZzWaX2dxBTgEkkjCWr/y81stqRLJe1PMGI/AnwUQNJU4AIzO6pBXtYgbNBw\nIeE4jtMqrU+BnQsc0CC84YjAzB4HNhMQZnYLcEtrpUjDhYTjRKqyXXSzvWfY415gHcdxnFKGgRfY\nyqfASpon6Q9x+ficGLaDpBslPSTphuJcYEmnS/qLpD9JemMh/EBJc+O5b1ZdbsdxnH5ZmfjXwWyJ\ndRIG9JjZK8zs4Bg2C7jRzPYGboq/kbQPcBywDzAD+I6k2jzg84GTzWw6ML2ZTxPHcZwtwgDWSXQK\nW2oxXf2Cj6OBS+LxJcCx8fgY4DIzW2Nm84CHCYtMpgDbFFYkXlpI4ziO0x7WJP51MFtqJPFrSXdK\n+kgMm2xmi+LxImByPJ4KzC+knQ/s0iB8QQx3HMdpHwNwy9EpbAnD9avM7AlJOwE3SvpT8aSZmaRB\nm+fb29u74binp4eenp7BytpxnA6lr6+Pvr6+wc+4w1VJKVQuJMzsifj/KUlXAwcDi2orC6Mq6ckY\nfQGwWyH5roQRxIJ4XAxf0Oh6RSHhdC/tnFbqU1o7j/oO4xlnDJILlmEgJCpVN0maIGmbeLwV8EbC\nMvNrgZNitJOAa+LxtcDxksZI2gOYDsyJy9SflXRINGSfUEjjOI7THoaBTaLqkcRk4Oo4QWkU8GMz\nu0HSncAVkk4G5gHvBjCzByRdQdhIYy0w08xqqqiZwMXAeGC2mf2y4rI7juM0Z1W7C1A9lQoJM3uE\n4PK2PvwfwBtK0pwJnNkg/C7gZYNdRsdxnJYZBuomX3HtOI7TKh2uSkrBhYTjOE6rdPj01hRcSDiO\n47TKMFA3aaNduPORZN1UH8dxqkESZjagrT8lGW9ObG+uH/j12oWPJBzHcVplGNgktpTvJsdxnO5j\nVeJfHZLGSbpd0r2S/iipN4b3SpofvWbf08iRqaTdJN0s6f6Y9pNVVQ98JOE4jtM6re9Mt1LS4Wa2\nXNIo4DZJ1xN83Z1jZuc0Sb4G+JSZ3Stpa+AuSTea2YOtlaY5LiQcx3FaZQDqJjNbHg/HAKPZuFd1\nU9tF9ECxMB4/J+lBghPUSoSEq5scx3FaZQBeYCWNkHQvwRP2DYWtEE6RdJ+ki4obspXkMQ14BXD7\nQKtShgsJx3GcVhnApkNmtt7M9ic4LD1E0r6EzdX2IHiqeAI4u+zSUdV0JXCqmT03OBXaHFc3OY7j\ntEqZTeL5Pljel5SFmS2VdDMww8w2CAVJFwI/b5RG0mjgZ8CPzKxSZ6e+TsJxnGHHoK2T2CuxvXl4\n0+tJmgSsNbMlksYDvwLOAu6ONgckfQo4yMzeW3ddEXb0fNrMPjWQOqTgIwnHcZxWad0L7BTgEkkj\nCWr/y81stqRLJe1PMGI/AnwUQNJU4AIzOwp4FfB+4A+S7on5nV6VZ2wfSTiOM+wYtJHElMT25glf\nce04jjP8GAYrrl1IOI7jtIp7gXUcx3FKGQZeYF1IOI7jtIoLCcdxHKcUt0k4juM4pQyDkYS75XAc\nx3FKcSHhOI7jlOJCwnEcxynFbRKO4zgt0/2WaxcSjuM4LdP9luvK1U2Stpd0paQHJT0g6RBJO0i6\nUdJDkm4obqwh6XRJf5H0J0lvLIQfKGluPPfNqsvtOI7TP2sS/zqXLWGT+CYw28xeAuwH/AmYBdxo\nZnsDN8XfSNoHOA7YB5gBfCe6xYWwGcfJZjYdmN5og3DHcZwty4rEv86lUiEhaTvgNWb2fQAzW2tm\nS4GjCf7Qif+PjcfHAJeZ2Rozmwc8TNixaQqwTWF7v0sLaRzHcdqEjyQGyh7AU5J+IOluSRdI2gqY\nbGaLYpxFwOR4PBWYX0g/H9ilQfiCGO44jtNGBrB/aYdQtZAYBRwAfMfMDgCeJ6qWasQNIHwTCMdx\nOpDuH0lUPbtpPjDfzO6Iv68ETgcWStrZzBZGVdKT8fwCYLdC+l1jHgvicTF8QaML9vb2bjju6emh\np6dn4LVwHKej6evro6+vr4KcWxslSBoH3AKMJbTDV5pZr6Re4MPAUzFqwx3nok32XGAkcKGZfb2l\ngqSUteqd3CT9FviwmT0Ub8CEeOppM/u6pFnA9mY2KxqufwIcTFAn/RrYy8xM0u3AJ4E5wHXAt+pv\nnu9M5zhOCoO2Mx33JcZ++WbXkzTBzJZLGgXcBpxKmLCzzMzOaXLdkcCfgTcQOst3AO8xswdbqUd/\nbIl1EqcAP5Y0Bvgr8EGC9LtC0snAPODdAGb2gKQrgAcIInpmodWfCVwMjCfMlqpkP1fHcZx0Wp+5\nZGbL4+EYYDQb1e79Ca+DgYfj5B4k/ZQw6aczhYSZ3Qcc1ODUG0rinwmc2SD8LuBlg1s6x3GcgdC6\nUVrSCOBuYE/gPDObI+nNwCmSTgTuBD5jZkvqku4CPFb4PR84pOWC9IOvuHYcx2mZMqP0vfGvHDNb\nD+wflwpcLWlfwnqwL8coXwHOBk6uT9pycVvAhYTjOE7LlI0kXhr/alxSEg/MbKmkm4EZZnZ2LVzS\nhcDPGySpn+CzG5suERhU3Aus4zhOy7Q2BVbSpJo7IknjgSOAByXtXIj2NmBug4veSfA6MS3aeo8D\nrh2sGtXjIwnHcZyWadkmMQW4JM5UGgFcbmazJV0qaX+CSukR4KMAkqYCF5jZUWa2VtIngF8RJgFd\nVNXMJtgCU2C3JD4F1nGcFAZvCuzPEmO/Y8DXaxc+knAcx2mZznbel4ILCcdxnJbpbJcbKbiQcBzH\naZnOdt6XggsJx3Gclun+kYRPgR3iVOOUrP10Y726sU7QvfUaHNxVuNNmuvUD7cZ6dWOdoHvrNTi4\nq3DHcRynlM4eJaTgQsJxHKdlun8KbNctpmt3GRzH6QwGZzHdlrteu+gqIeE4juMMLm64dhzHcUpx\nIeE4juOU0vFCQtJXJN0n6R5Jv5I0pXDudEl/kfQnSW9sZzlzkfS/JT0Y63ZV3JiE6B54RazvPZK+\n0+6yplJWp3iuk5/VuyTdL2mdpAMK4Z38rBrWKZ7r2GdVRFKvpPmF5zOj3WUakphZR/8B2xSOTwHO\nj8f7ELaGGg1MAx4GRrS7vBn1OqJWXuAs4Kx4PA2Y2+7yDXKdOv1ZvRjYG7gZOKAQ3snPqqxOHf2s\n6ur4JeDT7S7HUP/r+JGEmS0r/NwaWB+PjwEuM7M1FjYMf5iwgXhHYGY3WtjeEOB2YNd2lmcwaFKn\nTn9WfzKzh9pdjsGkSZ06+lk1oCNnHG1JOl5IAEj6qqS/A+8FvhiDp7Lpln7zCRuIdyIfAmYXfu8R\nh8d9kl7drkINkGKduulZ1dMNz6pItz2rU6L686LaTnHOpnTEYjpJNwI7Nzj1eTP7uZl9AfiCpFkE\nlVNvSVZDar5vf/WKcb4ArDazn8RzjwO7mdkzUVd8jaR960ZUbaPFOjWi455VAzr+WSUypJ5VkSZ1\n/AJwPvDl+PsrwNnAyVuoaB1DRwgJMzsiMepPgOsIQqJ+s/BdY9iQob96SfoAcCTw+kKa1cDqeHy3\npL8C04G7qytpOq3UiS54ViVpOvpZlTDkn1WR1DpKuhDIEYzDho5XN0maXvh5DFDb6/Va4HhJYyTt\nQfg452zp8rVKnGnxOeAYM1tZCJ8U98VF0osI9fpbe0qZR1md6PBnVccGHXcnP6s6inr7rnlWxZmQ\nwNuAue0qy1CmI0YS/fA1Sf9EMFjPAz4GYGYPSLoCeIDghWumxSkNHcK3gTHAjZIA/svMZgKHAWdI\nWkOo80fNbEn7iplFwzp1+rOS9DbgW8Ak4DpJ95jZm+ngZ1VWp05/VnV8XdL+BHXZI8BH21yeIYm7\n5XAcx3FK6Xh1k+M4jlMdLiQcx3GcUlxIOI7jOKW4kHAcx3FKcSHhOI7jlOJCwnEcxynFhYTjOI5T\nigsJx3EcpxQXEk4lSHqu4vx/F/9vJ+njLaTvlfSZzDTjozfXzdxLx8157in8nRbDx0i6peaew3E6\njW5wy+EMTSpdym9mr4qHE4GZBI+eWVm0cNkPAT8DdpZ0PLAK6DOzB4DlZvaKBuVcLekm4DiCA0rH\n6Sh8JOFsMSR9WtLc+HdqDJsWtzT9nqQ/KmxBO66Q5l/jNpm3Srqs1vsvjFTOAvaMvfevS9pd0txC\n+s9K+lI8/oKkP0u6FfinQpz3S7o95vHvksq+i/cC/0nYZOcS4LvAQQlVvwZ4X+JtcpwhhQsJZ4sg\n6UDgA4QG9pXAR6JzNYC9gPPM7KXAEuAdMc1BwNuB/YA3AweycQRQ+38a8Fcze4WZncbmO41Z4frH\nAS8nuCo/CDBJLwHeDRwaRwLradCgSxoDvMjM/k5wR/8W4P2EkQXA+Dp107sKye8nTZg4zpDD1U3O\nluLVwFVmtgJA0lXAawiupx8xsz/EeHcR9k4GeBVwTW1fBkmN/P2nbj9Zu/5KYKWka2Pa1xGEz53R\n1DAeWNgg/SSCAMPM1gKX1p1f0UjdFOOvk7Ra0lZm9nxieR1nSOBCwtlSGJs26GLjaGBVIXwdoaEu\nS9Mfa9l0hDy+cNwoLwGXmNnn+8l3BTCunzjNGAus7DeW4wwxXN3kbCluA46NM4S2Ao4FbqV5w/87\n4K2SxkraGjiqQZxlwDaF34uAF0jaQdJYglrIgN/G64+TtE0h/CbgnZJ2AojpXlh/ETN7BhgZ1U5Z\nSNoRWGxm63LTOk678ZGEUxUTJD1W+H0OcDEbdzG7wMzukzSNzWcaGYCZ3RnVQn8gNP5zgaV1cZ6W\n9LtorJ5tZqdJ+nK8zgLC5jiY2T2SLgfuA56slcPMHpT0v4AbosF6DWG21N8b1OkGgorspgbnxku6\np/D7+sLo5HDgFw3SOM6QxzcdcoY0NT2+pAnALcBHzOzeNpXlFcCnzOzEzHQ/A04zs4erKZnjVIeP\nJJyhzvck7UOwB1zcLgEBG0YjN0saYWbrU9JIGk0wvruAcDoSH0k4juM4pbjh2nEcxynFhYTjOI5T\nigsJx3EcpxQXEo7jOE4pLiQcx3GcUlxIOI7jOKW4kHAcx3FKcSHhOI7jlPL/A5RLdvYKyZI1AAAA\nAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f35d0f08a20>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"j = find_index(y_woa, 38.0)\n",
"i0 = find_index(x_woa, -30.0);\n",
"i1 = find_index(x_woa, -5.0)\n",
"q = woa_s.variables['s_an'][0,:,j,i0:i1]\n",
"print('j,i0,i1,shape(q)',j,i0,i1,q.shape)\n",
"plt.pcolormesh(x_woa[i0:i1],z_woa[:],q);\n",
"plt.gca().invert_yaxis();\n",
"cb = plt.colorbar(); cb.set_label('S (ppt)')\n",
"plt.xlabel('Longitude ($\\degree$E)');\n",
"plt.ylabel('Depth (m)')\n",
"plt.title('38$\\degree$N, 1995-2004 annual clim. (WOA 2013)')\n",
"plt.clim(34.9, 36.7); plt.xlim(-30., -5.); plt.ylim(6500, 0);"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## Model results"
]
},
{
"cell_type": "code",
"execution_count": 80,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"xh,yh,zt,zw,time,nv,xq,yq,thetao,so,age,uo,vo,uh,vh,thetao_xyave,so_xyave,time_bounds,average_T1,average_T2,average_DT,"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEfCAYAAAC04jrjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmcXUW177+/dKczMEMQEkCZggqIAR7DA5FGQaMoOFwF\nB0DFMVfkOV1A75VGFMH34OLIVUAJKgiXSZCgIKYRvEIYAoYZhAgJJkwShoydrPdH1emuPr33ObtP\n9+nT5/T6fj77c+rUtKvOUGvXqlWrZGY4juM4ThbjGt0Ax3EcZ/TiQsJxHMfJxYWE4ziOk4sLCcdx\nHCcXFxKO4zhOLi4kHMdxnFxcSDiO4zi5uJBwHMdxcnEh4TiDQNJCSW9t4P3vlfTmIdbR0D44zYUL\niQpI+qWkf0haJukhSccmadtKmiPp+ZjnB5LaYlq7pFMk/ZukWUmZhZKWSpqcxH1S0tyR7Vkfkjok\nnR/b9qKk+ZJmluX5vKQ7JK2U9PNB1D09lvlFWfymkq6U9HK874cG2eYPx/a8JOmp+D3sH9MWSlol\nabOyMvMlrZP06qLtzMHi1RDMbFcz+9NQq6GBfaiGpG5JK+L3+5KkBxrdprGMC4nKfAfYzsw2Ag4D\nviVp95j2Y2ApsCUwAzgQKAmEmcA1ZvZdYLKkDZM6xwHH17PRktoHkb0deAJ4s5ltCPw7cKmk1yR5\nFgOnAj8bZFN+BMxj4ID0I2Al8CrgI8A5knYuUqGkLwH/CXwrlt8m1vfuJNtjwIeSMm8AJmW0o1o7\nncZgwL+a2Qbxen2jGzSWcSFRATO7z8xWplHADjG8LXCJma02s6XA74DSQPcssFUMbwSsSMr/P+Ar\nkjYaTFskfVXSZWVx35d0dgwvjDOXvwIvSWqTdIKkRXGG8KCkt2T0cbmZnWJmT8T31wKPA3skea40\ns98Azw2ivUcC/wRuBJTErwe8D/iPeO8/A78BjipQ50bAKcAsM7vKzFaY2Vozu9bMTiw1F/glcHRS\n9BjgwrQd1dpZsI+vj0+9/4xqoHcnaR+XdHXy/hFJlybvn5S0WwwvlPQVSffE2dV5kraQdF387m6Q\ntHGS960xvIOk50oPLpKmSXqmoDpq93i/FyT9WtKEAr+xbSRdIelpSc9K+kHBz6Jw/9JbF+iDMxKY\nmV8VLsKM4RVgHXAHMDnGfxq4gPCEuhWwADg8Kfcx4EvAgUnc48BbgcuBU2PcJ4G5BdqxJfAysFF8\n306Yyewe3y8E7optmQC8ljBD2DKmvxrYvsB9tiAItZ0y0r4F/LxAHRsCDwHTgC7gF0na7sArZfm/\nBFxdoN6ZwBpgXIU8pc/4QeB1QBvwZOz/OuDVRdpZpf63AOOBR4ET43dxEPBi6XMDtgP+GcPT4vfz\nRHy/PfB8WZ3/A2we8y4F7gTeGL/LG4FvpPdPyn4SuC/+Dn8PfLdAHxYCt8bf1CbA/cBnKv3G4ud4\nD3BmvNcEYP+YL++zmJ7cr1D/Yv65wNPAM8AtJP8hv0b+8plEFcxsFrA+cABwJbA6Jt0M7Er4MzwJ\n3G7habtU7gIzO8vMbiqvEvgGcJykKYNoxxLgT8AHYtRM4Fkzm5/U+30zW2xmq4C1hD/gLpLGm9kT\nZvZYpXtIGg/8CrjAzB7OakbB5p4KnGdmT2WUWZ/wmaW8BGxQoN7NCH1eVyDvLwiziUMIg+DiQbaz\nGvsC65nZ6WbWY2Zzgd8S1Vxm9jhhRrc78GbCAP6UpNcSVJPl6wo/MLNnYltuBm41s3vid3klYaAe\ngJmdRxig5xEE/NcLtL30W1liZv8ErgFm5PzGnom/sb2BqcBXLczgVlmYBVb6LD6c3G8w/TuBIGSn\nAT8FrpG0fYF+OXXAhUQBLPBnYGvgc5JEUC9dDkwGpgCbSjqjYH33Ef5EJzK4wWk28NEY/ihBhZLy\nZHKPR4H/Q3hCXirpYklT8yqWNI4wsK4EPp+XLaPcR5IFxmslzSA8yZ+dU+ZlwhN8ykYEQVGN54Ap\nsa2VMEJfPkKOqqlAO6sxjeTzjvydPjUjwE1AJ+EB46Z4HUgQGuUPD0uT8Iqy9ysJwjWP84BdCAPx\nmmLNZ0nZ/Ur1l//GSov52wB/zxHQeZ/FtOR94f6Z2Twze8XM1pjZhcCfgXdW7ZFTF1xIDI7xBFXB\npoQ/zQ/jD/l5guppMD/kk4FP0X9QqcZvgN0k7QocSnjqT+kncMzsYjM7AHhNTMsUYlHonU9QB7zf\nzNbm3H+AQDOzX1nfAuOhhEFwW+AJSf8Avgy8X9IdscjDQLukHZNq3gjcm3PPlL8Aq4D3VstoYY3l\nMeAdwBUZWTqrtLMai4Ft4mdX4jXAouT9TQTVywFAN31C40AGColyyoVW5sOEpPUJgu484BRJmxRr\nfm7deb+xJ4FXK1rwlfEU2Z9F1uytt+k1tNNpAC4kcpC0uaQjJa0fF4HfDhwJ3GhmzxF0w5+LaRsT\nnljvKVq/mf0NuIQySydJFyjHzNTCIvplwEXAbWa2KCtfrGcnSW+RNIEwsK4kqKCyOIegvz8sTv/L\n62qTNJGgb26Li5xZgwUE9cD2hIF/BvBfwLXA22MfXiEM2t+UNFnSmwiWSb8o0P9lBFXdjyQdHsuP\nl/SOnFncsQT9/YqMtJ9UamcBbgOWA/8W29AJvAv4dZKnJCQmRjXLLQQVzqbAfAZH3qD6PWCemX06\ntv+/Bllvv7or/MZuA/4BnB4/94mS9kvSqn0WxRoibSTp7bH+dkkfIQjZ39XQL2cYcCGRjwGfJTxB\nPQ98FzjezH4b099HeEp9BniEMBB/cZD3+CZBXZU+yW1NGEzymE1YC6lm0z+BYML7DOHPPQU4qTyT\ngqnrpwmD5ZJEdZTuXfgPwiBwAkEFsYIc3XfUVz8dr6UE9dKKKFhLzCIsfj5NsET6rJmVbOEr9t/M\nziIsdP97LP9ErO/KjLyPmdldadQg25lLVOu8m77fwA+Bo9K1HDN7hKBGuzm+fxH4G/BnM6umZrSy\n8ID8kg4H3gZ8LkZ9CdhDg9x3klH/gN9YVDO9G9iR8Jk/CXwwpq2mymeRc8+s+3cQ1opKC9f/SjAI\neXSQfXKGCVX/rY4eFDZ5nU2wtDjPzAqtATQLkjoIT5i75al8FDaDPQBsYWYvj2T76k2R/jv1p5V/\nY87gaRohEdUbDwEHE3SdtwMfSp5AW564YHsWsL6ZfbLR7XFaD/+NOeUMZmduo9kbeNTMFgJI+jVw\nOOGJp+VR2IS2lLAWMrNKdqdG4lP0fRlJBuxcaR1otFBrH/w35mTRTEJiK/qb2S0C9mlQW0acuOBb\nyQzSGQaiVVSRPRujllr74L8xJ4tmWrhuDr2Y4zhOC9FMM4nFhL0JJbahv006klyQOI5TCDMb0l6N\nwY43Q71fo2gmIXEHMF3StoTNO0eQePosYWdWqCHtbVtGfF56Vvlq8dU+2bxdBmXluy6CrtRVXVv/\n9AHhau2eGF/XS+KSsE0Ir2urtL9aek9b3yS1fe3ATbqnfhv+IxrRtvX0xZfCSm2bkvSKcUXS12ak\n9wwivacsXxLXNRu6jqlQZxpfLX1lTnopPt3Nkta1KiMu61559Wd8bl1/ga7dM/JmkfX7TEnjz2nc\nM13/fX+1862C+f59WO7WGJpGSJhZj6TPE3zgtAHnZ1o2VepRlmBIwxMK1pOXXm3gr6XecWS3u1pb\n8+5TqiunfGnwT4VAOuD35avckbVtfY1enZm+ilUT2ge0pa2n2ugfyBI8eQxKCFUTSJUG2UkEByOV\n6iwqpCZkxEGfQK8mRPL6VxIi6UPCyoz0tEw7fQ8X1YRclhBJfyo9NNGoU53xjW7ACNBUX5eZXQdc\nVzFTpSfpwcwkqj2d55FVVxYFZxL9/qBpuazZQV69ad44+FgySCxfr08IrJ4QMqyiozdubZXOrM24\naVZcyspxy3ihLWONtK30UvmRta3qVKKP9qSuUr1ta5O4RDBlCZ+OZBDtFS7pIFsqPomwlxryn/SL\nzlSy6k/Dtc4kVpW9QrZASgVHR9KeVRnpKVm//wllec5uHa1wUw2gNdK6faz2pJ8lEAYjGGp5eq+h\nrs49yG53nmCo1pdYbnVSviQYAJYTDs1b3U9I9FXQE8N5gqOacCixa2cbL8R7tfUbxAcO/u1VdRwD\nacsQDGm4LZletLdl5y2xakJf3IRVYV6UNTs58CD6bIPKn55LVBMSWYN83u82K64UzqqzCBmyt3N7\nBg70eWVK989RZ7YakxrdgBGgdYVEi9CZ6SC6+dmrc3L1TE1G55uorrNvQjpfUz3PWMXVTc1IJXVP\nnjqpFK72xJZXV4m8h+ihzjqyVEtpvvVy8mbE9cS8yyf3PQO9lJjUl2YQq5JHx3RWUYpPZwxZs4c8\ndVHRmUbeTKAvrrK6KUvFVKSurFlLWzLT6Ji8ekCZCcmqS0mNlaqwJqzqU2H1romkT/ep6qZc918e\nzlqzqDZTyVIt5d0/S92U3quauqm0JrNpEjcxK2Nr0HoD6EDGQh8dx3Hqgs8kmpEs3WnRhem89CwG\noyfOK1e0fJYlUhqXPKlZRv9TS6WXNgwziHT2sDzRrK6I6wTpTCJ9+i/NKnr6zSQG/zOqNqPIm4nU\nsmaRPXvIm2n0DKgzK++E5FF8RZoeZx0dbUl6sqbRsba0plFlplHN+mowJrZZi9x5C98VTHz7lUvj\n0t/cZvE1PcE9b9bRAtQ6gEbX+zcRPr124DIz64ppxxE8G68FrjWzEzLKb0zfYVMGfMLMbq2xORVp\nPSGRNbWtpk4aqrqpFoFRTYhUERJ51kkl09TUBDUdkEsL06mQyFY3ZVs3leoqqjbKoyenfGlwTutP\nB+lSfP+49iTvQCGS3qta/W29/cteTC+pllIhmqWuyqt/VVsoly6cp0KkFuurrEX0lEyBU83sNyUt\nX+0hpxROBUNxQ7Smo9aZhJmtlHSQmS2X1A7cIuk6wtEBhxE8Ia+RtHlOFd8D5pjZv8TydTMPaD0h\n0Uj803ScMcVQ1E1mtjwGO2JVpTNsvlM6htbMnikvJ2kj4AAzOybm6QGWDaEpFWkm302O4zijikkF\nrywkjZN0N8Hz7vVmNg/YCXizpFsldUv6XxlFtwOekfRzSXdJOldS3cwF/dnXcRynRvIG0LvjVYl4\n2t+MODO4UtIuscpNzGxfSXsBlxKO2S2/7R7A583sdklnAycSjvYddlpPSGyUEVdUj1ptzSKv/GD8\nKRU0zbWkTLrxbdWEsFawuq1vzSBdU0jXEkqkOvvSIvXL/Rau+x5CSuUHsxhd6yJ0X3qf0rq0ftCe\nsQ5RpK56UMtmvlrXbHr7n27wS9eX1mYswqcbA9urrFlkuGVRle5ZRlfyfHe1vxIDrySRdVOENJ48\nddNe8Soxu0IdZrZM0lzCGR6LCGfAEwXAOkmblR2ruwhYZGa3x/eXEYREXXB1k+M4To20F7zKkTQl\nWighaRJwCOEAtauAt8T4nYCO8nPXzWwJ8GRMh3BaZ9YhU8NC680kHMdxRoghLFxPBWbHY5nHAZeY\n2RxJ44GfSVpA8I15NICkacC5ZnZoLH8c8Kt4LvzfgI/X3pTKNM0Z10WQZLagQoY8ddCEKukZ9CTT\n9VUTwoQs9YGU5e8oxLcPSE8pxeeZoGbtiF6RLI2lqqOstizv3QdR3IFfSn32KQw0IR3yLuncfRbF\n2tJvF3UN+ygG079+bcnYsZ1SMoHNMnVN4/upkLL2KeTts+itKLv+qg4sS8+8D5bFpzuw39L4MUfS\nsJwncU/BvG/Ez5NwHMfJZxQIhnrgO64dx3GcXNwLrOM4jpOLzySakPt23QGo7pm0mmfTrHL9dfsD\n/R1lmZKGew1cq8i7V5a7ilT/XUpfkdyrmh+lrHtNZnlvuCPRv5d8Eg3VC2s130h5Ov9q98/ynZSm\nd8T4WtcMSu3qWNVX/+RX+sxKleUFdTC+k6qdHFfNt9IrGekrM9JTE9Rq6w8p1Q61mlD2Ws698bXc\ni9C/VLlvk9JyA2gGY6GPjuM4dWF80RG0if1XuZBwHMepkSpHvffRxEKi5Uxg/7f9cUB8lopjY17o\nDW/B0t7wh7kIgB8zqzcuSx2VqpayzFKrmbjmUVKDVFMHpSqmtF+lch0ZcQCv4mkAduTR3rhpPNUb\n3izaMHYk+oy0/lK9qbpn0qq++ktur9vzdtyW4tPtQS8m4VJ8WmZZRt7nctKz1DGvZIQTFc2aJP2l\nGH4u+VOvYGA4jevJSF9eJX1NjenNysmjbJwZLhPY5QV9r05+pXlNYH3HdRn7cwf7c0ejm+E4ThPQ\n3l7samaavPmO4ziNY3zeAn4L4ULCcRynVsbACNpyXfzLHw8KgSwXGxv3RXVs2acIn7HpQKe+d9Dn\nxv2p5dMAWL1yoIdVgJ41GWar4xMTy/YME9L2NL1P6zy5I2il03WESUl4ctRap2sGm/Bkb3hvbgNg\nN/r8k7whCU/508shMCdpTGquWFqeSPX4iX7/xRj/fKLTfyrJ+lJ8zdLjQ5+uPi+9JyOu1fTzTgvR\nciPoQMZAFx3HcerEGBhBx0AXHcdx6sTQjnpvCuoqJCT9DDgUeNrM3hDjNgUuAV4DLAQ+aGYvxLST\ngE8Q9oV+wcyuj/F7AhcQ9n7OMbPjc29asuxMd4muH1/7rF5Z/eyGveE7XrtnX8IWA6t8+dEpsZ4+\nM771t+yzwdxg41XlRVh/3Eu94cmJ8mSDqJDJM8HNMjHdgL661o/hNG47FvaGD+BmADa5J1HYPJI0\nrKRa+nMS91hfcHGwkO1nlVrNBDTNuyYjvWcQ6Wm844x6xsBjdr1NYH9OOG0p5UTgBjPbCbgxvkfS\nzsARwM6xzI8lleyKzwGONbPpwHRJ5XUOOw+xW+/lOI6TyYSCVxNTVyFhZjcD/yyLPoy+0/xmA++J\n4cOBi81sjZktJMwJ9pE0FdggHhIOcGFSxnEcp3HUeDSdpImSbpN0t6R7JXUlacdJeiDGn5F3a0lt\nkuZLumY4u1ROIyZLW5hZSb+ylD4FzzT629ksArYiaCgWJfGLY3w2JdVSlrop5R99wXX3JtsmZw3M\n2tuqbfs2TK7euM/SaePJfaqj0u7rWs5FhuqO80qkjvwyd3rnPb1kVZ80dU3Za3m4UlytDGddjjOi\n1DiCmtlKSQeZ2XJJ7cAtkq4DJhMepHczszWSNq9QzfHA/ZAcWF8HGqpRMzOTNLz79S/vCq/twC6d\nsGvnsFZfidSzrOM4o4fu7m66u7uHv+IhLFybWckivIPgddyAzwLfMbM1Mc8zWWUlbQ28E/g28KXa\nW1GdRgiJpZK2NLMlUZUUl0pZDGyT5NuaMINYHMNp/OLc2t/fFV4n5uZwHGeM0dnZSWdnZ+/7U045\nZXgqHsIIKmkccBewA/BDM5snaSfgzZJOI3gY+4qZZfkJ+k/gq8CGGWnDSiN8N10NHBPDxwBXJfFH\nSuqQtB0wHZhnZkuAFyXtExeyj0rKOI7jNI4a1yQAzGydmc0gPPjuI2mXmHsTM9uXIAQuLS8n6V0E\ni9H5QN2dBtbbBPZi4EBgiqQngW8ApwOXSjqWaAILYGb3S7qUoGPrAWZZn4vaWQQT2EkEE9jf5d50\nbnxNj4zKmlWkuvmXq3SkO74mSxerb+0T4Iu23HBAPU+nivb0jMPS+siUJC41uy1pF19K4jYfaHq7\nweS+DAfwp95waU1jv9f9T2/clssSN6kbxdd0zSKZMo8ve4XsZYxq6SlZu6fzyvuOaqepyBlBu58P\nVxHMbJmkuQSrzkXAFTH+dknrJG1mZqnf4/2AwyS9kzC6bSjpQjM7uuZ+VKCuQsLMPpSTdHBO/tOA\n0zLi7wTeMIxNcxzHGTo5BiKdU8NV4pS/9U+XNAXoMbMXJE0CDiE8QL8MvAW4KaqeOsoEBGb2NeBr\nsZ4DCSqpuggIGBNbQRzHcepE7SPoVGC2pDaC2v8SM5sjaTzwM0kLgNXA0QCSpgHnmtmhGXXV9bCO\nljt0iB1jfyZlZEhNYVMVVOL4z67IqHe/GEjPEs7Si6RxKzPSy9tQIn0aKbU71dFkqc6SNvOmvuCr\nvvQE0Hd4EsARXNIb3ve66MzwvKR84t/w+VCcF5O+ZDnYy3PQV21HdcmcI9WmZZXP2qVdHu/Uj0k5\n4dJ3kX5/RWjVQ4fsAwXz/nfzHjrkMwnHcerCaBMMdWEMjKBjoIuO4zh1YgyMoGOgi47jOHViDHiB\nbb01iSkV+pPq9ssPMI+2UxXXJFLnIKnZbJaivkgbSmT5pchTvmdpNdN1jtfG10P6ojq+0uen9Vub\n/jsAX73+B30ZfpKUL3mEzfMq0lP2Cv0OJcpKX5EcYFQ6tKjamkZK1keWPt1MzoifNIjNlD1JW0vn\nEY9vGxgHyXGVaQPSe00oey1PX69Kelb5rHulv98073oZ6Wm49FtJt2BtlIQ3yyiT3N/i56LEqrrf\nAVZnJeG/jd6xZdjWJD5ZMO95vibR9GQJh960uOVArx6ZtjiO0ySMAc8OLiQcx3FqxdVNzYUkY1KF\n/qQisWwTTLYbrbL60yl6qm5Kb1m6R822mtWMQLNIFTLRYDGd2L4uCUcn65t9q8/91dfG9e1f3J35\nQH/PslmOC9sSfVRbjR5vsyh5sV2b/Puy2rI8UTItT4w0V8T4ND1tf6mu/vUP7F/qZTfNm9W+nn7p\n7bll8u6fVVd6//7t7xiQ3r/dbQP6lFXX2nVJm1Ymn288xz09t33d2ox79bRlht+2fZ/X6t9zWGYb\nRwPDpm76csG8Z7q6yYlYHM/VlD8Hx3EGxRgYQcdAFx3HcerEGFA3uZBwHMeplTEwgrbemkRpgaCa\nuqfMKqH3+I+K9VfPU/o48/OWDD7z1hl6qqQXPcctXadIFlMmRf39jCQ5PTG8ZEKbnnWV9Ueo9udI\n06u5eS0al8bnfQyryvLVWn+1uKHev1pd6TJPtbrS8lnLQ4MpX6meIryrL2ifqLGOEWDY1iS+WTDv\nN3xNwikjlb2+PuE4LUreMcEthAsJx3GcWhkDI+gY6KLjOE6dGAMj6BjoouM4Tp1w6ybHcRwnlzEw\ngo5rdAMcx3GalvaCVxmSJkq6TdLdku6V1JWkHSfpgRh/RkbZbSTNlXRfzPOFenStxBiQg47jOHWi\nRnWTma2UdJCZLZfUDtwi6TqCU+PDgN3MbI2kzTOKrwG+aGZ3S1ofuFPSDWb2QI29qEjrCgmjb69E\n1laQIm69K7I0CacHPG5YnrGMktvuWg/lLLpPIscZ94rY1r8kbf5L0ubSZ5b34y9tv8hrRlFX6CmD\nsckvuq2nUWbHrbPtqDjpKLIgCY/ifRLDxhC8wJr17s7qIPxzDPgs8B2z4ODHbKBXOTNbAiyJ4Zcl\nPQBMA+oiJFpa3WTrwjVs9VnfNRLlHMcZ5dSobgKQNE7S3YQnzuvNbB6wE/BmSbdK6pb0vyrdXtK2\nwO7AbcPQm0xadybhOI5Tb4Zg3WRm64AZkjYCrpS0C2FM3sTM9pW0F3ApsH1W+ahqugw43sxezsoz\nHLiQcBzHqZWcEbR7QbiKYGbLJM0lOMhZBFwR42+XtE7SZmb2XFpG0njgcuCXZnZVze0vgAsJx3Gc\nWskZQTt3D1eJU37dP13SFKDHzF6QNIlw4PDphJNq3gLcJGknoCNDQAg4H7jfzM4epp7kUtc1iTxT\nLUmbSrpB0sOSrpe0cVLmJEmPSHpQ0tuS+D0lLYhp36tnux3HcQrRVvAayFTgj5LuAeYR1iTmAD8D\ntpe0ALgYOBpA0jRJ18ay+wMfBQ6SND9eMwfeYnioqxdYSVsCW6amWoSz0T4OPGtm35V0AkEHd6Kk\nnYGLgL2ArYA/ANPNzCTNAz5vZvMkzQG+b2a/K7ufpeYl1T2y9jHYj0HKtm4yq2zd1Feu3tZNKRkn\n1+VZZLl109AYi8YJ6dP0tn1Be2SkG1KcYfMCe1PBvAc2rxfYus4kzGyJmd0dwy8TTLS2ItgBz47Z\nZtN7qCaHAxeb2RozWwg8CuwjaSqwQVz9B7gwKeM4jtMYap9JNA0jtiZRZqq1hZmVHqmXAlvE8DTg\n1qTYIoJQWRPDJRbHeMdxnMYxBlZ1R6SLUdV0OcFU6yUl+p+oShrGSXpXb6i7u5POzs7hq7puDEZ3\nU8suwJyNdb0s7guWvol+GrDkZ9IzINCfwWjOhkSVn65lfXbNymD+pllfQD3+5snn25OoK1+YNDDr\nKKC7u5vu7u7hr3gMCIm6n0wXTbV+C1xXWomX9CDQaWZLoipprpm9TtKJAGZ2esz3O+Bk4O8xz+tj\n/IeAA83ss2X3ylyTyG5X//eN3OgmZStvzaZn5P3rEO9Wy+CZ9U8YMWmQQ7V/pwuJ2soXJefkwynJ\n+tyAvcKjh2Fbk7i7YN4ZviaRSQVTrauBY2L4GOCqJP5ISR2StgOmA/PiNvQXJe0T6zwqKeM4jtMY\nfE1iyJRMtf4qaX6MO4lgD3yppGOBhcAHAczsfkmXAvcTHolmWd9UZxZwAcEsZ065ZZPjOM6I4+qm\n5sLVTUVwdVPz4eqm4WbY1E2PFcy7ffOqm8aAHHQcx6kTY2AEHQNddBzHqRNjYAQdA110HMepE2Ng\nBB0DXXQcx6kP1uSWS0VwIeE4jlMja8fACDoGuug4jlMfXEg4juM4uaya0FEw5+q6tqOeuJBwHMep\nkbVtrb8o4ULCcRynRtY2u8+NAtTVd5PjOE4r00NboascSRMl3Sbp7nhqZ1eSdpykB2L8GVn3lTQz\nnt75SDy4rW74TMJxHKdG1tY4hJrZSkkHmdlySe3ALZKuAyYTDmXbzczWSNq8vKykNuCHwMEEP/+3\nS7razB6ouSMVcCHhOI5TI0NRN5nZ8hjsIDjDMuCzwHfMbE3Mk+UBa2/g0Xh6J5J+TTjV04WEk/J8\nEk6dulU7ULqaY756/ySKOgas1YHfaPxJt6KzwbRPyW/xWWDK60ewPY1lKEJC0jjgLmAH4IdmNk/S\nTsCbJZ0GrAS+YmZ3lBXdCngyeb8I2KfmhlRhNP6jhkRRb67N6vzWbDcApO7GNsRxchjN3l+Hm1Vk\nm8DO617JvO6VFcua2TpghqSNgCsl7UIYkzcxs30l7QVcCmxfXnTIDR8EhYWEpPWAlWa2to7tcRzH\naRry1iTeX8lDAAAdzUlEQVT27FyfPTvX733/o1OW5dZhZsskzQVmEmYFV8T42yWtk7SZmT2XFFkM\nbJO83yaWqwu51k2S2iR9WNK1kp4GHgKWxFX3/ytpx3o1ynEcpxlYS1uhqxxJUyRtHMOTgEMIawpX\nAW+J8TsBHWUCAuAOYLqkbSV1AEcQTvWsC5VmEn8EbiScJHdfaQYhaTPgIOAMSVeZ2S/q1TjHcZzR\nzBDWJKYCs6Ol0jjgEjObI2k88DNJCwjbtI8GkDQNONfMDjWzHkmfB35POBz1/HpZNkGFk+kkdZhZ\nxb3kksaXVuFHA5KsGU/aG8zJdH1lupN3vnBdvFwjaPWF60n9cpiN/oXr4TqZbp7tWijv3rq39U6m\nSwWEpE2AV5Mc6W1md40mAeE4jjPS1LpPopmo2kNJpwIfAx4D1iVJB9WpTY7jOE3BWHDLUUQMHgHs\nUE315DiOM9ZYnWMC20oUERL3ApsAS+vcFsdxnKYiyy9Tq1FESJwG3CXpPmBVjDMzO6x+zXIcxxn9\n+JpE4ELgDMKMorQm0XwmRI7jOMPMaF+TiHsw3gUcAEwDVgALgGvN7L4idRQREq+Y2fdrbqXjOE6L\nMpqFhKRTgHcD3cCtwDPARGAn4PQoQL5kZn+tVE8RIXGzpO8QdvSV1E2Y2V1VGjgRuAmYEO9zmZl1\nSdoUuAR4DbAQ+KCZvRDLnAR8AlgLfMHMro/xewIXxA7OMbPjC7TbcRynrozyNYl5ZnZyTtqZkrag\nv3uPTIoIiT0I6qV9y+IrmsBW8Jf+fuAGM/tuPCzjROBESTsTLKl2Jng5/IOk6XF33DnAsdFL4hxJ\nM83sdwXa7jiOUzdWM6HRTcjFzK4FkPQBM/vvNC2Jq2qQVFVImFnnEBqZ5S/9MODAGD+bMBU6keAP\n/eK4QW+hpEeBfST9HdjAzObFMhcC7wFcSDiO01BGs7op4WvAfxeIyyRXSEg6CvhVdGeblb4jMNXM\nbq5QR5a/9C3MrCS9lgJbxPA0gt6sxCLCjGIN/T0cLo7xjuM4DWU0q5skvQN4J7CVpO8DJbcgG1Dd\nf08vlWYSmwHzJd1F8Dr4LGFNYEfgzcBzhBlALhn+0nctSzdJbik1rJQOgFmRxGX9HlIfSrWY8Q3V\nH1G132gr+TtqRdLfzKSytNHvu2m4GOUmsE8BdxK0NHcm8S8BXyxaSSXfTWdL+iHBbe3+wG6EkecB\n4Cgze6LoTRJ/6W8Hlkra0syWSJoKPB2zlftI35owg1gcw2n84rx7dXV19YY7Ozvp7Ows2syGUcmR\nX36Zzsx46VdDbI3jtB7d3d10d3cPe72jWd1kZvcA9ygMCuOA1xFU/g8OxoNGrhfYoSJpCtBjZi9E\nU6vfA6cDncBzZnaGpBOBjc2stHB9EeH81q2APwA7xtnGbcAXgHnAtcD3sxaum9UL7HDSJyR8JuHU\nk/yZhNkRI9uUGhguL7Dn2McK5f2cLmiYF1hJhwL/RfC/B+Gku8+Y2Zwi5es5V8rzl34rcKmkY4km\nsABmdr+kS4H7CSPYrGTEn0UwgZ1EMIH1RWvHcRrOaJ5JJJwFHGRmjwJI2gGYE6+q1E1ImNkCgvls\nefzzwME5ZU4juAEpj78TeMNwt9FxHGcorBrFJrAJL5YEROQx4MWihXOPL3Ucx3EqM4TjSydKuk3S\n3ZLuldQV47skLZI0P14zs+4r6Yux3AJJF0mqJK3ujPvLPibpY8BvgTskvU/S+6r1sch5EhMJG+C2\nTfKbmX2zWlnHcZxWplZ1U4XNxgacZWZn5ZWVtBVwHPB6M1sl6RLgSMK+sywmEgyESvvTSu453h3f\nX1GprUXUTb8BXiCYUK0skN9xHGdMMJR9EjmbjaFvP0Ml2oHJktYCk6lg8WlWcHW9wo2qsZWZvX0o\nN3Ecx2lFhrJPImez8TuA4yQdTdif9uWSb7sSZrZY0pnAEwQzxt+b2R8q3GcH4GzgfxME0f8AXzSz\nx/LKpBTp4f9I2q2ap0DHcZyxRp66aWH33/l7998rls3YbLwLwU9dSZV/KnAmcGxaTtImBPdG2wLL\ngP+W9BEzy9skdRHwQ6C0/nAEcDGwT8UGRiq55VgQg23AxyU9Tv9Dh3YrcgPHcZxWJU9IbNO5Pdt0\nbt/7/k+n3JJbR7LZeKaZnVmKl3QecE1GkYOBx83suZjvCmA/IE9ITDKzXyTvfynpq7kNKqPSTKK0\nqGEM1JGN7R1rjuM4wKoaz7jO2Gx8COGMhy3NbEnM9l7CAUHl/B3YN5ZbSRAa8zLylbguHsNwcXx/\nRIzbFHq3JeRSyS3HwtiZX5jZUWUd/AVwVFY5x3GcscIQ1iTyNhtfKGkG4UH8ceAzAJKmAeea2aFx\n7eIywnpGT3z9aYV7HRHr+3RO/PYDSiQU6WE/p3zRXGvPAuUcx3FamiGYwOZtNj46J/9TwKHJ+y6g\nq+C9tq2ljSVyN9NJ+pqkl4A3SHqpdBHsba8eyk0dx3FagVo3040Ekg6okr6hpKqeLCqpm04DTpN0\nuplVdAnuOI4zFhnN50kA75f0XeA6wj63Zwj+73YkOFrdFvhytUqKqJtOkvR+4E3AOuAWM7uytjY7\njuO0DqP5PAkz+z+SNiN4zPgAYR2kdNzDT8ws3+QqoUgPf0zY7HExwcrps5IOMbNZNbXcGSFW5MRn\nuegufEhVgfqHi3rXP1rpqZ5lVFD+m9mxIa1oNKPdC2w0k/0plRe2K1JESBwE7Fw6xlTSBQR33s4o\nxOwjAEg/aHBLnLFEM5whUQ9W12gC20wUERKPAq8mnP1ADD+am9txHGeMMMrXJIaFIkJiQ+ABSfMI\nNrV7A7dLuoaw8/qwejbQcRxntDKa1ySGiyI9/EaFNN957TjOmGU0r0lI2ht40sz+Ed8fQ1jEXgh0\nVdtpXaKqkDCzbknbEs6b/oOkyUC7mRU+2chxHKcVGc1CAvgJ8FYASW8GTgc+D+xOWMj+lyKVFDl0\n6NPAp4BNCVZOWxM8Fb61llY7juO0CqN8TWJcMls4gmD2ejlwuaR7CldSIM+/EvZIvAhgZg8Drxpk\nYx3HcVqOtbQXuhpEm6TxMXwwMDdJK9yoIhlXxSPygF7fTb4W4TjOmGeUm8BeDNwk6VlgOXAzgKTp\nhNNGC1FESNwk6euEo/IOAWaR7ePccRxnTDGa1U1m9m1JfwS2BK4v7XUjbIo+rmg9RYTEiYSTkRYQ\n3NbOAc4bXHMdx3Faj9FuAmtmf8mIe3gwdRSxblor6SrgKjN7ejCVO47jtDKj3LppWKjkKlySuqI+\n6yHgIUnPSjpZpQUKx3GcMcxodhU+XFSybvoisD+wl5ltYmabEHZb7x/THMdxxjS1CglJEyXdJulu\nSfdK6orxXZIWSZofr5lZ95W0saTLJD0g6X5J+9arj5WExNHAh83s8VKEmT0GfCSmFUJSW+zsNfH9\nppJukPSwpOslbZzkPUnSI5IelPS2JH5PSQti2vcG00HHcZx6sYoJha5yzGwlcJCZzQBmADMl7UOw\nHD3LzHaP1+9ybv09YI6ZvR7YjeD+uy5UEhLtZvZMeWSMG8xqzfEEr7Els9kTgRvMbCfgxvgeSTsT\nNnzsDMwEfpyotc4BjjWz6cD0POnqOI4zkgxF3WRmy2OwAxhP3xhZUZ0vaSPgADP7Waynx8yWDU+P\nBlJJSFQ6ZKDQAQSStgbeSbCGKnX8MGB2DM8G3hPDhwMXm9kaM1tI8DS7j6SpwAZmNi/muzAp4ziO\n0zCGIiQkjZN0N7CUYKJaGuOOk3SPpPNTTUvCdsAzkn4u6S5J50Z3SXWh0oxgt3imdRaTCtb/n8BX\nCZ5kS2xhZktjeCmwRQxPA25N8i0CtiIIpEVJ/OIY7wwb1Q66qeVQIqeP4TxIKP3LNssBRa1L3j6J\nld3zWNk9LzOtRNy3MCPODK6UtAtBa/LNmOVU4EzCFoSUdmAP4PNmdrukswkamUrOWGum0hnXQ1qS\nl/Qu4Gkzmy+pM+ceJmlYd293dXX1hjs7O+nszLx1y2NWeK8M4Rhcp7lwATEYuru76e7uHvZ68/ZJ\njO/cj/Gd+/W+X3bKj3LrMLNlkuYCM83szFK8pPPI3ri8CFhkZrfH95cR1fb1oJ47QfYDDpP0TmAi\nsKGkXwBLJW1pZkuiKqm092IxsE1SfmvCh7E4htP4xXk3TYWE4zgODHxgPOWUU4al3lrNWyVNAXrM\n7AVJk4BDgNNLY2PM9l7CJuZ+xLHzSUk7xY1xBwP31daD6hRx8FcTZvY1M9vGzLYDjgT+aGZHAVcD\nx8RsxwBXxfDVwJGSOiRtB0wH5sUP7EVJ+8SF7KOSMo7jOA1jCGsSU4E/Rm+s8whrEnOA70r6a4w/\nkLjdQNI0Sdcm5Y8DfhXz7QacVq8+juSe8pJa6XTgUknHEg6/+CCAmd0v6VKCJVQPMMvMSmVmARcQ\n1kLmVDALcxzHGTFWra7NwZ+ZLSCsK5THZ24vMLOngEOT9/cAe9V080GivnG4+ZFkrdSfkaL6moQv\nXA+Neixcj641CbOTG92EQSEJMxuS5whJtv4rA3YJZPLyepsP+X6NYnR7p3IcxxnFrO1pbpcbRXAh\n4TiOUyMuJBzHcZxceta4kHAcx3FyWLe29YfQ1u+h4zhOvXB1k+M4jpPLytYfQlu/h47jOPVidFki\n1wUXEo7jOLXiQsJxHMfJxYWE4ziOk8sYcEbgQsJxHKdW1ja6AfXHhYTjOE6tuLrJcQbLGPjXOE6J\nlY1uQP1xIeFg9m+F80rfTsp9PSN9eA5zcSp7Vh0Nn3OzeX6tC2Pgmahuhw45juO0PD0FrzIkTZR0\nm6S7Jd0rqSvGd0laJGl+vGbm3VpSW8yTdcTpsOEzCcdxnFqpcSZhZislHWRmyyW1A7dIuo5wONtZ\nZnZWgWqOJxzStkFtrSiGzyQcx3FqZU3BKwMzWx6DHcB4+k7vrHo4kaStgXcC5xXJPxRcSDiO49TK\n2oJXBpLGSbobWEo443peTDpO0j2Szpe0cc6d/xP4KrBumHqSiwsJx3GcWqlxTQLAzNaZ2Qxga2Af\nSbsA5wDbATOAfwBnlpeT9C7gaTObT51nEeBrEo7jOLWTZwL7UDc83F2oCjNbJmkuMNPMeoWCpPOA\nrEXp/YDDJL0TmAhsKOlCMzt6ME0vigsJx3GcWslbuN6hM1wlftvfZFnSFKDHzF6QNAk4BDhd0pZm\ntiRmey+woLxqM/sa8LVYz4HAV+olIMCFhOM4Tu3Uvk9iKjBbUhtB7X+Jmc2RdKGkGYRF7MeBzwBI\nmgaca2aHZtRlGXHDhgsJx3GcWqndBHYBsEdGfOaMwMyeAgYICDO7CbiptlYUw4WE4zhOrbgXWMdx\nHCeXMeAFtu4msJIWSvpr3D4+L8ZtKukGSQ9Luj61BZZ0kqRHJD0o6W1J/J6SFsS079W73Y7jOFVZ\nWfBqYkZin4QBnWa2u5ntHeNOBG4ws52AG+N7JO0MHAHsDMwEfiypZAd8DnCsmU0HplfyaeI4jjMi\nDGGfRLMwUpvpyjd8HAbMjuHZwHti+HDgYjNbY2YLgUcJm0ymAhskOxIvTMo4juM0hiG45WgWRmom\n8QdJd0j6VIzbwsyWxvBSYIsYngYsSsouArbKiF8c4x3HcRrHENxyNAsjsXC9v5n9Q9LmwA2SHkwT\nzcwkDZudb1dXV2+4s7OTzs7O4aracZwmpbu7m+7u7uGvuMlVSUWQWV33YfS/mXQy8DLwKcI6xZKo\nSpprZq+TdCKAmZ0e8/8OOBn4e8zz+hj/IeBAM/tsWf02kv1x8ql0KI7ZySNyaE7WoTij4bCeagzl\nMJ9Gfa7NhiTMbEh+jyQZnys43pwz9Ps1irqqmyRNlrRBDK8HvI2wzfxq4JiY7Rjgqhi+GjhSUoek\n7YDpwLy4Tf1FSfvEheyjkjKO4ziNYQysSdRb3bQFcGU0UGoHfmVm10u6A7hU0rHAQuCDAGZ2v6RL\nCQdp9ACzkqnBLOACYBIwx8x+V+e2O47jVGZVoxtQf+oqJMzscYLL2/L454GDc8qcBpyWEX8n8Ibh\nbqPjOE7NjIE1Cd9x7TiOUytNrkoqggsJx3GcWmly89YiuJBwHMepFVc3OY7jOLm4kHAcx3FyGQNr\nEiPlu8lxHKf1WFXwKkPSREm3Sbpb0r2SumJ8l6RF0Wv2/CxHppK2kTRX0n2x7Bfq1T3wmYTjOE7t\n1H4y3UpJB5nZckntwC2SriP4ujvLzM6qUHwN8EUzu1vS+sCdkm4wswdqa01lXEg4juPUyhDUTWa2\nPAY7gPH0nVVd0X1H9ECxJIZflvQAwQlqXYSEq5scx3FqZQheYCWNk3Q3wRP29clRCMdJukfS+emB\nbDl1bAvsDtw21K7k4ULCcRynVoZw6JCZrTOzGcDWhHNzdiEcrrYdwVPFP4Az824dVU2XAceb2cvD\n06GBuLrJcRynVvLWJF7phuXdhaows2WS5gIzzaxXKEg6D7gmq4yk8cDlwC/NrK7OTl1IOI7j1Ere\nmkRHZ7hKPNvfhbukKUCPmb0gaRJwCHC6pC3jmgPAewles8vLCjgfuN/Mzh5iD6riQsJxHKdWavcC\nOxWYLamNoPa/xMzmSLpQ0gzCIvbjwGcAJE0DzjWzQ4H9gY8Cf5U0P9Z3Ur08Y4/ooUP1xg8dcooy\nlMN5hutQoFY4vKdZGbZDh6YWHG/+0byHDvlMwnEcp1bGwI5rFxKO4zi14l5gHcdxnFzcwZ/jOI6T\niwsJx3EcJxdfk3Acx3FyGQMzCXfL4TiO4+TiQsJxHMfJxYWE4ziOk4uvSTiO49RM669cu5BwHMep\nmdZfua67uknSxpIuk/SApPsl7SNpU0k3SHpY0vXpwRqSTpL0iKQHJb0tid9T0oKY9r16t9txHKc6\nawpezctIrEl8D5hjZq8HdgMeBE4EbjCznYAb43sk7QwcAewMzAR+HN3iQjiM41gzmw5Mzzog3HEc\nZ2RZUfBqXuoqJCRtBBxgZj8DMLMeM1sGHAbMjtlmA++J4cOBi81sjZktBB4lnNg0FdggOd7vwqSM\n4zhOg/CZxFDZDnhG0s8l3SXpXEnrAVuY2dKYZymwRQxPAxYl5RcBW2XEL47xjuM4DWQI55c2CfUW\nEu3AHsCPzWwP4BWiaqlEPADCD4FwHKcJaf2ZRL2tmxYBi8zs9vj+MuAkYEnpmL6oSno6pi8GtknK\nbx3rWBzDafzirBt2dXX1hjs7O+ns7Bx6L5yWo1EH/vhBQ42hu7ub7u7uOtRc2yxB0kTgJmACYRy+\nzMy6JHUBnwSeiVkzT5yLa7JnA23AeWZ2Rk0NKdLWep/kJulPwCfN7OH4AUyOSc+Z2RmSTgQ2NrMT\n48L1RcDeBHXSH4Adzcwk3QZ8AZgHXAt8v/zD85PpHMcpwrCdTMc9BXO/ccD9JE02s+WS2oFbgOMJ\nBjsvmdlZFe7bBjwEHEx4WL4d+JCZPVBLP6oxEvskjgN+JakD+BvwcYL0u1TSscBC4IMAZna/pEuB\n+wkielYy6s8CLgAmEayl6nKeq+M4TnFqt1wys+Ux2AGMp0/tXk147Q08Go17kPRrgtFPcwoJM7sH\n2Csj6eCc/KcBp2XE3wm8YXhb5ziOMxRqX5SWNA64C9gB+KGZzZP0DuA4SUcDdwBfNrMXyopuBTyZ\nvF8E7FNzQ6rgO64dx3FqJm9R+u545WNm64AZcavAlZJ2IewH+2bMcipwJnBsedGam1sDLiQcx3Fq\nJm8msWu8SszOyQdmtkzSXGCmmZ1Zipd0HnBNRpFyA59t6L9FYFhxL7CO4zg1U5sJrKQpJXdEkiYB\nhwAPSNoyyfZeYEHGTe8geJ3YNq71HgFcPVw9KsdnEo7jODVT85rEVGB2tFQaB1xiZnMkXShpBkGl\n9DjwGQBJ04BzzexQM+uR9Hng9wQjoPPrZdkEI2ACO5K4CazjOEUYPhPYywvmfv+Q79cofCbhOI5T\nM83tvK8ILiQcx3FqprldbhTBhYTjOE7NNLfzviK4kHAcx6mZ1p9JuAnsKKc+TskaTyv2qxX7BK3b\nr+HBXYU7DaZV/6Ct2K9W7BO0br+GB3cV7jiO4+TS3LOEIriQcBzHqZnWN4Ftuc10jW6D4zjNwfBs\nphu5+zWKlhISjuM4zvDiC9eO4zhOLi4kHMdxnFyaXkhIOlXSPZLmS/q9pKlJ2kmSHpH0oKS3NbKd\ng0XS/5X0QOzbFfFgEqJ74BWxv/Ml/bjRbS1KXp9iWjN/Vx+QdJ+ktZL2SOKb+bvK7FNMa9rvKkVS\nl6RFyfczs9FtGpWYWVNfwAZJ+DjgnBjemXA01HhgW+BRYFyj2zuIfh1Sai9wOnB6DG8LLGh0+4a5\nT83+Xb0O2AmYC+yRxDfzd5XXp6b+rsr6eDLwpUa3Y7RfTT+TMLOXkrfrA+ti+HDgYjNbY+HA8EcJ\nB4g3BWZ2g4XjDQFuA7ZuZHuGgwp9avbv6kEze7jR7RhOKvSpqb+rDJrS4mgkaXohASDp25KeAD4M\nfCNGT6P/kX6LCAeINyOfAOYk77eL0+NuSW9qVKOGSNqnVvquymmF7yql1b6r46L68/zSSXFOf5pi\nM52kG4AtM5K+ZmbXmNnXga9LOpGgcurKqWpU2ftW61fM83VgtZldFNOeArYxs39GXfFVknYpm1E1\njBr7lEXTfVcZNP13VZBR9V2lVOjj14FzgG/G96cCZwLHjlDTmoamEBJmdkjBrBcB1xKERPlh4VvH\nuFFDtX5J+hjwTuCtSZnVwOoYvkvS34DpwF31a2lxaukTLfBd5ZRp6u8qh1H/XaUU7aOk84DBCMYx\nQ9OrmyRNT94eDpTOer0aOFJSh6TtCH/OeSPdvlqJlhZfBQ43s5VJ/JR4Li6Stif067HGtHJw5PWJ\nJv+uyujVcTfzd1VGqrdvme8qtYQE3gssaFRbRjNNMZOownckvZawYL0Q+CyAmd0v6VLgfoIXrlkW\nTRqahB8AHcANkgD+YmazgAOBUyStIfT5M2b2QuOaOSgy+9Ts35Wk9wLfB6YA10qab2bvoIm/q7w+\nNft3VcYZkmYQ1GWPA59pcHtGJe6Ww3Ecx8ml6dVNjuM4Tv1wIeE4juPk4kLCcRzHycWFhOM4jpOL\nCwnHcRwnFxcSjuM4Ti4uJBzHcZxcXEg4juM4ubiQcOqCpJfrXP+f4+tGkj5XQ/kuSV8eZJlJ0Zvr\nAPfS8XCe+cl1QozvkHRTyT2H4zQbreCWwxmd1HUrv5ntH4ObALMIHj0HVUUNt/0EcDmwpaQjgVVA\nt5ndDyw3s90z2rla0o3AEQQHlI7TVPhMwhkxJH1J0oJ4HR/jto1Hmv5U0r0KR9BOTMr8Rzwm82ZJ\nF5ee/pOZyunADvHp/QxJr5G0ICn/FUknx/DXJT0k6WbgtUmej0q6LdbxX5Ly/hcfBn5DOGRnNvAT\nYK8CXb8K+EjBj8lxRhUuJJwRQdKewMcIA+y+wKeiczWAHYEfmtmuwAvA+2OZvYD3AbsB7wD2pG8G\nUHo9Afibme1uZicw8KQxS+5/BPBGgqvyvQCT9Hrgg8B+cSawjowBXVIHsL2ZPUFwR/8u4KOEmQXA\npDJ10weS4vdRTJg4zqjD1U3OSPEm4AozWwEg6QrgAILr6cfN7K8x352Es5MB9geuKp3LICnL33/R\n4ydL918JrJR0dSz7FoLwuSMuNUwClmSUn0IQYJhZD3BhWfqKLHVTzL9W0mpJ65nZKwXb6zijAhcS\nzkhh9B/QRd9sYFUSv5YwUOeVqUYP/WfIk5JwVl0CZpvZ16rUuwKYWCVPJSYAK6vmcpxRhqubnJHi\nFuA90UJoPeA9wM1UHvj/DLxb0gRJ6wOHZuR5Cdggeb8UeJWkTSVNIKiFDPhTvP9ESRsk8TcC/yJp\nc4BY7tXlNzGzfwJtUe00KCRtBjxrZmsHW9ZxGo3PJJx6MVnSk8n7s4AL6DvF7Fwzu0fStgy0NDIA\nM7sjqoX+Shj8FwDLyvI8J+nPcbF6jpmdIOmb8T6LCYfjYGbzJV0C3AM8XWqHmT0g6d+B6+OC9RqC\ntdQTGX26nqAiuzEjbZKk+cn765LZyUHAbzPKOM6oxw8dckY1JT2+pMnATcCnzOzuBrVld+CLZnb0\nIMtdDpxgZo/Wp2WOUz98JuGMdn4qaWfCesAFjRIQ0DsbmStpnJmtK1JG0njC4rsLCKcp8ZmE4ziO\nk4svXDuO4zi5uJBwHMdxcnEh4TiO4+TiQsJxHMfJxYWE4ziOk4sLCcdxHCcXFxKO4zhOLi4kHMdx\nnFz+PxF/SAG5J95QAAAAAElFTkSuQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f35d0ee7e10>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"m1 = netCDF4.Dataset('/archive/aja/ulm_201505_awg_v20151106_mom6sis2_2016.01.16/CM4_c96L32_am4g7_2000_lowmix_hycom5/gfdl.ncrc3-intel15-prod-openmp/pp/ocean_annual_z/av/annual_20yr/ocean_annual_z.0021-0040.ann.nc')\n",
"for v in m1.variables:\n",
" print(v, end=',')\n",
"x = m1.variables['xh']\n",
"y = m1.variables['yh']\n",
"z = m1.variables['zt']\n",
"j = find_index(y, 38.0)\n",
"i0 = find_index(x, -30.0);\n",
"i1 = find_index(x, -5.0)\n",
"q = m1.variables['so'][0,:,j,i0:i1]\n",
"plt.pcolormesh(x[i0:i1],z[:],q);\n",
"plt.gca().invert_yaxis();\n",
"cb = plt.colorbar(); cb.set_label('S (ppt)')\n",
"plt.xlabel('Longitude ($\\degree$E)');\n",
"plt.ylabel('Depth (m)')\n",
"plt.title('38$\\degree$N, yrs 21-40, CM4_lowmix_hycom5')\n",
"plt.clim(34.9, 36.7); plt.xlim(-30., -5.); plt.ylim(6500, 0);"
]
},
{
"cell_type": "code",
"execution_count": 82,
"metadata": {
"collapsed": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"xh,yh,zt,zw,time,nv,xq,yq,thetao,so,age,uo,vo,uh,vh,thetao_xyave,so_xyave,time_bounds,average_T1,average_T2,average_DT,"
]
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAEfCAYAAAC04jrjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztnXmcXFWZ97+/pLOybxMSQNmCCoosgzA4SqPgRFFwGQUX\nQGXcMiKv2wA6M3Z0RPB9YXCDUUAJLgiDgChBQaRRnJGwBAwQRJQoCRA2CUvI/rx/nHPTpyv31r1V\n3dXVVf18P5/69K2z3XPqVp+nznOe8zwyMxzHcRwnj3Ht7oDjOI4zenEh4TiO4xTiQsJxHMcpxIWE\n4ziOU4gLCcdxHKcQFxKO4zhOIS4kHMdxnEJcSDiO4ziFuJBwnDYgabGk17bx/ndJenW77u90Di4k\nRgBJ35P0sKTlkn4v6YQkb2dJ8yQ9Gct8TdL4mNcjaY6kf5E0O6mzWNIySVOTtH+SdMPIjmwASRMl\nXRD79rSkBZJm1ZT5qKRbJa2U9J0G2p4Z63y3Jn1rSVdIejbe950N9vldsT/PSHooPodXxrzFklZJ\n2qamzgJJ6yW9oGo/C7D4agtm9lIz+9VwtyvpeEmHDHe7TvtwITEyfAnYxcy2AI4E/kPSvjHvHGAZ\nsD2wD3AIkAmEWcBPzOzLwFRJmydtjgNOamWnJfU0ULwH+AvwajPbHPhX4FJJL0zKLAW+AHy7wa58\nA5jPxpPqN4CVwN8A7wbOlbRnlQYlfQL4T+A/Yv2dYntvSor9CXhnUudlwJScfpT1s+uR9EFJbx54\nO+i908G4kBgBzOxuM1uZJgG7xeudgUvMbLWZLQN+BmQT3ePADvF6C+D5pP7/Az4laYtG+iLp05Iu\nq0n7qqSz4/XiuHL5HfCMpPGSTpa0JK4Q7pX0mpwxrjCzOWb2l/j+auABYL+kzBVm9mPgiQb6ewzw\nV+B6QEn6JsBbgX+L9/4N8GPg2AptbgHMAWab2ZVm9ryZrTOzq83slKy7wPeA45KqxwMXpf0o62fF\nMb5EUr+kv0Y10JuSvPdJuip5/wdJlybvH5S0d7xeLOlTku6Mq6vzJU2TdE18dtdJ2jIp+9p4vZuk\nJ7IfLpJmSHosTx0l6ei48speq+IK9gJgd+D/AKcBawjPw+lwXEiMEJLOkfQcsAh4CJgXs84GjpE0\nRdIOwOsJggIz+y2wVfzV+wszW5M0eSvQD3yqwa58F5iVCZe4WjgamJuUOSb2Y0vCP/4/A38bVwiv\nAxZXGO80YA/g7rzsKh2NK6c5wMdz6uwBrDWz+5O0O4G9KjT9d8Bk4IqScr8FNpf04qgCPJogOBrp\nZ10kTQB+Qnjm2wEnAt+XtEcs0g+8KpadAUwADorvdwU2MbPfxbJGEJyHET6fNxG+Z6fEtscBH0vK\nGoCZ/RE4GfiepCnAd4Dv5KmjzOwSM9vMzDYDZgB/BH4Qx70+Flsfirr30G7AhcQIYWazgU0J//BX\nAKtj1q+BlwJPAw8Ct8Rf21m9C83sLDO7sbZJ4N+BEyVt20A/HgF+Bbw9Js0CHjezBUm7XzWzpWa2\nClgHTAL2kjTBzP5iZn+qd4848X0fuNDM7svrRsXufgE438weyqmzKeEzS3kG2KxCu9sQxry+tGQQ\nqscBhwP3EFRmjfSzjIMIE/3pZrbWzG4AfkpUc5nZA4QV3b7Aq4GfAw9JehFBNVk7kX/NzB6Lffk1\n8FszuzM+yyuAfcnBzM4H7ieoy6YBn63XaUnjgIuBG8zsPOD9hJXj2QRV4yRXN3UHLiRGEAv8BtgR\n+IgkEX5B/giYCmwLbC3pjIrt3U2YUE6hsclpLvCeeP0eggol5cHkHvcTVAh9wDJJF0uaXtRwnDy+\nS9gr+GhRsZx6705UGFdL2gd4LWHSyavzLLB5TdoWBEFRxhPAtrGv9TDCWN5NgaqpQj/LmEHyeUf+\nzICaEeBGoJfwA+PG+DqEIDRqfzwsS66fr3m/kiBcizifsBL7Ws2qNY8vApsQVyZm9i0zy1ZmZmbf\nNLMrS9pwOgAXEu1hArArsDVhw/TrZrbGzJ4ELgTe0EBbnwM+wOBJpYwfA3tLeilwBOFXf8oggWNm\nF5vZq4AXxrxcIRaF3gUE1cbbzGxdwf03Emhm9v1MjWFmRxAmwZ2Bv0h6GPgk8DZJt8Yq9wE9knZP\nmnk5cFfBPVP+F1gFvKWsYNxj+RNB/XZ5TpHekn6WsRTYKX52GS8EliTvbwQOJQiJfgaExiFsLCRq\nqRVauT8mJG1KEHTnA3MkbVXYYNh/ORr4x9pnbGZzc1a9TgfjQqLFSNpO0jGSNo2bwP9A0Plfb2ZP\nEJboH4l5WxJ+sd5Ztf2oT76EGksnSReqwMw0bqJfRtAl32xmS/LKxXb2kPQaSZMIE+tKggoqj3OB\nFwNHRvVGbVvjJU0mWEKNlzQp6vrz+BZBkL6cYPX1X8DVwD/EMTxHmLQ/L2mqpL8n6OC/W2H8ywmq\num9IOirWnyDp9QWruBOA15jZ8zl536zXzwrcDKwA/iX2oRd4I/DDpEwmJCZHNdJNBDXh1sACGqNo\npfMVYL6ZfTD2/79yKwe119eAt8Tvr9PluJBoPQZ8mKBSeBL4MnCSmf005r+V8Cv1MeAPhIn44w3e\n4/MEdVX6K3FHwmRSxFzCXkiZTf8kggnvY8DDBJXYqbWFFExdP0iYLB9JVEfp2YV/I0yIJxPUXM9T\noPuOFkePxtcygnrp+ZqJaTbBJPVRwobyh81sUcyrO34zOwv4BEF//ijBfHc2OZvZZvYnM7s9TWqw\nn4VEtc6bGPgOfB04Nt3LMbM/ENRov47vnyZsGP+mwuaw1VxvVF7SUQSDhI/EpE8A+yn/3MmRBIOG\nm1L1YOlAnY5FnWSAoHA462xgPGGjsJLufqwhaSLhF+beRSofhcNgi4BpZvbsSPav1VQZv+M41egY\nIRHVEr8nmPctBW4B3pn8cnQqEjdszwI2NbN/and/HMcZvTRyorbdvAK438wWA0j6IXAU4dewUxGF\nQ2jLCHshs0qKO00SV2p5Z0QM2LPePpDjjCY6SUjswGBTwSXAgW3qS8cSN3zrmUE6w0C0iqpyZsNx\nRjWdtHHdGXoxx3GcLqKTVhJLCWcKMnZisC05klyQOI5TCTNr9ODjIBqdb4Z6v3bRSULiVmCmpJ0J\nvo+OJvHQmfEje/1GFdfFYa6jyCR/Y/LKrk3S1iUfXVZ23aD88bn18urX4+d98zms7+/q9i9tvyce\nYRhfcJShKL0eeZ9F0WeZN9Y8ftn3P7ym7+BKZXuSPo/PGd8kVtXNn5jk530+Zfl598rLv6DvIT7S\nt00st3pD/lRWJPdavVHalk8PGJf1PBovHmKA1Jj2yfh3eZKWnkjJ3Eg+V5Kfpj1XP7/vfujbmWIm\nJ9ebxL/jS/IBzm3fb7rBZxeb5z8qlvvXYblbe+gYIWFmayV9lOC7ZjxwQZ5lUz2BUDZxl5EnGNLr\novbzBEKRQMkrl9duUV+yqalYSKwtvFdRX2v7UyWtLH81E3m2jsq+TMil+auZWDd/PFM2XJcJgYH6\nA59TOuHn5Wc8x3IeZgYwIAxgsECYEq+nMnAub8XmG8KCsM2kxwHYZHziVmpScpNskt06SWtESKyq\nKVdUNs1/HNglXmcfVTr89CszKSeti5nQ7g6MAB31KM3sGuCaRuutipNI2QRYNIFW/dU/nnUNrVbq\nkQmD9YzLFThFQqaZ+zciMMsEXlXWMIEVTC3ML1vxpM/q+UQI5AuJdRvVy1udwMCqIhUMq5NZup4Q\nWcd4Vg2a0Tfu68SkbMag6574+fYM3H/QL/G8r2jeL/WiOnlCYFVO2TTtPuAFNfnp40kff96/Sm2f\nv9g9WuGOmkCbpOvGmKkD0n/WnpIJJ4+yOuk/fjpxD0wi6SSwcVtlqqusnT16tx80YWVtpeMrU43l\nUST4mpnwm6mzQ+9udfPTNssExtoGyuaR1sl77mVCJmPP3srOeAcxSEiMj9fpJJ83Caf56YSfkap1\n0gk/q180sWdf66T93r0IrhPTtopWErXt1F53GVPKi3Q8XSckuo3dendsdxdawo4lQqIT2bN3O4rd\nWnUuvbuXlxmruLqpA9kseoqekuiBn49qjbJf12Wrh0Z+nad68rw2ivYU8tLy8lOd94rk98yAaiT/\n/mVkv5SrqthCnbxf39V/Ppav2ur/us9Lb2RPIb1/9rnmrR7SshMHtVk/v5GVbPas1yZaq55Ncgqm\njyf9epTtGazNScvLX5eTlpbNW70U3b+L6boJNIexMEbHcZyW4CuJDmSbaC+4KvklPSlnszCPMn12\nWf2ifYK8enkbnGnZMhPTol+yq2P6qoKVRNUVQiP7DM3tAwxtpVG0MZ3XfpkJ7XCayGb3TVcveZZO\naZsT81Y66cezruC6HkV7Amtz0vKso4pWElnZ1CVkWjazlEpXGl28qmh2Ao0u828k2IP1AJeZWV/M\nO5HgkXgdcLWZnZxTf0sGgkQZ8P4Y7njY6TohkZob1tKsdVNZ/XXJxF3VHDadWMpMXAf3MeQXqTAG\n1B0bhXOo2+5wUXXyb0QFUy4QqlsvZdd5ZyvS63J1VfVzGnkmsIPOUawYMIedlJ1/yDsbAQOTcNkk\nn5eWpqf5z5bkp/fKa7/MBLeLhUSzKwkzWynpUDNbEePM3yTpGoLL/yMJHozXSNquoImvAPPM7B9j\n/Tyl5LDQdUKinQyX+avjOJ3BUNRNZpb9epgYm8piz3wpCx9rZo/V1pO0BfAqMzs+llnL4OOVw0on\n+W5yHMcZVUyp+MpD0jhJdxC8Ml9rZvOBPYBXS/qtpH5Jf5tTdRfgMUnfkXS7pPMkFR86GiK+knAc\nx2mSogn0jviqh5mtB/aJK4MrJO0Vm9zKzA6SdABwKSE8bu1t9wM+ama3SDobOIUQknfY6Toh0Yy7\njap69LwDdMNNnq4+vVeeTn6wb6D8vYih3L+MxkxUN/6syzam89LL9hTyNpur5GebyKX1VyX56wZc\naGQbzj1FewLP1fytvS7bk8g7MZ23P1F0/7wT1WV7Gnkb32mf00eVt/HdxVrYInXTAfGVMbdOG2a2\nXNINhPguSwix24kCYL2kbWrC4S4BlpjZLfH9ZQQh0RJc3eQ4jtMkPRVftUjaNlooIWkKcDghgNqV\nwGti+h7AxNp46Wb2CPBgzIcQrTMvwNWw0HUrCcdxnJFiCBvX04G5MSzzOOASM5snaQLwbUkLCf46\njwOQNAM4z8yOiPVPBL4f47n/EXhf812pT8fEuK6CJPtf2weo9YwazgzkudQuIk/dUXZ2ocxLbF6d\nIhox1y1Tp6SU2fnXS6u9b9798+o3YqJaVj/PVXeq+pm0Kqh+esrUJWUmokUmpmXnDPJOLOepm4pc\ndWfmqE/npBX1r0zdk2fCWtS/dTlp5OTntVnbbkrqtfa+9s85koYlnsSdFcu+HI8n4TiOU8woEAyt\nwE9cO47jOIW4F1jHcRynEF9JdCALeRnQmGfSPJ16M3sGVWIgVHUbUWaWm7p1yDzfptdp2mA/RKtj\n/QIT0nWZ76CCzyqae44v6F6WriKdeJ5OP69smc69aE8gz0Q0r/2i/DIT0DydfN51me+joqA/VcOP\npjTy+VXdUymzCk/aXJOUzb42PTX/Ht06mXbdBJrDWBij4zhOS5hQdQbtYP9VLiQcx3GapHbFVIgL\nidHDHQQT2Ik56pgis9BUdTODhwD4IwPhuFIvrxlpUKGBQDUD34RU3ZN3ojdd46ftTqp4YjpVEaUB\nljJX6dusenxD2tTnBk4EK2s+VWeUmTCWhaLMSyszsSxrs0xdVaQOyjMRzVO9lMV4bkRdM9Tx591/\nqM+n7PMrqh/rrSmZ1FJt5JoKh/S7Vd00oYtPk2f4iesa3s1lvJvL2t0Nx3E6gJ6eaq9OpsO77ziO\n0z4m5McO6ypcSDiO4zTLGJhBu26IDzEDGOwZdUU88jKV53Pr5LGMv9lw/SybAQP6/lqy/YH0num9\n8sxly7y1FpnglrkTydrN3FNAsg8BA15Gi3TeGUU67TwaMRGtuidRtifQ7J5EmduMMrcUZTr/vLRm\n3GaUjb/sXs08M8r3InKbyqkzFnT1QBfOoBszBoboOI7TIsbADDoGhug4jtMixsCKqaVCQtK3gSOA\nR83sZTFta+AS4IXAYuAdZvZUzDsVeD9hAfwxM7s2pu8PXAhMJgT/Pqnono+zDQCToooIYAUhsl+q\n4mkkaNDveREAW/LUhrTduH/DdaaOSvNTs9aU1HQ2j8zcNj29nee9Nm0/z/R20uYD+VuuHVB9KRt2\nOvz0Ok/dUJbfDtJ/zrU56UPtZ9mYiyaHrGxaP1UXlX3tsvxG1H3NUDC+UjVRrDehQFuamcNOmdxs\nxzqMMfAzu9UmsN8hRFtKOQW4zsz2AK6P75G0J3A0sGesc46kzLXuucAJZjYTmCmpts1h55uctOHl\nOI6Ty6SKrw6mpULCzH4N/LUm+UgGovnNBd4cr48CLjazNWa2GLgfOFDSdGCzGCQc4KKkjuM4Tvto\nMjSdpMmSbpZ0h6S7JPUleSdKWhTTzyi6taTxkhZI+slwDqmWdiyWppnZsni9DJgWr2cAv03KLQF2\nANbE64ylMT2XhatftlHapIlB9ZKegk5ZVSLqFz4Z2pw4eUCF89DUGQPtJ2qsZ6LqaXo8uR3yB+pl\nFlBpX1KrqSw/VTGtzulfqjpLLammET7anXhwQ9pOWw9cT9869Gvak8s3pCm1dMqz3iEnv4iqJ36L\njLuasd7Js3Qqc/DXrIPBsv7l5TfrLLDq6fFGTn83QvYVLJol4vdmSmL0NyUda4f/gq5EkzOoma2U\ndKiZrZDUA9wk6RpgKuGH9N5mtkbSdnWaOQm4BxLdegtoq0bNzEzSsEYjWfmF/7fhuufVB9NzyMHD\n2Xxdnmnts3K6ibKJu8w7S7MTf1XK9iaeK8kfZQKiv7+f/v7+4W94CBvXZpbZzE8keC4x4MPAl8xs\nTSzzWF5dSTsCbwC+CHyi+V6U0w4hsUzS9mb2SFQlPRrTlwI7JeV2JKwglsbrNH1pUeOT/+1Tw9xd\nx3E6nd7eXnp7eze8nzNnzvA0PIQZVNI44HZgN+DrZjZf0h7AqyWdRlgPf8rMbs2p/p/Ap4HNm+9B\nNdrhu+kq4Ph4fTxwZZJ+jKSJknYBZgLzzewR4GlJB8aN7GOTOo7jOO2jyT0JADNbb2b7EH74Hihp\nr1h6KzM7iCAELq2tJ+mNBIvRBUDL42a32gT2YuAQYFtJDwL/DpwOXCrpBKIJLICZ3SPpUoKObS0w\n28wyVdRsggnsFIIJ7M+K7rn8wu3DRaqb3TT8eXTbJC010UsDzR+5cZurrwzCOjVqXdSTNJYtrRPd\n7B8n7zXwZpOkYqaR2mYgadzOA2v3rbYLZrQrnp26cUeAdWvD+nZ8TxIoKLnecmqon+5JvIjfb7jO\nTHd32XrxQJ2tU9PdoOcoOtldFEypNj8tlxfgaHCgpfobHXkBmvICNaVtpftAWSAlGAimlAVP2qj/\nOV1J01Sm5inzElt2Irrqnk6zeyZ5JtDjS/LzZol0z+eh5DrvRH83UzCD9j8ZXlUws+WSbiBYdS4B\nLo/pt0haL2kbM0vdPRwMHCnpDYSZbHNJF5nZcU2Pow4tFRJm9s6CrMMKyp8GnJaTfhuw8Y604zhO\nOynYe+mdHl4Zc/44OF/StsBaM3tK0hTgcMIP6GeB1wA3RtXTxBoBgZl9BvhMbOcQgkqqJQICxsRR\nEMdxnBbR/Aw6HZgraTxB7X+Jmc2TNAH4tqSFBOXFcQCSZgDnmdkROW0Nq/FPLRrQ6HQ+koy3xPGk\nS+gp8W/6QLdMrpNtccuxE9DsePHsxnmDeKogPf21kfVl04K+bBH/5pmlQjAIhuIoLpkaLd3q3z29\nDqqfHWcMqKO2ZSBA0ZRoTpuqgIriceeRF6M7L+jS0Nsfvhjijdy/rFwjsdXz7l/1vkXqurK2MtPq\n1ONAero/uy6L8Z4fSGswn/jGOQNv/nl0zTOSMLMh6fMlmb29Ytn/Zsj3axe+knAcpzWMMsHQEsbA\nDDoGhug4jtMixsAMOgaG6DiO0yLcC2wHkjlnzfNymaalewK3Us0b1KLkOk/1m5oFpvsX6X0n56Tl\nmSOWua0oenJZemJiS2JlwbSwQbJkp4GNiiXTkk2LzES3qH9VvZim/zwTcvJTE+SeErVET86HnZcG\njIv2qj0T8k2Es+vxPYkeP6et8eMa3ycoSwvpQ9ufqGpCnLaTt79Q5AU5zwtxI6T3aukx4NFC982g\nGzEGhliNvA3rDXlx/00HjExfHMfpEMaAS3QXEo7jOM3i6qYO5PGctDU5aVNy0spITVxTE9W0/e1y\nyqahtTPVS09BfqYlKDMMKTKmy760qVnt4uR6y5z81M/k5Jq/kP8tKfvmlKmbBqmwksHk/dNNyLlZ\nT34H1vcEddrqMnVZUf/L/umLTI+rtl9VXVeUlqfOyyPtZ973v4iqWra8uNy1vK6B+3Yq3TeDbsQY\nGOLIYveEv2q52y3HcdrOGJhBx8AQHcdxWoSrmxzHcZxCxsAM2n1DzPYCinT+Gc14qUzrPNxAvUFe\nPPMUxGnEvDIFclQ2F+1ZrI0DfzzZdHk8UVBn6v8ityCZC5G8fYTa63o0sidRVLbqPfPyG7l/Ub16\nFO1NNLN/U9aPss+kmfbLIgw2Uv/ZgnTfk+gKxsAQ24M9PXCtjvTY4jhOKaMsAl8rcCHhOI7TLGNg\nBh0DQ3Qcx2kRY2AGHQNDdBzHaRFu3eQ4juMUMgZm0HHt7oDjOE7H0lPxVYOkyZJulnSHpLsk9SV5\nJ0paFNPPyKm7k6QbJN0dy3ysFUPLGANy0HEcp0U0qW4ys5WSDjWzFZJ6gJskXQNMBY4E9jazNZK2\ny6m+Bvi4md0haVPgNknXmdminLJDpvuExJrkb3ZUIM8fTXVP0AUkNq404oPjyfg3PbyRd5Ajz9c5\nlJ+jyPMlnhjaW0x/JjlH8czUnPsWHZQoo8y5UQ6piXD2T1fVR1KV/LwzB81+8xvxWVWPMt9KZR9/\n0fiGSjO+m4Z65qKTGYIXWDNbES8nEp64AR8GvmRma2KZx3LqPQI8Eq+flbQImMHgYAbDRlerm+zp\nwecVhtzeQwOvhurZwMtxnC6iSXUTgKRxku4AlgHXmtl8YA/g1ZJ+K6lf0t/Wu72knYF9gZuHYTS5\ndN9KwnEcZ6QYwirOzNYD+0jaArhC0l6EOXkrMztI0gHApcCuefWjquky4CQzezavzHDgQsJxHKdZ\nCmbQ/oXhVQUzWy7pBmAWsAS4PKbfImm9pG3M7Im0jqQJwI+A75nZlU33vwIuJBzHcZqlYAbt3Te8\nMub8cHC+pG2BtWb2lKQpwOHA6QRPWK8BbpS0BzAxR0AIuAC4x8zOHqaRFNLSPYkiUy1JW0u6TtJ9\nkq6VtGVS51RJf5B0r6TXJen7S1oY877Syn47juNUYnzF18ZMB34p6U5gPmFPYh7wbWBXSQuBi4Hj\nACTNkHR1rPtK4D3AoZIWxNes1gwQZC3cTZW0PbB9aqoFvBl4H/C4mX1Z0skEHdwpkvYEfgAcAOwA\n/AKYaWYmaT7wUTObL2ke8FUz+1nN/YyegfFYtBpRBVFo6xsdW751U9nHKS2LV22ybtqQnobmc+um\nyrh1U6CCdVMWgGs0IgkzG5LrTUlmN1YsewhDvl+7aOlKwsweMbM74vWzBBOtHQh2wHNjsbkEwQFw\nFHCxma0xs8XA/cCBkqYDm8Xdf4CLkjqO4zjtofmVRMcwYnsSNaZa08ws+0m9DJgWr2cAv02qLSEI\nlTXxOmNpTHccx2kfY2BXd0SGGFVNPyKYaj2jJMBCVCUNn85rXd+Gy/7+Xnp7e/PLjaozC6m6KXsk\n6WnAlOwwXtEJpqpR78v0GUPV9zQQlceSslnQpEHDa0SFVTaWoeqbymhC3VZK3vdgFJF+lFsWlmor\n/f399Pf3D3/DLiSGTmKq9d3EVGuZpO3N7JGoSno0pi8Fdkqq70hYQSyN12n60twbju/bcJnJh7z9\nhqEGAjJr5JR1Wm/aRmnSkzklwWxmTtllOSUdx6lHb+/gH4xz5swZnobHgJBotXVTkanWVcDx8fp4\n4Mok/RhJEyXtAswE5sdj6E9LOjC2eWxSx3Ecpz34nsSQyUy1fidpQUw7lWAPfKmkE4DFwDsAzOwe\nSZcC9xAUDrNtwPxqNnAhYe09r9ayyXEcZ8QZAyuJlprAjjRFJrD5ZQe/b+fHIN2evBv41pntnVO2\nP1510Z7EkPvSSPu+JzHsFOxJbOyabvQwbCawf6pYdtfONYEdA3LQcRynRYyBGXQMDNFxHKdFjIEZ\ndAwM0XEcp0WMgRl0DAzRcRynNViHWy5VwYWE4zhOk6wbAzPoGBii4zhOa3Ah4TiO4xSyatLEiiVX\nt7QfrcSFhOM4TpOsG9/9mxIuJBzHcZpkXaf73KhAS303OY7jdDNrGV/pVYukyZJulnRHjNrZl+Sd\nKGlRTD8j776SZsXonX+Igdtahq8kHMdxmmRdk1Ooma2UdKiZrZDUA9wk6RpCmMgjgb3NbI2k7Wrr\nShoPfB04jOAN+xZJV5nZoqYHUgcXEo7jOE0yFHWTma2IlxMJTr8M+DDwJbPgec4s1wPWK4D7Y/RO\nJP2QENXThcTYoMhxXy1pjO08p35V26llqF+JoTq4a4WDvKJ436OZshjlFOSPFEk/1iYOCB8Htt04\nDkq3MhQhIWkccDuwG/B1M5svaQ/g1ZJOA1YCnzKzW2uq7gA8mLxfAhzYdEdK6DohUc/z66ByHer8\n1qwXAOny9nbEcQoYzd5fh5tV5JvAzu9fyfz+lXXrmtl6YB9JWwBXSNqLMCdvZWYHSToAuBTYtbbq\nkDveAJWFhKRNgJVmtq6F/XEcx+kYivYk9u/dlP17N93w/htzlhe2YWbLJd0AzCKsCi6P6bdIWi9p\nGzN7IqlSG8Fzp1ivJRRaN0kaL+ldkq6W9Cjwe+CRuOv+fyXt3qpOOY7jdALrGF/pVYukbSVtGa+n\nAIcT9hSuBF4T0/cAJtYICIBbgZmSdpY0ETiaENWzJdRbSfwSuJ4QSe7ubAUhaRvgUOAMSVea2Xdb\n1TnHcZzhgJCPAAAdRklEQVTRzBD2JKYDc6Ol0jjgEjObJ2kC8G1JCwnHtI8DkDQDOM/MjjCztZI+\nCvycEBz1glZZNkGdyHSSJppZ3bPkkiZku/CjAUnWiZH2BkemG8Bsvzp10j0J37iuj29cDz9pPwZH\nzjMb/RvXwxWZbr69tFLZV+iu7otMlwoISVsBLyAJ6W1mt48mAeE4jjPSNHtOopMoHaGkLwDvBf4E\nrE+yDm1RnxzHcTqCseCWo4oYPBrYrUz15DiOM9ZYXWAC201UERJ3AVsBy1rcF8dxnI4izy9Tt1FF\nSJwG3C7pbmBVTDMzO7J13XIcxxn9+J5E4CLgDMKKItuT6DwTIsdxnGFmtO9JxDMYbwReBcwAngcW\nAleb2d1V2qgiJJ4zs6823UvHcZwuZTQLCUlzgDcB/cBvgceAycAewOlRgHzCzH5Xr50qQuLXkr5E\nONGXqZsws3zj/oEOTgZuBCbF+1xmZn2StgYuAV4ILAbeYWZPxTqnAu8H1gEfM7NrY/r+wIVxgPPM\n7KQK/XYcx2kpo3xPYr6Zfa4g70xJ0xjs3iOXKkJiP4J66aCa9LomsHX8pb8NuM7MvhyDZZwCnCJp\nT4Il1Z4EL4e/kDQzno47FzghekmcJ2mWmf2sQt8dx3FaxmomtbsLhZjZ1QCS3m5m/53mJWmlBkml\nQsIyt6PNdTLPX/qRwCExfS5hKXQKwR/6xfGA3mJJ9wMHSvozsJmZzY91LgLeDLiQcBynrYxmdVPC\nZ4D/rpCWS6GQkHQs8P3ozjYvf3dgupn9uk4bef7Sp5lZJr2WAdPi9QyC3ixjCWFFsYbBHg6XxnTH\ncZy2MprVTZJeD7wB2EHSV4HMLchm5PvyyaXeSmIbYIGCY6FbCeFEJgO7A68GniCsAArJ8Zf+0pp8\nk+SWUsPK8/Fvke+mqt+NTvFx5JQrBJr141X1XpvXvB/9vpuGi1FuAvsQcBtBS3Nbkv4M8PGqjdTz\n3XS2pK8T3Na+EtibMAMtAo41s79UvUniL/0fgGWStjezRyRNBx6NxWp9pO9IWEEsjddp+tKie/X1\n9W247u3tpbe3t2o320Y9R37Fdd6amy5dOMTeOE730d/fT39//7C3O5rVTWZ2J3CnpO8TPM2+mKDy\nv7cRDxqFXmCHiqRtgbVm9lQ0tfo5cDrQCzxhZmdIOgXY0syyjesfEOK37gD8Atg9rjZuBj4GzAeu\nBr6at3HdqV5gh5MBIeEribHD6FpJFP2AGU0MlxfYc+29lcp+RBe2zQuspCOA/yL434MQ6e5DZjav\nSv1WrpWK/KX/FrhU0glEE1gAM7tH0qXAPYRv9exkxp9NMIGdQjCB9U1rx3HazmheSSScBRxqZvcD\nSNoNmBdfpbRMSJjZQoL5bG36k8BhBXVOI7gBqU2/DXjZcPfRcRxnKKwaxSawCU9nAiLyJ+DpqpUL\nw5c6juM49RlC+NLJkm6WdIekuyT1xfQ+SUskLYivWXn3lfTxWG+hpB9Iqietbovny94r6b3AT4Fb\nJb1VUqlusEo8icmEA3A7J+XNzD5fVtdxHKebaVbdVOewsQFnmdlZRXUl7QCcCLzEzFZJugQ4hnDu\nLI/JBAOh7Hxa5p7jTfH95XmVMqqom34MPEUwoVpZobzjOM6YYCjnJAoOG8PAeYZ69ABTJa0DplLH\n4tOs4u56nRuVsYOZ/cNQbuI4jtONDOWcRMFh49cDJ0o6jnA+7ZOZb7sMM1sq6UzgL4RjCT83s1/U\nuc9uwNnA3xEE0f8AHzezPxXVSakywv+RtHeZp0DHcZyxRpG6aXH/n/lz/5/r1s05bLwXwU9dpsr/\nAnAmcEJaT9JWBPdGOwPLgf+W9G4z+37BrX4AfB3I9h+OBi4GDqzbwUg9txwL4+V44H2SHmBw0KG9\nq9zAcRynWykSEjv17spOvbtueP+rOTcVtpEcNp5lZmdm6ZLOB36SU+Uw4AEzeyKWuxw4GCgSElPM\n7LvJ++9J+nRhh2qot5LINjWMjXVkY/vEmuM4DrCqyRjXOYeNDyfEeNjezB6Jxd5CCBBUy5+Bg2K9\nlQShMT+nXMY1MQzDxfH90TFta9hwLKGQem45FsfBfNfMjq0Z4HeBY/PqOY7jjBWGsCdRdNj4Ikn7\nEH6IPwB8CEDSDOA8Mzsi7l1cRtjPWBv/fqvOvY6O7X2wIH3XjWokVBnhIKd80Vxr/wr1HMdxupoh\nmMAWHTY+rqD8Q8ARyfs+oK/ivXZupo8ZhYfpJH1G0jPAyyQ9k70I9rZXDeWmjuM43UCzh+lGAkmv\nKsnfXFKpJ4t66qbTgNMknW5mdV2CO47jjEVGczwJ4G2SvgxcQzjn9hjB/93uBEerOwOfLGukirrp\nVElvA/4eWA/cZGZXNNdnx3Gc7mE0x5Mws/8jaRuCx4y3E/ZBsnAP3zSzYpOrhCojPIdw2ONigpXT\nhyUdbmazm+q5M0KkLsHzXEWXuQx/viTfaS1Dde9dOfBYg+S5kH+U8ON07DHavcBGM9lvUX9juy5V\nhMShwJ5ZGFOFgAX3NHtDp7VkJ/Clc9vbEWdM0QkxJFrB6iZNYDuJKkLifuAFhNgPxOv7C0s7juOM\nEUb5nsSwUEVIbA4skjSfYFP7CuAWST8hnLw+spUddBzHGa2M5j2J4aLKCP+9Tp6fvHYcZ8wymvck\nJL0CeNDMHo7vjydsYi8G+spOWmeUCgkz65e0MyHe9C8kTQV6zKxyZCPHcZxuZDQLCeCbwGsBJL0a\nOB34KLAvYSP7H6s0UiXo0AeBDwBbE6ycdiR4KnxtM712HMfpFkb5nsS4ZLVwNMHs9UfAjyTdWbmR\nCmX+mXBG4mkAM7sP+JsGO+s4jtN1rKOn0qtNjJeU2SwfBtyQ5FXuVJWCq2KIPGCD7ybfi3AcZ8wz\nyk1gLwZulPQ4sAL4NYCkmYRoo5WoIiRulPRZQqi8w4HZ5Ps4dxzHGVOMZnWTmX1R0i+B7YFrs7Nu\nhEPRJ1Ztp4qQOIUQGWkhwW3tPOD8xrrrOI7TfYx2E1gz+9+ctPsaaaOKddM6SVcCV5rZo4007jiO\n082McuumYaGeq3BJ6ov6rN8Dv5f0uKTPKdugcBzHGcOMZlfhw0U966aPA68EDjCzrcxsK8Jp61fG\nPMdxnDFNs0JC0mRJN0u6Q9Jdkvpiep+kJZIWxNesvPtK2lLSZZIWSbpH0kGtGmM9IXEc8C4zeyBL\nMLM/Ae+OeZWQND4O9ifx/daSrpN0n6RrJW2ZlD1V0h8k3SvpdUn6/pIWxryvNDJAx3GcVrGKSZVe\ntZjZSuBQM9sH2AeYJelAguXoWWa2b3z9rODWXwHmmdlLgL0J7r9bQj0h0WNmj9UmxrRGdmtOIniN\nzcxmTwGuM7M9gOvjeyTtSTjwsScwCzgnUWudC5xgZjOBmUXS1XEcZyQZirrJzFbEy4kEH+zZHFlX\nnS9pC+BVZvbt2M5aM1s+PCPamHpCop5D+krO6iXtCLyBYA2VDfxIYG68ngu8OV4fBVxsZmvMbDHB\n0+yBkqYDm5nZ/FjuoqSO4zhO2xiKkJA0TtIdwDKCiWo2x50o6U5JF6SaloRdgMckfUfS7ZLOi+6S\nWkK9FcHeMaZ1HlMqtv+fwKcJnmQzppnZsni9DJgWr2cAv03KLQF2IAikJUn60pjuDBuNBLhpVTCb\nsUzR5z8azSv9+acUnZNY2T+flf3zc/My4rmFfeLK4ApJexG0Jp+PRb4AnEk4gpDSA+wHfNTMbpF0\nNkEjU88Za9PUi3E9pC15SW8EHjWzBZJ6C+5hkob19HZfX9+G697eXnp7c2/d9Zh9pHLZEAbXGV2M\nRgHRufT399Pf3z/s7Radk5jQezATeg/e8H75nG8UtmFmyyXdAMwyszOzdEnnk39weQmwxMxuie8v\nI6rtW0Erv4kHA0dKegMwGdhc0neBZZK2N7NHoiopO3uxFNgpqb8j4cNYGq/T9KVFN02FhOM4Dmz8\ng3HOnDnD0m6z5q2StgXWmtlTkqYAhwOnZ3NjLPYWwiHmQcS580FJe8SDcYcBdzc3gnKqOPhrCjP7\njJntZGa7AMcAvzSzY4GrgONjseOBK+P1VcAxkiZK2gWYCcyPH9jTkg6MG9nHJnUcx3HaxhD2JKYD\nv4zeWOcT9iTmAV+W9LuYfgjxuIGkGZKuTuqfCHw/ltsbOK1VYxzJNW2mVjoduFTSCYTgF+8AMLN7\nJF1KsIRaC8w2s6zObOBCwl7IvDpmYY7jOCPGqtXNOfgzs4WEfYXa9NzjBWb2EHBE8v5O4ICmbt4g\nGpiHOx9J1k3jGSka25PwjcvhJ2/jurP2JMw+2+4uNIQkzGxIniMk2abPbXRKIJdnN9luyPdrF531\nTXQcxxlFrFvb2S43quBCwnEcp0lcSDiO4ziFrF3jQsJxHMcpYP267p9Cu3+EjuM4rcLVTY7jOE4h\nK7t/Cu3+ETqO47SKRtyedSguJBzHcZrFhYTjOI5TiAsJx3Ecp5Ax4IDAhYTjOE6zrGt3B1qPCwnH\ncZxmcXWT4zTLGPjvaSlrae7fsxWfu08Thaxsdwdajz99B7N/qVxW+mJSb2PPn2m+MzTqeVYdDZ9z\np3l+bQlj4LdQy4IOOY7jdD1rK75qkDRZ0s2S7pB0l6S+mN4naYmkBfE1q+jWksbHMnkhTocNX0k4\njuM0S5MrCTNbKelQM1shqQe4SdI1hOBsZ5nZWRWaOYkQpG2z5npRDV9JOI7jNMuaiq8czGxFvJwI\nTGAgemdpcCJJOwJvAM6vUn4ouJBwHMdplnUVXzlIGifpDmAZIcb1/Jh1oqQ7JV0gacuCO/8n8Glg\n/TCNpBAXEo7jOM3S5J4EgJmtN7N9gB2BAyXtBZwL7ALsAzwMnFlbT9IbgUfNbAEtXkWA70k4juM0\nT5EJ7O/74b7+Sk2Y2XJJNwCzzGyDUJB0PpC3KX0wcKSkNwCTgc0lXWRmxzXS9aq4kHAcx2mWoo3r\n3XrDK+OncwZlS9oWWGtmT0maAhwOnC5pezN7JBZ7C7Cwtmkz+wzwmdjOIcCnWiUgwIWE4zhO8zR/\nTmI6MFfSeILa/xIzmyfpIkn7EDaxHwA+BCBpBnCemR2R05blpA0bLiQcx3GapXkT2IXAfjnpuSsC\nM3sI2EhAmNmNwI3N9aIaLiQcx3Gaxb3AOo7jOIWMAS+wLTeBlbRY0u/i8fH5MW1rSddJuk/Stakt\nsKRTJf1B0r2SXpek7y9pYcz7Sqv77TiOU8rKiq8OZiTOSRjQa2b7mtkrYtopwHVmtgdwfXyPpD2B\no4E9gVnAOZIyO+BzgRPMbCYws55PE8dxnBFhCOckOoWROkxXe+DjSGBuvJ4LvDleHwVcbGZrzGwx\ncD/hkMl0YLPkROJFSR3HcZz2MAS3HJ3CSK0kfiHpVkkfiGnTzGxZvF4GTIvXM4AlSd0lwA456Utj\nuuM4TvsYgluOTmEkNq5faWYPS9oOuE7SvWmmmZmkYbPz7evr23Dd29tLb2/vcDXtOGOUWn1J59m7\n9Pf309/fP/wNd7gqqQoya+k5jME3kz4HPAt8gLBP8UhUJd1gZi+WdAqAmZ0ey/8M+Bzw51jmJTH9\nncAhZvbhmvZtJMfjFCPNKcwz+1zd/OHC7HMbpY3EfYdKXr+r0q7PtdOQhJkNye+RJOMjFeebc4d+\nv3bRUnWTpKmSNovXmwCvIxwzvwo4PhY7HrgyXl8FHCNpoqRdgJnA/HhM/WlJB8aN7GOTOo7jOO1h\nDOxJtHrdOA24Ihoo9QDfN7NrJd0KXCrpBGAx8A4AM7tH0qWEQBprgdnJ0mA2cCEwBZhnZj9rcd8d\nx3Hqs6rdHWg9LRUSZvYAweVtbfqTwGEFdU4DTstJvw142XD30XEcp2nGwJ5E5+1AOY7jjBY6XJVU\nBRcSjuM4zdLh5q1VcCHhOI7TLK5uchzHcQpxIeE4juMUMgb2JEbKd5PjOE73sariqwZJkyXdLOkO\nSXdJ6ovpfZKWRK/ZC/IcmUraSdINku6OdT/WquGBryQcx3Gap/nIdCslHWpmKyT1ADdJuobg6+4s\nMzurTvU1wMfN7A5JmwK3SbrOzBY115v6uJBwHMdpliGom8xsRbycCExgIFZ1Xfcd0QPFI/H6WUmL\nCE5QWyIkXN3kOI7TLEPwAitpnKQ7CJ6wr01CIZwo6U5JF6QB2Qra2BnYF7h5qEMpwoWE4zhOswwh\n6JCZrTezfYAdCXFz9iIEV9uF4KniYeDMoltHVdNlwElm9uzwDGhjXN3kOI7TLEV7Es/1w4r+Sk2Y\n2XJJNwCzzGyDUJB0PvCTvDqSJgA/Ar5nZi11dupCwnEcp1mK9iQm9oZXxuODXbhL2hZYa2ZPSZoC\nHA6cLmn7uOcA8BaC1+zaugIuAO4xs7OHOIJSXEg4juM0S/NeYKcDcyWNJ6j9LzGzeZIukrQPYRP7\nAeBDAJJmAOeZ2RHAK4H3AL+TtCC2d2qrPGOPaNChVuNBh5yqDCU4z3AFBeqG4D2dyrAFHZpecb55\nuHODDvlKwnEcp1nGwIlrFxKO4zjN4l5gHcdxnELcwZ/jOI5TiAsJx3EcpxDfk3Acx3EKGQMrCXfL\n4TiO4xTiQsJxHMcpxIWE4ziOU4jvSTiO4zRN9+9cu5BwHMdpmu7fuW65uknSlpIuk7RI0j2SDpS0\ntaTrJN0n6do0sIakUyX9QdK9kl6XpO8vaWHM+0qr++04jlPOmoqvzmUk9iS+Aswzs5cAewP3AqcA\n15nZHsD18T2S9gSOBvYEZgHnRLe4EIJxnGBmM4GZeQHCHcdxRpbnK746l5YKCUlbAK8ys28DmNla\nM1sOHAnMjcXmAm+O10cBF5vZGjNbDNxPiNg0HdgsCe93UVLHcRynTfhKYqjsAjwm6TuSbpd0nqRN\ngGlmtiyWWQZMi9czgCVJ/SXADjnpS2O64zhOGxlC/NIOodVCogfYDzjHzPYDniOqljJiAAgPAuE4\nTgfS/SuJVls3LQGWmNkt8f1lwKnAI1mYvqhKejTmLwV2SurvGNtYGq/T9KV5N+zr69tw3dvbS29v\n79BH4XQd7Qr444GG2kN/fz/9/f0taLm5VYKkycCNwCTCPHyZmfVJ6gP+CXgsFs2NOBf3ZM8GxgPn\nm9kZTXWkSl9bHclN0q+AfzKz++IHMDVmPWFmZ0g6BdjSzE6JG9c/AF5BUCf9AtjdzEzSzcDHgPnA\n1cBXaz88j0znOE4Vhi0yHXdWLP3yje4naaqZrZDUA9wEnEQw2HnGzM6qc9/xwO+Bwwg/lm8B3mlm\ni5oZRxkjcU7iROD7kiYCfwTeR5B+l0o6AVgMvAPAzO6RdClwD0FEz05m/dnAhcAUgrVUS+K5Oo7j\nVKd5yyUzWxEvJwITGFC7lwmvVwD3R+MeJP2QYPTTmULCzO4EDsjJOqyg/GnAaTnptwEvG97eOY7j\nDIXmN6UljQNuB3YDvm5m8yW9HjhR0nHArcAnzeypmqo7AA8m75cABzbdkRL8xLXjOE7TFG1K3xFf\nxZjZemCfeFTgCkl7Ec6DfT4W+QJwJnBCbdWmu9sELiQcx3Gapmgl8dL4yphbUA7MbLmkG4BZZnZm\nli7pfOAnOVVqDXx2YvARgWHFvcA6juM0TXMmsJK2zdwRSZoCHA4skrR9UuwtwMKcm95K8Dqxc9zr\nPRq4arhGVIuvJBzHcZqm6T2J6cDcaKk0DrjEzOZJukjSPgSV0gPAhwAkzQDOM7MjzGytpI8CPycY\nAV3QKssmGAET2JHETWAdx6nC8JnA/qhi6bcN+X7twlcSjuM4TdPZzvuq4ELCcRynaTrb5UYVXEg4\njuM0TWc776uCCwnHcZym6f6VhJvAjnJa45Ss/XTjuLpxTNC94xoe3FW402a69R+0G8fVjWOC7h3X\n8OCuwh3HcZxCOnuVUAUXEo7jOE3T/SawXXeYrt19cBynMxiew3Qjd7920VVCwnEcxxlefOPacRzH\nKcSFhOM4jlNIxwsJSV+QdKekBZJ+Lml6kneqpD9IulfS69rZz0aR9H8lLYpjuzwGJiG6B34+jneB\npHPa3deqFI0p5nXys3q7pLslrZO0X5Leyc8qd0wxr2OfVYqkPklLkuczq919GpWYWUe/gM2S6xOB\nc+P1noTQUBOAnYH7gXHt7m8D4zo86y9wOnB6vN4ZWNju/g3zmDr9Wb0Y2AO4AdgvSe/kZ1U0po5+\nVjVj/BzwiXb3Y7S/On4lYWbPJG83BdbH66OAi81sjYWA4fcTAoh3BGZ2nYXwhgA3Azu2sz/DQZ0x\ndfqzutfM7mt3P4aTOmPq6GeVQ0daHI0kHS8kACR9UdJfgHcB/x6TZzA4pN8SQgDxTuT9wLzk/S5x\nedwv6e/b1akhko6pm55VLd3wrFK67VmdGNWfF2SR4pzBdMRhOknXAdvnZH3GzH5iZp8FPivpFILK\nqa+gqVFl71s2rljms8BqM/tBzHsI2MnM/hp1xVdK2qtmRdU2mhxTHh33rHLo+GdVkVH1rFLqjPGz\nwLnA5+P7LwBnAieMUNc6ho4QEmZ2eMWiPwCuJgiJ2mDhO8a0UUPZuCS9F3gD8Nqkzmpgdby+XdIf\ngZnA7a3raXWaGRNd8KwK6nT0sypg1D+rlKpjlHQ+0IhgHDN0vLpJ0szk7VFAFuv1KuAYSRMl7UL4\n55w/0v1rlmhp8WngKDNbmaRvG+PiImlXwrj+1J5eNkbRmOjwZ1XDBh13Jz+rGlK9fdc8q9QSEngL\nsLBdfRnNdMRKooQvSXoRYcN6MfBhADO7R9KlwD0EL1yzLZo0dAhfAyYC10kC+F8zmw0cAsyRtIYw\n5g+Z2VPt62ZD5I6p05+VpLcAXwW2Ba6WtMDMXk8HP6uiMXX6s6rhDEn7ENRlDwAfanN/RiXulsNx\nHMcppOPVTY7jOE7rcCHhOI7jFOJCwnEcxynEhYTjOI5TiAsJx3EcpxAXEo7jOE4hLiQcx3GcQlxI\nOI7jOIW4kHBagqRnW9z+b+LfLSR9pIn6fZI+2WCdKdGb60bupWNwngXJ6+SYPlHSjZl7DsfpNLrB\nLYczOmnpUX4ze2W83AqYTfDo2VATTdz2/cCPgO0lHQOsAvrN7B5ghZntm9PP1ZKuB44mOKB0nI7C\nVxLOiCHpE5IWxtdJMW3nGNL0W5LuUghBOzmp828xTOavJV2c/fpPViqnA7vFX+9nSHqhpIVJ/U9J\n+ly8/qyk30v6NfCipMx7JN0c2/gvSUX/F+8CfkwIsjMX+CZwQIWhXwm8u+LH5DijChcSzoggaX/g\nvYQJ9iDgA9G5GsDuwNfN7KXAU8DbYp0DgLcCewOvB/ZnYAWQ/T0Z+KOZ7WtmJ7NxpDFL7n808HKC\nq/IDAJP0EuAdwMFxJbCenAld0kRgVzP7C8Ed/RuB9xBWFgBTatRNb0+q3001YeI4ow5XNzkjxd8D\nl5vZ8wCSLgdeRXA9/YCZ/S6Wu40QOxnglcCVWVwGSXn+/quGn8zuvxJYKemqWPc1BOFza9xqmAI8\nklN/W4IAw8zWAhfV5D+fp26K5ddJWi1pEzN7rmJ/HWdU4ELCGSmMwRO6GFgNrErS1xEm6qI6Zaxl\n8Ap5SnKd15aAuWb2mZJ2nwcml5SpxyRgZWkpxxlluLrJGSluAt4cLYQ2Ad4M/Jr6E/9vgDdJmiRp\nU+CInDLPAJsl75cBfyNpa0mTCGohA34V7z9Z0mZJ+vXAP0raDiDWe0HtTczsr8D4qHZqCEnbAI+b\n2bpG6zpOu/GVhNMqpkp6MHl/FnAhA1HMzjOzOyXtzMaWRgZgZrdGtdDvCJP/QmB5TZknJP0mblbP\nM7OTJX0+3mcpITgOZrZA0iXAncCjWT/MbJGkfwWujRvWawjWUn/JGdO1BBXZ9Tl5UyQtSN5fk6xO\nDgV+mlPHcUY9HnTIGdVkenxJU4EbgQ+Y2R1t6su+wMfN7LgG6/0IONnM7m9NzxyndfhKwhntfEvS\nnoT9gAvbJSBgw2rkBknjzGx9lTqSJhA2311AOB2JryQcx3GcQnzj2nEcxynEhYTjOI5TiAsJx3Ec\npxAXEo7jOE4hLiQcx3GcQlxIOI7jOIW4kHAcx3EKcSHhOI7jFPL/Aa0aFXJapoOOAAAAAElFTkSu\nQmCC\n",
"text/plain": [
"<matplotlib.figure.Figure at 0x7f35d10312e8>"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"m1 = netCDF4.Dataset('/archive/Niki.Zadeh/ulm_201505_awg_v20151106_mom6sis2_2015.11.30/CM4_c96L32_am4g7_2000_sis2_lowmix_intel15_c3/gfdl.ncrc3-intel15-prod-openmp/pp/ocean_annual_z/av/annual_20yr/ocean_annual_z.0021-0040.ann.nc')\n",
"for v in m1.variables:\n",
" print(v, end=',')\n",
"x = m1.variables['xh']\n",
"y = m1.variables['yh']\n",
"z = m1.variables['zt']\n",
"j = find_index(y, 38.0)\n",
"i0 = find_index(x, -30.0);\n",
"i1 = find_index(x, -5.0)\n",
"q = m1.variables['so'][0,:,j,i0:i1]\n",
"plt.pcolormesh(x[i0:i1],z[:],q);\n",
"plt.gca().invert_yaxis();\n",
"cb = plt.colorbar(); cb.set_label('S (ppt)')\n",
"plt.xlabel('Longitude ($\\degree$E)');\n",
"plt.ylabel('Depth (m)')\n",
"plt.title('38$\\degree$N, yrs 21-40, CM4_lowmix z*')\n",
"plt.clim(34.9, 36.7); plt.xlim(-30.,-5.); plt.ylim(6500, 0);"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {
"collapsed": true
},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"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.4.3"
}
},
"nbformat": 4,
"nbformat_minor": 0
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment