NOAA tides from R.ipynb

{
 "metadata": {
  "name": "NOAA tides from R"
 },
 "nbformat": 3,
 "nbformat_minor": 0,
 "worksheets": [
  {
   "cells": [
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%load_ext rmagic "
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 1
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "X = np.array([0,1,2,3,4])\n",
      "Y = np.array([3,5,4,6,7])"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 2
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R -i X,Y -o XYcoef\n",
      "XYlm = lm(Y~X)\n",
      "XYcoef = coef(XYlm)\n",
      "print(summary(XYlm))\n",
      "par(mfrow=c(2,2))\n",
      "plot(XYlm)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "display_data",
       "text": [
        "\n",
        "Call:\n",
        "lm(formula = Y ~ X)\n",
        "\n",
        "Residuals:\n",
        "   1    2    3    4    5 \n",
        "-0.2  0.9 -1.0  0.1  0.2 \n",
        "\n",
        "Coefficients:\n",
        "            Estimate Std. Error t value Pr(>|t|)  \n",
        "(Intercept)   3.2000     0.6164   5.191   0.0139 *\n",
        "X             0.9000     0.2517   3.576   0.0374 *\n",
        "---\n",
        "Signif. codes:  0 \u2018***\u2019 0.001 \u2018**\u2019 0.01 \u2018*\u2019 0.05 \u2018.\u2019 0.1 \u2018 \u2019 1 \n",
        "\n",
        "Residual standard error: 0.7958 on 3 degrees of freedom\n",
        "Multiple R-squared:  0.81,\tAdjusted R-squared: 0.7467 \n",
        "F-statistic: 12.79 on 1 and 3 DF,  p-value: 0.03739 \n",
        "\n"
       ]
      },
      {
       "output_type": "display_data",
       "png": "iVBORw0KGgoAAAANSUhEUgAAAeAAAAHgCAIAAADytinCAAAgAElEQVR4nOzdd1gUV9cA8LNLWfrS\npYggqFgQCyqIIIiiKIqg2BAVY4k1SsQSxBJsRIwmQRALVnw1YJSgaCQRCyhGUUEJICJNQHoRkc58\nf0yyH8KCsOzsDsv5PXnyDOPsvWfYu4fZO3fuZRAEAQghhOiHKewAEEIIcYcJGiGEaAoTNEII0RQm\naIQQoilM0AghRFOYoBFCiKYwQSOEEE1hgkYIIZrCBI0QQjSFCRohhGgKEzRCCNEUJmiEEKIpTNAI\nIURTmKARQoimMEEjhBBNYYJGCCGawgSNEEI0hQkaIYRoChM0QgjRFCZohBCiKUzQCCFEU5igEUKI\npjBBI4QQTWGCRgghmsIEjRBCNIUJGiGEaAoTNEII0RQmaIQQoilM0AghRFOYoBFCiKYwQSOEEE1h\ngkYIIZrCBI0QQjSFCRohhGiqJybo3r17M/4jLy9vb2+fl5fX8ZfHxcWNGjWqxc7y8nJFRcVOhREf\nH29kZNSpl3yRqqoqoxk5Obnm0YqLizc0NDTfaB8PJ4W6i/z8fAaD8dNPP3H23Lhxw9ramqLquLZ2\ngiCCg4PNzc0VFBQGDhz4zTfffPjwofVrO3iYSOqJCRoAbt68WVZWVlpa+vz58w8fPmzfvr3jr+3b\nt6+3tzd1sXVRVFRU2X9ycnJoHi0SIgaD4e3tnZOTI6wAdu7cuW3btg0bNrx69er06dPp6emjR4/+\n9OkTb4eJpB6aoOXl5RUVFZWUlPr37+/q6pqenk7uj46OHjFihKysrJ2dXW5uLgA0NDSsXr1aSUlJ\nVVV1z549AJCRkbFz507yeD8/Px0dHR0dndOnT5N7Hj9+bGZm1nr75MmTffv2lZaWNjMze/36dfNg\nWlfBYWdnd+LECXLb19d3/vz57Rzc/NQ4ONFOnjy5sbHRwMBg3Lhx5EZVVVXr8+V6UkgkSUpKLl++\n/Jtvvmn9T1evXjU0NGSz2bNnzy4qKgKAxMREa2vrvXv3Ghsbp6SkjBs3zsPDQ1VV1cLCIjY2dvTo\n0fLy8hs3biRf3k5r50hLS/vxxx8jIyPnzZunq6trbm4eHh6uqKj4448/8nCYyCJ6Hm1t7ejoaHI7\nNzfX0dFx3759BEEUFxerqKiEh4eXlpauXbvW2tqaIIiQkBBDQ8OMjIznz5+zWKy0tLSnT5+amJgQ\nBPHgwQMlJaX79+/n5OTY2Niw2WyCIGJjY01NTcnCOdvZ2dmSkpL3798vKipyc3NbuXIlQRAvXrwY\nMmQI1yo4oR47dszBwYHcNjc3v3btWjsHEwShoqLy9OnT5ns40RIEISYmVl9fz9nger5cTwqJnvfv\n37NYrMrKSh0dnfDwcIIgrl+/bmVlRRBEeno6m82OjIwsKSlxc3ObN28eQRCvXr1is9lubm6JiYnJ\nyclMJvPixYslJSUmJibq6uqZmZmxsbEAUFhY2H5r5wgKCpo4cWKLqIKCgsh22NnDRJW4sP9ACIed\nnZ24uDhBEB8+fDAzM9u6dSv81wc3Y8YMADh06JCKikpjYyMA1NfXFxYWjh49OicnR0FBoaysjCzk\nt99++/rrr8ePHw8Ae/bsmTZtWlvVqampvXnzpk+fPlVVVaqqqu/evWtxQIsqOPtnzpy5efPmmpqa\n8vLypKQkOzu769evt3UwydraWlz837c1ICBgwIABbUXF9Xw7flJIBMjJyR09enTt2rUTJkzg7Pz9\n998dHR1tbW0B4ODBg9ra2uQHobq6OjAwkMVipaSkaGpquri4AMCkSZPKy8t1/1NRUdG7d+/2Wzsp\nMzOzT58+LXbq6upmZ2fzcJio6qFdHGfOnImPj09ISHj69GlxcXFwcDAAvHv3LjIyUk9PT09Pb+DA\ngZKSkoWFhbNmzfL09Fy5cqWWllZgYCDZUkn5+fkGBgbktr6+futaCIIgN8TFxU+dOjV69OiZM2cm\nJSW1OKydKjQ1NY2MjO7duxceHj5z5kwpKal2Dib973//i/+Pg4NDO78Eruf7xZNCIsbBwcHExGT3\n7t2cPfn5+Xp6euS2mpqapKQk2cuho6PDYrHI/XJycuSGuLi4hoYGZxu+1No59PT0WifZrKwsAwMD\nPz8/soPu9OnT7RzG4wl3Kz00QWtqapKJadSoUbNnz37x4gW509bWNjMzMzMz8+3bt3/++aeGhkZ6\nerqNjU18fPzff/99/fr1oKAgTiFaWlpv374ltzMyMjj7OQMkOLdfQkNDIyIibt++/ddffy1YsKBF\nMO1UAQBOTk4RERFhYWHz58//4sFkVHr/4XyK2voltD7ftk4KibCff/751KlTCQkJ5I8aGhpZWVnk\ndlFRUW1traqqKvyXf7+o/dbOYWlp+ejRo+TkZPLHs2fPPn36NDAw0N7efv369eXl5eXl5V999VU7\nh/F8vt1ID03QzWloaJDfwuzt7aOjo2/evFlcXLxt27aNGzcyGIzw8PAFCxYUFBQ0NjbW1tZKS0tz\nXujs7Hz8+PHo6Oi8vLydO3cyGAwAYLPZCQkJ8fHxJSUl/v7+5JElJSVycnLS0tKFhYV+fn7V1dXN\nA2inCgBwdHS8du1aQkLCxIkTv3jwF338+JGzwfV8uZ4UEm19+vTZuXPnvn37yB9nzJhx9erVO3fu\nlJWVeXh4ODk5dTA1k9pv7Rz9+/f/9ttvJ0+eHBoampWVVVlZaWZmVlRUtHbtWh4OE1nC7gQXguY3\nCQmCiIiIUFdXr6ioIAji9u3bQ4YMkZGRmTBhwtu3bwmCqKysdHR0lJWVVVZWXrVqVV1dXfPbbn5+\nfr1799bW1j59+rS2tjZBEE1NTevXr5eTkxs6dGhoaCh5k7C8vNzW1lZZWdnc3Pz69eu9evU6f/48\n57ZJ6ypaBDx48OCvv/6a3G7/4PZvEs6dO1deXv7jx4+cjdbny/WkkOghbxJyfqyvrx82bBh5k5Ag\niCtXrgwYMEBeXt7R0bGgoIAgiFevXhkaGpL/mpyczNnevn37rl27yG0DA4M3b96039qba2pqOn/+\nvKmpqaysrJaW1oYNG0xNTQ8cOMDbYSKJQfzXT4oQQsL16dOn9PT0Lz7A1cHDRAAmaIQQoinsg0YI\nIZrCBI0QQjRFuwdVCgsLw8LChB0FEhwmkzl37tzWT9yIBmzPPQ1/2zO1Cbqurk5SUrJTL7l161ZC\nQgJ1s2ohugkLC9PR0ZkyZYqwA/kygiDIGf46PvoQ23NPw9/2TEmCTkpK2rJlS2xsrIyMTGVlpYWF\nxeHDh9t55riFUaNGzZkzh4rAEA2lpaUJO4QvWLNmTUBAQHx8/Lx585hMJpPJDAoK4kyD9UXYnnsU\n/rZnShL06tWr9+zZY2pqymKxampq4uLi1q1bFxkZSUVdCFHtyZMnALB9+/bQ0FBjY+O8vDwXF5d7\n9+4JOy4k+ihJ0Ewm09LSkvwaKCUlZWFh0aknkRCiIWlpaWNjYwDQ0tJqPQUKQlSgJG+amJhMnTp1\n7Nixenp6WVlZjx8/Jls2Qt1Renq6k5NTdnZ2cHCwq6vrwYMHtbS0hB0U6hEoSdC+vr7R0dExMTEJ\nCQnKysqenp7jxo2joiKEBKCoqCgzMzM1NVVZWRkAJCQkOKsoIEQpShI0g8EYP348OacwQt2dmJiY\ngYEBZ35Ld3f3to6sq6t79epV8z1paWmamprUxodEl4C6hnV1dTkTGDaXlJS0YcOG5nvevn07ceLE\npUuXCiYwhHjQVnt+9+5daGho8z23b9/u16/fmjVrBBUaEikCStBRUVFc9w8ePPjPP/9svmfFihVN\nTU0CCQohHrXVng0MDHx8fJrvKSkpwfaMeCagR717yPIHqIfA9owEA0e/cREVFfXkyRNjY+OpU6cy\nGIz4+HgfHx8FBQUrK6uFCxcKOzqEkHC0SAXv3r17+vSpkpKStbU1RUtbUJKgDx061Hqnh4cHFXXx\nnbe3d1ZWlrOzc2RkZEREhL+/f1xc3C+//KKiojJr1ixM0D1Qt27PiI+ap4I+ffp4eXk5OzvHxsYG\nBASEhIRQkaMp6eJQUlLy9PSUkJCQaoaKiviuqanp9u3bp06dmmpicuTIkdTU1A8fPixfvpzJZHp4\neKxatUrYAaJOKy0tff36dX5+/v79+1NTU3koofu2Z8RfzVPBgQMHrly5sn79el9fX21t7djYWCpq\npOQKetmyZZcuXVq7dm23e4CwsbFRWlqa8c8/MGYMvHnDZrNramoyMzPPnj27detWzurFqBtZtmyZ\nm5vbrVu3DAwMVq5cycMj2t23PSP+evnyJScV/PDDD+RaugDQu3dvcuFzvqPqJuFff/3VHVuzhISE\nlrp64ezZH7/7Lnnr1k+fPqmrq//888/Z2dkeHh67du0SdoCo0+rr6x0cHLKysjZv3szzgIpu2p4R\nfzVPBVZWVr6+vk1NTZmZmdeuXaPoWTxscy0FGRvfLS/3f/TI79WrS4mJABAUFCTsoBDvxMTEvLy8\nxowZEx0dXVFRIexwUDfWPBXU19fv27fP1tZWTk7ul19+4VxN8xcm6M+9fSvx11+TIyMnM5mwbh2k\npECHZ5VE9OTj43P79m03N7eQkBD8W4v4RUJCYvfu3VTXggm6GYKAtWvhyBFgMgEAFi+Gs2cxQXdf\nXl5enG1yJEZ2dvaoUaOEFxFCnYMJupngYDAxgaFD//1xzBhwd4faWmCxhBoW4tHw4cOFHQJCXYIJ\n+j+FhXD8ONy589nOGTPg+nVwdhZSTKhLnFu9cXFxcUKJBCHeYIL+z6ZNsG9fy4vlRYtg7VpM0N1a\naGjohQsXmpqaCILIyspKTEwUdkQIdZSA5uKgu+vXQUoKrKxa7tfWhtpayM8XRkyIP44fP75t2zYt\nLa1NmzbZ29sLOxyEOgETNMDHj3DgAHB7nBcAYOFCuHxZsAEhfmIymebm5iwWy8bGhv4L1CLUHCZo\ngJ074dtvgc3m/q+zZsHVq4INCPGTlJRUSEhIXV2dv79/VVWVsMNBqBN6fIKOjYWsrPZ6mWVkwMAA\nXr4UYEyIny5fvjxy5Mjdu3cXFBTgs6Coe+nZNwnr62HrVvj11y8ctmQJnD/fZh8IojcZGZl+/foB\ngLe3t7BjQahzenaCPnQI5s+HLy4ZZ2UF330HDQ2AszF0Q82fTJGRkXnw4IEQg0GoU6jNOARBlJeX\nKyoqUjSbdZckJ8OdOxAZ+eUjGQywtYXISJg2jfqwEJ/FxMQAQENDw507d548eUJpXS9fvty0aVPz\nPcnJySYmJpRWikQYJQl6zZo1AQEB8fHx8+bNYzKZTCYzKCjIjFbPTBMEbNwI/v7/PtX9RYsXw44d\nmKC7I87czTNnzjx9+jSldRkbG+Mam4iPKEnQ5HXK9u3bQ0NDjY2N8/LyXFxceJiHl0JBQWBuDoaG\nHT2+Xz8oLITyclBUpDIsxH+cVVwLCwtxFAfqXijs4pCWljY2NgYALS2txsZG6irqtNxcOHMG7t7t\n3KvmzoWQEFi5kpqYEFV69+5NbvTp02fz5s3CDQahTqEkQaenpzs5OWVnZwcHB7u6uh48eFBLS4uK\nini0aRMcPgySkp171fz54OyMCbrbcXV1FXYICPGIkgRdVFSUmZmZmpqqrKwMABISEidOnKCiIl6E\nhYGqKpiadvqFbDaoqMDr153oGEFCZWRkBAC1tbVlZWUsFqu2trZ///4UrR2HEBUoeVBFTEzMwMBg\n6tSppqamAODu7s5u6zk9ASsvhwMH4MABHl++aBFcvMjXgBCFEhMTExMTzczMwsLCsrKyrl27ZmBg\nIOygEOoEAQ3s1dXVzcrKar0/Ly/vxo0bzfekpKRoa2tTFcf27bBrF8jL8/jyKVNg717YvbujYz8Q\nDeTn51tYWACApaUlPkmIuhcBJeioqCiu+yUlJZWUlJrvYbFYTIrSX0wMlJR0aaicuDiYm0N0NJd5\n7xBdMZnMwMBAKyurqKgoaWlpYYeDUCdQm6Dr6uokJSUBoK2vlqqqqnPmzGm+JzIykpJxozU1sGUL\nH6Y9WroUDh/GBN2NXLhw4YcffoiIiBg8ePC5c+eEHQ5CnUDJtWpSUtL06dNVVFQMDAwUFRWnT5+e\nmppKRUWd8MMPsGwZaGh0tRwjI0hNhY8f+RETotbevXsBICAgQF5e3sTERFpa+ujRo8IOCqFOoOQK\nevXq1Xv27DE1NWWxWDU1NXFxcevWrYvsyEPVFPnnH3jyBHbu5E9ps2dDWBjg4C3aGzBgAPw3lgMh\nShEEceHChUmTJvG3WEquoJlMpqWlJYvFAgApKSkLCwtxIU4z1NQE69bBTz8Bv+YDWbgQLl3iT1GI\nSnPnzgUAGxuboUOHWlhYpKamkk9OIcRHt2/fTkxMZDAYixcv5vsDH5TkTRMTk6lTp44dO1ZPTy8r\nK+vx48fC/GAEBsLkydC/P98KVFcHcXF49w50dPhWJqLMsmXL3Nzcbt26ZWBgsHLlyq5MOUDryb+Q\nYMXGxtbW1lpbW0+ZMoW6Wii5gvb19fX09JSQkEhISGAymZ6engd4HnrcRTk5EBICfH/A19UVB0R3\nF/X19Q4ODllZWZs3b+bt/vOaNWsAID4+fuDAgebm5kZGRo8fP+Z3mKh7ePv2bVhYGACMHTvW2tqa\n6uoouYJmMBjjx48fP348FYV3ztq14OvL/3mcZ8yAiRNh2zY+F4soICYm5uXlNWbMmOjo6IqKCh5K\n6AaTfyGKlZaWXr9+feHChQYGBoJ83InLFXRpaenr16/z8/P3798v/NEXXXHlChgYwOjR/C9ZSgqM\njeHpU/6XjPjNx8dHTU3N3d09OTk5KCiI53LoO/kXokxDQ8OlS5eKi4uVlJRcXV0Ffy+NS4JetmxZ\nSkrK7t27JSQkVnbfuYHKyuDIEdi7l6ryFy8GHFTbHQwaNGjjxo1kY26+ukrHkZN/ZWRkBAcHAwDt\nJv9CFLh582ZaWpq4uLizs7OqqiqDwRATExN8GFz+IJB9doGBgYGBgdevXxd8TPyxZQvs2AEyMlSV\nP3YsbNoEtbXAYlFVBeKHBw8eLF26VF1d3d7eftCgQbNnz+5sCbSe/Avx1fPnz+vq6szMzKysrGRl\nZQFAQkJCiPFwuYLuep+d8N25AzU1YGdHbS329nDzJrVVoC7bvn37o0ePDA0N169ff4inlX/pO/kX\n4pOMjAxyOgpDQ8MxY8YAAJmdhY7LFbSPj8/t27fd3NxCQkK60mcnNNXVsHMnhIdTXtGiRbBxIzg5\nUV4R6oLGxkZ1dXUAYLPZ/PrUtTX5V1pamq+vb/M9MTExgwcP5kuliO9KS0tjYmJmzJihpqZGdlvR\nJC9zfJagvby8ONvktUZ2djZv3XbCtHcvrFoFKiqUV9SnD3z6BEVFoKZGeV0ioba2Njk5WU1NjcIJ\nC1sZOHDgmjVrsrOzvby8+NV33NbkX1paWi1u2xQWFsrJyfGlUsQvdXV1UVFRlpaWMjIyVlZWDAaD\ntu/RZwl6+PDhwoqDb168gIQE2LdPQNUtWACXL8P69QKqrjv7559/VqxYYWpqmpmZaWBgwFtvAw82\nbNjw999/S0hIaGlp7dixoytFfXHyLxkZmRZreKuqquKisXxXXl4eEBDw7t27iRMnzp49u4PPDREE\nERsb27dvX1VV1WHDhpEXy5w1henpswTt7Ozc4p/j4uIEGEyXNTXB5s1w7JjganRygpkzMUF3xI4d\nO86dO9e/f38AWLRo0YsXL0aMGCGAel1dXZ8/f96VWz1JSUlbtmyJjY2VkZGprKy0sLA4fPgwOdEH\nErz6+npHR8evv/7aycnpl19+yczM9PDwaP8laWlpTU1NAwYM6NWrl5qamri4uKampmCi7SIufdCh\noaEXLlxoamoiCCIrKysxMVHwYfHo2DGYNImfT3V/kYICaGtDSgoMHCi4SruLujrIzIS0NHjzBt68\n2XbvXv+5c0FGBqKiTExMMjIyBJOgR48ebWFhYW1tTebovZ0feUm7yb96thcvXowcOXLBggUA4O/v\nb2Vl1VaCfv/+fWFh4bBhw8TFxTU0NKDtrz60xSVBHz9+3Nvb++zZs/Pnz799+7bgY+JRejqEhEAb\nnYMUWroUzpyBH34QdL208uEDvHkD//wDSUmQng6ZmdDUBEpKoK//73+TJl1qbMyxtZ01a9anT59+\n//13gU3N7ODg4ODg0JUSyMm/yO/Rwp/8q8djMpnt9xpVVVUlJSWNHj26oaGBvNuhp6cnoOD4jUs7\nYzKZ5ubmly5dsrGx8ff3F3xMPNq4Efz8QPCDySdMAC8vaGwUQtVC0SIXZ2VBY+O/uXjwYJg0CfT1\nQVu79fDwPb6+69at8/f3b2xs3Lx5c58+fQQTr6OjYxdLoNfkXz3e8OHDt2zZcu7cuQEDBpw5c4az\n4kd9fX1CQsKwYcM+ffrUq1cvANDp/tOZcUnQUlJSISEhdXV1/v7+VVVVgo+JFxcvwuDBIJSPDZMJ\nNjbw119A5aRWfJeRkfH999/n5eWZmZl999133NeCqqiAtLT/z8XZ2dDQAJqaMGQIeVEM+vrQuzdI\nSnakRjk5ubNnz/L3LATD19c3Ojo6JiYmISFBWVnZ09Nz3Lhxwg6q5xIXF79+/frJkycvXbrk5OQ0\nderU1NTUXr16MZlMCQkJcXFxNREaVcUlQV++fDkvL8/S0vLYsWNdXGSTc9ebWkVF4OcHd+9SXlFb\nliyB77/vRgm6trZ24cKFAQEBQ4YMuXDhgoeHh//+/fDq1b+JOD0d3r8HKan/76Mgc7GODgj1qSph\nodHkXwgAAGRlZTdu3JibmyshIdHY2FhRUaGnpycpKTls2DBhh8ZnXBK0jIxMv379AMDb25u3QgV9\n1/u77+D770GI64EOGAC5uVBZyft64YKVmJhoamo6fOBA8PD4KiVl5N9/Q1oa6OtDv34wdiwsWgQG\nBkDv4Ueox6qsrCwrK+vTp09OTs6AAQPExMRGUzEhGj1wSdDNn0yRkZF58OBBZwsV6F3vu3ehrk74\nV6+zZ8Nvv4Gbm5DD6Bg2m/3+/XvYuxcGDqzZs2fdtGkx3ehucIeRi13V1taWlZWxWKza2tr+/fvH\nxsYKOy7Ei9ra2vz8fF1d3czMTCUlJQAgn7wXbVzm4oiJiYmJibl3796OHTssLS15KVRgS159/Ajb\nt4OfHyWFd0q3msK/X79+gz59SgwJ+R+b7ezism7dOmFHRInExMTExEQzM7OwsLCsrKxr1651u1FW\niCCI7OxsgiDevXtXXl4OAEOHDu3du7ew4xIQ7jcJyY2ZM2eePn2ah0IFd9d7927YsAHoMHONkhIo\nKsLbt9AtUkB9/a6yspjdu6urqw8dOjRQpAdx5+fnW1hYAIClpWUX76mgDkpNTb148SKLxVqyZAnP\nj/UXFBQoKyuXlZVlZ2fr6OiQ/a49DffJksiNwsJC3kZxCOiud1wcZGSAoJ4Y/rLFiyE4GLpFCvjp\nJ5gzx8LFxULYgQgAk8kMDAy0srKKioriPlgF8dXLly/XrFnj7e1dU1Mze/bskJCQTg2pJHOOuLj4\nkydPJk2apK6uTs511TNxSdCcrw99+vTZzNNqfh2/693Y2NhiVrDKykpylAxBEHV1dSwWi/vGp0+s\nb78lLl2qq61t8xgBb0ydShw4ULd1K0tKihbxtLXx5g0rMpL44w+a/OpYFE+ofeHChR9++CEiImLw\n4MECezqmJwsKCvrpp5/IW1n19fWXLl3aunXrF19VV1dXW1srLy8fHR1tYmKipqY2Y8YM6oOlOy59\n0K7/cXFx4dcT67q6ulz3p6SkbPvckydPjIyMmpqaqqqqyCfouW/cutU0b14Vm93eMQLeYDKrbG3T\nHj+mSzxcNz5+TPvrr6bDh6uqq2kRT1UV1U+sqKur//jjj5cvX/7hhx9UVVUprQsBQGNjI2fmE3IY\nXDsHEwRRWVkJAA8fPnz//j0A2NnZidJA5q4imhkyZMiQIUP69eunoqKipaWloqJiZmZG8ENaWloH\nj1y+fPlXX331hYPevCEmTSIaG7saFt89f04sXy7sINp15gzh5SXsID6zf//+P/74g7ry79+/r6+v\nb2ZmtmfPnitXrlBXEVcdas+i5eHDh9bW1s+ePXv06JGpqWlqairXw6qrqwmCiI+Pf/jwoWADpBZ/\n2/NnV9B8v+tdV1dHbvDz7jlBwLp14OcHTC6X/0I2YgSkpMCnT8KOow0FBRAUBM1m/e4Jur6iCuoU\nc3PzQ4cOnT9//urVq+fPn+//+eRltbW1BEFkZmaS426HDRtmbm4upEi7AS45jrzrLS4ubmlpmZeX\nx0OhSUlJ06dPV1FRMTAwUFRUnD59Oj9XBz97FkxN6Tt73MyZ8Pvvwg6iDZs2wcGDPW0RRSpWVEHt\nMzEx+emnn3x9fTmPpxEE0dTUVFFRcfny5bq6Oj09vS7OYNVDcEnQ5F3v5ORkf39/3u56r169esuW\nLXl5ee/evcvPz9+2bRvfRtoWFMDp07B9O39Ko8KiRXDpkrCD4CY8HOTkYOxYYcchaFSsqII6pamp\n6eLFi8XFxWw2e8mSJVTfFhYlXEZxdP2uN4XTM7q7g69vB2fnEY5evYAgICcHaDWWvrISfHzgjz+E\nHYcQBAYGnj17li8rqnzRp0+fHj161HxPbm5uT77l9fvvv/fv33/w4MGurq7CjqVb+ixv7t2718vL\nKyAgQF5enly55+jRo7t37+5soVQ9qPL778Bmg5kZH4qi1MKFcOkS8DRCkSpeXuDpCQoKwo5DCC5f\nvrxs2TIxMTEAOHfu3JIlS6irq6Cg4NmzZ833FBcX98B+lb///ruhoWHcuHEzZ84Udizd22cJmuww\nImcw6ApKHlSprIQffuge14COjjBxIo0S9MOHkJcH06cLOw7h2Lx582+//RYcHCwvL3/9+nVKE3Tf\nvn1bjPklRxNSVyOtpKSkJCcnOzk59YRZMgTjswQ9d+5cALCxsSkqKmKz2adPn269SmFH8HF6xvr6\n+qioqOrq6ikREdLbt3ePa0ApKRgyBJ4/h5ZmtB0AACAASURBVJEjhR0KQF0dfPcdhIQIOw6hsbS0\nXLJkiZ2d3eXLl4Udi2gqLi7+448/XFxcBg4cKNrTBggel5uEy5YtS0lJ2b17t4SERIs15AWstrZ2\n2rRpDx8+rLt79/G1a8n6+kIMpnOWLAGaPLT2ww/g6goaGsKOQ5hmzJhx9OjRWbNm8XM0UY9XV1d3\n7ty58vJyVVVVV1dXJg2HvXZ/XH6n9fX1Dg4OWVlZmzdvFu63s/DwcFtbW28vr7kvXmiEhvr6+jY2\nNt69e3c9/VfRHjcOnj2D+nohh5GUBLGxsGKFkMMQqvnz5wPAiBEjwsLCut59hwDgxo0b6enpkpKS\nCxcuVFRUFHY4oozL4AoxMTEvL68xY8ZER0dXVFQIPiaO4uJiHR0dKCuDbds0R4woPnKkoKDg6dOn\n9UJPfB0xcSJERoK9vdACIAjYtAkOHQIGQ2gxCBV505t8/IrcQ+GqET3Ao0ePJCQkRo8ePWXKFPJh\nblw8l2pcrqB9fHzU1NTc3d2Tk5ODgoIEHxPHxIkTg4KCCgDqJk3y9vaePn26lpaWcHtdOuGrr0C4\nS/AdPw4WFjBkiDBjECrOTe/mhB1U95OWlnbnzh0AGD58ODkFkkSPXPlMKLj8ARw0aNCgQYOqqqqE\nngoHDBjg6enp5uZWV1fn4OCwont9VdfVhYoKKC4GoUzQk5MDly/DX38JoWraIG9683ajGxUWFj56\n9MjR0bF37959+/YFABkZGWEH1eNwSdAPHjxYunSpurq6vb39oEGDZs+eLfiwOGxsbGxsbIQYQJfM\nnw+//gpr1wqh6rVrwdcXevY3UFzyigc1NTWRkZG2trYKCgq2trbQbAUPJHhcujhoPrmMoqJiYGCg\nsKPoGGdn+O03IdT722+grw+iu5JmB+GSVx1HEMSdO3fev38vLi4+fvx4aWlpKSmpHviIDd1wSdA4\nuQzfKCiAlhYkJwu00rIyOHwY9u0TaKU01vXJv0Tbq1evUlJSGAzGwIEDe/XqJS4ujgMz6INLgsbJ\nZfhp8WK4cEGgNW7dCjt2AHYX/qfrk39xEARRVlZGEAS/YhOirKws8ql0VVVVcrk/bW1tHMtMN1ze\nj8DAwBEjRhgZGWlpaZ08eVLwMYmUSZPg7l1od1EJfrpzBz59Ajs7AVXXHXh6er5582bLli3Z2dm8\nTf61Zs0aAIiPjx84cKC5ubmRkdHjx4/5HaaAlJeXP3z4EAAkJSXJmZo1NTVxtBxtffbGlJeXX7x4\nsV+/fitWrHj16lVqaurmzZt/+eUXnksnCKK8vFxRUZHRU4fiApMJVlZw9y5MmkR5XdXVsHMnhIdT\nXlG3sm7duufPn3dlZNiTJ08AYPv27aGhocbGxnl5eS4uLvfu3eNbiNSrra199uzZmDFjGhsb9fX1\nc3JyUlJS9PT0FBQUAKCxsfHBgwdXr1718/MTdqToMy3n4mCz2devXz9z5kxGRsaYMWPaWkuwfWvW\nrAkICIiPj583bx6TyWQymUFBQWb0n4WOIkuWwP79gkjQe/fC11+DigrlFXUro0ePtrCwsLa2JnP0\n3r17eStHWlqanJRRS0ur/XX26IMgiISEBD09PTExMVVVVXFxcRUVlStXrhw/ftzKyiogIMDMzGzL\nli3d6fmvHuazBJ2fnx8ZGVlbW9unT5+srCyeh9eIwBUHPw0aBO/eQWUlyMtTWMuLF5CQgPcGW3Nw\ncOji4h3p6elOTk7Z2dnBwcGurq4HDx6k/72ZzMxMFoulpqbGZDJlZWUlJCTk/2t+hw8fjoqKkpKS\nIgjCxsZm1apV5PNf27ZtE27MqLXPErSysjIAsFisQYMGdX3wY3e84qAKOd7OzY2q8hsbYdMmOHOG\nqvK7M0dHR852XFwcDyUUFRVlZmampqaSHxAJCYkTJ07wLT6+Kikp+fjxo66ubmlp6YABA8TFxVtM\nxd7U1CQpKUl+uhkMRt++fQsKChS6xSSRPRIlNwe64xUHtVxcwMWFwgTt5wf29sBTf5TICw0NvXDh\nQlNTE0EQWVlZnHk5Ok5MTMzAwIAzhtrd3Z3fMXZVTU1Nbm6ugYEBZwGXkW1MdctkMhUVFe/du2dt\nbf369eukpCT9bjRJZM/zWYJ+8uQJOZ1rdnY2Z17XlJSUzhbaja44BERZGWRkICMD+vblf+Hp6XDt\nGkRF8b9kkXD8+HFvb++zZ8/Onz//9u3bfClTV1c3KyuLL0V1RWNjY2Zmpr6+fl5eHvkltSNLFx0/\nfnzTpk3e3t4KCgpnz54l15rpTs9/9SSfJejs7Gy+FEr/Kw4hWLwYLl4ELy/+l7xxI/j5gZgY/0sW\nCUwm09zc/NKlSzY2Nv7+/nwpM0rYfw7z8vJUVVVLS0sLCgoMDAw6dRXcq1ev4OBg6mJDfPRZglal\nbFofmlxxCJO9PRw8CNu383nyz4sXYcgQ4MuSjyJKSkoqJCSkrq7O39+/qqqqK0XV1dVJSkoCQFuP\njCclJW3YsKH5nuTk5NH8e+a+tLSUxWIxmcykpKTx48draGho9OylGESegAaot3XF8c8//2zcuLH5\nnuTkZHJKQ1EjIQGjRkFsLJib863MwkI4dgzu3OFbgaLo8uXLeXl5lpaWx44d27VrFw8lJCUlbdmy\nJTY2VkZGprKy0sLC4vDhw1ynlh48ePCff/7ZfM+KFSu6vupFbW1tdXU1m81+/PixmZmZsrLyJAGM\n2kQ0QG2C/uIVx5AhQ6ho0DS1eDGcOsXPBL15M+zbBywW3woURTIyMuSjzN7e3ryVsHr16j179pia\nmrJYrJqamri4uHXr1kVGRvI1TC4aGxsrKysVFRWfPHmio6OjqKg4bdo0qitFtELJo/dJSUnTp09X\nUVExMDBQVFScPn06rgUHADBqFPzzD1RX86e027dBQgKsrPhTmigiZ+jv37+/qqqqtra2qqrq2LFj\neSiHyWRaWlqyWCwAkJKSIqde4newn/nw4QMAxMfHp6WlAYClpaWenh6lNSJ6oqSdCeuKoxuYPh2u\nX4e5c7tazseP4O0NN2/yIyaRRY6oW7Ro0ddff21mZhYbG3v8+HEeyjExMZk6derYsWP19PSysrIe\nP37ckcESPPj06ZOMjExmZuabN29sbW1NTEyoqAV1I5RcQQv+iqPbWLQILl7kQzm7dsHGjcBm86Eo\nUdf16UZ9fX09PT0lJCQSEhKYTKanp+eBAwf4GGFDQ0NjY2N5efnvv/9eX1+vp6dHzpSPECV5U2BX\nHN2PlhY0NEB+PnTl5ntsLGRmwo8/8i8sUUZON2plZRUVFcXbdKMMBmP8+PHjx48HAA8PDwsLC74E\nRj47AwBXr16dNGmSsrLyggUL+FIyEhmUXEFTfcXRvS1cCP/7H+8vr6+HrVsBZx3rsAsXLnRxutHm\ncnJyuh4S+VBJeHh4enq6mJjY3LlzyUe6EGqBkitoiq44RISTE9jawrff8vjyQ4dg/nzo4Y/Od4a6\nuvqP/Pu2MX369C6W8OjRIwAwNzdvPkkIQlxR3jXMlysOkSItDYaGEB8Pw4d3+rUpKXDnDuDt1s7w\n8/Nrfm+Qh7k4mnN1deXthUlJSenp6dOnTzfn4zhLJOooT9Bdv+IQQUuWwPnznU7QBAEbN4K/P+C6\nRJ0RGhp6//59FeFNk62mptbU1DR48ODBgwcLKwbUTVH+Uef5ikOUWVpCbCx0dn70oCAYOxYMDamJ\nSWTp6+sLt4e3qKgI1/pDvMHRb8LAYMDkyXD7NnT860VuLpw5A3fvUhmWaGpqaho1atSECRPIsZ4+\nPj7CjgihjsIELSRubrBtWycS9KZNcPgwSEpSGZNomtv1x4IQEhJM0ELSty+UlkJZGSgpffngsDBQ\nVQVTU+rDEkHN74LwtqIKQsKCCVp45s6FX3+FVau+cFh5ORw4AH/9JZCYRFDXV1RBSFjw3oXwzJsH\noaFfPmz7dti1i9oFZ0Xa8ePHt23bpqWltWnTJnt7e2GHg1AnYIIWHgUFUFeH9lcUi4mB0lLASSa7\ngFxRhcVi2djYkJPDIdRdYIIWqsWLoZ3Fh2pqYPNm+OknAQYkgvi4ogpCAoYJWqgmT4a//oK2Fijw\n8YEVK6BXL8HGJGouX748cuTI3bt3FxQU8LaiCkLCIqAVVRB3YmJgaQn378OECS3/6Z9/4OlTwITS\nZREREXPmzAEAb2/vixcv8jZnfwelpKQcOXKk+Z6YmJhBgwZRVyMSbZQk6I6v4YbAzQ18fVsm6MZG\n+OYbOHWKzyvM9jAhISEhISFPnz799ddfAYAgiNevXy9cuJC6GvX19bdu3dp8T0VFhZSUFHU1ItGG\nK6oI25AhkJ4OHz+CnNz/7zx2DCZPhr59hReWKLCxsTE2Nj5y5Ii7uzu5pxfF/UWSkpL6+vrN98jL\ny4vsGpuIepQkaHJFFQaDAbiiSkfMng3XrsGiRf/++O4dXLmCA5+7TlVVVVVVlZzK7sOHDwRBsHEN\nGtStUHKTkFxR5fvvvz937py3t/e0adNwRZX2uLh8NoX/unXw88+Af9K67NmzZ6NGjSotLb13756h\noeHw4cPDwsKEHRRCnUBJFvD19Y2Ojo6JiUlISFBWVvb09Bw3bhwVFYkINTWQkoLsbOjTB379FQYM\ngGHDhB2TKFi/fv3FixeVlZX37NkTERFhYGBgY2OD0+SjboTyFVVQRyQOH548c+bDIUP2JSTI/v23\nsMMREUwm09DQsLq6urS0dOTIkQAgISEh7KAQ6gQBjYPW1dUVTEXd0YMHDzxjY2cwGLuqqrYzmUmZ\nmcKOSER8/PixqakpMjJy9OjRAFBfX48JGnUvAurojIqK4ro/Nzc3+PNH6V6+fNm3B4xeiI6O9vf3\nl5WVdXFx+e233/YfOiR19KhUZaXdDz9cvXoVl97gizlz5owYMaKoqCgiIiI1NXXjxo04FwfqXgSU\noA0MDLjul5OTMzExab7n/fv3Ojo6AglKmF69enXy5MnCwsKjR49KSUlVVVXB/v0gJfXpjz9YLJaw\noxMRnp6eU6dOVVNT09HRSUxMXLRoEc4NjboXIQ8VYLPZkyZNar4nNze3J4wbXbNmTVxc3IoVK3x8\nfHR0dFasWPHtt9/W1dX9/PPP165dE3Z0IoLBYJBdzwBgZGRkZGQk3HgQ6ixKEvShQ4da7/Tw8KCi\nrm4qJCRk1qxZDx8+dHFxCQsLu3jx4q+//iopKRkeHq6uri7s6BBCtEBJglZSUlq9erWvr6+YmBgV\n5YsASUlJV1dXZWVl8sljPT29Fo8IIxoiCKK8vFxRUZGBj+AjgaAkQS9btuzSpUtr167FBwjb4ujo\niANyu4s1a9YEBATEx8fPmzePyWQymcygoCAzMzNhx4VEH1UJ9C98UhmJiidPngDA9u3bQ0NDjY2N\n8/LyXFxc7t27x0NR9fX1UVFR1dXV48aNU1NTaz6YZ+LEiXyOG3V/lI+Dxq5nJBqkpaXJGQu0tLQa\nGxt5KKG2ttbOzi46OjojI8Pe3j45OZkczOPp6Xnjxg1+x4tEAeUJOicnh+oqEKJUenq6k5NTRkYG\nOWb/4MGDWlpaPJQTHh4+ZcqUvXv3uru7BwUF+fr6rlmz5vXr187OznZ2dvyOGokCyvuImy96j1B3\nVFRUlJmZmZqaqqysDAASEhInTpzgrRzOGH8dHZ2ioqLmg3mmTJnCz6CRSKA8Qbu6ulJdBUKUEhMT\nMzAw4DxsxZldurMmTpy4du1aOzs7BQWFffv2TZ8+vcVgHoRawFEWCHWarq5uVlZWZ19laGjo6enp\n4uJSV1fn4OCwcuVKBoOBg3lQOzBBI9Rpbc0tU1FR0WLloIyMDA0NDc6PNjY2NjY21AaHRAgdE3Rc\nXJxc8/WfOiY9PT0rK6sr67+VlJSoqKjw/PLi4mJVVVWeX15VVSUmJtat4x8yZIimpmZnX5iYmMh5\nILu7aGtumfLy8vT09OZ7yIlWQkNDO1X+P//8U1paSuljBF18u7+o6+35i6g+hZqaGl1d3RZrmH0R\nf9szgyAIfpXFFzk5OTdv3uThhdevX//w4UOfPn14q7e6uvrZs2cWFha8vRwA/vzzT1tbW55f/vr1\na2lp6a7EHxcXZ2lpyXMAXY/f0NCQhwAYDMb8+fPl5eV5rprOeGvPwcHBMjIyampqVIRE6uLb/UVd\nbM8dQfUpZGdnKygozJgxo1Ov4m97pl2C5pmfn5+Wltbs2bN5e3l+fr67u/ulS5d4DmDChAl3797l\n+eVHjx7V0NBwdnbm7eUFBQUbNmy4fPkyzwF0MX5/f391dfU5c+bwXAJtCX5ume3bt9vb25ubm1NX\nRRff7i/qYnvuCKpP4erVqzk5Od988w11VXyRgCbsR6j7UlJS8vT0lJCQkGpG2EGhHoGOfdAI0QrO\nLYOEBRscQl+Gc8sgocAuDoQ6AeeWQYIkOjcJs7OzpaWleb7x3djY+OrVq+HDh/McwLNnz1os39Up\n3T3+d+/esVgskV9tYP78+V25E9tBqampmpqalI5s6eLb/UVdbM8dQfUpFBUVVVdXUzoQ5YtEJ0Ej\nJADBwcE4ewESGEzQCCFEU9gHjRBCNIUJGiGEaAoTNEII0RQmaIQQoilM0AghRFOYoBFCiKa6ZYKu\nrq6ePXv2xIkTR40a9fjx4+b/1NjYqKysPGrUqFGjRq1Zs4bSMOrq6oyNjTMzMwUfQDu1CCaAY8eO\nTZgwYdiwYbGxsYKvXVQFBQW1njmPIAgPD48JEyZYWVklJSXxVnI7hXT9LWuncL4ET3X8zVH3FvCI\n6IZOnTq1ZcsWgiAePHgwYcKE5v+Ulpbm6OgomDD27NnDYrEyMjIEH0A7tQgggMTERBMTk7q6upcv\nX5qamgq4dlE1adIkSUlJX1/fFvvv378/efLkpqam+/fv29vb81Z4O4V0/S1rp3C+BN9+OXxscpS+\nBbzpllfQZmZmq1evBgCCIBQVFZv/U2pqanZ2trOz8/Tp0+Pj46mLISUl5cmTJ6ampi32CyaAdmoR\nQAARERFz5syRkJAYOnRoeHi4gGsXVbdu3Tp48GDr/TExMebm5gwGY8yYMS2+L3ZcO4V0/S1rp3C+\nBE91/ByUvgW86ZYJesiQIXp6eitWrJg8eXKLJZaVlZU3btwYGhq6b9++efPmEdQ8J9nU1LRhw4af\nfvqJwWC0+CfBBNBOLQIIoKCgIDk5ecqUKVZWVs+fP+9gYKh94uLiYmJirfcXFxfr6ekBgJSUlLS0\ndF1dHQ+Ft1NI19+ydgrnS/BUx89B6VvAY0gCq4mPPnz4IC0tffLkSU9Pz4kTJ759+5aTKE1NTcmr\n2mHDhjGZzPLyciUlJb4HcPz4cVtbW66LlQkmgHZqEUAA8vLyubm5169fLygoGDlyZE5ODovFEljt\noiQ4OPjGjRu9e/fmumgLSVlZmVxBvKamprKyUlJSkofy2ymk629ZO4V3JfgOViGAJsevs+BBt7yC\n3r9///HjxwGAxWK1+Gt25MgR8ktKXl4eg8Fo0QHCL0+ePPnjjz/s7Oxevny5ZMmS5vcJBRNAO7UI\nIABLS0sFBQUJCQl5eXkmk9n8mkUwpy8yXF1dL1++3E52BgArKyvya3VcXNy4ceN4K7+dQrr+lrVT\neFeC72A5Amhy/DoLXgiyw5tfCgsLp06dam5uPnz48IiICIIg0tLS+vTpQxBEeXn53LlzLSwsxo8f\n/+jRI6ojsbKyIm8SCjiA1rUIMoCmpqatW7eOGTNm2LBhISEhAq5dhPn5+XHuUHF+pU1NTZs3b7az\ns5s0aVJKSgpvJbcuhI9vWTuF8yV4quNvjrq3gDc4mx1CCNFUt+ziQAihngATNEII0RQmaIQQoilM\n0AghRFOYoBFCiKYwQSOEEE1hgkYIIZrCBI0QQjSFCRohhGgKEzRCCNEUJmiEEKIpTNAIIURTmKAR\nQoimMEEjhBBNYYJGCCGawgTdaTk5ObKystb/Wb9+fUhIyN69e4uLi0NCQgCAs9E+IyOjTtV75coV\nLy8vHoNGou7rr7+2trYeMGCArq6utbX1rFmzbty4wd8Gw2nYZIP/4vGtW2xSUtKMGTPGjx9vaWm5\ncuXKkpISfgUjqp+ObrkmodD17dv33r17LXampKSEhITMnTuXbDpz584VRmiohyIXgQsMDMzJySGz\n540bN/hbBadh89a2CwsLZ86cefnyZRMTE4Igzp496+joeP/+fSaTl8vEFsFcuXKFh0LoD6+g+eDK\nlSu7d+8+cuTI06dPQ0NDORv19fXLli0zNzcfP378gwcPAKCoqGjq1KmTJ09etGhRbW1t80KmTJkS\nFxcHALdv3164cOGnT58cHBwmT548efLk5o0vODjYx8cHAGpqakaNGgUArWtJTk6eOXPm9OnTZ82a\nVVpaKshfBaKVZ8+eubq6jh8/3s/PD7g1laqqqvnz50+dOnX8+PFRUVEAEBYWtmzZsqFDh758+bLF\nwZyGTTb4ysrKefPm2dnZjR8//p9//mmrxXKcP3/e1dXVxMQEABgMxtKlS8XFxZ8+fdq6Sbcu6saN\nG4sXL166dKmlpSV5Li2CIasQvc8CXkHzIiMjw9ramtwmF6wEAHd395KSkjlz5gwdOpTcOHnypKSk\n5MOHDwsLC83Nzd++fevj4+Pg4LB69eqYmJhbt241L3Px4sW//vrrqFGjzp8///XXX+fl5bm6us6d\nO/fRo0c+Pj7Ozs5tBXP27NkWtdy+fXvo0KF79uwJDw8vKipSVlam7DeBaK2kpOTmzZtFRUU2Njbr\n169v3VQCAwONjIy8vLyys7PHjx9PLn+cmJj4/Pnz1gdzWjiZNH/55RcjI6MdO3ZER0fHxMSwWKz2\nW2xGRsb48eOb7xk4cGB6enrrsLk2/sTExGfPnpWWllpZWa1fv75FMCTR+yxgguZFiy6O7Oxsroe9\nevXq9evXS5cuBQB1dfWqqqrU1NTly5cDgJmZGYvFan6wo6PjgQMHduzYkZKSYmlpWVBQcOfOnQcP\nHtTV1TU1NbUuvKGhoa1a3Nzcfvzxx2nTpunq6pIr0qOeacKECQwGQ11dnfyxdVNJSUlxcXEBgD59\n+gDAp0+fAMDa2lpCQqL1wS0KT0xMXLVqFQBYWlpaWlrm5+e332L19fVfv35Nbjc2NoqJiSUnJ3/1\n1VecnZwmLScn17oo8lxUVFTaOV/R+yxgFwc/cVbgJTcMDQ3Nzc3Pnj37888/29vby8rKDho06O7d\nuwAQGxvbootDVlZ2zJgxHh4erq6uDAYjICDA2Nj46NGj8+fPb76wL5PJLCwsBIA7d+6Qe1rXEhIS\nMnny5Fu3bunr658+fVow545oSEJCovmPrZvKwIEDHz58CACZmZkEQcjIyHBe1fpgaNbCyQNiYmIA\n4O7du9u2bWurxXIsWrQoODj42bNnAHDhwgVHR0cxMTETE5PWTZprUS3OpUUwbZ1gd/8s4BU03/Tq\n1ev169cXL16cNm0aubF8+fLly5dPnDixqqpq8+bNALBly5ZFixZduXJFU1Ozf//+LUpYvHixnZ1d\nbm4uANjZ2Xl7e0dFRfXr1y8rK4vsngYAa2vroKCgmTNnGhkZkR+n1rWMHDly6dKlioqKTCbzxIkT\nAv0tIBpr3VRWrVq1dOlSW1vbmpqas2fPtn8wp4WTX/7c3d2XLl06ceLE6urqoKCgsrIyri2WQ11d\nPSwsbMuWLaWlpZWVlYaGhgBQWlraukm31fibaxFMWzF3988Cg+vfOoQQotq7d++0tLTExMSEHQh9\nYYJGCCGawj5ohBCiKUzQCCFEU5igEUKIpjBBI4QQTWGCRgghmsIEjRBCNIUJGiGEaAoTNEII0RQm\naIQQoilM0AghRFOYoBFCiKYwQSOEEE1hgkYIIZrCBI0QQjSFCRohhGgKEzRCCNEUJmiEEKIpTNAI\nIURTmKARQoimMEEjhBBNYYJGCCGawgSNEEI0hQkaIYRoChM0QgjRFCZohBCiKUzQCCFEU5igEUKI\npjBBI4QQTWGCRgghmsIEjRBCNIUJGiGEaAoTNEII0RQmaIQQoilM0AghRFOYoBFCiKYwQSOEEE1h\ngkYIIZrCBI0QQjSFCRohhGgKEzRCCNEUJmiEEKIpTNAIIURTmKARQoimMEEjhBBNYYJGCCGawgSN\nEEI0hQkaIYRoChM0QgjRFCZohBCiKUzQXLx8+XLy5MlsNltFRcXBweHNmzcdf218fLyRkVFHjvz4\n8SODwSguLuY1TC7ExcUbGhri4uJGjRrFx2KREPXu3ZvxH3l5eXt7+7y8vI6/nGtjKC8vV1RU7FQY\nHW/YHURF+xc9mKBbamxstLe3HzNmzIsXLxITEw0NDZ2cnAiCEHZcndC3b19vb29hR4H45ubNm2Vl\nZaWlpc+fP//w4cP27ds7/lpsDN0aJuiWcnNzc3Jytm7dqq+vr6mp6ePj07dv34qKCgAICQnp37+/\niorK6tWra2trAeDkyZN9+/aVlpY2MzN7/fp1i6Kio6NHjBghKytrZ2eXm5vbwQCuXr1qaGjIZrNn\nz55dVFRE7uxI1ZMnT25sbDQwMEhMTNy5c2dbpaWkpFhYWBw6dEhbW7tv375RUVF8+K0hKsnLyysq\nKiopKfXv39/V1TU9PZ3c37qBNTQ0rF69WklJSVVVdc+ePQCQkZHBaQx+fn46Ojo6OjqnT58m9zx+\n/NjMzKz1djsNu3UVHHZ2didOnCC3fX1958+f387BbWl9Uq2L5XpYYmKitbX13r17jY2N2zqFs2fP\n6unp6enpnTt3Tk9Pr60aaYRAn6uvrzcyMpoyZcqtW7c+ffrE2f/69WsVFZVHjx69ffvWxMTk1KlT\n2dnZkpKS9+/fLyoqcnNzW7lyJUEQL168GDJkCEEQxcXFKioq4eHhpaWla9eutba2blFRZWUlABQV\nFTXfmZ6ezmazIyMjS0pK3Nzc5s2b1/GqCYIQExOrr69/+vSpiYlJW6UlJyfLysoeOHCgqqpqy5Yt\nY8eOpfC3ibpMW1s7Ojqa3M7NzXV0nuFSEAAAIABJREFUdNy3bx/RRgMLCQkxNDTMyMh4/vw5i8VK\nS0vjNIYHDx4oKSndv38/JyfHxsaGzWYTBBEbG2tqakoWztluv2G3roIT6rFjxxwcHMhtc3Pza9eu\ntXMw1/bP9aRaF8v1sFevXrHZbDc3t8TERK6n8PLlS1VV1b///js3N9fCwkJXV7etGukDEzQXNTU1\nAQEBdnZ2KioqU6ZMefLkCUEQe/bsWb9+PXlAfHz8/fv3q6urs7KyCIL4+PGjh4cHmf447fjs2bOz\nZ88mj6+urpaRkWloaGheC9cGeuTIkSVLlpDbhYWFEhISDQ0NHayaaJWguZaWnJysoKBQX19PEMSr\nV68MDQ35+stDfKatrS0rK8tmsxUUFADAzMyMbEhcG1hISIi+vv7ff//d1NRUVFRUW1vLaQwbNmzY\ntm0befzDhw/bSdDtN+zWVXBCzcvLk5OTq66ufv/+vaKiYnV1dTsHc23/XE+qdbFcD3v16pWkpGRN\nTU1bp+Dp6bl582byVb///juZoL/4ORUucaFevtNRXV0dQRCrV68mOxMuXbpkaWkZHR2dk5PTv39/\n8phhw4YBQENDw6lTp27dusVms1kslry8fPNy3r17FxkZyfkaJSkpWVhYeOXKlR07dgDA4cOH586d\n27r2/Px8zkvU1NQkJSWLioo6W3X7pQGAhoaGuLg4AJD/RzR35syZ0aNHA0BxcfGCBQuCg4OXLFnC\ntYHNmjXrw4cPK1euLCgoWLt27aZNmziF5OfnT5o0idzW19dvXQvx340WcXHxdlpXO1VoamoaGRnd\nu3cvOzt75syZUlJS7RzMFdeTal0s18MAQEdHh8VitXUKubm55K8RAPr06dN+je3HKTD4+WwpNDT0\n9OnTd+7cAQAWi+Xm5nbu3LkXL1706tUrJyeHPCY2NjYtLU1cXDwiIuLPP/9UVlYODg6+ceNG83I0\nNTVtbW1/++03AGhsbHzx4oWGhsb69evXr19PHvDx48fWtWtoaCQkJJDb5BWHqqpqZ6tuv7Ty8nIG\ng9HF3xISJE1NTTKD6OnpzZ49+8WLF0uWLOHawNLS0mxsbJYtW5adnT1nzhxFRUVOt7KWltbbt2/J\n7YyMDE7hDQ0N5AanjYWGhrbTutLT01tUsW7dOs6/Ojk5RUREvH379ptvvvniwVzPtPVJtS6W62El\nJSWcCw6up6CpqZmdnU0e8O7du/ZrpAm8SdiSra3t8+fPd+3a9ebNm9evXx87duzZs2cTJkyYPXv2\nhQsX/v777/T09I0bNxYXF5eUlMjJyUlLSxcWFvr5+VVXVzcvx97ePjo6+ubNm8XFxdu2bdu4cSPX\ntFhRUVH+n9ra2hkzZly9evXOnTtlZWUeHh5OTk7i4uKdqrp53udaGnW/OiQAGhoaZHLh2sDCw8MX\nLFhQUFDQ2NhYW1srLS3NeaGzs/Px48ejo6Pz8vJ27txJtkY2m52QkBAfH19SUuLv708e2X7DbqcK\nAHB0dLx27VpCQsLEiRO/eDC0av9tfWpaFPvFDxfXU3B2dj5z5kxcXNz79+9//PFH8sgOfk6FRth9\nLHSUmppqb2/fq1cveXn5sWPH3rp1i9x/7ty5vn37KigoLF26tLa2try83NbWVllZ2dzc/Pr16716\n9Tp//jynq44giNu3bw8ZMkRGRmbChAlv375tUQvZB9fcyZMnCYK4cuXKgAED5OXlHR0dCwoKOl41\nQRBz586Vl5e/d+8e2e3ItbTk5GROv3PzbURPzW8SEgQRERGhrq5eUVFBcGtglZWVjo6OsrKyysrK\nq1atqqur4/RBEwTh5+fXu3dvbW3t06dPa2trEwTR1NS0fv16OTm5oUOHhoaGkn3Q7Tfs1lW0CHjw\n4MFff/01ud3OwW21/7Y+Nc2L5XpY8xsqbX1AAgMDNTU1BwwYEBgY2MHPqXAxiG41whchhHiTkpJS\nUFBgZWUFAJGRkQcOHLh7966wg/oCqro4Ghoa6uvrOT+SQ3cRQkhYysrKFixYUFhYWF1d7efnN23a\nNGFH9GWUJOirV69qa2sPGDDg+PHj5B5TU1MqKkIIoQ4aO3bsN998M2LEiP79+2tqaq5du1bYEX0Z\nJV0cxsbGUVFRMjIyU6ZMOXDggIWFxfDhw+Pj4/leEUIIiTBK7ulraGioqqoCQGBgoJubW0xMTMdf\nW1hYGBYWRkVUiJ6YTObcuXPJpzBED7bnnoa/7ZmSBK2np/fVV199++23RkZGM2bMmD17dmlpaQdf\ne+vWrYSEBGtrayoCQzQUFhamo6MzZcoUYQfyZQRBkPPAdXwkVuv2LC4uLikp+enTJ0pCFBR2dvaw\nc+fejxqVNmUKwRTcaN0BERFy+fnJjo7VKioCq7RT+NueKUnQR48evXv3bmNjIwDs2LEjLCzs9u3b\nHX/5qFGj5syZQ0Vg9BEdHe3v7y8rK+vi4kIO7eyx0tLShB3CF6xZsyYgICA+Pn7evHlMJpPJZAYF\nBXEeAPmiFu2ZHBEsIyNDTbAC5O7e68SJ4f37w8iRgqt0zhy4d8/gwAGwtQUPD8HV22H8bc+UJGhJ\nSUnOHxAGg+Hk5OTk5NTWwWVlZc1/rKqqaj2aXfS8evXq5MmThYWFR48enThx4v/+97+LFy9KSEis\nXr26W1xL9ihPnjwBgO3bt4eGhhobG+fl5bm4uNy7d4+30sTExKSkpPgZn7AwmbBq1b/b16+DpSV0\nco5pHllbg7U1kI/PVFVBSQn899y26BHQc2W6urpZWVmt9yclJe3evbv5nufPn48YMWLp0qWCCUxY\n1qxZExcXt2LFCh8fn4iIiFu3bv3vf/+rq6tbunSpqqqqiYmJsANELUlLS5PzWGppaZHfDnlTWVlZ\nUlLCmfxBRMjJgaMjzJsHK1aAYJ5WJS/jKith7VpQV4ft24HbBCPdnYASdFuTDg8ePDgkJKT5nhUr\nVjQ1NQkkqDZFRUU9efLE2Nh46tSpDAYjPj7ex8dHQUHByspq4cKFXS39wwcoLw8JCZk1evRDLy+X\n7dvnSUkdHj2afeAAuLq6u7vfuHEDEzStpKenOzk5ZWdnBwcHu7q6Hjx4UEtLi+fSZGRkxMTE+Bge\nLUyYABYW8PPPEBYGzs6Cq1dDA65fhydPYNcuCAwEWVnBVS0QAkrQBgYGgqmo677//vvs7GxnZ+fI\nyMiIiAh/f/+4uLhffvlFRUVl1qxZLRN0fT2Ul0NFBVRUQFnZ/2833+BMjsFgAEEAmw1stmRxseu5\nc8ps9sIxY2ry83N1dNQGDwYvr/fOzmw2W/BnjdpRVFSUmZmZmpqqrKwMABISEpz543kgOl0cLUhI\n/H+ncGwsKCjAkCECqnrMGLhw4d/t/fthxAiwswNaTanBK5w65zNNTU1//PHHo0ePGImJU21tD2zd\nWnPkyPLa2uJ9+zzu3l2lpARTp0KzJyRBQgIUFYHNBkXFf/9js0FLi8zC//4oJ9e6IkcAx/+2s7Ky\n5syZs3TpUpusrMhDh37CJU5oRkxMzMDAgHOR4e7u3pXSRLOLowU1Ndi6FaSkYM8eQfc8LFkCAQHg\n4wO7doGNjUCrpgAlCfrQoUOtd3rQ8pZrC42NjbKysoxnz+Dbb8HJic1i1bLZqY2NZ1+92nrihEa/\nfqCoyPcuNl1d3du3b9+4cSNz3rzTL16IKyvzt3zEd+3cU9mwYUPzPW/evLG2tm5+T0U0uzha6NcP\nfvsN4uLgxg345huBVq2tDfv2wcePQE7GVFkJ1dWgri7QGPiHCQA//vhjSUkJHwtVUlLy9PSUkJCQ\naoaP5VNHQkJCT1m5aP784mPHwvr2vaGuzv7qq58fP86urfU4enSXnx9FN0CUlJQWLVo0Zds28dpa\nSE2logrER+3cU/nzc7a2ti3Ssch2cbQ2atS/2TkvD9avh/fvBVe1nByQk+7n5MCsWbB4MTx7Jrja\n+UccAJKSkioqKq5evcqZOoODxWLduXOns+1p2bJlly5dWrt2bXecfTiQxbpmZnbF29vQ0PDSpUsA\nEBQUJLjqv/sODh6EU6cEVyPqvK7cU+kRXRwtaGnBrFmwZAkYGMDhwyDIcbSDBkFMDDx6BCdPgrEx\nSEgIrmp++P8EumLFihUrVgBAaWlpUVERm80+ffq0s7Mzb3/t//rrL77FKEiXL4tLS885cUJoz8mY\nmoKXF2Rni/DQzh6uJ3Rx5OXlhYSEsFis+fPnKykpxcfH+xw/rqCnZ6Wjs5B87LCiAgR5M9zcHMzN\n/91euxaUlGDlym7xEeNyhbts2TI3N7dbt24ZGBisXLmS5wH53U9WFvz0E9y5I+QwPDzgyBE4ckTI\nYaD/8Peeish3caSmpi5atMjd3b2mpsbOzi48PPyzcVAsFgDAgQMQHw8eHvDfMomCQ44FXLYMFiyA\nr74SdO2dxOUh+vr6egcHh6ysrM2bNwt9SLLgNDbCsmUQECD8oZRTpsDz51BUJOQwREVpaenr16/z\n8/P379+fylP/Pn/vqVRWVnJWxhNJJ06cOHz48Pz5893c3DZt2nThwoXly5czmUwPD49VnCcPfXwg\nIABu3QJB9h+SxMXB2Rn+/PPf7Pz2LRw5AsXFgg6jY7hcQYuJiXl5eY0ZMyY6OrqiokLwMQmHry9M\nnCjQWQXasX49HD0K338v7DhEQde/EfL3norId3FUVVUpKSmR20pKSq9fv3758uXZs2e3bt362Xqs\n+vrw38KA8M8/8NtvsHIlCH7BVl1d0NGBpUuBxYLTp4FmsypyuYL28fFRU1Nzd3dPTk7+v/bOPJDK\n7I3jX8reYkyrdtoHTSWVbFEoRotUUxKppGG0ykiTUCkt0zalpsYkJVpERkqbhEH9IiXaUJmKVglX\n3N8f19wuLnGX933v7Xz+up1773O+9L6P8zzvOc9D6fMxGklLw+XLWLWKbh3/MXUqLl5Eg6ZtBAEQ\nSUSYkJAgqifeUp/isLOzW716NedoT2Bg4JQpU3bs2FFYWLhixYq1a9fy/86gQRg2DIsWYeZMUFzn\nj7OgjonB1q3gVLC6dg1FRZRqaJw615yPjw/3NSfvVlhYqKurS7UoivnwAe7uOHkSFFZN/AKyspg3\nD8HBzKzXJVkwLSKU+l0cpqamLBbL09NTXl4+ICBAS0vry+s8WVlYW8PaGsXFtXs8Ll3C99+DyjMB\nvXp9FuPhgepqLF5MQ4q8LnUc9Pfff0+XDjpZsQKenujenW4ddZkzBwYGcHODVK+2KCAwMDA+Pt7R\n0TEiIoIJEaHUpzgAWFpaWlpaCvLNjh1rX3z6hLlz0a4d3NwwerQItX2ZMWMwZgyKi5GbCwBsNrKz\noa1NqYb/qOOgpzWocpKRkUGhGDqIjER5OWxt6dbRADk5zJqFv/6CiwvdUiQVZkaEUp/iEA3m5jA3\nR34+8vMB4NMnPH0KKsOOjh1r/1pUVeHAAaSmwtoabm6ULur5PiSMjIwMDQ2tqalhs9kFBQXZ2dlU\nCqKUoiJs347z5+nW0QgLFsDUFM7OFNVvlDqYGRFKfYpDlPTuXeuUP36ElxfevoWdHWbPpjSslJfH\nzp2orER0NF68gJoaXr9G69bUPE7kc+cHBwf7+fmFhITMnDmzRZ1QJAw2GwsWYMcOvsWMGIGyMqys\nEBmJH3+kW4pEwsyI8GtIcYiedu0QHo4PHxAVhX//RZ8+KCmBigp1hxIVFMBti/PgAby80LEj5s2D\nmNtr8HksJisrq6+vr6CgYGpqyvx2RIKzYwdGjMCIEXTraJKff8bevRBD5/Wvh8jISBsbG2traysr\nK0dHR7rlkBSHELRpA3t79OkDAFlZMDPDzJmIjaVahp4eLl1CUBAePgRnX9CzZ/xv0upqXL4Md3eB\np+LjoBUVFSMiIlgs1p49e8rKygQ2zWgyM3HmDHhylAylfXuMGkXDJShFBAcHe3l5qaurL1++3MrK\nim450n9QhSJMTZGcDF9fcH+Zd++CxaJOQM+eWLy4dutXRARGjcKyZbhzp85nXrxAenqdAsUthI+D\nDg8PHzZsmK+v74sXLxrdtyjRlJdj0SIcOiQZud3lyz/v5ye0HKZFhMrKyp0ktvol4xg4EK6uta+v\nXoWJCRwckJZGtYylS5GcjB9+wMOHiYmJWVlZtePq6li4UBjDfDyUsrJy3759Afj5+QljujlkZWXN\nq3scvqCgQJ9b1kRMeHlh0aLaQIn5dO6MgQORmAgjI7qlSCRMiwhJikNcuLrC1RW3b9f2ky0vR0wM\nzMzw7bdUzN6qFcaOBVCzZ0+VEEvmevBx0Lz7kJSVlRMTE0U1WUN0dHTqPbcRe0/Cs2dRXIy5c8U4\nhcjx9MTPPxMHLRjh4eFFRUWGhoZ79+5lQkRIdnGIF+6GZTk5FBVh1izIyiIoCFpatMoSED4OOikp\nCcCnT58uXryYRn2wIFaeP4efH3P31TVGnz5o3x43boA0k205VEaEAN68ecP7z8rKSrm6NYjJLg6K\naN0aS5ZgyRKUloLzC09JQVwcJk6Enh51x4ZVVbFvn8Df5uOgufHXpEmTDh06JLBpxsFmY/58bN0K\nVVW6pbQcb2/4+SE8nG4dkgeVEWFWVtby5ct5R3Jycur1aCcpDqpp27b2hZ4ePnxAWBi8vBAdzbS6\nSHzh46ADAwM5L16+fMmEnJ3I2LcP2towNKRbh0AMHgxON6z+/emWImFQGRHq6OhcuHCBd6Rhyo6k\nOGijVSuMH4/x4z+P7NmDkydhYgIbGzDyWBOfdX73/9DV1Q3lNjOXdG7fxvHj8PenW4cQeHlh82a6\nRUgenPLNbdq0mTRpEhOOxZJdHAzip58QE4Phw8FtMnnlCnNK2YHvCtre3p56HeKFxcJPP0nMvrrG\nIN2wBIJpESFJcTALFRVYWYG7Qf75c+zYgeJiODuDpxc7XdRxWFpaWgAqKyvfvHmjoKBQWVnZr1+/\nlJQUmrSJDh8f2Nujb1+6dQgN6YbVcrr/V6ewZ8+eK1eupFcMSIqD4cyciZkz8ekTOJVpWSzY2UFP\nD2PHQk+P+hVenfk4AeCcOXNcXFxGjRqVkpLSsM+35HHxIgoLpSQ5YGGBDRtQXPy5KiPhSzAtIiS7\nOCSA1q1rd0/Ly+PwYVy7hhMn8P49LC3x4QMePMCQIZCRoUAInxz08+fPDQwMWrdubWhoWMSkdIwg\nvHoFLy/s3k23DtHB6YZFaAZaWlpaWlr9+vXr0KFDt27dOnToMJriysL8ICkOCaN9e1hbY9s2cCpc\nf/yIkBBMmICpU/Hw4Re++/Ah7OwwfTrS02tHLl6EsTHs7dG8x3v8iyXt27cvJydnz549SpQVixIT\nLi7YsAEdOtCtQ3SQbljNJjs7Ozs7e9SoUVFRUQUFBadPn9bU1KRbFKnFIeF06oTffsO5czh4sLbL\nx8WLsLHBtm343//qf3jfPvj44PffP6cls7JqN2I1r/0pHwcdGhp6//59T0/PwsLCv/76S+AfhH4O\nHULPnnV21UgB3G5YhObBtIiQ7OKQEr75BgoKAGBmhoMH0bMnQkPx4cO7d+/ev39f+5lnz6CpiQ4d\nwH06bW6OLVuwfTt++aU5k9Rx0AEBAQB+//33tm3bDh8+XElJabfkRtN5eTh0CBs30q1DDMyZg8hI\nVFTQrUMyYFpE+LWkOBITsXEjvpJYoWNHTJuGbdvQppbacXV1PH6MV6+golI7cucO5OTQrl0zU9h1\nHhL2798f/+3lEBXV1dUsFovqG6OqCvPnY//+2j9xUgbphtUSQkNDN23aFBsbO3jwYCZEhF/LLo6d\nO2FmBldXVFVh9mzY2jK3M4ZIadWqlSz3HPmiRVi7FjIyWLYMCQl49w5du8LJCVVVWLasOdbqrKCn\nT58OwNTUVFtb28DAIC8vT0dHRwCJ//zzj6Gh4cSJE1NTU/v169ezZ8+IiAgB7AiOnx/s7DB4MKWT\nUsmCBfjzT3z6RLcORsPMiPCrSHGUleHNG7i6IjYWERGorMSkSbCxQWSkMMWRJY++fREWhiNHoKuL\nceNgawtDQxw/jlOnYGzcHAN8tvU5Ozs7OjrGxcVpamouXLjwypUrLVXl5ua2d+/ely9fGhoa5ubm\nqqqqjh07luP9qeDKFWRlISqKouloQVkZ1takG1bTiDwiZLPZb9++VVVVlRFii9VXkeKIjwe3q7eq\nKhYuxMKFyM3FsWPYtAnDh9c2rSd8CT4PCauqqmxsbAoKClauXClY5c9WrVrp6uqOGzeue/fuvXr1\nUlNTo+6KfPMGK1fiwAFqdinSibs76YbVNKKKCBcvXgzg1q1bAwcO1NfX19LSSk1NFVjVV7GL4/Rp\nTJlSf3DAAPj6Ii0Nc+YgNBSjRsHXFwUFdOiTGPisoFu1auXj46Onp3ft2rV3nOM0LURLS2vy5Mnv\n37/v27fvypUr1dTUqIvpfv4Z/v6Q+hASPN2wrK3plsJohI8IOSWWVq9eHRkZqaOjU1RUNGvWLAHs\ncJD+gyosFh4/bvTgrqwsDAxgYICKCsTEwMMDlZX48UfY2n5+kkb4Dz4r6MDAwI4dOy5dujQnJ+fg\nwYMCGA0ODl6wYIGvr298fHzPnj3Ly8spejgTFoZ27T7HVlIP6YbVDISPCDkoKSlxFuDq6urV1dUC\n25H+FEdiIqe3yBdQVISdHaKicPAg3rzBxIlwcEBCAgkKeeGzgh40aNCgQYPKysoWCtpNq1WrVtzu\nnEuWLBFcXYvIz8e+fahb7FHKId2wmoHwEeGjR4+mTJlSWFh45MgRe3v7zZs3q6urC6xH+ndxnD4N\nZ+cWfF5dHR4e8PDAnTsIDYWXF8aNg7Mz+vUTm0SJgY+DTkxMdHJy6tSpk5WV1aBBg2xtbYWfplev\nXgX8kk1FRUX1KppmZmZqaGi0eILqasyfj717Id1rk4aQblhfIjAwMD4+3tHRMSIiQrCIsLi4OD8/\nPy8vT01NDYCcnNz+/fsF1iPlKQ42G5mZzTwmV5/vvkNgIKqrcfky/P2Rnw9bW8yeLVUngVsIHwe9\nevXq5OTkVatWubu7W1paisRBX+KWW62LkpJSvX4TV69eFeTy3bgRlpYS2nZMKEg3rC8hkohQU1OT\ne0x86dKljX2SxWLdvn2bd6SkpKR9+/b1rElziiM9XdhLsVUrjBuHcePw9i2io+HoCFlZzJmDyZNR\nt3nY1wAfB11dXc15pte+fXsV4dL2LBZLXl4eQGM1EL755ptx48bxjhw/frzFicLUVCQl4e+/BZUp\n4ZBuWE1CZUT45MmTyMhI3pHHjx/3rfu4TMpTHHz3bwiGqiocHODgULs/LygIQ4di4cKvai3Cx0EP\nHDhw8eLFhYWFPj4+guXa7t696+npmZKSoqysXFpaamBgsG3btv5i6tVUWoqff8bp09R1gWQapBtW\nk1AZEWpqanL7A3B49epVvQWHlKc4EhMRECBim5z9eb/+iuRk7N+PrCxMnoy5c9Gli4gn+hIXL170\n8/Pr0aOHhYXFnDlzkpKS4uLiOnfuPG/evDbiOSfJx6ktWbJk6NChWlpa6urqBw4cEMCoq6urp6dn\nUVHRkydPnj9/7uXl5ebmJrTURliyBL/8gm7dxGVfIiDdsBpHtBEh54UwVfGkOcWRk4P+/SGmPz+c\n/XnBwbh8GRoaWLQIEybg8GFQ2CUnKyuLs9AcNmzY8ePHg4KCJkyYoKKiYmVlVVlZKY4Z+aygZ8+e\nffPmTTkh0j2ysrKGhoac01aKioqcWmKCa2yC48chIyOykEpyId2wGodpEaE0pziioqi4GTn78+zs\nUFSEyEhYW6NHDzg4wMxM3MfTzM3N582bx2KxnJ2dy8vLz5w5o6ysbGBgUFBQkJSUZGZmJvIZ+fjN\nESNGGBgYmJiYcHx0QMsDluHDh0+YMGH06NG9e/cuKChITU0V7ATXFygsxPbtuHhR9JYlEdINqxE8\nPDz++ecfOTk5dXX1NWvWCGDB1dXV399/5MiRCgoKFRUVGRkZbm5u58+fF0yPNKc4LlwAZdtq0WB/\nnrc3TE0xf774mtvduXOnT58+ioqKMjIy1dXV3FWsoqJilXhqjPBx0DY2NjY2NsIYDQoKunbtWlJS\nUmZmppqamre395gxY4QxyIfqajg7Y88ecvqoFtINqxHs7e0ZFRFKbYrjyRN8+y1oKejKuz/Pz692\nf569fW3bKtHRtWtXJyenqqqqZcuW3b17183NbcmSJYWFhTExMR4eHqKdiwOf62zy5MlCGpWRkTEy\nMjIS6+bcbdtgbPxVPc/9MpxuWOvW0a2DWTAtIpTaFEd0NCZNolMAd3/eu3c4cwaOjlBWhr09JkwQ\nVbNXQ0NDQ0NDzmtjY+OTJ08GBgZ269YtKipKyMcbjUF1k1rRkJ6OCxdw7hzdOhjG1KnYsQMrVqBt\nW7qlMAimRYRSm+KIjcXRo3SLAAC0b1+7P6+w8PP+PEdHDB0q2nlsbW1FsimoCSTQQX/8iCVLcPz4\n17uvrjG43bBWrKBbCoNgWkQonSmO168BQFWVbh116dkTq1Zh1SrcuIGQEHh4wMoKjo7o3JluZc1F\nAn3cihVYurS2XSOhHvb2iIxEeTndOgiNIp3lRs+eZXRVxeHDsWMHzp+HhgZcXDBxIg4fxsePdMv6\nMnUcNDPb1NchOhqlpZg2jW4dTIXTDevwYbp1EBpFOjuqREeDJ1Jhs9nZ2dn379+nUREf6tXP++EH\n5tfPq+Ogmdmm/jP//ouAADCgaxGjId2wmI0UpjjKyvD6Nf7bY/727VszM7M9e/Zs2LBh8uTJ3NM9\nDKJrV3h44OJFrFqFhASMGgUvLzx8SLcsPvBJcTCtTX0tbDbmz8f27ahbeoZQH243rK8eZkaEUpji\nuHAB5ubcf+3cuXPRokV79+79888/jY2Njxw5QqO0L8DZn5ecjHHj4O8PCwvs34/SUrplfYaPg2ZU\nm3oWi3XixIk///yzdMMGDB9ePWrU5cuX3d3d6VXFdJjaDevx48cHDx48e/asMIXzmw8zI0IpTHFE\nRYFnM0N+fj63ROWwYcMePXrE5ZKDAAAgAElEQVREk6xmw9mfFxKCiAgoKmLWLMyYgZgYJoShfBx0\naGjo/fv3PT09CwsL6W1TX15ebmFhcf/+fdXCwvubN2dPnfrixYv09HQxHdqRHrjdsJjEpUuX5syZ\nIycnl5KSMnnyZGGakrQIpkWE0pbi+PQJDx/yFtc3MDAICQkBUFNTc/jwYe7GYQmAsz8vJgZbtuDu\nXZiZwcMDt27RqIiPg+7UqdPWrVvDw8M3bdrUgdZS2dHR0dbW1r8sXz4lIaFdWNjmbdvU1dUFrur7\ndeHhgZ076RZRh6CgoDNnzjjY2Kxfv75v377Xr1+nZl5GRYSQvhRHg4Y+Tk5OAIyNjY2NjYcMGWJh\nYUGTMiHo0QOrVuHqVTg44M8/YWyMTZvw8iX1QijqqCIYb9686dq1K96/h4/Pt3p67wQqrfeV0q0b\nevVCcjL09emWUktlebna0aOIjkZMjLq6+tu3b6mZNzQ0dNOmTbGxsYMHD6Y3IuQgbQdVoqIwZw7v\ngIyMjL+/P11yRMzw4Rg+HJWVOH8eHh549w4zZ8LOjrIT7XxW0Jz6uQMGDHB3d9+yZQs1Ovhibm5+\n4MCBpxUVZQYGvr6+wp84+Lrw9sbGjXSL+I+MjEOFhWcTEytjYu4/eXL69GkDAwNqZmZORMhBqlIc\nbDYyMqCrS7cOMaOggB9+wLFjOHIEFRWYPh0uLkhKouAxDx8HLcL6uV+kurr6UV1KS0vZ//3YGhoa\nAQEBP/30k7W1tZaWlqOjIwBVVdV9+/aJVZWUwO2GRS9lZfDwwNat3RMT7+npjRs/fs2aNcHBwZz+\nfhSQmJioqak5bty4gICAkydPinWukpKS/XW5d+9ejx492Gw2i8V6/vw5m81+/fr13bt3eUck+MWt\nW6yRI5+/eMEUPeJ+wWKxFyxgnTz5fO5cdng4y8bmeWgou7qa9zOirdwvlo4qzScvL++3337jHcnM\nzBw0aBD3n2PGjDlz5oxYNUgzq1dj3To6u2ElJMDfHytW4IcfWgMrV65cuXIlxRLE0VGlMaqrq7/5\n5hveEQUFBe4zbVlZWQBKSkrcLoWy/5UrkNQXZ87Ayop+GdS/6NsXu3bh40fZ5GRUV/O+JeIdSmw2\ne968eQ8fPmT/R2VlZXBwsLu7+549eyoqKtjUMn/+/Hnz5lE8qTQzeTL73j0a5i0pYTs4sF1c2G/e\nNP3BDRs2nDt3TnxCRo8eXVNTM3fuXDabbWZmJr6J+ML3eq6urqZYhrgwMGBXVdEtglmYm5uvX79e\nVNb4rKDDw8OdnZ05zzH++uuvuXPnivIPAoFiON2wDh6kdNLISOzZA39/MGCLFZURYXOQnnKjubnQ\n0BBVJU8CX/jkoFeuXDl16tTS0lIAMTExlEsiiJSRI/HkCSjb1/X0KWxtcecO4uKY4J0B7Nu3T8ge\nm6JFeg6qREWBPLcXM3z++hkaGs6dO9fS0jKcxtwlQYRQ0w2rpga7diEmBlu3YsgQ8c7VEpgWEUrP\nLo7z50HO9IoZ/uVGf/jhh927d0+dOjUvL49iQQTRY26OmzdRXCzGKe7dg5UVAJw7xyjvDOZFhFJy\nUOXpU7RvD2VlunVIOXwc9MyZMwEMHTo0KipKS0uLckkEMeDujl27xGKZxYKvL1auxJ498PBgYEbS\n0NBw/vz5lpaWT548oVsLIDUpjpgYkt+ggDq3U0BAgI+PD6fEDGdE4N7yBGbB6YZVWiribljp6Vi1\nCo6OWLtW3B3vheGHH37o3r371KlTKysr6dYiLSmOs2cRGkq3COmnjoPmuGOyapZCRN4Nq6wM3t54\n+RJHj6JLF9HYFA+8ESH1u7AbIg27ON68QXU1qDpq9DVTx0FPnz4dwDTSr0QqsbeHgQF++kkEZQR4\njp+IQpm4YGZEKA21OP7+m9ENrqSIOg6as3aurKx88+aNgoJCZWVlv379UlJSaNJGECncblguLoIb\nefUKy5dDVRXR0czvnCDyiJDNZr99+1ZVVVVGiHyONKQ4oqJQ9wAwQUzUcdCchcacOXNcXFxGjRqV\nkpISHBxMkzCCGFiwAKamcHYW8FHe4cM4dIghx0+ag6giwsWLF//++++3bt2aMWOGrKysrKzswYMH\nR40aJZg1iU9xlJejpATdutGt46tAclpeEYRH4G5YT57A1haPHuHcOUnxzhBdy6u0tDQAq1evjoyM\nzMnJuXDhgpeXl8CqJH4XR0ICxo+nWwQjiImJWbx48bp160pKSsQ0BRUtr/Lz84U3QhANLe2GVVOD\nHTswbx7WroWvLyQqNhdtyyslJSUdHR0A6urqwrSDkfgUx+nTmDKFbhH0c/DgwVOnTrm7uw8fPnzS\npEllZWXimIWPg/7111+FbHmVX5e5c+cSH80U2rfH6NHN7YZ17x4mTgSAuDjo6IhVl/gQPiJ89OjR\nlClTHj9+zOl/unnzZmFqekj2QZVPn5CbC556k18tYWFh+/btGzRokLW1taWl5bVr18QxC59c5KJF\ni27evCknJyewURsbm7y8vJEjR3Kepdy+fdvR0fHKlSsCGySIkmXLMHPmF57Cs1jYsAE3b+L336Gh\nQZUyscCJCI2NjS9duiRYRFhcXJyfn5+Xl8epYS0nJ7d//36B9Uj2Lo7r1yUoxyVWuPVFAdTU1Ijp\n/5SPgx4xYoSBgYGJiQnHRwcEBLTUaEZGxubNm9PT03/77bc+ffpYWlqeO3dOBGIJIqFzZwwciKtX\nYWzM/wNpafDyYv7xk2bi7e0dHR0tTMurVq1aaWpqctMjS5cuFUaPZKc4Tp/Gjz/SLYIRzJ0719HR\n0d3dvbCw8MKFC6tWrRLHLHwctI2NjY2NjTBG5eXlfXx87t275+rqOm7cuE8M6F5OqIOnJ9zd+Tjo\nDx+wfDkqK3HsGDp3pkOZ6HFzcxMyImxIr169CgoKBPuuBO/i4DS4EnfVLQlhzpw53bp1i4uLU1NT\ni42NFdMfXT4Omrf1X0ZGhsCmBw4c+PfffwcHB/fo0UNgIwSx0KcPVFVx4waGD/88eOECAgKwcqWU\nnUEQPiJsyKVLl/iOP3jw4I8//uAdSU9Pr/dkUoJTHP/7H3R0pCCoEhWmpqampqZinYKPg46MjAwN\nDa2pqWGz2QUFBdxTWAIgKyvr6urq6uoqhEKCWPh33rxnM2Z4a2iYmJiscHSUX7YMnTtLxPGTliJ8\nRMiFxWLJy8sDaGw3iLq6up2dHe9Ibm5uvSZ1EpziiIoi+zcoho+DDg4O9vPzCwkJmTlzZnx8vEim\naSwkvHv3roeHB+9ITk7OiBEjRDIpoTGqqqqmrVkTra7+9/btFzdsKNTT6xseDqrabFOM8BHh3bt3\nPT09U1JSlJWVS0tLDQwMtm3bxvfUuLKy8nDeoATo0KFDvSZ1EpziuHoVa9bQLeLrgo+DlpWV1dfX\nP3bsmKmp6Z49e0QyTWMh4eDBgy9cuMA7smDBAhF3XSQ04Pbt2yNGjPjW1hZWVhZTpozr0SNBSr0z\nRBERurq6+vv7jxw5UkFBoaKiIiMjw83N7fz584LpkdQUR14eevaESFP5hC/CZx+0oqJiREQEi8Xa\ns2ePkLuvWSwW54UwBwQIIqddu3YlJSUwNERaGmvTpgpJ9BfNJjg42MvLS11dffny5VacrgItRFZW\n1tDQUEFBAYCioiJnV7XAeiQ1xREdTQpAUw8fBx0eHj5s2DBfX98XL16sXbtWAKN37961trb+9ttv\nNTU1VVVVra2tSWcWRtG3b182m+3t7X00IWHGjBkuwpRPYjyciFBBQcHU1PTBgwcCWBg+fPiECRPW\nrVv3119/+fn5TZw4UUeIYzuSelAlPh4TJtAt4quDz0JAWVm5b9++APz8/AQzKtqQkCAOjhw58vff\nf+fl5fn5+Wlra9MtR4wIHxEGBQVdu3YtKSkpMzNTTU3N29t7zJgxAuuRyBTHv/9CRYU0uKIesZQb\n5YSEnGOEwoeEBHEgIyNjZWUlWMgvWYSHhxcVFRkaGu7du1ewiFBGRsbIyMjIyAjAihUrDITL10tk\niiM6GpMm0S3ia0Qs5UY5IeHo0aN79+5dUFCQmpoqTEhIIAiD8BEhL0+fPhXSgkTu4jh7FiEhdIv4\nGhFLudGgoCBvb285ObnMzExZWVlvb++NGzcKLZVAaBmiKjfKi7XQp3gkr9zo27dgsfDtt3TrkEAe\nPoSdHaZPR3p67cjFizA2hr19Mzs68t9mJ2RxGd6QkECgC3E0oLC3txfSguSlOOLiaosaEppBdXX1\n543C+/bBxwfduuHnn3H0KABkZaF/f5SXY9iw5ljjs4IODQ0VstwogcAcmNaAQvJ2cZw5QzbYfYHi\nYpw6hbw8AKWlpaWlpbXjz55BUxMdOoD7dNrcHFu2YPt2/PJLcwzzWUF36tRp69atIpFNINCO8BGh\naJGwXRwVFXj+HL160a2Dqbx/D2trqKrCxASjRwNQVVVtz62XoK6Ox4+hrg4VldqRO3fQpw8UFZtZ\n0oSPg961axdvJChMLQ4CgXZCQ0M3bdokTLlR0SJhKY6LF2FmRrcIJnHpEs6eRUYGHB0xbx7atUNi\nYqMfXrSotmbvsmVISMC7d+jaFU5OqKrCsmXNmY1/saSrV69+S54JEKQCpkWEEraLIyoKdavlfHVU\nViIpCSUlmDEDAPLyMG0aNm9uVuflvn0RFlZ/sCUdD/jMoaGhwekcQSBIAUyLCCUpxVFdjXv3oKVF\ntw76WL0a167BwAC2trUjixZROT8fB11TU6Orqzt27FjO6ZLAwEAqBREIooVpEaEkpTiSkyH0xkRJ\norgYZ87g6lWwWAgPh4wM1q+nVxEfBz19+nTqdRAIYoJpEaEkpTiioj6vHKWVykqkp+P779GmDW7e\nBJuN9evRsyfdsmrh46B5t+IL01GFQGACTIsIJSnFkZKCoCC6RYiNBw/g6YniYpiYYMgQALCwoFtT\nfcTbUYVAoB0qI8I7d+4sWbKEdyQnJ0dXV5d3RGJSHJmZ0NGBLJ+jEpJKdTUyMnD6NCwsMHYsunTB\nrl3o1o1uWU1BUUcVAoEuqIwIv/vuuy82oJCYFMfp01J1PqWiAlZWGDwYkybV7qNo0wZ1u5ExEIo6\nqhAIdMG0iFBiUhxXrmD1arpFCEFJCaKicPIkjIzwyy9QVMTFi1RrqK5GYiJOncKuXYIZ4OOgRdhR\n5YsUFBSEh4fzjmRmZvYix5YIooNpEaFkpDgePYK6ukQ2uHr3rrbx8e+/o0MH/PknunShTcyLF0hP\nR1WVwAb4OGjh6+c2H1VV1XpNNq9evcrpLUQgiASmRYSSkeKQuAbexcU4ehTx8ejWDQcOAMCvv1Kt\nITMTR49qp6crcgshqatj4UJ4eQlsko+Djo2N5bSO9/PzCwsLE75CYxO0b99+3LhxvCPHjx8nTWMJ\nIoTKiLA5SEaK49w5nDpFt4hm8Pw52rWDsjJyc9G7N06dAvXRyfv3aNcObDYiIzFlyrft239Xd9Ep\nDHUcdERERERERHp6+vHjxwGw2ezc3NzZs2eLajICgXqojAibgwSkOJ4/h6Iiox+gVVXhjz9w9iza\nt8fmzVBWBvVt6V+8QFgY4uMxYgQCAiAjg4AAALh8WYST1HHQpqamOjo627dvX7p0KWekc+fOIpyM\nQKAeKiPC5iABKY6zZxna4KqsDG/fols3vHsHJSUcP07DX5Hy8tpadElJ6NIFUVFoukSiqir27RN4\ntjoOukOHDh06dOAULnj//j2bzf5cN49AkDSYGRFKQIojOhp//EG3iLqkp2PPHty7h1WrMGUKOnSA\noyOlAthsXLqEP//E06eIjETHjtScsazjoG/cuOHi4nL+/PmsrKwff/xRUVFx+/btk6VpLyTha4KZ\nESHTUxzv36O8HAxpysVJ7wLIyYGbG+oe+aGUhw9x/Tr8/dGnD5XT1nHQ7u7uYWFhampq/v7+sbGx\nmpqapqamxEETJBRmRoRMT3GcO0d/g6uqKpw7h9BQtGqFo0chIwMHB6o1lJTgyBEcPw5XVzg4oG9f\nGraF1Gt5JSsrO2DAgPLy8tevXw8bNqx9+/ZykrgRkkAAANy4cUNXV/f169dXrlwZMGDA999/HxUV\nJYxBNpv95s0bNpstjBGmN41lwga769eRkYFt23DsWDM7j4ietWuhooL4eBr+NvBQx0F/+PChpqbm\n/PnzI0aMAFBVVUUcNEFyqRcR3rp1y9/fXwA7ixcvBnDr1q2BAwfq6+traWmlpqYKrIrRKQ4WC8+e\ngZbVfXIyZszAuXMAYGKCdevQvTulAoqK4Of3ObO8Zw8WLKhNsNBHnRSHnZ3d0KFDi4uLY2Nj8/Ly\nlixZYmVlRZcyAkFI6kWEAARbcKSlpQFYvXp1ZGSkjo5OUVHRrFmzrly5IpgqZqY4Xr58uXXrVtXU\n1PHt2n3/6VPr5rQLESHTpqFHDwQEoF8/SuflcvYs/voLc+fC25seAY3w+b+hsrLS1dV1zJgxHTp0\n6NatW05OztSpU6fQHuwQCIIi2ohQSUlJR0cHgLq6enV1tcB2xL2L4/bt256enh8/fiwtLZ04caK/\nv7/Ml7IEnz590tfXt7Cw8O3V65CS0qlff5WVle3evfuiZnQP0dLSys7OjoiIyMvL8/Hx4X2rpKTk\n0qVLjVYT/PdfnD6NxYsB4MSJZv50oqSwEIcOYeJE6OnB2ho8RbWYw2cHferUqdOnT9d7OyEhISQk\nhLkRGYHQOKKKCB89ejRlypTCwsIjR47Y29tv3rxZXV1dYFViTXF8+PBh1qxZx44d09LS+vTp08SJ\nE48fPz5z5symv5Wdnd29e3dtbW2lI0d+SkoyMDDQ09MLCwtLT0+fNWuWWTOaxvL1wiUlJREREXze\nKirCr7/i33/x3+4aqqmqgoMDyssxdy64x7IZSWsA8vLyu3fvNjY23rJlS0/GtBIgEITE29t7woQJ\nHTt27NGjR3Z29pw5cwSrDV1cXJyfn5+Xl8fpzCInJ7d//36BVX1Ocbx6hUuXAEBGBubmaNeuWSNN\nEhsbO3bsWC0tLQCtW7cODw+vrKwsKytzdnZ+9+5dWVmZr6+vqalpvZGuXbtWVlYCeLF7t52RUWlp\nadeuXSMjI7OyshYuXDho0CB5efk//viD25imuLjYwcGhurq6c+fOnC+eOHEiOzt7xowZXl5e1dXV\nnM9v3749PT09MjLSyspq5syZFRUVABYuXDht4MCz330XwWK1Cgt7sG7d9OnT3d3dS0tL58+f/+7d\nu48fP+7du7d///6LFi3Kyclp3bp1QECAkZGRwL/wz7x/j9RUmJtDTg47d6JjRxHYFDOtAQQFBaWl\npf3vf/+LjIwsLCxUVlbW0dEZOnTo0KFD+/fvLytExW4WiyUvLy86tQRCC5CRkRn23/pIS0tLS9Dm\np61atdLU1NTU1OT8c6lw677PKY5Xr/DoUe3omze17viLI03y6NGjPjwbdTkudevWrVpaWj4+PoWF\nhUZGRvn5+fv27eMdefz48YcPH06ePLlp06bBgwfPmjVr0aJFGRkZ8+bNMzIyCgsLi46OLi4u5jro\nwMBAGxsbV1fXpKSkuLg47nTx8fHa2tr+/v6czy9duvTVq1d2dnYPHjywHzVq+j//JI8YEXj48LTo\naOTnZ4eG3rhx4/Xr18bGxu7u7jt37tTS0lqzZs21a9eSkpKSk5Pl5eWvX7/+8uVLfX39hw8fCvM7\nx5Mn2LIF2dlwd68dkQTvDI6DbtOmjampqampKWfo48ePWVlZJ0+enDdvnqam5s2bN1vqZO/evevp\n6ZmSkqKsrFxaWmpgYLBt27b+/fuLXj6BQAe9evUqKCgQ7LufUxz9+2PVqjrvNWekSXr37p2cnMz9\n5z///JOWlnbv3r1Zs2YB4MTHHz9+rDdSXl4+f/78jRs3tmnTZvjw4YsWLYqIiJg6dWp6erqJicnE\niRN79eo1cuRIrtm8vLz58+cDGDVqFG/tSUdHx61bt3I///btW854m3XrLqalJerpsQoLuaXQxo4d\nKyMjw23mm52dzUl5GxoaGhoa/vzzz7m5uU5OTgA6depUVlamoqLS/N9DffLyMGMGduwQ3AJN8Fkd\nKygo+Pr6FhcXx8TEZGZmCrAEdnV19fT0LCoqevLkyfPnz728vNzc3EShlkBgBJc4aYcGZGVl6dYl\nJiaGE9pzKS0tLSwsFJOwiRMnJiQk3Lx5E0BVVdW6devatWs3cODA69evA8jPz2ez2crKyg1HWrVq\ntXbt2uvXr4eHhz9//lxeXt7e3v7HH3+cPXt2XFychobGoUOHuLMMGjTo8uXLAFJSUjgpDg4RERHm\n5uafP//oEbumBsDvffro/Pzz7tDQmTNncneR13tgO2DAgKSkJACXL1/28vIaMGCAvr5+SEjIjh07\nrKysWuydP3zAvn0wNkZmJgCYmUFfv+W/Tvrhs5mGE3/9+uuvGhoaghmVlZU1NDTkPDtWVFQ0MDCg\netcOgSBOuOmOeujo6NTrqvXXX3/VK58r1l0c7du3P3bs2PLlyz98+FBeXj5p0iQHB4ePHz86OTmN\nHz++oqIiJCQEwKJFi+qNAJCRkencufOvv/7q5uZ24sSJyZMnZ2RkODk5XbhwQVZWljft7unpOWfO\nnBMnTnTt2rUfz8a4YcOGOTk5qaqqyrJY+zt37pSenpubGxYWZmlp6efnd+nSpb59+xYUFPBtPLZ0\n6VInJyczM7Py8vKDBw9qaGjMnz/fzMysrKxs5cqVLfstsFiYMgU2Nvj7bwiz7mYAMnyPRRUVFbVt\n27Zt27aCGV2xYkV2dvbo0aN79+5dUFCQmpqqo6PTzG7KnB5uBw8eFGxqgsSxcePGYcOGWTCvobJI\n4DhoTqjOpaamRphHO0zn9Ws4O8PLCzxZESp4/BiXLsHZmdJJGyDa65n/wlaYXUQAgoKCOJn+zMxM\nNTU1b2/vMWPGCGOQQKCRLVu2NBxcsWKFYNaYeVBFBFRV4eJFWFpCTQ0NNuyKl8eP4eODigosX07p\nvOJHLJkHGRkZIyOj5uyMYbFYt2/f5h0pKSlhQkUbAoHLN9984+rqGhQUJJLUhASUGxWAhASsWUNb\n2YrXr+HlBW1temYXJxSlhht76v3o0aN6W0qfPXvWneIz+ARCkzg7Ox87duynn34SyaMURtfiEIz0\ndBw9iqgoUFnNNSUF69fDzQ2WlhBdiymmQZGDbuyp98CBAznVILk0fKhCINBOQkKCqExJT4qDzUZe\nHgYMwIgRGDGC0ql//BHt2mH7dtpqd1CFeJ9UsFgszovGnnoTCJKFwKlnLkwvN9pMHj6EtXVt8TnK\nyM0FZ1PDsWMIDpZ67wwxOei7d+9aW1t/++23mpqaqqqq1tbWeXl54piIQKCYp0+fCmlBGlIcV67A\nyQmbN8PDg6IZ797F9Onw9cWnTxTNyAzEkuJwdXX19/cfOXKkgoJCRUVFRkaGm5vb+fPnm/n1jIyM\nNi3vBfno0aOCggJhLv1Xr15xzzUJQElJSYcOHQT+ellZmZC3Lu36v/vuu65du7b0i9nZ2cOYXbCG\nF+uW1zyrdz2zWKySkhKBDyLW4/Xr19wT2MLTgmtAVxdN7oUV/nrmovb+/fiEhEQjo387dGjRuUq+\nfPr0qaysTFQ7EV6/fm1hYcH7cEK017NYHLQwB1XMzMwqKyvfvHnT0kmTkpLev38vcLGn8vLyGzdu\nGAjRvP3ChQvjx48X+Ou5ublKSkrC6M/IyDA0NBRYgPD6P378KIAAExMTfck55WVvb9+izze8nu/f\nv3/jxo1u3bqJRE9KSoq2trYAC5qGVFdXX7lypYnadYaPHvUtKflTT6851oS8ngGosFgDX7y40aNH\nMbC7omI8m43iYoGtcSkuLn769OnQoUOFNwUgJSVFU1OTt9eliK9nthhYvny5hYWFr69vSEjIunXr\nJkyYsGrVKnFMxMvOnTtPnDgh8Nf//fdfzjlUgTExMRHm67t27YqMjBT468+fP58xY4YwAoTUv3v3\n7oiICGEsfCVcv37d29tbVNYcHR0fP34sElOVlZXm5uYiMcUW+nquh5AXJy8pKSleXl6isjZv3ryH\nDx+KylpDxLKCJgdVCAQCQXhoPqhCIBAIhMaQ3oIABAKBIOEQB00gEAgMpZWvry/dGkRD165d+/bt\nK3BVbyUlpX79+nXp0kVgAdra2sIUmerSpYuQ+vv370+j/q5du2pqagpVVf3rQFVVVVNTU5gNkbxo\naGhoaGiIpLhHq1atBg8eLGShNC5CXs/1EPLi5IXz+xdmRykvffr00dTUFF9xFf7lRgkEAoFAOyTF\nQSAQCAyFOGgCgUBgKMRBEwgEAkMhDppAIBAYCnHQBAKBwFCIgyYQCASGQhw0gUAgMBSJdNDl5eW2\ntrZmZma6urqpqam8b1VXV6upqenq6urq6i5evFisMlgslo6OTn5+PvUCmpiFGgF79+4dO3bskCFD\nUlJSqJ9dgmCz2StWrBg7dqyxsfHdu3eb+VZLTXE4ePAg3wbkLbXGYrFcXFzMzc11dHROnToljKkm\nblUBrHHlNbzvBLDW0mu1aWGN3RHCIr5CeeLjjz/+8PT0ZLPZiYmJY8eO5X3rwYMHkydPpkaGv7+/\ngoJCvWKP1AhoYhYKBGRnZw8fPpzFYmVlZY0cOZLi2SWLq1evmpub19TUXL161crKqplvtdQUm80e\nN26cvLx8UFCQ8MLOnTvn5OTEZrNfvnzZqVOnqqoqgU01casKYI0D3/tOAGstvVabMNXEHSEkErmC\nHjVqlKurKwA2m62qqsr7Vl5eXmFh4bRp06ytrW/duiU+Dffu3UtLSxs5cmS9cWoENDELBQJiY2Pt\n7Ozk5OS0tbWjo6Mpnl2ySEpK0tfXl5GR0dPTq7eEbOKtlpoCEBcXt3nzZpEI69Kly7JlywCoqKgo\nKyuzv3TYuAlTTdyqAlhD4/edANZaeq02YaqJO0JIJNJBf/fdd717916wYIG5ufnSpUt531JTU1uy\nZElkZOT69es5BezFIa+UO+8AAAbZSURBVKCmpsbDw+O3337jdI2hXkATs1Ag4MWLFzk5ORYWFsbG\nxjdv3mymsK+TkpISTgNvRUVFJSUlbhvlpt9qqSkArVu3blFFiCasDRkyREtL686dOxYWFl5eXnJy\ncgKbauJWFcBaE/edANZaeq02YaqJO0JIxFIPWty8f/9eSUnpwIED3t7eZmZmDx8+5P6HjRw5kvPX\ndciQIbKysm/fvv3mm29ELiA4OHj8+PEaGhoN36JGQBOzUCCgbdu2z549i4mJefHixbBhw54+faqg\noEDZ7JKFmpoap/1gRUVFaWmpvLx8c95qqSnRCmOz2evXrz9z5sz27dub0wSuCVNN3KoCWGvivhPA\nWkuv1SZMNXFHCIlErqA3bNgQHBwMQEFBod46Yvv27ZxAr6ioSEZGpplRVUtJS0s7d+6cpaVlVlbW\n3LlzeZ9XUCOgiVkoEGBoaNiuXTs5Obm2bdvKysryLj2o+fElCGNjY044nJGRUa+vUBNvtdSUaIVF\nRUWlpqYmJSU1s0VnE6aauFUFsNbEfSeAtZZeq02YauKOEBYR5rMp4+XLlxMmTNDX1//+++9jY2PZ\nbPaDBw969uzJZrPfvn07ffp0AwMDIyOj5ORkcSsxNjbmPKygWEDDWagUUFNTs2rVKj09vSFDhnD6\nENL1+2c+NTU1K1eutLS0HDdu3L1799g8v6uGbwlsisOuXbua/5CwCWsuLi49evQY8h8VFRUCm2p4\nqwr/Y7J57jthrLX0Wm36v7LeHSEqSLlRAoFAYCgSmeIgEAiErwHioAkEAoGhEAdNIBAIDIU4aAKB\nQGAoxEETCAQCQyEOmkAgEBgKcdAEAoHAUIiDJhAIBIZCHDSBQCAwFOKgCQQCgaEQB00gEAgMhTho\nAoFAYCjEQRMIBAJDIQ6aQCAQGApx0AQCgcBQiINuMU+fPlVRUTH5D3d394iIiICAgJKSkoiICADc\nF02jpaXVonlPnDjh4+MjoGgCoYVkZGRYWlrSreJrRyJ7EtJOnz59rly5Um/w3r17ERER06dP5zjo\n6dOn0yGNQCBID2QFLQJOnDjh6+u7ffv29PT0yMhI7ouqqipnZ2d9fX0jI6PExEQAxcXFEyZMMDc3\nnzNnTmVlJa8RCwuLjIwMAPHx8bNnz/748aONjY25ubm5ufmJEye4Hzty5EhgYCCAiooKXV1dAA1n\nycnJmTRpkrW19dSpU1+/fk3lr4IgxTS80updtA0/EBUV5ezsrK2tfePGjXrX87t372xtba2trV1c\nXPT09PjaJ5AVtCA8fvzYxMSE85rTdxLA0qVLX716ZWdnp62tzXlx4MABeXn569evv3z5Ul9f/+HD\nh4GBgTY2Nq6urklJSXFxcbw2HRwcjh8/rqure/jwYRcXl6KiInt7++nTpycnJwcGBk6bNq0xMSEh\nIfVmiY+P19bW9vf3j46OLi4uVlNTE9tvgvAV0fBKq3fRNvwAgOzs7Js3bxYUFNS7nrdt26arq/vL\nL7/ExcVdvHiRr326f2L6IQ5aEOqlOAoLC/l+7Pbt27m5uU5OTgA6depUVlaWl5c3f/58AKNGjarX\nmH3y5MkbN25cs2bNvXv3DA0NX7x4cfHixcTERBaLVVNT09D4p0+fGpvF0dFx69atEydO7NWrF6ex\nPIEgPA2vtHoX7YkTJ+p9AICJiYmcnFybNm3qXc+3b99evnw5AH19/cbsq6io0PXDMgSS4hAl3A68\nnBcDBgzQ19cPCQnZsWOHlZWViorKoEGDLl++DCAlJaVeikNFRUVPT2/FihX29vYyMjK///67jo7O\n7t27Z86cydvYV1ZW9uXLlwA4iw6+s0RERJibm8fFxWloaBw6dIian50g9TS80updtA0/AEBOTg5A\nw+t54MCB165dA3D9+vXG7NP3szIF4qBFRufOnXNzc8PCwrgv5s+fX1BQYGZmZmFhMWjQIACenp4x\nMTGmpqb79u3r169fPQsODg6HDx92cHAAYGlpGRMTY2trGx8fX1BQwMn0ATAxMcnMzJw0aVJaWpqy\nsjKAhrMMGzZs8eLFhoaGsbGxtra2lP4WCFJEamqqwX8kJiY2vNJQ96Ll+wEODa/nlStXXr9+fcKE\nCQkJCUpKSk1//atFhnd1RiAQCNQQHx8vKys7fvz45OTkTZs2nTlzhm5FTIQ4aAKBQAPPnj1zcHBQ\nUVFhsVibN2/W0dGhWxETIQ6aQCAQGArJQRMIBAJDIQ6aQCAQGApx0AQCgcBQiIMmEAgEhkIcNIFA\nIDAU4qAJBAKBoRAHTSAQCAyFOGgCgUBgKMRBEwgEAkMhDppAIBAYCnHQBAKBwFCIgyYQCASGQhw0\ngUAgMJT/A/2yWM6mr0aQAAAAAElFTkSuQmCC\n"
      }
     ],
     "prompt_number": 3
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "XYcoef"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "pyout",
       "prompt_number": 4,
       "text": [
        "array([ 3.2,  0.9])"
       ]
      }
     ],
     "prompt_number": 4
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "station = \"8518750\"\n",
      "startdate = \"20121030\"\n",
      "stopdate = \"20121031\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 5
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R -i station,startdate,stopdate -o urltotal\n",
      "url1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\n",
      "url2 = \"SixMin_Verified_Water_Level.ascii?\"\n",
      "url3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" #return stationId\n",
      "url4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" #return datum\n",
      "url5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" #return record date-time\n",
      "url6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\" #return water level value\n",
      "url7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" #return quality flag\n",
      "url8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" #return quality flag\n",
      "url9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" #return quality flag\n",
      "url10 = \"WATERLEVEL_6MIN_VFD_PX.T\" #return quality flag\n",
      "url11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" # station ID goes here\n",
      "url12 = \"%22\"\n",
      "url13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"#we want MLLW as the datum\n",
      "url14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\n",
      "url15 = \"%22\"\n",
      "url16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\n",
      "url17 = \"%22\"\n",
      "urla = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11)\n",
      "urlb = paste(station,url12,url13,url14,startdate,url15,url16,stopdate,url17,sep=\"\")\n",
      "urltotal = paste(urla,urlb)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 8
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "print urltotal[0]"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "stream",
       "stream": "stdout",
       "text": [
        "http://opendap.co-ops.nos.noaa.gov/dods/IOOS/ SixMin_Verified_Water_Level.ascii? WATERLEVEL_6MIN_VFD_PX._STATION_ID, WATERLEVEL_6MIN_VFD_PX._DATUM, WATERLEVEL_6MIN_VFD_PX.DATE_TIME, WATERLEVEL_6MIN_VFD_PX.WL_VALUE, WATERLEVEL_6MIN_VFD_PX.I, WATERLEVEL_6MIN_VFD_PX.F, WATERLEVEL_6MIN_VFD_PX.R, WATERLEVEL_6MIN_VFD_PX.T &WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22 8518750%22&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%2220121030%22&WATERLEVEL_6MIN_VFD_PX._END_DATE=%2220121031%22\n"
       ]
      }
     ],
     "prompt_number": 11
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "#download the data\n",
      "cat(\"Contacting server...\\n\"); flush.console()\n",
      "dat = getURL(urltotal) #use RCurl to retrieve text into a vector 'dat'\n",
      "cat(\"Data returned...\\n\"); flush.console()\n",
      "Sys.sleep(2) #If you access data in a loop, be courteous and give the server\n",
      "#a short break between requests\n",
      "#cleanup\n",
      "rm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\n",
      "rm(url15,url16,url17)\n",
      "\n",
      "con = textConnection(dat) #create text Connection to dat vector\n",
      "all.lines = readLines(con) #read lines of text into separate slots in a vector\n",
      "close(con) #close connection to dat vector\n",
      "\n",
      "if (length(grep('^Error',all.lines))>0) { #check for error in retrieval\n",
      "\tcat(\"There was an error...\\n\")\n",
      "\tcat(dat,\"\\n\") #print contents of dat to show error\n",
      "\tflush.console()\n",
      "} else { #retrieval was successful, parse the text\n",
      "\n",
      "\t#The column headers are typically preceded by a line of dashes\n",
      "\theaderlines = grep(\"^--------\",all.lines) #find index of headers (-1)\n",
      "\n",
      "\t#read column header names into a vector\n",
      "\tcon = textConnection(dat)\n",
      "\theaders = scan(con, skip = headerlines, nlines = 1, sep = \",\",\n",
      "\t\t\twhat = \"character\", strip.white = TRUE)\n",
      "\tclose(con)\n",
      "\n",
      "\t#read rest of the data into a data frame 'df'\n",
      "\tcon = textConnection(dat)\n",
      "\tdf = read.table(con, skip = headerlines+1, sep = \",\", header = FALSE,\n",
      "\t\t\tquote = \"\\\"\", col.names = headers, strip.white = TRUE,\n",
      "\t\t\tstringsAsFactors = FALSE)\n",
      "\tclose(con)\n",
      "\n",
      "\t###########################################################################\n",
      "\t#The following operations will need to be altered if you change the \n",
      "\t#fields or data type being returned by the OPeNDAP server\n",
      "\n",
      "\t#Convert the time column to POSIX time (seconds since 1970-01-01 00:00:00)\n",
      "\tdf[,3] = as.POSIXct(strptime(df[,3],format = \"%b %d %Y %I:%M%p\",\n",
      "\t\t\t\t\ttz = \"GMT\"))\n",
      "\n",
      "\t#Give the columns shorter names\n",
      "\tnames(df) = c(\"stationId\",\"datum\",\"TimeUTC\",\"TideHT\",\"Flag.Inferred\",\n",
      "\t\t\t\"Flag.Flat.Tol\",\"Flag.Rate.Tol\",\"Flag.Temp.Tol\")\n",
      "\n",
      "\t#Uncomment this if you want to plot the data\n",
      "\tplot(df$TimeUTC, df$TideHT, type = \"l\",\n",
      "\t\t\txlab = \"Date\",ylab = \"Tide Height, meters\")\n",
      "\n",
      "\t#Save data automatically to a .csv file. \n",
      "\tfilename = paste(\"Station_\",station,\"_\",startdate,\"-\",enddate,\n",
      "\t\t\t\".csv\",sep = \"\") #make file name\n",
      "\twrite.csv(df,filename,row.names = FALSE, quote = FALSE) #write file to disk\n",
      "\tcat(\"Saved to \",getwd(),\"/\",filename,\"\\n\",sep = \"\")\n",
      "\n",
      "\t#Alternate file save method lets user specify file name at run time\n",
      "\t#Uncomment this if you wish to use it instead of the automated file \n",
      "\t#output above\n",
      "#\twrite.csv(df,file.choose(),row.names = FALSE, quote = FALSE)\n",
      "\n",
      "\t#cleanup\n",
      "\trm(dat,con,all.lines,startdate,enddate,filename,headerlines, headers)\n",
      "\n",
      "} #end of if-else statement\n",
      "# coops_tide_ht_retrieval.R\n",
      "# This script will download a set of verified tide height data from a NOAA \n",
      "# CO-OPS DODS/OPeNDAP server and parse it into a data frame to be saved to disk \n",
      "#\n",
      "# Station list: http://opendap.co-ops.nos.noaa.gov/stations/index.jsp\n",
      "# OPeNDAP server gateway: http://opendap.co-ops.nos.noaa.gov/dods/\n",
      "# The gateway has links to several other data types available including \n",
      "# water temperature, air temperature, wind etc. \n",
      "# See http://docs.opendap.org/index.php/UserGuideOPeNDAPMessages for info on\n",
      "# structuring OPeNDAP queries.\n",
      "#\n",
      "# Six-minute water level data can only be retrieved 1 month at a time. Other\n",
      "# fields such as water temperature could return up to 1 year at a time. This \n",
      "# script only deals with 6-minute tide height data, but could serve as a guide \n",
      "# for accessing other data types from NOAA CO-OPS. \n",
      "# The tide height is reported in meters.\n",
      "# Author: Luke Miller Feb 2011\n",
      "###############################################################################\n",
      "require(RCurl)\n",
      "\n",
      "noquote(print(\"Enter NOAA tide station number (i.e. Monterey = 9413450): \"))\n",
      "station = scan(\"\",what = character(),nlines = 1)\n",
      "noquote(print(\"Enter start date (format: 20080123 = Jan 23, 2008): \"))\n",
      "startdate = scan(\"\",what = character(),nlines = 1) #get one line of values\n",
      "noquote(print(\"Enter end date (format: 20081231 = Dec 31, 2008): \"))\n",
      "enddate = scan(\"\",what = character(),nlines = 1)\n",
      "\n",
      "#OPeNDAP query for 6-minute verified water level looks like this (on 1 line):\n",
      "#http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\n",
      "#SixMin_Verified_Water_Level.ascii?\n",
      "#WATERLEVEL_6MIN_VFD_PX._STATION_ID,\n",
      "#WATERLEVEL_6MIN_VFD_PX._DATUM,\n",
      "#WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\n",
      "#WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\n",
      "#WATERLEVEL_6MIN_VFD_PX.I,\n",
      "#WATERLEVEL_6MIN_VFD_PX.F,\n",
      "#WATERLEVEL_6MIN_VFD_PX.R,\n",
      "#WATERLEVEL_6MIN_VFD_PX.T\n",
      "#&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%229449880%22\n",
      "#&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\n",
      "#&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%2220080801%22\n",
      "#&WATERLEVEL_6MIN_VFD_PX._END_DATE=%2220080808%22\n",
      "\n",
      "#The parts of the url to be assembled:\n",
      "url1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\n",
      "url2 = \"SixMin_Verified_Water_Level.ascii?\"\n",
      "url3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" #return stationId\n",
      "url4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" #return datum\n",
      "url5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" #return record date-time\n",
      "url6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\" #return water level value\n",
      "url7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" #return quality flag\n",
      "url8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" #return quality flag\n",
      "url9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" #return quality flag\n",
      "url10 = \"WATERLEVEL_6MIN_VFD_PX.T\" #return quality flag\n",
      "#The remaining parts of the url specify how to filter the data on the server \n",
      "#to only retrieve the desired station and date range. Values must be enclosed\n",
      "#in ascii double-quotes, which are represented by the code %22\n",
      "url11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" #...insert stationId here...\n",
      "url12 = \"%22\"\n",
      "url13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"#we want MLLW as the datum\n",
      "url14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" #...insert start date here...\n",
      "url15 = \"%22\"\n",
      "url16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" #...insert end date here...\n",
      "url17 = \"%22\"\n",
      "\n",
      "#Assemble the URL\n",
      "urltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\n",
      "\t\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\n",
      "\n",
      "#Download the data\n",
      "cat(\"Contacting server...\\n\"); flush.console()\n",
      "dat = getURL(urltotal) #use RCurl to retrieve text into a vector 'dat'\n",
      "cat(\"Data returned...\\n\"); flush.console()\n",
      "Sys.sleep(2) #If you access data in a loop, be courteous and give the server\n",
      "\t\t\t #a short break between requests\n",
      "#cleanup\n",
      "rm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\n",
      "rm(url15,url16,url17)\n",
      "\n",
      "con = textConnection(dat) #create text Connection to dat vector\n",
      "all.lines = readLines(con) #read lines of text into separate slots in a vector\n",
      "close(con) #close connection to dat vector\n",
      "\n",
      "if (length(grep('^Error',all.lines))>0) { #check for error in retrieval\n",
      "\tcat(\"There was an error...\\n\")\n",
      "\tcat(dat,\"\\n\") #print contents of dat to show error\n",
      "\tflush.console()\n",
      "} else { #retrieval was successful, parse the text\n",
      "\n",
      "\t#The column headers are typically preceded by a line of dashes\n",
      "\theaderlines = grep(\"^--------\",all.lines) #find index of headers (-1)\n",
      "\n",
      "\t#read column header names into a vector\n",
      "\tcon = textConnection(dat)\n",
      "\theaders = scan(con, skip = headerlines, nlines = 1, sep = \",\",\n",
      "\t\t\twhat = \"character\", strip.white = TRUE)\n",
      "\tclose(con)\n",
      "\n",
      "\t#read rest of the data into a data frame 'df'\n",
      "\tcon = textConnection(dat)\n",
      "\tdf = read.table(con, skip = headerlines+1, sep = \",\", header = FALSE,\n",
      "\t\t\tquote = \"\\\"\", col.names = headers, strip.white = TRUE,\n",
      "\t\t\tstringsAsFactors = FALSE)\n",
      "\tclose(con)\n",
      "\n",
      "\t###########################################################################\n",
      "\t#The following operations will need to be altered if you change the \n",
      "\t#fields or data type being returned by the OPeNDAP server\n",
      "\n",
      "\t#Convert the time column to POSIX time (seconds since 1970-01-01 00:00:00)\n",
      "\tdf[,3] = as.POSIXct(strptime(df[,3],format = \"%b %d %Y %I:%M%p\",\n",
      "\t\t\t\t\ttz = \"GMT\"))\n",
      "\n",
      "\t#Give the columns shorter names\n",
      "\tnames(df) = c(\"stationId\",\"datum\",\"TimeUTC\",\"TideHT\",\"Flag.Inferred\",\n",
      "\t\t\t\"Flag.Flat.Tol\",\"Flag.Rate.Tol\",\"Flag.Temp.Tol\")\n",
      "\n",
      "\t#Uncomment this if you want to plot the data\n",
      "#\tplot(df$TimeUTC, df$TideHT, type = \"l\",\n",
      "#\t\t\txlab = \"Date\",ylab = \"Tide Height, meters\")\n",
      "\n",
      "\t#Save data automatically to a .csv file. \n",
      "\tfilename = paste(\"Station_\",station,\"_\",startdate,\"-\",enddate,\n",
      "\t\t\t\".csv\",sep = \"\") #make file name\n",
      "\twrite.csv(df,filename,row.names = FALSE, quote = FALSE) #write file to disk\n",
      "\tcat(\"Saved to \",getwd(),\"/\",filename,\"\\n\",sep = \"\")\n",
      "\tflush.console()\n",
      "\n",
      "\t#Alternate file save method lets user specify file name at run time\n",
      "\t#Uncomment this if you wish to use it instead of the automated file \n",
      "\t#output above\n",
      "#\twrite.csv(df,file.choose(),row.names = FALSE, quote = FALSE)\n",
      "\n",
      "\t#cleanup\n",
      "\trm(dat,con,all.lines,startdate,enddate,filename,headerlines, headers)\n",
      "\n",
      "} #end of if-else statement"
     ],
     "language": "python",
     "metadata": {},
     "outputs": []
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R\n",
      "require(RCurl)"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "output_type": "display_data",
       "text": [
        "Loading required package: RCurl\n"
       ]
      }
     ],
     "prompt_number": 6
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R\n",
      "station <- c(\"9413450\")\n",
      "startdate <- c(\"20080123\")\n",
      "enddate <- c(\"20080124\")"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 11
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R\n",
      "url1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\n",
      "url2 = \"SixMin_Verified_Water_Level.ascii?\"\n",
      "url3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \n",
      "url4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \n",
      "url5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \n",
      "url6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\n",
      "url7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \n",
      "url8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \n",
      "url9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \n",
      "url10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \n",
      "url11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \n",
      "url12 = \"%22\"\n",
      "url13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\""
     ],
     "language": "python",
     "metadata": {},
     "outputs": [],
     "prompt_number": 13
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R\n",
      "station <- c(\"9413450\")\n",
      "startdate <- c(\"20080123\")\n",
      "enddate = <- c(\"20080124\")\n",
      "\n",
      "\n",
      "url1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\n",
      "url2 = \"SixMin_Verified_Water_Level.ascii?\"\n",
      "url3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \n",
      "url4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \n",
      "url5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \n",
      "url6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\n",
      "url7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \n",
      "url8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \n",
      "url9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \n",
      "url10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \n",
      "url11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \n",
      "url12 = \"%22\"\n",
      "url13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"\n",
      "url14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\n",
      "url15 = \"%22\"\n",
      "url16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\n",
      "url17 = \"%22\"\n",
      "####### DON'T CHANGE ANYTHING ABOVE THIS LINE ############\n",
      "##############################################\n",
      "\n",
      "#Assemble the URL\n",
      "urltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\n",
      "\t\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\n",
      "\n",
      "#Download the data\n",
      "#cat(\"Contacting server...\\n\"); flush.console()\n",
      "dat = getURL(urltotal) #use RCurl to retrieve text into a vector 'dat'\n",
      "#cat(\"Data returned...\\n\"); flush.console()\n",
      "Sys.sleep(2) #If you access data in a loop, be courteous and give the server\n",
      "#a short break between requests\n",
      "#cleanup\n",
      "rm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\n",
      "rm(url15,url16,url17)\n",
      "\n",
      "con = textConnection(dat) #create text Connection to dat vector\n",
      "all.lines = readLines(con) #read lines of text into separate slots in a vector\n",
      "close(con) #close connection to dat vector"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "ename": "RInterpreterError",
       "evalue": "Failed to parse and evaluate line u'station <- c(\"9413450\")\\nstartdate <- c(\"20080123\")\\nenddate = <- c(\"20080124\")\\n\\n\\nurl1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\\nurl2 = \"SixMin_Verified_Water_Level.ascii?\"\\nurl3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \\nurl4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \\nurl5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \\nurl6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\\nurl7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \\nurl8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \\nurl9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \\nurl10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \\nurl11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \\nurl12 = \"%22\"\\nurl13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"\\nurl14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\\nurl15 = \"%22\"\\nurl16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\\nurl17 = \"%22\"\\n####### DON\\'T CHANGE ANYTHING ABOVE THIS LINE ############\\n##############################################\\n\\n#Assemble the URL\\nurltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\\n\\t\\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\\n\\n#Download the data\\n#cat(\"Contacting server...\\\\n\"); flush.console()\\ndat = getURL(urltotal) #use RCurl to retrieve text into a vector \\'dat\\'\\n#cat(\"Data returned...\\\\n\"); flush.console()\\nSys.sleep(2) #If you access data in a loop, be courteous and give the server\\n#a short break between requests\\n#cleanup\\nrm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\\nrm(url15,url16,url17)\\n\\ncon = textConnection(dat) #create text Connection to dat vector\\nall.lines = readLines(con) #read lines of text into separate slots in a vector\\nclose(con) #close connection to dat vector'.\nR error message: u'Error while parsing the string.'",
       "output_type": "pyerr",
       "traceback": [
        "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mRInterpreterError\u001b[0m                         Traceback (most recent call last)",
        "\u001b[1;32m<ipython-input-17-9de0d21a54a2>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mget_ipython\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_run_cached_cell_magic\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34mu'R'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34mu''\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/core/interactiveshell.pyc\u001b[0m in \u001b[0;36m_run_cached_cell_magic\u001b[1;34m(self, magic_name, line)\u001b[0m\n\u001b[0;32m   2552\u001b[0m         \u001b[0mcell\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_current_cell_magic_body\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2553\u001b[0m         \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_current_cell_magic_body\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mNone\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2554\u001b[1;33m         \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrun_cell_magic\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmagic_name\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mline\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mcell\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   2555\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2556\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0mrun_cell\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mraw_cell\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstore_history\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mFalse\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msilent\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mFalse\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mshell_futures\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/core/interactiveshell.pyc\u001b[0m in \u001b[0;36mrun_cell_magic\u001b[1;34m(self, magic_name, line, cell)\u001b[0m\n\u001b[0;32m   2132\u001b[0m             \u001b[0mmagic_arg_s\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mvar_expand\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mline\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstack_depth\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2133\u001b[0m             \u001b[1;32mwith\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mbuiltin_trap\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2134\u001b[1;33m                 \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfn\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmagic_arg_s\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mcell\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   2135\u001b[0m             \u001b[1;32mreturn\u001b[0m \u001b[0mresult\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2136\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/extensions/rmagic.py\u001b[0m in \u001b[0;36mR\u001b[1;34m(self, line, cell, local_ns)\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/core/magic.pyc\u001b[0m in \u001b[0;36m<lambda>\u001b[1;34m(f, *a, **k)\u001b[0m\n\u001b[0;32m    189\u001b[0m     \u001b[1;31m# but it's overkill for just that one bit of state.\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    190\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0mmagic_deco\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0marg\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 191\u001b[1;33m         \u001b[0mcall\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;32mlambda\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mk\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mk\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    192\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    193\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mcallable\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0marg\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/extensions/rmagic.py\u001b[0m in \u001b[0;36mR\u001b[1;34m(self, line, cell, local_ns)\u001b[0m\n\u001b[0;32m    550\u001b[0m                 \u001b[0mreturn_output\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mFalse\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    551\u001b[0m         \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 552\u001b[1;33m             \u001b[0mtext_result\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0meval\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcode\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    553\u001b[0m             \u001b[0mtext_output\u001b[0m \u001b[1;33m+=\u001b[0m \u001b[0mtext_result\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    554\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/extensions/rmagic.py\u001b[0m in \u001b[0;36meval\u001b[1;34m(self, line)\u001b[0m\n\u001b[0;32m    159\u001b[0m         \u001b[1;32mexcept\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mri\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mRRuntimeError\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mexception\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    160\u001b[0m             \u001b[0mwarning_or_other_msg\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mflush\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;31m# otherwise next return seems to have copy of error\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 161\u001b[1;33m             \u001b[1;32mraise\u001b[0m \u001b[0mRInterpreterError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mline\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstr_to_unicode\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mstr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexception\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mwarning_or_other_msg\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    162\u001b[0m         \u001b[0mtext_output\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mflush\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    163\u001b[0m         \u001b[0mri\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mset_writeconsole\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mold_writeconsole\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;31mRInterpreterError\u001b[0m: Failed to parse and evaluate line u'station <- c(\"9413450\")\\nstartdate <- c(\"20080123\")\\nenddate = <- c(\"20080124\")\\n\\n\\nurl1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\\nurl2 = \"SixMin_Verified_Water_Level.ascii?\"\\nurl3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \\nurl4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \\nurl5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \\nurl6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\\nurl7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \\nurl8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \\nurl9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \\nurl10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \\nurl11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \\nurl12 = \"%22\"\\nurl13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"\\nurl14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\\nurl15 = \"%22\"\\nurl16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\\nurl17 = \"%22\"\\n####### DON\\'T CHANGE ANYTHING ABOVE THIS LINE ############\\n##############################################\\n\\n#Assemble the URL\\nurltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\\n\\t\\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\\n\\n#Download the data\\n#cat(\"Contacting server...\\\\n\"); flush.console()\\ndat = getURL(urltotal) #use RCurl to retrieve text into a vector \\'dat\\'\\n#cat(\"Data returned...\\\\n\"); flush.console()\\nSys.sleep(2) #If you access data in a loop, be courteous and give the server\\n#a short break between requests\\n#cleanup\\nrm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\\nrm(url15,url16,url17)\\n\\ncon = textConnection(dat) #create text Connection to dat vector\\nall.lines = readLines(con) #read lines of text into separate slots in a vector\\nclose(con) #close connection to dat vector'.\nR error message: u'Error while parsing the string.'"
       ]
      }
     ],
     "prompt_number": 17
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [
      "%%R\n",
      "station <- c(\"9413450\")\n",
      "startdate <- c(\"20080123\")\n",
      "enddate = <- c(\"20080124\")\n",
      "\n",
      "\n",
      "url1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\n",
      "url2 = \"SixMin_Verified_Water_Level.ascii?\"\n",
      "url3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \n",
      "url4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \n",
      "url5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \n",
      "url6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\n",
      "url7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \n",
      "url8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \n",
      "url9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \n",
      "url10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \n",
      "url11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \n",
      "url12 = \"%22\"\n",
      "url13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"\n",
      "url14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\n",
      "url15 = \"%22\"\n",
      "url16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\n",
      "url17 = \"%22\"\n",
      "####### DON'T CHANGE ANYTHING ABOVE THIS LINE ############\n",
      "##############################################\n",
      "\n",
      "#Assemble the URL\n",
      "urltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\n",
      "\t\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\n",
      "\n",
      "#Download the data\n",
      "#cat(\"Contacting server...\\n\"); flush.console()\n",
      "dat = getURL(urltotal) #use RCurl to retrieve text into a vector 'dat'\n",
      "#cat(\"Data returned...\\n\"); flush.console()\n",
      "Sys.sleep(2) #If you access data in a loop, be courteous and give the server\n",
      "#a short break between requests\n",
      "#cleanup\n",
      "rm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\n",
      "rm(url15,url16,url17)\n",
      "\n",
      "con = textConnection(dat) #create text Connection to dat vector\n",
      "all.lines = readLines(con) #read lines of text into separate slots in a vector\n",
      "close(con) #close connection to dat vector\n",
      "\n",
      "if (length(grep('^Error',all.lines))>0) { #check for error in retrieval\n",
      "\tcat(\"There was an error...\\n\")\n",
      "\tcat(dat,\"\\n\") #print contents of dat to show error\n",
      "\tflush.console()\n",
      "} else { #retrieval was successful, parse the text\n",
      "\n",
      "\t#The column headers are typically preceded by a line of dashes\n",
      "\theaderlines = grep(\"^--------\",all.lines) #find index of headers (-1)\n",
      "\n",
      "\t#read column header names into a vector\n",
      "\tcon = textConnection(dat)\n",
      "\theaders = scan(con, skip = headerlines, nlines = 1, sep = \",\",\n",
      "\t\t\twhat = \"character\", strip.white = TRUE)\n",
      "\tclose(con)\n",
      "\n",
      "\t#read rest of the data into a data frame 'df'\n",
      "\tcon = textConnection(dat)\n",
      "\tdf = read.table(con, skip = headerlines+1, sep = \",\", header = FALSE,\n",
      "\t\t\tquote = \"\\\"\", col.names = headers, strip.white = TRUE,\n",
      "\t\t\tstringsAsFactors = FALSE)\n",
      "\tclose(con)\n",
      "\n",
      "\t###########################################################################\n",
      "\t#The following operations will need to be altered if you change the \n",
      "\t#fields or data type being returned by the OPeNDAP server\n",
      "\n",
      "\t#Convert the time column to POSIX time (seconds since 1970-01-01 00:00:00)\n",
      "\tdf[,3] = as.POSIXct(strptime(df[,3],format = \"%b %d %Y %I:%M%p\",\n",
      "\t\t\t\t\ttz = \"GMT\"))\n",
      "\n",
      "\t#Give the columns shorter names\n",
      "\tnames(df) = c(\"stationId\",\"datum\",\"TimeUTC\",\"TideHT\",\"Flag.Inferred\",\n",
      "\t\t\t\"Flag.Flat.Tol\",\"Flag.Rate.Tol\",\"Flag.Temp.Tol\")\n",
      "\n",
      "\t#Uncomment this if you want to plot the data\n",
      "#\tplot(df$TimeUTC, df$TideHT, type = \"l\",\n",
      "#\t\t\txlab = \"Date\",ylab = \"Tide Height, meters\")\n",
      "\n",
      "\t#Save data automatically to a .csv file. \n",
      "\tfilename = paste(\"Station_\",station,\"_\",startdate,\"-\",enddate,\n",
      "\t\t\t\".csv\",sep = \"\") #make file name\n",
      "\twrite.csv(df,filename,row.names = FALSE, quote = FALSE) #write file to disk\n",
      "\tcat(\"Saved to \",getwd(),\"/\",filename,\"\\n\",sep = \"\")\n",
      "#\tflush.console()\n",
      "\n",
      "\t#Alternate file save method lets user specify file name at run time\n",
      "\t#Uncomment this if you wish to use it instead of the automated file \n",
      "\t#output above\n",
      "#\twrite.csv(df,file.choose(),row.names = FALSE, quote = FALSE)\n",
      "\n",
      "\t#cleanup\n",
      "\trm(dat,con,all.lines,startdate,enddate,filename,headerlines, headers)\n",
      "\n",
      "} #end of if-else statement\n"
     ],
     "language": "python",
     "metadata": {},
     "outputs": [
      {
       "ename": "RInterpreterError",
       "evalue": "Failed to parse and evaluate line u'station <- c(\"9413450\")\\nstartdate <- c(\"20080123\")\\nenddate = <- c(\"20080124\")\\n\\n\\nurl1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\\nurl2 = \"SixMin_Verified_Water_Level.ascii?\"\\nurl3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \\nurl4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \\nurl5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \\nurl6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\\nurl7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \\nurl8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \\nurl9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \\nurl10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \\nurl11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \\nurl12 = \"%22\"\\nurl13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"\\nurl14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\\nurl15 = \"%22\"\\nurl16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\\nurl17 = \"%22\"\\n####### DON\\'T CHANGE ANYTHING ABOVE THIS LINE ############\\n##############################################\\n\\n#Assemble the URL\\nurltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\\n\\t\\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\\n\\n#Download the data\\n#cat(\"Contacting server...\\\\n\"); flush.console()\\ndat = getURL(urltotal) #use RCurl to retrieve text into a vector \\'dat\\'\\n#cat(\"Data returned...\\\\n\"); flush.console()\\nSys.sleep(2) #If you access data in a loop, be courteous and give the server\\n#a short break between requests\\n#cleanup\\nrm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\\nrm(url15,url16,url17)\\n\\ncon = textConnection(dat) #create text Connection to dat vector\\nall.lines = readLines(con) #read lines of text into separate slots in a vector\\nclose(con) #close connection to dat vector\\n\\nif (length(grep(\\'^Error\\',all.lines))>0) { #check for error in retrieval\\n\\tcat(\"There was an error...\\\\n\")\\n\\tcat(dat,\"\\\\n\") #print contents of dat to show error\\n\\tflush.console()\\n} else { #retrieval was successful, parse the text\\n\\n\\t#The column headers are typically preceded by a line of dashes\\n\\theaderlines = grep(\"^--------\",all.lines) #find index of headers (-1)\\n\\n\\t#read column header names into a vector\\n\\tcon = textConnection(dat)\\n\\theaders = scan(con, skip = headerlines, nlines = 1, sep = \",\",\\n\\t\\t\\twhat = \"character\", strip.white = TRUE)\\n\\tclose(con)\\n\\n\\t#read rest of the data into a data frame \\'df\\'\\n\\tcon = textConnection(dat)\\n\\tdf = read.table(con, skip = headerlines+1, sep = \",\", header = FALSE,\\n\\t\\t\\tquote = \"\\\\\"\", col.names = headers, strip.white = TRUE,\\n\\t\\t\\tstringsAsFactors = FALSE)\\n\\tclose(con)\\n\\n\\t###########################################################################\\n\\t#The following operations will need to be altered if you change the \\n\\t#fields or data type being returned by the OPeNDAP server\\n\\n\\t#Convert the time column to POSIX time (seconds since 1970-01-01 00:00:00)\\n\\tdf[,3] = as.POSIXct(strptime(df[,3],format = \"%b %d %Y %I:%M%p\",\\n\\t\\t\\t\\t\\ttz = \"GMT\"))\\n\\n\\t#Give the columns shorter names\\n\\tnames(df) = c(\"stationId\",\"datum\",\"TimeUTC\",\"TideHT\",\"Flag.Inferred\",\\n\\t\\t\\t\"Flag.Flat.Tol\",\"Flag.Rate.Tol\",\"Flag.Temp.Tol\")\\n\\n\\t#Uncomment this if you want to plot the data\\n#\\tplot(df$TimeUTC, df$TideHT, type = \"l\",\\n#\\t\\t\\txlab = \"Date\",ylab = \"Tide Height, meters\")\\n\\n\\t#Save data automatically to a .csv file. \\n\\tfilename = paste(\"Station_\",station,\"_\",startdate,\"-\",enddate,\\n\\t\\t\\t\".csv\",sep = \"\") #make file name\\n\\twrite.csv(df,filename,row.names = FALSE, quote = FALSE) #write file to disk\\n\\tcat(\"Saved to \",getwd(),\"/\",filename,\"\\\\n\",sep = \"\")\\n#\\tflush.console()\\n\\n\\t#Alternate file save method lets user specify file name at run time\\n\\t#Uncomment this if you wish to use it instead of the automated file \\n\\t#output above\\n#\\twrite.csv(df,file.choose(),row.names = FALSE, quote = FALSE)\\n\\n\\t#cleanup\\n\\trm(dat,con,all.lines,startdate,enddate,filename,headerlines, headers)\\n\\n} #end of if-else statement\\n'.\nR error message: u'Error while parsing the string.'",
       "output_type": "pyerr",
       "traceback": [
        "\u001b[1;31m---------------------------------------------------------------------------\u001b[0m\n\u001b[1;31mRInterpreterError\u001b[0m                         Traceback (most recent call last)",
        "\u001b[1;32m<ipython-input-9-9de0d21a54a2>\u001b[0m in \u001b[0;36m<module>\u001b[1;34m()\u001b[0m\n\u001b[1;32m----> 1\u001b[1;33m \u001b[0mget_ipython\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_run_cached_cell_magic\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;34mu'R'\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;34mu''\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/core/interactiveshell.pyc\u001b[0m in \u001b[0;36m_run_cached_cell_magic\u001b[1;34m(self, magic_name, line)\u001b[0m\n\u001b[0;32m   2552\u001b[0m         \u001b[0mcell\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_current_cell_magic_body\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2553\u001b[0m         \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0m_current_cell_magic_body\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mNone\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2554\u001b[1;33m         \u001b[1;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mrun_cell_magic\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmagic_name\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mline\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mcell\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   2555\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2556\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0mrun_cell\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mself\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mraw_cell\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstore_history\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mFalse\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0msilent\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mFalse\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mshell_futures\u001b[0m\u001b[1;33m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/core/interactiveshell.pyc\u001b[0m in \u001b[0;36mrun_cell_magic\u001b[1;34m(self, magic_name, line, cell)\u001b[0m\n\u001b[0;32m   2132\u001b[0m             \u001b[0mmagic_arg_s\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mvar_expand\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mline\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstack_depth\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2133\u001b[0m             \u001b[1;32mwith\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mbuiltin_trap\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m-> 2134\u001b[1;33m                 \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mfn\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mmagic_arg_s\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mcell\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m   2135\u001b[0m             \u001b[1;32mreturn\u001b[0m \u001b[0mresult\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m   2136\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/extensions/rmagic.py\u001b[0m in \u001b[0;36mR\u001b[1;34m(self, line, cell, local_ns)\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/core/magic.pyc\u001b[0m in \u001b[0;36m<lambda>\u001b[1;34m(f, *a, **k)\u001b[0m\n\u001b[0;32m    189\u001b[0m     \u001b[1;31m# but it's overkill for just that one bit of state.\u001b[0m\u001b[1;33m\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    190\u001b[0m     \u001b[1;32mdef\u001b[0m \u001b[0mmagic_deco\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0marg\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 191\u001b[1;33m         \u001b[0mcall\u001b[0m \u001b[1;33m=\u001b[0m \u001b[1;32mlambda\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m*\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mk\u001b[0m\u001b[1;33m:\u001b[0m \u001b[0mf\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m*\u001b[0m\u001b[0ma\u001b[0m\u001b[1;33m,\u001b[0m \u001b[1;33m**\u001b[0m\u001b[0mk\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    192\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    193\u001b[0m         \u001b[1;32mif\u001b[0m \u001b[0mcallable\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0marg\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/extensions/rmagic.py\u001b[0m in \u001b[0;36mR\u001b[1;34m(self, line, cell, local_ns)\u001b[0m\n\u001b[0;32m    550\u001b[0m                 \u001b[0mreturn_output\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mFalse\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    551\u001b[0m         \u001b[1;32melse\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 552\u001b[1;33m             \u001b[0mtext_result\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mresult\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0meval\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mcode\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    553\u001b[0m             \u001b[0mtext_output\u001b[0m \u001b[1;33m+=\u001b[0m \u001b[0mtext_result\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    554\u001b[0m \u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;32m/home/rsignell/epd-7.2-1/lib/python2.7/site-packages/IPython/extensions/rmagic.py\u001b[0m in \u001b[0;36meval\u001b[1;34m(self, line)\u001b[0m\n\u001b[0;32m    159\u001b[0m         \u001b[1;32mexcept\u001b[0m \u001b[1;33m(\u001b[0m\u001b[0mri\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mRRuntimeError\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mValueError\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;32mas\u001b[0m \u001b[0mexception\u001b[0m\u001b[1;33m:\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    160\u001b[0m             \u001b[0mwarning_or_other_msg\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mflush\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m \u001b[1;31m# otherwise next return seems to have copy of error\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[1;32m--> 161\u001b[1;33m             \u001b[1;32mraise\u001b[0m \u001b[0mRInterpreterError\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mline\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mstr_to_unicode\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mstr\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mexception\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m,\u001b[0m \u001b[0mwarning_or_other_msg\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0m\u001b[0;32m    162\u001b[0m         \u001b[0mtext_output\u001b[0m \u001b[1;33m=\u001b[0m \u001b[0mself\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mflush\u001b[0m\u001b[1;33m(\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n\u001b[0;32m    163\u001b[0m         \u001b[0mri\u001b[0m\u001b[1;33m.\u001b[0m\u001b[0mset_writeconsole\u001b[0m\u001b[1;33m(\u001b[0m\u001b[0mold_writeconsole\u001b[0m\u001b[1;33m)\u001b[0m\u001b[1;33m\u001b[0m\u001b[0m\n",
        "\u001b[1;31mRInterpreterError\u001b[0m: Failed to parse and evaluate line u'station <- c(\"9413450\")\\nstartdate <- c(\"20080123\")\\nenddate = <- c(\"20080124\")\\n\\n\\nurl1 = \"http://opendap.co-ops.nos.noaa.gov/dods/IOOS/\"\\nurl2 = \"SixMin_Verified_Water_Level.ascii?\"\\nurl3 = \"WATERLEVEL_6MIN_VFD_PX._STATION_ID,\" \\nurl4 = \"WATERLEVEL_6MIN_VFD_PX._DATUM,\" \\nurl5 = \"WATERLEVEL_6MIN_VFD_PX.DATE_TIME,\" \\nurl6 = \"WATERLEVEL_6MIN_VFD_PX.WL_VALUE,\"\\nurl7 = \"WATERLEVEL_6MIN_VFD_PX.I,\" \\nurl8 = \"WATERLEVEL_6MIN_VFD_PX.F,\" \\nurl9 = \"WATERLEVEL_6MIN_VFD_PX.R,\" \\nurl10 = \"WATERLEVEL_6MIN_VFD_PX.T\" \\nurl11 = \"&WATERLEVEL_6MIN_VFD_PX._STATION_ID=%22\" \\nurl12 = \"%22\"\\nurl13 = \"&WATERLEVEL_6MIN_VFD_PX._DATUM=%22MLLW%22\"\\nurl14 = \"&WATERLEVEL_6MIN_VFD_PX._BEGIN_DATE=%22\" # start date gets put in here\\nurl15 = \"%22\"\\nurl16 = \"&WATERLEVEL_6MIN_VFD_PX._END_DATE=%22\" # end date gets put in here\\nurl17 = \"%22\"\\n####### DON\\'T CHANGE ANYTHING ABOVE THIS LINE ############\\n##############################################\\n\\n#Assemble the URL\\nurltotal = paste(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,\\n\\t\\tstation,url12,url13,url14,startdate,url15,url16,enddate,url17,sep =\"\")\\n\\n#Download the data\\n#cat(\"Contacting server...\\\\n\"); flush.console()\\ndat = getURL(urltotal) #use RCurl to retrieve text into a vector \\'dat\\'\\n#cat(\"Data returned...\\\\n\"); flush.console()\\nSys.sleep(2) #If you access data in a loop, be courteous and give the server\\n#a short break between requests\\n#cleanup\\nrm(url1,url2,url3,url4,url5,url6,url7,url8,url9,url10,url11,url12,url13,url14)\\nrm(url15,url16,url17)\\n\\ncon = textConnection(dat) #create text Connection to dat vector\\nall.lines = readLines(con) #read lines of text into separate slots in a vector\\nclose(con) #close connection to dat vector\\n\\nif (length(grep(\\'^Error\\',all.lines))>0) { #check for error in retrieval\\n\\tcat(\"There was an error...\\\\n\")\\n\\tcat(dat,\"\\\\n\") #print contents of dat to show error\\n\\tflush.console()\\n} else { #retrieval was successful, parse the text\\n\\n\\t#The column headers are typically preceded by a line of dashes\\n\\theaderlines = grep(\"^--------\",all.lines) #find index of headers (-1)\\n\\n\\t#read column header names into a vector\\n\\tcon = textConnection(dat)\\n\\theaders = scan(con, skip = headerlines, nlines = 1, sep = \",\",\\n\\t\\t\\twhat = \"character\", strip.white = TRUE)\\n\\tclose(con)\\n\\n\\t#read rest of the data into a data frame \\'df\\'\\n\\tcon = textConnection(dat)\\n\\tdf = read.table(con, skip = headerlines+1, sep = \",\", header = FALSE,\\n\\t\\t\\tquote = \"\\\\\"\", col.names = headers, strip.white = TRUE,\\n\\t\\t\\tstringsAsFactors = FALSE)\\n\\tclose(con)\\n\\n\\t###########################################################################\\n\\t#The following operations will need to be altered if you change the \\n\\t#fields or data type being returned by the OPeNDAP server\\n\\n\\t#Convert the time column to POSIX time (seconds since 1970-01-01 00:00:00)\\n\\tdf[,3] = as.POSIXct(strptime(df[,3],format = \"%b %d %Y %I:%M%p\",\\n\\t\\t\\t\\t\\ttz = \"GMT\"))\\n\\n\\t#Give the columns shorter names\\n\\tnames(df) = c(\"stationId\",\"datum\",\"TimeUTC\",\"TideHT\",\"Flag.Inferred\",\\n\\t\\t\\t\"Flag.Flat.Tol\",\"Flag.Rate.Tol\",\"Flag.Temp.Tol\")\\n\\n\\t#Uncomment this if you want to plot the data\\n#\\tplot(df$TimeUTC, df$TideHT, type = \"l\",\\n#\\t\\t\\txlab = \"Date\",ylab = \"Tide Height, meters\")\\n\\n\\t#Save data automatically to a .csv file. \\n\\tfilename = paste(\"Station_\",station,\"_\",startdate,\"-\",enddate,\\n\\t\\t\\t\".csv\",sep = \"\") #make file name\\n\\twrite.csv(df,filename,row.names = FALSE, quote = FALSE) #write file to disk\\n\\tcat(\"Saved to \",getwd(),\"/\",filename,\"\\\\n\",sep = \"\")\\n#\\tflush.console()\\n\\n\\t#Alternate file save method lets user specify file name at run time\\n\\t#Uncomment this if you wish to use it instead of the automated file \\n\\t#output above\\n#\\twrite.csv(df,file.choose(),row.names = FALSE, quote = FALSE)\\n\\n\\t#cleanup\\n\\trm(dat,con,all.lines,startdate,enddate,filename,headerlines, headers)\\n\\n} #end of if-else statement\\n'.\nR error message: u'Error while parsing the string.'"
       ]
      }
     ],
     "prompt_number": 9
    },
    {
     "cell_type": "code",
     "collapsed": false,
     "input": [],
     "language": "python",
     "metadata": {},
     "outputs": []
    }
   ],
   "metadata": {}
  }
 ]
}