{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Dfs3 - Basic"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import mikeio"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<mikeio.Dataset>\n",
       "dims: (time:1, z:17, y:112, x:91)\n",
       "time: 2001-12-28 00:00:00 (time-invariant)\n",
       "geometry: Grid3D(nz=17, ny=112, nx=91)\n",
       "items:\n",
       "  0:  Diss. oxygen (mg/l) <Concentration 3> (mg per liter)"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds = mikeio.read(\"../tests/testdata/dissolved_oxygen.dfs3\")\n",
    "ds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<mikeio.Grid3D>\n",
       "x: [0, 150, ..., 1.35e+04] (nx=91, dx=150)\n",
       "y: [0, 150, ..., 1.665e+04] (ny=112, dy=150)\n",
       "z: [0, 1, ..., 16] (nz=17, dz=1)\n",
       "origin: (10.37, 55.42), orientation: 18.125\n",
       "projection: PROJCS[\"UTM-32\",GEOGCS[\"Unused\",DATUM[\"UTM Projections\",SPHEROID[\"WGS 1984\",6378137,298.257223563]],PRIMEM[\"Greenwich\",0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",500000],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",9],PARAMETER[\"Scale_Factor\",0.9996],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1]]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds.geometry"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<mikeio.DataArray>\n",
       "name: Diss. oxygen (mg/l)\n",
       "dims: (time:1, z:17, y:112, x:91)\n",
       "time: 2001-12-28 00:00:00 (time-invariant)\n",
       "geometry: Grid3D(nz=17, ny=112, nx=91)"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "do  = ds[0]\n",
    "do"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA9g0lEQVR4nO29ebhdRZX3//kmYXIAgaAdCXRAQV+kHSCmGRRQUHAiqCChRYLSTUvjPAHSbzs1/cMJleaV7rTMjQwiCtogM0RmAw5MolEQIkiIINAOIbl3/f7YdXLq3Oy6p/a9Z5+z773r8zz7uXvXrqpdZ7hn1apVay2ZGY7jOI6Tw7RBD8BxHMeZOLjQcBzHcbJxoeE4juNk40LDcRzHycaFhuM4jpONCw3HcRwnGxcajuM4EwBJW0i6RtI9ku6S9MFQvomkKyT9MvzduNZxuJ+G4zhO85E0C5hlZrdLejZwG7AfcCjwmJkdL+loYGMzO6qucbim4TiOMwEws4fN7PZw/hRwD7A5MB84I1Q7g0KQ1MaU0zRmzpxpc+bMqa3/X9z269r6dpzJxLY7bj3oIXTltttuW2Fmm42nj71f80z7/WNDec/72cq7gL9ERYvMbNHIepLmAIuB7YEHzOw50b3Hzay2JaoZdXXcVObMmcOSJUtq6/910w6orW/HmUxcseRbgx5CVyT9Zrx9rHhsiFsum51Vd51Zv/qLmc3tMqZnAd8GPmRmT0oa7xArMeWERt3MeNEL2xdDw6NXjj/sWOOb3l41tGnlK4gajvqOnxOXT5/e7mf9dUYfy8gxpGiNOaduXCfxViiuE7+O6D1IPit+/+K2cf1U/xadrxPem+HEINctf+8sfn5qoTeukziP+1Fq7KvbY9PqaNYajzn6rtiM6aOPKyb+jqW+k9G5orHw9Kr2+cqny8c11B6vxfWnDMaQdfktyETSOhQC42wzuzAUPyJplpk9HOwey3vysARu03Acx6kRA4axrGM0VKgUpwD3mNkJ0a2LgYXhfCFwUR2vo4VrGj1m+BnrrjnX0OhfAstRK2OxHs9Oo747ZqeJ2bKt09Y6SD23in0ro27H659WPms1orFMT8zKUxpLPEGOL1Ia3rRYMykbcPkzbXri/Uq9jxnlqc/eEtqI1onaxp9lgjX9JN7T5HcvMY2MP0tbN/ruxWPZIPHdjzSjWEva5+X/0u5zRvvBly35dPkgJjDDKVW7GrsC7wLukPSTUPZJ4HjgfEmHAQ8Ata6Ru9BwHMepEcNY1YPlKTO7HkjNNPcc9wMycaHRa+J15Xh2UcVYFc+s4/7imaeimV/H2nP0nPh7mjPLTJSrTKvIeD2m2J6Q0IxSGsu0hE0jngnH2kA8nBnt2W/Hs1ZHdYbydrOsRVWjY877nkGqrTreg7W1tpRNLMsGEz+/Q7lJaISRlhv/RnZoI3GdqDwe++t3/lxicAWX3/R/R73fNAwY6rL0NJGozaYh6VRJyyXdOaL8/ZLuDR6NX4jKj5G0NNzbOyrfUdId4d6JYV0PSetJOi+U3xK2oDmO4zSOXtg0mkKdmsbpwEnAma0CSa+hcER5qZmtlPTcUL4dsAB4CfB84EpJ25rZEHAycDhwM3AJsA9wKXAY8LiZvVDSAuDzwIE1vp4sLvvRpwY9hEq8ft5ny2/EE/3hsc2QO1pFk80O5SLqO94Rllzbz9FSIjq0sxldbAEdu6ESO4pS9RPl49EucrCUzaI1/ooaRUffGZ+BzYjtRFGdabEmHD+r/RkMxxpIbG/q+IzXHtfrdvnXNedX3PjPidE3BwOGJpE/XG2ahpktBh4bUXwEcLyZrQx1WlvD5gPnmtlKM7sPWArMC9vHNjSzm6xYgzmTtrdj7AV5AbBnSwtxHMdpEsOZx0Sg3zaNbYFXSzqOwuvxY2b2IwpX+JujestC2apwPrKc8PdBADNbLekJYFNgxciHSjqcQlthyy237OXrmfBcfuu/lJa/fqdIAxnj1KJDQ0lpC7FtJn5Q4pkdu60U2S6U8NOIZ9+RJqNY9SnbbZXjG5JTHpN6H6v+WiR21HVqSms3u/zm8s86Jp7Fp0jaVzK0yQ4NdkZi7PHH12UH4uteddya8yuuP3bUuoPCMLdpjIMZwMbATsDHKbaJifIdATZKOV3udRaaLTKzuWY2d7PNxhURwHEcpxJmsCrzmAj0W9NYBlwYlppuVTE9nBnKt4jqzQYeCuWzS8qJ2iyTNAPYiLWXw5wu7HLgl0vLb8yYlZaRo6HEGkjK7tFBcnZfrqWU7vaCETvR1tZGOrzsldi9VdUWUPe0LGE/ydEqyohtBK/btbvWEWOpXVUdmhGldWKzhzr6qTSEhiKGJscLAfqvaXwXeC2ApG2BdSmWky4GFoQdUVsB2wC3mtnDwFOSdgoaySG0vR1jL8j9gattqkVfdByn8RgwbHnHRKA2TUPSOcAewExJy4BPAacCp4ZtuE8DC8MP/V2SzgfupthNf2TYOQWF8fx0YAOKXVOXhvJTgLMkLaXQMBbU9VomG7u8o1y7iCdDOx+UqNOqmliHT2koKQ0ky+4RkdzRExkGklOHxE6tlrd60qaS8nVI7ZLKmYrVYPVM2afGyhU3lO9MSmogOZPpintVrOS97DBfTZBASJNJ06hNaJjZQYlbByfqHwccV1K+hCL878jyv1Czu7zjOM54KZz7XGg4k5yxzuZSWkyOjSTWRqziP1mWTSO2gUTag1r+BYn4WJ19ZMRsSsXNqoGB+AXlaAsV7UAdTcvsG7FvyATRLloYsGqiDXoUXGg4juPUiCGGJlFAcRcaU5Abz//omvPU7qmYUvtF98l8B7seED0nMfu+oaI20kGOB3vKq7llYKjq+Z7SLhLP7JXWcdltn+lJP2Ml9oeI/SQ6J9O9WY5paR2a4Ms7wzaxxx/jQsNxHKdG3KbhTCpuPK+a1tHSMDo0ijirWyLgbkoDieu86u1f6vr46yvuzsrapaQSm0YO44laU3GMl/14sNpFispL9WN8z1LfpSuv++SY+usvYshtGo7jOE4OReY+FxrOJCRH61gzy4u1i+GS+2SuaqcqJcpT2khKA4lJRvQtm/1WtT/kZEPMsW9U1ZIGzJWLq8V72mv3f+teqcvnMdHMA2biaeuebTEHSacCbwaWm9n2oexlwH8AzwLuB95pZk/25IElTB7x5ziO01CGUdaRwekU6SFivgEcbWZ/A3yHIq5fbbimMcWpYseI6dAoEvaNHFK2jlRcohwN5Ppvf6y0TjeP6b3nfnrU+yPHlbU+X1W7iGiqHWM85Ngg9nzN/7dW2VXXTgTbRTmFIbw383MzW1yScO5FwOJwfgVwGVBbekPXNBzHcWqlMITnHBRhl5ZEx+EZD7gT2DecH0Bn8Nee45rGFCQn9lQvSGoROVjCONKRna68aY7WUcZlSz5dWr73jt29rnNyeFf2pJ6iXHXNMYMeQk+paAhfYWZzKz7iPcCJkv6FIpDr02WVJJ2Y0deTZjZqOkQXGo7jODUzVKP13sx+Drwe1kQPf1Oi6nyg246RowEXGk4nVT3Cy0jZGZLaRUUP8tS+fJsW79rK0DreFmkdF+ZrHTGxB/ber2yf5+T/tlTk3pStwzMWTzoMscrq+6mV9FwzWy5pGsUP/n8kqn7FzM5I3Gv1tXG357nQmOJUde4rDetQktCoqFwuEXKWrZJ1OgIpJgRI4nf31W9tC5AffmeMAiQRIDC5nTcjvElHqS9PTTp6aQhPpJx4lqQjQ5ULgdNKx2H2VUnTgQ+Y2VdSdbqNwYWG4zhOjRjq2fLUKCknvpbZfkjSfKBUaOTgQsNZQ6x1pCgzoneGsi6nQ3MYR4oyTYs1mfg00joyLPq7zf/imvPFF41/W3tqO29W+ttovBpql8dLYQMJge70jIZ5hN8g6STgPOCPrUIzuz2ncW2vRNKpkpaHLH0j731MkkmaGZUdI2mppHsl7R2V7yjpjnDvxJD2lZAa9rxQfkvJ3mXHcZyBY0aVLbf9YBfgJcBngS+Ho3vgt0CdmsbpwEnAmXGhpC2A1wEPRGXbUaRrfQnwfOBKSduGlK8nA4cDNwOXUHhDXgocBjxuZi+UtAD4PHBgja9nQpMTDDBFhybRxVCbo3X0jMSuXEtZ3RPstm+hdSy+uPeOtJenAizu/Lk1550G9bbRRkPt1xE7Hqa2BveC5LgS257jMOlOOYUhvDdhRHqBmb1mPO1rE21mtpgid/dIvgJ8gs79NPOBc81spZndBywF5kmaBWxoZjeFXOJnAvtFbVo7AS4A9lS3XzTHcZwBMMS0rKMfSHqepFMkXRqut5N0WG77vto0JO0L/NbMfjri931zCk2ixbJQtiqcjyxvtXkQwMxWS3oC2BRYUc/oJzjJcB19HcXa5CgFHRPxyHaRsG/kvKay3Vm7v+kL5ZUTY7zukk90f1CCy29qR3noSGSk8oiFWl2udXTQEXql5E1I/SZVDIdi09vne+7xb6V1cj6D1Jblq686unvjCYShpiVhOp1ih1VLTfwFhX3jlJzGfVtEk/QMikGW6etl76iNUj5am7JnH95yy3/00Udzhus4jtMzmqRpADPN7HzCrMTMVgNDozdp009N4wXAVkBLy5gN3C5pHoUGEcdLmQ08FMpnl5QTtVkmaQawEeXLYZjZImARwNy5c6fMRvjYsS1FHGYjZfdQ2S6lfq4EJn02En4gvX5O3HdUZ483lGsm115aTQOJ7QJ77XZcdCf6EYnDzw8nYqbH4erLbsd+IvHvU45zoXV8CbrXT5DjEPnaPY9vX3T5Hb36iuZrJQYMNysJ0x8lbUr4JCXtBDyR27hvr8TM7jCz55rZHDObQ/Gjv4OZ/Y4iXsqCsCNqK2Ab4FYzexh4StJOwV5xCHBR6PJiYGE43x+4Otg9HMdxGoQYyjz6xEcofj9fIOkGClvxB3Ib16ZplHkumlnpmpmZ3SXpfOBuYDVwZNg5BXAExRrcBhS7pi4N5acAZ0laSqFhLKjppTSeDg0hI1xHR4rVRJiNrrutmiCfq+ZJqnPI0fuxxz6fL61y7Q+O6t5Pyi6QSs7U7XNIzezjPnKmjj3SLNWRTCmjz9Y4E2OMtZKm2kIMGrV7CrgL2J0ipLqAe6mgQNQmNEbxXGzdnzPi+jjguJJ6S4DtS8r/QhEG2HEcp7GYqWnLUzeZ2Q4UwgMASbcDO+Q0do/wyUC8JJ3hbZ3ymI5jM6U05W5zw86ETF2HMi56rjnUrD29Zu+2BnLNZW2tI54td+4Ii2wUvdIYyuqmdl2lvgPxMzs2sCW2sMV1xqqxpF5ndP6a10fv7+UZWl0f6aPjXhJJf0Wx63QDSa+g/clsCDwjtx8XGo7jODVS5NNoxJbbvYFDKTYUnRCVPwVkp0Z0oTEJSM24U7P+ODpsss/El3zQloxa7BI1ahip8cZaR4drRuQFrpTtIjXejtn4GHe5VUx4NRASWoeaYGcrRY3QNEJY9DMkvd3Mvj3WflxoOI7j1Eix5Xbwmoakg83sv4E5kj4y8r6ZnVDSbC1caEwCpkXewsPRJo1xxYHqlgujYx080UdnQKgx12mqdlF5XLHtKbZXdGgXcXnsEV4ek6qj+zKtIsMHo9OOUe75nfWb1y/fncFP2ivRoNhTzwx/nzWeTlxoOI7j1EwTQqOb2X+Gv5/pVnc0XGhMAq77Xnl01lfv184ZkUyf2oPJYU7K1sRmmsp9DpqscSU0ik6v7nKNYtqqcu0ix6bRqbSt/S53+EhE3uGKXlRn1+UaSBapnVdV/TRaZPzmjnlnVs0UodF75OcinQq8GVhuZtuHspdTpHhdn8LP7Z/M7NaStieOPk7LcvBzoeE4jlMzPbRpnM7aKSe+AHzGzC6V9MZwvUdJ29t6MQAXGpOYH363uwaSIkt76AUds/Ie990jkuNK2WBiG0WHptE+j+1QKdsFCS0li1b9qkvpHTvuehNvahA0yWejiHLbm+UpM1tcknDOKHwtoIjB9xAlhN1TXZH072b2/tR9FxqO4zg1UoQRyRYaMyUtia4XhYCro/Eh4DJJX6JYyNul8iA72XW0my40piCxBpKjdcSs8R2I18GT69fl50ntokH77HO0iyx7RYftoly7iH0zUraLXvggdKz5J3NrVOy0X0pHhud7yl4Sax2dDdqnsYd+76mkaawws7kVH3AE8GEz+7akd1DE5durYh/ZDN6k7ziOM8kZRlnHGFkIXBjOvwXM68mgE7imMcVJ2T12mx9pIGWaQVwWtUvt0urwdM7QLqraN2r1ncrRLnI8tpOvNcPDu6NBuS9Ft7p0+F0kYkPNmBbVSfTTo/c6uZOq21R2HFrHIOjl7qkED1FErb0WeC3wy3H2N+pgXWg4juPUTK8M4WUpJ4B/AL4WktH9BTh8lPbTgePNrHy2WPC10cbgQsMZGwmv5KRzeEXbRY6dpNeTt55F6I1ntqnXmmHHiEnOljviQ7V/mFozbeu4n9BQEnU6n1NenDXGuJvU+1FFu6iBVATiXtDLHOGjpJzYMbP9kKQdJSmVtM7MTh+tDxcajuM4NWLA6gYELIz4MXCRpG8Bf2wVmtmF6SZtXGg4tZG0S+TMxBMMIu5bzuy4w8M6adNI9FnVHyN+bqRddNgsymKETU9oFBXtFeOxEQzavjAoGpaEaRPg9xT2jxZG25g+KnWmey1zd/8i8BbgaeBXwLvN7A/h3jHAYcAQ8AEzuyyU70g73eslwAfNzCStR+EVuSPFG3Cgmd1f1+txHMcZE9a75aleYGbvHk/7OsXf6cA+I8quALY3s5cCvwCOAZC0HUWO75eENl8PBhuAkykMO9uEo9XnYcDjZvZC4CtAYjO2MxYWX/TxNQdizWHhiMtSmDKOaVpzoJqPBLKgFcVHrzBbc2h4eM0RlzNM+0gxLTpSj5Lax7RpxTE9OpLvNdHnq65HFtMqHmPtL+O9GDRG7VtuKyFpW0lXSbozXL9U0j/ntq9NaJjZYuCxEWWXm9nqcHkzRQYpgPnAuWa20szuA5YC8yTNAjY0s5uC0eZMYL+oTcst/gJgTyVzYjqO4wyO4aBtdDv6xH9RTNhXAZjZzygm7VkM0qbxHuC8cL45hRBpsSyUrQrnI8tbbR4EMLPVkp4ANgVW1Dhmp/W9jtfE47zkOXkzUtTtEJ6wpbT+V1NDLIsYC3le2h3Z91KZ+DL6jG0mpHwz4pfXqpPIldGM7KNtqmgETfLByMFoRhKmiGeY2a0j5tirU5VHMhChIelYikGe3SoqqWajlI/Wpux5hxP2Lm+55ZaVxuoUS1Ut1oQdiXXU6Ms3nqCD8Q7ApiZeqvycHMGSqp/zg9hlraDqb1XWD3KP1ifG+oNfuV2H92l5lTrDiBhi9XCjDOErJL2A8G5I2h94OLdx34WGpIUUBvI9o33Cy4AtomqzKbwcl9FeworL4zbLglPLRoxYDmsRAn4tApg7d25zAhw5jjMl6Je9IpMjKX4PXyzpt8B9wDtzG/dVaEjaBzgK2N3M/hTduhj4pqQTgOdTGLxvDY4oT0naCbgFOAT496jNQuAmYH/g6pSzitN7Uqlkx6OFK+rJUs6DE+ETTi0JxZPNaAumDUeh0cfzrAlAlpYwsV5Sd6xZy1Nm9mtgL0nPBKaZ2VNV2te55bbM3f0YYD3girCedrOZvdfM7pJ0PnA3xbLVkWY2FLo6gvaW20vDAUUkx7MkLaXQMLINOY7jOP2iaTYNSZtS/B6/CjBJ1wOfNbPf57SvTWgk3N1PGaX+ccBxJeVLgO1Lyv8CHDCeMTrVKQtw+Oq3fmnsHXas4UensbE3MrR35AKaCFpHRIcjXmQVl8UaVqJxHEgwPp8en0f9tE5TgQkTM/7xhGAfV58VflMr//4mNkBc+4P+JWdqktAAzgUWA28P1++k2JSUFU7dPcIdx3FqxBBDzTKEb2Jmn4uu/1XSfrmNXWg44+aH3/nYmvNXva2a1qEO9SKhUVB+kdzdW0P4kq502Cvap8Mz2s+ctjquEjWI26Z+W2ItItIu4v57Yd+oumMqrj+8TvmHMy21mTPxqNJJecZrS07me2V0GwcNM4RfI2kBcH643h/4n9zGLjQcx3FqxBpmCAf+EfgIcFa4ng78UdJHADOzDZMtcaHh9JjrL4y0jrd31zp+GNWP2f0tURKolCNhYv99MhFUTInWkWyXfE7snxJpSZFW0NF0OuXlkY0i6dyXsF0kKakyHi0iRyt4esPppXXW/d+2/Wba0+3zZJKnLmT9/qbq5Lx3NWANEhpm9uzR7kt6iZndlbrfqIU2x3GcyUdeCJEcbUTSqZKWt+JGhbLzJP0kHPdL+sk4B3zWaDdd03Bq4/pvl2sROVz3vfYurd3f9IX2jTq1jpLQImu1y9A6Oogn34lEVB3+Kamo6mOc3uVoF8nUrwlNoLO8fbryOe1OV6/fLl/9h3alDaJNndNWpRJOlQ95ItNDTeN04CSKOHyhbzuwdS7py8AT43zGqIN1oeE4jlMjZjA03BuhYWaLJc0puxcCtr6DzjwZY3rMaDddaDiN57r/+cSa8yxbR4KuWkdiVm7R/1BS60g9M64U/7cpmupH3u+K/DRK/S7GQw3aRTzGlfu0J7h/fPhZa87XW95Wt6avap+v93jHdrJSrrvkE+U3Aru/MdJCM96jWIPtJxV2T82UtCS6XhTCIOXwauARM/tlpcFVxIWG4zhOjRiVlqdWmNncMT7qIOCcMbaNeXq0my40nAlFlo2iYj9V+qscNTZuG19Fto6OXVjTUkaNRHmlXUe90S46Hh8Na6Oz25ty7jjxP9ecb33Fe9acr3pwvTXn6z7Z7nTxxWPTAGJNpMP2FTOtF2raeKg/V0YI2vo2ikym3eruUFL8BPAbM1ttZjuN1t6FhuM4Ts30IZTqXsDPzWxZ15rwdWAH4GcU04Htw/mmkt5rZpeP1tiFhtN4dpsf2TESHuG9mMeld0yNvfcOe0hHebl20VFuiVeY1DpKyvq4qf7G8z8aXbXP739XVByf95jY9tU0erV7qiwQrJmdQhGwNXdp6n7gsJYvRki3/XHgc8CFgAsNx3GcQVHsnuqN9E4EgsXMDq3QzYtj5z0zu1vSK8zs1zkZs11oOI1kTYZA0lpEnVFue7YErYRq1JGtL+F5HdePUsV2pJDtBSk7htMzGpbp515JJ1NEuwU4EPiFpPUIecNHw4WG4zhOzTQpjAhwKPBPwIcopinXAx+jEBiv6dbYhYZTG7se8OX2RTTVuuGCck/xVF6OyjumUrP1nPo10mmhiGb3Fe0eOa9prPm3O0h1MfDdSBMLQ40SGmb2Z+DL4RjJ/3Zr77GnHMdxasYyj34g6c2SfizpMUlPhpTaT+a2rzPd66nAm4HlZrZ9KNuEIkPUHAoL/jvM7PFw7xjgMGAI+ICZXRbKd6Sd7vUS4INmZmH97UyKfcm/Bw40s/vrej1OHh3aRYKO6LeJ8FExHdpFxzp/bzPCJenVzDqVgTBD68iKc+U0EwPrURiRHvFVCp+OO8yqW1vq1DROB/YZUXY0cJWZbQNcFa5bW74WAC8Jbb4uqeX+dDJwOLBNOFp9HgY8bmYvBL4CfL62V+I4jjMOzJR19IkHgTvHIjCg3hzhZYG15lPsMQY4A7gWOCqUn2tmK4H7JC0F5km6H9jQzG4CkHQmsB9waWjz6dDXBcBJkjTWN8LpDb3y2M5iAGvr4/m/ViJWVk7e86TXeK93UqVIvNeDiuU00WjYr9IngEskXQesbBWa2Qk5jfttCH+emT0MYGYPS3puKN8cuDmqtyyUrQrnI8tbbR4Mfa2W9ASwKbBi5EMlHU6hrbDlllv27MU4juN0o2LsqX5wHIXBe31g3aqNm7J7quwdtVHKR2uzdmERJXIRwNy5c5sl86c4nfaKRDTZBDlaTa3/qzl9p9wx4pn7UGzYGceuqljr6LLwHI8lGWPKd0n1BqNpDjCbmNnrx9q437unHpE0CyD8XR7KlwFbRPVmAw+F8tkl5R1tQrCujYDHahu54zjOGDHLO/rElZImjNC4GFgYzhcCF0XlCyStJ2krCoP3rWEp6ylJO4UEI4eMaNPqa3/gardnNAC1D4sOpPbRo/5jTOVHsm3ZUaVu7suo+qzxMK19mLTW0TmG6POIy6dpzZH1njoZCBvOO/rEkcAPJP25aVtu1wqsBRwPnC/pMOAB4AAAM7tL0vnA3cBq4Egza6WjOYL2lttLwwFwCnBWMJo/RrH7ynEcp3k0aDprZs/uXitNnbunSgNrAXsm6h9HYaAZWb6EInTvyPK/EISOM8FIxWCqSq8nZn2cRSfDTXWUV7NvlJKwXcQ5x+PseynGmu/CofDTaICKJumvzOx3463jHuGO4zh10wyX8Et6Uacpu6ecScKN57XzKOxyYNs7PJ5o3fCtduypOJptTN0Ts174kKTGWLnvHP+NHK/xEpJZ+aaXqzqLL3KNoh4Gr2kAL+tiuxDQ1bbhQsOpjViApPjhd9s/Uh0CJMdgPo6lrbIf/HgsHYmfurQbE/FSUVw8DgFS/pyobryu4IKiv/TLCXMUzGx691rd8eUpx3GcOmn5aeQcXZB0qqTlku4cUf5+SfdKuktSIll6b3BNw5m4pLSRMWognZrOmLooHl912SpH6xhKNE68B2vGkHDWa4BddkrRQ2eA04GTKIK1AiDpNRRhlV5qZiujSBu14JqG4zhO3fTIEG5mi1nbifkI4PgQuw8zW75Wwx7imobTU3bdvzyRUkwq1Pn1kU3hVW/r3k+yzx5rIL2iajDHOBjgHm/ovuJQFg4kZfAet5OlU4181W6mpCXR9aIQBmk0tgVeLek44C/Ax8zsR6M1kLQx8Hzgz8D9ZpZtdXGh4TiOUzMVdtStMLO5FbufAWwM7AS8ksKBeuuRETIkbUThDX4QRaDCRymCFj5P0s3A183smpyHOVOMeCtsB2OciCdn+Tkz+6h6R3KmihPhpLkgRwOpk8R7kNQ6EkO89tJPdH3UbvuubZOxyI7xw++Up9l1asYE9YYIWQZcGITErZKGgZkUQiHmAgpbyKvN7A/xjZDs7l1B2Jwy2sNcaDiO49RNvSuj3wVeC1wraVsKLWKtFBFm9rpUB2Z2G3BbzsNcaEwRdj6orV10zHlylIEqX/jEzLpqqPOqpJoOwopRVfPq1U4mD/XRYHr0RUzE9DsVODVsw30aWFgWvFXSDqMO0ez2nDF0FRqSLs7o5zEzOzTngY7jOFOOHgmNUWL6HZzRPLEuXXRNoa10JUfT+D/A349yX8D/y3mY02ziGfIN32p7c+96QPRd62KnyNEorv92e229o+9xkPKM7psVI5VsaQLs6nJqpjlJmN5pZg91rzY6OULjWDO7brQKkj4z3oE4juNMVnoR66wHnBK22l4L/AC43sxWV+2kq9Aws/N7UcdpEHF20cSXuRfaRaxRxOzyjqjvEt+CUUmNvZtn9Kh9jv8/OmUnyvIOT9g9Yg/1uL7HipqANEBomNkbJK1PYRN5K/AlSQ9QCJAfmNkDOf1kG8IlzQWOBf46tFMxDntpxbE7juNMKRqiabTyEP0gHIRMqW8ATgq5NOZ166PK7qmzgY8Dd9CImI1OFW46J2GjSFFh9t0R9vyCtnbRoVHEpJb2o6A2Sn3DUnGaxrM7q4v/RlbfcbjyDL+L1NjjSLvOJKIZNo21MLP7gK8DX5e0bk6bKrGnHjWzi83sPjP7TesYy0AlfThEY7xT0jmS1pe0iaQrJP0y/N04qn+MpKUhiuPeUfmOku4I904MecQdx3GaQ27cqT5pI62c4COOB4HzJG3drX0VTeNTkr4BXAWsbBWa2YUVB7w58AFgOzP7c8gNvgDYDrjKzI6XdDRwNHCUpO3C/ZdQxEq5UtK2IYf4ycDhwM0UGaf2oZ1D3EmQ3Bk11v5S2kWGRpEiVSelgWRN5BLpTrs9J9l3yi7RkY+1u9+Ke2pPARqyPBU4AXgI+CbFt3UB8FfAvRQ+H3uM1riK0Hg38GJgHdrLUwZUEhrRczeQtAp4BsULOIb2YM+gsPAfRRHy99wQwfE+SUuBeZLuBzY0s5sAJJ0J7IcLDcdxGkZyuXUw7GNmfxtdL5J0s5l9VtInuzWuIjReZmZ/U318nZjZbyV9CXiAIsLi5WZ2uaTnmdnDoc7DUUz4zSk0iRbLQtmqcD6yfC0kHU6hkbDllluO9yVMaWKtoowbz29rMbEXepLxuHJ3zOi710llrit9fCrHWSqtxXB5FVG+M+r6C127mFI0S9MYlvQOilhUAPtH97qOtIpN4+awVDQugq1iPrAVxXLTMyWN5s1Y9u9to5SvXWi2yMzmmtnczTbbrOqQHcdxxows/+gT7wTeBSwHHgnnB0vaAHhft8ZVNI1XAQsl3Udh0xjrltu9gPvM7FEASRcCuwCPSJoVtIxZFC8ICg1ii6j9bIrlrGXhfGS5U4GkfSNao4/r9ISK9ocOv4eqO6wS2kVn7okK+ycibUE542rWDNMZFA3aPWVmvwbekrh9fbf2VYTGPhXqjsYDwE6SnkGxPLUnsAT4I7AQOD78vSjUvxj4pqQTKDSTbYBbzWwo7ALYCbgFOAT49x6N0XEcp3c0aPIQfDPeD8whkgFmtm9O+5yAhbeb2Q6jba9t1cl5oJndIukC4HZgNfBjYBHwLIrkIYdRCJYDQv27wg6ru0P9I8POKSjSHJ4ObEBhAJ+SRvCd/y723m6fxr4ZKTpya0STodg2Uen5MRmTq9QE7Oazq2k3SftJh9Yx/hzZsY2iwzdjbN05U4SmOPcFvgucAnyPMfjcZQUslPSzUe4L2KjKQ83sUxQhfWNWUmgdZfWPA44rKV8CbF/l2Y7jOH3FGrd76i9mduJYG+cIjRdn1BnqXsWpCyXe/WSGvgRVtIvxUMfybo5W5TgDo1maxtckfQq4nE6fu97k0xir17fjOI4TaJbQ+BuKHVOvpdPnrmf5NJymU9WnIUFKM7nxvAqz+ETcpaq7lXY6+IQ15zf/90fyn+84DaRXNg1JpwJvBpab2fah7NPAP9DOCf5JM7tklG7eCmxtZk+PZQxV/DQcx3GcwXI65TtZv2JmLw/HaAID4KfAc8Y6ANc0pjipGEz9shEks9xFuNbhTHh6l+51saQ54+zmecDPJf2ITptGb7bcOo7jOOOg2u6pmZKWRNeLzGxRRrv3STqEwufto2b2+Ch1R+5crYQLjUlAvOup6o6pXnDTN6tpJbHmEJOT5c61DmdCkq9prDCzuRV7Pxn4XHjK54AvA+8ZWUmSrCCZvrtVZ7SHuU3DcRynRkS9safM7BEzGzKzYeC/gFT2vWskvV9SR9RWSetKeq2kMyiicYyKaxpOKbGXd1VNohuxhvC3h0RaR+KfJqWBuNbhTBhq3HLbitkXLt8K3Jmoug+FBnJOCCXyB2B9YDqFz8ZXzOwn3Z7nQmMSE2+VHcSyVWVynP7iIIFR/b99V1uA3HKWCxCnQfQwgq2kcyjyDs2UtIzCPrGHpJcXT+J+4B9Lh1HkB2+ldl0HmAn82cz+UGUMLjQcx3HqpkdhRMzsoJLiU8bQzyrg4a4VS3ChMdlIzGiytI6UM14P6FiGqkjHttyOG5TemHdo+1m3nu5ahzN4GhawcFy40HAcx6kbFxrORGCXd4w9ZHkvmLcw0i4ynqmEvSJue8sZkRE9smPE9WPt4pXvbtf50WmudTgDwJhUQsO33DqO49RMk9K9SnqbpF9KekLSkyGZ3ZO57V3TmAQkNYqYcaRArbK1tUO7SDy/gwztImVfSe2SirWLuJ9XvifSOk51rcPpI83SNL4AvMXM7hlLYxcajuM4NdOwJEyPjFVgwICEhqTnAN+gyLpnFA4n9wLnUeStvR94Ryt+iqRjgMMokj19wMwuC+U70k73egnwwW4u8FOJVFKlZHrWBN12UlW1XdROxhhc63D6RvNsGksknUeR9jUOWHhhTuNB2TS+BvzAzF4MvAy4BzgauMrMtgGuCtdI2g5YALyEwqPx65Kmh35OBg4HtglHWchgx3GcgaEKR5/YEPgT8HrgLeF4c27jvmsakjYEdgMOBQiJQJ6WNJ/C0xHgDOBa4ChgPnCuma0E7pO0FJgn6X5gQzO7KfR7JrAfcGmfXkojibWLqhpFDqX+Fhnf9mQwwpz6Ff+bYs0h1ihStGwgsSHS/TucntIgTcPM3j2e9oPQNLamyDB1mqQfS/qGpGcCz2vFTwl/nxvqbw48GLVfFso2D+cjy9dC0uGSlkha8uijj5ZVcRzHqY2G7Z7aVtJVku4M1y+V9M+57Qdh05gB7AC838xukfQ1wlJUgrJ5po1SvnZhEY9+EcDcuXMbJPN7RB16bWQaUvQAaxd2fX6vvcoB5h5WzbM8awihklZXHo7j5NGsX53/Aj4O/CeAmf1M0jeBf81pPAhNYxmwzMxuCdcXUAiRRyTNgiJqI7A8qr9F1H428FAon11S7jiO0xxCEqaco088w8xuHVGWPWXqu9Aws98BD0p6USjaE7gbuJh2LPeFwEXh/GJggaT1QjjfbYBbwxLWU5J2kiTgkKiNMwaSarJZ++iCqX3kULV+VbL6D7tbbHr7cJyeYplHf1gh6QWtJ0ranwrBCwflp/F+4GxJ6wK/Bt5NIcDOl3QY8ABwAICZ3SXpfArBsho40syGQj9H0N5yeylT3AjuOE4zaVjAwiMplutfLOm3wH3AwbmNByI0QqKPspSGeybqHwccV1K+hMLXwykjEQV2PLS+/OPSDOK2ffxn6oiWW/Le3OqxqZy6aJDQMLNfA3uFDUjTzOypKu3dI9xxHKdmepiE6VQKn4rlZrb9iHsfA74IbGZmK0bp4yMjrgGeAG7zzH1Of4j+IbL8LmLiGFOJ8l6R+sd1j3CnVoyeJWGiWI4/CTgzLpS0BfA6iqX9bswNx/fC9ZuAHwHvlfQtM/vCaI09yq3jOE6NiN75aZjZYuCxkltfAT5B3kLYpsAOZvZRM/sohQDZjMjpejRc03C6E0W8LdMYOvJgdNyo+pzy4hz7Scc/XDwenxY5TaBGm4akfYHfmtlPpax/ui2Bp6PrVcBfm9mfJa1MtFmDCw3HcZyaUX4c1ZmSlkTXi4Jzcnm/0jOAYyniSOXyTeBmSS0XhbcA5wTD+N3dGrvQmCqkvrOpmFAd2sPa3uE53Y1H66i8Oyt62LSh9vlwE6LuOlObaj4YK8ysbGdpihcAWwEtLWM2cLukecEnbu3hmH1O0iXAqyj+M98bdqICvLPbA11oOI7j1ExdfhpmdgftOH2EQK5zu+ye+hJwmpl9bSzPdKExCchat8+ZcSe+2J3aQ1GpQxOZ1pvpfK98P4YzvtW+Y8rpJ70KESLpHIpo4DMlLQM+ZWanVOzm58AiSTOA04BzzOyJ3MYuNBzHceqmR5qGmR3U5f6cjD6+AXwjhHJ6N/AzSTcA/2Vm13Rr70JjqhP7WAwl6sQawLS1y5L/DzXbE6qq/K5dOAOhj2HPcwmJ7F4cjhXAT4GPSPpHM1swWlsXGo7jOHXTIKEh6QSKHVNXA/8WRbz9vKR7u7V3oTFVSNkrorXWTl+HWAVZW2VIRYIdj12iw06SihNVEdcunEHTcu5rEHcC/2xmfyq5N69bY3d9chzHqRkNW9bRr+HEAkPSdEmfAsgxiLumMRlIzO5vOuejpeU7HxTlDk9oFyk/jbGqEkmv8Y5K5fVjbSgnhpXjNIr+5srIYU9JbwcOowgpchpwXW5jFxqO4zg108esfF0xs7+TdCBwB/An4CAzuyG3vQuNKci0RGLHVPymjkl8K59GfD9ul/rnSOUUT7iQxxrF9KfbN4bWieJgJRZX3Y7hNI4GaRqStgE+CHwb+D/AuyT9OGHjWAsXGlOdeEkq9YNva19oqPzXPraZxz/8sQNgTqiRWBCtXi8V6yQxXsdpGA0zhH8PeJ+ZXRlSZX+EIjT6S3IaD8wQHowvP5b0/XC9iaQrJP0y/N04qnuMpKWS7pW0d1S+o6Q7wr0TlRni0XEcp28YxeQs5+gP88zsSgAr+DKwX27jQWoaHwTuATYM10cDV5nZ8ZKODtdHSdoOWEAhBZ8PXClp25An/GTgcOBm4BJgHzxPeFemrU58OVOG8FgDCBEA1WGpjpaMEltxbzmr+5LRvENPKH1mBz4tcCYgTbJpABtKOoMiYOEwcD3F73EWA9E0JM2myBb1jah4PnBGOD+DtuSbD5xrZivN7D5gKTBP0ixgQzO7ycyMIpPVfjiO4zSIXiZh6hGnARcDs4DNKZarTsttPChN46sUWaaeHZU9z8weBjCzhyW1IjduTqFJtFgWylaF85HlayHpcAqNhC233LIHw5/YxJqGpYzScXHp9tvuU/6bzy7f8pvi1tPb2sgr393WOnK0Czd+O42lv0tPOWxmZrGQOF3Sh3Mb913TkNRKin5bbpOSMhulfO1Cs0VmNtfM5m622WaZj3Ucx+kNDdM0Vkg6ONiVp0s6mCL+VBaD0DR2BfaV9EZgfYr1tf8GHpE0K2gZs4Dlof4yYIuo/WzgoVA+u6R86lH1yxZ5nioOa57jw9dl1n/TN6tpFyl+dJprDs4kolGKBu8BTqLIK27AjaEsi75rGmZ2jJnNDiF8FwBXm9nBFGtsC0O1hUArFeHFwAJJ60naCtgGuDUsZT0laaewa+qQqI3jOE5jaJKmYWYPmNm+ZraZmT3XzPYzs9/ktm+Sn8bxwPmSDgMeAA4AMLO7JJ1Pkbt2NXBk2DkFcARwOrABxa4p3zmVQ4czXvk3NdZASrWOaMfUjef1RrtwnEmJAUPNUjXGw0CFhpldC1wbzn8P7JmodxxwXEn5EmD7+kboOI4zfnqlRUg6FWjZhbcPZZ+j2GU6TLGsf6iZ1bZU3yRNw6mRPd7whfZFyiXbfSAcpx56t3vqdAp7xJlR2RfN7P8CSPoA8C/Ae3v1wJF4aHTHcZya6ZVNw8wWA4+NKHsyunwmYzC7S9oht64LDac7ig+BhIk1h+M4o2AVDpgpaUl0HJ7zCEnHSXoQeCeFplGVI3Ir+vKU4zhOjQhQviF8hZnNrfoMMzsWOFbSMcD7gE9VbP8PuXVd05gEaLh99ArXJBynd8gs6+gB3wTePupYpF0lPTOcHyzpBEl/nfsAFxqO4zh1Um15qjIhP0aLfYGfd2lyMvAnSS+jCOf0GzoN66Piy1OTjF0ObKdyjf0nUhrDdZd8ou4hOc4Up3expySdA+xBYftYRrEM9UZJL6LYcvsbuu+cWm1mJmk+8DUzO0XSwi5t1uBCw3Ecp2Z65adhZgeVFJ9SsZungu3jYGA3SdOBdXIb+/KU4zhO3TQrCdOBwErgMDP7HUV08C/mNnZNw3Ecp06s0u6pfvAUxbLUkKRtgRcD5+Q2dqExCYhtFzsf9OVRaha4HcNx+kyjZAaLgVeHlNpXAUsotI935jT25SnHcZya6eOW28zh2J+AtwH/bmZvpUinnYULDcdxnLpplk1Dknam0Cz+J5RNz23sy1OO4zh1YhSbYZvDh4BjgO+E1BNbA9fkNnahMQmIfTOSDtzu2e04A0H0dempK2Z2HXBddP1r4AO57V1oOI7j1M3w4FUNSV81sw9J+h4lpnkz2zenHxcakwzPouc4DaM5y1Nnhb9fGk8nfTeES9pC0jWS7pF0l6QPhvJNJF0h6Zfh78ZRm2MkLZV0r6S9o/IdJd0R7p0YcoU7juM0iibsnjKz28Lf6yjSZ99tZte1jtx+BqFprAY+ama3S3o2cJukK4BDgavM7HhJRwNHA0dJ2g5YQLEl7PnAlZK2DXnCTwYOB24GLgH2YYrkCd/lHZE/hotKx2k2DbBphEn1pyhCpwuYJmk1xbbbz+b203dNw8weNrPbw/lTwD0UbuzzgTNCtTOA/cL5fOBcM1tpZvcBS4F5kmYBG5rZTWZmFFEa98NxHKdRZG63rV+wfAjYFXilmW1qZhsDfwvsKunDuZ0M1KYhaQ7wCuAW4Hlm9jAUgkXSc0O1zSk0iRbLQtmqcD6yvOw5h1NoJGy55ZY9fAX95VVvay9FKhL3N1zwsdL6HXnBo/q77dsOM7P44o/3bHyO45RgQDPCiBwCvM7MVrQKzOzXkg4GLge+ktPJwJz7JD0L+DbwoRE5bteqWlJmo5SvXWi2yMzmmtnczTbbrPpgHcdxxkETbBrAOrHAaGFmj1Ihyu1ANA1J61AIjLPN7MJQ/IikWUHLmAUsD+XLgC2i5rOBh0L57JLyScWr92trBbGd//qUdrHP59sXGfsCdpsf9R/v8Bhuf4HjYGutOnHZ1Vcd3fU5jjOlaYBNA3h6jPc6GMTuKVHEf7/HzE6Ibl0MtBKBLAQuisoXSFpP0lbANsCtYSnrKUk7hT4Pido4juM0A6OYhOUc9fIySU+WHE8Bf5PbySA0jV2BdwF3SPpJKPskcDxwvqTDgAeAAwCCm/v5FFvEVgNHhp1TAEcApwMbUOyamhQ7p2LtIuaH38nQLiKuvbS+aLZ77f5vtfXtOJOLnmbuOxV4M7DczLYPZV8E3kKhLfwKeLeZ/WGtUZhlx5cadQzWDLWpb8ydO9eWLFky6GEAnUtDKRZfVG6o3v2NbSN3nBWsTkHhOFMNSbeZ2dzx9LHR+n9lu2xxSFbdHyz94qjPk7Qb8L/AmZHQeD1wtZmtlvR5ADM7ajxjHg2Pcus4jlMnBgwN5x3dujJbDDw2ouxyM1sdLm+m09bbczyMSJ/I0SpapLSLGFunbeS+LqO+4ziDwsCy44jMlBQvhSwys0UVHvYe4LwK9SvjQsNxHKdu8s0AK8a6HCbpWAq779ljaZ+LC40aiZ3oWqSc6VLG7xQ52ojjOA2gtXuqRiQtpDCQ72k1G6pdaDiO49RNjb/jkvYBjgJ2D2lca8WFRo/Z/S2Rs1xUbuEi1ih++N2Pl547jjPJ6N2W23OAPShsH8soAhAeA6wHXBEcgG82s/f25IEluNBwHMepEzMYGupeL6srO6ik+JSedJ6JC40esPub4sCA0a6m77n24DgOTQkj0hNcaDiO49SNCw0nJg70d+3/uHbhOE5MX+JK9Q0XGo7jOHViYPnOfY3HhUYPUDMSrDiO01QyQoRMFFxoOI7j1IkZDLvQcCKuuby2gJKO40wG3BDuOI7j5GKuaTiO4zh59C4JUxNwoeE4jlMnfQhY2E8mfBImSftIulfSUklHD3o8juM4MQbY0FDWMRGY0EJD0nTg/wFvALYDDpK03WBH5TiOE2EhCVPOMQGY0EIDmAcsNbNfm9nTwLnA/AGPyXEcpwMbtqxjIjDRhcbmwIPR9bJQ1oGkwyUtkbTk0Ucf7dvgHMdxgEmlaUx0Q7hKytYS1yHH7iIASY9K+g0wE1hR7/DGjI9tbDR5bNDs8fnYyvnr8XbwFI9fdqVdMDOzelM/gzVMdKGxDNgiup4NPDRaAzPbDEDSkrHm4q0bH9vYaPLYoNnj87HVh5ntM+gx9JKJvjz1I2AbSVtJWhdYAFw84DE5juNMWia0pmFmqyW9D7gMmA6camZ3DXhYjuM4k5YJLTQAzOwS4JIxNF3U67H0EB/b2Gjy2KDZ4/OxOVnIJpF7u+M4jlMvE92m4TiO4/QRFxqO4zhONlNOaAwiVpWkLSRdI+keSXdJ+mAo30TSFZJ+Gf5uHLU5JozxXkl7R+U7Sroj3DtRUpmvyljGOF3SjyV9v0ljk/QcSRdI+nl4/3Zu0Ng+HD7POyWdI2n9QY5N0qmSlku6Myrr2XgkrSfpvFB+i6Q54xzbF8Pn+jNJ35H0nEGMzamImU2Zg2KH1a+ArYF1gZ8C2/XhubOAHcL5s4FfUMTK+gJwdCg/Gvh8ON8ujG09YKsw5unh3q3AzhSOjZcCb+jRGD8CfBP4frhuxNiAM4C/D+frAs9pwtgoIg/cB2wQrs8HDh3k2IDdgB2AO6Oyno0H+CfgP8L5AuC8cY7t9cCMcP75QY3Nj4rfs0EPoK8vtviyXRZdHwMcM4BxXAS8DrgXmBXKZgH3lo2LYkvxzqHOz6Pyg4D/7MF4ZgNXAa+lLTQGPjZgQ4ofZo0ob8LYWiFsNqHYhfj98CM40LEBc0b8MPdsPK064XwGhfeyxjq2EffeCpw9qLH5kX9MteWprFhVdRLU5lcAtwDPM7OHAcLf54ZqqXFuHs5Hlo+XrwKfAOLgN00Y29bAo8BpYensG5Ke2YSxmdlvgS8BDwAPA0+Y2eVNGNsIejmeNW3MbDXwBLBpj8b5HgrNoYljcyKmmtDIilVV28OlZwHfBj5kZk+OVrWkzEYpH8+Y3gwsN7PbcpskxlDHezuDYknjZDN7BfBHiiWWgY8t2AbmUyyfPB94pqSDmzC2TMYynlrGKulYYDVwdtPG5qzNVBMalWNV9QpJ61AIjLPN7MJQ/IikWeH+LGB5l3EuC+cjy8fDrsC+ku6nCC3/Wkn/3ZCxLQOWmdkt4foCCiHShLHtBdxnZo+a2SrgQmCXhowtppfjWdNG0gxgI+Cx8QxO0kLgzcA7LawtNWVsTjlTTWgMJFZV2OFxCnCPmZ0Q3boYWBjOF1LYOlrlC8KOkK2AbYBbw/LCU5J2Cn0eErUZE2Z2jJnNNrM5FO/H1WZ2cEPG9jvgQUkvCkV7Anc3YWwUy1I7SXpG6HNP4J6GjC2ml+OJ+9qf4rsy5tm8pH2Ao4B9zexPI8Y80LE5ozBoo0q/D+CNFLuXfgUc26dnvopCVf4Z8JNwvJFizfUq4Jfh7yZRm2PDGO8l2k0DzAXuDPdOoofGPmAP2obwRowNeDmwJLx33wU2btDYPgP8PPR7FsVun4GNDTiHwr6yimLmfVgvxwOsD3wLWEqxi2nrcY5tKYUdovU/8R+DGJsf1Q4PI+I4juNkM9WWpxzHcZxx4ELDcRzHycaFhuM4jpONCw3HcRwnGxcajuM4TjYuNBzHcZxsXGg4UwJJcyT9WdJPKrQ5MITa/n6NQ3OcCYULDWcq8Ssze3luZTM7D/j7+objOBMPFxrOlCNoHT8PUXPvlHS2pL0k3RCSFc0b9Bgdp6m40HCmKi8Evga8FHgx8HcU4V4+BnxygONynEbjQsOZqtxnZneY2TBwF3CVFTF17qBIFuQ4TgkuNJypysrofDi6HqbI4+E4TgkuNBzHcZxsXGg4juM42bga7kw5zOx+YPvo+tDUPcdxOnFNw5kqDAEbVXXuA74OPF7XoBxnouFJmBzHcZxsXNNwHMdxsnGh4TiO42TjQsNxHMfJxoWG4ziOk83/D9CypREVJknKAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "do.isel(z=-1).plot();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [],
   "source": [
    "dst = mikeio.read(\"../tests/testdata/dissolved_oxygen.dfs3\", layers=\"top\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<mikeio.Dataset>\n",
       "dims: (time:1, y:112, x:91)\n",
       "time: 2001-12-28 00:00:00 (time-invariant)\n",
       "geometry: Grid2D (ny=112, nx=91)\n",
       "items:\n",
       "  0:  Diss. oxygen (mg/l) <Concentration 3> (mg per liter)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dst"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAY0AAAEKCAYAAADuEgmxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAA9g0lEQVR4nO29ebhdRZX3//kmYXIAgaAdCXRAQV+kHSCmGRRQUHAiqCChRYLSTUvjPAHSbzs1/cMJleaV7rTMjQwiCtogM0RmAw5MolEQIkiIINAOIbl3/f7YdXLq3Oy6p/a9Z5+z773r8zz7uXvXrqpdZ7hn1apVay2ZGY7jOI6Tw7RBD8BxHMeZOLjQcBzHcbJxoeE4juNk40LDcRzHycaFhuM4jpONCw3HcRwnGxcajuM4EwBJW0i6RtI9ku6S9MFQvomkKyT9MvzduNZxuJ+G4zhO85E0C5hlZrdLejZwG7AfcCjwmJkdL+loYGMzO6qucbim4TiOMwEws4fN7PZw/hRwD7A5MB84I1Q7g0KQ1MaU0zRmzpxpc+bMqa3/X9z269r6dpzJxLY7bj3oIXTltttuW2Fmm42nj71f80z7/WNDec/72cq7gL9ERYvMbNHIepLmAIuB7YEHzOw50b3Hzay2JaoZdXXcVObMmcOSJUtq6/910w6orW/HmUxcseRbgx5CVyT9Zrx9rHhsiFsum51Vd51Zv/qLmc3tMqZnAd8GPmRmT0oa7xArMeWERt3MeNEL2xdDw6NXjj/sWOOb3l41tGnlK4gajvqOnxOXT5/e7mf9dUYfy8gxpGiNOaduXCfxViiuE7+O6D1IPit+/+K2cf1U/xadrxPem+HEINctf+8sfn5qoTeukziP+1Fq7KvbY9PqaNYajzn6rtiM6aOPKyb+jqW+k9G5orHw9Kr2+cqny8c11B6vxfWnDMaQdfktyETSOhQC42wzuzAUPyJplpk9HOwey3vysARu03Acx6kRA4axrGM0VKgUpwD3mNkJ0a2LgYXhfCFwUR2vo4VrGj1m+BnrrjnX0OhfAstRK2OxHs9Oo747ZqeJ2bKt09Y6SD23in0ro27H659WPms1orFMT8zKUxpLPEGOL1Ia3rRYMykbcPkzbXri/Uq9jxnlqc/eEtqI1onaxp9lgjX9JN7T5HcvMY2MP0tbN/ruxWPZIPHdjzSjWEva5+X/0u5zRvvBly35dPkgJjDDKVW7GrsC7wLukPSTUPZJ4HjgfEmHAQ8Ata6Ru9BwHMepEcNY1YPlKTO7HkjNNPcc9wMycaHRa+J15Xh2UcVYFc+s4/7imaeimV/H2nP0nPh7mjPLTJSrTKvIeD2m2J6Q0IxSGsu0hE0jngnH2kA8nBnt2W/Hs1ZHdYbydrOsRVWjY877nkGqrTreg7W1tpRNLMsGEz+/Q7lJaISRlhv/RnZoI3GdqDwe++t3/lxicAWX3/R/R73fNAwY6rL0NJGozaYh6VRJyyXdOaL8/ZLuDR6NX4jKj5G0NNzbOyrfUdId4d6JYV0PSetJOi+U3xK2oDmO4zSOXtg0mkKdmsbpwEnAma0CSa+hcER5qZmtlPTcUL4dsAB4CfB84EpJ25rZEHAycDhwM3AJsA9wKXAY8LiZvVDSAuDzwIE1vp4sLvvRpwY9hEq8ft5ny2/EE/3hsc2QO1pFk80O5SLqO94Rllzbz9FSIjq0sxldbAEdu6ESO4pS9RPl49EucrCUzaI1/ooaRUffGZ+BzYjtRFGdabEmHD+r/RkMxxpIbG/q+IzXHtfrdvnXNedX3PjPidE3BwOGJpE/XG2ahpktBh4bUXwEcLyZrQx1WlvD5gPnmtlKM7sPWArMC9vHNjSzm6xYgzmTtrdj7AV5AbBnSwtxHMdpEsOZx0Sg3zaNbYFXSzqOwuvxY2b2IwpX+JujestC2apwPrKc8PdBADNbLekJYFNgxciHSjqcQlthyy237OXrmfBcfuu/lJa/fqdIAxnj1KJDQ0lpC7FtJn5Q4pkdu60U2S6U8NOIZ9+RJqNY9SnbbZXjG5JTHpN6H6v+WiR21HVqSms3u/zm8s86Jp7Fp0jaVzK0yQ4NdkZi7PHH12UH4uteddya8yuuP3bUuoPCMLdpjIMZwMbATsDHKbaJifIdATZKOV3udRaaLTKzuWY2d7PNxhURwHEcpxJmsCrzmAj0W9NYBlwYlppuVTE9nBnKt4jqzQYeCuWzS8qJ2iyTNAPYiLWXw5wu7HLgl0vLb8yYlZaRo6HEGkjK7tFBcnZfrqWU7vaCETvR1tZGOrzsldi9VdUWUPe0LGE/ydEqyohtBK/btbvWEWOpXVUdmhGldWKzhzr6qTSEhiKGJscLAfqvaXwXeC2ApG2BdSmWky4GFoQdUVsB2wC3mtnDwFOSdgoaySG0vR1jL8j9gattqkVfdByn8RgwbHnHRKA2TUPSOcAewExJy4BPAacCp4ZtuE8DC8MP/V2SzgfupthNf2TYOQWF8fx0YAOKXVOXhvJTgLMkLaXQMBbU9VomG7u8o1y7iCdDOx+UqNOqmliHT2koKQ0ky+4RkdzRExkGklOHxE6tlrd60qaS8nVI7ZLKmYrVYPVM2afGyhU3lO9MSmogOZPpintVrOS97DBfTZBASJNJ06hNaJjZQYlbByfqHwccV1K+hCL878jyv1Czu7zjOM54KZz7XGg4k5yxzuZSWkyOjSTWRqziP1mWTSO2gUTag1r+BYn4WJ19ZMRsSsXNqoGB+AXlaAsV7UAdTcvsG7FvyATRLloYsGqiDXoUXGg4juPUiCGGJlFAcRcaU5Abz//omvPU7qmYUvtF98l8B7seED0nMfu+oaI20kGOB3vKq7llYKjq+Z7SLhLP7JXWcdltn+lJP2Ml9oeI/SQ6J9O9WY5paR2a4Ms7wzaxxx/jQsNxHKdG3KbhTCpuPK+a1tHSMDo0ijirWyLgbkoDieu86u1f6vr46yvuzsrapaQSm0YO44laU3GMl/14sNpFispL9WN8z1LfpSuv++SY+usvYshtGo7jOE4OReY+FxrOJCRH61gzy4u1i+GS+2SuaqcqJcpT2khKA4lJRvQtm/1WtT/kZEPMsW9U1ZIGzJWLq8V72mv3f+teqcvnMdHMA2biaeuebTEHSacCbwaWm9n2oexlwH8AzwLuB95pZk/25IElTB7x5ziO01CGUdaRwekU6SFivgEcbWZ/A3yHIq5fbbimMcWpYseI6dAoEvaNHFK2jlRcohwN5Ppvf6y0TjeP6b3nfnrU+yPHlbU+X1W7iGiqHWM85Ngg9nzN/7dW2VXXTgTbRTmFIbw383MzW1yScO5FwOJwfgVwGVBbekPXNBzHcWqlMITnHBRhl5ZEx+EZD7gT2DecH0Bn8Nee45rGFCQn9lQvSGoROVjCONKRna68aY7WUcZlSz5dWr73jt29rnNyeFf2pJ6iXHXNMYMeQk+paAhfYWZzKz7iPcCJkv6FIpDr02WVJJ2Y0deTZjZqOkQXGo7jODUzVKP13sx+Drwe1kQPf1Oi6nyg246RowEXGk4nVT3Cy0jZGZLaRUUP8tS+fJsW79rK0DreFmkdF+ZrHTGxB/ber2yf5+T/tlTk3pStwzMWTzoMscrq+6mV9FwzWy5pGsUP/n8kqn7FzM5I3Gv1tXG357nQmOJUde4rDetQktCoqFwuEXKWrZJ1OgIpJgRI4nf31W9tC5AffmeMAiQRIDC5nTcjvElHqS9PTTp6aQhPpJx4lqQjQ5ULgdNKx2H2VUnTgQ+Y2VdSdbqNwYWG4zhOjRjq2fLUKCknvpbZfkjSfKBUaOTgQsNZQ6x1pCgzoneGsi6nQ3MYR4oyTYs1mfg00joyLPq7zf/imvPFF41/W3tqO29W+ttovBpql8dLYQMJge70jIZ5hN8g6STgPOCPrUIzuz2ncW2vRNKpkpaHLH0j731MkkmaGZUdI2mppHsl7R2V7yjpjnDvxJD2lZAa9rxQfkvJ3mXHcZyBY0aVLbf9YBfgJcBngS+Ho3vgt0CdmsbpwEnAmXGhpC2A1wEPRGXbUaRrfQnwfOBKSduGlK8nA4cDNwOXUHhDXgocBjxuZi+UtAD4PHBgja9nQpMTDDBFhybRxVCbo3X0jMSuXEtZ3RPstm+hdSy+uPeOtJenAizu/Lk1550G9bbRRkPt1xE7Hqa2BveC5LgS257jMOlOOYUhvDdhRHqBmb1mPO1rE21mtpgid/dIvgJ8gs79NPOBc81spZndBywF5kmaBWxoZjeFXOJnAvtFbVo7AS4A9lS3XzTHcZwBMMS0rKMfSHqepFMkXRqut5N0WG77vto0JO0L/NbMfjri931zCk2ixbJQtiqcjyxvtXkQwMxWS3oC2BRYUc/oJzjJcB19HcXa5CgFHRPxyHaRsG/kvKay3Vm7v+kL5ZUTY7zukk90f1CCy29qR3noSGSk8oiFWl2udXTQEXql5E1I/SZVDIdi09vne+7xb6V1cj6D1Jblq686unvjCYShpiVhOp1ih1VLTfwFhX3jlJzGfVtEk/QMikGW6etl76iNUj5am7JnH95yy3/00Udzhus4jtMzmqRpADPN7HzCrMTMVgNDozdp009N4wXAVkBLy5gN3C5pHoUGEcdLmQ08FMpnl5QTtVkmaQawEeXLYZjZImARwNy5c6fMRvjYsS1FHGYjZfdQ2S6lfq4EJn02En4gvX5O3HdUZ483lGsm115aTQOJ7QJ77XZcdCf6EYnDzw8nYqbH4erLbsd+IvHvU45zoXV8CbrXT5DjEPnaPY9vX3T5Hb36iuZrJQYMNysJ0x8lbUr4JCXtBDyR27hvr8TM7jCz55rZHDObQ/Gjv4OZ/Y4iXsqCsCNqK2Ab4FYzexh4StJOwV5xCHBR6PJiYGE43x+4Otg9HMdxGoQYyjz6xEcofj9fIOkGClvxB3Ib16ZplHkumlnpmpmZ3SXpfOBuYDVwZNg5BXAExRrcBhS7pi4N5acAZ0laSqFhLKjppTSeDg0hI1xHR4rVRJiNrrutmiCfq+ZJqnPI0fuxxz6fL61y7Q+O6t5Pyi6QSs7U7XNIzezjPnKmjj3SLNWRTCmjz9Y4E2OMtZKm2kIMGrV7CrgL2J0ipLqAe6mgQNQmNEbxXGzdnzPi+jjguJJ6S4DtS8r/QhEG2HEcp7GYqWnLUzeZ2Q4UwgMASbcDO+Q0do/wyUC8JJ3hbZ3ymI5jM6U05W5zw86ETF2HMi56rjnUrD29Zu+2BnLNZW2tI54td+4Ii2wUvdIYyuqmdl2lvgPxMzs2sCW2sMV1xqqxpF5ndP6a10fv7+UZWl0f6aPjXhJJf0Wx63QDSa+g/clsCDwjtx8XGo7jODVS5NNoxJbbvYFDKTYUnRCVPwVkp0Z0oTEJSM24U7P+ODpsss/El3zQloxa7BI1ahip8cZaR4drRuQFrpTtIjXejtn4GHe5VUx4NRASWoeaYGcrRY3QNEJY9DMkvd3Mvj3WflxoOI7j1Eix5Xbwmoakg83sv4E5kj4y8r6ZnVDSbC1caEwCpkXewsPRJo1xxYHqlgujYx080UdnQKgx12mqdlF5XLHtKbZXdGgXcXnsEV4ek6qj+zKtIsMHo9OOUe75nfWb1y/fncFP2ivRoNhTzwx/nzWeTlxoOI7j1EwTQqOb2X+Gv5/pVnc0XGhMAq77Xnl01lfv184ZkUyf2oPJYU7K1sRmmsp9DpqscSU0ik6v7nKNYtqqcu0ix6bRqbSt/S53+EhE3uGKXlRn1+UaSBapnVdV/TRaZPzmjnlnVs0UodF75OcinQq8GVhuZtuHspdTpHhdn8LP7Z/M7NaStieOPk7LcvBzoeE4jlMzPbRpnM7aKSe+AHzGzC6V9MZwvUdJ29t6MQAXGpOYH363uwaSIkt76AUds/Ie990jkuNK2WBiG0WHptE+j+1QKdsFCS0li1b9qkvpHTvuehNvahA0yWejiHLbm+UpM1tcknDOKHwtoIjB9xAlhN1TXZH072b2/tR9FxqO4zg1UoQRyRYaMyUtia4XhYCro/Eh4DJJX6JYyNul8iA72XW0my40piCxBpKjdcSs8R2I18GT69fl50ntokH77HO0iyx7RYftoly7iH0zUraLXvggdKz5J3NrVOy0X0pHhud7yl4Sax2dDdqnsYd+76mkaawws7kVH3AE8GEz+7akd1DE5durYh/ZDN6k7ziOM8kZRlnHGFkIXBjOvwXM68mgE7imMcVJ2T12mx9pIGWaQVwWtUvt0urwdM7QLqraN2r1ncrRLnI8tpOvNcPDu6NBuS9Ft7p0+F0kYkPNmBbVSfTTo/c6uZOq21R2HFrHIOjl7qkED1FErb0WeC3wy3H2N+pgXWg4juPUTK8M4WUpJ4B/AL4WktH9BTh8lPbTgePNrHy2WPC10cbgQsMZGwmv5KRzeEXbRY6dpNeTt55F6I1ntqnXmmHHiEnOljviQ7V/mFozbeu4n9BQEnU6n1NenDXGuJvU+1FFu6iBVATiXtDLHOGjpJzYMbP9kKQdJSmVtM7MTh+tDxcajuM4NWLA6gYELIz4MXCRpG8Bf2wVmtmF6SZtXGg4tZG0S+TMxBMMIu5bzuy4w8M6adNI9FnVHyN+bqRddNgsymKETU9oFBXtFeOxEQzavjAoGpaEaRPg9xT2jxZG25g+KnWmey1zd/8i8BbgaeBXwLvN7A/h3jHAYcAQ8AEzuyyU70g73eslwAfNzCStR+EVuSPFG3Cgmd1f1+txHMcZE9a75aleYGbvHk/7OsXf6cA+I8quALY3s5cCvwCOAZC0HUWO75eENl8PBhuAkykMO9uEo9XnYcDjZvZC4CtAYjO2MxYWX/TxNQdizWHhiMtSmDKOaVpzoJqPBLKgFcVHrzBbc2h4eM0RlzNM+0gxLTpSj5Lax7RpxTE9OpLvNdHnq65HFtMqHmPtL+O9GDRG7VtuKyFpW0lXSbozXL9U0j/ntq9NaJjZYuCxEWWXm9nqcHkzRQYpgPnAuWa20szuA5YC8yTNAjY0s5uC0eZMYL+oTcst/gJgTyVzYjqO4wyO4aBtdDv6xH9RTNhXAZjZzygm7VkM0qbxHuC8cL45hRBpsSyUrQrnI8tbbR4EMLPVkp4ANgVW1Dhmp/W9jtfE47zkOXkzUtTtEJ6wpbT+V1NDLIsYC3le2h3Z91KZ+DL6jG0mpHwz4pfXqpPIldGM7KNtqmgETfLByMFoRhKmiGeY2a0j5tirU5VHMhChIelYikGe3SoqqWajlI/Wpux5hxP2Lm+55ZaVxuoUS1Ut1oQdiXXU6Ms3nqCD8Q7ApiZeqvycHMGSqp/zg9hlraDqb1XWD3KP1ifG+oNfuV2H92l5lTrDiBhi9XCjDOErJL2A8G5I2h94OLdx34WGpIUUBvI9o33Cy4AtomqzKbwcl9FeworL4zbLglPLRoxYDmsRAn4tApg7d25zAhw5jjMl6Je9IpMjKX4PXyzpt8B9wDtzG/dVaEjaBzgK2N3M/hTduhj4pqQTgOdTGLxvDY4oT0naCbgFOAT496jNQuAmYH/g6pSzitN7Uqlkx6OFK+rJUs6DE+ETTi0JxZPNaAumDUeh0cfzrAlAlpYwsV5Sd6xZy1Nm9mtgL0nPBKaZ2VNV2te55bbM3f0YYD3girCedrOZvdfM7pJ0PnA3xbLVkWY2FLo6gvaW20vDAUUkx7MkLaXQMLINOY7jOP2iaTYNSZtS/B6/CjBJ1wOfNbPf57SvTWgk3N1PGaX+ccBxJeVLgO1Lyv8CHDCeMTrVKQtw+Oq3fmnsHXas4UensbE3MrR35AKaCFpHRIcjXmQVl8UaVqJxHEgwPp8en0f9tE5TgQkTM/7xhGAfV58VflMr//4mNkBc+4P+JWdqktAAzgUWA28P1++k2JSUFU7dPcIdx3FqxBBDzTKEb2Jmn4uu/1XSfrmNXWg44+aH3/nYmvNXva2a1qEO9SKhUVB+kdzdW0P4kq502Cvap8Mz2s+ctjquEjWI26Z+W2ItItIu4v57Yd+oumMqrj+8TvmHMy21mTPxqNJJecZrS07me2V0GwcNM4RfI2kBcH643h/4n9zGLjQcx3FqxBpmCAf+EfgIcFa4ng78UdJHADOzDZMtcaHh9JjrL4y0jrd31zp+GNWP2f0tURKolCNhYv99MhFUTInWkWyXfE7snxJpSZFW0NF0OuXlkY0i6dyXsF0kKakyHi0iRyt4esPppXXW/d+2/Wba0+3zZJKnLmT9/qbq5Lx3NWANEhpm9uzR7kt6iZndlbrfqIU2x3GcyUdeCJEcbUTSqZKWt+JGhbLzJP0kHPdL+sk4B3zWaDdd03Bq4/pvl2sROVz3vfYurd3f9IX2jTq1jpLQImu1y9A6Oogn34lEVB3+Kamo6mOc3uVoF8nUrwlNoLO8fbryOe1OV6/fLl/9h3alDaJNndNWpRJOlQ95ItNDTeN04CSKOHyhbzuwdS7py8AT43zGqIN1oeE4jlMjZjA03BuhYWaLJc0puxcCtr6DzjwZY3rMaDddaDiN57r/+cSa8yxbR4KuWkdiVm7R/1BS60g9M64U/7cpmupH3u+K/DRK/S7GQw3aRTzGlfu0J7h/fPhZa87XW95Wt6avap+v93jHdrJSrrvkE+U3Aru/MdJCM96jWIPtJxV2T82UtCS6XhTCIOXwauARM/tlpcFVxIWG4zhOjRiVlqdWmNncMT7qIOCcMbaNeXq0my40nAlFlo2iYj9V+qscNTZuG19Fto6OXVjTUkaNRHmlXUe90S46Hh8Na6Oz25ty7jjxP9ecb33Fe9acr3pwvTXn6z7Z7nTxxWPTAGJNpMP2FTOtF2raeKg/V0YI2vo2ikym3eruUFL8BPAbM1ttZjuN1t6FhuM4Ts30IZTqXsDPzWxZ15rwdWAH4GcU04Htw/mmkt5rZpeP1tiFhtN4dpsf2TESHuG9mMeld0yNvfcOe0hHebl20VFuiVeY1DpKyvq4qf7G8z8aXbXP739XVByf95jY9tU0erV7qiwQrJmdQhGwNXdp6n7gsJYvRki3/XHgc8CFgAsNx3GcQVHsnuqN9E4EgsXMDq3QzYtj5z0zu1vSK8zs1zkZs11oOI1kTYZA0lpEnVFue7YErYRq1JGtL+F5HdePUsV2pJDtBSk7htMzGpbp515JJ1NEuwU4EPiFpPUIecNHw4WG4zhOzTQpjAhwKPBPwIcopinXAx+jEBiv6dbYhYZTG7se8OX2RTTVuuGCck/xVF6OyjumUrP1nPo10mmhiGb3Fe0eOa9prPm3O0h1MfDdSBMLQ40SGmb2Z+DL4RjJ/3Zr77GnHMdxasYyj34g6c2SfizpMUlPhpTaT+a2rzPd66nAm4HlZrZ9KNuEIkPUHAoL/jvM7PFw7xjgMGAI+ICZXRbKd6Sd7vUS4INmZmH97UyKfcm/Bw40s/vrej1OHh3aRYKO6LeJ8FExHdpFxzp/bzPCJenVzDqVgTBD68iKc+U0EwPrURiRHvFVCp+OO8yqW1vq1DROB/YZUXY0cJWZbQNcFa5bW74WAC8Jbb4uqeX+dDJwOLBNOFp9HgY8bmYvBL4CfL62V+I4jjMOzJR19IkHgTvHIjCg3hzhZYG15lPsMQY4A7gWOCqUn2tmK4H7JC0F5km6H9jQzG4CkHQmsB9waWjz6dDXBcBJkjTWN8LpDb3y2M5iAGvr4/m/ViJWVk7e86TXeK93UqVIvNeDiuU00WjYr9IngEskXQesbBWa2Qk5jfttCH+emT0MYGYPS3puKN8cuDmqtyyUrQrnI8tbbR4Mfa2W9ASwKbBi5EMlHU6hrbDlllv27MU4juN0o2LsqX5wHIXBe31g3aqNm7J7quwdtVHKR2uzdmERJXIRwNy5c5sl86c4nfaKRDTZBDlaTa3/qzl9p9wx4pn7UGzYGceuqljr6LLwHI8lGWPKd0n1BqNpDjCbmNnrx9q437unHpE0CyD8XR7KlwFbRPVmAw+F8tkl5R1tQrCujYDHahu54zjOGDHLO/rElZImjNC4GFgYzhcCF0XlCyStJ2krCoP3rWEp6ylJO4UEI4eMaNPqa3/gardnNAC1D4sOpPbRo/5jTOVHsm3ZUaVu7suo+qzxMK19mLTW0TmG6POIy6dpzZH1njoZCBvOO/rEkcAPJP25aVtu1wqsBRwPnC/pMOAB4AAAM7tL0vnA3cBq4Egza6WjOYL2lttLwwFwCnBWMJo/RrH7ynEcp3k0aDprZs/uXitNnbunSgNrAXsm6h9HYaAZWb6EInTvyPK/EISOM8FIxWCqSq8nZn2cRSfDTXWUV7NvlJKwXcQ5x+PseynGmu/CofDTaICKJumvzOx3463jHuGO4zh10wyX8Et6Uacpu6ecScKN57XzKOxyYNs7PJ5o3fCtduypOJptTN0Ts174kKTGWLnvHP+NHK/xEpJZ+aaXqzqLL3KNoh4Gr2kAL+tiuxDQ1bbhQsOpjViApPjhd9s/Uh0CJMdgPo6lrbIf/HgsHYmfurQbE/FSUVw8DgFS/pyobryu4IKiv/TLCXMUzGx691rd8eUpx3GcOmn5aeQcXZB0qqTlku4cUf5+SfdKuktSIll6b3BNw5m4pLSRMWognZrOmLooHl912SpH6xhKNE68B2vGkHDWa4BddkrRQ2eA04GTKIK1AiDpNRRhlV5qZiujSBu14JqG4zhO3fTIEG5mi1nbifkI4PgQuw8zW75Wwx7imobTU3bdvzyRUkwq1Pn1kU3hVW/r3k+yzx5rIL2iajDHOBjgHm/ovuJQFg4kZfAet5OlU4181W6mpCXR9aIQBmk0tgVeLek44C/Ax8zsR6M1kLQx8Hzgz8D9ZpZtdXGh4TiOUzMVdtStMLO5FbufAWwM7AS8ksKBeuuRETIkbUThDX4QRaDCRymCFj5P0s3A183smpyHOVOMeCtsB2OciCdn+Tkz+6h6R3KmihPhpLkgRwOpk8R7kNQ6EkO89tJPdH3UbvuubZOxyI7xw++Up9l1asYE9YYIWQZcGITErZKGgZkUQiHmAgpbyKvN7A/xjZDs7l1B2Jwy2sNcaDiO49RNvSuj3wVeC1wraVsKLWKtFBFm9rpUB2Z2G3BbzsNcaEwRdj6orV10zHlylIEqX/jEzLpqqPOqpJoOwopRVfPq1U4mD/XRYHr0RUzE9DsVODVsw30aWFgWvFXSDqMO0ez2nDF0FRqSLs7o5zEzOzTngY7jOFOOHgmNUWL6HZzRPLEuXXRNoa10JUfT+D/A349yX8D/y3mY02ziGfIN32p7c+96QPRd62KnyNEorv92e229o+9xkPKM7psVI5VsaQLs6nJqpjlJmN5pZg91rzY6OULjWDO7brQKkj4z3oE4juNMVnoR66wHnBK22l4L/AC43sxWV+2kq9Aws/N7UcdpEHF20cSXuRfaRaxRxOzyjqjvEt+CUUmNvZtn9Kh9jv8/OmUnyvIOT9g9Yg/1uL7HipqANEBomNkbJK1PYRN5K/AlSQ9QCJAfmNkDOf1kG8IlzQWOBf46tFMxDntpxbE7juNMKRqiabTyEP0gHIRMqW8ATgq5NOZ166PK7qmzgY8Dd9CImI1OFW46J2GjSFFh9t0R9vyCtnbRoVHEpJb2o6A2Sn3DUnGaxrM7q4v/RlbfcbjyDL+L1NjjSLvOJKIZNo21MLP7gK8DX5e0bk6bKrGnHjWzi83sPjP7TesYy0AlfThEY7xT0jmS1pe0iaQrJP0y/N04qn+MpKUhiuPeUfmOku4I904MecQdx3GaQ27cqT5pI62c4COOB4HzJG3drX0VTeNTkr4BXAWsbBWa2YUVB7w58AFgOzP7c8gNvgDYDrjKzI6XdDRwNHCUpO3C/ZdQxEq5UtK2IYf4ycDhwM0UGaf2oZ1D3EmQ3Bk11v5S2kWGRpEiVSelgWRN5BLpTrs9J9l3yi7RkY+1u9+Ke2pPARqyPBU4AXgI+CbFt3UB8FfAvRQ+H3uM1riK0Hg38GJgHdrLUwZUEhrRczeQtAp4BsULOIb2YM+gsPAfRRHy99wQwfE+SUuBeZLuBzY0s5sAJJ0J7IcLDcdxGkZyuXUw7GNmfxtdL5J0s5l9VtInuzWuIjReZmZ/U318nZjZbyV9CXiAIsLi5WZ2uaTnmdnDoc7DUUz4zSk0iRbLQtmqcD6yfC0kHU6hkbDllluO9yVMaWKtoowbz29rMbEXepLxuHJ3zOi710llrit9fCrHWSqtxXB5FVG+M+r6C127mFI0S9MYlvQOilhUAPtH97qOtIpN4+awVDQugq1iPrAVxXLTMyWN5s1Y9u9to5SvXWi2yMzmmtnczTbbrOqQHcdxxows/+gT7wTeBSwHHgnnB0vaAHhft8ZVNI1XAQsl3Udh0xjrltu9gPvM7FEASRcCuwCPSJoVtIxZFC8ICg1ii6j9bIrlrGXhfGS5U4GkfSNao4/r9ISK9ocOv4eqO6wS2kVn7okK+ycibUE542rWDNMZFA3aPWVmvwbekrh9fbf2VYTGPhXqjsYDwE6SnkGxPLUnsAT4I7AQOD78vSjUvxj4pqQTKDSTbYBbzWwo7ALYCbgFOAT49x6N0XEcp3c0aPIQfDPeD8whkgFmtm9O+5yAhbeb2Q6jba9t1cl5oJndIukC4HZgNfBjYBHwLIrkIYdRCJYDQv27wg6ru0P9I8POKSjSHJ4ObEBhAJ+SRvCd/y723m6fxr4ZKTpya0STodg2Uen5MRmTq9QE7Oazq2k3SftJh9Yx/hzZsY2iwzdjbN05U4SmOPcFvgucAnyPMfjcZQUslPSzUe4L2KjKQ83sUxQhfWNWUmgdZfWPA44rKV8CbF/l2Y7jOH3FGrd76i9mduJYG+cIjRdn1BnqXsWpCyXe/WSGvgRVtIvxUMfybo5W5TgDo1maxtckfQq4nE6fu97k0xir17fjOI4TaJbQ+BuKHVOvpdPnrmf5NJymU9WnIUFKM7nxvAqz+ETcpaq7lXY6+IQ15zf/90fyn+84DaRXNg1JpwJvBpab2fah7NPAP9DOCf5JM7tklG7eCmxtZk+PZQxV/DQcx3GcwXI65TtZv2JmLw/HaAID4KfAc8Y6ANc0pjipGEz9shEks9xFuNbhTHh6l+51saQ54+zmecDPJf2ITptGb7bcOo7jOOOg2u6pmZKWRNeLzGxRRrv3STqEwufto2b2+Ch1R+5crYQLjUlAvOup6o6pXnDTN6tpJbHmEJOT5c61DmdCkq9prDCzuRV7Pxn4XHjK54AvA+8ZWUmSrCCZvrtVZ7SHuU3DcRynRkS9safM7BEzGzKzYeC/gFT2vWskvV9SR9RWSetKeq2kMyiicYyKaxpOKbGXd1VNohuxhvC3h0RaR+KfJqWBuNbhTBhq3HLbitkXLt8K3Jmoug+FBnJOCCXyB2B9YDqFz8ZXzOwn3Z7nQmMSE2+VHcSyVWVynP7iIIFR/b99V1uA3HKWCxCnQfQwgq2kcyjyDs2UtIzCPrGHpJcXT+J+4B9Lh1HkB2+ldl0HmAn82cz+UGUMLjQcx3HqpkdhRMzsoJLiU8bQzyrg4a4VS3ChMdlIzGiytI6UM14P6FiGqkjHttyOG5TemHdo+1m3nu5ahzN4GhawcFy40HAcx6kbFxrORGCXd4w9ZHkvmLcw0i4ynqmEvSJue8sZkRE9smPE9WPt4pXvbtf50WmudTgDwJhUQsO33DqO49RMk9K9SnqbpF9KekLSkyGZ3ZO57V3TmAQkNYqYcaRArbK1tUO7SDy/gwztImVfSe2SirWLuJ9XvifSOk51rcPpI83SNL4AvMXM7hlLYxcajuM4NdOwJEyPjFVgwICEhqTnAN+gyLpnFA4n9wLnUeStvR94Ryt+iqRjgMMokj19wMwuC+U70k73egnwwW4u8FOJVFKlZHrWBN12UlW1XdROxhhc63D6RvNsGksknUeR9jUOWHhhTuNB2TS+BvzAzF4MvAy4BzgauMrMtgGuCtdI2g5YALyEwqPx65Kmh35OBg4HtglHWchgx3GcgaEKR5/YEPgT8HrgLeF4c27jvmsakjYEdgMOBQiJQJ6WNJ/C0xHgDOBa4ChgPnCuma0E7pO0FJgn6X5gQzO7KfR7JrAfcGmfXkojibWLqhpFDqX+Fhnf9mQwwpz6Ff+bYs0h1ihStGwgsSHS/TucntIgTcPM3j2e9oPQNLamyDB1mqQfS/qGpGcCz2vFTwl/nxvqbw48GLVfFso2D+cjy9dC0uGSlkha8uijj5ZVcRzHqY2G7Z7aVtJVku4M1y+V9M+57Qdh05gB7AC838xukfQ1wlJUgrJ5po1SvnZhEY9+EcDcuXMbJPN7RB16bWQaUvQAaxd2fX6vvcoB5h5WzbM8awihklZXHo7j5NGsX53/Aj4O/CeAmf1M0jeBf81pPAhNYxmwzMxuCdcXUAiRRyTNgiJqI7A8qr9F1H428FAon11S7jiO0xxCEqaco088w8xuHVGWPWXqu9Aws98BD0p6USjaE7gbuJh2LPeFwEXh/GJggaT1QjjfbYBbwxLWU5J2kiTgkKiNMwaSarJZ++iCqX3kULV+VbL6D7tbbHr7cJyeYplHf1gh6QWtJ0ranwrBCwflp/F+4GxJ6wK/Bt5NIcDOl3QY8ABwAICZ3SXpfArBsho40syGQj9H0N5yeylT3AjuOE4zaVjAwiMplutfLOm3wH3AwbmNByI0QqKPspSGeybqHwccV1K+hMLXwykjEQV2PLS+/OPSDOK2ffxn6oiWW/Le3OqxqZy6aJDQMLNfA3uFDUjTzOypKu3dI9xxHKdmepiE6VQKn4rlZrb9iHsfA74IbGZmK0bp4yMjrgGeAG7zzH1Of4j+IbL8LmLiGFOJ8l6R+sd1j3CnVoyeJWGiWI4/CTgzLpS0BfA6iqX9bswNx/fC9ZuAHwHvlfQtM/vCaI09yq3jOE6NiN75aZjZYuCxkltfAT5B3kLYpsAOZvZRM/sohQDZjMjpejRc03C6E0W8LdMYOvJgdNyo+pzy4hz7Scc/XDwenxY5TaBGm4akfYHfmtlPpax/ui2Bp6PrVcBfm9mfJa1MtFmDCw3HcZyaUX4c1ZmSlkTXi4Jzcnm/0jOAYyniSOXyTeBmSS0XhbcA5wTD+N3dGrvQmCqkvrOpmFAd2sPa3uE53Y1H66i8Oyt62LSh9vlwE6LuOlObaj4YK8ysbGdpihcAWwEtLWM2cLukecEnbu3hmH1O0iXAqyj+M98bdqICvLPbA11oOI7j1ExdfhpmdgftOH2EQK5zu+ye+hJwmpl9bSzPdKExCchat8+ZcSe+2J3aQ1GpQxOZ1pvpfK98P4YzvtW+Y8rpJ70KESLpHIpo4DMlLQM+ZWanVOzm58AiSTOA04BzzOyJ3MYuNBzHceqmR5qGmR3U5f6cjD6+AXwjhHJ6N/AzSTcA/2Vm13Rr70JjqhP7WAwl6sQawLS1y5L/DzXbE6qq/K5dOAOhj2HPcwmJ7F4cjhXAT4GPSPpHM1swWlsXGo7jOHXTIKEh6QSKHVNXA/8WRbz9vKR7u7V3oTFVSNkrorXWTl+HWAVZW2VIRYIdj12iw06SihNVEdcunEHTcu5rEHcC/2xmfyq5N69bY3d9chzHqRkNW9bRr+HEAkPSdEmfAsgxiLumMRlIzO5vOuejpeU7HxTlDk9oFyk/jbGqEkmv8Y5K5fVjbSgnhpXjNIr+5srIYU9JbwcOowgpchpwXW5jFxqO4zg108esfF0xs7+TdCBwB/An4CAzuyG3vQuNKci0RGLHVPymjkl8K59GfD9ul/rnSOUUT7iQxxrF9KfbN4bWieJgJRZX3Y7hNI4GaRqStgE+CHwb+D/AuyT9OGHjWAsXGlOdeEkq9YNva19oqPzXPraZxz/8sQNgTqiRWBCtXi8V6yQxXsdpGA0zhH8PeJ+ZXRlSZX+EIjT6S3IaD8wQHowvP5b0/XC9iaQrJP0y/N04qnuMpKWS7pW0d1S+o6Q7wr0TlRni0XEcp28YxeQs5+gP88zsSgAr+DKwX27jQWoaHwTuATYM10cDV5nZ8ZKODtdHSdoOWEAhBZ8PXClp25An/GTgcOBm4BJgHzxPeFemrU58OVOG8FgDCBEA1WGpjpaMEltxbzmr+5LRvENPKH1mBz4tcCYgTbJpABtKOoMiYOEwcD3F73EWA9E0JM2myBb1jah4PnBGOD+DtuSbD5xrZivN7D5gKTBP0ixgQzO7ycyMIpPVfjiO4zSIXiZh6hGnARcDs4DNKZarTsttPChN46sUWaaeHZU9z8weBjCzhyW1IjduTqFJtFgWylaF85HlayHpcAqNhC233LIHw5/YxJqGpYzScXHp9tvuU/6bzy7f8pvi1tPb2sgr393WOnK0Czd+O42lv0tPOWxmZrGQOF3Sh3Mb913TkNRKin5bbpOSMhulfO1Cs0VmNtfM5m622WaZj3Ucx+kNDdM0Vkg6ONiVp0s6mCL+VBaD0DR2BfaV9EZgfYr1tf8GHpE0K2gZs4Dlof4yYIuo/WzgoVA+u6R86lH1yxZ5nioOa57jw9dl1n/TN6tpFyl+dJprDs4kolGKBu8BTqLIK27AjaEsi75rGmZ2jJnNDiF8FwBXm9nBFGtsC0O1hUArFeHFwAJJ60naCtgGuDUsZT0laaewa+qQqI3jOE5jaJKmYWYPmNm+ZraZmT3XzPYzs9/ktm+Sn8bxwPmSDgMeAA4AMLO7JJ1Pkbt2NXBk2DkFcARwOrABxa4p3zmVQ4czXvk3NdZASrWOaMfUjef1RrtwnEmJAUPNUjXGw0CFhpldC1wbzn8P7JmodxxwXEn5EmD7+kboOI4zfnqlRUg6FWjZhbcPZZ+j2GU6TLGsf6iZ1bZU3yRNw6mRPd7whfZFyiXbfSAcpx56t3vqdAp7xJlR2RfN7P8CSPoA8C/Ae3v1wJF4aHTHcZya6ZVNw8wWA4+NKHsyunwmYzC7S9oht64LDac7ig+BhIk1h+M4o2AVDpgpaUl0HJ7zCEnHSXoQeCeFplGVI3Ir+vKU4zhOjQhQviF8hZnNrfoMMzsWOFbSMcD7gE9VbP8PuXVd05gEaLh99ArXJBynd8gs6+gB3wTePupYpF0lPTOcHyzpBEl/nfsAFxqO4zh1Um15qjIhP0aLfYGfd2lyMvAnSS+jCOf0GzoN66Piy1OTjF0ObKdyjf0nUhrDdZd8ou4hOc4Up3expySdA+xBYftYRrEM9UZJL6LYcvsbuu+cWm1mJmk+8DUzO0XSwi5t1uBCw3Ecp2Z65adhZgeVFJ9SsZungu3jYGA3SdOBdXIb+/KU4zhO3TQrCdOBwErgMDP7HUV08C/mNnZNw3Ecp06s0u6pfvAUxbLUkKRtgRcD5+Q2dqExCYhtFzsf9OVRaha4HcNx+kyjZAaLgVeHlNpXAUsotI935jT25SnHcZya6eOW28zh2J+AtwH/bmZvpUinnYULDcdxnLpplk1Dknam0Cz+J5RNz23sy1OO4zh1YhSbYZvDh4BjgO+E1BNbA9fkNnahMQmIfTOSDtzu2e04A0H0dempK2Z2HXBddP1r4AO57V1oOI7j1M3w4FUNSV81sw9J+h4lpnkz2zenHxcakwzPouc4DaM5y1Nnhb9fGk8nfTeES9pC0jWS7pF0l6QPhvJNJF0h6Zfh78ZRm2MkLZV0r6S9o/IdJd0R7p0YcoU7juM0iibsnjKz28Lf6yjSZ99tZte1jtx+BqFprAY+ama3S3o2cJukK4BDgavM7HhJRwNHA0dJ2g5YQLEl7PnAlZK2DXnCTwYOB24GLgH2YYrkCd/lHZE/hotKx2k2DbBphEn1pyhCpwuYJmk1xbbbz+b203dNw8weNrPbw/lTwD0UbuzzgTNCtTOA/cL5fOBcM1tpZvcBS4F5kmYBG5rZTWZmFFEa98NxHKdRZG63rV+wfAjYFXilmW1qZhsDfwvsKunDuZ0M1KYhaQ7wCuAW4Hlm9jAUgkXSc0O1zSk0iRbLQtmqcD6yvOw5h1NoJGy55ZY9fAX95VVvay9FKhL3N1zwsdL6HXnBo/q77dsOM7P44o/3bHyO45RgQDPCiBwCvM7MVrQKzOzXkg4GLge+ktPJwJz7JD0L+DbwoRE5bteqWlJmo5SvXWi2yMzmmtnczTbbrPpgHcdxxkETbBrAOrHAaGFmj1Ihyu1ANA1J61AIjLPN7MJQ/IikWUHLmAUsD+XLgC2i5rOBh0L57JLyScWr92trBbGd//qUdrHP59sXGfsCdpsf9R/v8Bhuf4HjYGutOnHZ1Vcd3fU5jjOlaYBNA3h6jPc6GMTuKVHEf7/HzE6Ibl0MtBKBLAQuisoXSFpP0lbANsCtYSnrKUk7hT4Pido4juM0A6OYhOUc9fIySU+WHE8Bf5PbySA0jV2BdwF3SPpJKPskcDxwvqTDgAeAAwCCm/v5FFvEVgNHhp1TAEcApwMbUOyamhQ7p2LtIuaH38nQLiKuvbS+aLZ77f5vtfXtOJOLnmbuOxV4M7DczLYPZV8E3kKhLfwKeLeZ/WGtUZhlx5cadQzWDLWpb8ydO9eWLFky6GEAnUtDKRZfVG6o3v2NbSN3nBWsTkHhOFMNSbeZ2dzx9LHR+n9lu2xxSFbdHyz94qjPk7Qb8L/AmZHQeD1wtZmtlvR5ADM7ajxjHg2Pcus4jlMnBgwN5x3dujJbDDw2ouxyM1sdLm+m09bbczyMSJ/I0SpapLSLGFunbeS+LqO+4ziDwsCy44jMlBQvhSwys0UVHvYe4LwK9SvjQsNxHKdu8s0AK8a6HCbpWAq779ljaZ+LC40aiZ3oWqSc6VLG7xQ52ojjOA2gtXuqRiQtpDCQ72k1G6pdaDiO49RNjb/jkvYBjgJ2D2lca8WFRo/Z/S2Rs1xUbuEi1ih++N2Pl547jjPJ6N2W23OAPShsH8soAhAeA6wHXBEcgG82s/f25IEluNBwHMepEzMYGupeL6srO6ik+JSedJ6JC40esPub4sCA0a6m77n24DgOTQkj0hNcaDiO49SNCw0nJg70d+3/uHbhOE5MX+JK9Q0XGo7jOHViYPnOfY3HhUYPUDMSrDiO01QyQoRMFFxoOI7j1IkZDLvQcCKuuby2gJKO40wG3BDuOI7j5GKuaTiO4zh59C4JUxNwoeE4jlMnfQhY2E8mfBImSftIulfSUklHD3o8juM4MQbY0FDWMRGY0EJD0nTg/wFvALYDDpK03WBH5TiOE2EhCVPOMQGY0EIDmAcsNbNfm9nTwLnA/AGPyXEcpwMbtqxjIjDRhcbmwIPR9bJQ1oGkwyUtkbTk0Ucf7dvgHMdxgEmlaUx0Q7hKytYS1yHH7iIASY9K+g0wE1hR7/DGjI9tbDR5bNDs8fnYyvnr8XbwFI9fdqVdMDOzelM/gzVMdKGxDNgiup4NPDRaAzPbDEDSkrHm4q0bH9vYaPLYoNnj87HVh5ntM+gx9JKJvjz1I2AbSVtJWhdYAFw84DE5juNMWia0pmFmqyW9D7gMmA6camZ3DXhYjuM4k5YJLTQAzOwS4JIxNF3U67H0EB/b2Gjy2KDZ4/OxOVnIJpF7u+M4jlMvE92m4TiO4/QRFxqO4zhONlNOaAwiVpWkLSRdI+keSXdJ+mAo30TSFZJ+Gf5uHLU5JozxXkl7R+U7Sroj3DtRUpmvyljGOF3SjyV9v0ljk/QcSRdI+nl4/3Zu0Ng+HD7POyWdI2n9QY5N0qmSlku6Myrr2XgkrSfpvFB+i6Q54xzbF8Pn+jNJ35H0nEGMzamImU2Zg2KH1a+ArYF1gZ8C2/XhubOAHcL5s4FfUMTK+gJwdCg/Gvh8ON8ujG09YKsw5unh3q3AzhSOjZcCb+jRGD8CfBP4frhuxNiAM4C/D+frAs9pwtgoIg/cB2wQrs8HDh3k2IDdgB2AO6Oyno0H+CfgP8L5AuC8cY7t9cCMcP75QY3Nj4rfs0EPoK8vtviyXRZdHwMcM4BxXAS8DrgXmBXKZgH3lo2LYkvxzqHOz6Pyg4D/7MF4ZgNXAa+lLTQGPjZgQ4ofZo0ob8LYWiFsNqHYhfj98CM40LEBc0b8MPdsPK064XwGhfeyxjq2EffeCpw9qLH5kX9MteWprFhVdRLU5lcAtwDPM7OHAcLf54ZqqXFuHs5Hlo+XrwKfAOLgN00Y29bAo8BpYensG5Ke2YSxmdlvgS8BDwAPA0+Y2eVNGNsIejmeNW3MbDXwBLBpj8b5HgrNoYljcyKmmtDIilVV28OlZwHfBj5kZk+OVrWkzEYpH8+Y3gwsN7PbcpskxlDHezuDYknjZDN7BfBHiiWWgY8t2AbmUyyfPB94pqSDmzC2TMYynlrGKulYYDVwdtPG5qzNVBMalWNV9QpJ61AIjLPN7MJQ/IikWeH+LGB5l3EuC+cjy8fDrsC+ku6nCC3/Wkn/3ZCxLQOWmdkt4foCCiHShLHtBdxnZo+a2SrgQmCXhowtppfjWdNG0gxgI+Cx8QxO0kLgzcA7LawtNWVsTjlTTWgMJFZV2OFxCnCPmZ0Q3boYWBjOF1LYOlrlC8KOkK2AbYBbw/LCU5J2Cn0eErUZE2Z2jJnNNrM5FO/H1WZ2cEPG9jvgQUkvCkV7Anc3YWwUy1I7SXpG6HNP4J6GjC2ml+OJ+9qf4rsy5tm8pH2Ao4B9zexPI8Y80LE5ozBoo0q/D+CNFLuXfgUc26dnvopCVf4Z8JNwvJFizfUq4Jfh7yZRm2PDGO8l2k0DzAXuDPdOoofGPmAP2obwRowNeDmwJLx33wU2btDYPgP8PPR7FsVun4GNDTiHwr6yimLmfVgvxwOsD3wLWEqxi2nrcY5tKYUdovU/8R+DGJsf1Q4PI+I4juNkM9WWpxzHcZxx4ELDcRzHycaFhuM4jpONCw3HcRwnGxcajuM4TjYuNBzHcZxsXGg4UwJJcyT9WdJPKrQ5MITa/n6NQ3OcCYULDWcq8Ssze3luZTM7D/j7+objOBMPFxrOlCNoHT8PUXPvlHS2pL0k3RCSFc0b9Bgdp6m40HCmKi8Evga8FHgx8HcU4V4+BnxygONynEbjQsOZqtxnZneY2TBwF3CVFTF17qBIFuQ4TgkuNJypysrofDi6HqbI4+E4TgkuNBzHcZxsXGg4juM42bga7kw5zOx+YPvo+tDUPcdxOnFNw5kqDAEbVXXuA74OPF7XoBxnouFJmBzHcZxsXNNwHMdxsnGh4TiO42TjQsNxHMfJxoWG4ziOk83/D9CypREVJknKAAAAAElFTkSuQmCC",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "dst[0].plot();"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<mikeio.Dataset>\n",
       "dims: (time:1, y:112, x:91)\n",
       "time: 2001-12-28 00:00:00 (time-invariant)\n",
       "geometry: Grid2D (ny=112, nx=91)\n",
       "items:\n",
       "  0:  Diss. oxygen (mg/l) <Concentration 3> (mg per liter)"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dsb = mikeio.read(\"../tests/testdata/dissolved_oxygen.dfs3\",layers=\"bottom\")\n",
    "dsb"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmUAAAJcCAYAAAC4+tNCAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAABUGklEQVR4nO3de5hkVXnv8d+vewa8otxU5OJgBBXwBuMEgiciCAyIDiaoQySgIYfoQaNRIzcjeEHBu2jQEEDAKBcRZFQGIXg7coBhQJS7jIAygjCIwKjATHe/54/aHWuqVnev3rWranfX9/M89XTVqr32Xnt3dffqtd/1LkeEAAAA0F9D/W4AAAAA6JQBAADUAp0yAACAGqBTBgAAUAN0ygAAAGqAThkAAEAN0CkDAACoATplwAxk+y7bj9r+g+3f2/6u7S0z677F9k9ays6w/dHutBYAkINOGTBzvTYiniJpM0n3SfpCn9sDAOgAnTJghouIxySdL2m78TLbT7N9lu1Vtn9l+wO2h2y/UNKXJe1SjLI9ZPswSW+W9P6i7NvFPl5o+4fFNjfZfl3T/s+wfbLtpUWdK2w/y/bnipG7W22/bKI22/4r29fYfrj4+ldF+Zts32F7g+L1PrZ/a3tT2/9u+9Mt+/m27XcXz3e0/VPbq21/w/a5zaN/tvezfX1xPv/P9oub3rvL9vts/7xo07m2n9DBtwUApo1OGTDD2X6SpDdJuqqp+AuSnibpuZJeKelgSW+NiFskvU3SlRHxlIh4ekScIulrkj5RlL3W9lxJ35Z0qaRnSHqnpK/Zfn7TMd4o6QOSNpH0uKQrJV1XvD5f0mcmaO9Gkr4r6SRJGxfbfdf2xhFxbrGfk2xvLOk0Sf8YEasknSnpQNtDxX42kbSHpLNtryfpQklnSNpI0tmSXt90zB0lnS7pn4pj/oekJbbXbzmfhZK2lvRiSW+Z6JoDQDfQKQNmrm/ZfkjSI5L2lPRJSbI9rEYn7aiIWB0Rd0n6tKS/n8a+d5b0FEknRMSaiPi+pO9IOrBpmwsj4tpipO5CSY9FxFkRMSrpXEkTjZS9RtLtEfHViBiJiLMl3SrptcX7h0vaXdIPJX07Ir4jSRGxTNLDanTEJGmxpB9GxH1Fe+dIOiki1kbEBZKWNR3zf0v6j4i4OiJGI+JMNTqSOzdtc1JE3BMRD6rRIX1p/uUCgM7RKQNmrv0j4umS1pf0Dkk/sv0sNUaq1pP0q6ZtfyVp82ns+9mS7o6IsUn2cV/T80cTr58yyb5/1VL2P/uOiIckfUPSDmp0JpudKemg4vlBkr7atM/fREQ0bXt30/PnSHpvcevyoaIzu2VRb9xvm57/aZL2A0BX0CkDZrhi5OcCSaOSXiHpAUlr1eiIjNtK0m/Gq6R20/L6Hklbjt8qTOyjE/e0tG2dfdt+qaR/UOMW5Ekt2/2XpEW2XyLphZK+VZTfK2lz227atnk26t2Sji9u144/nlSM0gFALdApA2Y4NyyStKGkW4rbh+dJOt72U20/R9J71OjQSI0RrS2KOCw1lT236fXVkv6oRvD/XNu7qXF78ZwKmnyxpG1t/53tObbfpMYkhe8UwfX/JeloSW9Vo6P1f8YrRsRKSdeoMUL2zYh4tHjrSjU6pe8o9rlI0oKmY/6npLfZ/sviej3Z9mtsP7WC8wGAStApA2aub9v+gxoxZcdLOiQibiree6canao7JP1E0tfVCHSXpO9LuknSb20/UJSdJmm74tbetyJijaTXSdpHjZG3kyUdHBG3dtroiPidpP0kvVfS7yS9X9J+EfGApI9LWhkRX4qIx9W4RflR29s07eJMSS/Sn29dqmjv30g6VNJDRb3vqBE3pohYrkZc2Rcl/V7SChHID6BmvG4IBgDUm+2/VmM0bV5LzFvrdldL+nJEfKVnjQOADjBSBmDGKFJ1vEvSqa0dMtuvLHKlzbF9iBppLS7pRzsBoAw6ZQBmhCLx7UNqrGDwucQmz5f0MzXSZrxX0gERcW+v2gegfmxvafsHtm8pkmC/qyjfyPZltm8vvm44Qf2Ftm+zvcL2kV1vL7cvAQDAbGR7M0mbRcR1xcSeayXtr0ZM6YMRcULR2dowIo5oqTss6Rdq5IEcn2R0YETc3K32MlIGAABmpYi4NyKuK56vlnSLGjkRF6kxaUjF1/0T1RdIWhERdxSTic4p6nXNnG7uvI422WSTmDdvXr+bkXTD73479UYYLAxkY1B56k0m8qKNn1VdO2aBa6+99oGI2LRXx9v7VU+O3z042qvD6dqfP36TpMeaik4plo9bh+15aqw0crWkZ46HN0TEvbafkdj15lo3CfVKSX9ZVbtTBq5TNm/ePC1fvrzfzUiad+aJ/Tkwf/hrK0ZL/mXq4A/aQOGz310dfA49VP6bs/wtR0y90QCx3bqCRlf97sFRLfveVj073vBmtz8WEfMn28b2UyR9U9K7I+KRdfNMT1wtUdbV3xoD1ykDAADdE5LGNGG2mp4rZm1/U9LXitVPJOk+25sVo2SbSbo/UXWl1l0ZZAs1ViTpGmLKAADArFQsvXaaGqudfKbprSWSDimeHyLpokT1ayRtY3vrYgWUxUW9rmGkDAAAVCg0OnFe517bVdLfS7rB9vVF2dGSTpB0nu1DJf1a0hskyfaz1ciDuG9EjNh+h6TvSRqWdHrTqildQadsNkveMk/cDo8uByB1svvcu/dVnkIv4oy6fc3RfWW/hb34fGW2LS+sRpITjc6pmzrXxGc/mZmJnxFUICJ+ook/rXsktr9H0r5Nry9WY73enuD2JQAAQA0wUgYAACrTCPRnanMZjJQBAADUACNlAACgUnVKiTGT0Cmru1R4YiroNhkU275dKiFjepA5I8g2VTF3ckGufsT65kY/555W8jplVs6OxEZXJb4NZb81ye98BxNaku0YTgXmJ34fZAbwp7aLjED8GE3cjElE9XssddDyvzfmffWEtrLhue2dhF8uPrr0MYBuoFMGAAAqEwqNJqfUYirElAEAANQAI2UAAKBSzL4sh5EyAACAGmCkrE4yg/qH5rQHrMZoXtSxUwHASdX9l5OcXJAK7E1WzjxIZnNTwcmtQcydhUKkApZTDcncW+6cg5b9la3XM9mTRPogN5i+ZCC6EwdIXo7UMRM/S6myoeHEzLeqs/xnfHhiLHMlgIq/96mg/uE5o21l23zjo21lrae13nojbdvcvP9xpds2CELSKCNlpTBSBgAAUAOMlAEAgEoRU1YOI2UAAAA1wEgZAACoTEjkKSuJTlmN3HXwEf1uwozwwgs/VLrumsfntpW1BcmnMo6nxpQTgd5ja9s3TE4uyJ38kDmWXTZOOplJvQPZv4crPGxHQfiZ39fsjPapw7Ze4+T+U+1IlCU+N6m25X5ekxLnlfq+tn1ecyevJM5h7vpr28pGE6sBjK4ZzjpGaqLDUPLatW83Ojb1D9123zqurYzgf1Sha7cvbZ9u+37bN7aUv9P2bbZvsv2JpvKjbK8o3tu7qXwn2zcU751kN3592V7f9rlF+dW253XrXAAAQL6xHj5mk27GlJ0haWFzge1XSVok6cURsb2kTxXl20laLGn7os7Jtsf/JfqSpMMkbVM8xvd5qKTfR8TzJH1W0oldPBcAAICu6lqnLCJ+LOnBluK3SzohIh4vtrm/KF8k6ZyIeDwi7pS0QtIC25tJ2iAiroyIkHSWpP2b6pxZPD9f0h7jo2gAAAAzTa9jyraV9L9sHy/pMUnvi4hrJG0u6aqm7VYWZWuL563lKr7eLUkRMWL7YUkbS3qg9aC2D1NjtE1bbbVVleeDPrjl9cdmbZeKPVsvEbtSpbVuj3kZS8TGpP59GGvPbalIjc2nYplaA3qSm6Ri1hKBQNmJfTNjo5LtTWzW7bjg5L+g5Q+aHWdWOuAvFSuWl0x6OBFT9YsD/q1kQ6TnndeeZFUZ4V2ja/NiwEZHMmPFEueavCaJ8099v1KGh9atm4oxm5PY//YXHddWdtOi9rJBEAqSx5bU65QYcyRtKGlnSf8q6bxidGuivzKT/fXJzpMeEadExPyImL/ppptOv9UAAABd1uuRspWSLihuRS5zY+rLJkX5lk3bbSHpnqJ8i0S5muqstD1H0tPUfrsUAAD0UkijDJSV0uuRsm9J2l2SbG8raT01bjcukbS4mFG5tRoB/csi4l5Jq23vXIyoHSzpomJfSyQdUjw/QNL3i84eAADAjNO1kTLbZ0vaTdImtldKOlbS6ZJOL9JkrJF0SNGRusn2eZJuljQi6fCIGI+websaMzmfKGlp8ZCk0yR91fYKNUbIFnfrXAAAQJ7Q7EtV0Std65RFxIETvHXQBNsfL+n4RPlySTskyh+T9IZO2oiZadfL8pLs3vL67mZJSSWQnDs3Ea2fKksYSQQ7jzzeXpYzHJwMak4FpqcSZaaC/3Olqqbakkqom8rF2jp/oZP51cm2dbC/lOT+Wg+cSiacF8Du4fayX77pmMzGlbfijR9oK2sNbF+zpv3PyXDis5/6zKXucaS+1+myxHVq3ywpZ7IGU/rRS2T0BwAAFbJG6c6WwoLkAAAANcBIGQAAqExIGmPaXSmMlAEAANQAI2WotV0uPbKtbCgzVCF3QkBZN++fN5EgNSEglfl9zpxEUHQqRj4VnN8SsJwKYE5mV88Mwk9KbZab0jnzv+juL5yWnYM6b2/JCRZT7z8VwJ8K9O9FUH+unGz123yjfSUAJ7Lh5wb6p+Rm78/N6J+7HSZHTFk5jJQBAADUACNlAACgMiFGyspipAwAAKAGGCkDAACVGsuNTcU66JRhxkn9sA/1ITg3NQkh5eb9Tyh9jBde+KFEafuEgLUtQfypIOk567XXS64YkDl+nsyGPpYZ/d/JBIOcfWXLDAjPbVrGqgHJTP01D+ova731RtrKWj+rDakVHjr4PlQoN/C/H7+DMPtw+xIAAKAGGCkDAACVIdC/PEbKAAAAaoCRMgAAUJmQNcqYTyl0ylBrV+7VHiTf7Uz9Vc8aSk0ISGX0T7nl9R8vdczUBIGRkfYA6+FE8H+u1AoBqcz0MdpJlHxJuYsSJALsU9+b7CzvqTkNrastzNKg/pRU1v/tL2ovW5v4bHYi9+cr5/uau7IAsw1RBTplAACgUnRSy2F8EQAAoAYYKQMAAJVh9mV5jJQBAADUACNlmHGu2PPEtrJOgv9zYh9yA4ertvP3jipVL3eCQHrFgHapSQKpFQIiETc9mkren8r8n4qozomvz43BHx5rKxua0155dG0i0D81ISBxDqnAcbf86ztbg/o7MXdO3mcp9XOYOzEjXZbZwIy23bz/ceV2NitZo7lLg2AdXDUAAIAaYKQMAABUJiSNMeZTClcNAACgBhgpAwAAlWL2ZTl0yjArVBn8nxvUX3VyxKHcrPEZcicIlF0xQMpfNWBobnuAfcrY2sTAfcYlTgXcDyWC+lsD7htlqe0q/mNS4fd1Nkhl+e/ECy74cNZ2uRMCsn4O6W+gS+iUAQCAykQw+7IsrhoAAEAN0CkDAACoAW5fYsbpdqLYKuv16hhl49FSsWdX7Z0XZ3bL648tdUxJet55H20rm7N+ewLR1pifdPLQ9v1nxw8lksKOjCaOkZkoNhVrRLLY7rr1bz5Yuu72Fx3XVtb6c1h1DNygGCPwrhRGygAAAGqAkTIAAFCZxoLkjPmUwVUDAACoAUbKAABAhUiJURadMtTaLpce2VZWdW7PKqWCznshJ6g2dzJAJ8H/uVa88QNZ2237zY+s83ookey1EyNr25PdprzgE6vbym59/1PbyiwSxc4kBPGjbuiUAQCAyrAgeXl0ygAAwKxl+3RJ+0m6PyJ2KMrOlfT8YpOnS3ooIl6aqHuXpNWSRiWNRMT8braVThkAAKjUaA9yO07DGZK+KOms8YKIeNP4c9uflvTwJPVfFREPdK11TeiUAQCAWSsifmx7Xuo925b0Rkm797RRE6BThlq7cq8T2so6yeg/W7VOMHDin9TUigF1Cv5P+cXf/ts6r19wwYcr3f+cue2rCIw+3h78nwrqHxpun3SQWpUAGDQh9zpP2Sa2lze9PiUiTsms+78k3RcRt0/wfki61I0lPP5jGvsthU4ZAACYyR7oINbrQElnT/L+rhFxj+1nSLrM9q0R8eOSx5oSnTIAAFCpsRmQp8z2HEl/I2mnibaJiHuKr/fbvlDSAkld65TV/6oBAABU79WSbo2Ilak3bT/Z9lPHn0vaS9KN3WwQnTIAAFCZ8bUve/WYiu2zJV0p6fm2V9o+tHhrsVpuXdp+tu2Li5fPlPQT2z+TtEzSdyPikqquUwq3LzHjXLHniW1lucH/qcD21gB4J7ZJ/dinAudz5WTg78RQYmmBqoP/F1xydFvZsoUfy6pb1q1/88Gs7TqZEDA0tz2AP3XtUoH+qc8OgP6KiAMnKH9LouweSfsWz++Q9JKuNq4FI2UAAAA1wEgZAACoTMh1Sx47YzBSBgAAUAOMlAEAgEqxIHk5dMowK6SC/3PtcumRFbakXbeD+nOP2Unwf66XL20P/r9mn+4G/6fkTgh4fmJCwNz1R9rKnvOJ9rq/fn97WSrMf9tvfmSd162rFADAODplAACgMhHS6AxIHltHXDUAAIAaYKQMAABUyH0J25gNGCkDAACoAUbKUBs7f++ovhw3Few+KBKx/0mpa5Rbty7B/ym35a4QoPYJAcOJjP4pI6Pr/u+77fkfadvmFwfUO/h/u28d11Y2Mjq8zuu5c0az9nXTovZ9YXYJEVNWFlcNAACgBrrWKbN9uu37bbetqG77fbbD9iZNZUfZXmH7Ntt7N5XvZPuG4r2T7Mb/7LbXt31uUX617XndOhcAAJCvTguSzyTdPJszJC1sLbS9paQ9Jf26qWw7NVZr376oc7Lt8bHxL0k6TNI2xWN8n4dK+n1EPE/SZyWVT1QFAADQZ13rlEXEjyU9mHjrs5Ler3XzLC6SdE5EPB4Rd0paIWmB7c0kbRARV0ZESDpL0v5Ndc4snp8vaY/xUTQAANAfIWsseveYTXoa6G/7dZJ+ExE/a+k/bS7pqqbXK4uytcXz1vLxOndLUkSM2H5Y0saSHkgc9zA1Rtu01VZbVXIuwJDbI93rvAhvMst/Mgd9u9yg/tzj7rT0mFL7isT1vW7fj5ba13SkVgjY/qLjsurOaZkQsDaG27Z53nl55zA01ME3IlfyI9ze5lZr17Zvs/567asjvPjbeZMrcjnxc5iS+uz8/LXtEziAfupZp8z2kyQdI2mv1NuJspikfLI67YURp0g6RZLmz5/fg99qAAAMrtkW69UrvbxqfyFpa0k/s32XpC0kXWf7WWqMgG3ZtO0Wku4pyrdIlKu5ju05kp6m9O1SAACA2utZpywiboiIZ0TEvIiYp0anaseI+K2kJZIWFzMqt1YjoH9ZRNwrabXtnYt4sYMlXVTscomkQ4rnB0j6fhF3BgAAMON07fal7bMl7SZpE9srJR0bEaelto2Im2yfJ+lmSSOSDo+I8UyEb1djJucTJS0tHpJ0mqSv2l6hxgjZ4i6dCrqgk0SxV+398Ur3V2epOJiU3LiaVnVaCiX3XFvtePEHsrarOvYslQT1RUuObStrC0ROfKuS554oG03835kbZ1b2M5K9v8QprB1pjzNbb257nFlK6ppUHdT9ku9Ul7T3Z/u1JwUeVCFpjOSxpXStUxYRB07x/ryW18dLOj6x3XJJOyTKH5P0hs5aCQAAUA8sswQAACpkjdZoFH4mYXwRAACgBhgpAwAAlSGmrDw6ZahUbsB9JwG7Cy45uq1s2cLywf9lA8zrlEm67DnUyejY1OeQOs/W5KwTSU0ISH0Pr39N+YDtnM+EU4H5o+1FKZG4RqnrljpGar2T7GS0qckJLbenUhMJUmWjY/X5Y52aiDB3TuY3o0Vq0gDB/5guOmUAAKBSxJSVU59/WQAAAAYYI2UAAKAyESamrCSuGgAAQA0wUoau60VAfCr4P5VifKjirOaDIicIvxdSgeMjo+X/t0wFur/0u+0B26ng/xd/+4OpPbaXZHzmUue1zQdXt5XddtwGibpT7n5a26WC+lPhQa1tzt59Yv+ptuUumpd9Xpl1W4P/ywb+S4Md/D/KSFkpXDUAAIAaYKQMAABUJlSvdXVnEkbKAAAAaoCRMgAAUCETU1YSnTJUKjeoPxXE+9jauW1lQ0Pt2drX7yDwNtW+Ogf/pwLAu529vy5B/b0wlsqGn7jmL1pybFtZKjN9anWBst+vX3yoPag/mVk/UZadqT8lI6hfqvbnJjeoP7duJ8H/QD/RKQMAAJVprH1Jz7gMxhcBAABqgE4ZAABADXD7EgAAVGqUMZ9S6JShUqOJOILhREBwKhA3FUyc2m4kEWCdG8Cf2m694fITB9B/qc9NrlQQ/lji85U7+SHVljVr1/01mzpmMqN91hEn+llKlGUG/6d+bnKucSffh05UHdTfur/W758krTd3pNqDAgU6ZQAAoDIhE+hfEuOLAAAANcBIGQAAqNQYYz6lcNUAAABqgJEyVOrafY7P2m7BJUe3lT0hETz7WCLINvWvRCqrd+6CuGtGh9d5PZxYRSA1WSElFXQ8lGhcbtu6nb2/Fzo5hyqDx6u+lqns/SnDLds5MWlgbctncFoyA/OTZandZX/WuxvY34us/DmrAaRWR0hP1qjvyiC9FpGe9IWpMVIGAABQA4yUAQCASjH7shxGygAAAGqAkTIAAFCZRp4yxnzKoFOGvli28GNZ26UmBKSkAnZHE78UcjL/p4J4I5FfPZmFvYNY37oE9XdyXnU5h9x2pLZLZe8fTgR7pwPA248x1rK/fl2j3KD+fgTwp65bThD+dJT92UytcJDef+L3RqLsRUuObSu74XUfmn7DMCvRKQMAAJUazZxhjnUxvggAAFADjJQBAIDKhJh9WRadMtRaKvYsFWeWijUZViq5Z7nB4WSy11Qc21giji0Ro5SbeDTZlpK/7HJjY1LbjSTOq2o5sUydxIqlrls6likRUzY8mnXcsdT3v2V/kfg8pD6/Y7kxUDPsj18nsWKpz2buz0PZeLTk/lOfkUTS6dz9vfS7/zbtdo27/jUfKV0X9cPtSwAAgBpgpAwAAFSIlBhlcdUAAABqgJEyAABQqWQcLqZEpwwzTm7i2ZcvbZ8QMOz2YNw1o8PrvE4GBHcQJJ4KYB8ZbR+k7iQ4Pcd6c/OC1XMnBFQtGexeMvFqJzO/cidhpCZwrBkZbitrDewfTXzvq5ZqW4oTTUm1rmxC2U6SvaYMJYLpx0bbr3knyv581SVxMmY2OmUAAKAyEdIondRSiCkDAACzlu3Tbd9v+8amsuNs/8b29cVj3wnqLrR9m+0Vto/sdlsZKQMAAJWq2ezLMyR9UdJZLeWfjYhPTVTJ9rCkf5e0p6SVkq6xvSQibu5WQ2t11QAAAKoUET+W9GCJqgskrYiIOyJijaRzJC2qtHEtGCnDwGsN7E4F4ad0EtibGzidbEnJw+YGmKemA5QN9O6X1GSF0cxM+sOZgf7J1RsS17gt6D7zczOU2Y7U53BoqP38U+eV+r6W/VxXHdQ/J3MVhdR2axMTLrotNbmkF5NmUqsB9DvLf8i9XmZpE9vLm16fEhGnZNR7h+2DJS2X9N6I+H3L+5tLurvp9UpJf9lZUyfHSBkAAJjJHoiI+U2PnA7ZlyT9haSXSrpX0qcT2yTn4pdv5tQYKQMAAJWqe56yiLhv/Lnt/5T0ncRmKyVt2fR6C0n3dLNdjJQBAICBYnuzppevl3RjYrNrJG1je2vb60laLGlJN9vFSBkAAKhMqLMEzlWzfbak3dSIPVsp6VhJu9l+qRrNvUvSPxXbPlvSqRGxb0SM2H6HpO9JGpZ0ekTc1M220ikDWqQyuqeC/1PB1GWzgUv1Dqavum3J1QsystDnBqan/iAMZwa/p6TaljshJEvivFKB8+nzT+0vr27uMbotFRCfuzrGTJM61zVr2ycm5K7AgalFxIGJ4tMm2PYeSfs2vb5Y0sVdalobOmUAAKBSNctTNmNw1QAAAGqAThkAAEANcPsSAABUJ3qePHbWoFOGWeuafT7WVvbypUe3lQ215AJM5ddJBYmnfumkAqdTgb35v7Cqy7ieq1/B1GWzyyfrpbL3J76HKamg/sfWzM06RumJHsnPUmISQmZgfi8yyZfVSdvGUqsoZP4cdjIJJ0cn50VQP8bRKQMAAJUJ1T95bF0RUwYAAFADjJQBAIBKEVNWTtdGymyfbvt+2zc2lX3S9q22f277QttPb3rvKNsrbN9me++m8p1s31C8d5LdiBawvb7tc4vyq23P69a5AAAAdFs3R8rOkPRFSWc1lV0m6ahi6YITJR0l6Qjb26mxptT2kp4t6b9tbxsRo2qs5H6YpKvUyKq7UNJSSYdK+n1EPM/2YkknSnpTF88HA6I18F+SUsnmhzMDe1PB6rlBwam4jOys7jWWO5mgdbvcSQ7JCReZgf5r1rb/WhwbrfC//tS5Z65wkDqv1KoEudn7uy33c1715JLUcUczJwSU3X+u5GoTpfdWT3VbZmkm6dpIWUT8WNKDLWWXRsRI8fIqNVZcl6RFks6JiMcj4k5JKyQtKBYM3SAiroyIUKODt39TnTOL5+dL2mN8FA0AAGCm6WdM2T9IOrd4vrkanbRxK4uytcXz1vLxOndLUjHy9rCkjSU90Hog24epMdqmrbbaqrozAAAAbRgpK6cvsy9tHyNpRNLXxosSm8Uk5ZPVaS+MOCUi5kfE/E033XS6zQUAAOi6no+U2T5E0n6S9ihuSUqNEbAtmzbbQtI9RfkWifLmOittz5H0NLXcLgUAAL0VIqN/WT3tlNleKOkISa+MiD81vbVE0tdtf0aNQP9tJC2LiFHbq23vLOlqSQdL+kJTnUMkXSnpAEnfb+rkAUmpLP8LLlk3y3/qU5QM/s9MjthJEHNqMkEq2D2nLbPhp6PqgPA1a9tDrEdHq72B0DZZIZGVPhX8n4qQTa1KkLoi/cjo38kxU5/p3O917nFzV+Vor9c+kaKTVTXmzqk2e//1r/lIpftDf3WtU2b7bEm7SdrE9kpJx6ox23J9SZcVMflXRcTbIuIm2+dJulmN25qHFzMvJentaszkfKIasy6XFuWnSfqq7RVqjJAt7ta5AACAfGT0L6drnbKIODBRfNok2x8v6fhE+XJJOyTKH5P0hk7aCAAAUBdk9AcAANUJZl+WRacMaJGK5cmNM6tabvLYnCio3NsJsyH2LCUVezQnkXh1ZDQR35XK95qZ8LVtX4kYpeR2yTizRFnW3maeTuK2cqXixcq2I/lzmZmwGBjHguQAAAA1wEgZAACoDMsslcdIGQAAQA0wUgYAACrFSFk5dMow8JYtbE8o26o1wWyvDCWi7nOC/1PByblJNjvJL5Q7SSA1maLs/sdSyVgzZQdiJ9o7lAgSd+LEWr8XqckAnSTFTU8IKL27WstOKJv4vuYmBW49RmpfvUjOm+rUkCh29qNTBgAAKsMyS+URUwYAAFADjJQBAIBK9SLP3GzESBkAAEANMFIGZEhNBtj5e0dVeozcGIzUSgKtwfmpgOjc/1xzg5hHkwHrWVVrI3Wuc+eMtpWtHRnO2p8TQeFu+X4NJXaVumy5Gf1zpc61k89Jt2OGOmmbEj8jw4nVG3oRsN8q93s43Ie2VYkFycthpAwAAKAGGCkDAACVCRYkL42RMgAAgBpgpAwAAFSK2Zfl0CkDSrpq74+3lXUS/J+zsoAk7bT0mFL7zw0wzv1lOpwIaq/6lkXOCgGp80rdAkgF66eC+ucMt5eNJLLB566a0Nq+3Az82SswZF7z3IkDqbKR0fZr19q+TgLzU0H4qZUPcuv2I4A/pZOJGRhMdMoAAECFyOhfFjFlAAAANUCnDAAAoAa4fQkAACpFoH85dMqACqWC/6t27T7Ht5WVDf5P6SyTerVaA+BzAv8nkhvjkprAsF5mlv9U89oD4rOaka3qoPZU+1ITIlqP+/ja9j8nuasIrJeYXLHZ0x5pK3tkzfptZX94vL1sNHp/E4igflSBThkAAKhMiOSxZRFTBgAAUAOMlAEAgOpEZ6EGg4yRMgAAgBpgpAyYBVqD/6sM/J+OVGB3lbEl6SD5vIkJqbalgtPXnzvSVpYKdE8ZTWShb21LcgWCDoLE6xJg/oT11mZtN5T4fh0876q2sl8/vnFb2X2Pb9BW9sux9u0eXTu3rSz1mUhdu9REmhxV/8yVbUddjImYsjIYKQMAAKgBRsoAAEBlQuQpK4uRMgAAgBpgpAyYhVJxOyl1j/vIibUaGWv/3zI3zir1X2kncWbtkUzSY2uq+zWbe15Vx5lVmaA29Zn72q8XtJVd8eIL2so++sAL2spW/vFpbWWp7+E1+3Y3RqvbSZ1nlsFbkNz2SRmbPRIRH5hsAzplAAAAnVkk6YNTbHOkJDplAACgdwYwT9lnI+LMyTawveFUOyGmDAAAoAMR8Tnbw7b/ZbJtptoPnTIAAFCpCPfsURcRMarGbczSuH0JzAIvX3p0v5vQsbLB5HOGx0ofc2S0/f/S3OD/lNSEgCes117Wau3IcNb+q5ZKxpv7fahyMsGVe52QKG0v++CzKjtkT8z0BLAo5QrbX5R0rqQ/jhdGxHU5lemUAQAAVOOviq8fbioLSbvnVKZTBgAAKhMxuMljI+JVndQnpgwAAKACtp9p+zTbS4vX29k+NLc+nTIAAFCpsXDPHjVzhqTvSXp28foXkt6dW5nbl8AMs+CSckH9nWTvzw3qTt2ySAWOp36R5hwj95ZIal+puqlJAqng/+HE/kYTKwmkJgS0nn9yJYBEWSr4P3l9h/ImOqSC+gFUbpOIOM/2UZIUESO223/AJ0CnDAAAVGoAk8eO+6PtjdUI7pftnSU9nFuZThkAAEA13iNpiaS/sH2FpE0lvSG3Mp0yAABQqTrNvrR9uqT9JN0fETsUZZ+U9FpJayT9UtJbI+KhRN27JK2WNCppJCLmT3G4myS9UtLzJVnSbZpG/D6B/gAAYDY7Q9LClrLLJO0QES9WIxj/qEnqvyoiXprRIZOkKyNiJCJuiogbI2KtpCtzG8pIGVChXS49Mmu7dAbzdmWD+vsld0JAKnC+yv13UjcVdN/JJIHWuqkA/uEOViUA6iZUu+WPfmx7XkvZpU0vr5J0QCfHsP0sSZtLeqLtl0n/M7NqA0lPyt0PnTIAADCTbWJ7edPrUyLilGnU/wc1lkVKCUmXuvFf3X9Mst+9Jb1F0haSPtNUvlpS9n/XdMoAAEClejz58oHMW4ttbB8jaUTS1ybYZNeIuMf2MyRdZvvWiPhx60YRcaakM23/bUR8s0xbJDplAABgANk+RI0JAHtEpJN4RMQ9xdf7bV8oaYGktk6Z7YMi4r8kzbP9nsR+PtNalkKnDAAAVGcGrH1pe6GkIyS9MiL+NME2T5Y0FBGri+d7ad2Fxps9ufj6lE7aRacMKCk3qD/1y2nn70020adZvX+xDYrcFQKGhhKrF7TE8Nf9jxUw29g+W9JuasSerZR0rBqzLddX45akJF0VEW+z/WxJp0bEvpKeKenC4v05kr4eEZekjhER/1F8/VAnbaVTBgAAqlWjjP4RcWCi+LQJtr1H0r7F8zskvSTnGLZPmqIN/5yzHzplAAAAnbm2ip3QKQMAAOhAMftySra/EBHvnOh9OmUAAKBSxE5OaNfJ3qRTBgyIoUSQx1jNJxKk2txtnVyTodSEgHpf4iydrKTQ6tp9jq9sX8Bs07W1L22fbvt+2zc2lW1k+zLbtxdfN2x67yjbK2zfZnvvpvKdbN9QvHeSi2kQtte3fW5RfnXrEgoAAKA/Inr3mE26uSD5GWpfAPRISZdHxDaSLi9ey/Z2khZL2r6oc7Lt8QXiviTpMEnbFI/xfR4q6fcR8TxJn5V0YtfOBAAAoHOTjp13rVNWLEPwYEvxIknjwXBnStq/qfyciHg8Iu6UtELSAtubSdogIq4ssu2e1VJnfF/nS9pjfBQNAAD0R6gRU9arR13YHrb9ySk2+/xkb3ZzpCzlmRFxryQVX59RlG8u6e6m7VYWZZsXz1vL16kTESOSHpa0ceqgtg+zvdz28lWrVlV0KgAAAA0RMSppp8kGiCLijMn2UZdA/9QJxCTlk9VpL2ys6n6KJM2fP3+W3YHGbDGW+R9flTEU/Qikr7vcCRG5Wf6BgROaHTNcyvmppItsf0PSH8cLI+KCnMq97pTdZ3uziLi3uDV5f1G+UtKWTdttIemeonyLRHlznZW250h6mtpvlwIAAPTKRpJ+J2n3prKQVMtO2RJJh0g6ofh6UVP5121/RtKz1QjoXxYRo7ZX295Z0tWSDpb0hZZ9XSnpAEnfn2iVdwAA0DuD+tc4It7aSf1upsQ4W40O0/Ntr7R9qBqdsT1t3y5pz+K1IuImSedJulnSJZIOL+7NStLbJZ2qRvD/LyUtLcpPk7Sx7RWS3qNiJicAAEA/2N7W9uXj6cBsv9j2B3Lrd22kbIIFQCVpjwm2P15SW1bBiFguaYdE+WOS3tBJG4FOpBKFpuLCiD0aHKNjU39fO0nEmls3dx566jOce4yceMS6JydGFw3oSJmk/5T0r5L+Q5Ii4ue2vy7pozmVez37EgAAYLZ6UkQsaykbya1cl9mXAABgVqhX/rAee8D2X6gYK7R9gKR7cyvTKQMAAKjG4Wqk4HqB7d9IulPSm3Mr0ykDAADVGtCYsoi4Q9KrbT9Z0lBErJ5OfTplQIVyg/8xs6SC2teMDreVpQLsW1MDDA/l/bUa6yAIP1dqf2UTCl+zz8c6bQ4w49neWNKxkl4hKWz/RNKHI+J3OfUJ9AcAAKjGOZJWSfpbNXKorpJ0bm5lRsoAAEB1YqDT/mwUER9pev1R2/vnVmakDAAAoBo/sL3Y9lDxeKOk7+ZWZqQMAABUa0AD/SX9kxqrDH21eD0s6Y+23yMpImKDySrTKQNKumLPE9vKdr3siLayVPD/FXud0Fa24JKjs46bm629rEFds24ya8fygvqTGfJbAvtTt3VSKwH0Int/blA/QfxAnoh46mTv296+WFoyiduXAACgYu7hY0b56mRv0ikDAADojUl7kdy+BAAA1SIMYiKTXhlGygAAAGqAkTKgQqng/1zLFrYHU+cG/1ep6okEqYkDqcDxly/t/bnmyl2pISc4P10v75idyJ04QFA/KsFI2UTWTPYmnTIAAIAK2N4xUfywpF9FxEhE7DxZfTplAACgOiFpcDP6nyxpR0k/VyOof4fi+ca23xYRl05WmZgyAACAatwl6WURMT8idpL0Mkk3Snq1pE9MVZmRMgAAUKkBTkL9gubksBFxs+2XRcQdzgjYpVMGoDK5v4i7HdQ/1oOEkrmB+KNj696QGB4aa9um6sWbU0H9qez93V4dAhhAt9n+kqRzitdvkvQL2+tLWjtVZTplAACgWoM7UvYWSf9H0rvViCn7iaT3qdEhe9VUlemUAQAAVCAiHpX06eLR6g9T1SfQHwAAoAK297P9U9sP2n7E9mrbj+TWZ6QMAABUa3BTYnxO0t9IuiFi+tMd6JQBJe1y6ZFtZVUHbC9b+PG2siqz/M+GGVJPf83tbWUPfXebtrJU8H8q+L3qSQLXv+Yj67ze8eIPlN5Xblb+1HkB6Im7Jd1YpkMm0SkDAAAVq3iVsJnk/ZIutv0jSY+PF0bEZ3Iq0ykDAACoxvFqBPQ/QdJ6061MpwwAAFQnNMgpMTaKiL3KVmb2JQAAQDX+23bpThkjZai1XS87otL9jXV5RlAqELuT4P+dlh6TOkpbSZWB3b3Ihp9S9hxSQf0b/Vv7XYMHPtKeTDv3XKv8vl6370dL1ZPyV0JIZepftvBjpY8LTI8Hefbl4ZLeb/txNRLGWlJExAY5lemUAQAAVCAintpJfW5fAgCAakUPHzVg+1lVbEOnDAAAoDMXV7ENty9RG53Ej3U7VizXmtHhtrKqE5R2OzHobEg8+uBH1rSVbZRIMvvgd7ctfYzcRK5VumYf4sIwQ8z8XyPT9ZIpllOypCmXW6JTBgAA0IGIaP+PvAQ6ZQAAoFqDN1JWCWLKAAAAaoCRMgAAUJ3QIOcp6widMsw4uUH9V+51QlvZLpceWeqYqUSha8fyBpqHh9rH8a/Zuz1ge8EleYlBh0oGmOdet+jTbYfUdapSJ0H9AJDL9oaSni3pUUl3RcRYbl06ZQAAAB2w/TQ1svkfqMZC5KvUWJT8mbavknRyRPxgqv3QKQMAAJXqQ8aYfjtf0lmS/ldEPNT8hu2dJP297edGxGmT7YROGQAAQAciYs9J3rtW0rU5+6FTBgAAqjVgI2W2d5zs/Yi4Lmc/U3bKbC/J2M+DEfGWnAMC09FJpv7coP5UEH+r3KD+a/c5Pmu7VFB/KoB/eCg7PnQdqes2XKP7CTnXvGpzM69l7vc6105Lj5lym9zPDYDa+vQk74Wk3XN2kjNS9kJJ/zjJ+5b07zkHAwAAmIXeHBH3dLqTnE7ZMRHxo8k2sP2hThsCAAAwQ51WpML4oaRLJP0kIkamu5Mpx+kj4rwqtgEAAIPB0bvHlG2xT7d9v+0bm8o2sn2Z7duLrxtOUHeh7dtsr7A9YUxMROwjaTc1OmWvl3SV7QtsH2Z7q9zrlh08YXu+7QttX2f757ZvsP3z3PoAAAB9cIakhS1lR0q6PCK2kXR58XodtofVCM/aR9J2kg60vd1EB4mIxyLikoh4V0TMl/ReNe5IftH2spyGTmf25dck/aukGySViz4GJnHFnie2lZXNwD+RTjLzt0oFZ6famwrgH3L7McsG9aekjtnJpImquQ+TDnInF6QmBIyOtdcdU3sZAftAoUa/byLix7bntRQvUmNkS5LOVGOE64iWbRZIWhERd0iS7XOKejdnHvdOSSdLOtn2ejl1ptMpWxUROTMxAQAAemUT28ubXp8SEadMUeeZEXGvJEXEvbafkdhmc0l3N71eKekvJ9up7dVqTwjysKTltt873sGbyHQ6ZcfaPlWNYb7Hxwsj4oJp7AMAAMxmoV7nKXuguF1YtdRw31Rn9hlJ90j6elF/saRnSbpN0un68+hc0nQ6ZW+V9AJJc/Xn25chiU4ZAACYSe6zvVkxSraZpPsT26yUtGXT6y3U6HBNZmFENI+mnWL7qoj4sO32BJUtptMpe0lEvGga2wMAgEFUn1zVE1ki6RBJJxRfL0psc42kbWxvLek3aox6/d0U+x2z/UY11sKUpAOa3pvyqkynU3aV7e0iIivADajClXud0FbWSfB/TlB/KtD7qr0/ntWOVIB9Sj8C3auWe67dljuBIfeapyYEDA+1171m749l7Q9Af9k+W43bhpvYXinpWDU6Y+fZPlTSryW9odj22ZJOjYh9I2LE9jskfU/SsKTTI+KmKQ73ZkmfVyPAPyRdJekg20+U9I6p2jqdTtkrJB1i+041YsosKSLixdPYBwAAmOVq8j+bJCkiDpzgrT0S294jad+m1xdLungax7pD0msnePsnU9WfTqesNccHAAAACsWtzndKmqemPlZEvC6nfs6C5NdFxI4R8auptsk5YLH9v6ixnmaokffsrZKeJOlcNU7kLklvjIjfF9sfJelQSaOS/jkivleU76RGUrgnqtGTfVdE1Kh/DgDAABrcv8TfknSapG+rRE7XrAXJp8jcb0lPyz2g7c0l/bOk7SLiUdvnqRE8t50a2XVPKJYyOFLSEUX23MWStpf0bEn/bXvbiBiV9CVJh6lxz/ZiNUbzlua2BQAAoEKPRcRJZSvndMpekLHNaInjPtH2WjVGyO6RdJTS2XUXSTonIh6XdKftFZIW2L5L0gYRcaUk2T5L0v6iU4bCmtHhrO3KZmFPTULY9bLWhNBpnUwIyM1Mn3PM3CD5ukxgSJ17bts6WdEgNdEDABI+b/tYSZdq3Zyu1+VUnrJTNtltyzIi4je2P6XGbIdHJV0aEZfanii77uZqjISNW1mUrS2et5a3sX2YGiNq2mqr7HVBAQBAGYN7+/JFkv5e0u5aN6fr7jmVpxPoX4liJfZFkraW9JCkb9g+aLIqibKYpLy9sLHcwimSNH/+/MH9qAAAgG56vaTnRsSaMpV73imT9GpJd0bEKkmyfYGkv9LE2XUnyqi7snjeWg4AAPrEUa+UGD32M0lPV3qFgClNnUmzer+WtLPtJ9m2GnlCbtGfs+tK62bXXSJpse31i6mm20haVtzqXG1752I/ByudkRcAAKAXninpVtvfs71k/JFbuecjZRFxte3zJV0naUTST9W4tfgUJbLrRsRNxQzNm4vtDy9mXkrS2/XnlBhLRZD/QEgF2L986ZRLikkqH9TfC7lB8t0O/s/VSVB/2UkHnWTlzz3mbFhtAei7DibWzHDHdlK5H7cvFRHHqr3hjyuRXbfY/nhJbX9NI2K5pB0qbyAAAEAm246GH021zWT76cftSwAAMJtFDx/18APb77S9TooH2+vZ3t32mfpziNaE+jJSBgAAMIsslPQPks4u4t8fkvQENRYyv1TSZyPi+ql2QqcMlcpNnnrFnieWPkY/4sde8d/vz9wyL0YpFfM0OtY+cD081L5Kx09e/YnMtkwt9/uVG2fVSYxa2bq5cWFl4+4ATN+ghWZGxGOSTpZ0su25kjaR9GhEPDSd/dApAwAAqEhErJV0b5m6dMoAAEC1BmykrCoE+gMAANQAI2UAAKA6g53RvyN0ytAXu1x6ZNZ2I6Ptg7lOxGsvW/ixTps0bZ0EjicD0ROB/t3WyYQLAMC6bP+NpBMlPUONmV+WFBGxQU59OmUAAKBagztS9glJr42IW8pUJqYMAACgGveV7ZBJjJQBAABUZbntcyV9S43lIyVJEXFBTmU6ZQAAoFqDe/tyA0l/krRXU1lIolOGekhlXE9ZMzrcVjaU+MmeOzzaVtaLlQRy5Ga+T5mTOK+U1OoCVWb5BwCUExFv7aQ+MWUAAKBSjt496sT2trYvt31j8frFtj+QW59OGQAAQDX+U9JRktZKUkT8XNLi3Mp0ygAAAKrxpIhY1lI2kluZThkAAEA1HrD9FyqmOtg+QNNYnJxAf1QqN6g/lak/FdQ/PFQ+YKDO2eqHKg6E+OvL/3Wd1z/e45OV7h8ApqVmsV49dLikUyS9wPZvJN0p6aDcynTKAAAAKhARd0h6te0nSxqKiNXTqU+nDAAAVKeGsyJ7xfZ7Wl5L0sOSro2I66eqT0wZAABANeZLepukzYvHYZJ2k/SfttuTTLZgpAwAAFRrQEfKJG0saceI+IMk2T5W0vmS/lrStWosWD4hOmWo1JV7ndBW9vKlR2fVnTM8VnVzKlN1xvzWwPyqpfZP8D8AdN1WktY0vV4r6TkR8ajtxyeo8z/olAEAgGoN7kjZ1yVdZfui4vVrJZ1dBP7fPFVlOmUAAAAViIiP2L5Y0iskWdLbImJ58fabp6pPpwwAAFTGGujZl5+S9JWI+HyZ+sy+BAAAqMatkk6xfbXtt9l+2nQqM1KGvuhFUP+ulx3RVlaXLP+poHuC/wFgZouIUyWdavv5kt4q6ee2r5D0nxHxg6nqM1IGAACqFT181IztYUkvKB4PSPqZpPfYPmequoyUAQAAVMD2Z9SYcfl9SR+LiGXFWyfavm2q+nTKAABAdQZ4mSVJN0r6QET8KfHegqkqc/sSAACgGm7ukNkeLrL6KyIenqoyI2Xoi1Tm/1Rg/li49DE8QP+qDbUEVowp77oR/A+gKwbn12+rPWz/raRD1Vhy6SuSfpRbmU4ZAABABSLi72y/SdINkv4k6cCIuCK3PrcvAQBAtQZ09qXtbSS9S9I3Jd0l6e9tPym3Pp0yAACAanxb0gcj4p8kvVLS7ZKuya3M7UvURiqx6y6XHtmHllRrt8vfl7Vd6j+kVGxYa/xY1VLt/eEen+rqMQHMLgMU0ttqQUQ8IkkREZI+bXtJbmVGygAAAKqxge0Lba+yfZ/tb0p6NLcynTIAAFCtAY0pU2O25RJJm0naXI3bmV/JrUynDAAAoBqbRsRXImKkeJwh6Rm5lemUAQCA6vRylKx+I2UP2D6oSBo7bPsgNda/zEKgP/qi6gD+uiSKzQ3qz1U2qD+VADaVKDa1/1RQ/6u+/962sh/s/ulSbQOAWewfJH1R0mfV6DL+v6IsC50yAABQqZr8n9xzEfFrSa8rW5/blwAAADVApwwAAKAGuH0JAACqNaC3LztFpwyVevnSo7O2G8oMOBiL9oz2uVLHSAW7p4Lic1Qd1N+JnIkOueeZCurP3Y7gfwBYl+0dI+K6nG25fQkAACrl6N1jBnh77oZ0ygAAALokIv537rZ0ygAAQLVqkjzW9vNtX9/0eMT2u1u22c32w03bfLDsadve1faTi+cH2f6M7efk1iemDAAAzEoRcZukl0qS7WFJv5F0YWLT/xsR+1VwyC9Jeontl0h6v6TTJJ0l6ZU5lemUoS+u2PPErO1Smf/7kb2/TkH9dUbwP4CaLn8kSXtI+mVE/KqLxxiJiLC9SNLnI+I024fkVub2JQAAmMk2sb286XHYBNstlnT2BO/tYvtntpfa3r6Dtqy2fZSkgyR9txidm5tbmZEyAABQGRePHnogIuZPtoHt9dRY/uioxNvXSXpORPzB9r6SviVpm5JteZOkv5N0aET81vZWkrLzLjFSBgAAZrt9JF0XEfe1vhERj0TEH4rnF0uaa3uTksdZrcZty/9re1s14tkmGp1r05dOme2n2z7f9q22b7G9i+2NbF9m+/bi64ZN2x9le4Xt22zv3VS+k+0bivdOst3jzjkAAGhTk9mXTQ7UBJ0j288a7z/YXqBG3+h3+Se7jh9LWt/25pIul/RWSWfkVu7X7cvPS7okIg4ohhSfJOloSZdHxAm2j5R0pKQjbG+nxn3g7SU9W9J/2942IkbVmOVwmKSrJF0saaGkpb0/HYwbSwxaX7vP8W1lr/jv92fusfv/N7Rm+R+qOEI1d2JCZK5e0O2JDqnA/Nws/ympuqlz/eEenyp9DACYiO0nSdpT0j81lb1NkiLiy5IOkPR22yOSHpW0OCLK/qJ1RPzJ9qGSvhARn7B9fW7lnnfKbG8g6a8lvUWSImKNpDXFTIXdis3OlPRDSUdIWiTpnIh4XNKdtldIWmD7LkkbRMSVxX7PkrS/6JQBANBXdcq0HxF/krRxS9mXm55/UdIXKzqcbe8i6c2SDi3KhnMr9+P25XMlrZL0Fds/tX1qkWjtmRFxryQVX59RbL+5pLub6q8syjYvnreWt7F92PisjFWrVlV7NgAAAA3vVmMywYURcZPt50r6QW7lfnTK5kjaUdKXIuJlkv6oxq3KiaTu6cQk5e2FEadExPyImL/ppptOt70AAGA66hdT1hMR8aOIeF1EnFi8viMi/jm3fj9iylZKWhkRVxevz1ejU3af7c0i4l7bm0m6v2n7LZvqbyHpnqJ8i0Q5AABAz9j+XES82/a3legqRsTrcvbT805ZkbfjbtvPL5Y/2EPSzcXjEEknFF8vKqoskfR1259RI9B/G0nLImLU9mrbO0u6WtLBkr7Q49NBjY0lgsmHuhzo0I/VBlI6CczvhdSEEACYwb5afO1oxlK/Zl++U9LXipmXd6gxZXRI0nnFjIVfS3qDJBX3ZM9To9M2IunwYualJL1djammT1QjwJ8gfwAA+q0e/5/2TERcW3z9ke1Ni+fTDmLvS6csIq6XlMq+u8cE2x8vqS2vQkQsl7RDpY0DAACYhiLP2bGS3qFGzPtQkWLjCxHx4dz9kNEfAABUJxopMXr1qIl3S9pV0ssjYuOI2FDSX0ra1fa/5O6EThkAAEBnDpZ0YETcOV4QEXeosTD5wbk7YUFy9EUqo3svguSrDP6vur2p/eVm+Z9pql41AUDNDN6P+NyIeKC1MCJW2Z6buxNGygAAADqzpuR762CkDAAAVKpGsV698hLbjyTKLekJuTuhUwYAANCBiMhe33IydMqAFqnEpjMtBqrqJLmpWLxO/HCPjvIrAqi7mfUrszaIKQMAAKgBRsoAAEClBjCmrBKMlAEAANQAI2UAAKA6IWLKSqJThlqrOmC9SrmJXTtJMtvthLpVB/Cn/GD3T3f9GAAwG9ApAwAA1arv/9O1RkwZAABADdApAwAAqAFuXwIAgMpYpMQoi04ZKpUb/F51XUytF0H9AIDy6JQBAIBqMVJWCjFlAAAANcBIGQAAqJSDobIyGCkDAACoAUbKUGup4PThPkzrGVNekPxQnwIpqlz5YDTa/1cb9lhl+wcwy7HMUmmMlAEAANQAI2UAAKBS5Ckrh5EyAACAGmCkDAAAVIuRslLolKHWhofaA8w7yfzvxJh6Tqb7KgPpO9XttnQS1P+D3T9dYUsAYLDQKQMAAJWq0f+xMwoxZQAAADXASBkAAKgWI2WlMFIGAABQA4yUYVbICdaX8v4LSU0GSO9rsP8VJKgfAKpFpwwAAFQnCPQvi9uXAAAANcBIGQAAqBYjZaXQKUOtdZIoNlU3lRa1Tolhc+TGz1WJ+DEA6D46ZQAAoDIWMWVlEVMGAABQA4yUAQCAagVDZWUwUgYAAFADjJShL67Y88Ss7Xa97Ii2slSge+6EgJzg/+GKgyFSx8xNUFv2GFXvHwCmg19B5TBSBgAAUAOMlAEAgOqEyFNWEiNlAAAANcBIGQAAqJRTmboxJTplqFTV2eb7kb0+dyUAZTZtKDGO38lKBakg/pFYd9B7TqLFucH/ZO8HgP6gUwYAAKpFTFkpxJQBAADUAJ0yAACAGuD2JQAAqBTJY8uhU4ZKDfXgJzEZiN9B4HzrcHEqqD81pJw9uaiDCQEpqfOfOzSa2xoAQE3RKQMAANUJsSB5ScSUAQAA1AAjZQAAoFJ1iimzfZek1ZJGJY1ExPyW9y3p85L2lfQnSW+JiOt63U6JThkAAJj9XhURD0zw3j6StikefynpS8XXnqNThlobHSt/hz0V0uBE0H3rJIHhzH/xUgH3qaz5P97jk1n7S9nt8vdlHQMAamVm/ZpaJOmsiAhJV9l+uu3NIuLeXjekbzFltodt/9T2d4rXG9m+zPbtxdcNm7Y9yvYK27fZ3rupfCfbNxTvnVQMQQIAgMGxie3lTY/DWt4PSZfavjbxniRtLunuptcri7Ke62eg/7sk3dL0+khJl0fENpIuL17L9naSFkvaXtJCSSfbHi7qfEnSYfrzsOPC3jQdAACkWI2Ysl49JD0QEfObHqe0NGnXiNhRjduUh9v+60STW/VlrK8vnTLbW0h6jaRTm4oXSTqzeH6mpP2bys+JiMcj4k5JKyQtsL2ZpA0i4spiyPGspjoAAACKiHuKr/dLulDSgpZNVkrasun1FpLu6U3r1tWvkbLPSXq/1s2/+czx+7fF12cU5RMNK25ePG8tb2P7sPFhzVWrVlVyAgAAICGit49J2H6y7aeOP5e0l6QbWzZbIulgN+ws6eF+xJNJfQj0t72fpPsj4lrbu+VUSZTFJOXthY2hzFMkaf78+TMr/HCG6SSzfspoYn+5me9TcoL/UwH8qfndqYD7n7z6E2WblvTDPT7VVvaq77+3sv3/YPdPV7YvAKihZ0q6sAg5nyPp6xFxie23SVJEfFnSxWqkw1ihRkqMt/aprX2ZfbmrpNfZ3lfSEyRtYPu/JN03PtuhuDV5f7H9RMOKK4vnreUAAKCP6jJJPCLukPSSRPmXm56HpMN72a6J9Pz2ZUQcFRFbRMQ8NQL4vx8RB6kxfHhIsdkhki4qni+RtNj2+ra3ViOgf1kxtLja9s7FrMuDm+oAAADMKHXKU3aCpPNsHyrp15LeIEkRcZPt8yTdLGlE0uERMb768tslnSHpiZKWFg8AANBPNRkpm2n62imLiB9K+mHx/HeS9phgu+MlHZ8oXy5ph+61EAAAoDfqNFIGZBlLzPGoNPg/c65C1UH9uQjOB4DZiU4ZAACoVF0C/Weafmb0BwAAQIGRMgAAUJ2QNMZQWRl0ylBrqUSuqaStncSZtSaPTe3/ij1PzNoXAABl0SkDAADVYqCsFGLKAAAAaoCRMgAAUClmX5bDSBkAAEANMFKG2tjx4g9Uur+qk8wCADK1ZeVGDkbKAAAAaoCRMgAAUCliysphpAwAAKAGGCkDAADVCZGnrCQ6ZZi1COoHAMwkdMoAAEBlLMnMviyFmDIAAIAaoFMGAABQA9y+BAAA1RrrdwNmJjplqNQQyWkAACiFThkAAKgUgf7lEFMGAABQA4yUAQCA6pA8tjRGygAAAGqAkTJUyolA/52WHtNWdu0+x5c+Rid1AQDdFhIxZaUwUgYAAFADjJQBAIBKkR2pHEbKAAAAaoCRMgAAUC1iykqhU4auG2JuNAAAU6JTBgAAqhOSWfuyFGLKAAAAaoCRMgAAUC1iykphpAwAAKAGGClDpVLZ9l++9OisuqnVAMjeDwAzEANlpTBSBgAAUAN0ygAAAGqA25cAAKBSJtC/FEbKAAAAaoCRMgAAUC1GykphpAwAAKAGGCkDAADVCUkss1QKI2UAAAA1wEgZKrXz945qK7Pdh5YAAPrBCmZflsRIGQAAQA0wUgYAAKrFSFkpjJQBAADUACNlAACgWoyUlUKnDF23bOHH+t0EAABqj04ZAACoDnnKSiOmDAAAoAYYKQMAAJUiT1k5jJQBAADUACNlKG2XS49MlJK9HwCAMuiUAQCAanH7shRuXwIAANRAzztltre0/QPbt9i+yfa7ivKNbF9m+/bi64ZNdY6yvcL2bbb3birfyfYNxXsnmZWvAQDos2iMlPXqMYmJ+hwt2+xm+2Hb1xePD3bt0kyhHyNlI5LeGxEvlLSzpMNtbyfpSEmXR8Q2ki4vXqt4b7Gk7SUtlHSy7eFiX1+SdJikbYrHwl6eCAAAqLWJ+hyt/m9EvLR4fLi3TfyznseURcS9ku4tnq+2fYukzSUtkrRbsdmZkn4o6Yii/JyIeFzSnbZXSFpg+y5JG0TElZJk+yxJ+0ta2qtzGTQ7f++odV6nxiWv2vvjWfva8eIPtJWl9vfypUe3lV2zDysEAEBthWoTUzZJn+PmvjZsAn2NKbM9T9LLJF0t6ZnFxRu/iM8oNttc0t1N1VYWZZsXz1vLU8c5zPZy28tXrVpV6TkAAIC+2mT8b3zxOCy1UUufo9Uutn9me6nt7bvZ2Mn0bfal7adI+qakd0fEI5OEg6XeiEnK2wsjTpF0iiTNnz+/Ht13AABmq94us/RARMyfbIPWPkfL29dJek5E/MH2vpK+pUZIVM/1ZaTM9lw1Ls7XIuKCovg+25sV728m6f6ifKWkLZuqbyHpnqJ8i0Q5AACApAn7HP8jIh6JiD8Uzy+WNNf2Jj1upqT+zL60pNMk3RIRn2l6a4mkQ4rnh0i6qKl8se31bW+tRu91WXGLc7XtnYt9HtxUBwAA9IkjevaYtB0T9zmat3nWePYG2wvU6Bv9ruJLkqUfty93lfT3km6wfX1RdrSkEySdZ/tQSb+W9AZJioibbJ+nRlDeiKTDI2K0qPd2SWdIeqIaAf4E+Vdkp6XHtJXNbenCX7nXCVn7SgX1dyIV/J8ylrjDHdFeNtZSlvoZ/9l+H8lrHACgTibqc2wlSRHxZUkHSHq77RFJj0paHNGfmQr9mH35E028Fs8eE9Q5XtLxifLlknaornUAAKBj9Zl9OVmfY3ybL0r6Ym9aNDky+gMAANQAa18CAIDqhKSxeoyUzTSMlAEAANQAI2XQS77zb21ldnt//dp92sL62uQG9V+370eztuuHFy05tt9NAIAZbOo1KZHGSBkAAEAN0CkDAACoAW5fAgCAanH7shQ6ZbNYKjbKTv2gtKdw+dl+H55y/6kEs6l91Tl+LOWG132o300AAAwgOmUAAKBajJSVQkwZAABADTBSBgAAqkPy2NIYKQMAAKgBRspmoO0vOi5ru6FJl2D9s5+/duqg/uT+1f6f0DX7Tp1gFgAwm4UUY/1uxIzESBkAAEANMFIGAACqxezLUhgpAwAAqAFGygAAQHWYfVkanbKaS2XlT2XNv2nRcVn7e8l3/q2zBjW5Zp+PVbYvAAAGHZ0yAABQLWLKSiGmDAAAoAYYKQMAANVipKwURsoAAABqgJGyGtnuW8clSvPS8qey/KeC/3+230em1SYAANAbdMoAAECFgtuXJXH7EgAAoAYYKQMAANUJSWMsSF4GI2UAAAA1wEhZn7zwwg8lStuD+m95fSqjPwAANUZMWSmMlAEAANQAI2UAAKBajJSVwkgZAABADTBSBgAAKhTSGCNlZdAp6xO7/QN78/7H9b4hAACgFuiUAQCA6oQUQZ6yMogpAwAAqAFGygAAQLWIKSuFkTIAAIAaYKSsT9yevB8AgNmBPGWlMFIGAABQA3TKAAAAaoDblwAAoDoR0hgpMcqgU9YnNy06rt9NAAAANUKnDAAAVItA/1KIKQMAAKgBRsoAAEClgpiyUhgpAwAAqAFGygAAQIWCmLKSGCkDAACoAUbKAABAdUIsSF4SI2UAAAA1wEgZAACoVjD7sgxGygAAAGqAkTIAAFCZkBTElJXCSBkAAEANMFIGAACqE0FMWUkzfqTM9kLbt9leYfvIfrcHAACgjBndKbM9LOnfJe0jaTtJB9rerr+tAgAAmL6ZfvtygaQVEXGHJNk+R9IiSTf3tVUAAAwwAv3LmdEjZZI2l3R30+uVRdk6bB9me7nt5atWrepZ4wAAAHLN9E6ZE2Vt3fOIOCUi5kfE/E033bQHzQIAYIDFWO8es8hM75StlLRl0+stJN3Tp7YAAACU5oiZe9/X9hxJv5C0h6TfSLpG0t9FxE2T1Fkl6VfFy00kPdDtdtYc14BrMI7rwDUYx3WYXdfgORHRs9tEti9R4/r1ygMRsbCHx+uaGd0pkyTb+0r6nKRhSadHxPHTqLs8IuZ3q20zAdeAazCO68A1GMd14BqgP2b67EtFxMWSLu53OwAAADox02PKAAAAZoVB75Sd0u8G1ADXgGswjuvANRjHdeAaoA9mfEwZAADAbDDoI2UAAAC1QKcMAACgBgayU2Z7oe3bbK+wfWS/21Ml21va/oHtW2zfZPtdRflGti+zfXvxdcOmOkcV1+I223s3le9k+4bivZNsp1ZQqC3bw7Z/avs7xetBvAZPt32+7VuLz8Qug3YdbP9L8bNwo+2zbT9hEK6B7dNt32/7xqayys7b9vq2zy3Kr7Y9r6cnmGGCa/DJ4ufh57YvtP30pvdm3TXADBMRA/VQI5/ZLyU9V9J6kn4mabt+t6vC89tM0o7F86eqkVx3O0mfkHRkUX6kpBOL59sV12B9SVsX12a4eG+ZpF3UWM5qqaR9+n1+07wW75H0dUnfKV4P4jU4U9I/Fs/Xk/T0QboOaqyFe6ekJxavz5P0lkG4BpL+WtKOkm5sKqvsvCX9H0lfLp4vlnRuv8858xrsJWlO8fzE2X4NeMysxyCOlC2QtCIi7oiINZLOkbSoz22qTETcGxHXFc9XS7pFjT9Mi9T4A63i6/7F80WSzomIxyPiTkkrJC2wvZmkDSLiyogISWc11ak921tIeo2kU5uKB+0abKDGH6XTJCki1kTEQxqw66BGPsYnurECyJPUWIpt1l+DiPixpAdbiqs87+Z9nS9pj7qNHqauQURcGhEjxcur1FieT5ql1wAzyyB2yjaXdHfT65VF2axTDKW/TNLVkp4ZEfdKjY6bpGcUm010PTYvnreWzxSfk/R+Sc2r1Q7aNXiupFWSvlLcxj3V9pM1QNchIn4j6VOSfi3pXkkPR8SlGqBr0KLK8/6fOkUn52FJG3et5d3xD2qMfEmDew1QI4PYKUv9FzPr8oLYfoqkb0p6d0Q8MtmmibKYpLz2bO8n6f6IuDa3SqJsRl+Dwhw1bt18KSJeJumPatyymsisuw5FzNQiNW5HPVvSk20fNFmVRNmMvgaZypz3jL4mto+RNCLpa+NFic1m9TVA/Qxip2ylpC2bXm+hxu2MWcP2XDU6ZF+LiAuK4vuKYXgVX+8vyie6Hiv152H95vKZYFdJr7N9lxq3p3e3/V8arGsgNdq/MiKuLl6fr0YnbZCuw6sl3RkRqyJiraQLJP2VBusaNKvyvP+nTnFr+Glqv11aS7YPkbSfpDcXtySlAbsGqKdB7JRdI2kb21vbXk+N4MwlfW5TZYp4htMk3RIRn2l6a4mkQ4rnh0i6qKl8cTGLaGtJ20haVtzaWG1752KfBzfVqbWIOCoitoiIeWp8f78fEQdpgK6BJEXEbyXdbfv5RdEekm7WYF2HX0va2faTirbvoUac5SBdg2ZVnnfzvg5Q4+es9qNEthdKOkLS6yLiT01vDcw1QI31e6ZBPx6S9lVjVuIvJR3T7/ZUfG6vUGP4/OeSri8e+6oR53C5pNuLrxs11TmmuBa3qWlGmaT5km4s3vuiihUgZtJD0m768+zLgbsGkl4qaXnxefiWpA0H7TpI+pCkW4v2f1WN2XWz/hpIOluNOLq1aozoHFrleUt6gqRvqBEQv0zSc/t9zpnXYIUacWDjvx+/PJuvAY+Z9WCZJQAAgBoYxNuXAAAAtUOnDAAAoAbolAEAANQAnTIAAIAaoFMGAABQA3TKAAAAaoBOGYBK2J5n+1Hb10+jzptsr7D9nS42DQBmBDplAKr0y4h4ae7GEXGupH/sXnMAYOagUwagcsWo2a22T7V9o+2v2X617Sts3257Qb/bCAB1Q6cMQLc8T9LnJb1Y0gsk/Z0ay4C9T9LRfWwXANQSnTIA3XJnRNwQEWOSbpJ0eTTWdbtB0ry+tgwAaohOGYBuebzp+VjT6zFJc3rfHACoNzplAAAANUCnDAAAoAa4hQCgchFxl6Qdml6/ZaL3AAANjJQBqMqopKdNN3mspJMl/b5bjQKAmcKNyVAAAADoJ0bKAAAAaoBOGQAAQA3QKQMAAKgBOmUAAAA18P8B29/1vYq0chQAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 720x720 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "dsb[0].plot(figsize=(10,10))\n",
    "plt.title(\"Bottom oxygen\");"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "11.076560020446777"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dsb[0].to_numpy()[0,110,56]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "12.409002304077148"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dst[0].to_numpy()[0,110,56]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.05738005042076111"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dsb[0].to_numpy()[0,58,52]"
   ]
  }
 ],
 "metadata": {
  "interpreter": {
   "hash": "f4041ee05ab07c15354d6207e763f17a216c3f5ccf08906343c2b4fd3fa7a6fb"
  },
  "kernelspec": {
   "display_name": "Python 3.9.6 64-bit",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.10"
  },
  "orig_nbformat": 4
 },
 "nbformat": 4,
 "nbformat_minor": 2
}