{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "This is a notebook to experiment with calling a python wrapper for the teos-10 matlab library.\n", "\n", "The python module is gsw_calls.py, currently stored in /data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC, along with the matlab wrappers for the gsw routines. \n", "\n", "The function call_generic_gsw() links the appropriate matlab scripts to the current directory, saves all the input files to individual text files, then calls the matlab wrapper script which read the input files and saves the ouptut to a text file. Python then reads the output file, rehsapes and returns as a numpy array. It also removes the links and deletes the temporary text files. " ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "collapsed": true }, "outputs": [], "source": [ "import sys\n", "sys.path.append('/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/')\n", "import gsw_calls\n", "\n", "import netCDF4 as nc\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "%matplotlib inline" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "collapsed": true }, "outputs": [], "source": [ "F = nc.Dataset('/data/nsoontie/MEOPAR/NEMO-forcing/open_boundaries/west/SalishSea2_Masson_corrected.nc', 'r')" ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "collapsed": true }, "outputs": [], "source": [ "sal = F.variables['vosaline']\n", "temp = F.variables['votemper']\n", "dep = np.expand_dims(np.expand_dims(np.expand_dims(F.variables['deptht'][:],axis=0),axis=2),axis=3) \\\n", " + np.zeros(sal.shape)\n", "long = F.variables['nav_lon'][:] + np.zeros(sal[:].shape)\n", "lat = F.variables['nav_lat'][:] + np.zeros(sal[:].shape)\n" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mw_gsw_p_from_z('output_file','input0.txt','input1.txt');exit\n" ] } ], "source": [ "p = gsw_calls.generic_gsw_caller('mw_gsw_p_from_z.m', [-dep, lat])" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "(52, 40, 1, 870)" ] }, "execution_count": 5, "metadata": {}, "output_type": "execute_result" } ], "source": [ "p.shape" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAV4AAAEACAYAAAD7rx6dAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAGr1JREFUeJzt3X+wXWV97/H3JwlJiSimYJKaxACiGJhrI2qwN+1wrBgD\ntgl17nAV516Qar2XWhjtWBJuO4Fe50KcMYz9kWkLSAPFYsSLBEckyeCuF6/8MqQETkhPqwkQyOFX\nG2u9jfnxvX+sdczmuM/Za6+19jor63xeM2vY+9nP2t9nb06+5znPetbzKCIwM7PqTJnoBpiZTTZO\nvGZmFXPiNTOrmBOvmVnFnHjNzCrmxGtmVrHMiVfSFEnbJG1Kn8+StFnSLkn3STqxf800M2uOXnq8\nVwKDbc9XAVsj4gzgfmB1mQ0zM2uqTIlX0nzgAuCmtuKVwIb08QbgwnKbZmbWTFl7vDcAnwXab3Ob\nExHDABGxD5hdctvMzBqpa+KV9EFgOCK2Axqnqu89NjPLYFqGOkuBFZIuAI4HXivpNmCfpDkRMSxp\nLvBCp5MlOSGbWWYRMV4Hr6vXS7E/e/U9EXFKkXh5qJdFciSdC/x+RKyQ9Hng5YhYK+kqYFZErOpw\nTlxTWnOz+zbw3kkQ03GbHbdJn3VNhlwjqXDilRSfy1j3Dyme6PPI0uMdy/XARkmXAXuAi8ppkplZ\nMcdNdAO66CnxRsTfAX+XPn4FOK8fjTIzK6JIj7IKdW9fbqdMkpiO2+y4ExFzIuOW5fiJbkAXjU28\np06SmI7b7LiT6bOWqVFDDXl5PMLMqlT3HmXd22dm1jP3eM3MKlb3xFb39pmZ9cw9XjOzijnxAkvf\nVEUUM7OEp5OZmVWs7omt7u0zM+uZhxrMzCpW98RW9/aZmfXMPV4zs4rVPbFV075zKoliZgbUv8fb\nyy7DZmbHhOMzHmORNEXSNkmb0uezJG2WtEvSfZJObKu7WtKQpJ2SlmVpnxOvmTXOcRmPcVwJDLY9\nXwVsjYgzgPuB1QCSziTZBGIRcD6wXlLXHS2ceM2scaZlPDqRNB+4ALiprXglsCF9vAG4MH28Argj\nIg5FxG5gCFiSpX1mZo1yXNbMdqhj6Q3AZ4ET28rmRMQwQETskzQ7LZ8HfK+t3t60bFzVJN5frSSK\nmRkA03ImXkkfBIYjYrukgXHOLLR7etfmSZoBfAeYnta/MyKulbQG+ARHt3W/OiK+VaQxZmZlOG5q\n5/LvHE6OcSwFVki6gOT622sl3QbskzQnIoYlzeVo3tsLLGg7f35aNq5M27tLmhkRP5E0FfgucAXJ\nQPK/RsS6LudGfLFrCDNruiuq2979pyd2rwcwff/Y27tLOhf4/YhYIenzwMsRsVbSVcCsiFiVXly7\nnWTS7DxgC/CW6JJYM3XII+In6cMZ6Tkjb1r5fvRmZt0cN6P0t7we2CjpMmAPyUwGImJQ0kaSGRAH\ngcu7JV3I3uOdAnwfeDPw5xGxOh1quBTYDzxK8pthf4dz3eM1s0p7vPHGjHWfG7vH20+ZppNFxJGI\neAfJ+MWStHu9HjgtIhYD+4BxhxzMzCpTZD5ZBXoKHRE/ktQClo8a270RuGes86554ujjgXfCwLt6\na6SZNVOr1aLVapX/xjWfKNt1qEHSycDBiNgv6XjgPpLxjm0RsS+t82ng3RFxcYfzIx4tv+Fmdox5\nZ4VDDW/JWHdoYoYasvxe+CVgQzrOOwX4SkR8U9KtkhYDR4DdwCf710wzsx7UvMfbtXkRsQM4u0P5\nf+1Li8zMiip/VkOpav57wcwsh5pntpo3z8wsh5pntkqa99I7T6gijJnV2MlVBhvjluG6qPnvBTOz\nHGqe2WrePDOzHGqe2WrePDOzHGqe2WrePDOzHDydDB7F9wibTXbLqwxW8y5lzZtnZpaDZzWYmVWs\n5pmt5s0zM8uh5pnN27ubWfNMzXiMImmGpIckPSZpR7rhA5LWSHpW0rb0WN52zmpJQ5J2SlqWpXk1\n/71gZpZDzswWEQckvbd9j0lJ96Yvrxu9x6SkRSTbAC0i2Shiq6Ry9lwr6lHeWUUYM6uxSmc1/EL+\nU3vcY3IlcEdEHAJ2SxoClgAPjRfDQw1m1jw5hxog2WNS0mMkW5ptiYhH0pc+JWm7pJskjexjPA94\npu30vWnZuJx4zax5Cuy5lnGPyS8UbZ6ZWbOMkdlaTydHFhn3mNwLLGh7bX5alqd5ZmbHsDGGEQZO\nTY4R13731a932GPy/cD1kuaO7DEJfAgY2cJ3E3C7pBtIhhhOBx7u1rxKEu923lFFGDOzRP7M1tMe\nkxExKGkjMAgcBC7vNqMhU/MkzQC+A0xP698ZEddKmgV8BViYNuSiiNjf88c0Mytb/ulkPe8xGRHX\nAdf1EqfrxbWIOAC8Nx1sXgycL2kJsArYGhFnAPcDq3sJbGbWNzMyHhMk06yGMea1rQQ2pOUbgAtL\nb52ZWR4FZjVUIVPiHWNe25yIGAZIB51n96+ZZmY9qHnizRQ6Io4A75D0OuAuSWdx9G6On1Ub6/zB\na772s8dvGFjEGwbOzNFUM2uaVqtFq9Uq/41rviykMlyAe/UJ0h8BPwE+DgxExLCkucC3I2JRh/rx\n5nhidLGZTTL/yFld60giIjrdmpuZpIi/yFj3v1E4Xh5dhxoknTxye1zbvLadJPPXLk2rXQLc3ac2\nmpn1pgFDDWPNa3sQ2CjpMmAPyQo9ZmYTr+ZDDV0T7zjz2l4BzutHo8zMCimwOlkVfMuwmTVPzTNb\nJc37p8Hug+pm1nBVTmY61ocazMyOOTXPbDVvnplZDjXPbDVvnplZDh5qMDOrmGc1mJlVzD1e4IFK\nophZnVU5q6HmXUpvdmlmzZPzlmFJMyQ9JOkxSTskrUnLZ0naLGmXpPvadhlG0mpJQ5J2SlqWpXlO\nvGbWPDkTb68bP6Q7EF8ELALOB9ZL6rrojhOvmTXP1IxHBz1u/LACuCMiDkXEbmAIWNKteU68ZtY8\nBVYn63Hjh3nAM22n703LujbPzKxZCuynVnTjhyyqSbwPVhLFzOrsdyqMNUZma30vObKIiB9JagHL\ngWFJc9o2fnghrbYXWNB22vy0bFw970DRK0nBx/obw8zqL77UvU5pO1A8l7HuG1+9A4Wkk4GDEbE/\n3fjhPuB64FzglYhYK+kqYFZErEovrt0OnEMyxLAFeEt0SaweajCz5smf2Xra+CEiBiVtBAaBg8Dl\n3ZIuuMdrZhWpssd75OVsdaecNDF7rrnHa2aNc7jmma2a5m2vJIqZGdCAxCtpPnArMAc4AvxVRPxp\neivdJzh6de/qiPhW31pqZpbRgRnTM9b8aV/bMZYsvxcOAZ+JiO2STgC+L2lL+tq6iFjXv+aZmfXu\n8NR6L0+WZZfhfSR3cBARP5a0k6N3ZlQ+KG1m1s3hmq8L2dMtw5JOIVk44qG06FOStku6qX21HjOz\niXSIqZmOiZJ5CDodZrgTuDLt+a4H/jgiQtLngHXAb3c8+flrjj4+YQBeO5C7wWbWHK1Wi1arVfr7\nHq75hK1M83glTQO+AdwbEV/s8PpC4J6IeHuH14JpnsdrNtnFwe51yprHuydmd68ILNQLtZ7H+yVg\nsD3pSpqbjv8CfAh4ouzGmZnlUfcx3izTyZYCHwV2pEulBXA1cLGkxSRTzHYDn+xjO83MMjtA1ulk\nEyPLrIbv0nnJYM/ZNbNaqvsYb71bZ2aWwzE/1FCKQz+qJIyZ1dnrKovkxGtmVrGJnKObhROvmTWO\nx3jNzCpW96EG7zJsZo3zU6ZnOkaTNF/S/ZKelLRD0u+l5WskPStpW3osbztntaQhSTslLcvSPvd4\nzaxxCozx9rQao6RFJNsALSLZ6HKrpLrsudZ1000za7wqZzXkS205VmNcCdwREYeA3ZKGgCUcXUis\nIw81mFnjHGZqpmM8GVdjnAc803baXo4m6jE58ZpZ4xRNvKNXYwTWA6dFxGKSHvEXirTPY7xm1jhj\njfE+0XqZJ1uvjHtuuhrjncBtEXE3QES82FblRuCe9PFeYEHba/PJMLZazfbuDPY1hpnVX8SirnXK\nWhbyb+PCTHU/oq//XDxJtwIvRcRn2sp+thqjpE8D746IiyWdCdwOnEMyxLAFqMvFNTOz6uSdx9vr\naowRMShpI0nv8iBwebekC5UlXq/VYGbVyTudLM9qjBFxHXBdL3Hc4zWzxvEtw2ZmFav7LcNOvGbW\nOE68ZmYVc+IF4F+rCWNmBhxgxkQ3YVxd71zrsFrPFWn5LEmbJe2SdF/bLXRmZhOqjFuG+ynLLcMj\nq/WcBfwK8LuS3gasArZGxBnA/cDq/jXTzCy7uifeLLsMd1qtZz7JqjznptU2AC2SZGxmNqEatfVP\n22o9DwJzImIYkuQsaXbprTMzy6Ex83hHr9aTrMHwKuPcJndr2+NfTg8zm+xarRatVqv09637rIZM\ni+Skq/V8A7g3Ir6Ylu0EBiJiWNJc4NvRYRWMJEF/s+Rmm9mxJuL8rnXKWiTn6vijTHX/l/5n4Xh5\nZF2P90vA4EjSTW0CLk0fXwLcXWK7zMxyO8D0TMdE6TrUMM5qPWuBjZIuA/aQ7DtkZjbhjvkx3nFW\n6wE4r9zmmJkVV/cx3nr/WjAzy8GJF0juwTAzq0bd5/F6s0sza5zDTMt0jJZniQRJqyUNSdopaVmW\n9jnxmlnjFLhluKclEtI91y4CFgHnA+sldZ2e5jFeM2ucn+acKpZjiYQVwB0RcQjYLWkIWAI8NF4c\nJ14za5wyxngzLpEwD/he22l707JxOfGaWeMUncdbbImE7ipKvAerCWNmxtjTyV5sDfJSa3Dcc9Ml\nEu4EbouIkTtyhyXNaVsi4YW0fC+woO30+WnZ+DGyrNVQRPKb4mt9jWFm9Rfxoa51ylqr4TdjY6a6\n9+iin4sn6VbgpYj4TFvZWuCViFgr6SpgVkSsSi+u3Q6cQzLEsAV4S3RJrB5qMLPGyTvG2+sSCREx\nKGkjMEjyp/3l3ZIuOPGaWQPlHePNs0RCRFwHXNdLHCdeM2ucvNPJquKLa2bWOHW/Zdg9XjNrnGN+\nWUgzs2ONVyczM6uYE6+ZWcWceM3MKnaAGRPdhHF5IXQza5y693i7rscr6WZJw5IebytbI+lZSdvS\nY3l/m2lmll2B9XgrkWUh9FuAD3QoXxcRZ6fHt0pul5lZboeYmumYKFl2GX5A0sIOLxVayMLMrF/q\nPo+3yNY/n5K0XdJN7fsPmZlNtLoPNeT9tbAe+OOICEmfA9YBvz129fZlIRcBZ+YMa2ZN0mq1aLVa\npb9v3S+uZVqPNx1quCci3t7La+nrAX9TuKFmdmyL+GjXOmWtx3vigecz1d0/45cKx8sja49XtI3p\nSpqbbgoH8CHgibIbZmaW1+FD9R7j7do6SV8GBoCTJD0NrAHeK2kxcATYDXyyj200M+vJ4UP1HmrI\nMqvh4g7Ft/ShLWZmpSiSeCXdDPwGMDwyhCppDfAJju61dvXINFpJq4HLSO4UuzIiNneLUe/+uJlZ\nDocOFurx3gL8KXDrqPJ1EbGuvUDSIpJtgBaRbHS5VVJd9lzzQuhmVp0jh/Onth7vXVgJ3BERh4Dd\nkoaAJcBD48UoMo/XzKyeDk3NdvSm070L84Bn2ursTcvG5cRrZs3z79OyHdmtB06LiMXAPuALRZrn\nMV4za56xFkR8uAWPtHp+u4h4se3pjcA96eO9wIK21+anZeNy4jWz5hkr8Z49kBwj1l871jtkvXdh\nE3C7pBtIhhhOBx7u1jwnXjNrngJLgPdy70JEDEraCAySzCK4vNuMBsh4y3ARyS3DN/Y1hpnVX8TH\nu9Yp65ZhHsyY195TPF4e7vGaWfMcnugGjM+J18yap+a7jTnxmlnz/PtEN2B8Trxm1jzu8ZqZVcyJ\n18ysYk68ZmYVq/m6XE68ZtY8nk5mZlYxDzWYmVXM08mg9gMuZtYsNe/xdl2PV9LNkoYlPd5WNkvS\nZkm7JN3XtiiwmdnEO5TxmCBZFkK/BfjAqLJVwNaIOAO4H1hddsPMzHI71hNvRDwA/POo4pXAhvTx\nBuDCkttlZpbfwYzHBMk7xjs7IoYBImKfpNkltsnMrJhJMp2sy+KX32h7/Nb0MLPJrtVq0Wq1yn/j\nms9qyLQQerrV8T0R8fb0+U5gICKGJc0Fvh0Ri8Y4N5J94sxsMov4713rlLYQ+uqMC6FfNzELoWfd\nZfhV+w+R7DN0afr4EuDuEttkZlZMgTHeXmdySVotaUjSTknLsjQvy3SyLwP/F3irpKclfQy4Hni/\npF3A+9LnZmb1cDjj0VnmmVySzgQuAhYB5wPrJXXtQXcd442Ii8d46bxu55qZTYgCU8Ui4oF0eLXd\nSuDc9PEGoEWSjFcAd0TEIWC3pCFgCfDQeDF8y7CZNU/5c3THmsk1D/heW729adm4Kkq8Nb9/z8ya\nZaw5us+14PlWGREKbc/uHq+ZNc+BMcpPGkiOEY9dm/UdhyXNaZvJ9UJavhdY0FZvflo2rqyzGszM\njh3FbxnOOpNrE/BhSdMlnQqcDjzcrXnu8ZpZ8xS4HTidyTUAnCTpaWANycytr0q6DNhDMpOBiBiU\ntBEYTKNeHhlujsh0A0URyQ0Uf9LXGGZWfxG/17VOaTdQ/FbGvHbXxNxA4R6vmTVPza/nO/GaWfM4\n8ZqZVazmm9448ZpZ84w1nawmnHjNrHk81GBmVjEPNUDtvwUza5ZJsgOFmVl9eKjBzKxiTrxmZhWr\n+eimE6+ZNY+nk5mZVazJQw2SdgP7gSPAwYhY0rlmzb8FM2uWhg81HCHZ5v2fy2iMmVkpGj6dTHgx\ndTOrm5r/kV00aQawRdIjkj5RRoPMzAorvgNFXxXt8S6NiOclvYEkAe+MiAfKaJiZWW4Fx3g7Xb+S\nNAv4CrAQ2A1cFBH787x/ocQbEc+n/31R0l0k+8l3SLxb2h6fBry5SFgza4hWq0Wr1Sr/jYv3Zjtd\nv1oFbI2Iz0u6ClidlvUs99Y/kmYCUyLix5JeA2wGro2IzaPqBazNFcPMmiPiD7rWKW3rn8y7r3eO\nJ+mHwLsi4uW2sqeAc9t2Gm5FxNvytLFIj3cOcFfyIZkG3D466ZqZHaNGrl8dBv4yIm4C5kTEMEBE\n7JM0O++b5068EfFDYHHe883MqtdKj67ar19tlrSLn+9G594p2HeumVkDjXV1bWl6jLi2Y61R16++\nTnL9aljSnLahhhfyts5zcM2sgfLPJ5M0U9IJ6ePXAMuAHcAm4NK02iXA3Xlb54XQzayBCuWcjtev\nJD0KbJR0GbAHuChvAA81mFkD/b/cZ451/SoiXgHOK9Con3HiNbMGqvdf2U68ZtZA9V6swYnXzBrI\nPV4zs4q5x2tmVjH3eM3MKpZ/VkMVnHjNrIE81GBmVjEPNZiZVcw9Xur+JZhZ07jHa2ZWsXp39px4\nzayB3OM1M6uYp5OZmVXMPV4zs4rVe4y30A4UkpZLekrSP6TbHZuZ1cDBjMfEyJ14JU0B/gz4AHAW\n8BFJubY67o8fTpKYjtvsuJPps5Yp/9Y/VSjS410CDEXEnog4CNwBrCynWWXYPUliOm6z405EzImM\nW5Z693iLjPHOA55pe/4sSTI2M5tg9R7j9cU1M2ugek8nU0TkO1F6D3BNRCxPn68CIiLWjqqXL4CZ\nTUoRoSLnS9oNLMxYfU9EnFIkXh5FEu9UYBfwPuB54GHgIxGxs7zmmZk1T+6hhog4LOlTwGaSi3Q3\nO+mamXWXu8drZmb5FLqBYjz9vLlC0s2ShiU93lY2S9JmSbsk3SfpxLbXVksakrRT0rICcedLul/S\nk5J2SLqi37ElzZD0kKTH0phrqvq86XtNkbRN0qaq4kraLenv08/8cIVxT5T01fR9npR0Tr/jSnpr\n+jm3pf/dL+mKCuJ+WtITkh6XdLuk6RV9x1emP8eV/PuprYgo/SBJ6P9IMsB9HLAdeFuJ7/+rwGLg\n8baytcAfpI+vAq5PH58JPEYyrHJK2i7ljDsXWJw+PoFkjPtt/Y4NzEz/OxV4kGTaXt8/b/p+nwb+\nBthU4ff8A2DWqLIq4v418LH08TTgxKq+57Z/N88BC/oZF3hj+h1PT59/Bbikgp/js4DHgRnpz/Jm\n4M1Vfsd1OfrzpvAe4N6256uAq0qOsZBXJ96ngDnp47nAU51iA/cC55TUhq8D51UVG5gJPAq8u4qY\nwHxgCzDA0cRbRdwfAieNKutrXOB1wD91KK/s5wpYBvyffsclSbx7gFlpUttUxc8x8J+AG9ue/yHw\nWWBn1f92J/ro11BDp5sr5vUp1ojZETEMEBH7gNljtGVvGW2RdApJr/tBkh+avsVO/9x/DNgHbImI\nR/odM3UDyT+M9gsBVcQNYIukRyR9vKK4pwIvSbol/bP/ryTNrCBuu/8MfDl93Le4EfEc8AXg6fT8\n/RGxtZ8xU08Av5YOLcwELiDp3Vf5HddC38Z4a6BvVw0lnQDcCVwZET/uEKvU2BFxJCLeQdIDXSLp\nrH7HlPRBYDgitgPjzavsx/e8NCLOJvmH+buSfq1DnLLjTgPOBv48jf1vJD2ufscFQNJxwArgq2PE\nKS2upNeT3N6/kKT3+xpJH+1nTICIeIpkWGEL8E2SYYTDnaqWGbeO+pV49wJvans+Py3rp2FJcwAk\nzQVeaGvLgrLaImkaSdK9LSLurjJ2RPwIaAHLK4i5FFgh6QfA3wK/Luk2YF+/P2tEPJ/+90WS4Zwl\n9P/zPgs8ExGPps+/RpKIK/l/C5wPfD8iXkqf9zPuecAPIuKViDgM3AX8xz7HBCAibomId0XEAPAv\nJNdJqvqOa6NfifcR4HRJCyVNBz5MMo5UJvHqntgm4NL08SXA3W3lH06v2p4KnE5ys0deXwIGI+KL\nVcSWdPLIVV5JxwPvJxkT6+vnjYirI+JNEXEayf+/+yPivwD39DOupJnpXxRIeg3JuOcO+v95h4Fn\nJL01LXof8GS/47b5CMkvuBH9jPs08B5JvyBJJJ91sM8xAZD0hvS/bwJ+i2RoparvuD76NXhM0ivb\nBQwBq0p+7y+TXP09QPJD9DGSCwVb05ibgde31V9NckV0J7CsQNylJH8abSf5M2lb+jl/sV+xgf+Q\nxtlOckX4f6TlfYvZoQ3ncvTiWl/jkoy1jny/O0Z+dqr4vMAvk3QatgP/m2RWQxVxZwIvAq9tK+v3\n97wmPf9xYAPJ7KMqPut3SMZ6HwMGqv5ZrsvhGyjMzCrW5ItrZma15MRrZlYxJ14zs4o58ZqZVcyJ\n18ysYk68ZmYVc+I1M6uYE6+ZWcX+P/US2IX4Ha/nAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(p[0,:,0,:])\n", "plt.colorbar()" ] }, { "cell_type": "code", "execution_count": 7, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mw_gsw_SR_from_SP('output_file','input0.txt');exit\n" ] } ], "source": [ "sal_ref = gsw_calls.generic_gsw_caller('mw_gsw_SR_from_SP.m', [sal[:],])" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWEAAAEACAYAAABiV8coAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXuUZFd93/v5TtV0dfc8NBrJGkkIJMxDgGJHIrYg1sW0\nQNiD0AKcOATjlSDw9XKWzYWFswhC9r0wuU6WYC1wSIJvgsBEcCEGCwOCCPRYosEoV0IgjTVYg1CM\n9UCaGRBiRjPqZ1X/7h/n7Kpdu05Vne7qrpou/T5r1eqqU2fv7977nN59ev8eW2aG4ziOMxq2jLoB\njuM4T2d8EnYcxxkhPgk7juOMEJ+EHcdxRohPwo7jOCPEJ2HHcZwRUnoSlrRF0t2Sbsg/nyrpZkn3\nS7pJ0ikb10zHcZzxZDVPwu8A7os+XwXcambnA7cB71nPhjmO4zwdKDUJSzoHuBz4WHT4dcB1+fvr\ngNevb9Mcx3HGn7JPwn8KvAuIw+v2mNkRADM7DJyxzm1zHMcZe/pOwpJeAxwxs/2Aepzq8c+O4zir\npFrinEuA10q6HJgCdkj6FHBY0h4zOyLpTODHRYUl+eTsOE5pzKzXw15fdkl2rPzpD5nZeYPoDYpW\nk8BH0suBf21mr5X0AeCnZvZ+Se8GTjWzqwrK2AdL1D1VuhWdFM3+Xwcuzd9Xga25Rpm/OqvRifkG\ncEWuNYjOE32+D/WHn1/JdTdaJ6WsbhmdcP3LjFs33eN9yoW+lNXppju/hrJrYYrivu7Mv5ui1Z+1\n8mCX4zcA/zzXisdtEC4rMddIGngSlmR/UvLcP2bwSX9QBpkrrgE+J+mtwEPAG9anSY7jOIOxHn80\nhsWqJmEz+wbZAx9m9gRw2UY0ynEcZxAGebocNpupravivBFoPnsEmgDPd92x1R1VX88fke56Mcjy\n5rAZ20l4FBPiz49AE55ek9LTTdcn4bUxtssRjuM4m4HNNLFtprY6juOUwp+EHcdxRshmmtg2U1sd\nx3FK4U/CjuM4I8QnYcdxnBHiLmqO4zgjZDNNbJuprY7jOKXw5QjHcZwRspkmts3UVsdxnFL4k7Dj\nOM4I2UwT22Zqq+M4Tik205PwanZbdhzH2RRMlXylSKpJulPSPZIOSHpvfvwDkg5K2i/p85J2FulK\neqek70m6V9KnJU30a6tPwo7jjB1bS75SzGwRuNTMLgIuBF4t6WLgZuACM7sQeAB4T1pW0tnA/wG8\n2Mx+kWyl4Y392uqTsOM4Y0e15KsIM5vL39by08zMbjWzlfz4HcA5XYpXgG2SqsA08Fi/tvok7DjO\n2LG1Wu5VhKQtku4BDgO3mNldySlvBb6aljOzx4APAg8DjwJHzezWfm31SdhxnLGjWi33KsLMVvLl\niHOAl0h6UfhO0h8By2b2mbScpF3A64BzgbOB7ZLe1Let/U6QVAO+CUzk519vZvvyBevfpbUJ8dVm\n9rV+9TmO42w0WyvFx7/ZyF5lMLMnJX0d2AvcJ+lK4HLgFV2KXAb8MN9/E0l/BfwK0DFhx/SdhM1s\nUdKlZjYnqQLcLik8in/IzD5UqkeO4zhDottT7iuq7TPovz/W/r2k08medI9JmgJeBVwjaS/wLuBX\nc+NdEQ8DL5U0CSwCrwTSpYzOtvY7AYoXqkOby5R3HMcZJltray56FnCdpC1ky7WfNbMbJT1Athpw\niySAO8zs9yWdBVxrZleY2bclXQ/cAyznPz/aT7DUJJw36LvAc4CPmNldki4H3ibpXwDfAf61mR3r\nVY/jOM5QWGMYmpkdAF5ccPx5Xc4/BFwRfd4H7FuNZinDXLJQfXG+UP1nwM/nfnOHAV+WcBzn5GAQ\nH7Uhs6pm5AvVs8DeZC34WuDL3crdFL1/DvDc1Yg6jjO2zM7OMjs7u/4VnyQTbBlkZr1P6Fyovgm4\nBrjbzA7n57wT+GUz63DHkGQfLNGQQTLh/7jP91Wy6JgpBrs2ZXSmcq1BdJ7o832of5x0wvUfROd4\nn+9DXwbVmR+g7Gro9juxk1bY7aA5Eh7so7+T9nEbhMv6zDUAkjCzgWxNkrosHhSc+wAD6w1KmXux\n20L1JyVdCKyQXcvf27hmOo7jrIJN9CRcxkWt20L1v9yQFjmO4wzK2r0jhs4m+nvhOI5Tkk00s22i\npjqO45RkE81sm6ipjtOb5T7f+83+NKJL2PLJiN+XjuOMH5toZttETXUcxynJJprZNlFTHcdxSrKJ\nZrZN1FTHcZySuIuas5Es4xfOcXqyiX5BNlFTHcdxSuLeEY7jOCNkE81svsec4zjjxxpTWUqqSbpT\n0j2SDuTbuCHpA5IOStov6fOSdnaTzjcKvVvSDWWa6pOw4zjjR6XkKyHfuujSPH/6hcCrJV0M3Axc\nkOdPfwB4Tw/1dwD3lW2qT8KO44wfAyR1L9rOzcxuNbOV/PgdZBtcdCDpHLLNQD+2mqaeFNRH3YCT\niH5jcdJctE3MIPeb36ubgMm1Fy3azi055a3AX3Qp/qdkG4KeUlbPn4Qdxxk/1rgcAR3bub0k384N\nAEl/RLbJRcc29pJeAxwxs/1kmyCXShbvD1WO44wfXWa22UeyVxny7dy+DuwF7pN0JdlSwyu6FLkE\neG2+CfIUsEPSJ/vlXvdJ2HGc8aPLzDbz7OwV2Pf/tX9fsJ3bq4BrJO0lW2b41dx414GZXQ1cndfz\ncrId6PtufuGTsOM448fagzW6bef2ADAB3CIJ4A4z+31JZwHXmtkV3avsjU/Cm5h++XPXo/5wg2y0\nMWqZbDPJjdSpsz4bVp4srMdY9bqHBtl8d+SscWbrsZ1b4dahZnYI6JiAzewbwDfKaPY1zPVwXj5V\n0s2S7pd0k6TS1kDHcZwNZQAXtWHTdxLu4bx8FXCrmZ0P3EZv52XHcZzhUSv5Ogko5aJW5LwMvA64\nLj9+HfD6dW+d4zjOWhinJ2FoxkLfAxwGbsmdl/eY2REAMzsMnLFxzXQcx1kFm2gSLtWMPFzvojxp\nxRckXUD2NNx2WrfyN0XvnwM8d7WtdBxnLJmdnWV2dnb9Kx7XVJa58/IsmfPyEUl7zOyIpDOBH3cr\n98qCY7Fld1h/kDbam6DO8EJa6ww+bmXDo4PnQji23n0sSlK/Xjpbab/ucV82K8sMz3MhXJuN+t2Z\nmZlhZmam+Xnfvn3rU/FJ8pRbhjLeEacHz4fIefkgcANwZX7am4EvbVAbHcdxVseYLUd0c16+A/ic\npLcCDwFv2MB2Oo7jlGecliN6OC8/AVy2EY1yHMcZiAGyqA2bk+SB3HEcZx3ZRDPbUJq6u4f4VmAn\n7YaGHcDxVdT/oz7fb821pvO6A1PA/Cp07u7xXejDVK4X66TGoX58u+DYbuDJXGdHXmeq8wyydaGy\n3FtwbCfZmOymdX1CnwLnrpPOcv4zXJugt1aKdCAbo3B9wr0wiM7BLsenaL82aRtWc08DPJh8PoOW\n9btbmPczgEdXoVFk/Ay/F8EAuB5G4KEzTssRjuM4m45NNLNtoqY6juOUZBPNbJuoqY7jOCXx5QjH\ncZwRsom8I3yPOcdxxo817jHXI3XvByQdlLRf0ufzFA5F5fdK+r6kH0h6d5mmyqxryod1QZLdEX2O\nPRKCpXonsHsSdu6k9RfsmUDJvaDufLj398HCvxPYvQ2mtgHb8y93krkdlOD2H3b/boqsL0Fn5zaY\nOoU1/UW+s4dO6Etz3E6BrdvWRycMRbguwZtg92TWH2pr07k70YnvgdCXaWBH0NmWf7nK2OW7k/sg\n6MReHlPk1ybuxzZYPlZe52CXc2MvnCr5tYn/14zdG0pw70+764RrM1WFHdsSnVVwMNGIfx3axi1c\nm8A24KlVij3ef66RhJmV2iCzRx3WNun0OveldOhJmjazOUkV4Hbg7WRDc5uZrUi6BjAze09Sbgvw\nA7JMDY8BdwFvNLPv92qDPwk7jjN+DBC2XJS618xuzROZAdxBthNzysXAA2b2kJktA39BlvK3Jz4J\nO44zfgwwCXdJ3RvzVuCrBUWfQfv/7z/Kj/XEJ2HHccaPNa4JQ5a6N99J6BzgJZJeFL6T9EdkuzF/\nZr2a6t4RjuOMH11mttlvZ68y5Kl7v06Wuvc+SVcClwOv6FLkUeBZ0edzKBHA6JOw4zjjR5f942Ze\nlr0C+z7S/r2k08medI9FqXuvkbQXeBfwq/m+m0XcBTxX0rnAIeCNwG/1a+pQvCPi3A67t0G1mltz\ng8V9G3BK/oqt8C8EDvTXePRjncfOCJbpUFctqj/ohWX3khz5L+2fpyZhqpZ7J8RasUbozyp0Hk10\ndsSeI0EnHregF74v6VHwRDJu1Wqus53sX7V03ML7av6zpHX8yURnayX3HInHLbQ/9Cn8q7gdOFFO\nZ/kLiU5oZxj/+N4K4xl09gBHyunw5eRzJW9n6EfQOyU6thZuTD6H6xHGKtWNvUoWSmrcXHAsXId4\n3II7TmgDZP1bhVcJ1w7RO6KPx1Tz3Ge1e0dI+gWyPTPj1L3/TtIDwAQQ/EnuMLPfl3QWcK2ZXZGX\n3wt8OC/7cTO7pl8b/EnYcZzxY40zW4/Uvc/rcv4h4Iro89eA81ej6ZOw4zjjxyaa2TZRUx3Hccph\nnjvCcRxndDQ20cw2lKY+45W0jDrBaDRJu8FnG50GJshCPfvV//fRh2B4qdBpHAlGjNjAFAxD3eyd\nEXseKNCJDSWh3thAF4wbq9Dp6E8wjsRjFfclHrMKcBot80EPdsfGi9D+Wl53rFVkEILSYbg747Dl\nMO7dxiq8X0O479ZwziRwIfAw7YbSomu2FmKds/P3FYqvS7h2gdWE+8Y6Z9Bp8IuvWWwwg/Ih0vE5\n28gMlPF9EPejSOdZZON8kjFWk7Ckc4BPkl2eFeCjZvaf8sQWv0vrMl6dL0o7juOMlMXaRMkzlza0\nHWUo8/eiDvyhme2XtB34rqRb8u8+ZGYf2rjmOY7jrJ5GZfMsCpfZbfkwWQw1ZnZC0kFa8dAD+fM5\njuNsBI1NlNV9VbkjJJ1Http2Z37obXl+zY9JOmWd2+Y4jrMm6lRKvU4GSi9f50sR1wPvyJ+I/wz4\nt2Zmkv4E+BDwO0Vl3/dDmhFYM6fBTAljm+M448/s7Cyzs7PrXm9jEzl+lQpbllQFvgJ81cw+XPD9\nucCXzewXC74z+waZRXiB1rS/jcyKHyyxVWAS6rWWZbNeyR7UD9XOZpEac0zxosX7Otq37TsrrQ9P\nRjpVsr3bE08FyzUa1UyjUc0Ej1Z20aDCIhOc3TjUobPzjmgRP+ik1vfI6yL0JWg0KhUaVDjODhpU\nOK2LC8Op35xvfXiKbFU+9ooIetXMH3JpslOjToUT7KBOhdO76Jx+exQPvJjrBOv3aTTDYq2S9WOx\n1q6RFasxzxSL1NjTJe739DsjnVPITLkVWt4xk+3XZLE20aYBcJwdLDJBg2o5ndgLJw+F7qYRPxEd\nZRdL1Fhkgmd22VXgzNvzWN1wL4d7bRvY7s57a6ky0abRoMpPOa2nRodOeOXeL/VtnffWIhPN+gHm\nmeJxTu+p8Yzbo3sj9rDJvZfqtc7rXqdCg2rz+jzG2ZzNY101As8qERe+XmHLD1m5J71z9eOB9Qal\n7J+LPwfuiydgSWfm68UA/wT43no3znEcZy1spjXhMi5qlwC/DRzIEx0bcDXwJkkXkrmtPQj83ga2\n03EcpzThv4LNQBnviNspTn/sPsGO45yUbKY14c3TUsdxnJJspuWIoeQT/pb9I3ZxlGnmmsePs4Pj\n7GCJltEiXuxPaVChQoNFJjjODh7hmVzEfgB2cbR5XqzToMLP2MU8021Gi0ZTr1MrtGWJGg0qPM5p\nTcPGnijGc0fegwYV5phmjummgSrUHeovcoUJf6lD/6dpGeNiQ8oOjlNjkToV5pluGqiWIp3U6BPG\nKqZIJzam7OA4FRpNw2QwUMX1x08Xcf9irTmmAdqudTedcA9kmq1Y2EZbf1r61T468bhNMUeVRtNw\nmI1breOax9eo6OlpKb/fzoiMSs/mweb7iTziKu1PqlN0j4dxW2KCBhWmomsD8Fz+V5tOqPcou5hj\nupRGfM5RdrGLo23X63zub6sjvgfC/dbv/o11wu9CkfH0X/GJjmMp62WYu9teWOrcF+vgpjHMOY7j\nbBpOFh/gMvhGn47jjB3ZfwP9XymSapLulHSPpAN5jhwk/aak70lqSOpI+h6VP0XSX0o6KOlvJb2k\nX1v9SdhxnLFjrWvCZrYo6VIzm5NUAW6X9FWyjdZ+A/ivfar4MHCjmf2zPL5iup+mT8KO44wdSwO4\nqJlZMDSEMDIzs/sBJHVdP5a0E3iZmV2Z11MnC+vqiS9HOI4zdgySO0LSljwm4jBwi5ndVVL22cDj\nkj4h6W5JH813bO7JUJ6Eb+RydnGUHRwHaPMoKLIkdyNY0TMr9HZmuZRdHG2z8KY6x9nRVSe2Esff\nBU+Eudw6fJRd7OA4L6IVMh3rhPOWmCh0Eg9W/SIvhjj0M+hcmHt9hDaezuPUqbCU9z3TqUXn1GlQ\npZLH6oa1rqJxfYRnNq3kcX8qNJpeJnNMNUOel9p02j0uenmZPNbMdk6bTvCOAAq9CVKNuP6iX5pD\nkc4vRFtzB++I2Kukm6dFoMgLJPDXtPZJj3UmWGqO/zxT/IxdzTrC9QgEb420X934O57TphPKBO+I\nmFQr9p5Jicf4Qc4r/G6Jibbf0fi7+H4LHjTpOffxwjZPHIB/1bO360s3P+F7Z3/GvbNHC78LmNkK\ncFH+ZPtFSS8ys858CZ1UyTYJ/QMz+46k/wBcBby3XyHHcZyxotsfuAtmTueCmdObnz+976GudZjZ\nk5K+DuwFykzCPwIeMbPv5J+vB97dr5AvRziOM3aEp/N+rxRJp4e0vPlSwquA76enFWma2RHgEUnP\nzw+9khKTtz8JO44zdgzgJ3wWcJ2kLWQPqZ81sxslvR74T8DpwFck7TezV0s6C7jWzK7Iy78d+LSk\nrcAPgbf0E/RJ2HGcsWOpbTfS8pjZAbJ13fT4F4EvFhw/BFwRff4b4JdXo+mTsOM4Y8dmyh0xlEn4\nO/wSu/gZO8iSbmfeEVNNy2pKWLEJxFb44B3RoMJSnuPgODua56Y6rdwE7V4L1aalt1MnWNQXqTXj\n+n/KaW0W6VhnnqmmFTn1Wkj10vj+kAci+DUeZVfbv1KncrRpwZ5jqtmueNxiq3WrLy0PidCu0Jeg\nk/anllvg5yJvgmKd4msTvClCuToVqjQ6dMJYxJ4rcZ+ryfWP6w99i/uySI0K9baxz7wwMit+Om6h\ndbFWUQ6J1thNND0tgl5gmrlm+Xjc4jGrUO96/UPZkNsj9tqIPUzifCjxuLU02vuU5jCJc0IsUWvW\nF2uk7cw8cWpddeJrEbRCjos5ptvqGzabKWzZn4Qdxxk7PJWl4zjOCPHlCMdxnBHik7DjOM4I8Uk4\n4UHOY1cekgu0GaPSBfRWMG+nYSYYmGKj3tzcNGdPt5KGpzohoXu8RpQaZVKd2DDXoMLc0hRLCzXm\ndk536IQ2hYTuvYwygThxfBiLBhUWlyZo1KssTrcMPz/lONPMNY1eqUEzNpaloatxaPRiZDBbXJro\n6M8OTmOKuaahpSjMN9VKw6ODQaZpMCvQ2cWuyJCVGcxiY1gYryLjbGyIW2KCxkqFypZGU2dpZ2qY\naxnMwrgVXZugk4ZHN6gyv9IK/V9cqLG0MMHS7nKGudiQ1c8AOLeU6VSrjabO3O7pNp1Q7jg7mI+M\nnanBrMZS0xAd9zn8vgDUJhdZXKhxfHpHWz0TLDb7Hq5PYyXvz5ZGM0wbaBor0/uxSGdpYSLb+XxI\nxMbTk52+EXOSzpF0W54b84Ckt+fHT5V0s6T7Jd0Uokwcx3FGzVoj5kZBmbDlOvCHZnYB8I+BP5D0\nArLEFLea2fnAbcB7Nq6ZjuM45dlMk3CZ3ZYPk6V0w8xOSDoInAO8Dnh5ftp1wCzZxOw4jjNSxtZP\nWNJ5wIXAHcCePGEFZnZY0hnr3jrHcZw1MJZ+wpK2k6Vme0f+RJxu09x12+bH3/f/cIIFaixy6swv\nMD3zS2tsruM448Ts7Cyzs7PrXu/JstRQhlJb3ud7JX0F+KqZfTg/dhCYMbMjks4Evm7Wuc+0JDtl\n8RC7Jo7m3gTVZlLybmGxwUIbexUEa/ISE8zNTdOoV1icn2BlscYpZz7ePC/WCRbrYElv6myJV4WK\nPQqWViaYO5HpLC3UYGGCU85pbeO9Y+JE0wsjJMBeWpmgXs+9CaoNKltib4Iu3h65ZX9xfoKVRhUW\nJjjtvJa3x/Ytx5lmvhlyOs9004uiUq23aXVLiN5YqTQt1PXlSlMn7c80c01PiuNL22nUqzTqFSrV\nvB/VeodWqB+gXq8wf2Ka+nL2eWWxBvVKm86uiaNtydbDtYx1JiYX28YPoLFSaY5to17NvBQWamyp\n1At1dkycaHpwzK9MNa9lfblCdWuDSjW8sj4F6vUKjXrutRDusUa1pbNQ5ZTzDjfPn56Yb3kELE01\n+x/KVLc2mJhcaurEfWg03+f3WL3CltpiV50w3nNz08wdn8quI7ClUqc2tdTWn8WFWrPuoFNfrmR1\nA1Qb2ZhFvzvxmDdWKsydmG7dl7lO2p/FhVpTI4xtfbnCylPTUG1k/XlqGurCnk9f1mvL+6vt/yx1\n7r/X/z3yLe/L5hP+c+C+MAHn3ABcmb9/M/CldWyX4zjOmlnMd7np9zoZ6LscIekS4LeBA/m+SwZc\nDbwf+JyktwIPAW/YyIY6juOUZazWhM3sdui6wHLZ+jbHcRxncNa6JiypBnwTmCCbH683s32SfhN4\nH/BC4JfN7O6CsucAnwT2ACtkyd7/Yz/NzfPnwnEcpyRrnYTNbFHSpWY2J6kC3C7pq8AB4DeA/9qj\neIip2J87MnxX0s1mlm6P1MZQJuFjj59KY1eVpTwcd2llom1BP9Aylkw1jSWxoaCeGzDmT0xDvQIL\nVTgBx6qnNutIdWKDDEB1ayPRKjbMzB2fahp8WBAswLHJdp3GdFZnCGsuMmLEOs2yueEnGJeafakD\nC/Czn+xqtefUCvMTWQjo3Nx0e5m83bHBJGVpIQ9fzY2L1NXUSfszP5mFzi4u1DhxdEemEelQbTAx\nudjsU9aXxPCTG2GyA506nE7TQDV/YpqlcC3rgmpmJJ6fXGoam1pj1rqGwbBIXaxQ69qfSrXeupah\nXXVYqpJp5YajcE80+5Bfw/i6rFRrsADMw7HtLZ2l7TUq1QaNvD+cyNtDVmapaizl/QlGq2Yfwr0f\n7rE6rExWu+qEcVg6MQ0nWr+6K9SYn8wMYeEaLeWGueY1hKz/C3mhXCf+3QltbLuWebua/Zk0llKd\noFEXS1WLdKpt/Rkmg/gJm9lc/rZGNkeamd0PIKmrEa9LTMUz6Nyjrg1/EnYcZ+wYZE0431/uu8Bz\ngI+Y2V1rqOM8spiKO/ud65Ow4zhjR7fliEOzP+Dw7A96ljWzFeAiSTuBL0p6kZmV2fIe6Iyp6He+\nT8KO44wdS13cz06b+QecNvMPmp//Zt//6FqHmT0p6evAXkpsXQ/NmIrrgU+ZWSm33bJ+wo7jOJuG\nep6Gtd8rRdLpISOkpCngVXSu6fYK7iiKqeiJT8KO44wdWeRr/1cBZwFfl7SfbD33JjO7UdLrJT0C\nvBT4Su4xgaSzJH0lfx9iKl4h6R5Jd0va26+tw1mOOFojXhhp8zxoa03L6go0LfCBZnjniZZFnIWs\n/sCJ6NymJTlY3yGz3uY6sTU4WLihwDI8Dyx211mcn2izvlPNLNbdtNrCehdy63EdWM50Vo5ta+oc\nB5amljr7Xgeq2eVbmZxgCZpW+EDTAh/6n+ik/VmazD0pgvU91qlmP5cWJpp9ahJb+U+0rOlFOsc4\nlYnJxdYYx+ejTGehxkq1xvxiq1y7TtSXLjonyKz9K41qNmYn8vMCWwWTVVbqlczaH2sk3h0sA1sp\nvN/mg/dIvZKN2XFoRo9Xgk4t9yyot/chtD/cYw263tfz8e/K0Sptv1CVvL5qNbtGwasj1B3rLCTv\nqy2N4M3RbN8J2se5RvY7MVlt3QuhH817Re06oV8LDJUBXNQOAC8uOP5F4IsFxw8BV+Tve8VUdMXX\nhB3HGTs2UwIfn4Qdxxk7xjafsOM4zmZgrHJHOI7jbDa6uaidjAxnEj4OEBnnTtTaDQShFdVqd+MP\ntIUqs0y24N/hCp3rBINUm+EHMuNPprMyOdEyygTjSlNHLQNJMDQc7aITjGtpaOZUpFWNtFIDU2qU\n+VmrihW2MR/CdOM2BSpkx6tAvcZKNbr5UgNTPOYnoH1D2hpL2yutvh+lNW7h+lTa+9TSiX6GvhBp\nxndZtdYy6MSGn/g+qJEZwurV1rEiA1M9+S5cn0pWyUo1MuCeSHSqZPdQMGgV9SXcY3Uyw196H1RC\nG/Pyx4GnknELOqE/cf3QMhY+lfSlSCdwgvb7vpq3Lx63UE9sjAx9KBqzoDMV3TOpYS6MXejPpFqG\ny9gYuUjruleTOoaEL0c4juOMEF+OcBzHGSHuHeE4jjNCfBJ2HMcZIT4JO47jjJBFav1POkkYziR8\nLKjl1urjkG+o3N6S8McrtsCHJNHQbjlOLf6BSqQTrOKhbNMLgyTUk9aXQSOEKqeW+CKd49G5sU4j\nP28rWThnaoVP+5KWJy9br7XaFFvfQ7PDK4SNxsRW/lhjgdZ1aZ5bbR+31DtiMu1TdE5s5U+9GNIQ\n2ym1rO9F16eeawUrfFw/tI91rBMI9WzNyz1F5/iGtof+xARLf3xdYit/aHM1Kh/GNx23alQ+7kOD\n4vYX6cRtgsyjIe1vqhP3oUhnks77OnhZBIJ3RGh3qD++PqEf9eT7XjpDYDM9CfdN4CPp45KOSLo3\nOvZeST/KE1SUSlLhOI4zLBpUSr1OBspkUfsE8OsFxz9kZi/OX19b53Y5juOsmbWmshwFZXZb/pak\ncwu+6pVT03EcZ2RsJj/hQfIJv03SfkkfC0mQHcdxTgY203LEWv9c/Bnwb83MJP0J8CHgd7qe/an3\nZSGVU8A/moEXzqxR1nGccWJ2dpbZ2dl1r3etE6ykGvBNYIJsfrzezPZJOhX4LHAu8CDwBjNLTdtI\neifZXLgE3ydCAAAZ1UlEQVQCHADeYmad26DHZcysTMPOBb5sZr+4mu/y743PG+wCttHK95DGo6ev\nyfz72GIOLY+IYHGdp2VBh3ad+LxAqlFJykMrSXjsubBMf53Ywk1ef9Aq8phJ2xg8GaYSnVrUnjgx\ne/BQiPsTe5kUWflj74gwxkFnMvouToIe+h3qDz8DvXI6FOkEC/oJ2nNNxF4Xk7TfA2nOgyKrf3x+\nKB/GLM7PEcZoKtFJ+0LUn2pUR3z+dto9AlLviFinaMzi9ocy3XQCP6PdiyG+JqE/cftjnZDfZIrO\n351wH8XnxvkmQt2xXnpdqhTrNMDeRF8kYWYDLXVKslMWD5U691jtrA49SdNmNiepAtwOvB34p8BP\nzewDkt4NnGpmVyXlzga+BbzAzJYkfRb4H2b2yV5tKPskLKI1YElnmtnh/OM/Ab5Xsh7HcZwNp1Ev\nO7V1YmZz+dsa2RxpwOuAl+fHrwNmgas6Cmd/drdJWgGmgcf66fVtqaTPADPAaZIeBt4LXCrpQrJH\n7geB3+tXj+M4zrBopFunrQJJW4DvAs8BPmJmd0naY2ZHAMzssKQz0nJm9pikDwIPA3PAzWZ2az+9\nMt4RRf9EfKJfOcdxnFHRbRJe+dZfs3L7t3qWNbMV4CJJO4EvSLqA7Gm47bS0nKRdZE/M55KFQl0v\n6U1m9pleepvHj8NxHKckYTPdDl4yk70CH7imax1m9qSkWWAvcCQ8DUs6E/hxQZHLgB+a2RMAkv4K\n+BXgJJiET9C+4P8z2sMwY4NMHIKZtjAs/gfjRzCgpUYK6DSWpVrBeBFrxhppiG9I7N1LJ21zqD/W\nSvsSG/6C4SjJZd9mMEsNgMH4szU6FpMamNIQ5rg/sbFzseD7ZnJyOq9T0IrDXLvpxIa5dNyCEbMo\nhDs1Ai3QHjIc15OGhBeFLccht3G5WCeMV0hWHn8P7Qbbbv0JdaUGwLiecE276cR6T9F+X4d2x+OW\nGhhD2fgeLvrdiY268bjFYxN0YgNgMMSFUPFUZ8hhyyuNtU1tkk4Hls3smKQp4FXANcANwJXA+4E3\nA18qKP4w8FJJk2S9fiVwVz9NfxJ2HGf8WPua8FnAdfm68Bbgs2Z2o6Q7gM9JeivwEPAGAElnAdea\n2RVm9m1J1wP3kP0pugf4aD9Bn4Qdxxk/FtY2tZnZAeDFBcefIFtuSI8fAq6IPu8D9q1G0ydhx3HG\njyEvfwyCT8KO44wfPgk7juOMEJ+EE4J3RLA2H6Pd8h2H+YYwz/T72MqfbhGfWt9jz4kia3Vs6U+9\nI4jqDXWkFudU56ke38c/izwXikKKU1LPiKL+TCbHUo00AXfqNQKZ1btfcvKQzHsqKZt6YcTtjgkW\n+DQMO9VJr1GRlT+12sdacThxmpw89o4J9RUlj4/HK/YeSfsde2Mcj74LOqFvsX6ceJ2o7d10Yr0Q\n7h3rpN4ksYdK2sZYs8hzhejc8LsWPHCCTrgPYo34fS+dYbDc/5STBX8Sdhxn/EjdPE9ifBJ2HGf8\n8OUIx3GcEVK03HaS4pOw4zjjhz8JJwRjxXayv1AhbDkEtcTGlyrthoDUyJQay4LxICbopN/Fhrig\nVaMznzC0G43CK9UJxp2QezfNi5sa5uLw4mA4SA1/Rca5heQVSEOwY52YdEfcbv3pZQCs0G6Uia9f\namBKw4tTnTAm6a7bsU647unaXmoAjA1zqQ7R8Th0NhjJgnExBJkWhWEXGeaWk3PiHLpPJf1ZpjNk\nuSiMt8gwl+oETtB+HxTppGOT6oT3Rb8foXy8U3Yw/gadNGy5W3/C+2FPij4JO47jjBCfhB3HcUaI\nu6g5juOMEHdRcxzHGSG+HOE4jjNC3EUtIViMgwX5GMVhw/O0hyx3C/WNvRZSi3SsE0JjA8HCG1vg\ni0KKi8KJQ3hyL520T3ES9HT36FgnfqUJsYMVP04YnobfFiVbj63cRVpF/Yk9I9Jw1TQJeq8w7Fiz\nSGcr7Ts6F+mkYdi9vBaKdEKdYXfv1KukTnFSdyjW6qUTWI1OkUb8XTed2IMl9obpppNORvH9U6fz\n2sTXNPZeSb2XQvm0D6F87HVR5LWy0WyiJ+Et/U6Q9HFJRyTdGx07VdLNku6XdJOkUza2mY7jOKug\n6AGnnysoIKkm6U5J90g6IOm9+fHSc56kLZLulnRDmab2nYTJNvX89eTYVcCtZnY+cBvwnjJijuM4\nQ2GNk7CZLQKXmtlFwIXAqyVdzOrmvHcA95Vtat9J2My+RRZeEfM64Lr8/XXA68sKOo7jbDjLJV8F\nmNlc/rZGtsBilJzzJJ0DXA58rGxTyzwJF3GGmR3JG3wYOGON9TiO46w/jZKvAvLlhHuAw8AtZnYX\nsKfknPenwLvIJu5SrJdhrrfgLe9rGVjOnoHpmex4bMSapN34U6XYyJQal4LhheicYCgryokbjH/B\ncNPNyJTuTlwHjvbRCQay1Pi3HL0PbUkNZt10tuc6wUgS64S6F2g30MUhsnF/0vyv6biF8Vmk3Zgy\nFbU51S8yMqU66bgFnadofxpJQ7BjHZI6QrmicYvbFIxYqQEwNmJ1y1tMUq5IZ572HYbjMQ3h3XGb\nu+VHjkPeu+lA+y7I4fw47L6o/TFBJxhe42sTh7wX6dTJng3j/qV9COVDGHi4n7owOzvL7Oxs9xPW\nSjfviAdn4aHeema2AlwkaSfwBUkX0DnHdcx5kl4DHDGz/ZJmAJVp6lon4SOS9pjZEUlnAj/uefY/\nfl92wSbJBucna1R1HGesmJmZYWZmpvl5375V7ZHZnW7eEefMZK/AN7vrmdmTkmaBvZSb8y4BXivp\ncrJHlx2SPmlm/7JXU8suR4j2Wf0G4Mr8/ZuBL5Wsx3EcZ+NZ45qwpNOD54OkKeBVwEFKzHlmdrWZ\nPcvMfh54I3BbvwkYSjwJS/oMMAOcJulh4L3ANcBfSnor8BDwhn71OI7jDI21hy2fBVwnaQvZQ+pn\nzexGSXcAn0vnPElnAdea2RVda+xD30nYzN7U5avL1irqOI6zoawxWMPMDgAvLjj+BAVznpkdAjom\nYDP7BvCNMpoetuw4zvixiSLmhjMJH6XdMHeUVmhq7B1RlJy8W2L32LodW3hDPcdpWd9TL4zYAyO2\nmMekGtDbm6DIOyL2wFiNTrxjb9Dp5h0RdhWOk3qnu/pCu5W+27gVeUcET4LQj1gnTR4fdFLPj27j\nFntHxN4esU7quRJ7D8R1pjrBch+8XFJvj1B/3Kai34Y0bDnVCdc0eEfMJ2VD/8K9nXqUpInXu+nE\n17Vo3KA93DuuK+1P/N2J5PvUO6Jo9+n0fk3rJS8b93XYk6KnsnQcxxkhPdziTjZ8EnYcZ/zw5QjH\ncZwR4ssRjuM4I8R31nAcxxkhvhyRcJxWfoPjwE9pxaGHvATbaeVB2EbL6luUByHNHxF7EwSdo8nx\nrXm94fvgVZDmqSjSoItOlVYi7zTR+mT0Oc2HkfYl/rlMu9dCMDB006lHP1Od2DiRWqjX0p/UOh7n\nwujWF1apE9cfe9CkpH0p0gk/U51FsusfnpaCJ07YUKBbn4InS6wT6izywkg9RWr0zlER6irSifNf\nFI1bqKOefJeS6sT3WuztUqQT+hHn7dgavSf5HL6Pc2MMC5+EHcdxRoivCTuO44wQd1FzHMcZIb4c\n4TiOM0J8OSLhKC2DxAmyfPXQMvgEg0ww1IWWpQaMeCdaKE4aHnSO0grJjENLp6LyceL1XiHFdNEh\n15iPzov7FOt10ykKKy4yzB2lcwfkYGAsCvWNjSzdwpf79We1OoHYCNNNp8hgFsZtG50hxdBpOIpZ\nrQ7RsXhDgUCFYg0Snbg93QxZ4ft+4dHxeWV14vs5NZrGdaWGxvBdamQs0klDo0MdU7SSt8f1duvP\nMHEXNcdxnBHiyxGO4zgjZBNNwmvd6NNxHOfkZe07a9Qk3SnpHkkHJL03P36qpJsl3S/pprD7RkH5\nvZK+L+kHkt5dpqk+CTuOM34slnwlmNkicKmZXQRcCLxa0sXAVcCtZnY+cBvwnrRsvhvHfwZ+HbgA\n+C1JL+jXVJ+EHccZP+olXwWY2Vz+tka2ZGvA64Dr8uPXAa8vKHox8ICZPWRmy8Bf5OV6MtCasKQH\ngWPACrBsZhcXnhi8I8L7w7TCRGvADrIBCaGkcdLw1LoL7ZbzUGcg6PyELPk1UV0h8XYtOjdOhF20\nJX3QSUM8g8fBU7S8MOLk9NvzY3EIbpqkPtaJ9WKd0Ief0G75jr1K4nEr6k8avkyf/sTbhcfeCrUu\n/Ym10r500ynaHj72ikhDilOdorDfIm+CkAA9bu+2pI54TONQ3CK66cQJ96EzXDf2KOml0+36QHsy\n9aAT38/xcej0kAn1bqX4d6ebThxqXXSNyugM22VsAL38ifa7wHOAj5jZXWGnZQAzOyzpjIKizwAe\niT7/iGxi7smghrkVYMbMfjZgPY7jOOtHNxe15Vmoz/YsamYrwEWSdgJfkHQB2dNw22mDNjEw6CQs\nfEnDcZyTjW7eEZqBrTOtzwv7ulZhZk9KmgX2AkfC07CkM4EfFxR5FHhW9Pmc/FhPBp1ADbhF0l2S\nfnfAuhzHcdaHNa4JSzo9eD5ImgJeBRwEbgCuzE97M/ClAtW7gOdKOlfSBPDGvFxPBn0SvsTMDkn6\nObLJ+KCZfWvAOh3HcQZj7WvCZwHX5evCW4DPmtmNku4APifprcBDwBsAJJ0FXGtmV5hZQ9LbgJvz\nsh83s4P9BAeahM3sUP7zJ5K+QLYI3TkJ/9378gJAfQZOzGSfp2jlqYX2hf+wsB8v/qe5S8m/LzIu\nHCEz/MTGmAaZESsYzdKQ4qCTEtqQ6gTjUgj9DDrB0BjaGvcnGM1Wo9MAHo/qrNIyaELnuK21P2Es\ng2Eu6MQ7Qge9EJIdwlaLtHrphN2jg04wVsWGrFBnkU7Rml83nfQeCjqx8TQ2AsdaRRTpxPl+i0Lg\nY6PpWnSKjHKhjtQIDJ27lDe6vC9j/Ivv27h8arRLQ79LhA7Pzs4yOzvb/8TVssZgDTM7ALy44PgT\nwGUFxw8BV0SfvwacvxrNNU/CkqaBLWZ2QtI24NeA4gWW096X/ew2YTqO87RkZmaGmZmZ5ud9+7qv\n0Y4rgzwJ7yGzHFpez6fN7Ob1aZbjOM7TgzVPwmb292QRJY7jOM4a8QQ+juOMIZsnobBPwo7jjCGb\nJ43a8JK6BwvyPGSuC1MwvzP7HMYrtfJWaQ/v7BbuW2RFPg7wJNSnoL41006ttWVCcANFOsdDO/JG\nxjvyBtJQ7KL+9NNpjltUOE68nibQXmt/mjoR8efV6hSFxnbTCZ4qaVLw4CFTiT4H0hDmfjqidQ1i\nennixFq9dBZoxVBV6fTk6ea5UlYnDouO25kmoYk9aMro9PLASDcqCMR9Ch408bnd+jNU/EnYcRxn\nhKR/5U9efBJ2HGcM8Sdhx3GcEeJrwo7jOCPEn4Qdx3FGiD8Jt3OUaEx+RJbdbYos6G4nHJ8qTqRe\no2WBjbdXj7dwD/XH5QB4EHiCLMHCTrCdcHSq/Zw0r8OqdebJPD3m8/5MQX1PZ2j2cvSzSCflp9F7\nCwUfpZWwYSfUd2Y6aR6FXv0ppfNE0p+drf6k3iVr1iHXCRVMwfxUe76IcI1ij5KQmyDopA87cTvb\ndMgK2xTM524CsbdA8FyhdWpbX8OGAD11nqRZcX2q3QsieF4UJd1ftU7kYlHf2p4YP5wL2e9O7LWS\nXrswhql3RJtOznzuXRTKxJ5JNXrPd93uiQ3Hn4Qdx3FGiHtHOI7jjBBfjnAcxxkhm2c5wrcmchxn\nDFnb1hqSzpF0m6S/lXRA0tvz4/9Q0v+U9DeSviRpe0fh7LxTJP2lpIN5HS/p19LhPAnXITNgPZH/\nfBTYScuKtAeOb80W8UOLwoJ/2Ck5NgClRoa2sYw1ngR2k60PLYOd00rAXqOddDfkVKfjmj1BZmT8\ncd6HnWRGwK1guzOdCu2GiakuOmlfLP4r/kT+CltV7aAZ6227W7sxB53Qr267O3fVOU7rGi2TjdvO\n/NydrfDYuD+xoaynTtzBYMgMOjtpGmnrU5nOU0k93XTic9quz3LenydznSpNgya7M0NTMI4V7Xic\n6tTpoTOf64QbK+jsgBNbs3Kp8axIJzZAd+3PXK4XCuUG5xO0NmcP9cbGxjRsudvO2M3+hFc48bTM\nqPlUdGpRwvjUCNdtZ+yhsOYn4Trwh2a2P59ovyPpFuDa/Pi3JF0J/Bvg/yoo/2HgRjP7Z5KqwHQ/\nQX8SdhxnDFnbk7CZHTaz/fn7E8D3ybayf160ddutwD9Ny+a7M7/MzD6Rl6+b2ZPpeSk+CTuOM4Ys\nl3x1R9J5ZDnT7wD+VtJr86/eQLaTcsqzgcclfULS3ZI+mm8W2hOfhB3HGUPmS76KyZcirgfekT8R\n/w7wB5LuItvRcamgWJVsf7qPmNmLydaPrurXUveOcBxnDOn2lPv9/NWdfC33euBTZvYlADO7H/j1\n/PvnAa8pKPoj4BEz+07++Xrg3f1a6pOw4zhjSDdr4HPzV+BLRSf9OXCfmX04HJD0c/mu8luAPwb+\nS1rIzI5IekTS883sB8Argfv6tXSgSVjSXuA/kC1rfNzM3l985hHgAPAA2eAEr4U8NJZpqO/OLONT\ntI9fvxY2k8QHYh1o37d7vhVSmibD7rdyU6hzMCqce2AEa399KrMax9u6h2jgvjqPRgcORp+rwGlR\nf6ZgYSqzgi/Sbmnvxzy0QnqLdKBlst/ZMpaH8NsFOhN4kxRrbpue6oTPW/OT8tBypjKdWtQfuujE\nyfE7dEKY95H8y3CfBc3drTDchbyt4WeqEbwywn3ZoRO8cYLOjlZ5290ap1D/8iA64V/oKVoDvDXz\nXIjHq8hLJZDqtF2b2DviSZr3WKh4YaqVsD54MsX3dHzfLZJdx/CrMXTW5h0h6RLgt4EDku4h8+25\nGni+pD/IP/+Vmf23/PyzgGvNLGx7/3bg05K2Aj8E3tJPc5At77cA/5lstn8MuEvSl8ys97P+0HgA\neN6QNfcD5w1ZE1ieBWaGr3tsFk4bge5PZ2Hn00X3W0BfV9MN4LvAS0egu16szS/OzG6nOOPF14D/\nWHD+IeCK6PPfAL+8Gs1BDHMXAw+Y2UNmtgz8BfC6AepbZx4Ygeb+EWgC9dnR6D45It0nnk66t49A\nE+DuEemuF4N7RwyLQZYjngE8En3+EdnE7DiOM2I8d4TjOM4I2TxZ1GRm/c8qKii9FHifme3NP18F\nWGqck7Q2AcdxnpaYmQYpL+lB4NySpz9kZucNojcog0zCFeB+MsPcIeDbwG+Z2cGeBR3HcZwma16O\nMLOGpLcBN9NyUfMJ2HEcZxWs+UnYcRzHGZwNyx0haa+k70v6gaS+oXurrPvjko5Iujc6dqqkmyXd\nL+kmSadE371H0gN5js9fG0C3W67RDdOWVJN0p6R7cs33Dqu/eV1b8mQkNwxLV9KDed7WeyR9e4i6\nHblgN1pX0vPzft6d/zwm6e1D0H2npO9JulfSpyVNDGmM35Hfx0P5/dkUmNm6v8gm9/9Ftji+lcyB\n9gXrWP//Rpbd6N7o2PuBf5O/fzdwTf7+RcA9ZEsv5+Xt0hp1zwQuzN9vJ1sTf8FGawPT+c8KWUan\ni4fR37y+dwL/L3DDEMf5h8CpybFh6P434C35+ypwyrDGOfq9eQx45kbqAmfnYzyRf/4s8OYh3McX\nAPeSxdJVyJYynzPMMT4ZXxtTaRZq89Xo81XAu9dZ41zaJ+HvA3vy92cC3y/SBr4KvGSd2vBF4LJh\naZMliP4OWUTOhmuSpeu7hSwcL0zCw9D9e+C05NiG6pLFTf9dwfGh3VfArwF/vdG6ZJPwQ8Cp+QR3\nwzDuY+A3yUJ8w+c/Bt5FFss+1N/dk+m1UcsRRYEcz9ggrcAZZnYEssTMwBld2vLoerRF7blG92yk\ndr4kcA9wGLjFzO7aaM2cPyX7JYkNB8PQNeAWSXdJ+t+HpFuUC3Z6CLox/xz4TP5+w3TN7DHgg8DD\nefljZnbrRmrmfA94Wb78MA1cTvbUP8wxPukY53zCG2ZxVGeu0VRrXbXNbMXMLiJ7Mr1Y0gUbrSnp\nNcARy3YZ6OW3uRHjfIll+VgvJ8vh+rICnfXWTXPBPkX2JLbRugAoS/jyWuAvu+ism66kXWQpBs4l\neyreJum3N1ITwLK8Mu8n++/qRrKlhnSzsnXXPdnZqEn4UeBZ0edzaE8NthEckbQHQNKZZJu/hbY8\nc73aooJco8PStmyrlFlg7xA0LwFeK+mHwH8HXiHpU8Dhje6rZUlRMLOfkC35XMzG9zfNBft5skl5\nKNcWeDXwXTN7PP+8kbqXAT80syfMrAF8AfiVDdYEwMw+YWa/ZGYzwFEyu8qwxvikZKMm4buA50o6\nV9IE8Eaydaf1RLQ/od0AXJm/fzOtRKE3AG/Mrb/PJksm+u0BdDtyjW6ktqTTg7VY2VYpryJbQ9vQ\n/prZ1Wb2LDP7ebLrd5uZ/QvgyxupK2k6/08DSdvI1kkPsPH9PQI8Iun5+aFXAn+70boRv0X2xy6w\nkboPAy+VNClJtPLebnhfJf1c/vNZwG+QLb8Ma4xPTjZqsZnsae1+snRmV61z3Z8hsyIvkt1QbyEz\nMtyaa94M7IrOfw+ZZfUg8GsD6F5C9u/TfrJ/pe7O+7l7o7SBX8h19pNZlv8oP75hmgVteDktw9yG\n6pKtzYbxPRDunWH0F/iHZA8Q+4G/IvOOGIbuNPATYEd0bKPH+b15+XuB62gmWt7wvn6TbG34HmBm\n2PfyyfjyYA3HcZwRMs6GOcdxnJMen4Qdx3FGiE/CjuM4I8QnYcdxnBHik7DjOM4I8UnYcRxnhPgk\n7DiOM0J8EnYcxxkh/z+M/O0VazeAAQAAAABJRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(sal_ref[0,:,0,:])\n", "plt.colorbar()" ] }, { "cell_type": "code", "execution_count": 9, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mw_gsw_CT_from_pt('output_file','input0.txt','input1.txt');exit\n" ] } ], "source": [ "temp_cons = gsw_calls.generic_gsw_caller('mw_gsw_CT_from_pt.m', [sal_ref[:], temp[:],])" ] }, { "cell_type": "code", "execution_count": 10, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWEAAAEACAYAAABiV8coAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXt0nFd97/35eUYaS7IkW3Iko8ROHEIu+C0NWRBSUqgo\nlwZ6Xm5vFxC6ONx63qwDgSx6Tgv0dC3is2C95CzgtH0p7yFA00Ch0IZb6OKS5FCVcsuFOJAQOwlJ\nnNixI8VyJMuSPOMZ/94/9t7P7Hn0jGZ0HUv5fdZ6PM9t9ndfnvn50f799t6iqhiGYRitYUOrM2AY\nhvFMxoywYRhGCzEjbBiG0ULMCBuGYbQQM8KGYRgtxIywYRhGC2naCIvIBhG5W0Ru9sdbROQWEXlA\nRH4gIr0rl03DMIz1yULehK8B7o+OPwTcpqoXAD8EPrycGTMMw3gm0JQRFpGzgNcAn49Ovw640e/f\nCLx+ebNmGIax/mn2Tfh/An8GxMPrBlV1FEBVnwQGljlvhmEY656GRlhE/hAYVdV7AJnnVhv/bBiG\nsUDyTdxzOfBaEXkN0AF0i8iXgCdFZFBVR0VkGzCW9WURMeNsGEbTqOp8L3sN2Syik83f/piqnrMU\nvaUiC5nAR0R+D/gvqvpaEfkfwLiqXiciHwS2qOqHMr6jMNMg5WPA3oXkO8Vwxrlr/RbzK+DeJej8\ncYPr1wL/Dbh7iTp/0uD6LHAy0vkurst+obyvwfVjXmsEeCzj+q3AK5vQ+fMG148Co8DtdXTS/Cvw\nsozzf9Hge6O4d4WHgN80oZPmfwMvB/7zIr67CF7QB09cC2deW3u+w39OAMeXqNFf5/wT18I5Xnc/\ncGKJOoAebXyPiCzZCIuIfrTJe/+SuUZfRK6h+iP8nKr+TYbGMK6btg14SlVf5s9fAfwVrpfhC6p6\nXaM8NPMmXI+PA/8kIu/C/XLetIS0DMMwlo22RX5PRHYB7wZeAJSB74nIv6jqI9E9vcDfAq9S1SdE\nZKs/vwH4NO5/6kPAnSLybVXdN5/mgoywqv4b8G9+/yjwioV83zAMYzVYwtvlRcDtqloEEJEfAW8E\nPhHd81bg66r6BICqHvHnLwUeUtXH/He/iosim9cIr+MRc8PPEE2A57RI99wW6Z7TIt2dqy/ZPbz6\nmq3UXSY6mtwyuA94iR+M1onr59ueuud8oE9E/lVE7hSRt/nzZwIHovsO+nPzspTuiNOc4RZpnmyB\nbquM8LNbpNsCYwi05D+dnuHV12yl7jKx2O4IVd0nItfhHB7HgT1AJXVbHrgE+H2gC/iZiPxssXld\nx0bYMIxnKvUM2z4a9A0AqnoDcAOAiHyM2rdbcG+4R1T1BHDCd1n8NvAEsCO67yx/blF5NQzDWLPU\nexP+Lb8Fvp1xj4icoapPicgO4A3AZalbvg38vyKSAwrAi4BPAQ8A54nI2cBh4C3AlY3yakbYMIx1\nxxIN29dFpA/Xt/geVT0mIlcBqqrX+y6LH+BiXivA9ap6P4CIXA3cQjVErWHsrRlhwzDWHYvtEwZQ\n1ZdmnPts6vgT1EZMhPPfBy5YiJ4ZYcMw1h1LMcKrjRlhwzDWHXXCz05LzAgbhrHuWEuGbS3l1TAM\noymsO8IwDKOFrCXDtpbyahiG0RT2JmwYhtFC1pJhW0t5NQzDaAp7EzYMw2ghFqJmGIbRQuxN2DAM\no4WsJcO2lvJqGIbRFG3NWrbyimajKcwIG4ax7sivISPccHkjESmIyO0iskdE7hWRj/jzHxGRgyJy\nt9+uWPnsGoZhNKYt19x2OtDw/wtVLYrIy1R1xk9i/BMR+Z6//ClV/dTKZtEwDGNhNP0mfBrQVFZV\ndcbvFvx31B/LSmTKMAxjKbQVWp2D5mlqtWUR2SAie4AngVtV9U5/6WoRuUdEPi8ivSuWS8MwjIWQ\nb3I7DWjKCKvqKVV9Pm7huktF5LnAZ4BzVfVinHG2bgnDME4PlmCEReQa7/+6V0TeX09CRF4oIidF\n5I3Ruf0i8kvvQ7uj2aw2jV9raQS4ItUX/DngO/W/+dFo/6V+Mwzjmc7IyAgjIyPLn/Ai33JFZBfw\nbuAFuNiJ74nIv6jqI6n7NgAfB36QSuIUMKyqTy9bVkVkK3BSVSdFpAN4JfBxEdmmqk/6294I3Fc/\nlb9sNj+GYTyDGB4eZnh4ODnevXv38iS8+MiHi4DbVbUI4JezfyNz15N7H3AT8MLUeaHJHoZAM/9f\nPAu40Vv+DcDXVPW7IvJFEbkYZ/n3A1ctRNgwDGPFWHx/733AR0VkC1AEXgPcGd8gIkPA61X1ZSJy\naer7CtwqImEV5s8tOauqei9wScb5/9jou4ZhGC1hkdERfjn764BbgePAHtyy9jF/BXwwOo6jxC5X\n1cMicgbOGO9V1R/Pp3ma+AcNwzCWkTqWbeQYjEzN/1VVvQG4AUBEPgYcSN3yAuCrIiLAVuDVInJS\nVW9W1cM+jadE5JvApYAZYcMwnmHUsWzDfW4L7D489x4ROcMb0R3AG4DL4uuqem507w3Ad1T1ZhHp\nBDao6nER6QJeBTTs5DYjbBjG+mNpQ5K/LiJ9wEngPT4q7CpAVfX61L0a7Q8C3xQRxdnWL6vqLY3E\nzAgbhrH+WIJlU9U5MbSq+tk6974r2n8UuHihemaEDcNYf6why7aGsmoYhtEka8iyraGsGoZhNMka\nmsDHjLBhGOuPNWTZ1lBWDcMwmuQ0mbC9GcwIG4ax/lhDlm0NZdUwDKNJ1pBlW0NZNQzDaBLrjjAM\nw2gha8iyraGsGoZhNMnGVmegecwIG4ax/rDuCMMwjBayhizbGsqqYRhGk6why7aGsmoYhtEk1h1h\nGIbRQtaQZWu4KqiIFETkdhHZIyL3ishH/PktInKLiDwgIj8Qkd6Vz65hGEYT5JvcTgMaGmG/9PPL\nVPX5uAmLX+1XGP0QcJuqXgD8EPjwiubUMAyjWQpNbqcBDY0wgKrO+N0C7v8PBV4H3OjP3wi8ftlz\nZxiGsRiW8CYsItf4v/rvFZH3Z1x/rYj80vcO3CEil0fXrhCRfSLyoIh8MP3dLJoywiKyQUT2AE8C\nt6rqncCgqo4CqOqTwEAzaRmGYaw4izTCIrILeDduReWLgf8gIuembrtNVX/b9w68G/i8/+4G4NPA\nHwC7gCtF5MJmstoQVT0FPF9EenAL2e2idoE7Mo4jPhrtv9RvhmE80xkZGWFkZGT5E158dMRFwO2+\nGxYR+RHwRuAT4YaoZwBgE3DK718KPKSqj/nvfhXXY7BvPsEFdU37VUdHgCuAUREZVNVREdkGjNX/\n5l+mjjuA2YVIL448UF55mWcWfcDRVdDpBqaWnoy0gZ5cejoNWaVnep0xPDzM8PBwcrx7d8MV4ptj\n8U63+4CPisgWoAi8BrgzfZOIvB74f4AzgD/0p88EDkS3HcQZ5qVlVUS2AidVdVJEOoBXAh8Hbgbe\nAVwHvB34dqO0DMMwVoU6lm1kH4w8UP9rqrpPRK4DbgWOA3uASsZ93wK+JSK/i/tT/5XLnNUangXc\n6Ps7NgBfU9XvisjPgX8SkXcBjwFvWmwmDMMwlpU63RHDu9wW2P2dufeo6g3ADQAi8jFq327T9/5Y\nRM4VkT7gCWBHdPksf25eGhphVb0XuCTj/FHgFY2+bxiGseosYRY1ETlDVZ8SkR3AG4DLUtefraoP\n+/1LgHZVPSoidwLnicjZwGHgLcCVjfROk3BlwzCMZWRplu3r/s32JPAe7wu7ClBVvR74v0TkPwIl\nnCPgTbiLFRG5GrgF12vwBVXd20hMVOcJalgGRES5so5GqKhZYDo6vxU4sgCRbQ2ul/32FDARnV+o\n76e7iTycwP3xEpdnoz/fLF0Z5y4C9gODuB6qE8BDKZ0zaeKPn9T9aZ7jdc70aZ8E7lmizq6Mc8/B\nuS224tpgAniApfniXlDnfAew2WtM4N5RFvJ8LURnk9dIt/dCnwGY+yf1xbi26PDHE7hey5h+YHwB\nGv0Z58Lvos0f72fhec9Am/DnigiqKkvRERHVbzR57xtZst5SsTdhwzDWH2vIsq2hrBqGYTTJGrJs\nayirhmEYTWJTWRqGYbQQW2POMAyjhdibcIr/Gu1vLcORlOxGYKPCfnFTBIGLNniS5mg0RUY8arUD\n+A0usgBgC/B0kzpvmOfacZwHuew/N+I82c2WIeYddc6XcV7rEKnwlP98INJ54QJ0rq49bN92jNKT\nPU4nePjLuEiCE7gR8IvR+VDt4YbBaU6NdlWjYk7gyjWJq8ef+RsXOhI4PZnqZoUJcekfp9pGR6iN\nKtjE3CiD+Xh5nfNFn05oo/3U5n+hUT/1nrdQlglcezzJ4kdMvy11HP/ugsYBXL73R/ed7qO019Dr\n5RrKqmEYRpOsIcu2hrJqGIbRJGvIsq2hrBqGYTSJ9QkbhmG0kDVk2dZQVg3DMJrkNFk/rhlWxQif\nf8mvkv3ncj8DO0Y5TjdlclTIk6NMgRKj5w4yxgAPHLuA2SOb+Z3fH+FnD76sYfrPPv/Xc869iNvZ\nzAQTbGaGTir+75McFY48r5/97OTgoe3wZPOtdfH5P6853sQUQxxKdGbppEiBEu3kqHDopUMcKj2L\nyd9sW1CUxNnn107EP8gYF/AAOSpM0c0U3VTIMcFmAEYZZOzoAKW7ehbk5b/ouXtqjgsUOb/vAQqU\nOEI/s3RSJsc4W6mQY7Q0wOS+bS66ZCNNe/ovet5cnXMG9yflmaGDWTr9ficHHz+nGkGTFU1TT+eS\njPLwACUKSb2VaOdpNjN2bJDZg1uSuSo2nDPNqf1Zk3Zk6LyoVqfbPwdF2plgCyXamaGTsVMDjB8c\nhN8s7md20YurOjnKdDJLP0dq6mqCzYwf7af0UE91Po8FzFNx0Uv3zDnXff4Um5nwdbaJcbYyNdPN\n8YNb3Xwf+/2NZ+GOT0fW0OvlGsqqYRhGk6why7aGsmoYhtEka8iyraGsGoZhNIdadIRhGEbrqKwh\ny7Yqk7q/Sr/Ni/kpr+A2uplKHFlF2ilRoEKOIu0cYogxBjnEUOJ0aoYDbE/2X8K/81J+xLM4lDiY\nSt5ZNkMnRdoZY5D9nMMRtjLFJvJUKDbhTo3z9GJ+ysv4VwYYI0eZCbYwxabECVSknf3sZIwBjtDP\nFN3kqVBuIoCxMxoP+mJ+ynO5nyEOJXkI6Y+zlRk6OMwQo77eSrTTwQxT885A79gazf4ddAYYBUjS\nLlFglAGm6OYwQxxgO0fYSon2pB0bEfIO8Bq+mzgzQ3lmvKNp3NfTb3g2E2xJvtPNVFPlKVAEYDMT\nvI0vMkU3WxmnSDtTdDPBFmbp4AhbOcSzOMwQM3Q2TDdNu9fZwgQv4d/pZ5xupih5nbAdYohx+nmU\ncwCokKeTmaY1N/vx9FuY4OX8bzqZSb4fNMbp923/LMYYTL7bzVRTbfOsqG12sp+LuYdOZminxBTd\nzNLBKINMsDkpz2hKp5m2Abif5ze8Z7kmdT8x3fg+gI1drZ/UfUOjG0TkLBH5oYj8WkTuFZH3+fMf\nEZGDInK3365Y+ewahmE0plhob2rLQkSu8bbuXhF5f8b1C0TkpyJyQkT+NHVtv4j8UkT2iMgdzeS1\nmZf2MvCnqnqPiGwCfiEit/prn1LVTzUjZBiGsVpUcovrFBaRXcC7cYtYlYHvici/qOoj0W3jwPuA\n12ckcQoYVtVmpwVr/Casqk+q6j1+/ziwl+rqZC19jTcMw8iiQq6pLYOLgNtVtaiqFeBHwBvjG1T1\niKr+Amek0whN2NWYBd0sIufglhu83Z+6WkTuEZHPi0jvQtIyDMNYKcrkmtoyuA94iYhsEZFO4DUQ\nOZ0ao8CtInKniPynZr7QtA/Rd0XcBFyjqsdF5DPAf1dVFZGPAp/CvcbP4TfX/iNFDvAYh7h8OMcr\nh5tVNQxjPTMyMsLIyMiyp1upY9p+OnKSn42czLwGoKr7ROQ64Fbc+NM9uPXNm+VyVT0sImfgjPFe\nVf3xfF9oKjpCRPLAvwDfU9W/zrh+NvAdVX1exjW9T5/NrqcedsUZA44CPUAvbnn3fmAApgc2MFYY\nTLywIUoieNIr5Hg6w+P7Br6V7Cc6h3ATaw96rS5gCLQXRvt6OcB2JtjMGIM8yjmUKNDJTJJOlsf3\nzXwt2f+tow8it1NdXrzXb31O81hfOwdy2xmnnyP0c5ghDjFEhRztlACYSdYur+W9fCbZ33ngsJuA\nfgz3X2aPKwddMD20gfHC1hqv9RgDHGA7OSq0UyTvhwZn8T7+Jtm/8KnH4EFc2wAM4NqlF548o5dx\ntjLGAI/7ejvMEEDNH3dTbMrUeSd/n+zv+sXDbmfSl8fX2bGhdsZzrq72s3OOFz7n//KbL5Llj7gJ\ncFESzz56EDnky9PltvIQHOlxZTnAdkYZ4HAqEie0TY4Ks3Xa5218KdHZPnOQwiFfnoKrs2MDriwu\nksSVx0VnVHXCszZfpMSf8Pna8kzjnukuKPe6soTImFCe43QnvxEXE1ScN3ohftb6GWfrzDgFP9G+\ndsFUb/U5DuUZ89ESoTybmOJ4ExESX+ZdDe9ZruiIx3SgqXvPlrF59UTkY8ABVf1fGdc+AkzV84s1\nuh5otjvi74D7YwMsItui62/EvcYbhmG0nCX0CePfYhGRHbj1Tb4yj1RiwEWk0/cYICJdwKtowi42\n7I4QkcuBPwbuFZE9uD6PvwDeKiIX47yB+4GrGqVlGIaxGhTJDj9rkq+LSB9uYbT3qOoxEbkKUFW9\nXkQGgbuAbuCUiFwDPBc4A/imiCjOtn5ZVW9pJNbQCKvqT8ieIvn7TRfJMAxjFanXJ9wMqvrSjHOf\njfZHyXbWHccFLiyINTS4zzAMoznqdTWcjqyKEd71o4fdy/tdOEfJIZyDyTvkGAB2QNfQKXYOHGbn\n0GGK/TDaOZA458IwzU1MUaBEN1OJ42bXjx6uit2FW+U4OOZ6ffp9wLkgA7BtYJJtOyfRXhjv21Tj\n2ApDqQcYA5wDJTg2nveTB6s6vwB+DBzzx13VcjAAPQMldu18GAYfZrpnA4cLQxzxzq2n2ZwM1w69\nU6Xoz6ed3ztc1bkH5zALjrle3P/Bg66+uraPsWNgjOIgHOnsr3FohvlzoepwinUu/OFjVZ37cBHg\n45GOL8u2nZNsG5jkuUMPM9p3f+LQDMOxg1NmK0cydWra5+e+bSZx8956h2nPzhI9Q67tz9v+MOP0\n8zSbGWdrUpbgyOxghoLXif/s3PUTr1P0ZUk55vI7YdvQJNuGJjnv3IcZ7RxIhkpPsJlRBpMh9OCG\nDeep0MFsjVOtRudx3GrEkWOuZ0eJnmcdZuf5h3n6jA4O8HDNMOOQ1hTddDOVWWdzdB7yGt4xl++F\nbTsm2facSYpDeznQeVYyDN8Noa86AkveexyGGMehWTXP9DHf/l5HeqGnt8Suix6GoYd54oyH/e9k\nIClPcASGIfmdzFKgmOlAX03MCBuGYbSQZuZoOV0wI2wYxrpjKX3Cq83ayalhGEaTWHeEYRhGC0n3\nr5/OmBE2DGPdYX3CaX5CErUwOwoPTbvgiHNCdMQQzvs/VD0u9MOOvjF2DI5xrK+didzmmlV52ykx\nxCFX2T+JtLzOsUNw9AT0bYSeWCPalz7Y2n+crUMPUux10QWzdCSTZs/Q4Yd+ltyfNz9K6dwFs+Nw\nsgI9ITpip08/KlNXzynOGzzIs/sOMt63KZnIPF4FOgxnrpCDf4t0fDTByTFoy+M8/Tujcpzr62sA\nzuwb58z+cbY/60DivQ4TswdqdH6S0jlAbXTE417H15sMwbbtk2zrmmT70AEmCm4iexeFUZjzJ2Au\nDLlP62RERyTtPwRbDs2yZeAg5f6DTPV01ETHpEc6hcgCoNo+J3DRBCE6ouB1Io3CGOwYGmP7wBhT\nve1M5VzkQpH2mpXAA+fxm2p5gk4ZeMTX2zRJdASP4yJLjsKWgVm2nPsg0z0bmCp0++G+W6iQS6Ju\n6v7pHOptGniUmugIen35xqEwBOftPMj2wYNMdW5K6muCzZTJJe2fqRO3zaSvozA8OgzFH3flOXPn\nOGcOjvP0GQeSZzgMX674CXHyvo7OYX9LuwSsT9gwDKOFWJ+wYRhGCzEjbBiG0UKsT9gwDKOFlJpY\nuPd0wYywYRjrDuuOSLPXbU88AqM4x3UPcGwczpmGnmmcp3k62o7i5nuYdOPXe3rH0C6Y6OtgnP7E\ne90Obm6FwENO52BI4gQMPu50ajQqOK/vUXdc6HXRBdoFpY1u3org9U1I6ex/3H39JNBzAgYnoS/W\nCPs97lPGYevAcbZuOk6xi8STHY+zz1F2FRQ4AKO+3tqAjnE4J12OaVdP9LnybJmeZUvXLMXeMWY6\nQ7THpprJvbN0jj3kIkragO6NrkyEpcOP+/0TQBd0TZ+iq3eccu84Uz0dyZwbIYohRzmZTyBdbycP\nwdQ05PM+quSEvxbX2yTkj8KW3lk2980y1TvOTM5Feripygv1dcLcEWNw7KjX6MI97am2keT5Gqd7\nYIpSrp0ZOhONCjlm6Ewmla95DuK5I0J0xFGfdliLYdIdd/Weoqt3ku6BKTYXJqiQZ4pNVMhToj15\nnmt09kY66eiIPqornPnyFI5Boec4m/uPJ1EloV3cnBjVt8NEZy9VpsmOjij4PPh22XJ8ls29s8x0\nbWCisIVx+pPoCKiNTKgpzypi3RGGYRgtxELUDMMwWsha6o5Y0GrLhmEYa4ElLm90jYjc67f317nn\nb0TkIb/a/MXR+StEZJ+IPCgiH2wmr2aEDcNYdyzWCIvILtyq8S/ArZLxH0Tk3NQ9rwaerarPwS3r\n9r/8+Q3Ap4E/AHYBV4rIhY3yujrdEWPAuHMujeJ8F7N4R9MJaBuHjl7cMNZCRq6KbpMu2FKZpXLG\nBEXvoJmhEz//umPUaYzh5qgOboGOcRgMzpkw/DdczOOcQ2WnUShA/9ARKFDjoKnRGa+WpQxM4Xb6\njvpyBA182nmqzpppKHRBYdNxcme4k8ExUyFXXcEZODleLU8e6AB6xqAvlKMXt/hUeJ4qfn8TFE5A\nYdMs3YVZCj3FJM1EJ1WeJ064OmsDZk5ARzxUOrRJwWvmXZvki7ClOEu5MEuxsIHg+8lR9o6mlM4Y\njE26InaUYfYEDIY0QzlC21Rc3UkZeoolOrtKFAsbmC10JkOY2ynO1Sm6hjk6BuNeo2MaeoKTKS5D\nuL8IPfkS5UKJ9q4SpZxr97DictCpZOjMjsOxaejY6HywdFEdKh3qz2t0cYp8/yTl3AYKhWLyHLvy\nVOhgJltnDJh0Oj09vn42Of1k+DeuvvIV2FyeJddbSRyNwcGYnuS/pm3SjrngjBuLyoMrkxShq+sU\n+f5x6HRGr0Q7OSp+SH4HeSrVIfKrTL0VuZvgIuB2VS0CiMiPcAsZfyK653XAFwFU9XYR6fXrzu0E\nHlLVx/x3v+rv3TefYMM3YRE5S0R+KCK/jl/PRWSLiNwiIg+IyA9EpHfBxTUMw1gBltAdcR/wEm/f\nOoHXMHc9uTNxMTGBg/5cvfPz0kx3RBn4U1XdBfwO8F7/iv0h4DZVvQD4IfDhJtIyDMNYcRZrhFV1\nH3AdcCvwXWAP1b9h6yENrs9LM6stPwk86fePi8he4Czca/bv+dtuBEZwhtkwDKOl1IsTfnTkAI+O\nHMi8FlDVG4AbAETkY9S+3QI8Qe3b8Vn+XDtu/rz0+XlZUJ+wiJyD66z+OTDol35GVZ8UkYGFpGUY\nhrFS1IsT3jG8kx3DO5Pjf939szn3iMgZqvqUiOwA3gBclrrlZuC9wNdE5DJgQlVHReQIcJ6InA0c\nBt4CXNkor00bYRHZBNwEXOPfiDV1S/o44drfAMddrs7FOy8Mw3jGMzIywsjIyLKnu0Rn4NdFpA83\nGPY9qnpMRK4CVFWvV9XvishrROQ3OBfmO3EXKyJyNXALrqv3C6q6t55IoCkjLCJ5nAH+kqp+258e\nFZFB/z/ANmr9rDVcuw2OPQF34LziYaH1Y7iogmPTznvNJM4Du4lqlET4zON6pwvQ0TNDruC6aYq0\nu+95Zqedw/gY1dXoO/z+YPD6hiG+caRELtLogs7pU8wWZpKJvIsUqkN4gdnJqk7gmC9LTzz8etqn\nO+k/U89GZ9csM52zlGinQsXpROU5OunrCBe1ADAbD1sOdRanH2cKyJehozBDZ6HT94U5D3ZNvU1W\n2yM8FLNFaEuXJXwnjmLJOY1c2dVZ+AEU8T+GSOfYtNOZpRqc0jed0pmk2h45kmiTPE4DZigWCnTi\n2ifRCe1ThJOTcKzso1YATkTDsCepLlEfPP45dy5fgE5K5AplcoVKMhy3QClTJ0QsHAPKPgpjTln6\nosYoQCEP7YVTVPIlcjmnUKSdCvn6OsedztQJaMtBR96d43ikEe7Pg+SgM+/KQYFoSLH77GTG/Xai\nZzoZFj3tI0pyvo7qPQNAYSN0dro2z/naApJIiU5mahYvSDM8PMzw8HByvHv37sz7FspSjLCqvjTj\n3GdTx1fX+e73gQsWotdsnPDfAfer6l9H524G3uH33w58O/0lwzCMVlCkvantdKDhm7CIXA78MXCv\niOzBdTv8Bc6D+E8i8i7cy+2bVjKjhmEYzbKu5o5Q1Z8w54/ohFcsb3YMwzCWzlqaO2Lt/HdhGIbR\nJGaE0xShXHauxjLuczb6nAU3tDfeiqntRDWtQvEUUKRSyLlhpdURucwWq+kGrUSvDG1+eHJNutM4\n508474dltldKkCMZ6pnWKYe8e8o4h0bPfGUoVjXIQ/sJKHS6hDN1oroKxzNepyNLI+jEjk0gXzmV\nDL0NTpOkTn15glYYHj17AnrKdTSCTq6aX4GaUPjwb015TkRtHtqpDG1Z9RR/FpyWAPmNpyhQpER7\nMvy2Qr7qZCrPLU8bZLdFGFIePWP5IsApKvkK7bliUl95Kk4tlMenEXTaQlnitMN+vlZbgEKxRDm3\ngUqhukpxPZ2TRVd3M0BH0Q35nvMMBx2vmS9CJX+K9kqJXM6tupxL6cRtE9KaPeHqr2Mj1aHL8zwD\n7cUipUIhmdt57pCI8qobRZtP2DAMo4Wsqz5hwzCMtYZ1RxiGYbSQ0mkSftYMZoQNw1h3WJ+wYRhG\nC7E+4TTRj+9qAAAfvUlEQVRFN9Q2RCsE4uiFZBLvClWPb5jcO1zz+7ky5POnEk9vvKDryXL18GT0\n9TLec52Vvp88PPEuJzplcrloOOaJWp2QPim9mjSDZz9VhrBJhST9HDmnU65NM9ZpC/VWgY6QXsh7\nRl1RdmXKlYk88P4tIZqgL9RbaJ/w1WTnBHPyTpFqVImvt1ylkgzFzWqfEOkRJ3Oy4r39QSfkvxLp\nVKrHuTI1nnfCn55BJ4rGiSNyakTjtDcyp95yOdf+5KpaCVGkTvwcnEy3S/wcFOZq58oAp2rSr6cT\nonHKvq1Ia5QjnehZyJUh54SSdqnRiZ5pyj4KI5SrDG2hjk5EGsVov+wibyqVCvnc3KiImvKsItYn\nbBiG0ULMCBuGYbQQ6xM2DMNoIdYnbBiG0UIsRC1NykmSOH1i0s6kDCcWBXde0n390XHs7IsdTYnT\nJJ1++jPa8pVTlIicGBk68fGcssTHWTrBQeOdGpUgEA0lTZwxZJCuxHR5wkbVmeWyM9dhltZI9uOy\npNsjXZeF4MysDl9N0ojSDVpheHQ5/VDM1y654NCqMkenXOugTYqa1sgqQ859St61f3UYbjlTpxw5\nNMuhLI3KkCpHrlJJpsiaTyd2mp4seydzVnunnX/RCtjVOvM68W+p4nXK1bK0zanEufUWnH+5XM47\nZmt1kmHYq4h1RxiGYbQQ644wDMNoIUuJjhCRDwDvBk4B9wLvVNVSdL0H+Afcop454JOq+vf+2n7c\n2iOngJOqemkjPTPChmGsOxZrhEVkCHgfcKGqlkTka7gFO78Y3fZe4Neq+loR2Qo8ICL/oKplnPEd\nVtWnm9U0I2wYxrpjiXHCOaBLRE4BncCh1HUFuv1+NzDuDTC4WUqbXTYOMCNsGMY6pBi8kQtEVQ+J\nyCeBx3HTN9+iqrelbvs0cLOIHMItS/zmOAngVhGpANer6ucaaa6aET4ZbWkSb284iC9keZbLKQ95\ntJ8e5hvOZZIZpkFtREHkuU7rBK94W3S5RiuKTpjjiW9yNGdWnZ2kjhc+vUH9MmaUJ3y2RftJu1Qy\ntjidMCy3gU6slUlWWTLSgNQE6PE9ldqIhZO4KIxMjbTefOXIKM/JSkZ5suo+/Sws0EaE4dF1teoc\nS8VFebgfjRu6XjdyION5SM7HURdQU5YwjcBqR0DMR7034ZmRO5kZuavu90RkM/A64Gxc3+5NIvJW\nVf1KdNsfAHtU9fdF5Nk4o/s8VT0OXK6qh0XkDH9+r6r+eL68NnxtFpEviMioiPwqOvcRETkoInf7\n7YpG6RiGYawWWet7VMhRGL6MLddenWwZvAJ4RFWPqmoF+Abw4tQ97/TnUdWHgUeBC/3xYf/5FPBN\noKFjrpm+ixtwlj/Np1T1Er99v4l0DMMwVoWyj4RutGXwOHCZiGwUEQFeDuxN3fMYfpFjERkEzgce\nEZFOEdnkz3cBrwLua5TXZlZb/rGInJ1xSRp91zAMoxUsNk5YVe8QkZuAPbhembuB60XkKndZrwc+\nCvx91Dvw56p6VER2At8UEcXZ1i+r6i2NNJfSJ3y1iLwNuAv4L6o6uYS0DMMwlo2lREeo6m5gd+r0\nZ6Prh8noHVDVR4GLF6q3WCP8GeC/q6qKyEeBT+GCmzO59ikYBw4D24HNixQ1DGN9MTIywsjIyLKn\nu+6nsvSdzoHPAd+Z7/5rt8BDR+AO4BgwVk+4GY9+RBK5kBG1EGgjil7I8rrHerEHu4nohbSnOpQp\nmQw7TnshNIgmaIsPYo/7POnMmW8j/i5z660p0vNKNOHxrxvsMF99p+Y3kEoqamUhNDvHeCoCJ5/1\nxUoUpdIgreQzer6kgp9wv5xdljiyolH66UiMZp+79LNWqd3vIEVaZ4kMDw8zPDycHO/enX4BXRzF\n0vqbwEeI+oBFZJuqPukP30gTnc+GYRirRaW8doZANMypiHwFGAb6ReRx4CPAy0TkYtwQvf3AVSuY\nR8MwjAVRKa+j7ghVfWvG6RtWIC+GYRjLwroywoZhGGuN8kkzwgsSzde7ALXOjBSV3PwVHSfXFp9c\nhlLPN/S2rZ7GMi48mw/px1WQHoK7RNrwZUmz0PSbubdeU6YdqPO13WJ1YgdjnTZaM5OEN6iDectR\nbiqJeW/KcjK2ou5OVdbO++XayalhGEazWHeEYRhGCzmxdkzb2smpYRhGsyxDd9xqYUbYMIz1hxlh\nwzCMFmJGeH7CUOKa6IVlzkme6rLqIen8EjQq5Oc0bBswG+mFc4vXyGXO/hR0limwI9Gaj6WUI9ao\np5NPfS447Vy8EP38OulnrWmN1Jfm05kjuiCd2i800pmT/CIrsVHbtOUWnnYoS4UcuegHs+pzOcy7\ncsDphb0JG4ax/ljGcNCVxoywYRjrD+uOMAzDaCEnWp2B5jEjbBjG+mMNvQk3s8bcshIcPsFhEs/3\nm0meuUNN81WnyXxDIuP083hHQyNy1c9KvuoASvLit5DnJO1ovyZLTf4312ho5xy9tNMkrqcMD57m\nqjrNOLGSJLJuXSYPYfIsNJuWL6PWq6p58tVG/WvE13K1x/PVF1TzntTbAuolXY5yHcdsIP6ttIX2\nrjsHQC2VfD5J261RnZ/7zKW+m5Ql/Wxl3FvOVU1JWmexSw0tiazVx+utSJ5CRD4gIveJyK9E5Msi\n0p66/l9FZI9f5PheESn7VZoRkStEZJ+IPCgiH2wmq6tuhA3DMFacRRphERkC3gdcoqrPw/1385b4\nHlX9hKo+X1UvAT4MjKjqhIhsAD6NW/poF3CliFzYKKvWHWEYxvpjaSFqOaBLRE4BncChee69EvhH\nv38p8JCqPgYgIl8FXgfsm0/M3oQNw1h/VJrcUqjqIeCTwOPAE8CEqt6WJSEiHcAVwNf9qTOBA9Et\nB/25ebE3YcMw1h/1HHN7R2DfSN2v+b7d1wFnA5PATSLyVlX9Ssbt/yfwY1WdWEpWzQgbhrH+qBei\ntnPYbYFvzVlY9BXAI6p6FEBEvgG8GMgywm+h2hUB7s15R3R8lj83Ly3pjuig6ukN+4mnt8Y1n0Hw\nWuc2VIf5Znwn7eVvw3t766Ufe4Ajp3EYQjqfh7ymLHH6aY9ynUiDMDl9MzpJ2WLveL1ojFRZQ9qJ\nTo5ki9NOtnRZwhbIaLPgha8ZTjxP+7QRRXrkU2mSOhdRbZd83XqriYyJ6yOtlxF9M1fLlakmoiBX\n1Un0QnpZz1q9NqN2WG+WThjWnycVGZHWisuQdxEYWZEL6fvC6TavlY+vN4iMCJFKcRnS+6s+sfvi\noyMeBy4TkY0iIsDLgb3pm0SkF/g94NvR6TuB80TkbB9R8Rbg5kZZbWiEReQLIjIqIr+Kzm0RkVtE\n5AER+YHPkGEYxunBIo2wqt4B3ATsAX7pT18vIleJyP8d3fp64AeqOht9twJcDdwC/Br4qqrOMeBp\nmnkTvgEXchHzIeA2Vb0A+CEuTMMwDOP0YAlxwqq6W1UvUtXnqeo7VPWkqn5WVa+P7rkxaxFkVf2+\nql6gqs9R1Y83k9WGRlhVfww8nTr9OuBGv38j7n8FwzCM04OTTW6nAYt1zA2o6iiAqj4pIgPLmCfD\nMIyl8QycRU3nu3jt0zCBC5q7AOcy7PBb4hPYyFwHTJbjIZcxBHceh1wnKadZgaozpp5Dw38Gh0Zl\nHp0OUmVJO3+ytDKcXDUOppROcPzEWjVOsywHUIbzJAyLTRwldXRCueY4MmMHYx1nVjzMu5yhkzjJ\n4vJsxLVLo3Ikzr+oXF6jRqdc/WrcNg3rK9qCQ6s6BDfUXb623nJVnbZQlo3MdSjWGXofhsbHbZOl\nk89Dh//zuSNdV/HzlHL6OY180i61ztn8HKduPg/5ivtsK2S0S8ZzHTSgdp7nXOScrTd0eWRkhJGR\nkcxrS+IZMIHPqIgMquqoiGwDxua7+dozYP8RuB3XDXNskaKGYawvhoeHGR4eTo53754TMrY41uEE\nPuK3wM3AO/z+26kN0zAMw2gt66lPWES+AgwD/SLyOPAR4OPAP4vIu4DHgDetZCYNwzAWxHrqE84K\nw/C8YpnzYhiGsTysoe6I5XLMGYZhnD6YEU5RcJ7jTmAG50nuAbqJvON5qhESIYKhTgRAJQ+lQiHx\nJsce3g6vA64dOiKdxNO7kWxv7xydPEUKVe97pBN7+UP0RQfec11Iacyjo7lq9EWiU6jVCY00J5og\n9l6HdOPjaKvkqYmMKFGYE7WQjijoKGToZG2xp7zGA593OlG9dVKNVgj1lqQddNKe/lQdhkiPIu0U\nKFWjCYLORhc90laO2iddZ1kRK6lnzA3BrvX0l2ivfQ7yUbRHXBfhGQj7dTTi4fch/UQneg7aclDO\n4yI/4t9KVluk6iqOjAhtE3RqLIAfslwOBiykmX6eU21Szm2gkstRpECBYlSeCiXa/fNWMy/6ynOa\n9Pc2g70JG4ax/ii2OgPNY0bYMIz1h3VHGIZhtBDrjjAMw2gh6ylEzTAMY81h3REpNjrPcb5cO5dD\niJLo6QK6qPXCZkUXxB5f73EtUnD3etoiHUhFFGyMdOpFFuSq2sGrXKLgvPyRTsdG6DjhyhKiL3rw\n4+1jT3w9r3LkIS9SoER7pk6IJghl6cZHLWyMtno60T3BCx90KqkojFCeuM7auhropMqohdol4hMP\nfFSeEH2R9/VV0y5x2rH3PxWdUc5t8GWoeuBLtFd1yq6OOvwcEj35VNTKRuZGLkTPmOZqowpKvqKK\noTxdXqfo2rsmmmTTPPWVj86HchQKFH17FGmnnWJVJ9U+4KMxwjO8ibnPcqSjhWq7h3YJc1R0MEMx\nFYER8taRw71Jxr/JeaJkQhSJK0OhZr+DmWoU02piRtgwDKOFWJ+wYRhGC1lCiJqIfAB4N3AKuBd4\np6qWouthWaNH/KlvqOpH/bX9uAVCTwEnVfXSRnpmhA3DWH8ssjtCRIaA9wEXqmpJRL6GWyvui6lb\nf6Sqr81I4hQwrKrphTDqYkbYMIz1x9K6I3JAl4icwrllDmXcIxnnwvkFLaC8Oka4UHWUhKHEfUA/\nGU659HDP4JCLnB3FgnPIFYNDJnZgRDp5r5Pp/MvaIsdJuUCSfjKMtCul4x1ZPUGnC+csibXiLe3g\nKECxUOs4mVOefHUy7x6/taXTj9NNOX+Cc6tUKCRDSIMTaE69ReVJ6iurfdIOGt8+lTyUomHeoY1q\n6m0j5E94R1Y+0glpx/sZTq1yIQxZrzp/5uicgI4u6Jj2zsyuqM7S9RbyHz0XpY3hGWtPnGRhCHba\nEcwm6JmGkxWnmZShTnuH9phbjvZEI1MncgDWOJhj52lKJwzvrzr+ComTMcsJTMXVXfKnfJYjO+NZ\nC0OWS7TTTqlmv+T3izUewFVgkSFqqnpIRD6JW3V5BrhFVW/LuPV3ROQe3JL2f6aq94ckgFtFpAJc\nr6qfa6TZkiXvDcMwVpRFLvQpIptxa2ieDQwBm0QkPZPkL4Adqnox8GngW9G1y1X1EuA1wHtF5Hcb\nZdW6IwzDWH/U6xOeHoGZkfm++QrgEVU9CiAi3wBeDHwl3KCqx6P974nIZ0SkT1WPquphf/4pEfkm\ncCnw4/kEzQgbhrH+qNcn3D7stsCROcspPQ5cJiIbcR0zLwfujG8IS7v5/UsBUdWjItIJbFDV4yLS\nBbwKaLhekxlhwzDWH4sMUVPVO0TkJmAPzpTfDVwvIle5y3o98Eci8p/99Vngzf7rg8A3RURxtvXL\nqnpLI00zwoZhrD+WMGJOVXcz9w32s9H1vwX+NuN7jwIXL1RvSUa46cDkLuc57pmuivYA3Ruho9dd\nJ3yGLQzJTHt+N0Ip154Mv53jfe+t1ekD+noz0s+KXkhFLQQv9Qydzrsb6bR1Qc9ktSx9oSw9GWVI\ne/yjLUQtzNKZqdOd1gll6UmVaSPZkQwbq5EewUs+Q6fzjPfWlqd/0mvG7VKvzlJtUi7ATFe1XYDq\nflSenh43aXg+REakddIe+ajuyr5dQtRCkXbyybDlSCfv6r6vWC0bvdRGrqQ3r6E++sY9Y1GkQlKe\n9lod/2x34NNPP88ZUTgalcNFRFSHkucoz9Up+PSKOOMSP2PzRP0EjTjCJ0yuXkxrBE5QjZhIP8tp\nna7aKKK4nsJ+aCeb1L0+S30TXnBgsmEYxorzDJpFbcGByYZhGCvOGprAZ6kGNAQm3yki/2k5MmQY\nhrFkFhkn3AqW+iZ8uaoeFpEzcMZ4r6rOGxNnGIax4jxT+oSbDUy+9nGgDBPA5cBLcU6ZjuBkqucE\nytjKBZjxTqzEwdQfiXU5Z1ybL9xgP7WOrPm0orlgZwudiSMjcWhEjix6oe+Q0+lLl6WXuQ6N2Dnn\nHU7FrmpZqnMK1+q09ULfpJ/vt+CO65Yj7WDsgnKXc5jN0pk4TEreWTKnPJN+vtoCrk7r6QSHoy9b\ncDSVcu2JIxOqTpk5OhXc8N3YiZVun9gR1AXa5TUKBWbpTLYclbk63nkVVltO0kw7tNKOXz9kuVqO\nDmboZJaOmiHSNTrHqb5VZbVNhhNtpqtajqqjNHZmZej04pxmpHTqOBvLXeEZLiQOwPCcVfza0UUK\n7scSCNYgSyf9LEcO7Fk6axyAoTxBZzY4nDMYGRlhZGQk89qSOE3ecpth0UZ4IYHJ1/4WUIT9R6vL\ngxuGYQwPDzM8PJwc797dcGzDumMpb8KLCkw2DMMwqizaCC82MNkwDMOoYiPmDMNYh6wdz5wZYcMw\n1iFrxzO3Oka4322Dh/yqsTmcpzX2iPcxf2SBP57pavde6w5KtDtPfF+tVqJToNbDm9ao4/0v9tZG\nYIStRqcPzux1qzsnafQBA8z1+seecn+t3AUznR3M0kGRdma9Jz5LZ3A6lcegEW9xXYaoBR8ZMZNz\nGsFrPeMjC2qiSvrcJOjJUOHeSKcvpbWRGo/5VK8b5uvy78oEJBEM6fYBqsOTQ/r91NZbRkTBbKEz\nib6Y8R75TJ1KpJOP0umPyhK3jx9uXOyCqc5NST3FESVlcslzV6MTTxSzyacfl2dTpFGoRi2EcoRo\ngrCfo1zVGfDpFoFpqsOWQ77DigUZz7H7nXQyRXcymb9re9c27WE/fgbCb7PC3GiPvlTb9LqyTBW6\nk3YvUSDnh5EXKZCj7NePLiS6q4e9CRuGYbSQ2VZnoGnMCBuGsQ6xN2HDMIwWYn3ChmEYLcTehA3D\nMFqIvQnXMgRs9973sHx98MLHnuvYW56KntBe5yGfynUn3vApup3XdSjS2h55+WPPeNrDG+tFW9FH\nLcwk0QpR1EKsswPajlE78XgPtV7+4IkvUBNhUI681zN0cpxuaqIwYp2jvhxhroYCc6Mj4nJ5z7n2\nuqiFqZzzjoe6KkXRBTU64/4zjirpT5UpI7JkumcDM7mOZFL68Alk62z33w+z1G+kGk0Q2runVqfY\n67zwYY6CED0wRTd5KnN1Krg5HYJOaKO43lJzY2gUGREiCkK6U3RToFiNwoh1wD3TUPtshfaOozE2\nwVRP/Fx1UCGfRMaEKIZEZyDSKePmdKhE6fWnytHn6rPYC1O57kQjzFExhavDss9wjQa46IsC1YiP\n+HcTf/ZXfyfhd+iifAoUKFZ/l55Qn6vL4t+EReQDwLtx86XfC7xTVUsZ970Q+CnwZlX9hj93BfBX\nuBkqv6Cq1zXSs7mADcNYh8w2udUiIkPA+4BLVPV5uP/G35Jx3wbg48APUuc+DfwBsAu4UkQubJRT\nM8KGYaxDljShcA7oEpE80AkcyrjnfcBNwFh07lLgIVV9TFVPAl8FXtcop2aEDcNYh5xscqtFVQ8B\nnwQeB54AJlT1tvge/7b8elX9/3CrCwXOBA5Exwf9uXkxx5xhGOuQem+5v/RbNiKyGff2ejZuEeOb\nROStqvqV6La/Aj64TBldJSO8E+dUiNVih0k8DDYehhk56cLQ2ODEmKKbKbo5TjfsiNKdjnTyVIdI\nBydDPJQ5DCn2TrrpnurQWOdw2JToTLEJzk3pxBp5ah1XhWq6sSOr2Bscfy792PkTzvGcSCc4YoJD\nMywhnUq3xsnknXLOYdaROMyC0yzosTPSOUHVkRXaJq6nMCQ3cphN92xgquDaIAxNDfUFVPfjepsE\njvn94DwNjrksR2mvc5gF52WRdo4nbVJ1zNXohEnJx3G/xeBsTA8p97pagIm+at7j4d0TbE6O55Tn\nBNVnGGqHfMfOud4whHxD9Dz5ZwoS51bsmKvRKfv2L3rN8DyHYfJ+PzivZwvVZ3eCzc4BF7VHJXbM\nxW1znOrk8ZWoPIPUDMUOv5OZyCkX6q2TmUQ37z2XU2xKVuBePeo55p7rt8CX0je8AnhEVY8CiMg3\ngBcDsRF+AfBVERFgK/BqESnj3pxja3SWPzcv9iZsGMY6ZNEhao8Dl4nIRtx/ey8H7oxvUNXkvy4R\nuQH4jqreLCI54DwRORs4jHPoXdlI0IywYRjrkMWFqKnqHSJyE7DHJ3I3cL2IXOUu6/Xpr0TfrYjI\n1cAtVEPU9jbSNCNsGMY6ZPET+KjqbuYu1fbZOve+K3X8feCCheiZETYMYx1iw5YNwzBayDNk2HLT\nQ/T+D9xqdEVgL/Ao1VjpMJF0asLoeILyMDS2OuS2OpT4EEMu/VirAtyFC7Eukx2J0VOrWewl8fjG\nE5/Peg/5GAMQj325EOfp34cL1y5TOxR7Y6ThJw3XrmpkxKwfQhxP5p7ovCDSeQEuEuMuqhEZIf10\nhITfr0ZG1EZDTPkohgk2c4DtcFGk8xyf/t5IJ44mKdTqFP3k5EGj6CdBD9EEOSqM088og7VRGM/y\nz8GjVL3w8TDlqDzlqL5CGcLw23G2MkU33UwxTj9H6K/qlHE+6hPAQ/5cnsxolXgIeTU6wg0nnqKb\nUQYoUaCbKUYZdO0ToldO4KITpqmG7Mf15PXKvU7DTXzfkURGxFEkB9hOhRxDHMrW6fL1Nkk1Uib6\nnYSImFKuPRmeHKJIQnTEEfo5zBA5KmzngGubOBJnEhdREqIwQlRJFOURImLCEPI42iNMTH+Efvaz\nk05m2M4Bxhhkgs2sLmvnTXjRgzUWO0RvtRjZ1wLNh1dfE2BkT2t0fzLSmreNe0eOtkT30ZHHV11z\n5M7G96yIbgt+P8vLkkbMrSpLGTG3qCF6q8XIAy3QfGT1NQFG7mmN7k9GKo1vWgHua5ER3t8KI3zX\nqks63TVvhBc3Yq4VLKU7ImuI3qVLy45hGMZycHq85TaDOeYMw1iHrJ015kRVG9+V9UWRy4BrVfUK\nf/whXDDzdan7FidgGMYzElWVxnfVR0T24+Z+aIbHVPWcpegtlaUY4RzwAG5Y32HgDuDKZkaIGIZh\nGI5Fd0csdoieYRiGUWXRb8KGYRjG0lmxSd1F5AoR2SciD4rIss296dP+goiMisivonNbROQWEXlA\nRH4gIr3RtQ+LyEMisldEXrUE3bNE5Ici8msRuVdE3r/S2iJSEJHbRWSP1/zIapXXp7VBRO4WkZtX\nS1dE9ovIL32Z71hF3V4R+Wefzq9F5EUrrSsi5/ty3u0/J0Xk/aug+wERuU9EfiUiXxaR9lWq42v8\nc7wqv581gaou+4Yz7r/BdY63AfcAFy5j+r+LG4P3q+jcdcCf+/0PAh/3+8/FzYiUB87x+ZJF6m4D\nLvb7m3B94heutDbQ6T9zwM9xoYArXl6f3geAfwBuXsV6fgTYkjq3Grp/j1vUEaozRK9KPUe/m0O4\n5VBXTBe3VOkjQLs//hrw9lV4jncBv8KNw8vhujKfvZp1fDpuK5MoXAZ8Lzr+EPDBZdY4m1ojvA8Y\n9PvbgH1Z2sD3gBctUx6+hZsEelW0cetd3QW8cDU0cZNS3woMUzXCq6H7KNCfOreiurjB0w9nnF+1\n5wp4FfDvK62LM8KPAVu8gbt5NZ5j4I+Az0XHfwn8GW6w/Kr+dk+nbaW6Ixa11tISGVDVUQBVfZLq\nYt7pvDyxHHkRkXNwb+M/xz1AK6btuwT2AE8Ct6rqnSut6fmfuB9J7DhYDV0FbhWRO0XkT1ZJdydw\nRERu8F0D14tI5yroxryZ6goOK6arc9dRm1S3jtpKl/U+4CW++6ETeA3urX816/i0Yz0v9LliHkcR\n2YRbafUaVT2eobWs2qp6SlWfj3szvVREdq20poj8ITCqqvdQu5jhnOwtp67nclW9BPcjfa+IvCRD\nZ7l188AlwN967Wncm9hK6wIgIm3Aa4F/rqOzbLpSu47aEG5l4T9eSU0AVd2H63q4Ffgurqsha+z7\nMypaYKWM8KLWWloioyIyCCAi26jOa/UE7n/bZcmLuGWwbwK+pKrfXk1tVT0GjABXrILm5cBrReQR\n4B+B3xeRLwFPrnRZVfWw/3wK1+VzKStf3oPAAVUNszV8HWeUV6VtgVcDv1DVI/54JXWTddRUtQJ8\nE7eO2oqXVVVvUNUXqOowMIHzq6xWHZ+WrJQRvhO/1pKItOPWWrp5mTWE2je0m4F3+P23A9+Ozr/F\ne393AufhBpYslr8D7lfVv14NbRHZGrzFItIBvBLXh7ai5VXVv1DVHerW03oL8ENVfRvwnZXUFZFO\n/5cGItKF6ye9l5Uv7yhwQETO96deDvx6pXUjrsT9ZxdYSd1kHTUREVxZ719hTQBE5Az/uQN4A677\nZbXq+PRkpTqbcW9rD+BmdP3QMqf9FZwXuYh7oN6JczLc5jVvATZH938Y51ndC7xqCbqX4/58ugf3\np9Tdvpx9K6UN/JbXuQfnWf5v/vyKaWbk4feoOuZWVBfXNxvq997w7KxGeYHfxr1A3AN8g+qazCut\n2wk8BXRH51a6nj/iv/8r4EZcFNNqlPVHuL7hPcDwaj/Lp+NmgzUMwzBayHp2zBmGYZz2mBE2DMNo\nIWaEDcMwWogZYcMwjBZiRtgwDKOFmBE2DMNoIWaEDcMwWogZYcMwjBby/wNvdlkjwPsIkQAAAABJ\nRU5ErkJggg==\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(temp_cons[0,:,0,:])\n", "plt.colorbar()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## What wrappers exist?" ] }, { "cell_type": "code", "execution_count": 11, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/get_tides.m\r\n", "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/mw_gsw_CT_from_pt.m\r\n", "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/mw_gsw_CT_from_t.m\r\n", "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/mw_gsw_p_from_z.m\r\n", "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/mw_gsw_SA_from_SP.m\r\n", "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/mw_gsw_SR_from_SP.m\r\n", "/data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/startup.m\r\n" ] } ], "source": [ "!ls /data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/*.m" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Checking on the inputs for a wrapper" ] }, { "cell_type": "code", "execution_count": 12, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "function [] = mw_gsw_SA_from_SP(filename, SPfile, pfile, longfile, latfile)\r\n", " startup()\r\n", " SP = dlmread(SPfile, ',');\r\n", " p = dlmread(pfile, ',');\r\n", " long = dlmread(longfile, ',');\r\n", " lat = dlmread(latfile, ',');\r\n", "y = gsw_SA_from_SP(SP, p, long, lat);\r\n", " dlmwrite(filename, y, ',')\r\n" ] } ], "source": [ "!cat /data/nsoontie/MEOPAR/tools/I_ForcingFiles/OBC/mw_gsw_SA_from_SP.m" ] }, { "cell_type": "code", "execution_count": 13, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "mw_gsw_SA_from_SP('output_file','input0.txt','input1.txt','input2.txt','input3.txt');exit\n" ] } ], "source": [ "sal_abs = gsw_calls.generic_gsw_caller('mw_gsw_SA_from_SP.m', [sal[:], p[:], long[:], lat[:]])" ] }, { "cell_type": "code", "execution_count": 14, "metadata": { "collapsed": false }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 14, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAWEAAAEACAYAAABiV8coAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAIABJREFUeJztvXuQJFd95/v5TdV0dfe8H9ZoJFkS5o3WXom1Zda6rFsg\nzCC0Aq+9uoDjLkK+DkcsXAh8w4uQfa80ex0bQATYrI13F4FZwcJaIBtJsBJ6hNRgaVdCoBk0oEHI\nxnrPDEjDPHr6WdW/+0fmqfrVqayq7Kquqq7S7xOR0ZVZmed7zsns09nn9ziiqjiO4ziDYd2gK+A4\njvNixgdhx3GcAeKDsOM4zgDxQdhxHGeA+CDsOI4zQHwQdhzHGSC5B2ERWSciD4vIren+NhG5U0Qe\nE5E7RGRL76rpOI4zmqzkTfgDwKNm/2rgblV9JXAP8OHVrJjjOM6LgVyDsIicBVwKfMYcfhtwQ/r5\nBuDtq1s1x3Gc0Sfvm/CfAn8I2PC6Xap6BEBVDwOnrXLdHMdxRp62g7CIvBU4oqr7AWlxqsc/O47j\nrJBijnMuAi4XkUuBCWCTiHwBOCwiu1T1iIicDvwk62IR8cHZcZzcqGqrl722bBXR4/lPf1JVz+1G\nr1tkJQl8ROTXgf9bVS8XkY8BL6jqR0XkQ8A2Vb064xr9eI6y1+euRSNHM47dC1ycfi6m5Yefq6kT\na16efu6lTpHkr2HQuY1kwr5XOs3aklf3RJvvV3pvmunm0Zloc04e3bkuylgJm4CvA5dFxzdTuy95\n3qJa8UST418H3pXqdNNnlrfkGGtEpOtBWET0T3Ke+8d0P+h3Szf38CPAl0XkKuBJ4IrVqZLjOE53\ndPMS1G9WNAir6jeBb6afjwKX9KJSjuM43dDtfwj9ZJjquiLOfZFoArzcdUdW9xUD0Byk7mqxWlMo\n/WBkB+GXvEg04cU1KL3YdH0Q7oyRnY4YBV50DW6B94UzqgzTsz1MdXUcx8mFvwk7juMMkGEa2Iap\nro7jOLnwN2HHcZwB4oPwGqXXN+ZF1ZmOs4ZxFzXHcZwBMkwD2zDV1XEcJxc+HeE4jjNAhmlgG6a6\nOo7j5MLfhB3HcQbIMA1sw1TXoWE18rw6zlpgmN4oLcNUbx8rHMcZOYbJRW0lS947juMMBetzbjEi\nUhKRB0Vkn4gcEJFr0+MfE5GDIrJfRP5GRDZn6YrIB0Xk+yLyiIh8UUTG2tXVB2HHcUaOYs4tRlUX\ngItV9QLgfOAtInIhcCdwnqqeDzwOfDi+VkTOAP4v4LWq+kupxDvy1NVxHGekWJ93ZCs3HlLV2fRj\niWSMVFW925zyAPBbTUosABtEZBmYBJ5rVwUfhIeYUTIADpMhZa2wGvc+b78P2/0pdjEIi8g64LvA\nS4FPqepD0SlXAX8dX6eqz4nIx4GngFngzmjwzq5ruxNEpAR8CxhLz79JVfemcyW/R22p+2tU9Rvt\nynMcx+k16wvZx79VSbZWqOoycEE673uziLxGVR8FEJE/ApZU9UvxdSKyFXgbcA5wHLhJRN6Vda6l\n7SCsqgsicrGqzopIAbhfRG5Pv/6Eqn6iXRmO4zj9pNmb8BuK8Aaz/x+ONy9DVU+IyL3AHuBREbkS\nuJT6IiyXAD9OF0FGRP4W+DWgu0E4rUzDHEm6L3mudxzH6SfrS51dJyI7Sd50j4vIBPAm4CMisgf4\nQ+BfpMa7LJ4CXici48AC8EYgnspoIJd3hIisE5F9wGHgLjNH8r7UZeMzIrIlT1mO4zg9p1P3CNgN\n3Csi+4EHgTtU9Tbgz4GNwF0i8rCI/CWAiOwWka8DqOq3gZuAfcD3SF5SP92uqqKq7c6pnZzMkXyV\nxA3jp8Dzqqoi8ifAblX93Yxr9OM5yu5m4v9om++LafkTdGfM+Emb74upRrcGs37ptOs3jE43nGjz\nfbg33dIvnblVKCMPm5oc30zSjtVoy7MtvptItcLvT7dckmOsERFUtav/sEVE9eyc5z5F13rdsqLf\n4XSOZBrYE80FXw98rdl1d5jPLwVethJRZ2DYX7x+eWGMireHk4/p6Wmmp6dXv+AhepDavglnzJHc\nAXwEeFhVD6fnfBD4FVV9V8b1/ia8RnXa9Vsovx864a2uXzrd4G/CndHXN+GX5zz38eF4E94N3JD6\nzq0DblTV20Tk8yJyPrAMPAH8fu+q6TiOswKG6E04j4vaAeC1Gcf/TU9q5DiO0y0dekcMgiH6e+E4\njpOTIRrZhqiqjuM4ORmikW2Iquo4jpOTJmHLaxEfhB3HGT2GaGQboqo6juPkZIhGtiGqquM4Tk6G\naGQboqo6juPkxF3UHGft4Q/7i4ghutlDVFXHcZycuHeE4zjOABmikW2Iquo4jpOTIRrZfMl7x3FG\nj0LOLUJESiLyoIjsE5ED6VqaiMjHRORguojF36S51TNJF8F4WERuzVNVH4Qdxxk9OlxZI1266GJV\nvQA4H3iLiFwI3Amcp6rnA48DH26h/gHg0bxV9UHYedGQsbq5M6qM59wyyFpTU1XvTldhBngAOCvr\nWhE5i2Qx0M/kraoPwo7jjB4dTkdAyzU1A1cBtzdeCcCfkiwImnvduCGavnYcx8lJk5Ft+ulka0X6\nxntBOu97s4i8RlUfBRCRPyJZaahhGXsReStwRFX3i8gUOVej90HYcZzRo8nINvWSZAvs/V/Ni0jX\n1LwX2AM8KiJXkkw1vKHJJRcBl4vIpSSrQ20Skc+3WwDDpyMcxxk9OveO2CkiW9LPE8CbgB+KyB6S\naYbLU+NdA6p6jaqeraq/ALwDuCfPCkT+JtwDlsjXsd0aipa6vD5P+aEdvTZqLZEsJunGs/yU6f4X\nuNUzNEHtGej1s7bqdN4xzdbUfBwYA+4SEYAHVPXfishu4HpVvaxnVRWREvCttAJF4CZV3Ssi24Ab\ngXNIFvq8QlWPd1oRx3GcVaPDQbjFmpqZ6zer6iGgYQBW1W8C38yj2XY6ooXf3NXA3ar6SuAeWvvN\nOY7j9I9Szm0NkGtOOMtvDngbcEN6/Abg7ateO8dxnE7oMFhjEOQahJv4ze1S1SMAqnoYOK131XQc\nx1kBQzQI56pG5Df3VRE5j0Zn5KbOyXeYzy8FXrbSWjqOM5JMT08zPT29+gWPairL1G9umsRv7oiI\n7FLVIyJyOvCTZte9MeOYtYL36w9Sry28ZWoW62Gw8rerY7gvwXMhHFvttmV5k6yWznrq77tty7Cy\nROK5MApMTU0xNTVV3d+7d+/qFLxG3nLz0HY6oonf3EHgVuDK9LR3A7f0qI6O4zgrY8SmI5r5zT0A\nfFlErgKeBK7oYT0dx3HyM0rTES385o4Cl/SiUo7jOF3RJEPaWmSNvJA7juOsIkM0svWlqttbiK8H\nNlNvaNgEnFxB+c+0+X59qjWZlh2YAOZWoPNwi+9CGyZSvW50vp1xbDtwItXZlGrEOmcCz65AJ6s9\nm0nqup3a/QntWm2dpfRnuDfFSGelhrlHmhzfRO3+hGehG+PcwSbHJ6i/N3EdVvJMQxKGajmNmvW7\nWZj3Su9NFsGYGQyAqxEe3XdGaTrCcRxn6BiikW2Iquo4jpOTIRrZhqiqjuM4OfHpCMdxnAHi3hGO\n4zgDxN+E6/lF89l6CgRL9WZg+zhs3kztL9jPA23Wggo8+FTr74OFfzOwfQNMbAA2pl9uB47m01n/\n48ZjoQNDW4LO5g0wsYWO/iJPZOhU60DNW2EC2L4F1m+o6WQmPc2ps5nEA8O2ZT3pvdlAkkOvA53N\nkU7wFAheClXPlaCzIT2xBGSuYZDN9ug5CM9a7OWxeQNM2PuyAZZWkAn7tCbn2rYUSe+N/Q2z7g05\neOSFxmOvpnZf1pP02USpXqfBqb8FByON8AxA1G/h3gQ2AKdWINRvOhzZWuRP/xjwL0meyH8A3qOq\nJzKu3wP8GUlg22dV9aPtNH15I8dxRo8Ow5Zb5E+/EzhPVc8HHicjf3oaVfwXwJuB84B3isir2lXV\nB2HHcUaPLnJHZOVPV9W702ySAA8AZ2VceiHwuKo+qapLwF+T5F1viQ/CjuOMHh0u9AlN86dbrgJu\nz7j0TOonUZ9Jj7XEDXOO44weTUa26W8nWyui/Ok3i8hrVPVRABH5I2BJVb/U46o6juMMMU3Wj5t6\nfbIF9n6qeRFp/vR7SfKnPyoiVwKXAm9ocsmzwNlm/yxyRJGLatMFMVYFEVGb22H7BigWU2tusLhv\nALakW7COF0lMwc0C9Q3P/ufGY6cFy3Sw6Jei8rdQm/HJyZFIZyJYpkN9g1YofwO15Agr0Inbs8l6\njoR+K0Q6G8z3OTka6RSLqc5Gan1j+y18LqY/c1rHT3ymfn99IfUcsf0WPFaCRvglWoEVfukrkU6o\nZ2iL9fAI/Rn+Jd0FHMmnw9ei/VD/8JzZ5yB81wm3Rfuh/KBVoNZn1qtkJRnxY42gE+oe2hLcmEId\nIOnblayvfn37sUZEUFVZQalZZai28Ziqnns2dXoispPkTfd4mj/9DuAjwDLwceBfqGqG3wqISAF4\njGQdi0MkaWDeqaotRzF/E3YcZ/TofGRrlj/9cRK3tbtEBOABVf23IrIbuF5VL1PVioi8j8STIrio\ntX2N9EHYcZzRo8ORrUX+9EzXeFU9BFxm9r8BvHIlmj4IO44zcqhHzDmO4wyOyhCNbH2p6plvJjEi\nbKFmMBin3qhkjRjWwHRajvIfNzuhjGAMs+XGRoxwTs7w2F2xTmwosUYla9ywRrscRqYz/zHSiTXi\ntljjSQHYAWSaDurZbsOJ47bE98PqFFKt7eQKw60LW477J9yHDTT2W2BHPp31z5m2nA88Felk3TP7\nxlRprwGmLuPAGelnayi1xsbQtsBKwn2tzmk07zdroAvkvDc8Zz5vIDFQhmduI/XPsDWkBs4m6ec1\nxkgNwiJyFvB5ktuzDHxaVf9cRK4Ffo/arb4mnQ9xHMcZKAulsZxnLva0HnnI8/eiDPyBqu4XkY3A\nd0XkrvS7T6jqJ3pXPcdxnJVTKQzPpHCe1ZYPk4TvoaozInKQWiheV/58juM4vaAyRLksV5Q7QkTO\nJZltezA99D4R2S8inxGRLatcN8dxnI4oU8i1rQVyT1+nUxE3AR9I34j/Evj3qqoi8ifAJ4Dfzbr2\nur8neWceh6kdMJXD2OY4zugzPT3N9PT0qpdbGSLHr1xhyyJSBL4O3K6qn8z4/hzga6r6Sxnfqd5H\nkinahlJuILF6B+tuat0vlxLLZrlQe0k/VDqDBUrMMsFrFh5tqN+G7yzXdk6lW/B82E69Z0ERtASL\naVLvcmEdlWJyw44VtlKhwAJjnFE51KCz+QEziX+KxKMiWI3tVqy1ZaGUtKNSLFIpJGXPMUmFAjua\nuDBs+9ZcbSdY0q31PVjIo/4KGkkbSswxQZkCO5vo7Lx/prazQHJ/bNh1ao3XQqKxUFrHYikxjYd/\n9xYocZKNVCiyq0nc784Hjc4WElNu8JZJ70059SBIdMaq7QicZBMLjOXXsV44aXvsfbca9o3oGFtZ\npMQCY/x8k1UFTn8witUNz9pm0C21vgrP1WJhjDKFuoHhBXa01GjQsZ48W5L+Cn0FVJ8tqzHHBM+z\ns6XGmQ+aZyNohN/HLbCwofHZCm0J9+c5zuCMOjeLbM7OERe+WmHLT2q+N71z5Cdd63VL3j8XfwU8\nagdgETk9nS8G+FfA91e7co7jOJ0wTHPCeVzULgJ+BziQ5thU4BrgXSJyPonb2hPA7/ewno7jOLlZ\nIK+L2uDJ4x1xP9npj90n2HGcNckwzQkPT00dx3FyMkzTEX3JJ3yf/jM2cZJNnKRCgQKVdG8Ti4xV\nDSN2sj/ZT84NnyFxPTnGNp5jNxewH4CtHKteE3TCNT9jK3NM1hktKhSa3qRQl0VKVCjwPDuqho1d\nJg7UtmeWyWpbFtLr4jqH9hVSi1GoS2j/JDVjnDWkbOIkJRYoU2Au1VlgrFq/uHyrG/qumY41pmzi\nJAUqVcOkNVDV7kd9+Vk6s0wCMMls9VgznfAMBENiTSe7/yyL6flWx/bbBLMUqVQNlEm/lRqer0Bs\nOLP9dpJNnGaMSi/hiernsTTiKm5P/Hw1e8aDRlLnubrjL+Pv63QqFKr1mU2Nu2UKFNP+t23IqsMx\ntrKVY3X365U8VneOfQZOsok5JjP7P7Qj1pllklkm2cnzDde8l880HItZLcPcw/rqXOe+Vg4OjWHO\ncRxnaFgrPsB58IU+HccZOZI39fZbjIiURORBEdknIgfSHDmIyG+LyPdFpCIiDfmGzfVbROQrInJQ\nRH4gIr/arq7+Juw4zsjR6Zywqi6IyMWqOpsuV3S/iNwOHAB+E/gvbYr4JHCbqv7rNL5isp2mD8KO\n44wci124qKlqMDSEJJ6qqo8BSLq2URbp6syvV9Ur03LKJGFqLfHpCMdxRo5uckeIyLo0JuIwcJeq\nPpRT9iXA8yLyORF5WEQ+nS4W2pK+vAnfxqVs5VidJTtYUVt5KlhsmOwxtjLLBNNczFaO1Vl4t3Ks\nwTsieBJk6QRrcCDcmOCJMMskL7CDrRzjNdRCpoNOsOzXvCMa/wJb63VSr0YL9nOcwTG2somTnJ96\nfYT67eT5qnfELJNVL4zGtpQbyrc/Y51f5ECdTvAymWWCGTZRplD1QrDnBbLKtzoB22/BOwJo6U2Q\npZP1S3PI6Nj2BO8I61UShygXjUZ8b2L+jto66aE9BSqMsUiBMhWKnGRjqlMrI9yToBfuW55n/ke8\nwrRnrnpd8vxPVvdDMLFlMafOP/BSU9dKtd8XGav7HQ3fxwSPDVsGwKO8mlKUq/e9LWuyujS7j49M\n/4xHpo9lfhdQ1WXggvTN9mYReY2qNuZLaKRIsj7de1X1OyLyZ8DVwLXtLnIcxxkpmv3xOW9qJ+dN\n7azuf3Hvk03LUNUTInIvsAfIMwg/Azytqt9J928CPtTuIp+OcBxn5LA+zK22GBHZGdLyplMJbwJ+\nGJ+WpamqR4CnRST8C/NGcgze/ibsOM7I0YWf8G7gBhFZR/KSeqOq3iYibwf+HNgJfF1E9qvqW0Rk\nN3C9qoZl798PfFFE1gM/Bt7TTtAHYcdxRo7YlpEXVT1AMq8bH78ZuDnj+CHgMrP/PeBXVqLpg7Dj\nOCPHMOWO6Msg/B1+ma38jE0kSbeTfAsTdfkPLMHCavNG2GTSwRMhWIpPsql6bawT56cIFCONcH7Q\nCNbhoG0t0rHOHBPMRvkpkrLrLePQmEMg5IEIdTzG1rq6buMYT3AukHgtJHkw6vut1l81vdgDYyHV\nCG08xtY6D4ut/KxqzZ413gTZOvV9Zr0XQlvCfpFKQ79Zr4XQxwuUqvUPWQks1rMh3B/rMVOgXNee\nxAsj8VqI+816EwQt24bGvqv3eLE6k8xWr4/7rTbzWNOK73/t2kZPpjM4VKcT6mRzR4T7EetYLwzr\nxRLu0RgLABxhV1WjQLn6nIZFFObMvbPeE0HL/g5XKDDGYp33Tuy10S+GKWzZ34Qdxxk5PJWl4zjO\nAPHpCMdxnAHig7DjOM4A8UE44gnOZWsaKgvUrTicZTCzBoCANTBZo97s7CRnTNaShmfpNDOYtTIA\nBUPOImMsLI6xOF9idvNkg06o0xwTmQYzqxWwieNtgvags7C5Zvh5gZNMMls1qMQGTWssi40gNvF2\nqJvVse3ZxA4mmK0aveIw32KdTqVq9LL3BqgzGDXrN2uMCknDQ1nBOGSfgTg8umpkXE6NaeUCi/Ml\nFjfHhrmawSz0W3xvihk61ng2tzxBuVygWKywMF9icX6Mue2T1TLG0trEhrm4r6wBMCt8fXYxMcw1\n0wmGuXJqmLMGM6AaPg1QYrH6zNvnMfy+AJTGF1iYLzE3WW94CwY7a9AM/VxYV6k7JxgA7f2v9nmk\nszg/lqx83ieywvrXKm0j5kTkLBG5J82NeUBE3p8e3yYid4rIYyJyR4gycRzHGTSdRswNgjxhy2Xg\nD1T1POCfA+8VkVeRJKa4W1VfCdwDfLh31XQcx8nPMA3CeVZbPkyS0g1VnRGRg8BZwNuAX09PuwGY\nJhmYHcdxBsrI+gmLyLnA+cADwK40YQWqelhETlv12jmO43TASPoJi8hGktRsH0jfiONlmpsu2/z8\ndf+JGeYpscC2qV9kcuqXO6yu4zijxPT0NNPT06te7lqZashDriXv07WSvg7crqqfTI8dBKZU9YiI\nnA7cq9q4zrSI6JaFQ2wdO1Zdbt0uDZ8VFmutr4G6cOLZSSrlAgtzYyyfmmTLWbXlyLeOhWTrxarF\n2lrSoWblbeVRMLs4weJ8iUpqeWd+rE5n09hM1QsjhDgvLo9RLqfeBMUKhXX1HgX1OgUWl8eqluPy\nUoHlShHmx9hxbs3bY+O6k0wyRwifnWOShcUxKuXUm6BYrmplJd1eXE5Cbht0ZkpsOfdwXXsmma16\nUpxc3EilXKRSLlAoVppqVSjUeSnMzUxSXkr2lxdKUC403J+6sOX0XlqdsfHk/pfGaknBK8uFat9W\nykUW5xPPi3WFcqbOprGZqgfH3PIEszOJTnmpQHF9hUIxbEmbAuVyodq31WesUqzpzBfr+q00tlhd\njn52caLa/nBNcX2FsfHFqo5tQ6X6ucDiTOJNsK60UNXZeNbz1X6fHKsldZ+dnWT25ERyH4F1hTKl\nicW69izMl6plB53yUiEpG6BYgXKBjTtrCc4LxXK1zyvLBWZnJpmbqXlPZLVnIf0dsX1bXiqwfGoS\nipWkPacmoSxoLUd9U1Zryftr9P/Jde5/kP9v4Eve580n/FfAo2EATrkVuDL9/G7gllWsl+M4Tscs\npKvctNvWAm2nI0TkIuB3gAPpuksKXAN8FPiyiFwFPAlc0cuKOo7j5GWk5oRV9X5oOsFyyepWx3Ec\np3s6nRMWkRLwLWCMZHy8SVX3ishvA9cBrwZ+RVUfzrj2LODzwC5gmSTZ+39spzk8fy4cx3Fy0ukg\nrKoLInKxqs6KSAG4X0RuBw4Avwn8lxaXh5iK/akjw3dF5E5VjZdHqqMvg/Dx57dR2VpkcTI1ZKUG\nKWs0AIyxZKJqLAnGrWCYWZwvJcaCcgHmizADx8e3VctIdEpVHWuQCeQxzMyenKgafJgXmG/UqUzW\nhwEHIw7UjBhWp3ptavgJxqVqW8rAPPzsp1tr9dlWYG4sCR0NxrXqNZAYPyKtmk7NYBKMfpSlqnP8\n+fr2zI0nobML8yVmjm1KNMqFxIiTao2NL9RpBB2gZpApp3aOoBP1W2l8oWrEWwz3sixQTIzEc+OL\nrCuUWZxYrNOoGvxCW+aF5WIpU2dxY4mx8YXavQz1KsNikUQrNRwV15v7H/oK6u7LcrEE88AcHN9Y\n0xkbX2BsfJFK2h5m0vqQXLNYVBbT9gSjVbUN4R6GZ6wMy+PFqs7M+NZqvy9urBnaFmcmYaZYpzM3\nnhjCwj1aTA1z1XsISfvn04qnOjPUnrVwbd29TOtVbc+4spjqrCuUa+1I7+FiUROdGaBYTNozAyzR\nV7rxE1bVsCx8iWSMVFV9DEBEmhrxmsRUnEnjGnV1+Juw4zgjRzdzwun6ct8FXgp8SlUf6qCMc0li\nKh5sd64Pwo7jjBzNpiMOTf+Iw9M/anmtqi4DF4jIZuBmEXmNquZZ8h5ojKlod74Pwo7jjByLTdzP\ndkz9E3ZM/ZPq/vf2/o+mZajqCRG5F9hDjqXroRpTcRPwBVXN5bab10/YcRxnaCinKUPbbTEisjNk\nhBSRCeBNNM7ptgruyIqpaIkPwo7jjBxJ5Gv7LYPdwL0isp9kPvcOVb1NRN4uIk8DrwO+nnpMICK7\nReTr6ecQU/EGEdknIg+LyJ52dc0VttwNIqL8QGHrAhu3nqyFgi4YC3/AWF2Dtd9S9T6YqVnEmQdr\n4LU6VUtysL5D1SpudSxV67W1DM8BC8DObJ0QPl31PIBkoifSslStyvOpNbxMYkFeAE6vnbduyylK\nE4u18OnQ9jK1yaSglVrhGzRC+2OdqD1j4wu10NZgfbc64/VtqnWascDPmD5oohPquHxqsnb+ErA+\ntCXdxqOVeoOGvS/rm+hsLCchs+F5iS3060N7Ih37rIRnLNQtPG/m/rCxXA0BZqYIJ6EaPV6IdIrm\nd80+K6EtlfTcoBPaU0x1Qv1mUu+DQKwzXk6eq1C21al6R6SfN5pywv0N9Zuh9sxA4isQtye0wz4r\nVqdAcm/KoG+gLasVtvx2/e+5zr1Z3jnwsGWfE3YcZ+QYpgQ+Pgg7jjNyjGw+YcdxnGFgpHJHOI7j\nDBvNXNTWIv0ZhI8BlGq2hJlSbeK+zpCVWGSWx8eqoZF1mFBllkgMGQ2u0KlOMBZZQxEAkugUjU6M\nvS4YJuZDOzJ0bJ0sE0armD4UwYgD9QamitH6Wa2IZTYwF4VPV/sOakaZokC5VNOBeuOPva6c9lvd\ngrQlFjcWajrHqL83pPU0barpmJ+hLVbTPmXFtI6x4ccaAIMBKA3vrhp7YgNTmarRp06nkFy0bELB\nmYl0itQMbuWoLbFhkUgnPAeF9NrQFyeBU1G/WZ3xyBgXnpdw761+qLPVCczQ+NyPR+0J5VmNcsZm\nCfc31CM2zIW+qxpRjVHV/hplGQBjrR7j0xGO4zgDxKcjHMdxBoh7RziO4wwQH4Qdx3EGiA/CjuM4\nA2Sh3uq8punPIDxDYuENCbhPUrM625qEP17lDAs81Icq288Wq2MtzHGYb8HoYL4LVuMQqmwt8bFO\nuNEnzblWp0K990LRCFkruPUQsNdDzdod6mSt76G4sJVTHYu1XludeeB41KagY70JYp26NplzghXe\n1i++B4Rr0zoGK3/shRFb4W35UN/XzXSgFtJ8ikYvjFD34FUQyLLyh2vC55no/AK1/o3bU6S+Pfb8\nrHDiVjqBY9SHU1vvEasD2f0UiJ/rIvW/l8E7IvSPDUsukjz+sddFHLacpdMHhulNuG0CHxH5rIgc\nEZFHzLFrReSZNEFFriQVjuM4/aJCIde2FsiTRe1zwJszjn9CVV+bbt9Y5Xo5juN0TKepLAdBntWW\n7xORczK+GmjmIcdxnGYMk59wN/mE3yci+0XkMyEJsuM4zlpgmKYjOv1z8ZfAv1dVFZE/AT4B/G7T\ns79wXTLPx5x6AAAaN0lEQVSJPwH8syl49VSHso7jjBLT09NMT0+vermdDrAiUgK+BYyRjI83qepe\nEdkG3AicAzwBXKGqsWkbEfkgyVi4DBwA3qOqi/F5ddfkSeqeTkd8TVV/aSXfpd8rf6NJ4vUN1PI9\nhFh/a622ycOtBR4acyAEi2tI7B2wOva8QF3ScKMTCJbgBeq9MJZorZNl5a8mv4ZMj5m4jsGTYSLS\nKZn62MTsDUnQqfcysRby2ANjPj3f6oyb71omQaf+z3dWTge7H+uE+z5Dfa6JUO+J9Bp7ndWKcx9k\n6Wyk3lIfdMLzVmiiU45+Wu8Im3i9mU78HFidrD6z9Q/XNNMJ58beRfaehPaE3y/rpRB+X0jrFP/u\n2N/B2EMIar8vtj22v6x3RKxTAX0XbVmtpO5bFg7lOvd4aXeDnohMquqsiBSA+4H3A78FvKCqHxOR\nDwHbVPXq6LozgPuAV6nqoojcCPwPVf18qzrkfRMWzBywiJyuqofT3X8FfD9nOY7jOD2nUs47tDWi\nqrPpxxLJGKnA24BfT4/fAEwDVzdcnPyZ2iAiy8Ak8Fw7vbY1FZEvAVPADhF5CrgWuFhEzid55X4C\n+P125TiO4/SLSrx02goQkXXAd4GXAp9S1YdEZJeqHgFQ1cMiclp8nao+JyIfB54CZoE7VfXudnp5\nvCOy/on4XLvrHMdxBkWzQXj5vr9j+f77Wl6rqsvABSKyGfiqiJxH8jZcd1p8nYhsJXljPockFOom\nEXmXqn6pld7w+HE4juPkJCzY28CvTiVb4GMfaVqGqp4QkWlgD3AkvA2LyOnATzIuuQT4saoeBRCR\nvwV+DVgDg/AM9RP+P6M+DNMaZMI54yRGgdiYEQw6wTi1QKORAuoNeJY45NYaI6xOHOJbztDJCj22\ndYhDZOPv7PUh1LRMfd9AfWJs26ZQdmhTVntiA1McwmzbY42dpzLaUE0antGucK4Nc22mYw1zcb/F\n/dPMABiMVzZkOA6NtQazrH6zIbf2fKtjQ43jcGmoN9i2a09sALTlhOuCETb+3uqFMOy4nSVzXWxg\nDNfaZzjrd8cade1zbcO8bXtC+cEQZ1e/tjp9DlternQ2tInITmBJVY+LyATwJuAjwK3AlcBHgXcD\nt2Rc/hTwOhEJo9cbgYfaafqbsOM4o0fnc8K7gRvSeeF1wI2qepuIPAB8WUSuAp4ErgAQkd3A9ap6\nmap+W0RuAvaR/CnaB3y6naAPwo7jjB7znQ1tqnoAeG3G8aMk0w3x8UPAZWZ/L7B3JZo+CDuOM3r0\nefqjG3wQdhxn9PBB2HEcZ4D4IBwRvCOCtfk4jZZvmzCc6HtrhY9DULOs79ZzIstaHcKIs0Kmg/a8\nKSO2OAc2km2ttvUIoaRxe4OObUszi3hs4c9qz3h0DOot5Nb7wobJWkKS7lbJyUMS9Ino2tgLw9bb\nEizwC9QvQx/rWG+M2AvD3pesJOJWy3phBK+NEIIdyrTH4wTv1ksjywMjrsNJ813QCXWJdawXTLzk\nfbPnAGrh3lYn9iaJl7rPoxM/O/Z3zXrg2Ocg1gheEXna00uW2p+yVvA3YcdxRo/YzXMN44Ow4zij\nh09HOI7jDJCs6ak1ig/CjuOMHv4mHBGMFRtJ/kKFsOVghAvGl5ATNl7RN2ANbtbAZI0UVif+zubf\nDUaMEvU5VQNZqyDHOuF4yL1r8+La8N7w0+YutnmLyxlbVl3iMNJQ3jj1hpk4WMgaMW294/bY0Gi7\n0m8o02rY+xev7FuO9mOd0CdZeXGDTrjvNnQ9K8zXhhrHOpjjNnTWGv2gZjyLjVi2PdYwtxSdY3Po\nnjLfBR0b4lukMYw3vudhP9YJ2PBoG1IehxVnhWHbPsjqswlzTtAJ51sd2yZbdqwZ7/cLH4Qdx3EG\niA/CjuM4A8Rd1BzHcQaIu6g5juMMEJ+OcBzHGSDuohZhE4TPkYQtZ60WHFvg8yRcjy3SVieELQey\nQpatlr0+1gnhya104jbZJOg2fDXWsVusE8qIV1qGmmdE8CqJ+y22uFvviKz2WJ04XNUmQW8Vhh1b\n41vpWK+SWCcOw27ltZClYxO+Wys/1J4z60kQtOy+1YoTxAcWmnxXpvac2X6D5u2xxDrWA+WkKdP+\n3sRJ6rPaYp8daLw3C+bcObNvw+9tXeI22PsRdOIQ7X4wRG/C69qdICKfFZEjIvKIObZNRO4UkcdE\n5A4R2dLbajqO46yArBecrC1CREoi8qCI7BORAyJybXo895gnIutE5GERuTVPVdsOwiSLer45OnY1\ncLeqvhK4B/hwHjHHcZy+0OEgrKoLwMWqegFwPvAWEbmQlY15HwAezVvVtoOwqt5HEl5heRtwQ/r5\nBuDteQUdx3F6zlLOLQNVnU0/lkgmWJScY56InAVcCnwmb1XzvAlncZqqHkkrfBg4rcNyHMdxVp9K\nzi2DdDphH3AYuEtVHwJ25Rzz/hT4Q5KBOxerZZhrLXjXdTUDyxlTMDmVHLdGrHHqjT/FaINsI1YF\nOGa0yiTGhBASa40usRHDhsjGPRGvTlzO0DlFfUisNcxZ49+S+RwTt6VMfchwqLttj9Wxxr/weTy6\n3mrQRCfUL7THhrNOUB9+Ghsgmxl/mvVb0DlF/dtIHIJtdeL2WINerBPnn56n0QBojVjN8haTQ8eG\niIccybae1ohlQ6WzNGihE+6HXQU5nG/LDPUPWvHbXiX6bO+NbU+WTjA2BkJOYWt8C9fbcHMbmh4x\nPT3N9PR08xM6pZl3xBPT8GRrPVVdBi4Qkc3AV0XkPBrHuIYxT0TeChxR1f0iMgVInqp2OggfEZFd\nqnpERE4HftLy7H9+XXLDgiX/px2qOo4zUkxNTTE1NVXd37t3RWtkNqeZd8RZU8kW+FZzPVU9ISLT\nwB7yjXkXAZeLyKUkry6bROTzqvpvWlU173SEUD+q3wpcmX5+N3BLznIcx3F6T4dzwiKyM3g+iMgE\n8CbgIDnGPFW9RlXPVtVfAN4B3NNuAIYcb8Ii8iVgCtghIk8B1wIfAb4iIlcBTwJXtCvHcRynb3Tu\nl7wbuEFE1pG8pN6oqreJyAPAl+MxT0R2A9er6mVNS2xD20FYVd/V5KtLOhV1HMfpKR0Ga6jqAeC1\nGcePkjHmqeohoGEAVtVvAt/Mo+lhy47jjB5DFDHXn0H4GPWGuWPUQlOtd0RWcvJmid2tdTu2SI+T\neBME63vshWE9MKzF3BJrkKFjvQli7wQb9mu9MfLo2BV7g04z74iNNCb1zvJ/tAnLg5a1jId7EXtH\nrE/3g7eH1YmTxzfTadZv1jvCentYnSzPlVCO1Yh1godI8HKJvT1WomPDcWOdcE/DYgNz0bVL5jyr\nY712yk324/aEsmKvkrDytQ33DvvxYBTrxJ44Wd4R9liZxuc19kIivc62td+DoqeydBzHGSAt3OLW\nGj4IO44zevh0hOM4zgDx6QjHcZwB4itrOI7jDBCfjog4SS3x+EngBWpx6MGbYCPJvxDzwAYaE74H\nrKU1WJ+tN0HQORYdDzphC14FcZ6KWMfG4sc6RWoJw+NE6+NmP86HYSlHP5eo91oIBoZmOmXqk223\n0rEP5krbY8u3OlnYtkC2F0ZenXB+M40snQK13AE2d0TQWSC5/+FtKXiu2OXe4wTvWblKgtdASFIf\ne2FY74QJc9zqxG2weRji9oTjcXsqpoxyRlnWIyf+zt4b6+2SpYMpP0vHnhci0jype1v8TdhxnNHD\n54Qdx3EGiLuoOY7jDBCfjnAcxxkgPh0RcYyaoW2GJF891BvKgqEuGDBsmDHmWLuk4UHnGLWQzBBa\nuiG63oaStgoppokOqcacOc+2aYn6ENksHduWLJ1Tpj3xCsjWoBmH4GaFFLfSse0pUwtX7UTHGmHK\nNBoAw89mBs1wn8L14SmNDUeWWCe0p5kO5phdUCBQIFuDSMfWp5khK3xvQ+Vjw2Z8Xl6dCfM5NpqG\nzzbBeqyTpRHr2NQBWaHQ9li8EnPW+f3AXdQcx3EGiE9HOI7jDJAhGoQ7XejTcRxn7dL5yholEXlQ\nRPaJyAERuTY9vk1E7hSRx0TkjrD6Rsb1e0TkhyLyIxH5UJ6q+iDsOM7osZBzi1DVBeBiVb0AOB94\ni4hcCFwN3K2qrwTuAT4cX5uuxvEXwJuB84B3isir2lXVB2HHcUaPcs4tA1WdTT+WSKZsFXgbcEN6\n/Abg7RmXXgg8rqpPquoS8NfpdS3pak5YRJ4AjgPLwJKqXph5YvCOCJ8PUwsTLQGbSDokhJLapOFx\nGCnmWPiZZeX/KYlnAaas+VSjRC0k1ibCLpqf9gYFj4BYZynVmDHX2bBoTHusV0OrMNwsndAea/m2\nXiW237Lak1cntMcuF269FUK/xe2JteKQ4iydeHn42CsiWLdjj5LYyt9MJxASoMdeJbYM26fNlqVv\np2MT7kNjuG4zj5JYq9n9gfpk6kEnJNwnOk4Tndh7I6stsU54Zstk36M8Ov12GetCL32j/S7wUuBT\nqvpQWGkZQFUPi8hpGZeeCTxt9p8hGZhb0q1hbhmYUtWfdVmO4zjO6tGFi5qqLgMXiMhm4Ksich7J\n23DdaZ0r1NPtICz4lIbjOGuNZt4R5WmoTOcqQlVPiMg0sAc4Et6GReR04CcZlzwLnG32z0qPtaTb\nAVSBu0TkIRH5vS7LchzHWR2azgFPQeG62hYhIjuD54OITABvAg4CtwJXpqe9G7glQ/Uh4GUico6I\njAHvSK9rSbdvwhep6iER+TmSwfigqt7XZZmO4zjd0fmc8G7ghnReeB1wo6reJiIPAF8WkauAJ4Er\nAERkN3C9ql6mqhUReR9wZ3rtZ1X1YDvBrgZhVT2U/vypiHyVZBK6cRD+h+vSC4DyFMxMJfsT1Ixj\nUD/xHwwCWZP/Fpt31XKExPBjQ5aDwSEYzeIQ3GY6oQ6xzgK1cGKMTjA0Ysq3KzGvVKcCPG/KLFIz\naIZ22H6z561U5xi13MKkOnZF6KAXQrJtSGys1UonrB4d6hmMYtaQFcrM0sma82umE+oe8ldvSL+3\nxlNrBLZaWWTp2Hy/1pBon+2s/Mit5i6tTpZRLtQ1NgJDowE41ml1b2Id+9zaa+OQ5lY6TZienmZ6\nerr1SZ3QYbCGqh4AXptx/ChwScbxQ8BlZv8bwCtXotnxICwik8A6VZ0RkQ3AbwB7M0/ecV3ys9mA\n6TjOi5KpqSmmpqaq+3v3Zg8ho0w3b8K7SCyHmpbzRVW9c3Wq5TiO8+Kg40FYVf+RJKLEcRzH6RBP\n4OM4zggyPAmFfRB2HGcEGZ40av1L6h4syHOQuC5MwNzmZD/0V2zlLVIfRtosDDfLinwS4ASUJ6C8\nvj7xenxuqxDcVjonQ93TStoVeQNxKHYn7an2m7nYJniPE2h32p6qjmHOfN8qFDtLJyusPIT3xjoh\nOXmcFDx4yNiVje0q2c10stozQ+0eWLI8ceJwYhvKnqUzTy2Gqmiu30DtHmV54uTVsWHRtp5xEhrr\nGZNHp5UHRrxQQcAmkQ+eOvbcZu3pK/4m7DiOM0Div/JrFx+EHccZQfxN2HEcZ4D4nLDjOM4A8Tdh\nx3GcAeJvwvUcw/TJMyTZ3SZIgu42w8mJ7ETqJWoWWJsHISTMtuXb6wB4AjhKkmBhc7Idm6g/J87r\nsGKdORJPj7m0PRNQ3tUYmr1kfmbpxLxgPmu48FlqCRs2Q3lzohPnUWjVnlw6R6P2bK5vj31i8upk\n3p+jpoAJmJuozxcRzrMeJWEZ96ATv+wEnQom2+vRWsV1AuZSNwHrLRA8V2wbbZ8WyH4O6nROUC24\nPFHvBRE8L7KS7q9Yx7hYlNc3Lllvk+E30wlaViNTJ2VufeL9Ea6xOURKtB7vmj0TPcffhB3HcQaI\ne0c4juMMEJ+OcBzHGSDDMx3hSxM5jjOCdLbcsoicJSL3iMgPROSAiLw/Pf5PReR/isj3ROQWEdnY\ncHFy3hYR+YqIHEzL+NV2Ne3Pm3AZEgPW0fTnsyTGsmBF2gUn1yeT+DaZuF112RqAYkNUXV9ajRPA\ndpL5oSXQs2qGjBL1xKsh59J5hmSpqaW0PZuSAnR7olOg3jAx0UQn7iu1f8WPpltYqmoT1Vhv3V5b\nUTrohHZlre4cr4hcp3OSxFIX2rM9bROJcS6Ex8YrYAdDkzVExTp1/RYMmcEwt5mqkbY8keicisqx\nBrk4ibgNXa7qLKXtOZHqFKkaNNmeGJqC0SprZeVYJw6RrtOZS3XCgxV0NsHM+uQ6u6p0Mx1rgG7a\nntlUL1yUGpxnqC3ObsuNdQLNVsautids4cQdiVHzlDk1K5F7bIRrtjJ2X+j4TbgM/IGq7k8H2u+I\nyF3A9enx+0TkSuDfAf9vxvWfBG5T1X8tIkVgsp2gvwk7jjOCdPYmrKqHVXV/+nkG+CHJUvYvN0u3\n3Q38Vnxtujrz61X1c+n1ZVU9EZ8X44Ow4zgjyFLOrTkici5JzvQHgB+IyOXpV1eQrKQc8xLgeRH5\nnIg8LCKfThcLbYkPwo7jjCBzObds0qmIm4APpG/Evwu8V0QeIsmNt5hxWZFkfbpPqeprSeaPrm5X\nU/eOcBxnBGn2lvvDdGtOOpd7E/AFVb0FQFUfA96cfv9y4K0Zlz4DPK2q30n3bwI+1K6mPgg7jjOC\nNLMGvizdArdknfRXwKOq+slwQER+Ll1Vfh3wx8B/ji9S1SMi8rSIvEJVfwS8EXi0XU27GoRFZA/w\nZyTTGp9V1Y9mn3kEOAA8TtI5wWshDY1lEsrbE8v4BPX9166G1STxAasD9et2z9VCSuNk2K1mbqrL\nfzfTmaDqgRGs/eWJxGpsl3UP0cBtdZ41Bw+a/SKww7RnAubTkO8F6i3t7ZiDWkhvlk6oUGqFt8nd\nyyRhrMHDwIa02suq7Yl1wv769KQQWj6R6JRMe6AxUTjUJ8dv0Alh3keo3Z/Q8euB7bU6z0flxmHZ\nwSsjPJcNOsEbJ+hsql2v22v9NE9tkYL4mc6tE/6FnqDWwesTzwXbX1leKoFYp+7eWO+IE1SfsVBw\neNbWU/Nkss+0fe4WSO5j+NXoO515R4jIRcDvAAdEZB9JIPc1wCtE5L3p/t+q6n9Nz98NXK+qYdn7\n9wNfFJH1wI+B97TT7GbJ+3XAX5CM9s8BD4nILara+l2/bzwOvLzPmvuBc/usCZSngan+6x6fhh0D\n0H1hGja/WHTvA9q6mvaA7wKvG4DuatGZX5yq3k92xotvAP8x4/xDwGVm/3vAr6xEsxvD3IXA46r6\npKouAX8NvK2L8laZxweguX8AmkBlejC6Jwake/TFpHv/ADQBHh6Q7mrRvXdEv+hmOuJM4Gmz/wzJ\nwOw4jjNgPHeE4zjOABmeLGqiqu3PyrpQ5HXAdaq6J92/GtDYOCcinQk4jvOiRFWlm+tF5AngnJyn\nP6mq53aj1y3dDMIF4DESw9wh4NvAO1X14OpVz3EcZ7TpeDpCVSsi8j7gTmouaj4AO47jrICO34Qd\nx3Gc7ulZ7ggR2SMiPxSRH4lI29C9FZb9WRE5IiKPmGPbROROEXlMRO4QkS3muw+LyONpjs/f6EK3\nWa7RnmmLSElEHhSRfanmtf1qb1rWujQZya390hWRJ9K8rftE5Nt91G3IBdtrXRF5RdrOh9Ofx0Xk\n/X3Q/aCIfF9EHhGRL4rIWJ/6+APpc9yX35+hQFVXfSMZ3P+eZHJ8PYkD7atWsfz/jSS70SPm2EeB\nf5d+/hDwkfTza4B9JFMv56b1kg51TwfOTz9vJJkTf1WvtYHJ9GeBJKPThf1ob1reB4H/Btzax37+\nMbAtOtYP3f8KvCf9XAS29Kufze/Nc8DP91IXOCPt47F0/0bg3X14js8DHiGJpSuQTGW+tJ99vBa3\n3hSahNrcbvavBj60yhrnUD8I/xDYlX4+HfhhljZwO/Crq1SHm4FL+qVNkiD6OyQROT3XJEnXdxdJ\nOF4YhPuh+4/AjuhYT3VJ4qb/IeN4354r4DeAv+u1Lskg/CSwLR3gbu3Hcwz8NkmIb9j/Y+APSWLZ\n+/q7u5a2Xk1HZAVynNkjrcBpqnoEksTMwGlN6vLsatRF6nON7uqldjolsA84DNylqg/1WjPlT0l+\nSazhoB+6CtwlIg+JyP/ZJ92sXLCTfdC1/O/Al9LPPdNV1eeAjwNPpdcfV9W7e6mZ8n3g9en0wyRw\nKclbfz/7eM0xyvmEe2ZxlMZco7HWqmqr6rKqXkDyZnqhiJzXa00ReStwRJNVBlr5bfainy/SJB/r\npSQ5XF+fobPaunEu2FMkb2K91gVAkoQvlwNfaaKzaroispUkxcA5JG/FG0Tkd3qpCaBJXpmPkvx3\ndRvJVEMl69TV1F3r9GoQfhY42+yfRX1qsF5wRER2AYjI6SSLpYW6/Pxq1UUyco32S1uTpVKmgT19\n0LwIuFxEfgz8d+ANIvIF4HCv26pJUhRU9ackUz4X0vv2xrlg/4ZkUO7LvQXeAnxXVZ9P93upewnw\nY1U9qqoV4KvAr/VYEwBV/Zyq/rKqTgHHSOwq/erjNUmvBuGHgJeJyDkiMga8g2TeaTUR6t/QbgWu\nTD+/m1qi0FuBd6TW35eQJBP9dhe6DblGe6ktIjuDtViSpVLeRDKH1tP2quo1qnq2qv4Cyf27R1X/\nD+BrvdQVkcn0Pw1EZAPJPOkBet/eI8DTIvKK9NAbgR/0WtfwTpI/doFe6j4FvE5ExkVEqOW97Xlb\nReTn0p9nA79JMv3Srz5em/Rqspnkbe0xknRmV69y2V8isSIvkDxQ7yExMtydat4JbDXnf5jEsnoQ\n+I0udC8i+fdpP8m/Ug+n7dzeK23gF1Od/SSW5T9Kj/dMM6MOv07NMNdTXZK52dC/B8Kz04/2Av+U\n5AViP/C3JN4R/dCdBH4KbDLHet3P16bXPwLcQDXRcs/b+i2SueF9wFS/n+W1uHmwhuM4zgAZZcOc\n4zjOmscHYcdxnAHig7DjOM4A8UHYcRxngPgg7DiOM0B8EHYcxxkgPgg7juMMEB+EHcdxBsj/D5Qe\nz20KiwefAAAAAElFTkSuQmCC\n", "text/plain": [ "" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "plt.pcolormesh(sal_abs[0,:,0,:])\n", "plt.colorbar()" ] }, { "cell_type": "code", "execution_count": 15, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Output_derived.txt Output.txt\r\n" ] } ], "source": [ "!ls *.txt" ] }, { "cell_type": "code", "execution_count": 16, "metadata": { "collapsed": false }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ls: cannot access *.m: No such file or directory\r\n" ] } ], "source": [ "!ls *.m" ] }, { "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.5.2" } }, "nbformat": 4, "nbformat_minor": 0 }