dwhswenson/contact_map

View on GitHub
examples/advanced_matplotlib.ipynb

Summary

Maintainability
Test Coverage
{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "2b0b99b2",
   "metadata": {},
   "source": [
    "# Advanced matplotlib tricks\n",
    "\n",
    "The most common tricks to get the publication-ready figures for your contact map plots are described in [Customizing contact map plots](custom_plotting.ipynb). Here we will illustrate a few more advanced techniques. In general, we recommend getting familiar with [matplotlib](https://matplotlib.org/) through its own documentation, but some of the recipes below may be useful."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "c88acc64",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import matplotlib.pyplot as plt\n",
    "import mdtraj as md\n",
    "traj = md.load(\"5550217/kras.xtc\", top=\"5550217/kras.pdb\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "344ebc78",
   "metadata": {},
   "outputs": [],
   "source": [
    "from contact_map import ContactFrequency\n",
    "traj_contacts = ContactFrequency(traj)\n",
    "frame_contacts = ContactFrequency(traj[0])\n",
    "diff = traj_contacts - frame_contacts"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "4b8c5e7f",
   "metadata": {},
   "source": [
    "## Advanced color schemes\n",
    "\n",
    "In some cases, matplotlib's built-in color maps may not be sufficient. You may want to create custom color schemes, either as a matter of personal style or to create schemes that place more emphasis on certain values. To do this, the `diverging_cmap` keyword in Contact Map Explorer's plotting function can be useful, as can some matplotlib techniques."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "60de0fc0",
   "metadata": {},
   "source": [
    "### Customizing whether a color map is treated as diverging or sequential\n",
    "\n",
    "As discussed in [Changing the color map](custom_plotting.ipynb#Changing-the-color-map), Contact Map Explorer tries to be smart about how it treats diverging color maps: if the data includes negative values (as possible with a contact difference), then the color map spans the values from -1 to 1. On the other hand, if the data only includes positive values and if the color map is diverging, then only the upper half of of the color map is used.\n",
    "\n",
    "The diverging color maps that Contact Map Explorer recognizes are the ones listed as \"Diverging\" [in the matplotlib documentation](https://matplotlib.org/stable/tutorials/colors/colormaps.html). The sequential color maps recognized by Contact Map Explorer are the ones listed as \"Perceptually Uniform Sequential\", \"Sequential\", or \"Sequential (2)\". Other color maps are not recognized by Contact Map Explorer, and by default will be treated as sequential while raising a warning:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "bd215a01",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/Users/dwhs/Dropbox/pysrc/contact_map/contact_map/plot_utils.py:31: UserWarning: Unknown colormap: Treating as sequential.\n",
      "  warnings.warn(\"Unknown colormap: Treating as sequential.\")\n"
     ]
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYkAAAD8CAYAAACCRVh7AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABMHElEQVR4nO29eZwdVZnw/32ysjUJIaEJAUIgAUmCkMAAM0LeDIrEIAICgoqyDZgRGH1FBdR5oyMqyoAyoCJLWEQWWQUJIg4/DKhhC1sWIwmQkIWEEJI0SehOp5/fH6eqb93qqrp1763t3j7f/tSn61adOufculXnOed5nvMcUVUsFovFYgmiT94VsFgsFktxsULCYrFYLKFYIWGxWCyWUKyQsFgsFksoVkhYLBaLJRQrJCwWi8USihUSFovF0iSIyAwRWS0ic0POi4j8j4gsEpFXRGRipTytkLBYLJbm4RZgSsT5TwBjnO084JeVMrRCwmKxWJoEVZ0FrI1IcjxwmxpmA4NFZHhUnv2SrGCt9JW+2p9+tNNR0/X9HFnXSVeS1UqcfvShbx3fs9YyXYp2f/rRJ7JOjfK7hpHl96tUVlL4y+nn62e65yqlC7qm2nr09TRfUe9UVvcGWKOqw6q9aMoU0TVr4qV94QXmAR94Dl2vqtdXUdwI4C3P52XOsZVhFxRCSHTRVVfD2SiNSCdd7MxAVsX4rq20sIq2RMqsVI6Lt7yw49US9T0q1a3W39WtexL3rx7C6j+akSxiSaLPbVbvgL8c72f3ewHszPZl6VoY0n0Oyn8j77MWlM79Hb37nXTRSQejGQlQdo23HkF1TpEllZP0ZM0aeO6ZeGn79OMDVT2klnIcJOBYZGymQgiJ3kTchqvehjluHmFpkmpgvY1Akt8pKs+8hUMlFtXWlhQCv9D3NvBtrO3+3MIQ2hyth7vvntuZkcynZFf1CgZ/3t7PQb9rLfcyqQ5YYiiwNbPSlgF7eD7vDqyIusAKiSakUC8AydTHn0dW37GaBqVwjU/KuL14KDXWriDwNt4tDAG8v1l5w76IJYxlPADves75RxlR97aVlu5yKgmOwv1GCpqdkHgIuEBE7gIOA9araqiqCayQyJSiqEEs8anmt2qW3zVK/eM/744WvMeCPrcwpPv+jGACq5hVdt4vHNw8vPfUq0Yazcjusr3XBF3XECQkJETkTmAyMFRElgHTgf4AqnodMBOYCiwCNgFnVcyzCKHCRUQh/x83bvl519PS2BT5+alUN7++P44aMChPd+Qwn7kV7WJRx7yjGZcc1Xkv1GIvOOQg0WeeiJe23861lVEPhRpJ5P3iVCo/y5FAkRuSRqE33MOkn0l/Pv57uIglZWX6e/2uysdLC0PKbBRAmU3CW65/BBJWF1e91Mj2HRcFtMC+N4USEkUn7EVMozFq9sYtC5K6h0n/vkXNK07+Xt2//3wrLYGjDFc4RAkCr5dSJUERVC9/vg31/mRruK4aKyQSoKEeyBwJ8oxphHvXCHVMiyAPo6DRgnsuaJQRh7aQ+V+raCsTIO6+fxThV4M1GhkarqvGCglLIFmMjhqh8W0UQZYW/u8e1RhH3SvvdUFu0UGjDPeY11jt7rs2Da83ld/9tpI3VND3ywU7kqjMQAawB8Mr/rCNQNEalVrrU6TvkCeNeh+Seg79jXZYj9897zUkh6X15hlkePbbHdx8vILAtWn0dK+lIdsRO5KoQDsdDT1U9FKkhzNKt2upjjg9zyJ1ENIyZEPwbGl/Om9D7bdV+FVIfoO0N503rV+FNZFJzPG50sahSL8TYOZJbMm7EuEUQkg0On7/7aIIvEK9CAEU7mWNoIiT95K4zns9VBYKYbOtg7yNgkJmuARNtnP33c9trC0bSfiFiVdANPpkOqtuanK8D6W/Z9SoZPEdGv0e+UljZnlWZQd5MbnHw/KO8mzyT3bz4u1IBY0yKtUtrC6N/DxZF9heRiM/rC7N8B0stRPl7u3tsYcF6vPmUanDETT/wc3fy1jGd9siqp28F5dc3GjtSKJ3UK+3RKG8LRqYhvWVz5FqG1XvCCHIcOxPA+WjB//IoZWWsuffzdNfL+8EvChVVq1Uex+Semc129hNVVNRSIjIHsBtwK5AFyZ++dUiMgS4G9gLeBP4jKq+51xzKXAORj7+h6o+VkvlGkltU6+7XZzvaQVJZey9qZ5qjNP+EUOYqmg0I0NnavvThuXZSgtjGU+Hz3PJ/e81aofFcaqGaq9L7FlrAsN1J3CRqs4RkRbgBRF5HDgT+F9VvVxELgEuAS4WkbHAacA4YDfgTyKyr2r1sjKPFz6uAa8a0vQ0sVjSIExV5LcnhD2T3tDg/uvca91ywvIxQmluj2u89fPP5q5Ur8LSyDYJJ4zsSme/TUQWYFYyOh4TbRDgVuBJ4GLn+F2q2g68ISKLgEOBvyVd+TSot7fvYnv9+ZHmvQ9qgJr5t/Z/J/+IIMzN2js68Mdt8hq13fvpvYdhqqSg++4e80/Wy5q6BFOjq5u8iMhewATgGaDVjUOuqitFZBcn2Qhgtucyd3k8f17nYRbibhjCGoOgB8TaJvIjzXtXq/dNs+A3MscxTvvDiUfZEeIIIOg5hyJvt/O6n4FGHkm4iMgOwH3AV1V1g0jQKngmacCxHvHInXVZr3fyzj9eeQziuALWQ0MOkwtMo9zPvOZG+POCeB2guGnD5lQErQURNPFuLON5lyWRgihoDkW9ZN5Ra4aRhIj0xwiI36jq/c7hVSIy3BlFDAdWO8erXh6vEfE+6N6HyeumVy3+YXc1VPOSh13fCA1qHBptpnlecyPi5LWKNo5iGgDrmN9jhnM1z6ubJm69w8KJVzoWRJirbpAg83ps+b9fkFoMEhjJNPJIQsyQ4SZggape5Tn1EHAGcLnz/3ee43eIyFUYw/UY4NkkK50H/gcryAgHpRW20jByR+VZ7yinGQSEK7i9M3Yt1eG/hwBv8zQAHQFG5rijjLjzJNJ6DqNGI34PqahJgO65RCMrNIF300eALwCvishLzrFvYYTDb0XkHGApcAqAqs4Tkd8C8zGeUefX4tlULUFeGJCdrtJrQPMOn6ulFkEQxefkEe7QY2uqS6ORp+GyGrIYtdVahv8ettLCaX1eAOCprk9XjABbS5mraCtb4zqJexNm1A4iSCMQdl0av1vR50kUbvnScXy++/jbPB1bdZP0i+fPzz+SaAT1TJJ1dFUOT3BdIvnVS9yAey4N5/GSI1EB9xqFsJnhbnypsOcnTI0cdc7TIa1padEJY0T//LN4aQd9MvvlSwslJIpO0MQd6KnXrGdZxd7k3VSph1fvvYh64bMi6bARaRGme9+ZkTXb2IpC1FwPL0FRZasRLvUIiSevqpwOYPCnevka17WS1UsXp+GvR9XkXp8URRc4aUXrdL93I64rkBdhOvt3IxrGaqjmHU36uQ3TCPg/uwLCVX11+GyPED8MSVXY2E3ZkISgiMojritgvQater2bvNe6ut6gWavNoE6AyhOs8v5OjTKvIsz7Zy+O9KR6qqq6+8N5VEoHPZclrYew7xQVdbaVlu6R02hG9sgjrG71/qY2CmwVeH+UuO6kafU4KuUb1EPx9yj8utFq6uAnaK3fSsbaSm6EefXu6iXKm6vRXGCLoH4K91iKJxiCbHjVfKekvZvCnoEWhpQZyuczt4cw8y+I5ObnP+ctq656F9y7qVA2iSK8LFHEbSjr6Z0XrTFuZPxLadbaC+4NVNC151Gluok70veHD4njFRXSVtVkLzhob9EnfhAv7c6f6+U2iaBeYJFe0jTnHNTjQpgFjdBoRvX2askn6ryXoBFlnI5Emt54SV3bqHNNwkYSAzwjAxfvSN+rovTHqYLg3zWR39Kqm+JT9EYIap/VHEWRvzekU78gt2JvWdUKJjddvQ1bmO0JyoPVBdUz6Y5E3J58ks4S/u9bjUqyKM/xzozsntjqxTXExxGKQYImyM5R93dWa5NoCvyRKyFcj9qoxuCsqWQfqWdE5t+vh3B9ffixpMhK1VOLoPOT5bNe6bedz9zuEaXXGO6ONL2jB6/HUpRNqx6bS0Wsd1PjE+UjfQLXsJR7ABjFGdzHOZnVy+vPDj0N1ZM4nVncnll98qZSg14EitTjdolyBPh6n00AvNU1k7s5uftctd8jiV54mCoxaAJg2ByZIDWS9zqXIDVi0HV1220UtLP2y9OmKYREVsa1sAf6QS7s3l/Oi6nXw0tpNBPsBbaGlzKsTb74G62jmMZSHgWSnTNRr5ohDeNwXBftsLRhje9EJrGgy8z0WsgNddUxCQEepkrcmZHd74BfWAQJEG/5XmHiTxs0a9tfXt12m4Krmwrl3VRkKk3Xdxng8bN2ydPoW8ReayNgVYYlop7fap7tet+DSr9J3Hr68/FGUvAvkOTH7ynlXuccr827aaTo45fGS7vLv/dy76YiE9XT8/tM+xvmvBqZoguIpEeAfu8mL9WWk8V9C1Jd5PmsBNUjKFRFULo4pO3cEWYzCHofvaMM77NRjaE+ybkSRR5J9Mm7AlEUaVJUVF1cwdBKC/txPGAaQLexcs8lUQd/Pu5nb3kuq2jjJG6qu9y0iNtw+7+X9357cV/SAbQwgBZ2dEZ/RXfjdA2meQp0bx289djAEs6WFzhbXuAwppZdU+mdqIW03nn/8xL0XV1tQdj76k7AC8q7LhTjAhtny4HCqZv8w8Ki9ITj1qWeIXVS3zcN9Va93yvo2iL9vmE0Qh2TpN6Z00XBO0r17geNNsNGBEGdkIlMAozt0ess4lxXkyrowD1F/3hRvLS7ftWqm8oeyGqGfmkT9eJ4w4e0MIQBtNDiOR+3x1zJ6yNuQ+3vHSWh0knSD7/aPGttqJKwKxTl+cuLnRnJfhwEwEJnXbGwZ7MeW0PSbsvuiACC3784MZiCvJvKVW8JLjpUYO+mwo0k8iapXngjzFDOm0btpVryJ+rZCVL/+GdVhwkJ/wjDzcudaxG0XIAnnEdtI4ndRf/wlXhpd/tmAUcSIjID+CSwWlXHO8fuBvZzkgwG1qnqQSKyF7AAWOicm62q05KudJokpe5ptsYvjZX+0rpHVkAnR1G9vKLq4lcXeT2WWmkpG2X4RwtRrtJBoTwgIffqhOwNIjIFuBroC9yoqpf7zg8Cbgf2xLT//62qN0flGUfddAtwLXCbe0BVT/UUeiWw3pN+saoeFCPfhqMWdU+zkMY8lGYUprVQhPsQ9mxXO1ku7e9RSwfA2+sPsouFzYXw5xHmIZUESXg3iUhf4OfA0cAy4DkReUhV53uSnQ/MV9XjRGQYsFBEfqOqHWH5VhQSqjrLGSEEVUqAzwBHxf8qjUs1OvRqr0mDIjQ+edBo37kI9Q2zCVTzDFUrUGq5plaXW68XoDsfwsVr1PbbQcOM2FAaTbgjk5p/R9e7qX4OBRap6usAInIXcDzgFRIKtDht9w7AWiDSIlKv4fpIYJWqvuY5NkpEXgQ2AN9R1aeCLhSR84Dz6iy/kCRhwMuzHlmRdP3c+5eluqnSb9Zoqq96RhL1lJPENXG86NznI2hirFdw+K8LOpfk0rhVmIaHisjzns/Xq+r1zv4I4C3PuWXAYb7rrwUeAlYALcCpqtHjmHqFxGeBOz2fVwJ7quq7InIw8KCIjFPVDf4LnS92PRTLcF0LSTXuWTdwRaXW+5nHPavVQysNj54k8Rpwd3T2l/Ni6qMF//X1jBziXuu3twXNzo46545Oav2eWp1305oIw7UEZe/7fAzwEkb7sw/wuIg8FdRGu9Q8mU5E+gGfBu7uro1qu6q+6+y/ACwG9q21jGpISj8Yd+Kbd0KXt5dRTz2CZgrXQ9jku6JMUgyrR1BPsBrcSVKuT3uehE2Uq+RqmTVjGV82WWwRS1jEEtpYywfO5j7ncZ8h//eu9tmr5jkIOhdWnvc3cUcE/gjPUft+krDXaVe8rQLLgD08n3fHjBi8nAXcr4ZFwBvAh6IyrWck8THg76q6zD3gGELWqupWEdkbGAO8XkcZsUmqJ1aL/tU/AvDGeaomz6SNw0Hl1tOrq+f6INLQdXvpzw6FHpkF1S2pelbbCw9bJriFId3rMnhVNrUQ5TUUJ884Hk1B+fkbeH8cNvd8XNWhP4+6Sc4m8RwwRkRGAcuB04DP+dIsBT4KPCUirRgv1cg2Oo4L7J3AZIwubBkwXVVvcipwpy/5JOC/RKQTEyF9mqoWOyZCQngfrigviTCKqHLwUrS6hbkweo89w8zM6uES1vtNUmDHIam8vUbeJOvrvS9+T6F6ywkyNnvLizORzn+9W0/vaD9J4Z6Ed5OqdorIBcBjGBfYGao6T0SmOeevA74P3CIir2LUUxer6pqofO1kuhrJo1Evcq+4CPgN11CskVNUWUX9TaOiH3sp2ncIs/l4505AdJC+oDXSKzwXNU10+/Bw0d+fES/tyB8XcDKdpYS/9+NSq3qpUhl+qs23aC9uWvh7jFFui3HJ0kBbxN/IbSAH0MKqhJZPTVvw+vP3q4ChNMoPUxWFzZvwq5H9HZF6XWB1a22XZkHTCgnXaPmzCbOYlNA6QP6HwDX0DeWg7rg2N+yxLZ96q/aXoNoH7Y/7wccXBp8rYuOTBmHfM6vv34z3uVr7WLUePm5EWVcl6DbiB3Nl2ep3fsIEzUQmscGp8wgm9Ahv7v0+lZ6XMDVS1D0pgropLRpaSPzRCQwyYCB0bYVOx43M3HDzkGx+Hx4bA8e8FphFRa4YaH78b7SXeg5jGc91E+Zyy0u3dh+bpWaJ0GXv/hU4oLbCauDt1VcD4YFfbhl8A2euOzez+jQqD47Yn2EjFgDwkWfhGccxqu8AOORPOVaMUgN666Gmkeo7wBw//Ol86hMU7mIALYxgJKsC1p5we9luGP1V3M44+T4Az+hMWmlhP8wzukefqZFG3L040tl7qkz9c2Sf+3m8azIAgxkLIWtgBPX4/W6tWQt9tSvTxahEA9okLMXnremw3plrutOBMOI7+danyLw13fzvuy3sdkl4urGMp8NpRL1GbXffH5Lbu++OvN9lSfd1OzKSOczqHvnPYVYPIeQtwxUSb/IUOzOy2/PKTQ8lNVGUU0MYUbYMv9oqgJrsBQe0ij702Xhp977a2iQslsT41+951hS4ZySJhXZOkUmczixuz7TMznXQd5DZl6DpWB46fN5B3kbXO18oSGXjdbMtnTd5edVD/mvLG/anuvfeZUms8r2fK4UdqbSf1gJWRR5JFGplOu/kpy3v5FiRDNnyjtlWmrXmzdBT4d5dj2HjbNg4O9/6NTLexiyNAIVpkKWAcJ+1Ptsa4VBJQED8ZXzjUqtzgdd11s0jziqEUZMac6PgK9MVaiTh7U30G2oeoDdeND/idhPyqlW69B9W/rn0oj4Gh2ddm3DW3mP+bz8RBu6Tb11qoasd+gzMuxbFwatljrovKy6PVj+FEaX7h3iNc5SHn5t/0EgiTe+ptPK23k01YBrLtqYVDo3GkFPyrkH1eF/qTc/ChsdhzkXm3ORXcqyYj1O5l/f5OwAPLP0OA/ZMtzxXQESNHNw0y74bL0/v6owQHIcsqoH1u6h60wfFSwoyMkfNeYgirr0itVGHNVzHqEQvNlxfPuADAM66chtaL8y5Mpam4of9jbvft7aYvqD3VY8SEF3t0LXZ7PcbnFLlYpJG7z3FEUFthuthovefGC/tvjdYw3U362fCoKkwY0cTP/DsDadWuKIx+dqibQACe49znd77+HuMcTHvF9bSWLjCAeILCDddv4Ko5tJozAthh/BT4JFEYYXEIDPXhjPfMcLh7II8tEkTpVoYf09pvygComMp9PfEmYxj7CwC62fC4v+BCY+az41S7ySoRr1UKZ0leYo+T6KwQsLFGhuLRdr68rRwOx0db+Vbj0p0LDX3+PmPmc/1TORzG37tiH6PiiAgen1csgIr3AvlApsH3+mbnQg/lXsBuHLge1w58D1u3PFho/9thzlTot0Bz5YXsqpmUzNgz2ILOrduh/ypPgHRua60HyYgXBdYiO8CmxY950P0ItREjIiz5UHhRxIdS83/AXvCgi+Y/f1/XX0+Xe3m/5pbYJcvlY5Pf7sPlw0LvCRx3Jg0F7XvZA60wx0DpwHwBNdBxEsycdsWZmwy+xtnw/Y5u8e6jVBR1GAWw2LH+eH5T0aH7SjC6MHiwaqbasfb66tFOLjMGPYwAP+24biy4/55ClljhENlLthUWuAvbwEBsGV53jWojacONv+PfAFedTxKDnggv/qE4RXCkzgdqDzR7v2noP9uZn+fiLksbofJCofiUAAn01AKLySSwi8cLPWx7bi8a1AbR3q0dkUUDi4Lzirtx52F3XdI9ETHsNHDD/t3lnlC5UGvtkkktzJdKlS0SYjIDBFZLSJzPce+KyLLReQlZ5vqOXepiCwSkYUickxaFU+CsYznlsE35F0Ni6UHBzwQT4hdu90/uHa7fwDRgnv2EbB1vdn8I4h6BERSa6bXa5NIco37PHDtQ5W2PIhjuL4FmBJw/KeqepCzzQQQkbGYZU3HOdf8QkT6JlXZenANxF6evmcu96xblU+FLIXCjVBaRDbPM1tXe8+YZl9+b1++/N6+wRc6tC+GYROM6ippG1JRDM5J1SGP76KYsBxxtjyo2IVQ1VkislfM/I4H7lLVduANEVkEHAr8rfYqJkOQh4cJNWHjR1vKI5QWgXduMv9/8SVlm7Gm6x9kQ4hybd08z/zfsgL2uSbhClqSo+DzJOpxgb1ARF5x1FGOuw4jAK8n+jLnWA9E5DwReV5Enq+jDhZL4iw8C45iGkcxLbc6DDvHbBNHbNd9TLU6dci248y2w5GV0+ZNmNoqC/VPEiqzuuupMbccqFVI/BLYBzgIWAlc6RwP8pcI/Gqqer2qHpJ1HBIvnevKf9yfbfM279xkenFPfjivWlnyZr+bjddZXM+zNDluyeayOQxva7Q6xFWrunaKRiFMbZWF+ieOyqySEKh7+dIGt0n0QFVXqepWVe0CbsColMCMHDxBG9gdWFFfFdPj0MGTyn7cr36wa3cP7oG5a8rSulEqLc1PKy1sehE2JbQ2epJUclvtM9BsXpfp3kjSI5BUhZWjboqz5UFNbg0iMlxVVzofT4Ruhe5DwB0ichWwGzAGeLbuWqbEd3bbzKd9Isx9uPbfbhG8D8svM8el3xImXFJ68G4c2clxSzZnVVVLhqyijbf+J+9a9CSveQ2zjzD/81pTuxaKYEyvikaeJyEidwKTgaEisgyYDkwWkYMwX+1N4EsAqjpPRH4LzAc6gfNVs7fJ37jjw2XzIlRho/OAe/WzJ7zxHPgMf+7DNf39owHocBa6GnUVeGdEH9cYC51ZYrLlnfKJlfvdnF9dqsFVQewqJvy16/2U5CTRrR3J5WUJQIu96JBdT8JiwejxNzvqpe0OKx0v8qxk/6srUnLzTjIw5l8cZfJHCqsTKAw1rfUwdpDoHR+Jl3bCo3Y9idRw18/2LpFq6d24oS/6DjL/G0U4uHjr6Bo2P5jfM90kTq9r7ewznnXtcXb4nBYF6KuH0muEhBUO2ZGGyiMN3IllquUNbpTB2pv2JG5i/75n8f1O4/+Rp2D5/V7bMmf5JqZ39qxEPQICypcFzZLRjMyt7Mxp0nkSFksg/YcVV0D85VCzRbkVRq2r7hUE93EOl23tk3uYbTCusl/+VXklXj3RbN6w4VBaq6Lo9BoBEdP9Na/RhrVJWBKjEZZYddWOL+gsHhjxT4BxYKhHhx8UNj1oNbg4K8RZGprabBI7iv76sMrpAA75k7VJWBqYoguIzfNgm7ElteOnVzxnduo08gZ97yBB8PiHzP9WWnJ10XRXv/Mfg9LxuOHJXfL+TvVQhAi0RfZuskLCUliuGGhe2m+0JzMxapuxcNMgs67IF17MNnR8Ky0c/XfzfX406EbO3nBqpuV7CVqZ78Q9nQlBTiyzau0YtTawRRAueZefZ8iNOFghYSksSQkHF9O7N8Lh3yLWXUgDb0N01vpTObtgKqcd+FDstGE9b3+DX00P3Y1oEGWHqEagFEH4xEUptneTNVxbCsGKy82yrBtnl44lbaxzRyYAx3BpchnHoOjLhd7Nyd3L61bCG+vIGxzP3yhHxUTyB9VbxJKKhupqGv2sBUS9YUCs4bpSJazh2hJAEXTF9VJ04ZAkafTeG2lEQI2G6/1bRG+N8KjzcthT2Ruu7UjCUliKsqBNXA5jatnnuAJCFWbseDczdrw7t5XRkiCJ38r//eOORNK+b1HhxOsO/qnQtTXeVgkRmeKsCrpIRC4JSTPZWVF0noj8uWKediRh6Y28NR32+F56+cdxdw169R7/EBz99+YfdTQptY0kdhC9OebSBP/8t/AynFVA/wEcjYnI/RzwWVWd70kzGPgrMEVVl4rILqq6OqpMO5Kw5IJ/Kdnbh3yPR0bBI6OyKf+FG/dPJV+v7jiugBAxAfp2lRY+vrBxBUQjj4LyJiGbxKHAIlV9XVU7gLswq4V6+Rxwv6ouNeVGCwiw3k2WnPBPXjt97XRYm135JyxfkHie1aiX/GkaSa3mEmQvyMqO1Az2qjLi61KG+lbzvF5Vr3f2g1YG9U/T2xfoLyJPAi3A1ap6W1SBVkhYLHVSrXCIk64RRhNB3kxZl10UYVGPgb1Kz6U1ESqtOCuD9gMOBj4KbAv8TURmq2roUoZW3WRJjGrMW24QwLxYcXn9eVy/w5OJCwhLdRTFuaEgy5fGWRl0GfAHVd2oqmuAWcCBUZnakYQlMapp/DoWp1ePOOwW6PcRj8fGmP8tO08GkhcOVoiUaDAX2JpJaGnS54AxIjIKWA6chrFBePkdcK2I9AMGYNRRP43KtOJIQkRmiMhqEZnrOXaFiPxdRF4RkQccizkispeIbHbcq14SkfxXkrcUku0Pz7sGtXPMa2b7l9m1CQi/q2yjkoWhOgsB4bq3VvN9kv7uSYwkVLUTuAB4DFgA/NZZLXSaiExz0iwA/gC8glla+kZVnRuWJ8RwgRWRScD7wG2qOt459nHgCVXtFJEfO4VfLCJ7Ab9308XFusBaGoW49oJ6Ir56o+kWzT7h7dmn3ctPKv8MRyM1ucB+aDvRG2JGRZn0YgEn06nqLHx+J6r6R0dqAczG6L4slqZnw6PR5+O6wEbhjSr7n/2KtRpNpca22h55pbKSyKsR1FVFDsuRhOH6bMD76owSkRdF5M8icmTYRSJynog873PnsjQwqjD3lNIDPZqRPDhif978Jrz5zfK0/nkSjcKgCpqiv/2z2W5oeTKR8hZsvTmRfNIgqPEtiiG5kXAD/BVVSNRluBaRbwOdwG+cQyuBPVX1XRE5GHhQRMap6gb/tY5v7/VOPlbd1ASImJ7k204j8Zr6Arb9pLRbzyI/RabtXfP/vPcnJ5LffZyTSD5FJ8iVtTdNzkvIcJ0KNY8kROQM4JPA59UxbKhqu6q+6+y/ACzGTN6w9BK8L/n/jqV7ac8i6dUr8eqJpf1juSw8IbB+ZnkvzzVqJ0UBouZkQtAIpIUhqcSDKhwxRxENNZIQkSnAxcD/UdVNnuPDgLWqulVE9gbGAK8nUlNLw3Ca3AvAne0n173qWx4c8EBp/xFnEZ4wKqmf6qWRhGvStCU0Bb8h1qAocGegopAQkTuByZjp4MuA6cClmNf/cTFP8WxVnQZMAv5LRDqBrcA0Vc0w2IKlCDzJWQC8cmKFhBZLgQgSEFnN6C7yiLGikFDVzwYcvikk7X3AffVWytLYPHWueaFGX0tDjSTieCUdy2UVRxfNxvMfM/8P+VMy+QWtYBdmBM8Dr2CIClWeZP2KLCRsWA5L4ux7g9kazTgdx3aSp4BwG5KFZ5lt04ulc0cxLbVyt3aYzaVeN9dqG9ck3WrjEMdDK3EBpjG3HLBCwpI4eRrZaqGaRWNezViFtunFnvdzv5vNNmpCqeH85VnpBTc4bJbZXJJ2c43Ky+2xN7NbrSp0xdzyoGmEhNcDYO4pedem9+IN3Pfyp/KrRzVEra3sF3Zeo3YWbOdZ1nLTM+XnvA3nfilOp/jLoWbLg2YWDmUUeCTRNAH+REov9LjfwjFiFrp/jB8BMInTAZjF7bnUr5lxBUO/oeXHT3poJFRY3L7oFMm7aNsJ5l73HxaeZiwmIs58IsPxVMU7y90FmpJfg8NiKPLIu2mEBPhf6B+VnbPCIT3cRsuNMzT7CPM5qodeZEYzMte6d7WDDOh5PI6NJ0nh4JLGAk2WcoosJJpG3VSJOVPMZkmGznVmC5roc/jTZmtU8hZuD476p+59Edi63my1krUdxVI9RZ5M12uExMQ/mM2SDB8aPJIPDTYG34NlEgfLJLSjwkUeXLWIS6PGcgrDv6jRgyP2563p8Nb00rFjuJRjMGrR9sVmU4UTlz9Xdm2/weVB/6olazuKpbmoGCo8k0rY2E0Nhf+RKZLevlG5cceHAThn/XEcLJO45lDjTvSRZ+vP+6mD4cgXSp87lkL/Pcp/tyc/bP5PfqX+8noxNYXxHjNQ9OoR8dIe+0b2ocKbyiZhyY5dxbhfLnsnGe+TX2z3Kl/edED357GM55V1Rr/u7UVvngcHjMvXZpAUVww09+4b7S2cs/44AD6YD3OYlYhwcPEKCIABe/ZUvU680tk5OrlyLVVQ4G6yFRKWWPhHD65rYpSnTTV4BQQYA2yQimXbcdDoHlMu32gvzXPIejS2z39g1idz6LON+d/Vbgzk7tyRP1+2hBEZzB+MCn+RVWiMPFH6xky5NdV6BGGFhCUW48TYEDpoY97ibBrpiUwCTM/apXMdbPxreWC9+3czht5Pr3iOpd+GPX+QSfUaGn9gwh18K7+4I7X3Xo7Op5oG3JvW3W9hCItYQgtDQvMJO+edhd3YAkRQDXBnC2RzqjUJwgoJSyy8rpUD98mmzFGcAZQLiX6D4Y2vlaf79IqSodcKiNrwN/buSGLQ2CWB6dy0e2GkywjeL/udKjECM0twgyOM9uQTACziulDB450Z38ZaWhhSFinWH7ojScGRboTYPqBxY9hkLySs4dpSWNxGoY21ZS9oVztMHmgnR+ZFWCPuFyBRPX3Xu62DtlD7Uist3SMN93PU+tphz0vcgIK14smvJqPy6AED9KpdW2OlPf6tZdZwbbG4hDUeMgAO72sCEc/aGi0kln7b/LcjjGSIaqj9DW9UQxw16S8sCqtX8ATl7Y4qqqlHEtSffzXqpuyxQsISyFjGpzJ7Nykue8MMz3+yZ3Q6KxySJyoia7VGZlfQhAmGSoLIS5Rdw0vx7BfVqJuyp+JkOhGZISKrRWSu59gQEXlcRF5z/u/kOXepiCwSkYUickxaFbekx6P7wNz2uZzANZzANcyZApf17+Ky/mYh3l9s92rmdTqMqRzGVJ462HwesKfZKnH/bv/Ubdi2JIPbkIetAVFNI+ymjcovLotY0qCu0YIyINaWS+0q2SREZBLwPnCbqo53jv0Es0zp5SJyCbCTql4sImOBO4FDgd2APwH7qmqk35a1SeSHuyaBN9qoKmhHeKygPGIbtS82/585GSa9GJ3WYsmJGm0S2+sVw8bFSvvpFc9lbpOoOJJQ1VnQY7HZ44Fbnf1bgRM8x+9S1XZVfQNYhBEYloIyakJL2boEYHz2o4LJpS0grhz4XtlnVei/u9miBETHUvhm33YmcXp31F9LetSzEFCtCwn5r8t6QaKoutSOsUnE2fKgVptEq6quBFDVlSKyi3N8BDDbk26Zc8xSUIqnn4WL2o32cvM88zlqVAPlE/1+3DmwO4bUkkthiLP05qCp+Ud3bTbyeHayNkpHkVTZqn1R3SGRvNIg6QB/QfNGA1VJInKeiDwvIs9XW4iN5pouq66JPu9tlGfseHcqdVCFbcaaLUpAzJoATx9iNnf50T4DzTbqKtj+X8wG8Oo8KyCKQrOvNlcdzTmSWCUiw51RxHBgtXN8GbCHJ93uwIqgDFT1euB6qN4mYaO5pkvrhdHnvSEkzt5waiJljmYk85eaRtwffM6Pd72FoHReTxlvaI9t46l9LTUQZ+6B93dx5zTYkR24huuiUutI4iFwpsOa/7/zHD9NRAaKyChgDJBgqDJL0UhK9/+alhqLKAGhCu2LSqOGIFoYQgtDOImbmMik7vAenesSqWoitNJStkbAFQPbugP+NSJxRgXe0UMcT6RK+v68bBFB1FeXPqAD4m05UHEkISJ3ApOBoSKyDJgOXA78VkTOAZYCpwCo6jwR+S0wH+gEzq/k2VQLGx6Hlo+VPu8qaU6ZbwxO5V4A7tKT6Vxjjv1s97ayIHJ+wnp/X++zif/u2i5WuV/dbSHvr+gZYymqTCg1Kl7VVZRLqzddpRGB2/gs4pyy4x8ePB4KMvdjFW1lQi7qd+qtNNJ8h/rqUuyRhA3L0YS4Bt+iq1du2+lqvrD2K0Dl0YOLm871gHKN3Jb8STe+UUNQk3vqPv2H6Q93ird84Gnv3GDDcljqp+jCAcwiO6e/fVykcOhYav5veavnKKPRhMMjo+DYN8qPbXT8ALc/PPv6pEHSAqK3CB2lD6rb512NUBpSSPSG+PKNwupfQcsRZr/tadjlS9HpS6OCaAHhHT2MkZEksYaEd2Ejd92ErHhv/fcwmtoS65/MrvxGpNe832pjNyWO+/A8PHJbAI5bkn343CLz4Ij9ARi884LUl6OsJBS8qMLW9WY/as1mV0CUhEgyHjDehY2yFBAAp6+d3uPYbpdkW4dGo7eMJIxNon/elQgl6XkSmXLcks1WQARwwvIFnLA8fQERF68XT7/B4QKiq91svx5ytV03u8Hwe/e4s5Hd4/5Z0i7eNSKqpZ5rvYTNnK7GY6le76ZmnCdhsfRAtacBOsjoHHbtTYMeBuDfNhyXQu3KSSLKbZDKyrUzDNinfGnXrNVbWVNpJnRYJNcoN9hKo4ik5liElVNLoMLaEEhoJCEiU4Crgb7Ajap6eUi6f8JExzhVVe+NyrOhhcSWd8x/92X0f7Zkwz/OLe0HCYoo4eCdu+AVDmmrGpIIgx7U6IcZoZtZQFjqQ+lDVwKGaxHpC/wcOBozsfk5EXlIVecHpPsx8FicfBtaSBw4bLyzZ154KxwM/y6vA3DMbsdywvIFqZe37w3Ozg09z1Xr2urSO3TRliKTnU1EUE1kJHEosEhVXwcQkbswQVfn+9JdCNwHxIqh39A2ifnMLfTCOHnxS92bX+remQgIP4cxNVY6r3G60ewPro3lzW/mXZPGpEgzpaOIE2YksSiw9I+1YSY1P+/ZzvNkNAJ4y/O5R4BVERkBnAhcF7d2DT2ScJnIJOYwy8aDKQCzdWbNo4dG46Xf7A9kL4gbnWYZJUatt10dAvFHEmsiJtPFCbD6M+BiVd0qMV/AphASbjgIKxzyp9JzV1TBMJqR3Pcp8/wc+FDpeFc7vHIifPgB8/mJg0rnvr58U3YVtBSWRNa4TsZwHSfA6iHAXY6AGApMFZFOVX0wLNOmEBIWSz3MPgIOe6rUwVBKwmzRBXDMzBYmD7wZMLGxXF7TJakLvS3vWFtbs6MqqCbSFD8HjHGCqy4HTgM+V16WjnL3ReQW4PdRAgKaSEiMZmR3JNGi9lZ7I6pwuBg7xTPMDJgolz+HP20EgxsY0dsoG6N8G3djhMPdUtKp/+8p6atNrIDoDfShi3gBNaNQ1U4RuQDjtdQXmOEEXZ3mnI9th/DSNELCqpqKiREGM3sc97vKrv5VdbO368W7JoVL3AbZVS+M+y3BWmCCXYHj0rkueka6pdkQSGYkgarOxPfChQkHVT0zTp4N5d205R2zud4l808zm5/DmMrCs2DhWaVjG2eXJjpZ8uO188zmb0Dd+E8u7hoQaXHIwFL+IqVRRDVECYF6RkpWQPQ+lH6xtjxoqJGE29NTNUHf/KOH0os5k/1uLr+2WSJtNjr7BsylABMc0Mt3dtvMz1aYBY1mcXv38Rk73l3zanhb3oF+Q91P5WtfWLWOJT8SmyeRCoUXEq5bq3flMqhevdQoayz0Vvyqpk+veA5jhyunnuVSdx/WwttqVEVjZCR/X9c7VJSb50H/EeUjlHrUYZZkUSS3UUIcaq6ZiOwH3O05tDfw/4DBwLmAEySDbzl6sppoY62nTPP/rzWMCqxwqMz0fsaqPH2jlIWR2DyvOPfP37iNZiRf2+73QHmUV5f2xeb/gL0BX8yg3qLWKcpvZwmj2EIikZXpnFggy4HDgLOA91X1v6u4PrAS/qrZno+lWrwjUfv8lOhYGrxcbPtiGLhP9vVpEmpaNW5k37F66bZ3xEr77xsnNOzKdB8FFqvqkriz+MKY7Rgw95hi/psFZ6z3kiWcy/p3AfCdLT39MAY0SAiIrAlbT3zjnOryibMA2GjK7YfVzFD2pg267iimAfBE/CgTBSSxeRKpkNRIYgYwR1WvFZHvAmcCG4DngYtU9b2Aa84D3LgjBwfl61UtbJxtjc9ekl6dr5UWVrabvJb9F+z5g/C0RVI/AcxxOhQT/5BvPSy9nhpHEuP14m3uiZX2/E1jMx9J1C0kRGQAZur3OFVdJSKtwBrM/KTvA8NV9ewKeQRWwjsk9vdGejtJ96AaOabSCVwDwINcmHNNLL0RzwinpgZ8z74H6MXb3B8r7QWb9m1IddMnMKOIVQDufwARuQH4fa0Ze4fEVkCUU6tw8I9AGlk4gJkUhyMcxg+sfyEhS+8iSIWV9Ci9IgqqfbMpqwaSmEz3WeBO94OIDPecOxGSeWtvH/K97n07Ka6cV0+MH7Z4FW1lAkI7zBYlIJZ+u/yzaZjT5et9yoPnbXrRbKuuKU8nA+CxsWazAsJSLUGCwPuO1JpHdQhdDIi15UFdIwkR2Q6zCpLXy/0nInIQRt30pu9cVRzFtO4e82mvT+f0weZ4kW0TrlpsLGZBpCQbLm8Pxx0BTJFL6dr6o1gPqhtSHUojiKgV09w02gF4bBTf214JjUeREP/dVR7LZrsJwekacfRjsZQjUOCRRCKG67orEWKTyG5lqN7D+bIMgGu7do+97sOuYn8HSzGIowrKut2o2ybR50C9aOAfY6X96ge7NqRNIjVsw5QM7lyBffgchw2+D4jugbvqpFKa/H+Hc8VMmb9BjVuV21i88WJb6CjD0jtptHbDzLgu7kii0ELCEo83vgYjvlr6vO5R83+XLxkPsQ8WGqN/y8d+1DCrxvl7g8P67Gd2tpp/7jkrIHoXRRQASdRJCxxr1QqJBmUSJvjdCp7i6Ktg0VXB3l9hk6b8FElAQM8X74db7aNqaVYEVSskLAnjjYxaD7ftdDUAF/Q5iZ/r7t3HrT3IUg1ZPS+Zu6dmgtBV4Ka4uDVrUPJoXI9iGp/Z4UwApr9/dOzyVWGKvA3AY+zevYbDHGY12UtoSZusnpdmfC4V7EiiN5HHQ/wE1/Evm39R9XVtf4LFlAKLzfGtsWCxWOJRX+dQCm2TKG7NquD2Id+LNZGsWWmlhfOv6MP5V/Tpjr8Uh232MzPZ7Wx2SxS3DL6BWwaHrBblEOf9G83Ibk+7k7iJk7ip+9qxjGcs47ttbWG46fzlfU4eqaouQcSdkJoGxsOp8pYHhZ4nEUUrLfx0538HYED/lzn57ccSr5clWR4csT8nLF9QduzJD5v/k18pT3vFwDa+0d57BX9RCesx+4/H7VnX2gN3r0vSRpGAqrimOQy79zlYL+gfL4zEpR0D7DyJuKyijcveNWsZ2XAMtbH6Vz1XhEuTx1Y8glmbqsS6d/d39krCY/M86LdLC9/wLCl6/Q5PAnDe+5NTraMlmrBG1H887mS3Whvleq+PyjMuyQmoYns3NexIQhUWf9nsr3oBPvJs+fmFZ5n//rWuextd7TB8YKlH7n2gwxaeKQJ2eU1LEAX3bqqplz+izyF6fr9nKycEvr2lrx1JxOEopvH3L17H/r8OT7PhzcyqU0bRXEdX/RzO7bsegN23fZZp71cOfFWE73Ca3AucnGsdLOnh2ib89rCx1B/JtwjPb7XkZW+IQ8MIifbF7jrF8MH863jj+9HpD3YXoIkIYNcbOOhrLaxy/BOOen9a9/HRjGT5z8pf0Jc/Zf736dvGAQ9kVsVA7rYCoqkJc5Z4t4ITRZzGv9EEBBg32KLSMEJi4D6lWcHVrIrW1V4e6XThWemqoP4xs401j8Ix15ie0s6M4xlmpldgBbwvjH8NilFXlaf95kOXAvAYP0q9XpbGwZ1wuesuX+HjC5PJs+Bqo2xRULUjiZr4oszi1q5JNV0bFgL793e8B+xUe6UqMGgqnMq9vDDzZOdzMd1Ld+PIHr25P6gRDq4twF39bh6/sS9zL+aL733F7PRYhLh24hrAvTTz6pR2JFEjt+kkbnX2ReCOnb/pfPpJzXn+3w07cZFHgLQvhm+OWQPA1V1Da87XSzsr6evYilXhP/t1AbBR2/hp16BEyqiXoLAefkNxYy8ub/HSDGvEN7OA6Mq7EhEUzrvppanm/4GP+NNkWaP6OQzzRWZrSdVU9O9gVQCWONT6nHivqyWPSiOJAqw1UZPn0W5yiJ7b7/lYaf+rUxrLu0lE3sQsNrAV6FTVQ0RkCHA3sBdmZbrPqGrsgeoxM80P/TZtiBif+TT50/7wsQWV01WLa4dY/yAMOsEcU4UxMpL5S82DPmDP0kzRO/TY5CtRJVY4NC/XbvcP5mw2v+8MPRiAznXmXL/B1eVV79yGWvOoNJJoZKN2ETrrYdQ1knCExCGqusZz7CfAWlW9XEQuAXZS1Yuj8jl4vOhffgvbjPXnX3PVCkX7YvN/wN7QuQZ2H2Z77BZLE1LzSOKcvs/FSnvZ1j6ZjyTSmOZ3PHSbEm4FTqh0gXbC1rVmX8RsH8xPoWY5MXAfs61/EPoNhbe1jbe1jdGMZPWvzMxniyVr5p5iNj8b40WIsCSEAhrzLw/qNVwr8EfHpvArVb0eaFXVlQCqulJEdgm6UETOA84D2GM49B0Czx9VOl+Nm2ujMPhEeOcmGHq2+fyaLmFXcWdD21GFJVvG3xN8/NETjgFsLDQvUbaMJOwceQmAONQrJD6iqiscQfC4iPw97oWOQLkeYAcZpJPHHcHhfW4DjIdREjMvi8iwc6B9stkfsLcZVUDzqNYsjc9JKx8jbAJwVIA/qGw0rpSm0nX+fTevoP0kHTGi8qg/f2Vrgf2b6hISqrrC+b9aRB4ADgVWichwZxQxHFhdKZ+NbOAZZvKMxwW1GQWEy8B9zP/2xaV9iyVt4jaaUR2WsGtbGMIilkSWEdSgxyHM4B21n3bPP0lcdVNRqdkmISLbixhdiYhsD3wcmAs8BJzhJDsD+F29lWxWrICwZMmjU9p4dEpA431N+eeVV/VIEoh37QXX82gVbRUb4Goa6FrXdyiSEIiDqsbaKiEiU0RkoYgschyH/Oc/LyKvONtfReTASnnWY7huBZ4WkZeBZ4FHVPUPwOXA0SLyGnC089kSE1WztdJSZtT2LqpisVTLUUzjoN/BQZ4uW8dSs9180QdlaYd+IV6eWTTEaZRRPAGiKF2xtihEpC/wc+ATwFjgsyLi8xnlDeD/qOqHge/jqPyjqFndpKqvAz2kkKq+C3y01nx7M66rLBhbxRgZ6XxaUoh5FJbG5QmuY8awh51PxwHeMPHblKXtPwxLxlQSADE5FFjktM2IyF0Yb9NuX1FV/asn/Wxg90qZFneli4w4lsvyrkI3A/eBjtfNBsb76TU1et6OpfnWzdL4rG8/gvXtR+RdDYsPM47YGmsDhorI857tPE9WI4C3PJ+XOcfCOAd4tFL9ChG76cA94I/fhNYLsy/7Eb6TfaERuHYKr6vssnfa2H/YSGjS2DWWbLioPb3Alr2FIMN8Ky3dhvtaUSMA4rAmYjJdkMtBoCFDRP4VIyQq9hoKISRefisfAQFm4tCWlWZ/8In51CGIYefAuiFmf9AJZlThep383z7rCxMo0NI4nMq9gF2rI2mScIFNSN20DNjD83l3YIU/kYh8GLgR+IRjHoikEELCi2uo3eVL8K2+nfxwa7pV3P7w0vrJMDnVsqrFFVrugkuuc8N/9qvN48PSuymScAib7xBFEqHC63V/rSXEeRy6khESzwFjRGQUsBw4DficN4GI7AncD3xBVf8RJ9PCRYFNmsOYWrboz8bZ0M8xzLmqHfcWnCjX8CCVhzTug+ZdgtHtpQ2RifxS907uCzise6AUKBBKvuz+72exZIXXPTVuIw+1h/wu2vyGAGqKq7SLHKinSrwZ7tfq8MgyRGQq8DOgLzBDVX8gItMAVPU6EbkROImS7rqzUp2bXkg0E2748f99aiZbnXfl4KnJLcQylvFAc09ktCRHlJDwN+iu3h6MkIgSGH69f9xJgGHpMhQuNQuJz0i8jt7PdffGChXeiJzETTzNV4H6h4hZr7/gjhhOP/IaBjIcgMcuPJl9rom6Kj5WOFiqwTuajrOGw85O2lZaaGNt2bla8KufalEFFWN00kUXHTnXIZymd4Fd+u3yz/dxDjszsvuBrYcWhnT3jsL4d3m97nL8PMiF3M3J3M3JDP2EedDdLYpK50czsvult1hcfrHdq/xiu1d5aI9dejxDi1hSFo7DJShkRgdtdDj7+3E8+3E8YEaw7vN7Ejf1mLXtfm6lpewZdQUNmGfXzcMdEcfBFXSjPQIsa8zKdLFdYDOn6UcSp//gdPAt1bmHM5mo1p6zaweIo+ZZq3NqKiMug6ZC3AiylXpMzbo8pKU+vrzpAABaN/VsQKsZTZc/X7/rvq7FIxTe6F5loHIe3jK9gqrSe+0fPcQZjaSLVuMCmzlNLySC1nJ+jB/Vlef27Fk5kUORPEosllo4W14AYJH+lBU8FThKGM1I2lgb2cjGsXntySnMYRbQszGPE2E26Jj/Ov/nauwjaWGFRI6cwDVsYX3ZpDl3He2DKtiKln4bXnciGUx+pXT8Ca4Lvcb/siThtmexFIE1vBT6LFcSEFAuHMJ68l7vwkohyb1CxG8kjxNoMKh8P9mMLBSlM4NyaqPphUSQS2sl4eDywWrY/SNm/97Vx3Dy25Xd1PwPnBUQlqR5eOS2ABy3ZHMm5Z15kFkTe8aL5nPaves44caj0rYwpADG6PiYqXRWSBSKsJmn3jkHhzGVwYfM5KFvmKHEv204LjS/o5jGRkxwJf+chWJ4T1iaiX9bku1rO+1F1xBsRgJxGu8g0hAuQXml1TFL711WumhPId9k6JVCIsxO8Awzu0NyTz/zWDrfixYOLlHqJysgLI2OV01UTUPpFwpxFgHy2hYqrTAX5lkYNL+iXgGVbmfPjiQKS9CEHzck93Wfgx2PDr/23l3jqZ8slqR58XLzzE66JHt7V60rynnxT8Lzqof8kQy8+fivC3K9DSu/3gY+zc6eWptE/iy/DEZ4gr26K3HJgDZ2+VLpuPdhjRIQN+74MAedeBxP/sV89hq1LZY0eeNroE57kpWAuGKgeSe+0Z7MHAJ/gxvXjlermqv4KMqWvCsRSs1CQkT2AG4DdgW6gOtV9WoR+S5wLvCOk/RbqpprcKERIdHAZ158A3Bud2/k+QeWREaCXeGssddvp+NY9zzs7WqirJCwZMSomEuLJklSwsESRvOqmzqBi1R1jrPW9Qsi8rhz7qeq+t/1Vy8dPnyheehnf+1cuKrUG4kSEBseh5V3mP0v3myD6lnyoZUWvreDec0WbBrNieOHAr1jNBtkF8g6NE4aKF1s1Ww81WqhnuVLVwIrnf02EVlA9CpIhcF9oDa+VSGhhx2PLq0z/Qx2KVFLPqyijenvH929/52vOCfOya9OUNmwm4Thtx51U7GFSbHVTYnEbhKRvYAJwDPOoQtE5BURmSEigcthich57hJ8SdShGlTNNrbnZOwy1j1QimfkGrXv0GO7o7HG4QQSir5nsTh4J4r95vy3+c35b+dco8qNb1KNc5wYZWHlJxXQM3mULrbE2vKg7lDhIrID8GfgB6p6v4i0Amswcau+DwxX1bMr5JFZqPDOddB3kFtudNokej+qcEPLkwB85teTE1n97it91nB119D6M7JYLElSUxjvnWQv/Sj/L1ba+zinsUKFi0h/4D7gN6p6P4CqrvKcvwH4fV01TIh3birtVxIOLkn0fkxZkwH4vInuzbFcBlS3vvaf9oePLTD7P3ltKFfvU3fVLE2AX43iPlunDG7lzHXndqfb8g58b1ez+tllW5s++HODobmNEuJQj3eTADcBC1T1Ks/x4Y69AuBEyH+RAlVQJ1x7n4Hh6braYf7pZn/7PeDPM25g6o/Ni3bAl+ofVWx/uPlfjXBwcQUElFbUs/Q+utrhwVH/BMBlK7bl2XUmGF6/web8r+8xz9YRp5SHy+4/DHrBygANSrFtEvWMJD4CfAF4VUReco59C/isiByEUTe9CXwp6OKscLVpUcLBq3Ebf4/530oLr1xzbrdw8QqIU7mXuzm5LFDgMVwKwFs8zLss6U7fvtg26pbkMM/xc92fXeHgMuQUdy/3vpklJkoXW9mUdzVCaajlSzvXmf/+FyMI79eKUi+5ecbNNw4ncRP3Oe4mG2eXRhBJ0dUOP93xPQA++bmd2O/mZPO3NB+rf0XZxFFLKtRkLxgsI/QIvhwr7SN8J3ObREONP/8xzWyVcNVL2hEtILra4b37jHBISkAA3QICkhcQYHqTF7XvxEXt+QmIzfPMBkYQbpydTz3qYc4UswFM4vR8KxPBRCYxkUkAvP9UdNqudrN5uXHHhxl6Zjp1s9SPiQLbEWvLg4YaScRFNVo4xB1lgAnpAcGztmsxQFss9VDp2U6Dx8aY/8e8lm25DUhNvfxBMlz/mbNipX2MH9mRRBIkJSDACIewsB79GUR/jD/tJE5nEqf3WFO7mjkVteDvNVqam6wFhJdZE/Iru7lRlI5YWx40lJBwJ8EtvrC0r2rcW90gZGF4h+EiybxsD3Jh96JGs7idWdzOnj8oT5NF+I6udnjuX822MF6HJFEWnpVPuUniOh5Ywvnqi5PyrkIobuTYRqTo6qaGigJb1rD3mMgcPhuy2tFDIxHltZUVX7zZHS01bjyrP+iPmPcZs+96uBWRG3d8ONYaJ2ngrj1dRBp5BUhlK51szLsaoTSUkPAy/zTY/06zX6nhbzbBUDRm+4L8eu+3f43vLU5sYOO3XyzG/dbZKfDzkoeAkIbSNzQimtsoIQ4N+/OPvStabXTFwLaKKqjeQHrxZkq4v4MIjJHyYb+/h9d/WDEFhO1IhNOvn9miyOI5c8vxr15XazynJKmv/GKrmxpWSFTizF+0cOYvel8c/L8cWtq/d9djuHLI1zMt3y8UXNfNasnjpU/KVtVs9OlrtijqjUbgBtKMizccSRLB+/LE2iQyxLU9iMBQN6RgziGUs6b1YBj7rAnJ0LHlQE5fOz31MuefZv6PvavnOa8eeyzjefkdMxPYHU08PHJbAI5bUh5P/21tsw12QejIwIMurk2hqMKgvnp1sbXA6qammifhbay8AqPZmX8aSH+z/6HbSt853cXbDaMZyWtqXvB67rW/rnnMB7BY6qSmOQzbS4t+iImx0s5hVmNFgS0a3p7s6/+RXz2yZtR/wrHjzFT0J359XffxWgVEpVAifz0c/vlv7qd0vEra/lT+2VUvPDqljYl/SKVIiyUXtOCG66YaSXjxjyTcUca+1yUbgiMP7tj5m3x2zU/KjmXd6/be3yJ7LFkMC74A+/8671o0PTX18reT7XU0H4qV9lXmZD6SaFoh0ex4G+mkIs26wQ7v2vt7nL52erch8eXZS7pHFv7HxaqELJZuamrAt5VtdW9Gx0o7n7lW3WQJZvM82GZs8LmkQpGXRljG2O0aE7c/3DujdQm7SgstDOn+bCk+HUthwJ6lz13tcMuwuwE4e8Op3cc3vQjb2fAbmaJooQ3XVkhkwLf6dgLww6213+5tx5WPHr4ortdQMqESWmlhHJ8H4AmMXaNjqTnXfw/wCoNGdznsjXgFBLgz9U/tkW7UhBZI6bd1OxptrO0+Vuk58jo0jGU8833rZIzFePL5j8fJrzgoWyluEDarbioY3ofeDcXtH0HkoeKJioZraUzWOxPlB3liUKpC55ri2JcqNepH4ThscF1omgypSRU0QPppK4NjpV3Gu82jbhKRKcDVQF/gRlW9PK2yis7qX5n/cRZ98faIth1n/qvCrtLCY1Pz6wG9ZT2Kmo5BAQGKs+iA1NqbL+YooH7M4qXJqJsqtbvOstNXA1OBTcCZqjonuoKqiW9OBRcDewMDgJeBsRHp1W5m+0qfNfqVPmv0uX8tv01518tuvWebyKTc6+DdWmnRVlqqPldrWbWcc7bna2kv+yI6mAGxtqgyiNHuYoTDo5gIZYcDz1SqX1ojiUOBRar6OoCI3AUcD8xPqbyGYdU10Hph+PnZXV8E4PJrYePT5ljfIeHpLZakKVq016jRQ9IjixaGhOYZVZY/kGU1bEVZl8xIIk67ezxwmxqJMVtEBovIcFVdGZZpWkJiBPCW5/My4DBvAhE5DzjP+dhOcVZuHwqsSSvzKAEBpfUnHFVTqnWpEluXYGxdgmnIutTa0DvXjazpYngMU8c4bCMiz3s+X6+q1zv7FdvdkDQjgMyFRJBmU8s+mC92PYCIPJ+1MSYMW5dgbF2CsXUJxtYlPqo6JaGsKra7MdOUkVYU2GXAHp7PuwMrUirLYrFYLPHa3arb5rSExHPAGBEZJSIDgNOAh1Iqy2KxWCzx2t2HgC+K4XBgfZQ9AlJSN6lqp4hcgNG19QVmqOq8iEuujziXNbYuwdi6BGPrEoytS8aEtbsiMs05fx1mjeGpwCKMC2zF1ekLMZnOYrFYLMWkaVems1gsFkv9WCFhsVgsllByFxIiMkVEForIIhG5JOOy9xCR/09EFojIPBH5inP8uyKyXERecraAAAap1OdNEXnVKfN559gQEXlcRF5z/u+UQT3283z3l0Rkg4h8Nav7IiIzRGS1iMz1HAu9DyJyqfP8LBSRYzKoyxUi8ncReUVEHhCRwc7xvURks+f+JBZQKKQeob9HDvfkbk893hSRl5zjqd0TJ/+wdziX56UpSSMsR1rhO1Iofzgw0dlvAf4BjAW+C3w9h/vxJjDUd+wnwCXO/iXAj3P4jd7GTBTK5L5gQttOBOZWug/O7/UyMBAY5TxPfVOuy8eBfs7+jz112cubLoN7Evh75HFPfOevBP5f2vfEyT/sHc7leWnGLe+RRPc0clXtANxp5JmgqivVCW6lqm3AAszswyJxPHCrs38rcELG5X8UWKyqmS0coaqzwBNP2hB2H44H7lLVdlV9A+O1cWiadVHVP6pqp/NxNsbXPFVC7kkYmd8TFxER4DPAnUmVV6EuYe9wLs9LM5K3kAibIp45IrIXMAF4xjl0gaNOmJGFisdBgT+KyAtiwpYAtKrjx+z83yWjuricRvkLn8d9gfD7kPczdDYmYJrLKBF5UUT+LCJHZlB+0O+R5z05Elilqq95jmVyT3zvcFGfl4YjbyFR9RTxVCohsgNwH/BVVd0A/BLYBzgIE9Pkyoyq8hFVnQh8AjhfRJJZUahGxEzI+RRwj3Mor/sSRW7PkIh8G+gEfuMcWgnsqaoTgK8Bd4jIjilWIez3yPO9+izlnYpM7knAOxyaNOCYnQcQQd5CIvfwHSLSH/Nw/UZV7wdQ1VWqulVVu4AbyGg4qqornP+rgQeccleJyHCnrsOB1VnUxeETwBxVXeXUK5f74hB2H3J5hkTkDOCTwOfVUXY7Kox3nf0XMPrufdOqQ8Tvkdc96Qd8GrjbU8fU70nQO0zBnpdGJm8hkWv4Dkd/ehOwQFWv8hwf7kl2IhlEqBWR7UWkxd3HGEfnYu7HGU6yM4DfpV0XD2W9wjzui4ew+/AQcJqIDBSRUcAY4Nk0KyJmYZeLgU+p6ibP8WEi0tfZ39upy+sp1iPs98j8njh8DPi7qi7z1DHVexL2DlOg56Xhydtyjpki/g9MD+PbGZd9BGao+QrwkrNNBX4NvOocfwgYnkFd9sZ4XbwMzHPvBbAz8L/Aa87/IRndm+2Ad4FBnmOZ3BeMYFoJbMH0/M6Jug/At53nZyHwiQzqsgij13afmeuctCc5v93LwBzguJTrEfp7ZH1PnOO3ANN8aVO7J07+Ye9wLs9LM242LIfFYrFYQslb3WSxWCyWAmOFhMVisVhCsULCYrFYLKFYIWGxWCyWUKyQsFgsFksoVkhYLBaLJRQrJCwWi8USyv8PSP2fja4ksCsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "traj_contacts.residue_contacts.plot(cmap=\"gnuplot\");"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "5fbe6449",
   "metadata": {},
   "source": [
    "If you want to either force a sequential color map to use the only the upper half of the color space, or to force a diverging color map to use the full color space, use the `diverging_cmap` option in the color map. Setting `diverging_cmap` will also silence the warning for an unknown color map. This is particularly useful for user-defined custom color maps, which could be either diverging or sequential."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "51f4d4de",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAmwAAAD8CAYAAADKbPhPAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAACI80lEQVR4nO2deZwT9f3/X+8kuwsssKCAcqhQAQ8sbhUPFBEVFS+8K1aqEBRFpda236r1qiA/r1prlUOUCC0qHnjgUQ9URFAU0BVFVLBgRVBABGSB3U3y+f0x+Uw+mcwkc88k+3nymMcmM5/5fD6TDJ95530SYwwSiUQikUgkkvASCXoCEolEIpFIJJLCSIFNIpFIJBKJJORIgU0ikUgkEokk5EiBTSKRSCQSiSTkSIFNIpFIJBKJJORIgU0ikUgkEokk5EiBTSKRhA4iShDRBiL6zOA4EdE/iWgVES0jokP8nqNEIpEY4cUaJgU2iUQSRqYDGFLg+CkAemW20QAm+zAniUQiMct0uLyGSYFNIpGEDsbYfACbCzQ5E8C/mMIiAO2IqLM/s5NIJJLCeLGGxdycoF0qqYpVoaXp9pT5G1SNBgIQoQhSLO16v4D/10UGYxrtDwq35+n29Yn9keYYK3IcMD8Xfv8B0L0H3b6u7di6iTHWsei8OrRgaDT5f+LnpuUAdgl7pjLGplqYVlcA3wrv12b2rbfQh2u0bNuetenUNYih7cFvSIkkCAothB6w8evlptawaNt9GEvuNNUn27nR9zUsFAJbFVqiHwYEPQ1LxCIxJFNJ/WNR42O2xooqX5PYp94+M/24OS+/+jaiukU16nfV+zqmVapbVAMA6nfVq98ZAFRVVKlz558dP17omPoZs9zr9/KzmIeXvzHVsCkN6r+Hqabs9bW7GGP9HEzLiczrOm06dcX59zwd1PCWYYyByPhpWey4nfEAqH1q31vpx815ed1vMdJphkgk3JIzn6NYxpIxqPPmn52Z4+JnLF6715/DpHMONLWGseROVO1/gak+d338oO9rmDSJ2qSQcGJVcBEf5Eb9afvU21cM8cFvd35G5/strAEInbCm/Zxi0RgamhrQ0NSAWDSGqooqAIpAxvfVtK5RP7uqiiokU8k8YQ3Q/77NXL+d79sJFCFTmwusBbCX8L4bgHVudNwcKCacWBFezNSjJqKcPrXvzaIVDOzMT+/8IIQ1AKET1nKFLpYRsvhrvh/qvmQqLQjh2e9VK6wB+Z+x2Wv3t945WdgcY3kNkwJbCPBT2LEzlniOF5pDp23CRnWLalX4EgWt+l31eZ9l/c6s0FW/qx41rWtUYY8L2IU+g1g0VlC75q8gbU5Yc0lgmwPg4kyk1ZEAtjLGAjGHNnf8Fnasjqdt74Xm0O7xsJJOsxwBDFAELK3QHRX+L0cihGQqjWQqLQh7pL7WgzFWVLvm6/3lq7xmfQ0rvadhCLBjjpTkY+bzC+tnrGe+FPdzDZreOfy8ZCqJ9m12w8YtG9T9XIDjwhq/ftG0Wt2iGg1NDWo73nfgnxXBLWEMRPQEgEEAOhDRWgC3AqgAAMbYFACvADgVwCoAOwCMdGXgZoBdc6QkHzc1lX7C7wHRdKk9xrVo2v38PCKgKZlGZUVU3S8KcFozKBfM0mlFa8fbACzPXBos4V3DQiOw+fHAKTSGlfEDfzBKAAQnpGhNlRwuvGmFLb258tcbt2xATesabN2+1dBXUWsmFf3WuJk1DOZhIiASixZvaALG2IVFjjMAV7kymEt4/dBxy+8sPA9GSRCCSq6ZUnuMa9Jy2+XOM/u6siKKZCqNWDSSI9CJ18TP5X5wot8aPydUpmFyx/DoxRoWGoHNjwevn0ECzVGoK6XrdjJXI0GMBxZwM6h4jAtW/DX3YwOArdu35vStZwIVhTnuDxcGIS0X18ydJYnXD16j/r3QmIVL4+EPpXTNTuaqJ0xxuFDFTaFZsybA/eEZ4+cp58Sikby+9UyguSbSkAlpIiGdFhAigS1ItA9up4JHqQgtbmPnut0SnPw816gfMbBANFUCyNGQieZMrXDG/dGM/Nb0Pq/QmOhdNIlKzKP34HYqfJSK4OImdq/ZLeHJj/OK9ZMNLsjVruUGEWTP1Qpm6XRWY1ZIaBODIcNnog/LPPKRQQc6BP7g85hCvlV+U+qftfjZ6V2LqG0T3xtdNxfmxPOrW1Sr51RVVKnvAeQcCxqCslCb2STeEp6HnzcY+VYFQSl/1trPTu9aRMGMa96Mrpko38xKxP3WFE0dfw1kX4fnM6TsRRTbAkAKbCWKEyFLL0VIEJRiBKgW8bPj2jGt1k1PK8aDCLRaNr2cazzQQNW+xWJq31wbV92iuqBWzi98jBKVlDBOhSxttGJQBC0sOkX72fHAgEJmU96On6/1XdMKeADUYINIhJBWAx6y+d1EfzY9wpnaw39C8cSMUqRgItqg8Nsnq7kFPoTpGux+16JwpdWO8f1AVkDTtuHn8+PiedzE2tDUkGNS5T5vPNKU9xv459mcTaIsnD5Qfs7Jylhh+5zsEqbrsPNdawUrvdO5hkzvOD+fHxdNnaK/GlHWTMp93vh+3m9oPsuQTEOPUAhsKZYO/mGjg9/CWnPFrehdJ7jlRyfmT9Me40KVGOnJTZ5iwILWfCqaQPnfju065aQDMcLXHx1EiFQ0U6V9mB44An4Ka80ZI2HJT4HZbiJiLUaVCbhQJfq0cWFODFbgffDjRFBTefDjjU2pnHQghfD1hxDBtShRLwjvzAKEazu0Wg8vCdIPKWhh0e0IXbPX47Tqg3af1l+NmypFU6kofHG48FZVUaVq1biJk+dcE68rFo2pwhpvb/QZ+nlPKWudNImGAa7xyGo+vCVIP6QwCIt+JuO12q7YOYVMmmJiW9EXTSQbkJCtjsB91fj5Yk43xpgqrGU1cMafn//3lTSJlhT8ASpG64VNA+im5iRs1yZip/yWV31rzxFTeOi1ETVpAHIS3or7+H0mat3MzjMsAQcAmrdJNGTwh6hW8xEW3JxT2K5NxC2tlxd9i+k7jNpp86Zpm4raNTGYKLed8XfNNXChImzzEZAatiKE6oEoEMY5NTf0zKFcK6stLaXNo6aHVjMnUtO6BjWta3LGNVPj1QgrmkjzECgSMbVJ/COMAk0Y59TcMDKHcq2sXmkpvVQe4rl6/XCSqXTOuE4ifQuVu3KGWe2a1LCFCqsarLBq4rymWFqLcsLMPcE1aHpmSlG7xjVrolZNjCjl5lHeh5hcVyxTZeca+HjF2tn5PqWGLTxY0WKFLxeWfzSXazdfCYO3z9WUia+5Zk30ZdMT7MQxebCBmKvNq+twpMEN8X1Q9KcuEe1FRG8T0QoiWk5E12T270ZEbxDRyszf9sI5NxDRKiL6kohOtjOxsPpVGc2rkCYu6GvxkmJ5xcoJo4S1ekKQKIiJOdTE41otl975Na1rVD847tfGo0b5e6vRxWba2hLWCIjEIqY2vwhq/QqrX5XRvIr5oIXherwiXHnAvEMvSa74V9tGFMZ4DjXxmFbLZeTblkylkU5nU3fwiFH+3qpwZU7odPB9UsTcFgBmRk0C+CNj7AAARwK4iogOBHA9gDcZY70AvJl5j8yxYQD6ABgCYBIRWS4w6JUAUMhp3AxO/Z4k5YM2GpQHFBjdT2JBeDHhrSj4ifdKMpXE1u1bUb+rHg1NDXkF33nUaaExfYVCGXQQyPrllQBg5DRuVqCyO6/mINA0N/QiQbMJbvPvJzFwQEx4Kwp+WuGXiBCLRtT6pHwf70/bh6QwRQU2xth6xthHmdc/A1gBoCuAMwHMyDSbAeCszOszAcxijDUwxlZDqUR/uMvzto1RNvpihOKBGGKcfj56PlWl8JlrtWKi8KXdeFuxjqiez5u4XxxH1I7x9mICXi+w9h2YE9b8FNjKbf0ychA3YyKSGOPG56PtoxQ+cz2tmCiAiVvW9yxrNSzk76b1YRPvU95eTMDrBbZ83RRJsvgWAJaeiETUHcCvAHwAYA/G2HpAWRSJqFOmWVcAi4TT1mb2afsaDWA0AFShpeWJu4lWs6FnYjL7QCylAuhu4vSa7QrSYUKbg83oXuDCmjaKVOxH/KtFzM/mdQF4q99BmH3Y3Fy/Mv2pa1jrjp09mnVx9Hyw8nNoSV+2Qrhxvdo+SvEz1OZg0zNXZn3cmKEfmp65NXt+NkWI12XqrH8HwQUUmMG0IZaIWgOYDeD3jLFthZrq7Mv71hhjUxlj/Rhj/SpQaXYanuBWqabmKqxxgtSIWRnbyTyLaQH1hH+tn5p4XPRD48eBTDmqlvnaM/G1aGJ141ocE06TqDI1l9cvIHcNa9l2NzemaQs9DZsdYaFQmofmQNAaMS/zrxmda2Ri1wr/HO39IfqhaY9HIoSURoOW+yPC3WtxBUKoNWymBDYiqoCy2D3GGHs2s/sHIuqcOd4ZAE+7vhbAXsLp3QCsc2e63qKtx8jTKJjFqi+RFRNg2M2DQc/P6/xrRucmU0n07NYLHdt10jVZFnLyT6aSakJdvTZbt2/NqXCgN3axQBetwCiaaHkQhF56EB7kYAVCOAW25rJ+AcirychTKZhBa9IqhlUTYNACUSHCMDezgrITgVpPsN/RkERjUwqNTSlds6XReESkJtQ1ahOLRnSFNf7e6DzxPtTek2JyXv5ePAZADXIoN8xEiRKAaQBWMMb+LhyaA+CSzOtLALwg7B9GRFVE1ANALwAfujdldxEfaPzhyffV76zPaWcGvQdoIdNWsXPNHAsKUTskRkCWO2IVguoW1di4ZQN2CP5rgHFEqdFxsY0b37VokhWDIkT0kvjyv2LKEdMQIRqLmtr8otzXL9HHCID6AOX7opHsw9cMRg9RI9NWoffF2oeBrIYo6Jn4h1iBIJ1mqIxFEIsqm1HUp/b7L3aPOP2uxWAE8Z4UhTDtEGJ1DzHliPXBSztK9GgAvwVwPBHVZbZTAdwJ4EQiWgngxMx7MMaWA3gKwOcAXgVwFWMs5fbExQegFyWktGkYrPiw6fVlhmMHHG9uciFB1A7x137iplbPSl/8Wnn05knDT0crnZJTVnPUJVNJ1LSusX1dRgELHDGqVPSD05pc9YrYm4EARCNkavORUK5f2gegF9oAMRWDFR82M8KZHpu27bI0vzDABVz+12/c1OyZ7Sv3moHvftyBZCqNZCqd9xlY8WkkIrUfO9dlFKyQO0ZWmON+cFpfTcdfY4hNokWfDIyxBTD2wjvB4JwJACZYmkg0hu6de6jvN27ZkJMstBBWNQFGvmZajYQdh24n/m/vLHjL1rk9u/XCqrUrbZ1rlUJ+el4nD9ZLe+EWVvsS5/Lcw0+6FjRh9p43Qix1xZP36pXQ4mWytPe4Nk0JAMDsc5iASMiqGPi1fmXOw85GRbbjWguzWNEGFMpbpdVIWNUy2BVaGGPo0LaFrXN3NCTRqso/lwqjz8/rgAt9B373xrLSlziXvTtWFzR7WsHKPa+HWOpKmWfW102raSPKv7/F6ypHk2hoHKOSqaQlocMtPyQxNQKQ+8ASTX1mhTe7gQdOrsepsGYUbWg1WlZMV+EF2sjLIII8ROEnFo2humW1YyHLLcR5GUWSdmynBEPy4vFA9j6v15h1eXJesxAI0ZAJbH5CRKYFD7f8kMTUCDw2IjcFg2LuMyu42RVanFyPU2HNKNrQSDAzFna9VZzoRV0WEr69QBR8GGNIpZljIcstuHZM/EyI8r/fxqaUWjweyL3HRc0eT85rcRaBac/MEBqBzQpuP6iNhDGzWeG159ghyNJWha7fLDwDfxh97dxEKzTW76y3dT8WO8fJPS4KtGJ/QK6gVtO6BslkMuf715bUsiSAExCJhnexCwtuP6QLCWNWTURONGxOzneC0fVb1ziF08/OTbRCYzSin7qjGGbOsXufG/nO8ddcWEum0ogQ91XL1chl+7I8PMoirUfYsOPnY6bKgVPfOKtRoqKgxtM4aPvRSyobFHqfFxcu/Ag88MokahZtgtuO7TuhY/tOpr+fQoI5/yy1fmVW5yf+1ZbIEttxzaB4TVrh3co8FB+2iKmtuePEx8fKPqsaBqtRoqKglkoz3WNW+/UavehW0d/PS7w0iZpBzK8GAI1NaTQ2pU1/P4UEczEARutbZmV+Yl9AfoksfpxrBsVr0grvttLUhNiHLVQrpyiYFEqpYUfzxc/j45hpo9Uu8IepmTFExHPE19rrEE1rWsGkmDbGL4xSSwQtUPo1vjbNxsafNmD9pnWmooPNaM14MIMdtONxAZr/X9IGQYjmUm3lBL3+ikFEqIhFTW3lCn/IFEqnYeshAuQ8rMy2E4fhD9Ji5+rNTYzAM2ormta0x4qlcPATPaEpDMKkH3MQU2wQESorIqiqjBp+P1qh20wkMA9osIr+jw/lL///JAqMoqmUzy2t+dFgGUKoBbZwqG2Q/zDz0i9IKwwV8tUqlPRUD72+xAewF5npw2CG5Bo27sxu53NyOr4faOetNSEWmpMZYc0JouAsajzF/0tcIBMLyGvvT70gBbOELejAT8QHmpd+QVpByEjroRXwzJhHjfriD0YvMtOHwQzJP0fuzG5mTm6btv34HPJ96ArLHtp7rRBupPLQ3n/8vfj/SQw40Pq3RSL5qUBsTMTmFXhPaAQ2PS2NXw9hN/KfOXnIOSWoCgvaiEI3AjP0UmEUa+/EfGjVT0/b3ssgC6vkafRiueZRHlgAZIvPawMNtGZTSyZRomYtsGk1NX4JIm48SB0/5Bzgt+O93thWyyQVizQFtJpO477tXr/V8/SjVC0P6xl6GrE0Y4ggN2qUC2pcaCuUO87650oIsw9baAQ2IDhNkfbBZEcICFLL5XRsUStjJQrTrnO69rMGsuY7UTtnRggudu2FBEkn37E4bzPmaj/uj+qW1TnJngsFRIj1TEVEfzY7c/Y5x1roCDKXl5OHVZBaLqdji9evjZA1K8xad3PKd4zn2irt+E4E6kKCpJNIXu2cC+GXQJ1Ks7z1Iz8ogn/X2bmJiP5sUsNWhnATnvhg1Xsw82PliJEjv5nrNRL2rI6td54bn7fbJmirmlQv7hmjz3nr9q2oblGtatK4Txp/LwqXhTTadoU15Rdv89WwBQHXNABa02f+wzkMpkcvMDLbmb1epw95I6HPjc/bbRN0vhBZ3M/L7fvG6HOORSNqFQYx0EDcJwrFWmHNiqBeqjT71VUviu6wA45Ax3adcGTfo9T9bj54eXCFNrBCTBxcKoifjd8CbRCBDnp+acP/NAoDDj9W3Wc1YtTKdWgjSMU+eKCM1idN1H7ydmZ+oFiuJUqEioqoqU3iDuIDXXwQb61vRGNTClu2N7hu6uQRfNrAip0NpfmDVnTE95swBFys/uFn/LhtF34UKlVYiRi1cg3a6FFtVCg3c4rwigziOcV+nPCxbBHi0lQlqWGz4i9lBu1DePGKDwAAP/282bUx9MbTBlZ4NZ5XaDUxPbv1wvpN6xzlY7Oi3RHb2bknzJjCtfv0EtF++vxSW8mL7Qi6efU/M4l7Re2YNvUIP6ZnbtZLmsv7tOOfJ1N2FMeqv1Qx9ASNmupKAEBFzP3vIxtYkTuuF2N5jaiN2dGQRIuKqOmgg2L9mUFs64YfXbHEvHpJaHdvU4XWLSpMj6udu7XrzX2vJO7NTdehNWvnjpVr+tZLmpvt0/Il6U8yRJSkwOamsKYX1cg1C7FYLE+o8tLBX/Q/KgW0n4MVocXIt8vuZ+tGGTEzZlntOBu3bMBPP282FXmsxcq9ZKSF27p9a04/egKZXp488bXWJYCbVq0KbITmHXRgFjeFNdEkqn0wA4rTtihYeWkuKkX/RfGzsFpxQc/U7OSzdaOMmFGksNEYlRVRVMQiuue5HSyhp4njBee1wp/WLzE3Vx5T/yptctvHopTz/6KcKEmBzU30HvRibiqtYOGVsBZUpCdgX2OpjRLlmOnLbROzF/2aGatju0451QOszMPKXM0EXugFDIhCW6GKBlrE/wMwW/qcSvOhXcoYPeRVHyCmLfXjzfcTpN+QE42lXv4us325bWb2ot9i42jLPIm4ncbDqL3W90x7L4lCmzbvmhFiO4uTDLWGLbQ/h/3yTzKqfhCLxtCt017qPlXrZmNeepF4Yp+A8mA9su9Rlvt2g0ICljaRqmEKiFgMrWxoZdzCr5QqWp/GHbvqcfLpp6Nzhy5qG6N7xKpvmxOM8gdWt6jOC7ThaH0qLSfOzWjYzGzljl++SUbj8Afcrkwxei6Q2J2XNoGq2CegPFS3bA/m/34hAUvrx6TN6M8f+umMT16Qz2o/0qrwMbJm7Qi+/2kHGhpTaGjM/jLTT2Jr/t5xev/rfU/qd5XRnmnH0PpU2p8Dmdz8J1QaNj9TIBSCj792w7fqPjE4wQgjLZmR6U2bZX7Rsvdszdeqaa2Yn5aIKMwVGsNOomO3tIpBBDuIfmGvvfSSqfl4HVnKtWaiTxqPFAWU71J7z4lC2dbtW9GxXSf89PNm1Y/NktAmNWyhiVDjc2hRqWhPxKg7IwrNXc/sp9V2tGttrzSd03xihc7XK1WkRywaAWzEwrj1fft5z+RqXIE927cyNR+v08TwCFDukybOlQvekQjlJMcVhbJYNILGphQqYhG1qL29agfB//81IlQCm/igKmaycXtM8SFc07oGW7dvRVVFlZp8NJnMFbKM+rIqEIlzsGuadGpasyOEGrW1Mp+gBXM3qG5ZjfZtdlOF+2JBDIVwmgONa82095DePaXn8wZki8Pb8ackIlTEQrWk+E7WvOPPum8kuCRT6cwDK5N8NNOsWKF4oxQgZoQ5u6ZJp6Y1t/ys7KQ/CYNw7oRUmqEpmVIFe8B+Pj+rOfC0cK2Z3j1kJHhrfd64eddRoZEQf6ehXV29KOEEmNOCcW2ROAezD1En+cS8uuZiFLo2t/ysygntvWKkXbT6eZj1kzT6USC+5lo2ox8BotaNt+faOG7Wthp4QAAi0fAudn7iRQknwLwWjJfyiURIzRRv5gFqpm8jvLrmYpj9TIpR6sKXWbT3SaEyanY/P6tCtNbkzt/qlZ8CkKNx4+2JcoNw7AcehPc+KCqHElGCiDYQ0WfCvieJqC6zrSGiusz+7kS0Uzg2xcO528LqQ1SMsgsrVvNlaXHj2oLIiRYGtH59fmD0o4Dvr6qoUrVtWl81o3xt/Bg/n/dh6d4gJa2Hmc1Ud0RDiOhLIlpFRNfrHK8hoheJ6BMiWk5EIw36KZs1zJ6pKRxmWiNs58vK4Ma1+Z0PLSxofcX8wMjcqtVOc42b1l/NKF+bWL1BTLBr7f6g7ODFtuLX6cr6JWLmSTMdwIMA/sV3MMYuEAa9F4CoYviaMVZrot/Q4JaGKSicaubCLpCawa9r0I5jRhtmpp3d8fXgWjGteVasJaoVMkXhTOvbZhY303oQURTARAAnAlgLYDERzWGMfS40uwrA54yxM4ioI4AviegxxlijprvpKOM1rJgpL8zCGuBcMxd2gdQsflyHXgSm1XPcHF8P8bDWPKvNuSaeo7TN922zjis+iW6uXypFV1fG2HwAuhldSfmEfg3gCfOXEj6K+XD5oUUpdQ1V0AKfX+NrhZli35uZ6FUz371ZTa82IjQWjeWUp+L9iHnY9Prkx63el9EImdpMcDiAVYyx/2YWsFkAztS0YQDaZNah1lDWqbyLKfc1rJgWwQ8tSqlrqMIg8PkxB60wY+Z7KzYvs9+9OWEt917mZamUv7km0ELz48dt3ftkciuMa+uXiNOfw8cA+IExJmZM7UFEHxPRO0R0jNGJRDSaiJYQ0ZImGAqUgaNN41AIJ0JXkAKP3bGtCC5O+jaz30+saNis9leojVE7vbQxyVRS1aaJ5k9u7tSeW1VRpb7mJlOrqVIUM4XptB4d+P//zDZa011XAN8K79dm9ok8COAAAOsAfArgGsZYGtZwZQ3buS28VUrMmoWcCF1BCjxOxs6a0dwVOAulWwkaqxo2K/3Zaaf9TLhgR0RCdLMYOZp/LmO536VoDrV+jaYltkJrmCfrl9On34XI/WW6HsDejLEfiehQAM8TUR/G2DbtiYyxqQCmAkAbahf8XSxg17zmJMIvTDgJsHALs5GrTiMrncB9B1u1qM6rdmCE0xQseuhFOhvNVy+QRowsFfepQQcmE+cSESpjpkvcbGKM9SvUnc4+7TpxMoA6AMcD2BfAG0T0rt56UwBX1rBOPQ8K1RpmL0rPOEq0lLBy7V5dZzGNp9bUxz97Pz937juYTKXVsmJmtGlW5mhWq2a2rRiAIH5uevtsBR0QYKFOaKE1zJP1y7aGjYhiAM4B8KQ6G8YaGGM/Zl4vBfA1gN52xyiEXU1LofO05iKrYzhx/teaXf3UJGnH0vo+mSGZUrL+e41W66OXnsLLz05MLssLqzcmG0zfM1Y+W21boyTPhd7zz4trzbTttH9FrPpGEghRipjaTLAWwF7C+25QfomKjATwLFNYBWA1gP1NzzfANcyupqXYedwEJJqEzOLE+V9revLTmT3fn4kMjxlBpGT99xo9jQ9/LwouXiEml+VF1RUzI0zdM1Y/22Lt9bRren1kBdrcdoU0pLZ9I90xiXqyfjkxiQ4G8AVjbC3fQUQdM852IKJfAOgF4L8OxjDErjalWD40bkoSk45yX6BiD9j6XfWO5mXX1OZUQDGTI87UPGL+R0wCxln97fZRCG0Kj6qKKtTvrDeVWFlLodQcem0LpYzh96u2Db93tQKtaOoXhTp+jl2ikaipzQSLAfQioh5EVAlgGIA5mjb/A3ACABDRHgD2g7X1JrA1zK4mxawPkJhclPsBFXvAKg9we9n2tedZ6cepgFJoHGvaoOBMluK4du8NM3PXpvBgTPE75QlprYyvbVdsfL1+tdo1vTb8/tVrK258fGdRx2altaKfkSfrV9GnFBE9AWAQFHvtWgC3MsamZSagddQdCGAcESWhGFKuYIyF17mjAPwBZxRVpyVIU6if42of/uIc1m/S/oBwNo7Yt974esftYrePhqYGx0meRQHKblSp9jvR+rFp24lj8/38B4rYzpoPG5lO2VEMxliSiK4G8BqUHPQJxthyIroic3wKgPEAphPRp1BW0OsYY5t05tUs1zD+gNOLqtPitzlOxO9xtQ9/PoeqShulDgqMwdEzf2rx05cs9xznSZ5FAcqu4Kn9PsTvpVC1CrEd/3EitrH1ubrwXbi5fokUFdgYYxca7B+hs282gNnF+gwrTtJ7mPW5KnWMhCWrgobZ6g9m9wf1GfN8ZXrY8QV0Q/DTlp3iEaJ6AqE2WtRy7jUNJrVnpmCMvQLgFc2+KcLrdQBOMtFPs1nD7AoEbvsuhRk9k6Qd371Cn4nV/cX684piwpoVHzTta+tzyfflE4MKuAY497h6dkZoc0P4dec7cGv9Egk+5C4EaM1JYoZ3qw8vuxqgUhPqtL5PVnN3uf2ZBPHZVbeoRiyWXwqKYyadhxvz1gvC4AEGWt818Ry9mqNiol31OzXp2hMhQmXUdNCBxCW0WgYgN8u7VSFE75xyFOhyowqt5+5y+zPx+/NLpxnSjCFWQCvuh4CpZ+7UVjLQDsP9NLWmf0VjmL2HbZm4Q3wfu2O/CJiO7Tph+rZ/2D5ffGjWtK5Bx/aK8/wLPz3uqC8jZtdPt3VemND6P/F9bo8RZup31Rcsem/W97Fzhy455ww4/FjDttq+OrbrhKqKKnRs1ynnuxB96rQaPNFXk+/j5zj7TsnVSgfNicamFJ57bYWtc7W+P8lUGg0ZB/rXFn7jqC89XnxrZd6+UhPWgPy0D25fQ9g/k0iE8nzatHD/sUI0NGZ/0aXTSnTmj9t26bbV9tXYlAJjyl/xu9D61On5R2p97kQfTPE7tfw9KGq64lsAlJyGbXb9dOxM1SNKMVRGqhBB1gTzXP2/cXb1b033ddvM8bh1+M0AFEHt/nW34YnfzAUArHlpNf676Dt3J5+haY1+qpXHlz2C3/S91JMxw8CLm2dhfeQbjG53He7deBN2JHfg5s5/93zc6hbVuO/7WwAAO1M7cc3uf/VsLK65amhqQCwWQ6uMTxgvqs7bJFNJdOu0F9asX42+Rx6C9S+tQywaQ89uvdB9YE8s+PAd3f47tu+EjT9tULVfg4adiNdnvoSa1jU5Y4jj6L33QqNLcNckWo68+NZKpJIpRCIRkOaBU1FZgZfeXoXTj+tpqq9Pvt6Eg/ftoL5PptJ4ce6X6H/4PgCAllWtsPvuLV2+AqD9bvpBKQs/WYejD+6ie6wceP29b/Dzlu0AgHNP7YOnX/wUqVQKw86q9XTcdJrhuVc/z7xO4/zTf+npeFx7RQSkGUM6xdSi6vw4EWFXYwotq2L4qb4Be1a2AmMM23c1AQC27WzC7m1b5PXd2JRGZUVE1X6t/2knuu7eCsk0Q6XBXPR8AL3T6JoLAQ0KCioqRqQNtWP9MCDoaUgc0n1OV5xy4DF4rm4uvj+voO9kydJ9Tlf8vGMHfhz2U96xmtY1SCaTaiAC/wsgx/TI/9a0rkH9znpUVVShVYtqbNyyAR3bdcLGLRtyBD/RT44LbNUtq1G/M5sMV+wfKC6Q6fm7iaZUkXl4eWmBfEMqnffvwkZNHWXqc5xw7O2m+iwVOvU8iJ1/z9NBT0PigMS0N1HZpgaN235C/NITg56OJySmvQmkU4hfpu86lUylEcn4gnGhTWt65H+TqTSimZQgyVQalRVRNTVKZUU0R/Djf7nAlkozteKJtn9lX2GfQu1cRLOoXjqPSeccaGq9ibTvzqoG3WjikwR2PT/a9zWs5DRskvCy8ddbMHnXLMfF6J3QvXMPAMCa9atd73vh0OWY9886DJpbq3s8mcw3M+qVgOJ/uTlVNE9yLZnYVhS8Nv6kHK/fWZ83htb0WUiA086F/zUKoDCD1LBJSpkRI49HJEKOi9E7ZWdDEi2r3H80J26foghrt1xl2CZCudUFgHzhif/l5lQioDLz/17UxCntssEAynHlnGiECgpnWp827WvjvyY/jIKEV8MWGmcSMenq++d8EeBMrPP+OV+gekZrAMC8wXV4fv1jWD7WvRQXpQIXOpwWo3fCmvWrPRHW5g2uQ6oxZSisAfrXXb+r3tTnYTVgQwzyMBKy3Io+NQs1cx82rl1ITJga8EyskZgwFYmH/qO8HjcRiXETMXfR/5D428yAZ+YvXFBxWozeKZ4Ia+MmAowVFNYA/Ws383lYsdRpU4AUErLcikC1RIh92EKzcor+Nw3bdiEWjWHZ5d8WOCM89H92f9Rfovg2DJpbi7M6X4Q+D/jny7Hx/+3E0uFf+zZeMeYPWRb0FFxl3uA6AMDAV/vmHdt9VnvT/ehVJSiWLFdvX6HUHE5w6tdGmShRM1s5omoXUk1gjCFx94xgJ2SS+I2jEb/8FOX1LVchfstVGHzk3oj/abgv4yceeBaJOxO+jGWGxPjJQU/BVRLjJgIA4jePyT/28Oum+9GrSmA1Wa5eChW3vLLc8msTAxcKbUEQGoFNZNDcWiRTSfR9aK/ijSXo+JeWOHTmvkFPQxUkOh7QCV9dtxEzt030bewBhx+Lg3vWYvGwlVg8LD+KzS5cWNPTrM0bXIc2rVoZniuWsQLyU58YVSYAcoU77WttP1aSOxfCDS2ci5UOSpb4LVeBiBD/8yVBT6UkiI89B/Hr40FPIyt8tKpB4u9P4PnX/bX0/LhtF7btaETijmmu9akKazqatcS4iUCR/4tiKSu9xLhawYVHiQK5Al6hFCpiHdBiFGoT9qhcN5A+bC5wU+JW3M3uRuOonUFPJVD4A/+A+/YEAAy/q7D63U14VOVhq3rZ7uP6yTfizjETAGQFNUBfWJs/ZBkqKiuwZqhxJLFeyo9iCYG1PnDFfM+K9S2e77VZlMi9SgcS//h09WYsfu19xK84LeipBAp/4Mf/+JtAxudRlfEbzAXu6PHZms04qPtuAIoIa+MnA5Eo4qNOKNifNu1HsYTAovnUjO+Zmb71zvcK8mEMJ4RydV19o1IJZsbShwKeiTnOfH14nrA2ZtUwAIqjuqQ0MCuszRtch3QyjaPn9PF8Tk6FLD9z2SlBB83Xh42TuE+pJf/esvUBz8Qci2c9pyuszXxqCYCMs7qkJDAlrI2bCLA04jdd4cucnApAYRag/CaUGrYeE5SbrseN/TEQ+X5DYeOwWflanck9ZwGALw91AFg8bCXqNynamUKO8X6w+sbNSIy5BSeMuDTwuVilkAnUzPFmDVHZmzvNEL/2AgDAFy88j6P65vsNhQ0jjc7wXysZC/x4sCfumAY0KclWiznGez6X+55Eq+77YcenCwOfix2KCmsGxyQIexq2cApsHD0nb4k+ekJjUPSYsBvGT3gQ1cP8T+/B/dcOm9UL49f/wVJi3nmD6xCJRQzvOz+EtVIrUSZCAKJU3tozK+g5eUv0cWIGdBsucCe+WBrI+Ik7piF+wyjMer7OclLexLiJAEX0Awx8ENbs1GQNFwQKscTWLFbXy8ePdbU/Xjro9lkTkPhgMuYPWYab116t69R98umnuzp2KRGEEHnYrF7quFaENW66DlJYA8JfjqsQRITKWIWpTWKer74zLn9mB142qG7VJrz/6Xokxk/GrOfqDB26v/9ph6vjlwpBCZF8XMvCWsZ0HZSwBtgsBRUyZJSoTbi2ZPhy80LP/CHLwCbnCk4jPnb3JuUO7jcNuxHxI8Yg/tm5GN/tQd2H7SlXHQcAWD52nS+52aTPnHUGvtsPSx59WppBHaLkYZNRohwe7TfjiQ+snTd+MhKTX1LfL5j+pKvz4s7ttT07oP8vO2PYny/DsLNrDR9Ca39QUhb5kZdN+svZI/GvhWhTe7Q0g7oBmdwCINQCG9eUzOzzUpGWWVbdthA0Jldw6v/s/q7OK2/MtcZpJH5/yh8BAH0e6OJLbrbqjsFVGdCS2KpouK5ZfXEg43PhlVc/0GPFtd/jvoeeNKwvynPKSWHNHFGKmNqaA1xTcsmFR1g674CzzkJ8TPZHavzG0a7OS0urIola++2vJDX3JS9bRX79yaB47tUVAICZTy8JZHwuvO5sMNa6J+59HIl7HwfWfalbY5TnlJPCWnlQditn/Ijm7TdSO6170FNQidf8AQBwf49/BTL+U/ffBaBwmaqGbbt0c9jNG1ynRoNyYe36yeZqzDnBSsWDsEEkNWxu0P+XnYOeQmDErxsR9BRUzh5yAABg+PnBlLyt3E8Zt2Dlg2QjkGzUzWOnRoNmhLXP1mz2ZJ4iZpLphp2SNokSUYKINhDRZ8K+vxLRd0RUl9lOFY7dQESriOhLIjrZq4mbpaZ1DR5f9kjQ05AEQCFB8R//uReAvoB7/49/RUWrClS0qsjRrPG0H2awK3hZ9WGzO45XgmGEIqY2PynlNSyZSmPhJ82vzJ2ksKC45AulMlD8uhG6Qu7TL32q5FkTNGs85YdZ7AheVoUZu8JdqQuFdjGzck4HMERn/32MsdrM9goAENGBAIYB6JM5ZxIR+fpzWlsWadVfvsdd5zzo5xQkBmgz//tF3ag1mD9kWU6N2kPuP1G37dLhX2P25y/h6Dl9HKVk8St4wO44XsyPQIhFYqY2n5mOElrDxNJI/5r0PH7Zq4Ofw0sMELP++0nirumKX6NQo3bZ07ON29+ZwM+btzlOx+KHFsnuGF7OLcwatqIrJ2NsPhF1N9nfmQBmMcYaAKwmolUADgfwvv0pWkMb5dfxLy3xCer8Gl5SAL3M/17S+E/Cn5v+iLP+MzzPB00vGrRu1Bq0aluN+ccE47NS6hDgu/bMDKW2holRfvGx5/g1rKQI2qz/XpOYNAf9zxgINNTn+aAZpYxJ3DUdqGiB+MVH+zDDMiTkedic3IFXE9GyjLmBV8DuCkCs2L42sy8PIhpNREuIaEkTGh1MQ1JqnFt3Inp28z7lR+XvGE66QnH4nje4rqgZsHZad7WsVlDoFYj3qm/XzyEKpUm0AK6tYTu3ee8fJAkP0x9fhB0FggHcIH7lUKxZ8xMAxR/NjBkwft2IwEprAfo+bG6aL632ZdWnjiz8CwK7K+dkAPsCqAWwHsC9mf16V6H7aTHGpjLG+jHG+lWg0uY0irNw6HL1oXPX7DvR+E/ytSi5JJ/ZtW8UjKx1k1Orz1e1a3OP09eczR+yTPVpCxqtqdJN02WhvowEM6vjl5DA5uoa1rKtNf8gsyRun6I+cD5euRGJSXN8L0ouyWfEb44sGl3rBqcc00PVrj06fpJhu8T4yapfW5DomQvdNB8Wqjdqdj7FBzG5BYCtlZMx9gNjLMUYSwN4GIrJAFB+je4lNO0GIFCP2bPmH6c+dK4793pU/o7hhT+/px6vbhGeNBjNhVg0hmWXf1u8ocsYpeYY+GpfNf1KKeNEG+eGYEggREz+C5pSWcOGX3eZ+sD5Va+OiF85FL+q7aYeT6ebp/N10DDGkLh7hq9jFkrNEb95jJp+pVRxookr9WS9ZrG1chKRGHd+NgAefTUHwDAiqiKiHgB6AfjQ2RSdMW3l/TnvY9EYjh+mmMn2fKYDWkyvRCwaQywawyv1TwcxxWZHMpXEbWPyw9C9xMs8avf/+FfP+rZCGKoklIqGrVTWsHeWrM15zxjDD5mawYlH3sD0aW+oZh/GGP7zrnEKG4l7EBEi3fbzbTwv86g9/dKnnvVthbAIXWEOOjCT1uMJKA63+xHRWiIaBeBuIvqUiJYBOA7AtQDAGFsO4CkAnwN4FcBVjLGUZ7MHkPggG1E1b3AdVlz7fc7xdr/NjfRLppIYe4KSH2yv3fbEj8N+QjKVRDKVxKnV53s51WaPGKU5u/aNAGeS69MmzssOLaMt3ZhSyRPWKNGwr2Hvf7pefZ0YN1FJhJrh27mvaK8Fhx2Q0aRUtED8spNyHiKnHGOcJFriDDFKE1DMokGh9WnTzs0KkUjwP6DCRJgFNjNRohfq7J5WoP0EAOYTVjlETJQ7aG4tMDf3uF40INdELD7+s7xjEu9o2tGI5WPXYeMKxdciqOoBvNTU3OOWYNDcWjTtkEEvrkChjRIN9RomJso1Ew3IHxbxS47xdmKSXNKpbHmuHUrEexAVBHipqUfHT8qOn/b0N4UkJJRuWnUDOrbrhI1bgne+lGRZOHQ5mnY0AUkELqyJY/NqBu26tlOPde/co2BlBD2u3XOci7MrXbgPm8QZjU0pVFbIahBhIXH7lKxAtCObmiiock983MS4iarwhqqsL/bOhmTh6ggazh5yoKvzK3nCYZnVpexWVyms5ePU5OcUnoB20NxaDJpbi9323T2vDReeAODIvkfh8vFj1fdeseOFlXj+lJk51Q6sCmsAUL+r3sVZZSnFgJhS8WELM1JYy8eJyc8pPAFt/Jar1A0tWue148LTlu0N+Oq7rVlhyiM6H3cq+o+8CP1HXpRT7cCKsAYAkYh3EkopBsWEOa1H2WnYJPn0f3b/QMaduuUu9D4vv7JP34f2ytsnatwWLXsPi5a9h0GozWvnJqdWn4/Gf5KaVOaa1Rfj19dcl1PhYPz6P+Dmzn/3dB5GeCUIekVYE+dKSp/4jaN9H3P2K8uxdck83WPxP1+Svy+j+WrXugrtWleht8cauFOO6YHEpDkAgAOuHApAKVTf+OUSVcic9Xwdhp1V6+k8CuGlMOgJZZw4VxJyFg5dHuj4N3ZXUltte2IFIjHrt5o4f1EDx9+7QeXvsr8A7+/xr7xyVEEJa6VJySXOlYSYxO1TAh3/jBOVH7p7nXg6YOOeFeefY77UeW+X+JVDEc8Ia4BSf1QsSRWksFaqSA2bJBCc1MJ0St2oNdjyTR0AYGjHCzEQ+cEfxRDnr/V5+/GFOlzQ9lLPU1ksHrYSh83qpfu+mL9bLBoLRaoNjtfzIaIg6oRKyhSntTCdkLhrOtCgaLhP7L8P0F+/FFQhxPlr/d06DjgJgJKmxcuIw8Qd0xC/YZTuezO+bl7Pzyq+zCc8l5uH/KkrMeS2meNtn7vlmy3474T30WbPNi7OSCEWjWH3M2sx7UPvK1aIwhoAXLokm/qlmL+bFeHIzRJURvghPJZK4lxJ+fPJ15vsn9xQjwPOPgeodD9lD2MMGxe8jo0LXsf7n35f/AQHiMIaAJx3zW/V12Z83cwKR26WnyqE18IawVxKj9DmYZM0X24dfrPtcwfNrUX8iDE4dOa+Ls5IgQse+1x3hOt9F6PNbjUFj+vVAhX3FRPMjAINtP2YGdtvuA+bNIlKwsDB+3awfW78lqvQ/5edEb/e/QTf4sP+i+efdb3/QjQl00XbFKoFWkwwS6eZYaBBsXP9EvpKGblySlR2n9Uey8dmq/A48RPj2rkDexzkdFp58HkFkRpkwYfvFDzOhUkuPPGkzNrjIqKgZRRooO2n0NhuYk0IlD5skuBIPPx6Nk8a4MhHjGvnft7Z5HheWsR5+Z0aZPe2LYq24QIlF6BEAVNPsyQKWpEIGQYaFNNKeaG1slr8XZmIyS0ApMOJROXHYT+hD7qo7we/3Q9J2BMCuHbu89XuJicOUlizgp7gVqytHkH7wdkp/i6RBEH8spNy3o+8+UrbfXHtXJuWFY7mpIULa0HlcLOCKLgVEqaKCVpB+sHZGTeogAIzyNVVYkgYHOY7d8gKkIWEtXmD6zBj6UOmTId+YkYzpkWcv5FGjrdx61r1PjereeB40IFbpamIaAgRfUlEq4joeoM2g4iojoiWE1Fh9aekWREGZ/mGxmwFgmLCWmLcRLy3bH3oTIN2PkfxGrTni1ovt65VT5NmNwecWz5sXqxfUmCTeEr3OV0dnb9+k2KiLSasAYpP25PbHsGT2x7xPOmul5gxffI2bgnVeoKl1TxwBPeCDogoCmAigFMAHAjgQiI6UNOmHYBJAIYyxvoAkMWAJa6SmPamo/OrKpUkyGaENUDxaXvp7VWeJ931mmIaOX7cLaFaT4iylQPOrDm0SNderV/hUUVIAmH+kGVqvdWnvnwU1XtX47SWv3at/wcG3IszMMxRH2aEtUFzaxGLxnBu9QjlvcdJd81ixiRaPpCbJtHDAaxijP0XAIhoFoAzoRRl5/wGwLOMsf8BAGNMljlpZiTGT1brrc7/aC22b2/EqQN/4Vr/3Q7o6bgPs8Ja/JarsubD43vptg2CsKX2KBE8Wb+kwNbM4cIaAPx6v5Gu93/Gbv4Ia0B4hCI9vzM3fdHCLARaENg6ENES4f1UxphYf6grgG+F92sBaMOCewOoIKJ5ANoAuJ8x9i9rM5aUMlxYA4CBh3Rzvf+TjtrH0flWhDUgHGZcIF9Ic1No0wtmCBMW5lVoDfNk/ZICmyS0WBHW9I4HFZhQLCrUrf79CEiIRWNAqng7gFsKTC92mxhj/Yp0p0XrlBIDcCiAEwC0BPA+ES1ijH1ldhISiVcUEtaKRYomxk0MNDChWGSoG337obnz2B+w0BrmyfolfdjKGDN+XF4Xht99VnvL50x990FHwlpzwQ8Nm7UxCEQRU5sJ1gIQi852A7BOp82rjLF6xtgmAPMBHGxhwpIQY8aPy+ui8ImHX7d8zgfLv3ckrDUX/NCw2Uly61LQgSfrl9SwlTFmBJq2Xdt6Oocfh/1kqf1z9f/G942KJtmJsFYuwlzQaT2sEkXUra4WA+hFRD0AfAdgGBSfD5EXADxIRDEAlVBMDve5NQFJsJgSZqpaeTsHTaqQYrz09irs2rFLObeAsFbs2spJkCs5Hzh3purJ+lVUYCOiBIDTAWxgjB2U2XcPgDMANAL4GsBIxtgWIuoOYAWALzOnL2KMBVcQTlKUPg90Kd7IR86uVkqn6AUNaIW1zh26qFGkYcJNIctNYa1YuhCxnZ1x3SyKzBhLEtHVAF4DEAWQYIwtJ6IrMsenMMZWENGrAJYBSAN4hDGWl/hPrmHlS/xPw4OeQg6nH5cJUjgtP2G4nrDW0JhSI0nDhJtClpvCmhn/Nydzd2sNc3P9EjFjm5gOYIhm3xsADmKM9QXwFYAbhGNfM8ZqM5tc6CSmmDe4rqAJV0+zVkhYWzh0ec55fuB2/je7/Rmdx1N3eFYxgdzLYQQAjLFXGGO9GWP7MsYmZPZNYYxNEdrcwxg7kDF2EGPsHwZdTYdcwyQeU8iEa6RZMxLWErdPKdqnF7jt82Wnv0LnmFk/HAuILlU6cHH9UikqsDHG5gPYrNn3OmOMr+qLoNhnJRLb7HPE3tjniL11j9nxWTt6Th8AwBPHTnM6NdMU01rZ6c/ueUFBiJja/ESuYRJfaKNfu9ROdYP4TcrvhAGj/NUgGgk7dgU5IrJ8btDmUzL5LwjcWDnjAP4jvO9BRB8T0TtEdIzRSUQ0moiWENGSJjS6MA2JlnmD6zBm1TDMG1yH6hbVeHHzLBz8+n45beYPWRbQ7HLpMWE39Jiwm+6xFU+/gq9uXWCr30+eXupkWrbQE5hKyQ/NCWFe7ArgeA3buW2zUTOJTRLjJmLmU0sw86klSIybiHSa4fX3vkFi+vxsm/GTA5xhLvFrL9Dd3/bQY9Hn3HNt9dmppnjtTy/QCk1BC1G+4VLiXK9wZMMhohsBJAE8ltm1HsDejLEfiehQAM8TUR/G2DbtuZl8JVMBoA21C1ctjjJh0NxaPLzfM8BxwMsDFqLNr/dX9gv+YWIetrCyZ+VeOPuY39o6d9Gy91yejbeI/mNhKrFlBgIhSuHzxymEW2tYp54HyTXMZXgiWQB4dPwHmH77JOHoQKWNkIctrLRo1QJH9NnT1rntWle5PBtv0fqYha3MlhnCLJza1rAR0SVQHHkvYplvhTHWwBj7MfN6KRRn3t5uTFRiD/7w/+H5JRg0tzYU0ZPXrL4YAHBwz1rDNqtv3KyaQnkggh1KrUSVqImrqqiypZkLUtALo0nUCLmGhR/RZ6nD0ScifstVoYignPm0ki912w5961DividVU6gaiGCDUitRpfUxY8y6AFSKQp5f2FrZiWgIgOsAHMsY2yHs7whgM2MsRUS/ANALwH9dmanENjdtuwadzzscAxEObdr9PZRkzp+sqjNs02PCbugBfROpFcIgoNqloanB1nlmfOm8MdGG0typi1zDSoszfz8KL9z/KHBCOGTn4ecr+VLbtqrUPW5kHrVKGIRTJ9hRVhUT8LxOE1LSGjYiegLA+wD2I6K1RDQKwINQSim8QUqleR71MBDAMiL6BMAzAK5gjEnnjoBZtHQh/vPwpOINJWWPXpSomxo5N6NEXZyTXMNKnN3aVKFlnyODnoYkBGjXj+akkSu6UjPGLtTZrRt6xxibDWC200lJ3OWUDwbglMuuDIWGzahkFDePFtK6hZXx6/+Amzv/3da5Wp81r8ta6Y3tVt1TxRc3HOZOEbmGlT4zHv8A6f9+DJzzq0DnUSjic9uORkONW5iZ9Xwdhp1Va+tcUdulp/ny8scZH09PgHM0bngVbCFcXSWu82K/eaEJLjAyUX6yqs5TYY37sp1bdyKWXa5UUujZrZcrfbeKZbOtW9VWmRGMvPJJMxrbiaBozoMtxCuiJJSMuOjIUAQYFPKh81JY44Li9McXIXH3DADAjgZ3ftBFo7mBQlY0VhZyKlqak5OxnQqJYY50lwJbMyAMjvfVLaoLHr9m9cVqMIKbLLv825zrn137Bg55pAcA4OCXursyxh873q6+diLs6J1bWqWpCBGKmtokEisE7XyfThcXOHgggpsk7p6Rc+0jfnMkRv6fsk4+Ndud8c4/45c5721XCSiQwy3MfmF5hDitR8kKbDwz/phVw4KeSqjhxd3/uCYe6Dzqd9Xn7RMFqft7/EsNRnCTvg8p9Xc7HtBJ3ccFoNm1b7gyxtQtd7nSjx6lI6zxdSy8v07DRGLcRHWb+ZT7D/pyITFhqlrg/bFn/M+pyIlE9O9ZUZjigQhuEv/zJcqLVjXqPi78jPiNOz59s19Z7ko/RpSSsGbWBzeoaypZgY2b1i644noc2CNbt6175x4BzShcvH/OF5g3uA4N25RCxFP2fzLgGeXjZwTnjys3qcKrlprWNbr7zdI5vY+j88sHKqm0HkEimtUav/gAP+9sUt/vdMnUVcokJkxVhKFUk7IBuPCcQwKeVT6+RXHu/FkVXPVIptK2u27TrrXtcyX+UlqZOTXwB/7nyNZLXbN+dUCzCRf9n91fdfC//8e/4prd/xr0lFSqW1TratzcZP6QZUgns4tYIR++rdu3OhrrjN2klpcj/dPMY/Swb1lV0suyK8RvHI3EuInqZ/T0S58aarmCIJ1mns4nMX4ywLLrVzH/vVjU/o+gk46SPzhFwqwRLLufujevvTroKQTKwqHL8wqph0lYA/TNo26z5d+Kmn/Q3FpUtKowfZ4XfnTNBQJAFDG1SYyZ9Vxd0FMIjMTtU1RTscj5p//S4Ixg8Fp43GvwqQAyQn3Ems+nF750knBQdivn+G4PBj2FQDnp9cMBAOcvOQmRWPGvVzQHhqWuqB67z2qvvn5x8yx0n9NVfc9N4kuHf60Kq20vPEA9zgvBm8ELP7rmg1kPtvD+gg0Dw86uDXoKgTHiL5cDAIbfcAVgQrAXTYFhqiuqJfHw6+rr19/7BolpbyIx7U11HzeJJ+5MIDFuIr594yX1GC8EbxYvfOmaFTLoQOIH8wbX4eUBCwEAG7dsMJXKQzQHhiX1hx4/DvtJfX3GbsOwZuh36vvPVysm8U/GZhfF85echK+eec3yOImtufnUFg9bmaOtnLmttErF+ApBRolKbJMYN1GtF1pZETWVxkM0BYYh7YcR8ctOUl+fdNQ+iI86AfFRJ6j72rRUrAAHnHYKAEVgrek3yNZYz726Qn2duGNajrby+df1/XglWWTQgcQ3Br/dD1VtW9g694HX71Nf17SuwcKh2eihulFriqbmCIrbZo4HAPzixv5ot087tNunHTZu2YDR7a6z3Fe85g857w+b1QtvTn9EfX/4hF87m2yZIzVsEieMvPlKIGrehUFk8Yof1NfJVBqJ26eo7xN3TUfirumm0nMEwSdfbwIArHjuWaCqGpUVUZx7qnnLgMjZQ7LWhfgNo9Dql0er7zd/9rGziTYDwiywSe/WMkDUACVTSfR/dn9b/Yw96Vr19dbtW3NMibXTuqMe3vue2eHW4TcD8C7qND55nFrXNNWYwvwhyzDw1b6oblGNNjNb4vvzNrk2ll7etVLJxUaMQEwKYxJriBogIkL8xtG2+jnsgD3U17FoJMeUGL9uhO35+cHB+3YA4E3U6Y41XwKoVd4whsT4yYjfPAbpNMP0xFzELz3RtbF4glzXqw9IAEgNW1lw5qJjcdqCo9Fmzzau9tuxXTZ32cKhy7H6xmxJxTkbnwAA9H65O3q/3N3VccNGjwnZIvQH3Lenajqu31WPs2sH655jpjoBbxOLxlTtZVVFVV47vX2xaMyzCgh2YWBgzNwmkXAuvnEMRtx0JVDZ0tV+G5tS6uvE7VOU7b5seqM33v9GOZZ429Vxw4ZYiD7+x9+opuNIhFDZtr3RaUX/n/LjjDFVc8mYsuW3zT83lOsAKQXrzWxBIAW2MmDr9q2o31WPQ2fu62q/++7dU3199Jw++PUJA9X3Qzsq5Rm/Om0NvjptjavjhhlRSKpuUY3/fP5uzjF+vGP7TjkCbyHat8kKhJ07dMkbh4/Ft1g0liPE8XGdCHBuCX98wS62SSScWDSCSIQQv97d5N5i6ab4TVcoGrfduqj7TuyvpLOIx49zddwwIwpJ6TRD489b847zNo1N6Ryh14gmIX3SrqYUdmXO0Qpk6bQi2Clj5I/pRIBzV/gLb9SBFNgkhnz61Sc5D/J773682ScmFk2T9bvqc4IfxGPrN63Dxi0b8oQp8TVv35hsUM9btXaleoxr3up31aOhqQENTQ1qSpT6XfXq+Vx4E4vIc8wKYkZlsSzBAJZKm9okEq9p26oy/0H+v8+adWJi0SwZiVBO4IP2eFVlFJUV0RxhSnzN24opTlpVxdAqk0eQiFTtWyRCIMq2Vd7zovG5Y2u/MzPCmJ651Y4Qp6Qmkho2SQkiCgUAkE6mcewlhf0dyt08WgjuayYKOslUUt3E9yJcQ6rXF9/Pz9PzZ2toasjZ55a/m61+0iY3icRjRKFAhaWx/qcdhueUu3m0GNzXTBR2RCd7PYd7riE16osf4/1qz1cEIBLeuyMN2a+JKgU2ics4Ladkl6vX/F/B483JPKqHkVaqmLZK1MAB+toycb+Ino9bofZeIn3YJGZwUkrJKfMefdrwWHMyjxqh9//Tij+bVgOnFf7yzy08H3+DFcyaQ0NqEiWiBBFtIKLPhH27EdEbRLQy87e9cOwGIlpFRF8S0cleTby58uKOJ5RtwLs47IAjcPPaq/GnB64HkJuWw006d+ii5ic7bFavgm15MEJzRU+DxvcXO0+vnRmBy4/KEaYw6b/mt7wm17Bw8co7/8W//t9D2FrfiFnP1WH5N0qORTEth5s0NKbw3Ksr1Pxk8RtGGbblgQjNGT0tWjGhqZAGrti5YSo5BpS+hm06gCGafdcDeJMx1gvAm5n3IKIDAQwD0CdzziQimSXTKcsu/1Z93WboAWgzVMmzs3jFBxjf7UH8beydAIDrh97kyfgvDXkXEYqYSpvBgxGaI6WQesNz0szc5i/TIdewQEncPUN9/f07/wEA1FRXYtjZteizjyIrH7qfuSAdqzz29xmqUqRY2gweiNBckak3wk1Rr2LG2Hwi6q7ZfSaAQZnXMwDMA3BdZv8sxlgDgNVEtArA4QDed2m+zZJDHumBJBRhoJDQ5Jam5fZZE3DTsBsBKDneIpsiGNH293ntFg9biQe73wMAeGfGG1izfrUr45c6VvOmWWkvtvUyP5udvnlaj7Ah17DgGfl/2Rq9RkKTW5qWulWbUNtTyWuWGDcRoAjOPvkA3baJO6Zh0Mjz0bl9K7SsCleanCCxkjfNblsvc7M56TvMMqvdO3QPxth6AGCMrSci/tOoK4BFQru1mX0SB/itueHCWt2oNYh8F9EtWTVvcB2wCRi56UoAQOTgCP742G/w1Jvz0WPCbmp0Y3PEy+/Li+CCYuOYhgEsFT6BzQC5hvmIn5obVVi7azpAEcOSVTxh77yp/1bqlu59ELB5HeLXXqBGNjZXvPq+vAguKDaOtRPDrWV0O+hA70p1V3AiGk1ES4hoSRMai3Z889qrnc6tLKic1hKV0/QTTPKKBzOWPuR4nHmD67Dlmy26wtr0bf/Af2e/hUFza9Vt4Kt9sfj4z7Du0/UAgIUXLc87T6JPOZlSw+jDZhFba9jObZv1mqjMeq7OhamVB4kpL+vvzwhQ7y1b73yMcROBhnpDYe2511ag/eHHI37LVcp28xjELzkGqFf86abfM0P3PEk+YRZwyg27GrYfiKhz5pdpZwAbMvvXAthLaNcNwDq9DhhjUwFMBYA21K7oEj6+24M2p1peNI7aaXiMm0svOfRyW31Xt6jGvLPqUL+pXtf0On/IMqSTaV3zKDeh8XJWtdO625pDOVHMrMiPl5U2MqS1GnVwdQ3r1POgghc+7OxaxxMuF+JXnKa/P2MqPapvZ1v9ptMM0+9KAE27DM2uifGTAZbOO85NaLycVdhLWflFIdMiP1ZO2kiehy2s2NWwzQFwSeb1JQBeEPYPI6IqIuoBoBeAD51NUWIVu8ltXx6wEIC+n9y8wXVo27WtoQ9dVUUVjux7lJrdXywc7yWxaEzVLPIi8GHBbGSokbBWKBWIX2WpLI3DSiqth1zDQordxLbTb58EQN9HLjFuoqJ1q2ype5wxYMv2BjWzv1g43msYY6p2kReBDwuFtGd6iXO1GP1f92sNsLfehDetR9HVmIiegOKc24GI1gK4FcCdAJ4iolEA/gfgfABgjC0noqcAfA4gCeAqxljx2hYm+Oq6jVi3VMkqP/jtfmVhRhpw+LG4ve39qGrbAq+f87RaxFxEq6UZ/iclJH3m36YZ9vvPZf8Po3pdg41bNhi20eb74kKPXtoOfqyQ1qx+Vz0WLXtPfX/qW0dhK7YatneLZCqpCpF6n18pU+ge9+v+tzIOQzjNnWFYwxJ/fwLYrphNR958ZVmYkX7ctgsv/GMaEK3AocPOUQuYa9FqaVb/8DN67GFc93jBR9/h2H7dUFlRODhX7JcLPHopO8Ti8kaas0iE0K51Np/hxTfYK0BvByJShUijz7BUMbrP/br/7YwT5v+aFIZfu22oHeuHAUFPIzDqRq0J1IQ4a3kCe15ziKFmDchq3cQIUklxtIJxKTIPLy9ljPUr1q62by1768W3TPW5e/fdTfVZKnTqeRA7/x7jhKzlTOKu6YGaEBfUrcNXc54z1KwBuVo3MYpUUhwvozn9YtI5B5pab1p27s16jDTnfrXijpN9X8NkpYMQEKSwlvhgMrr8sZ+usLZ42EpUd6jOORYGYe3lnU/lvF8+VtfFKBQYJdI1g19mT9dg2QLPxTZJ+RCksPb+p+vx1Ysv6Atrd0wDKlrkHQtaWHtl/n/z9iX+NjOAmZjDibAWBoWQVXjy32JbEJSMwCYWzZYUhk2OoW7UGrDJxp/XvMF1mDe4Dr+4sb9h2o76TfUY9HytrTnwqgvzhyyzdX4h6v+X6/M18Og+ro8RBkpSKydrieoSIt+90JOY/BISd01HYvJLhduNm4gVzz2rGwmaGDcRaNqFEdfFbc2BV11IjJ9s6/xCtG5dmb+zxpukwUFTipq5MFc6KBkJiD+8Xql/GqdWnx/wbNzlxc2zsDX2I4a3LZyF2yw0JoladAem6T/w5w2uQ0WrCgBQozq1xwHFDFoPe9GLY0+6FgB0hUGn/Hq/kTnvfxz2k+tjhAEvE+N6hRRK9OEPrv+8uxqnHGMvKCisvP7eN9ixfSfOOml/V/qLjzm9aJvEuIlAJKpGdeYdQ/GqBoU47IA9lD4M0oI4YeAh3fL2xS87yfVxwkA5mFPDRMkIbJxyE9YA4Izdhvk2FhfG9AQ1QNGIfX/vRxjWx94vU0kWrcAl+rPxY3ZTe9hJBaInAJpJPQKzLvcMpZTWIxDKTVgDgJOO8reckyqQ6Qlr4yej99CzMaC2i69zKle0lQkA5ccH3283tYfdVCBaAbCYQGj5B2QzS5wrCRFcONO+N0rNMW9wHVbdttBVYa2mdY3lc7Rm1OVj1+H9c74wPB5W9Aq5a4u8F0rtUUiQspO3zUlRerO4mTiXiIZkCrCvIqLrC7Q7jIhSRHSepclKQo0Y3andZxRgcMBZZ7kqrCVT1u33WjNq4m8zkZgw1fB4mNFWJhCLvIt/jYQvI+HHbt42O0XpgxLAvFi/Sk5gEx/c4mtJLmcsGQQgX2gzEtZ43rT4EYoJwC1/wa3braf20JpR+zzQBf2f3d/wuCQcMCilqcxsxcgUXJ8I4BQABwK4MFOYXa/dXQBec/dqvCHnwT1has57SZbpjynVwdT8aQK6wlomb1r/XypJd90yzcei1h+RWjNq/E/DEb9xtOFxSXjgiXOd+rB5tX6VnEl0yOtHqPm9xId4KTP0N+fisgfPd9U0+mK/eRiE2px9hTRr2JF7vNR8pyT2cNVPzl3H+sMBrGKM/RcAiGgWlMLsn2vajQUwG8Bhbg3sJRdff6n6WnyIlzLf/ViP5V9uctU0OuKiIwEcmbffbOqOMJu1JO7hvo+caxo5T9avktOw2dHYhJ05j8/2zI+tc4fC5oFiZlK/mTe4Dge/vl/Q0/CNoCOfjYQ1u/NiaWZqg5LEdomwaaWXrgC+Fd7nFWEnoq4AzgbgX1p6h9jR2ISdrrtXe+bH1tCYQkOjsROlGwEGbpIYNxGJ6fODnoZvBB1kVKhslg8UWsM8Wb9KTsMGAB3bdcLGLRvKqwajRzxR+0rBhLhAeIQ1zu39bsMZ8C8QI0jCqskUAyKsYCHH2qYiSSfNFGH/B4DrGGOpUtKoNDalUFkRVT+rcqnD6AWP3ak8y8xq1sJAt97+BmEESVj/34kBEdbPNd200BrmyfpVkgIbL7kkhbXiGAljQQhp1S2qccUXF+De7gkASvDAfx6ehFMuuxLrn1HKNe5xVj8M65JfXkbiP7YCDtyLEjVThL0fgFmZxa4DgFOJKMkYe96tSXgBL7kkBbXiFBLG/BbU0mmGJ579CABw0XmHAlACCFr2ORI7ly9Ch6NOAAAMPnJvX+cl0SdgYdKT9av89POSUHL/j3/FywMW4oxLf6dq9+65/Xe4p/sj+MvmsdjjrH7Y4yzlxwovQu8GMjDFX9KMmdpMsBhALyLqQUSVAIZBKcyuwhjrwRjrzhjrDuAZAFeGXViTlCZPv/Qppt8+CQ2fL0LD54tU7V7kF7/CsLNrceY1I7Fp4RsAskXo3UAGpfiPS5UOPFm/SlZgq25RnRcBKSnMvMF1qk8br3TgF9fs/lcAQFXbFqp278V+85BMJbHgw3cwaG4tBs2txeC3++HJKXe6Nm65BKaUAowxpJJpU5uJvpIAroYSPbUCwFOZwuxXEFF+Aq4SI51muhGQksIkxk1Ufdr8/OzOP/2XyotoBRCtULV7Iy46EkSE3dsqJbBG3nwlKvc/wrVxyyUwpWQwGSFaTF7zav0qWYFNmkOtM2huLdZvytXKikJboVJWdpk/ZFmOcFhMgEqmkrjgivyUNVaFS56mROIvFoIOivfF2CuMsd6MsX0ZYxMy+6YwxvKcdBljIxhjz7h8OZ4hzaH2iN9yFaoqo+p7UWgrVsrKDonxk3ME6/iNowsKUUSExi8+yO/HonDJ05RIShcv1q+SENjeP+cL9aF/2Ve5ueU6d+iCc+tOVN+HuRB4mLj3zt/j3jt/n+PL9smSVerrju3cqW137qLBABRhsaptC1Pn6PnXWfW5M6rkIPEQWfxdl8SEqepD/99PLs473tCYwvTHF2Xbh7gQeJiI7Htorh/btk3qy8Yms+U5CjP8ekU4i99ylaJdM4Geb51Vfzu9Kg4S71HysIW3+HtJBB30f3Z/zBtch9MWHI2He2eF0EFza7Ee6zC7Niuk9XlAliQxw4v95uXtO7hfT7X+6LSV9+N3ff+CNetXAwBmLH0Ilxx6uen+3z/nCzRs25WzT5onyx8Xgw7KhviNo5EYNxEjbroyR7MmPsRH/Cabcyz+p+G+zq9UUXK1CbTtoL58Z8laDDikK1pWKY+495atBwAc1bezqb4TE6YCqaacfdI82TwIaeArgJBr2LifGjeHmTGD1o1a4+2kyhgak40KHNrxQlVYA2BJWAOAY144CABw2oKj1ULzpUDdqDV55lTpK2kOBnMBByaDDkoerZ+aWTNo4q7pHs6qfBGLxp/Yfx9VWAMUQc2ssAYAI/9yGQBgxE1XApFokdbhInHX9ByTqvSTLB9sC2xEtB8R1QnbNiL6PRH9lYi+E/afaneMhqYGAIombcXTr5g6p3Zad7vDhZ6r77oW84csy6ml6ZeAKgot1S2q8cDr9+m2Wzr8a8wbXIe5xy0BoAjZpWSerJ3WvaTmGyqYuz5sXuP1GsZ/qcdvuQptDz3W9Hnx60bYGa4k+GLtlpxamn4Jp9oAj3SaYfGKH/Lb3ZlAYtxEPDpeifSMRKjkzJPx60aU3JzDRFmaRBljXwJK7aNMPazvADwHYCSA+xhjf3MyMa1WY0z7vzjprix48Lr7MBC5dTT9ElBFH7L6XfUYe9K1uu0GPtMXLw9YGLpkvE6o7lCd837p8K9x6Mx9A5pNeGGAqQjQsODlGqbVapx32kG251lO7N+tHfYXamn6JZxqfcgiEcJhB+yR127En0di+u2TQpeM1xEVub7DiTsTiF8fD2gy4SfMJlG3fNhOAPA1Y+wbJ5Ln/T/+FQBwx9sPAlMUc5qMBg0nf3rgevxtbH76jVisJNwiLXHYrF457/c+sLifZKEqAWKFjmLVBPhxXipKbNuzWy+sWruy6Fx8w93EuX7jyhr29Euf4ufvvweAPJ81SbhY/s1P6LNP+5x95Wiuj9+gSURe3c7UeUaVAtJppt7XZqoJ8DbatjsakmhVVX7PCy9xy4dtGIAnhPdXE9EyIkoQUXu9E4hoNK/B1YRGAEqurmt2/yu+P0+J9qnfVV92UZ92azTGojHMH7IMvV/ujt4vd9dt46f/Xsuzt+vu37p9a1lp1/To+JeWRdsUEsLEHyHFqgnw48lUMq9tqIQ1AIC5CNGQRok6WsN2btsMQMnXFb9UiVqPRKjsIj6d1GhkjClpMhJv6x7323fv04++ydsXi0bKS7umQ3zsOabaGQli4o8QMz9ueBtt21AKaxRuk6hjgS2TxXcogKczuyYD2BeKqWE9gHv1zmOMTWWM9WOM9atAZd5xboY6/OFeecdKme6de9g6b+5xS5BOpvHVaWvw1WlrdNv46b/36on56Qkk5YWdHxeMMVNbmHBjDWvZdrfcgxkz1Ig/XOTJnINiZ4FC7MV4dPwkgKURjx+ne9xv371TTjrQ1/Ek/hK2dcYN3NCwnQLgI8bYDwDAGPuBMZZijKUBPAzgcDudcjNUuZlEzWpFxIcl9+cLi+Zq/pBluHefh1HTuiboqUhMoid82dX2GsFKNw+b62sYN0OVm0nUrFZEK5iHrVB7YvxkzP77NCRTpeNz2dzRCmBeCGRKHjbnlQ68wg2B7UIIpgQiEmOnzwbwmZPOn/ryUQDllxD3mtUXF3xgcvPXvMF1iMQihsKaaB4Vo0edMvxPWb+HZZd/i8ppWTNgOpnG1meXY+v2ra6NJ/EWPdOrWXOsaRhcK03lM56tYfM/Wgug/JLhznx6ScEHpmg2SoybCJC+qVE0j4rRo05Z/cPP2X7vnoHElJeRmPJytgFLY8+BJyMWDXVmK4mA1gxZzCxpz2xpzhxakiZRImoF4EQAzwq77yaiT4loGYDjAOiHExahZzdFw9b1/5Qft0EkxK1ukY0OdKJN0mrLDuxxEGJk7GzOqwxwzdrAV/vqtps3uA6PPJCtJ/vUcY/anqOWmX+bpr7u+9BeaBy1U30/aG4tzmh1oWtjScoDBnMpPcIUmODVGrajQfm/veoVpVxSUMlwuTbTqSaJC2eJcRPx886mgg8tscqAqlkTIkNzjn33hfr+qJHurSk99mijvo7/+RLErzgN8StOy+675SqceuwvXBtPIvEDRzYRxtgOALtr9v3W0Ywy8KStQebEEs2xTrRJomA2aG4tPsdn+Ly78Y/2jVs24OwRF2De2jpDzRoX5o5/4xAkofT/4HX6udEk5YNRVGmxaFM3xypEmIQxM3i1hrXM1LsMOh8WN8k61SRx4YxryS4671DDtpUVUazdlF07dTVrGUFupPD57N+tnaM5SsKPUVSpmWhTt8YqRpjTeoRWH+z2w6dUqG5RjQN7HIQL/u9kXWFNTJw7aG6tL5/Tqeecqb6ORWNYdvm3no8psUYo/r+Urg+b6wRlMgmadJrh551NWLN2K9as3aovrGVMn/FbrvLlc1q/eUfOe8YYEnfP8HxciXma6/8Xq4QwrrZ5cNhbB2HS1OcAAP2OPQA0JonFw1aiZp92eHnpQnzeJ18DN29wHTLKNM8DEETtSpdfdlUNRslUEn0f2svTsSXG+CmY2RmrHCOzJPkkZrwLrPtKedNuD8THnI7EHdOAqlbA9s2GwQV+BR+I2pXtu3JrghIR4n++xNPxJfr4KZjZHSvMwqMU2Hyme+ce2PjTBnx+6moctiuTsmSW8mDkkbG90THvPL8jRcWH9SO3TfJlTElpw4BmoT1rzuxsSKKqIooRvx2ASOSYnGN5CVo1+BkpKj50e3WR0eyS8kAKbD4jFlQ3y6zlCTx5z2t4bvqT6j4vfJYk7hIWvzLfYCyMEaASF2lpM9npgrp1OOnKEeiWya/phc+SxH3c/p5C/71TuDVsofVhCwNu56kS6dmtFybP+2fRMeYNrsO40/+O56Y/qUaPAiHxWZIUxIvvKMzfOyux4u/NAa9M1Dsakvjw8x9M9Z8YNxEH79cR3TpUqxGkYX4oSrK4/T2F/XsPex42qWErgJcPx1VrV6Lu3eJJdLsc2hXfLFC0chu3bPBsPhKJHrFoDLCQ3D6dlhq2MOHVA7JVVQztaloUbwgArXdDdUYzV1kR9WQ+Eoke5eZTW5Iatqe+fNRT7ZcfxKIxzNz7Ubx14kcF2239ZkvZVXsoNx5f9oju/mL3KM/zd2Tfo9T2Na1rCpYvq2ldk9PvsQOONzWWOCe+uQ4r6VqivjH/o7Ul/yBhjGHBS+/h0dunFG/csKPsKj6UEws/MU5KX+g+5f+Pt2xvUNsmU2kkU2nsbNBXdiRT6Zw+N23bVXQcvTl5+f8nzBq2khLYYtEYnl41A5VtK0NtGjJDMpVE/SXbDZPicngggpu8uHlW3r6Z2ybmvL9t5njXxy1XftP3UgDFhSbtcS6IL1r2nuqbtnX71oJ+jmI+wFg0hncWvAXAujbYbHsr/TJIk2ghGGN4Z+laNDamQ28aKgYRIX75KboJcbUUC0aww+vv5RZuf/71L3Lef/L1JtfHLFeOPjiblL6QIKQ9xoXwdq2rVN+0WDSCWDRi6Oso5gNkjKFDW0VLa/X/g9n29qoSkMnNf0pKYEumkri09nc4q3N5FVQuBJvsvibk4aufzttXk8zmDq0btQYnPXu++n7quw+6PodyRCvcFHsPZIU4K4IRb2vnR0sylTR1ni0NHANSTSlTW3OEiHB0bRcMPnLvoKfiK4nJL7neZ5/9OuS8b9U6Wzovcdd0LJ2VLVzxwfLvXR+/XClU/sko4a3RsWJjeJl2w4kGTmrYXGLe4Doc9saBmLrlrpz959adGNCM3GX+kGV55qol76xwfZw5j8/O23fGbsPU17XTuqMho6oGgNHHXO36HCQKVoUuv1wBbOVgA0M6nTa1NUcS4ybiX48vwuxXlucdm/74ogBm5C6J8ZP1zVVbfnB9rK67V+e8P+mofdTX8etGAKls7rUj+uzp+vgSBatCl1+uAE402FJgc4Ge3Xrhkftvwtz+72N0u+tyjvXr/CtPxvTbT+6UC4/CqL+OwQNv/t2wjWf+Rxpu2naN52NIFMSatYC9urWh8OmUlQ4M2dGQRLRXP8R/exTOPTW/3F7L6pY6ZznHV1+5Tt2xct02LF5RODjKjzmd+Xv3zbASffT+P1utXRsun87wmkRDsMoXZunwr/Hz9z+jbp81OK7iYN02R18yHANR2BesFGiI78JDqQfUwvfVLapx5eizsXiWUvXgj2viSDJF83F/j395OpcFH77jaf+SLNqgkvqdxkEmRpqvsPh0Nlf/NCMSdyaAxp1AVTUquhvXRa5f9h5wpv76ViqMHD0ERIQdHbI/QNJpBnTprb5/7JmlYIxh+Pn9PJ3L7m1NRrBKHKMXUBI1CDIx0nyFxqczQO2ZGUKvYTt05r4AFDPdw72fMWzH62sC7plIV16/AQPf7YeB7/ZDdYtqdO7QpfhJDuAP3VVrs+k+Fh+fLVH1n+M+xP09/uW5sCaxxqzlCcyun277/FBox1yAMWkS1RK/Pq78vW4EfnvBYQXb8hqbgHsm0kf/8RQS/1qIdJqhodFb30H+0G2lcTiPX5KtiDD0tL6eC2sSayyoW4cX31qJF98qnmZKj3Bpx8qb0D4pThh8Mm7GXcUbAnmRlr/6YhBm4w3Hc+gxYTek31B+LX72x2/RY4K/6TU6tu+EeiFicNK+MzFoda36vme3XlizfnVotCvNlWF94o7ONxOcwM2mYU/xIjVsChu27sRL9ydMt9dGW9Ye3tOVecSvvQA/btuF6fc/hfi1F7jSpxUamlI5EYNP3zNVLU21oyGJlpXR8GhXmikDap0pIvS+P7196TQLfXoXnjg3rIRWw/bm3NcAKLUzN075xNK5A6efkfN+6fCvcf4V9iJLG37ehYafd2HXjkbMG1yHy8ePxbCxF9vqyyra9A7aOqKr1q6UwpoHLB9rnJcoKOp31YdfWGNAMpkytZU7nWoUn7T4LVdh39POtHz+R09kI7kTdybwzYbttucSi0aAdAqJcRPx1XdbsWbDz7b7soo2vYNYR7RVVUwKax6R+NvMoKeQR9iFNQVSU4EU24IgVBq2G769AgBw8sgrcvaf3/MSS/1oNW6HztwXT+MxW3NavOIDAMDxGw/BE7Wv4MJ3RvlWgF1LqOtIlhF9HrD/i9PKd8TbmjmnkIbN6Hz/75fmm2ON88SzH2PnZ+/l7Dv20G6W+xE1btysapea6kpcNPYCPHbnFCyYNtOX4utGhL6WZJkQ/9Nw2+ea/Y54O7PtjTRshc4P4n4J893pSMNGRGuI6FMiqiOiJZl9uxHRG0S0MvO3vdn+7un+CO7prmSNHzS3Fu32aedkeipPbtfPRG+F9ZvWofdAxZ9u3uA6VLeoxuJhis2fZ5v3Gims+cM//nMvTj79dJx8+unqvoVD81Mx6OFVPrVCGrbQBCIwlJwPm9tr2LCzawFkNElV1YUbm2TOm1857qOqMgq02wOAkl4knWZI3DENQDbbvB9IYc0flnyxAd//tEN9nzBTkSKDlaS0VtobadgKnR/I/RLeIFGQswRztAZAP8bYJmHf3QA2M8buJKLrAbRnjF1n1AcAxDq2ZKna3XP2BaXFKgaPWq1q2wLHvHCQFKIkZc88vLyUMVbUU/wXXXuz/zfmAVN9XnjzEFN9eo1ba1ik9Z6sqtd56vsgtViFUKNWoxUY+ZfLpAAlaRZMOudAU+tN270PYEf8+VFTfc4d29/3NcwLH7YzAczIvJ4B4KxiJ1BEmcagubWuata84NCZ+6L3wH3RsG0X5h63BNUtqj2pRiAJF2NWDct5H0Y/t+BhYOm0qS3kWF7D+FLqpmbNC+LXxxVNW6oJj46fpGjaPKhEIAkXM59akvM+jD5uYaGcE+cyAK8T0VIiGp3ZtwdjbD0AZP520juRiEYT0RIiWpKqb8L7jz2uHqud1t3htLylyy1tcNRQJVnvywMWokIWAih7JvfMrb+68i9vBzQTb9BLLWI53UhpJs51ZQ1jTfWoPngAgEym/RAT/915QIe9AADTb5+EkVecFvCMJF4z/Ne5iqC9BwwMaCbeoGcpLMd0I04FtqMZY4cAOAXAVURk+i5gjE1ljPVjjPVr31CDj279Tj1mJ9O731T+jqHNnm0AAHOPW1KktaTcaPfb/CSoWgGHV6UwEnzMCkRiO/662F+r6Jn2rZr7WWnWEnVlDevQrRuO7K/4uFrN8h4E8SuHApVKJOuj4ycFPBuJ3/zv9Rd192uFHN1SYwZtC8Hbmv1rFbOpRYr3ky0YH8YoUUcCG2NsXebvBgDPATgcwA9E1BkAMn8L1ykB0JhswtNTslGcW7dvdTIt3zh05r5os2eb0PrbSexTrASY3neuFXCqKqrUYut6fRntL9SvNlBB/KsXFepvUl5z2rUwadjcWsOICPt0ag0gk0ajBIhfHwcqW4bW307ijELCj9F3rhVEGIMaCarX1qyApQ1Q0PtrFBHqt6YsxDEH9gU2Iqomojb8NYCTAHwGYA4AnofjEgAvOJ1kmOGVGCTlxQ3fXIEbvslNL1M5LVvvsXpGa8NzuZAkRnXajeS0InAFHQDDGFz1YSOiIUT0JRGtyjj/a49fRETLMtt7RGSptpNcw5ynDJGElyefz81fmpjycvb1Q/8peC4Xknhkp5OSUlaFuqBxy4fNi/XLyU/BPQAsIKJPAHwI4GXG2KsA7gRwIhGtBHBi5n2zYd7gOsSiMTUQwa+UHxL36NmtF068dDROvHS0um/xsJX4M/1ZfT/3BWMzuJuCk9O+/Bbi3NKwEVEUwEQopsoDAVxIRAdqmq0GcCxjrC+A8QCmWpyuXMM0JMZNRGLcRDCWDUbwM+2HxB12NCSx47P31feJO6bhsJP7Zxts+rbg+W4KT077CosgZwWv1i/b9hLG2H8B5EmEjLEfAZxgt99SZunwr4HvFZ+20649GvVI4p0FbwU9LYlFVq1diVWzFwIA4q8qCUwPm9ULt8+6TW3T/9n9A5lbuGFuRoAeDmBVZp0BEc2CEr35uToaY2KG2kUALGWolWtYLok7s6W0Hh0/CSNuuhIA0EEWUi85WlXFcMBZZ6nv4zeMyjkev3E0JPq4JCB6sn6VhrOFQw7uWevLONynDVCiR2PRmJpcV1JarFu9GetWbw56GiWFxaCDDjzCMrNpnyBdAYhqgLWZfUaMAlDYzlOibNvR6Ms43KeNByNMv32SomnLJNiVlBYtW1YEPYWSxIJJtNAa5sn6FYoEYpF9IsA33vX/yao67zrXcOjMfdH4T8J7cz7Gu2d+hsHPH4p6hLsGpCSfm4bdGPQUAkEMXIhFY6iqqALMWsQyaT1MsqlI0km9n7m6nRPRcVAWvAFmB3ebTRu9C5Rq26rSs761cJ+2xKQ5wKZv8ej/exgjbrjMt/El7lHbs0PQU/AdMXBBiXC1dr7FgIJCa5gn61coNGzpb7wNg18+dh3WjfOv4HHl75iaXPflAYppbdjYi30rGi9xzoDDj8WAw48NehqBkkwl0dDUYOEMBpZKm9pMsBbAXsL7bgDyshUTUV8AjwA4M2PKDIQOHb1LRZT420wk/vmMZ/3rEb9yqJpgd/rtStoPP4vGS5zzYzP3PVTSb1hPcutSWg9P1q9QCGwcNjnrrH/prVe61m+fB7rgpeNmFG/oIl1uaaOaR+cNrkNNpxrUdAp/fjmJwoIP38GCD9/xZSwzOdSqW1SrBeCt9msFbZCClaAFxgCWYqY2EywG0IuIehBRJYBhUKI3VYhobwDPAvgtY8x5wU0X4I76K9e5p3GL/2k4+pzgv/Iw/rvzVPNoYtxENDaFP7+cJMvuPvkems2hZiedj9WUHlpBKsCcaZ6sX6ES2GhMEjRGeUA8cpu1ZI6dO3RRXy8fuw5Lh3+tBAFk6H3bABx2wBGG5/OHG38oDjj8WAz9zbmW5qCFl7ECgAvfGYUL3xmVN1dJ+VEshxuH32tmisA3NDVY1HYFVPw9lTa1Fe2KsSSAqwG8BmAFgKcYY8uJ6Aoi4vlWbgGwO4BJJBRvD5L4mNMBAL26WP9x1tCYTSic+NvMnCCA5bNnY2t9YV82/nDjD8Yft+3Cdz86c8dQS1kBWDBtZt48JeVJoYS5HH6faXOqGWFHbgpC2HIjD5tX65ej4u9u0YbasX7BuZ94TucOXfDWbz8CALRoVYmD7t0rJ0eXWWpa15RMUuHmDBfW9JLYav3D6nfVo7pFdd79oG1bSPjSHi/W3ipmi7/v3aEH+9MZ40z1ec30i0NR/N0tOvU8iJ1/z9NBT8MzGhpTeOyBJ4FIFCOu+bWan8sqyVS6ZBILN2e4XCAKTHr+YZEIIZ1muveDtr2R8KU9VqitXcwWf2/X/QB27I3/MtXnnNGH+76GhSLowG2O7HsUliz/EIB1LYPbDzsAWL9pHS5+9UwAQFWbFoi8UY/5x1hXBkhhrTRIppKoblFteB/xe6y6ZTVi0Rgamhos33eikGfGjOnFfa2FMYZUY7DJe8uFLdsbUFNdaevB5cUDr6oyinOvOA/JVBrTZ76H+MVH2+pHCmulAZEiiOndRvz+SqXTIKZozqzec6KQp2fG1MOL+zqf4MpOmaFs/vf0frm7+nrRsvdQ3bIa1S2t+fwASjkhI5yYSBev+ACLV3yABR++gxmvvm6rxqRVHyaJPzzw+n14edtTOd9b/a76vO8xr5xUUvnbrZPim1rTugaxaAxH9j0qrxwVv1/4PcDNo9UtFKGvWP1dLkRyPClZ5a4PW7MjkXhbfd2udRVSNkt4GRlNnJpHa6orFb+oH9farjEZprJkEoXFK37AqwvW4NUFa3K+u0gkt/SU1vQZyfzdlTGRJ1Np9b7Ysr0hR/DhpaeyZlRlfzqttDdTf5cLkUB5FnY3Q9lo2BovzZWKu3ZUHoJmtVKdO3TB+k3rCpoqN6/aZH+CAj0m7IYkjDURRpoQO2ZUifeMPela3cLvxTRa/Ptcu+FbpeZoUqk7+vX/VhU9h/fNBUOj+1ych5lSWU5Ju5c4t9lx4UXH5Lzf2ZhCm5bmf1M3NKZQVRk1NFe2qIg6mh8nfu0FBY8X0lDYNaVKvOOwA/bQFYCKabT4d9miMgoiQoSy332rKn3RQqtV40JhIc2rOI9ipbKcQrDna+cXZaNhW7N+dc77z1d/hs9Xf2b6fDPaK7+iBiWlxcmnn45unfbK0c5yjVYhTaqRVqx75x4AcrVgvIi8Fj1BUUQ8x+ged0vb5nKUaLOjpeYh18Zi4tNUEa2DX1GDktLi+592YFdjCrsaUznaWa7RMtJmFdKK7cxo3fQ0dCLavo2KzHP0NLRua9vI5L8gKBuB7bADjsipaHDDt1fghm+vMGzf++XumLltovp+1dr8igTaB640SUqM+OnnzXka0IamBkNBC8hqf7UasMUrPsjZL6JNAcL7Nwp0EDHS0LqmbWOu5mFrdmytb8ypavDEsx8btk0k3sbzr3+Rs09PqyE+cKU5UmJERSyCilgkTwNKZKzN4loxrfYLUMzn4jGONgUIT7uhF+Sgh56GNsw+Z25TNgLb4hUf5FQ0uGOvKbhjrymG7ffbozumL5+J59c/Ztimfld9zgNXmiRLn1fq3Y/ku/Dxwarw5Za2yqgfbTQop5DvpZ+kU8zUJsmnproyp6rBhef8yrBtpKolNv+wGXMX/a9gn5EI6T5QJaXJf95dXbyRRd7/8BvEohFVAHNLY1VMWyYeD5NLmoXSVL5TNgKbFr1M9Tz/WecOXTBn8QL8JnIxzup8Ud65Pbv10s2V5omjtsRXhra90PU+r+lyq/paGyxQCKf3k9i/Gz8mnM6HMYZUU9LUJimONlM9z3/W0JhCetsmHFDbA4OP3Fv33B0Nybx8ac3VUbucGDKgu+t9njF4v5z32mABI5zeT2LfbvyYcOv+dqnSgSeUrQSi52+2ftM6HDvgeHR4MIp5X+xE/IgxuufqmUeBABKRSkoCUbtmRkgzSpTLj4lmT/5e21arUdPmbLOaxsOVtB8ZHzaJO2h9zqoqo9i0bRdeeqkOiFWi/y87G56rZx5tTqYjiXlEh/9igQbicT1zp5h3TdyXn2sNEEtravO1mTWRmp23WSzWEvWdstWwcUStQSwawzsL3sIdTzyOH4f9lNe2kHlUUh7UPNbGM19EMwKPnpDG71EuhCVTSVRVVOUIdrwdb1vMH82q8OXWjxHpw+YuuWYjhg5tWyD9/WrELztJt30xE6mktHn0kTc89UUsXq0gX0jLCme57RjLTQUitjXSqIntrQhgrv0YMWkOlSZRF9jzmQ7q68ppLVE5rSWaHswe5w/E3nd1zDs38cFkPLjmwZxABEl5cdhbB6EiFnXdF/G2meNtn2vkI6mdI2+nNbkWCmrwG55nycwm0SfxyBvZ11NexqNTXlbf8wdi/A/6Zv33P12P/61cmxeMICkPEjPeBYhc90X85Gv76apEwUo7r/wAhlwhLGjzYiliW2Ajor2I6G0iWkFEy4nomsz+vxLRd5naWHVEdKp70y3M9+fl33hPHK0EHsSiMaz8y/e65+0+qz1GfjQWLSuqMO7dezydoyQ4Fh//me494pRbh9/sep+lSilp2MK4hsUvPTHn/X5HHQJAEYanPzhb95zEw68j8fDrWLFgKSKRCDavXe/5PCX+E7/kmLz7ww0O3rdD8UbNiHL1YUsC+CNj7CMiagNgKRHxn4f3Mcb+5nx69kmPblIm8sY/ACjaiC63tMlr99V1G3HD/ifivcHL8eKmeT7OUOI3sWgMD7z5d7w16wOcefdRGN72qqCnZBq9tB1+lJuyRJohVVqFwUO9ho28/FQ8+q8FOPrgLiAixH93Xl6bxN+fQGRPJW/fhecdhqpKd5LjSsIHYwyLV2zAHh2q8XHdWpx10v5BT8kSftQMdYPwzSiLbYGNMbYewPrM65+JaAWArm5NzCn8QdZv74OwGMYJdHvf1RGbFqawftM6v6YmCYhkKomxJ/wByVQSjx34OCqLn+KIQgKVVWFLr62VgvB+wACwEqp0EPY1jIhQ2bJlwTbxP1yITZloUimslTdEhMMO6AQiwpurvgJ8ENiMhCo7wpbZmqG8/2JtvCKEMqSKKz5sRNQdwK8AfJDZdTURLSOiBBG1NzhnNBEtIaIlTWjUa+KIeYPrMG9wHS68+i+6x9eN+1lNjPvOgrd003hoOeyAI9yepsRnuBBzX9c7fBvL6rFCmE29Yad/x2lrWOnmYXO6hu3cttn1OSXGTUTjl4v1j/3zGTUpboe2LdChbYu8NB5atta7v85K/IULML86sb+v45ndXwyzqTfsmB3dS1tDJjf/cSywEVFrALMB/J4xtg3AZAD7AqiF8uv1Xr3zGGNTGWP9GGP9KlzWdSwculx9PfDVvrpt9r6tvZoYF4ApDds9XR/C1HcfxLpxP1ue0/lXXITzr8jP+SYJhuvOvT7oKdjCS62Z875Ls9KBG2tYy7a7uTqnxO2K7238Zv3UQyPHnpuTFBcormGbfe/D+GD590j88xlbc/pmw3Zb50nc51e98gPnSgEvNWaupfUo1yhRIqqAstA9xhh7FgAYYz8wxlKMsTSAhwEc7nya5mn8p/JJDppbi0Fzaw3b2Xk4DZpbi9HHXI2fNioCm1gKy4gntz8CALh++zhcv32c5TEl3iBqkw7uWYvHlz2ivn//nC9w+fixuHz82CCmVrqUYC3RMK5hiUlzAADxW4x9LO08nOK3XIUj+uwJNComVLEMViHmvPkVAODNxFOWx5R4g1abtG1HIxZ+klU6JCZMxVffbfV7WmVBePVrzqJECcA0ACsYY38X9ovZHM8GCjiQucy8wXVY9MonOHpOH93j84csw5hVw3DYWwfh8WWPgE02LsxdiD4PKOZTsRSWERe0vhQAcOjMfXHozH0tjyVxj/lDlmHOxifQsV0nvHPaJ+r+ufH3MeaoP6rv+z+7Px66+QE8dPMDQUyzZGEA0um0qS0MhHENS4ybCPz4HeI35ddBToyfjJlPLUFixrtY+Mk6JCa/ZMsMFP/TcADIKYNViKEn9FbOuz5ueSyJeyTGT8Yb73+DxqYUHp3wUM6xZ6a9hCMO2lN9H79xNHp3rfF7iqWPWWmtBDVsRwP4LYDjNeHvdxPRp0S0DMBxAK51Y6LFmDe4DoCxCXTe4Dqkk2lM7jkLH5/4BUYsHouqygpV08bLWPGC8Qf2OAg1rWtykpUuHf61x1ch8ZKBr/bF0I4XYuOWDTlCfce/tFSrFUjswxjLyQ1XaAsJoVrDEuOUHJB6ZtDEuIkAS2P4r/th5MUD8OX7HwMaPx9eyooXjf95ZxOSqbQq1CXuTHh9CRIPid88Bif23weVFdE8gT4+9pycigUSuxAiJrcgcBIlugD6cuYr9qejz8Khyw21ZkBWWDMygS4cuhwVqFD7SKaSwKjchwYvZcULxn++WvlRfWTfowAAi5a9hxbtC0dsFWL+kGXK3xEv4uP952F27RtFzpAEBZscA40JjVBRQjCk06WT1sOvNSxx+xRdjVlOGy6s6ZhBE7dPASLZhzQRIX7FaXnteCkrXjS+TcsKAMCW7Q1o17oKiDnzFU6Mn4xDLjwfdR+uwojfHOmoL4m3JCa/hPiY04OeRmlS7lGiXjP9b7caHps3uA6RWMRQWJs/ZBkOG3xQQYGvEIuWvYdFy94DkDWF2mHgq30x8NW+uGnYjZ4Ka3Wj1gAAlo8NT5qSm9deDQDo3rlHION3bNcJALDiWv3EyVyYBpSKF+++8JEv8yo3GIAUS5vamhMVvQ4teDwxbiJAEX1hbfxkYLcuRQW+QrRrrVR44aZQu8RvHoPanh08FdYSd01X/v5tpmdjWGXWc3UAgJ0NwfyIa2xSfgQl7n3csE1i/GT19fufrgc2fev5vMqVEFtES0Nge7h34agmPTMoT+uRTqZR+Tt9Pw+xlBXHTCBBmKmd1h2AM+HSbcZ3U+qDrVm/OpDxN27ZAAD44VN9gU28f+JHjDE0qxvxXP2/7U+unGCl5cPmF7+94LCibQqZQeNXDtU9RyxjJWI2mCCMxK8bofx1KFy6ybCzawEALascpr2xSWVFJvq3Pr/+NUe8f/r/srNhdLEeL729yvbcyg0lSjS8lQ5KQmArhJ5mjZtIjY5z9MoUVWbqjXbv3AO9X+6ec8xMrjYziBodiX8UuhecMn3bPzzru3RgUmCzga5mrYCJVD3PoEwRj0PY2ZBEIvG2ur9YnjaziNociX8Uuhec8txrK/Dcays861/iDsH8ZLDIvMF1uOXWSzHuNiX1wlFDf4U7dhunW8Nx/pBlQNL+w3nxCiVv5pr1qwGNm4ib1RDmD1mGhTNmYsn6jz0xkZ5bd2Lo/OQO7HGQ6htYbvzf3v8vkHGrW1S7XszeLgxAupmZO82QGDcR6F4LrKlTdnTYC4eefJRhDUcuENl9QNdUK75qLatiiMePU/e7WQkhMX4yqvsehZ31Oz0xkU5/fFHo/OR+3tmk+gWWG6cd3zuQcdNp5noxe6eEudJBSQhsg+bWYv7cJRiEWmXHXIZbkS+szRtc50hY8wtucrthj7s9G+O9wcuLN/KZSfvOxJNv3InJPWcFMn7ig8mIH2HeVGAFbnb1m7AIa0AmSjTZFPQ0QkdW8Dq6aFszmrWgsWJus8uF5xU3I/vN0/dMReX+SrWb4b/u5/v473+6Hv1/2bl4QxuoZlefCZuwBgRTDsssJSGwcS776jxcdOVNhgJZ2AU1P3miNhvoxj8XURvz/jlfAFByjvnJBVdcnxW8fcYrYS2N0omM9BYmNWwF+PeTi9H05YeFzZwhFtT85LE7p+S855+LqJFJTJiK+I2jfZ1X4xe8cpn/AptXwpp7JZ3Kg/CKayUmsD3c+xndh/1tM8frmkdLDTcLducIaVCENFEb47egxuck+heWC43pBt39XhRgF/vk+QG9zGsWi8ZgWh7NBB1I9FGCD/Q1R598vcnQRFoq2CkIboQooImIGhm/hbX4LVepGtBygqWNBTY3v1Ntf34UeC83YbTkgw4A4IbNtwQ9BdtM3XIXAOD59Y/h8c8fdr1/UUjj6S3M4LgQuAHlqAVNMX2ByY4gVd2iumgbUVALURLakqt0ECaWvvZe0FOwxexXFNeLuYv+h3c//s71/iMRyhHSeIqLYnj1oI7fclXZaUEL/X+0Kkyl0yxPyNYiCmqhMz+arCNakrVEg4Zra96b83GwE3HAUyueR03rGjTVJ/Hr/Ua61u9lX52Xt4/7WdW0rsH753yhmkUB4JX6p3Pazj1uiWtzKXdaRosLWWYp5pPmt5BmbSyGFEuZ2iSKv5qqsSnRvFlbf9yKZCqNiooIBh7SzbV+//3kYt393NcqmUojMWGquv8/7+amDHp0/CTX5lLuRGPu+a9pBWwtfgtpVsdTcqyZ+xcEJWUS1fLYpNuL5mgLK5d9dR52Nu3CpdfcjkHba3HhfqNc67u6RbXi62fgK7Z1+9Y8k+ip1ed7YsJrDpxbPSLoKYQCJehA3j9mqdjP15ryrvLvJxcrmpnVHyMWPQoDD+nqWt/pNEPTlx/CyHwMALFoJMckesoxPVw33zUXzji+V9BTCBVhvoVKWmDjwtott16K+XNLSyN01V1/w3mvnYhBa2sBmNdk8AoGeolxJ//0/3DA+ae6Mr8uh3YF5iqvY9EYbvjmCjUBrkSih0zrYY2chLrdawObhx2a1izHsN8PR6sqJWrSiqCU+NtM3cS4z7z8GbYtfcedCbbeTX3JGMOTz38CIJsEVyIxIswCW0mbRDk8PxugaK4WDg1fSounV83A06tmqBUYtnyzBavWrrTcT58HuhhWMRjT/i8AFD+xqkxdQSuIQmPvuzrm7G9Owtrw5bIGny2YTJxrG56jDYr2KnH7FOO2AfHO0rVZM25DPVrZzPxvVMXgvNMOUo7fchUQtZ7vTBQa43+4MGf/sLNrm5WwNuOJD4o3khgQ3uJUZSGwiY7sD/d+xnbdUC85v+clOL/nJQCU+bbZs43lPhYOXY6nvnwUT335KADF9Mk1blwQ5AQRBVouzOzzUtBTKEkUDVvK1CbJRXRk/+0FhzmqHeoVxx6q+KjFb7kKqGxpq4/E7VMw/6O1ABTTJ68ZmuPPB/8jQMuNSy48IugplCzhFddK3CRaCtSNWoMt32zJ23/ozH0t93X0nD749ZxsYEL9rnr0eaCLGlk4+O1+qKqoUtN4SOyxeNhKHDYr69cxf8gyrJ7wPi459HIAwLLLv0Xfh/YKanohhkntWRmSuGs60JC7psSvj9vqSxREIxFC/E/D1ajCkTdfiTLLwhAYiTumIX6D4hedGD8Z+595Fo7qq+RxS9w9A/E/XxLk9EJNmE2iUmCzyKW3XolHbjMfgVQ7rTvmDa5TtYAnDD4Zb+I1y+PGojF079wjx4y6eNhK1G/KLqRhS/NQqojCGqBUprjk1cvV94c80gNJOP+cq1tUo6FJyeFm9L3xQJCa1jUAlIARTk3rmpz3Ruf6hszDFnpWrtuKXl1qLJ0Tv24EEuMmIn7LVdiwdaftsRlj2NmYUk2piTumAU271ONKRJ/t7iUCXFgD8itTjPy/i10ZI51m6vdVyIeRB4MkU2nEolmjnvZ9oXMlClJgs0gxYY0/RI00a2/OtS6sAcoDXevzxgWLPZ/poFvIXuIeq2/cjB4TFEfmucctQVXbFo7NzlbKSukJZh3bdSoosPktvDMwNMkfDKGmmLCWTCkCdywa0dWsdaqxZwoFlIe66PemaoAeecOwkL3EHRL3PYn4tRcAyKQ8iVY4NjtbLSulFc4ak8UFNr+FNUIIc8MJeObDRkRDiOhLIlpFRNd7NY7fsMmFZVz+AK2d1h2AYqZ87VFvHYhvOO5qT/uXQBXWAMUH0S0fQbMJimPRmGfJjN2CMYZ0OmVqM0OxNYQU/pk5voyIDnHrWsp1/UpMLuyfGYtG1Ido/LoRABRTZcuDjvJsTm323NOzviUKXFgDMsl/XfIRtJKguCSqDriYONeL9cuTJwARRQFMBHAigLUAFhPRHMbY516M5yc0prAG4fwrLsIht3VF/4t+A0DRctyxl7cC2zW7/9XT/iX5dGzXyZWC78W0YMlU0tC0uWb9ap0zzKHXpxsmVLfSephcQ04B0CuzHQFgcuavH2OXJPExhSOgv9mwHQCw6P2vUf/JAgCKluPCc37l2ZzOP/2XnvUt0aexKeVKwXcz2igi0jVttqy0P75ef26ZT93Qr3m1fnmlYTscwCrG2H8ZY40AZgE406OxAqFymr5pYMEzb+OUv1yJPX65J9rt087fSUl8ww1hzSxGQpQT4aqqospUf2ZKZXFcLk1lZg05E8C/mMIiAO2IyI0K2WW/fiWmvKy7f892LbFnu5ao/+JjoLo9UOVeFQ9JeHBDWLOCniDlRLjSU9bp9WemVJZOR26o2DxZv7yysXQFINZaWQuN5EhEowFwvWzDPLz8mUdz8QajwgSbAEwDAHQAsGkewpcTzgEdoFxhOdE8r2lXwaPadvuYabodW1+bh5fNVjBvQURituupjLGpwvuia4hBm64A1pucgxFmxs5bwyadc2BJrWGTXv+/Yk06ANg06Zy7fZiNb5Tb//dyux7Am2sytYZt/Hr5a5POOdCNNcyT9csrgU1P/MwRczMXNhUAiGgJY6yfR3MJBHlNpYG8JvdgjA1xsbuia4jJNl6NLdewEqTcrqncrgcI9ppcXMM8Wb+8MomuBSAmquoGYJ1HY0kkkvLDzBri1Toj1y+JROIET9YvrwS2xQB6EVEPIqoEMAzAHI/Gkkgk5YeZNWQOgIsz0VZHAtjKGHNqDjU7tkQikRjhyfrliUmUMZYkoqsBvAYgCiDBGCvkzDW1wLFSRV5TaSCvKYQYrSFEdEXm+BQArwA4FcAqADsAjDTqz42xi5xW8p+5DvKawk+5XQ9QBtfk1fpFJZEbRSKRSCQSiaQZUxbF3yUSiUQikUjKGSmwSSQSiUQikYScwAW2cikBQ0RriOhTIqrjuVmIaDcieoOIVmb+tg96noUgogQRbSCiz4R9htdARDdkvrcviejkYGZdGINr+isRfZf5ruqI6FThWKiviYj2IqK3iWgFES0nomsy+0v6eypV5PoVHuT6VTLXJNcwuzDGAtugOON9DeAXACoBfALgwCDn5OBa1gDooNl3N4DrM6+vB3BX0PMscg0DARwC4LNi1wDgwMz3VQWgR+Z7jAZ9DSav6a8A/qTTNvTXBKAzgEMyr9sA+Coz75L+nkpxk+tXuDa5fpXMNck1zOYWtIat3EvAnAlgRub1DABnBTeV4jDG5gPYrNltdA1nApjFGGtgjK2GEulyuB/ztILBNRkR+mtijK1njH2Uef0zgBVQsmOX9PdUosj1K0TI9atkrkmuYTYJWmAzKs1QijAArxPRUlJK1gDAHiyTVyXzt1Ngs7OP0TWU+nd3NREty5gcuOq9pK6JiLoD+BWAD1C+31OYKafPVq5fpfXdlfz6Bcg1zCpBC2xelZYJgqMZY4cAOAXAVUQ0MOgJeUwpf3eTAewLoBZK3bZ7M/tL5pqIqDWA2QB+zxjbVqipzr5QXlMJUk6frVy/Sue7K/n1C5BrmB2CFtjKpgQMY2xd5u8GAM9BUdn+QESdASDzd0NwM7SN0TWU7HfHGPuBMZZijKUBPIyser0kromIKqAsdI8xxp7N7C6776kEKJvPVq5fpfPdlfr6Bcg1zC5BC2xlUQKGiKqJqA1/DeAkAJ9BuZZLMs0uAfBCMDN0hNE1zAEwjIiqiKgHgF4APgxgfpbhi0KGs6F8V0AJXBMREYBpAFYwxv4uHCq776kEkOtX+Cm7/xelvH4Bcg1zRNBRD1BKM3wFJfLjxqDnY/MafgEliuUTAMv5dQDYHcCbAFZm/u4W9FyLXMcTUFTsTVB+1YwqdA0Absx8b18COCXo+Vu4pn8D+BTAMiiLQedSuSYAA6CYA5YBqMtsp5b691Sqm1y/wrPJ9atkrkmuYTY3WZpKIpFIJBKJJOQEbRKVSCQSiUQikRRBCmwSiUQikUgkIUcKbBKJRCKRSCQhRwpsEolEIpFIJCFHCmwSiUQikUgkIUcKbBKJRCKRSCQhRwpsEolEIpFIJCHn/wPBNNjTqtx2bAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x288 with 4 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, axs = plt.subplots(1, 2, figsize=(10, 4))\n",
    "# force a diverging color map to use the full color space, not just upper half\n",
    "# (as if it is not diverging)\n",
    "traj_contacts.residue_contacts.plot_axes(ax=axs[0], cmap='PRGn', diverging_cmap=False);\n",
    "\n",
    "# force sequential color map to use only the upper half (as if it is diverging)\n",
    "traj_contacts.residue_contacts.plot_axes(ax=axs[1], cmap='Blues', diverging_cmap=True);"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b6b6acd1",
   "metadata": {},
   "source": [
    "### \"Clipping\" at high or low values\n",
    "\n",
    "You might be interested in somehow marking which values are very high or very low. This can be done by creating a custom color map. Details on this can be found [in the matplotlib documentation on colormap manipulation](https://matplotlib.org/stable/tutorials/colors/colormap-manipulation.html). The basic idea we'll implement here is to use the `'PRGn'` colormap, but to make values below -0.9 show up as red, and values above 0.9 show up as blue. We do this by making a color map based on 200 colors; the first 10 (-1.0 to -0.9) are red, then we use 180 representing the PRGn map (-0.9 to 0.9), and finally the last 10 (0.9 to 1.0) are blue. Note that colors in this approach are actually discrete, so you need enough colors from the PRGn map to make it a reasonable model for continuous behavior. For truly continuous color maps, see matplotlib documentation on `LinearSegmentedColormap`.\n",
    "\n",
    "This is very similar in principle to one of the matplotlib examples on [\"Creating listed colormaps\"](https://matplotlib.org/stable/tutorials/colors/colormap-manipulation.html#creating-listed-colormaps)."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "dacabaa9",
   "metadata": {},
   "outputs": [],
   "source": [
    "from matplotlib import cm\n",
    "from matplotlib.colors import ListedColormap\n",
    "import numpy as np\n",
    "\n",
    "PRGn = cm.get_cmap('PRGn', 180)\n",
    "red = np.array([1.0, 0.0, 0.0, 1.0])\n",
    "blue = np.array([0.0, 0.0, 1.0, 1.0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "9d223a64",
   "metadata": {},
   "outputs": [],
   "source": [
    "# custom color map of 200 colors with bottom 10 red; top 10 blue; rest is normal PRGb\n",
    "new_colors = np.array([red] * 10 + list(PRGn(np.linspace(0, 1, 180))) + [blue] * 10)\n",
    "custom_cmap = ListedColormap(new_colors)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "c2559a2e",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAD8CAYAAABKKbKtAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABCNUlEQVR4nO2deZwcVdW/n9OzZSYLIUwSQhJMwIgG1JDkFxSR5ZU1gpHXVw0om0tERQVUZFFAEOEFQUVZDIuACAFfRBZRNpFIWBOIEAgxCQTIQkISSDKZtbvP74+q6lR3eqnurq6l5z7zqc9U37pVdbq6u7517z33HFFVDAaDwWDwm0TYBhgMBoOhPjECYzAYDIaaYATGYDAYDDXBCIzBYDAYaoIRGIPBYDDUBCMwBoPBYKgJRmAMBoOhThCRG0VknYgsKrBdRORKEVkmIi+KyGTXtsNFZIm97Uw/7DECYzAYDPXDTcDhRbYfAUywl1nANQAi0gBcZW+fCBwjIhOrNcYIjMFgMNQJqjoX2FikygzgFrV4GhgqIqOAacAyVX1NVXuBOXbdqmis9gB+0N7errvuuitI2JaEgBNIoT++d4MhCriDmVTxO3zh+RfWq+rwcvcTOVxhvcfaC14Gul0Fs1V1dhmnGw285Xq90i7LV75PGcfNSyQEZtf37cq8efPCNiMUWlvbAOjq6vSlnqH+SGsagITUf4dDWtMF36fX61Dp9co9d7nHaWtte6OsE2ZYD8z3WFe6VXVqZeexDpCnTIuUV0UkBKY/41UwjLD0X/qDsDjk3uCd17k3/5QmAWiQxsx2Z/9Kr1c/uc4rgbGu12OA1UBzgfKq6BdX1GAwGAwA3Ascb3uTfQzYpKprgOeACSIyXkSagZl23aowLRiDwRAKTqsD8rceirUonJZLbt2UJrfbVqltcWzRiMjtwIFAu4isBM4DmgBU9VrgAWA6sAzoBE6ytyVF5BTgQaABuFFVX67WHiMwBoOhZni5UXu9kbu7wfq0F4Amac6q40VcinW9lWtT1FDVY0psV+DbBbY9gCVAvmEExmAw1IzcG3WlYyW5dR1hqaTF4j5WvkH9Quc0lI8RGIPBEBhebtpuEXK8J7d2duQ9hl/dYV5tM5SHERiDwRA4xbrO3OWO96TXW38lYydGWGqHERiDwVA3FBMLd0ulkIOBac34ixEYg8EQGNXcwKu9+ecb2Lfm0xQfkzFiUznREBiNr1ugwWDwTjljMLn1B7YNArafdFzK3dmp42wrNIaTew8K7H40+EXYZ4y3uo/U1hS/iYbAiHlKMBgM2eTeE/JFs/D6YOruFsutX8jlOTAEpMFbELSqY7cETDQEpo7I189rxNNg8Ea58cDc3V7O63wtmmJClE9Ygu1RERIJb+dK1dgSvzEC4zOhNLENhjqhUCyyXHKFp9TvLtKeZQKSqM9w6kZgDAZD4ATRQih0Dr/CyfiFYATGYDAYfMNLS8WvKAC55A7sFzpfYIiQaGoI59w1puQVFZGxIvKYiCwWkZdF5Ht2+TAReVhEltr/d3Ttc5ad13mJiBxWyzdgMBjiTaluMHe9SsPLtLa2kdZ0Zsk9Zu7YTaDYXWRelrjh5dNKAt9X1Q8BHwO+bedqPhN4VFUnAI/ar7G3zQT2xMoNfbWd79lgMPRzcm/wxcjXgqlUBLq6OrPEJB+OG3TQCJBoEE9L3CgpMKq6RlWft9e3AIux0mvOAG62q90MfNZenwHMUdUeVX0dKyz0NJ/tNhgMMaTS7i2nG83tFVZKbMoVo66uTlpb2zLxzwKjjlswZY3BiMg4YG/gGWCknagGVV0jIiPsaqOBp127OTmfc481C5gFMHbs2NzNBoPBkCFXlGoxXpLWdGiZY+MoHl7w/CmJyCDgLuBUVd1crGqesu3mB6nqbFWdqqpT24e3ezXDYDAYPLVgEpIgpclMeuVShDetQJBEwtMSNzy1YESkCUtc/qiqf7aL14rIKLv1MgpYZ5cXyvlsMBgMhhxEoKEpfuLhBS9eZALcACxW1Stcm+4FTrDXTwDucZXPFJEWERkPTACe9c/k6FHoaaqcpyeDoV5xxjWckCx+4HUsp0EaIzXnJS/9fAzmE8BxwEsistAuOxu4BLhTRL4KvAl8HkBVXxaRO4FXsDzQvq2qcYtwUBaFvuiR/2L3U7ykzDWUT6FwLR2dVo96g/ce+X6G+Nb9JSKHA78GGoDrVfWSnO0/BL5kv2wEPgQMV9WNIrIC2IIVkSapqlOrtafkHVBVnyD/uArApwrscxFwURV2+UaYN5AoxyJb3/c27U07h21GKJhwPrWh0HV9o3s5ALu0vI8EtQ0oWetZ+q2tbf47ApQR7LLoYazpIFcBh2ANVTwnIveq6itOHVW9DLjMrn8UcJqqbnQd5iBVXV+1MTaRe8TuTG1LjdqSGAB4bwk47oXuL0CYN5Ao37z8FBfnM2trCGcegSHavG/AhMDOVStxcR4WHVdmsNI4+/Eb9zFUzDRgmaq+BiAic7CmjbxSoP4xwO1+nLgQkROYam5SQboYltM6qfdumFoKS1AxqyDaDwRxJWpxvyqlUBpnX747/gW7HA285Xq9Etgn7ylF2rAmwp/iKlbgIRFR4HeqOrtag+L/yYdEpSErDOURxLUzn0/tqAdxqTkiNHiPRdYuIvNdr2e7hMDTFBGbo4B5Od1jn1DV1facxodF5FVVnevVsHyYT7+GRPnJuF6eLP0gyp9T3OlJd2fWnS7vWhPk5xlCF9n6IoPv5UwRmUlO95iqrrb/rxORu7G63KoSmEj/ouLu5lvIlbKamErFjlXOcY24GIKgJTEgswRFtdGWveL+rVV7n/LJTfk5YIKIjBeRZiwRuXe7c4nsABzAtqkliMhAERnsrAOHAouqelNEvAVTrzdBP7/8xiPKP+p9rCwM+ss1repeJf4EslTVpIicAjyI5aZ8oz1t5GR7+7V21aOBh1R1q2v3kcDd1rRHGoHbVPXv1dpUn3dwg8Ej/eUGGBb1fG0Htg3yxbFIwHPK5FKo6gPAAzll1+a8vgm4KafsNeCjvhjhwgiMoV9TyxtgWGM7URpTqgcBLzRe6ZfXqojQ3FSft+L6fFc1ppIfTb45Oob6Jqwba5Ru6FGypVIapLFg5ktfxlIFGmIYyNILRmAqoNiPZlPS8vpraxhER8oKkbFjY3towlLo6asmM5INhjycPO4HAPx82RkMaxxRonYwlNvKK1TPHy8y8a2LLGrUjcBEpSm+Q+OwzPqOjeGnISg0+NjRuRmTZtTCiUQwINEWie+QF6LyfXcoZE+f9nLZ8h8DMDAxJGizClLo2jneYIE6GAk0xDBbpRfqRmCC/LEV+3G7m8yKtR5Fb7go2hQWcQxxEyVxgcL2NEkzTQ21jT/mJ2H8LgTTRWZwUezHnb2tPr80QVLrJ/XcPvSo3bjjTp/20iTxEZgwEBEjMIby6U5bYxwDEm2hehTlO+e7yfWR6MIrhdfrVen1VdKk7ZamuRH6TwONrO1dCcDQxvZAJ1zGBRFoMl5khnIZkGjLrEfNoyhq4lKt00Gl17dBGs1YVA1JSIKRzWPCNqMqvOQPqs6bzLRgDAaDwVADxAzyh0OUJox5xe2F4tivpBHXeEwU3k/UvJCMy3R9ktJkJlWyu0UfJ7yEY6r2t2TclEMgSjdAr7i9ULbZH733Ecdra4gfDdJYdx6L+VJDV/N7EhEaGurz91jykxeRG4EjgXWqupdddgewh11lKPCeqk4SkXHAYmCJve1pVT3Zb6MN0SKOLU2DoRJyW/5+fef78xjMTcBvgVucAlX9orMuIpcDm1z1l6vqJJ/sM8QAIyyG/kItvusi0n+9yFR1rt0y2Q6xYjt/Afgvn+0yGAyGSJPrXVYpItDgT8rkyFGtHH8SWKuqS11l40XkBRF5XEQ+WWhHEZklIvNFZP76d9ZXaYbBYDAEi19dZWK7KXtZ4ka17bJjyE67uQbYVVU3iMgU4C8isqeqbs7d0c4jPRtg8pTJhfJGGwyGGJKV7RHLs9JMZC2A0H8H+QshIo3AfwNTnDJV7QF67PUFIrIc+AAwv0o7Y4OfA95m8DwYHDdacwOsnNyo3e7vrGrxSYr9HRNNOT8HA6+q6kqnQESGAxtVNSUiuwETgNeqtDFWFJrlW8kPy/wYgyHztF2f3eCBUMgVOa3punNTduOPmzL9N+GYiNwOHAi0i8hK4DxVvQGYSXb3GMD+wAUikgRSwMmqutFfk+OHEYpoY+Jj1Q6/5opEFb/ywTQk6jNgkRcvsmMKlJ+Yp+wu4K7qzTIYDPVAf+kWqyoWmUCjT11kInI48GugAbheVS/J2X4gcA/wul30Z1W9wMu+lVCf7bIQqeenNYPBK06KcKi/MED5fuNVe5FJ9S0YEWkArgIOAVYCz4nIvar6Sk7Vf6nqkRXuWxZ1fxfs014enT2fR2fX3s8gpUl6tZte7Qbg7Cd+wtlP/KTm583lviv+Ffg5DQY3XV2dmcULaU1X1AroSXdvd5yNyXVsTK4r+1hg3S/6tJe0pjPruSQk4esDpNNF5mUpwTRgmaq+pqq9wBxghkczqtm3ILFvwTg305bWFtIp6wva19uXWXdzz2WPM+OHB1R8rrnvPALA/sMPzpSlNMk/r1sIwF4n7pwpH9Dcxr677V3xuaph7fRXsKYoFebq42/jW7ccG4xBdca5T/2UTy2xfnsHnDiJp++wHvKamhuYcvQexXYNlLSm+dfNL2ZeNzVbP/d9j90rLJOycLcEnOyvfZos6c3Xne7MBM58L7k+Kx3A1vRmzn7/pQBcu+IXBY/hCFPu+NvqnjcAGNMyniFtQ4EAWmACjd7HYNpFxP20PNue8gEwGnjLtW0lsE+eY3xcRP4NrAZ+oKovl7FvWYhq+FNQJk+ZrPPmzQvbDEM/Zf0rlh/KO6s3M3LsUACG7TE0PINiyPpXNtLYZN0kh07YoWA9p+tsa2dH3hwrXsZs3C7RaU2TYpso5WbQLBSY0i0w7gjo7v3KbaW0tbYtUNWpZe0E7LzHzvqlq4/3VPeKgy8reA4R+TxwmKp+zX59HDBNVb/jqjMESKtqh4hMB36tqhO87FsJsW/BGAzVMuxDQwFonzisysRR4eB+og+clPWvfeKwsnarZvwid75Ngm2Cktv6KXQed8slfDdqKacFU4yVwFjX6zFYrZQM7knvqvqAiFwtIu1e9q2EsK9sFu6nD+1TpKn/TUzQPqtFufmtLeyw25BM+cXzLYeO70/4Ps07NIViW71Si+i4QRJmnhWxY2gp3npCSnU3ReX6B2mHgF9uys8BE0RkPLAKaypJVj+4iOyMFd5LRWQa1jj8BuC9UvtWQqQExv304RaXvi1JmgZHytSa4bxvt7gAnDX1zDDM8UTn210AtOzQTENrvP35Re0bpoTfdRwHSl2n95ZuKtplFlWCdK8W8SdlsqomReQU4EEsV+MbVfVlETnZ3n4t8D/AN+25il3ATLXGSfLuW61N/eOubTAYDBHGr4mWqvoA8EBO2bWu9d9ipV/xtG+1xEJg+kvrJa607dwatgm+0f1ej7WSVh6/898AHPrNaSFalJ+NyXWMHjwOgK3vdZBoiUbXUi7JvpSnevkG2ivBr5ZHoF1k4tsYTOQwd+6Icf6UywE4/a5ZDBk3OGRr+h8tO27rpo2isDgMaxyRGc+oxa3w6Y1z+diw/SvaV1TQtNV15nXw36+B9qiM4ZSDIDQ31ue4amQ/je53tk2g+uPym8IzJGDOffI0zn3ytIListfFh2574e3h0GAom0rFBexxmQasxVASpwXjZYkbkW3BDBi+zY3wy7ud5NlLJe6U6upYdNZD215E7PuW7klHtqvGK21HTQBg0qETefI794RsjaE/4KMXWeSIrMC4MR498SDu4gLQed+25KzpHntiXsTfl1vYF9y9xLdoAsajLiCkfqMpR/uXEyCnffBcTvvguYGe86iBMzPrP/vEldyx4lZEBVFh39+UDgO0tndlyTqGykm0JCIvLpAtgFWLS4pM16uKGnEJAMGKpuxliRuxaMG4n9CWPv4WEw4YW2KPwjhPZZvf2sLgXQdlyq946afVGVkB922dk1n/8bzvArA1tQXAU/fM8xsWcMQoKw5T76a+aE3ATBG5LjxDcd6av5a3/vMOEJ14Zf0BEaG5sT6zqcZCYNxPaNWIC8CcN/4AwBfHfTmrPCpRA9oaBpWuZHPEqG2tnEiJC5DsTtI4MBZfr4I8dsPzABz01cm8/ODr7HnY+JAt8oBL2L2GkOl514oYvMcnx9HZNbJoXVExrRqfEf9CxUSOeN8BKiBXWAy1Ie7iApawOMRCXIBJlx/OwjP+DngPIdPQbD3AdXZ1laybKy5Pb5wLVOd11t8RwZeZ/FGk5LsSkRtFZJ2ILHKVnS8iq0Rkob1Md207S0SWicgSETmsVob7SWtrG1cffxtXH39b2KYYDFXhiEspLj3sWi49zJrg3TiwsegDwZO3Lcoam3HzsWH7G3GpGivhmJclbnh5zLwJK7TALTnlv1TVrIQLIjIRK0jansAuwCMi8gFVjcyMjXyeMRte38Bxw3vCMskQM9zh4qNMcmuSpjar6zSdTGd1A//o798EigepbGu1IjQs+ddwM55WQwSJ5QRRL5T8lajqXBEZ5/F4M4A5qtoDvC4iy7AypT1VuYn+kq//2Ap1Uj/hTgy1Jcri0rFyK0sGWUnGpgz9eEZAcscYS42jJLcmeXf1JgBadiw+LmOonnoVmGre1Ski8qLdhbajXZYvK9rofDuLyCwRmS8i89e/s74KMwyG8Fj+xCo6Ux10pjrCNgWAQWMGsvzKDpZfWZ09jQMbadmxOSt0jqE2iAhNiSZPS9yo9FHsGuBCQO3/lwNfwXLpziXvo5Kd5nM2WBktK7TDH3Jcah9ZawUU/Z9Tf8j0L3wcgNuOvj4EwwxRZ/f98j4/hcoXzj2kov1Ehf89/BoAznjwZD9NMhRBqN8WTEUCo6prnXURuQ64335Zk6xotaYv0UuTKyvewSMtn4X3bp/OG91Lt6sfZK4IQ/Rxxiq8eGFFGRU1whIK9TsGU9G7EpFRrpdHA46H2b3ATBFpsTOjTQCerc7E2nPJs5cW3LZw/UIWrl8IwM7fmGItX89OiX3G3LM5Y+7ZtTTREGFeengZLz28LGwzIsuTty2yPNEMBUlIwtMSN0q2YETkduBAoF1EVgLnAQeKyCSs7q8VwDcA7OxpdwKvAEng22F5kN2x4lZg+3kv+Z42Lz7wIn5S4OlzxpjPZ9abB1mtnDcvz/ZZuHT/n1dvsCE2tLW2Zn1/othNVgm1SlPe15v0/Zj1hNhuyvWIFy+yY/IU31Ck/kXARdUY5QeFJlTm68bw2rWRKyyG/klXZzc9m/qA6EVQqIZEY6LfRC2PFFK/YzD1+a4MBoMhJghCozR6WkoeS+Rwe5L7MhE5M8/2L9nevy+KyJMi8lHXthUi8pI9eX6+H+8tug79NaK11Qqf4WQDNBg84fI0VNG6ark49HX2Zc3o705bvxGvIWcK8ckTPlLV/vWOX11kItIAXAUcguVw9ZyI3Kuqr7iqvQ4coKrvisgRWJ68+7i2H6Sqvs0b6XcCY4QlfGLpdVXg99+3JUnT4NI/o3eT1m92ZecbfHjIFD8tq4o7L3iY3b9rBVidMvTjWduqFRaHMLt/0poO3QYv+GTfNGCZqr4GICJzsCa/ZwRGVZ901X8ay9O3ZkT7qhvqks6urliIy+M3LeTxmxYWreNFXAB2bGxnx8b2SIkLWHNm9uj4CHt0ZLcyJl16eN74YwvuXhKccT4QB+8rZx6MRy+ydmeCur3Mch3K80R3m68Cf3O9VuAhEVmQc9yK6XctGEM4xLHVsu8JEzPrfzj7rxx/0ZGADxkenZt2RByHBo0ZuF1ZoaCZfmXLNLiQsubBrFfVqQW2eZ7oLiIHYQnMfq7iT6jqahEZATwsIq+q6lyvhuUj2tJuqBvi0mpxSG5N0iTNNInlmn7czz/tX4bHBjyLy31X/Kv68/mIk0a6VHl3ujMzhmMojuBbqBhPE91F5CPA9cAMVd3glKvqavv/OuBurC63qjACY4g9F+zzS9+P2TiwkSu/eDNXfvHmTOsrDI46/ZP89thb+e2xt4Zmg5tCKaQHDs1OlDcg0ebbGE5/wKeJls8BE0RkvIg0Y0W2v9ddQUR2Bf4MHKeq/3GVDxSRwc46cCjbJtBXjOkiM8Sec585rSbH/e4dJ1j/Cbfldcpt0U+St2rLCoaFbQTxDOPkV7h+VU2KyCnAg1ht5Bvtye8n29uvBc4FdgKuFhGApN3lNhK42y5rBG5TVW/JhYpgBMYQG95buom2EdZTca3dhOe+8wj7Dz8YgI7UZgY1DKnp+eLOsMYRYZsARN9bLC8+TrRU1QeAB3LKrnWtfw34Wp79XgM+mlteLUZgDLFh6IQdAjuXIy6AERdDTenXoWIMBkNt6El305IYUPH+zrhMHLrQDMVJ1OlwuBEYg6EE61/ZSPtE/0cYqhEX2CYs913xL446/ZN+mGQIgXpuwdSnbBpij6ggmu3WP3vWnZz+2Bmc/tgZgdry0O3PBHq+cjHiEn9EEp6WuGFaMIZIkm++yazZXwC+ELgtx154RODnNPQvJO8cyfhjBMZgMBhCRJC6HYOpz3dliD1tra2hTnB0897STVUf41ef+z2/+tzvfbDGUHeI6SIzGALlvbc322vhh8Wv1j36nsseZ+pnoxXk0hAlhIaoBKbzGS8pk28EjgTWqepedtllwFFAL7AcOElV3xORccBiwAm5+rSqnlwLww31TT3lW5nxwwOq2r8n3Q1U73VmiCYC2DPo6w4vba6bgMNzyh4G9lLVjwD/Ac5ybVuuqpPsxYiLwVAlLYkBpcUlT2h9Q3wQEp6WuFHSYjtc88acsodUNWm/rHnSGoPBYKhnvMlL/Fo5fkjiV8hOWjNeRF4QkcdFxDjoG/ot3e900/1OdzAnc6UAeGnzAl7avCCY8xp8QOq2BVPVIL+InAMkgT/aRWuAXVV1g4hMAf4iInuq6uY8+84CZgGMHTs2d7OhTtnr4kNZdNZDgBX59vynLwTgxjv+xspfPZ2p50yy9CX/SkgMGF56zOSJP7wIwPy/LODUu07y5bxj2t7ny3EMwWCNwcRPPLxQscCIyAlYg/+fUlUFUNUeoMdeXyAiy4EPAPNz91fV2cBsgMlTJsf3LmIoC0dcwIoge8HHzwPI/HeIs7CUw4a33wXwTVzASs9siBP160VWkWyKyOHAj4DPqGqnq3y4iBVUR0R2AyYAr/lhqKH+OOmv3wzbhKp5+cHXM+utraUTbOV2mc344QFVe5kZCpPWNGnNn4UzSoiIpyVueHFTvh04EGgXkZXAeVheYy1YeZthmzvy/sAFIpLE8mk5WVU35j2wod+yMbkOgJumX4vmTxkeG/Y8bHxmvaurdIpgL91mBv+IQ34YgViOr3ihpMCo6jF5im8oUPcu4K5qjTLUN05yqmm/OpJnTr0vZGsMhrCJp4eYF8xMfkPgvP7kagCe/d79sW/BFMPpMvPSsqlnFtxtzbuecvQeIVsSXep1kL8+35Uh0ozfdxfG77tL3Q/kd3V1Rkpclj+xiuVPrMqK8daZ6qj5eft6U/T1mlmghRAf3ZRF5HARWSIiy0TkzDzbRUSutLe/KCKTve5bCaYFYzB4pNwWieMA4B6nCYK+LdYc6KbB2T/v3fcbDUBnV1embM1Tm9h9v0E1tedjX5xY0+PHHsGXhGO2g9VVwCHASuA5EblXVV9xVTsCy/lqArAPcA2wj8d9y6ZuWzCLH1nB4kdWhG2GIQf30/OUyz8doiXlU26LZM/DxgcuLmAJiyMuvZv66N3UV7CuIzq15PGbFvL4TQtrfp44Ix7/SjANWKaqr6lqLzAHmJFTZwZwi1o8DQwVkVEe9y2bum3BfOjgcYD11Om+KZh+8eDRPkWarB+H++n5udPNAH+taRnSDEBr64Csa5+PlCZpEP9vCY/ucQ8ABzDJ92PXA/mytxahXUTc8wpn23MKAUYDb7m2rcRqpbjJV2e0x33Lpm4FxiFXSIywBI8jLg5P3rYIgH2P3SsMc3wl9wEmLApFPnBelxIXoCbiAttPojVkoyj2XHUvrFfVqQW25VOp3AMXquNl37Kpe4HJ5YV7lgKw94wJIVtS56Sg0OTkehAWhyiIC8At59wPwHE/d3U7FvkMvBDWGFJ/RNO+OLysBNxxt8YAqz3Wafawb9n0O4ExwhIM6UQ6kwa2T3tppgXwHgImZQfrdj9Vi0rdep69t3RTVmKz235ixY899Jh9aJ84LFPekbLC+g1qGJIZz+rs6soWFocqx42NsASEgqZ8+V4/B0wQkfHAKmAmcGxOnXuBU0RkDlYX2CZVXSMi73jYt2z6ncAYgsE9g7pJmsue75Kvu6ZexQW2z5p57IVH5K03qGFIZv2Sz1xrr51An/YC8OTNr3DAiZN8semxG54H4KCvZjxZSfekSbRs7xv00DXPcug3p/ly3v5IGV1kxY6RFJFTgAexHi9uVNWXReRke/u1wAPAdGAZ0AmcVGzfam0SP95YtUyeMlnnzZsXthkGn9E+67uVOwZTKb/49HX84K9fz7zOtHLSjVlP68mtVnmiLRGLUCHlMPedRwDYf/jBmbLk1iSNA4N7Vtz3N5Zz0ZPfuSdT1rOhh5adWgKzIYq0tbYtKDI+UpDJkybr4/+Y66nukJ0GV3SOsDAtGEPN8EtYHNziAq5WTk5XUJA326BxC4tD0O934UP21IjvuAoTkuVo4ASYfO8/mxm2x9BA7YslEXjQrwX1+0s0BIbTkhASaLf1Q2loDS78eJ/20iTNLoOsf90bu7OCS/7h7L8C1kD4ukUbABix106B2VkvdN63dLuylh2zu0GdluPat96rucCkNZ05n3s9Xz23bVFBUb8G+SOHERiDwWAImTptwBiBMVRP1oB8a+F6taIjtTk7yZbdeNrw1mZGu1owbk8r03KpLU5rYfguQ0rUhJ50d6ZVkdUS9UCf9tLguo11pztpayge+sadHyYhifBbNgqajH7OmkowAmOIPTs0DMtbPmbvkXSltwIwIFE6GZjBP5ybtdvFuhAticpz5CRyol0VEpdiXWf5yovVrwX12oKJVmekwVABhW4EKsobXct5o2u5p+M44zKG+NAgjZ6EoFyxCLQ1o0BavS0xw7RgDGVTq5hVtWBC456e65puM4NDkC0YxbRgDAbuvuQx7r7kMRpporW1jdbWNl64ZynPvvsEz777BAD3r/4z96/+c2g29qS7eeyG5zOTBBMtibwTA/PheJkZDMGOx1heZF6WuFHyMVREbgSOBNap6l522TDgDmAcsAL4gqq+a287C/gqlrPod1X1wZpYbqg5fVuSWTlFjj7zIGtFt8Xg2huwUktYTB/12cDsy0djTxPpVGUDpnnDrRgMtUb9mckfRbzI9E3A4TllZwKPquoE4FH7NSIyESuGzZ72PlfbiWwMMSQ3YZVDsZAtCQlm9vxj6/I/tzQOaORTs6byqVmFJzune9Is2foSS7a+RHc6GsEqDf0bTamnJW6UbMGo6lwRGZdTPAM40F6/Gfgn8CO7fI6q9gCvi8gyrEQ2T/lkr8EAwEEjDsusJ7cmaWprArzFK0u0JNij5cOAFUBz5fNrAdhp7JDMxMygvYgM/Rj1LZpy5Kh0pHakqq4BsCNxjrDLRwNPu+o5yWy2Q0RmAbMAxo4dm69KQfLFQjLUns0rtgAwZNxgT/X/uPwmAL60+4k1ssiicWBjyWCaj86eT6LBEgx38EYVZfRk++vrShuf7kyTGGgExhAQ9akvvnuReU5aY2dhmw1WsMtyTmKEpTLuenMOn9t1Zsl6r3QsBGDH5mGMat41U966U3mzKGstLF5w4mMdeNzeJFsKpw4GsmKa1XM8M0O0UEwLJpe1IjLKbr2MAtbZ5V4S3hhCwou4AEwcNClveaExmaBIaxqsqPSePMPc+WMaWhtoqDZBiiFy1EdXZlkZLWNFpZ/MvcAJ9voJwD2u8pki0mInrpkAPFudiYa40p3u9HUQvdwbSV9niRYL1g3q3eR63k2uz+RUAbK6y6LKBfv8kgv2+WXYZoRK/MUFUEgl056WuOHFTfl2rAH9dhFZCZwHXALcKSJfBd4EPg9gJ7e5E3gFSALfVlXff6o9G3oA+n3+CS9sTK5jWKM1xqB9yoX7/QqAc585rexjLeu0wrS/v22ip/p/+vFjAMy86JCyY0wVwuucFvDWzZWQRHYcM5tUIpkV4yqKVPIZGqKH1vEgv0k41o8IOjGVH1z/6u/42ge/UfH+jjuz2+vMYKgFlSYc+8jEj+r9tz7gqe77powxCccM0SRO4nLHilsB+PoeJ5f0EEv3uKLj5rRw4iospz92BgBXHHRpVnnvJqvbr3mHpsBtMtSKeM7S90IddGAaosqWNzvY8mYHba2tmXWvfHHcl/niuC97ntdCM9biE7/49HWZdccTLUg+ePtUPnj79g+qnes66VxnJofWE6qQTqunpRpEZJiIPCwiS+3/O+apM1ZEHhORxSLysoh8z7XtfBFZJSIL7WV6qXPG55G2AGfMPRuAS/f/eciWRJ/bfvI32nexAjoe+s1pNT/f4F2t0OmdXV1A8RwdWaTYLg1yKfwe7HWnZ/Yicn4za/YX8pYPnbBDwJYYgqBa8fCIE4HlEhE50379o5w6SeD7qvq8iAwGFojIw6pq58nml6r6C68njL3AGGHxzrEXHhG2Cd4oIi5Oa+K6JddWNTZjqB+qdVUOO+GYKqT6AvEQKxSBxWWLrgGcSfRbRGQx1mT5V6gA00VmMBgMYaJlRVNuF5H5rmVWGWfKisACjChW2Q4RtjfwjKv4FBF5UURuzNfFlkvsWzCG/sWcN/4AEGjrxY/8N07LK7e7zRm0H7rzELsrsXh9w/ZU2/KIwlyadNpzC2Z9MS8yEXkE2DnPpnPKsUdEBgF3Aaeq6ma7+BrgQqzgAxcClwNfKXac2AuM9lk/QGmSrDL3a0O47HbO/rx20dzKD+DMpGqwBv+Dxo/kaoWEwvEGc4tLsfqG+kPxbwxGVQ8utE1ECkVgya3XhCUuf1TVTHInVV3rqnMdcH8pe8KX7ipJN6ZIN2bP5TTikp9VvSs496mfcu5TPw30vFWJC1hjMibKi6FesSdaBpBwrFAElgwiIsANwGJVvSJn2yjXy6OBRaVOGPsWTFxS90aB0c3juODj54VqQ0+6m5bEgFBtCIMxp34MgJW/erpETYMXWlvbgG2J7+JN9S7IHskbgUVEdgGuV9XpwCeA44CXRGShvd/ZqvoAcKmITMJqdK0ASvZT183duU97M+FI6iMAXn3SH8UF4CtfjIkHX0yoD2GxUDsWWe3PoxuAT+UpXw1Mt9efIH9UfFT1uHLPWTcC4451ZcTFUEvSmub/XXEUAAu+/9dMuagw7ddHAvDs9+7nxAdOBuD3n76G8z/2k+ANNcSGgFowgVM3AmMw1Jonb7O6nPc9dq8sYXF47alVPHPqfQBsSK7l95++JrMtyIeetlYrb0+u44AhotRxsMu6EpiwJ0wZKicOYzP7HrtXZj3fTXz8vrtk1p0I1g57XXwoAIvOeqiWJm5nkyH6aHBjMIFTVwJjhCW+eBGXLW92ZMLPBIk7cZlDuTfxWguLabXEGK3fLrLY3ZGdeS8OH/rZwXzoZ9u7fveku+lJd7P8iVVZ5b2b+jKT2wzRZLdz9s9bPnKP4Vmv+7Q3O0lYjeilJ+t17ncwCnR2dRlxiSkKpJIpT0vciF0Lxj3HJa1pFv/4kbz1nCfi3fcbnVVuwpxHn0LzZtYueQd30Mw55zwMwOd/dhADEm2Z8j8uvwmAL+1+YsU2uCfr5iZLM/OsDL6i9RuuPxYCU8jtuJIuMdOVEF9yu8eO+/mn89arRlgcch9kEulEv5vsmdyaBKBxQGO/e+9BU69dZBULjIjsAdzhKtoNOBcYCnwdeMcudybpVEyukBx5y5cAuP/4P5Z9LCMslbFw07PsPWQfIDuMSVyyZKY1zQNr/gLAkbv8d946qa4UDa3b30kT0v/EBeKVoC7OqEI6FUg05cCp+BukqkuASQAi0gCsAu4GTqLMnAHlUomwGKpj0g7T8maWjMtNKCGJgsLiIANM15cXnAyiudlDU13WGEE+kTYUx7RgivMpYLmqvmGFsqmcI35/DAAvPP0f3v7dAsDMzDf4w7PvPgHAtB33y7tdSRNDv5fAyRUWh55NlsOF0w3tJ0HcAzpTHbQ1BO+laI3BmBZMMWYCt7tenyIixwPzsbKjvZu7g53HYBbA2LFjM+V/O8k+zEnb6iYkYXKRh4yosPbl9QCM2GungvUy/fYRbNk0zR1prczIv93EtauOtp39FxaHIB4wQxEXrC6yZF/8PMS8UPWnJiLNwGeAP9lF1wC7Y3WfrcHKGbAdqjpbVaeq6tT24e15j53uSWea441DGmgcYprexWhtbcsEAfQbFWXEXjsVFRewhCWK4gKw78yPsu/Mj4ZthsGQhTPR0ssSN/y4ExwBPO/kCqgkZ0Ah3E1x00VWmnoKAOg3okJ3p+XgkdI+01oxRAczyF+UY3B1jzkJbeyXnnIGGAwGQ3/GjMHkQUTagEPIzgtQds4AL1z7ytUAnDzxW2Y8xgOTLj0cgIVn/L3iY+QLkZLLukUbsrrNgkz1u6zzFQDe3zYxU9a3JUnXBqulMmTc4Ey5inL3JY8BcPSZB9XcNoPBK6rx7P7yQlUCo6qdwE45ZWXnDChEZ6oDsAbfTt7jW5nyuAiLOylSSq3B7yC6ZjpSm/nDh6/2XL+1tS1v95oXkRi5Z3uW+/ILm58BLLfmWuMWFoemwY00DR6cp7YRFkNEqeNB/kh3RGd5dcRwfN990w6yz39QwxA+fMQQz/W7ujo5YdS3Abh5zVVlnStXhIIQFoOhnrAG+U0XmaHOcNIbdKY7+ORjHylZP8jur0pZ0/smo5p33a68b0uSpsHm626IJiYWmSFW7Pr9jwOw4ufzMmVb13Zm4nmle9L0dVhjWYN2GsLXPlh6qCzKwuLwbu/GvAJjxMUQWQLyIhORYVjhvcZhjY9/ocAcxRXAFiAFJFV1ajn7uzG+v3VEd7qT7nQnaU2z4hfzWPGLeSRaEpnFHSwy0ZKgZacWWnZqCdFi/5k4aFLYJhgMZeF0kXlZquRM4FFVnQA8ar8uxEGqOskRlwr2B4zA1BUDEm0MSLSRkERmqYbrX/0dJ4z6dmZ8xmAw1AA74VgAEy1nADfb6zcDn631/kZgakgtZ9aXojPVQWeqgytm3FDxMWZOOIab11zFzWuuCu19GAz1jiqk+lKeFqBdROa7llllnGqkM0fR/j+iQD0FHhKRBTnH97p/BtMxXUPCnFnveOC9teSNsvft2WBlcGwbtq1LzUQJMBhqRVleZOtzuq2yEJFHgJ3zbDqnDIM+oaqrRWQE8LCIvKqq+bMAlqBuWjDORExDNuc/8ANEJeMB5oWmQU00DWoy4XkMvnD18bdx9fG3+XY8x/sR4N3k+qxtKU166jlIaTIzN83N+r63Wd/3tj+GesTKB6OeltLH0oNVda88yz3AWhEZBVbEFWBdgWOstv+vw0rB4sw98LS/m9i3YK7/1l0AtA4aC5eGbEwE2WG3IXnzuBSjUDj2WnPuUz8F4IKPn5dV/tA1zwJw6Dez59hcsM8vrf2eOS0A6wyV8q1bjvX1eO4Hnx0bswPlNkijp9Z2oXlp7U35Hv5rT0DzYO4FTgAusf/fk1tBRAYCCVXdYq8fClzgdf9cYi8wJ11lxV43wQv9Y8ubVgSF3BTFteYbU07IW75+9YbtypJbk/zkiVO3K7/rzTl8bteZfptmMNSO4PLBXALcKSJfBd4EPg8gIrsA16vqdGAkcLed16sRuE1V/15s/2LE/q688lkrM/OKxWs54MRJ221f/sQqdt9vdMBWRZ9iccYGjgxnQH9087i85cdeeMR2ZYVSAhhxMcQNBU/dX1WfR3UDVnLI3PLVwHR7/TUgb06LQvsXI5YC48QoW/XEu0w4wEpW9r59RuWt+/WXv8c/9vu/wGxz445FFjXee30T50+3slqP3eN9nH7PV0O2qHw2Jq0u4GGNJZ1ZDHVKvkyXQcb98wVVUr3bjwfVAzH5BLLzfTd3DwDgM/86icUHPFJ0v8dm3VX2GER/YIfdhvDLV62uVUewwfrBjjv7EwC8eflTmfKX/vYax71kBRytJkKznxhhMeRzRImNsNg4g/z1SGw+iYbWbdEune6RxT8uLi4O+WJovf/cA1l2wT/9MzAP7765EYB35q9l9JThAPRpLy2JATU9b7nkpop1C4vD+EPbWXhENITFEC+uf/V3AIx8YCJHnf7JkK2JIApqEo4Fz7q+1QCMaNql4mMUGme47pSzKj6mVwYMt4TkW+NO5E8rbsoqiyq92s0A2X4MZlBDdnTmKHf/GaJFJs7dB8O1wyFft1q4aN1mtIzSVd6OEU27bCcujltytfzX8MOzXqe6UqS6UnxlzHd9Ob6bR9L3Io0JpNG63Kd98FxO++C5vNb1qu/nqpYBCW8D/F1dnUZc+glOgr96IVriYnWRaSrtaYkbkW7BGAwGQ92jSqrHJBwLhH/fvwyAjx75/rzbv3b153w5T27XmTPGc+PKK305vpuurk5Id2deO4PrBkMciEsGWa9ErYtMAwrXHwZVCUy+vAGV5Axwk09YkluTBec9+MVJf/0mAL//9DU1Ob4zsN+1rovWEa2Z8kzYi15rBr0TpiKsGcUGQy5/W3MPk3eaAsDI5jHbNqSIZabZKIkLYE20NAJTkINU1R0QyMkZcImInGm//lHRI6SKi0itxQVqJyy5tI5ozXK5lqTl4SYt1n8jLIaoccSoGfk3xFBcokq9uinXQsrLzhnwwr8TDGkfnFWW3FqfE4/AEpaG1ga61nUhTYI0WeKS1jRb3uzIhGoxGKLMXhcfmre83pwCak4dD/JXKzD58gZ4yhkgIrOcnAZjx6xj8/otPPt/27yqgmi1hE3riFY6Vm6lY+VWwGq6D951UOAxwAyGSlh01kN5yy9fennAlsQbtd2UvSxxo9q7+HZ5A7zuqKqzgdkAU0V0cPsQrt3pa0z7H2uQPaXJ2M3IrYRBYwYClpu0ezKpwRBXzppaMpOuwU0aUr3Gi2w73HkDRMTJG7BWREap6hqvOQPSk/ema948bnSV9QdxcdPQ2mBExtBvcZxdIjcAHwBK/Q7yV/xpishAERnsrGPlDVjEtpwB4DFngMHCiIshquz7mxns+5vtB/s3r9jC5hVbsso2vba57OMnJNEvxQUA203ZdJFlkzdvgIg8R5k5Awze2fJmR2aMZn3f28brzFBzOlMdPHXKvQBZgWPTPWmu+NxsAM5f8P1M+ZCx2Q47hlLUb6iYigWmUN6ASnIGGPLjdmd2GDh2WygXIy6GIGhrGMScFX8A4IvjvpwpT7QksoTFwfGKNHjDChVTn27K/WugowhbUu8BMLhhaKh2uHGExT024+5GSPekQ0tvbOhfLP2SFRmceeHaUZeYmfy1JdWTZvOKLQwZF17TOkrCkktDa0PGlXnQmIFon/200xyiUYZ+xY/n+R8E1mChdZxwLBKPvw0tiVDFBazJYb2b+hj02T1CtaMQg8YMZNCYgVmTM92tmShGZjbUD0cNnMlRA0066lqRTqc9LdUgIsNE5GERWWr/3zFPnT1EZKFr2Swip9rbzheRVa5t00udMxIC48Y9k/27E87huxPOCeS8zTs00bxDEzdf+dNAzlcpTqgZZ3zG4Td73xaSRYb+wH1b53Df1jlhm+ELmfh/UUE1EIFhWxivCcCj9uscU3SJqk5S1UnAFKATuNtV5ZfOdlV9oNQJI9FF5sY9i/3KpRf5csweO5KxO5Nk76Y+GpotfXUPon9u15lsSlr9zTs0Dit57NbWtkxeFPf6xuQ6utLW+ujmcdW/CReOve7Ame4IzT3p7shlzTQY/KLaaMhRc4dWIJ0OZKLlDOBAe/1m4J8UjxP5KWC5qr5R6QkjJzC1IN/NtlgIci/C4uBOuuVeDyJffOuI1ox4simBJq0nnOZ2/wdn+ktkBUO8KCY27ge+aKOkqm+deCErjJcdgaUYM4Hbc8pOEZHjgfnA90tFyo+WlAfIu8n1pSuViZNGOEhaEgNoSQxg6C5D+fy4E/n8uBNZteAd389jxMVQK8rtslK21S/WGvEqLmlNh9ptpqr0Jfs8LUC7E8PRXma5jyUij4jIojxLgZDY+RGRZuAzwJ9cxdcAuwOTgDVAyaBz/eKuMfo70wBY9ZtnM2U7NraTUstzw6+b59bOjpKKvap3he9dZmD9mO6z17vf6c4Su1I/NC9PevF5GjREjftX/xmAJ197gZ/vd+F228vtshLXr6w73cmARFvWb9l5eNyxsb3gMdyhabJc/13lft8fCqGUJbLrVXVqwWOpHlxom4iUE8brCOB5VV3rOnZmXUSuA+4vZWy/bcEYDAZDJFAlnU55WqqknDBex5DTPWaLksPRWKHBitIvWjDLf/3PvOXOIPyghiEVH9s9oO7lSaw1UftutAHDB5TV2vBS17ReDJVy5C7/nfW/Wty/swH278ndqvHyey70W3WXB9UtrBDUGMwl5AnjJSK7ANer6nT7dRtwCPCNnP0vFZFJtskr8mzfjn4hMAMK3NSrERYHpxntlSAG/w2GKLG2dyVgOc80S/6HsXKiKZdyOOlMdWTOFwvsFkztT5M/jJcdFX+663UnsFOeeseVe85+ITDOeMS6jtVZM/b3+dVRADxz6n35dsti9HemccB/fQSA246+PlPe1lA4OZjzoxnYNiivK7PB0J9okmZPLYdS5BMX9/6xERYby005YnNzfKJfCIxzQ8+NFeBFWBxahgzgqWet2fIXj73EU1Il50tfyJXZYAiSM+aeDcCl+/880PMuuultAD41a0yg562WoHLUqCp9qfpMM90vBKYQTuiL3BnK7omZzvp7b7zLNT+z5iS5I8rm4jTPG6Qx7/wb04IxhEXQwuJw4NcnhXLeagluQmYwXWRh0K8FplDoC0cY1ve9zccutESocUBTUWFxKNZlBqYFY+h/RH0OVdjZNAMc5A+caH/yIeGM2bR3dbL6smcA6OzqKlj/4vmXACYXuSHajPjaZADevm5+5MKlVEs+kfAaUib0a2HHIqtH+pXA7PyNKbz9uwWZ1zsdbw3aa1rZeOtLefcpJiwAd6y4ldl3WfONXnprWZYDgMEQFXb9/sdJNFiJwIK+oc595xEA9h9ecA5g1eR7T6ELh0cCjEUWOBULjIiMBW4BdgbSwGxV/bWInA98HXDilZztJepmELjFxc0V552a9brhsLEAdBQ5lvM0mOzuY+j7rKjXj//jRWv6kcEQMd68/KnQzl1LYakLNLBYZIFTTQsmiRXs7HkRGQwsEJGH7W2/VNVfVG9ebdlwy4uA9XR34uVfzZR3/GVJ0f3aWlvZ5YfWHJqXz/+7iVxsiAX3rLTCSk1qn8RZf/tfANPijgCKkjReZNnYUTmdyJxbRGQxMNovw4JkSHt5yc46u7pY32c5CBhxMcSFGWM+n1l/4E67RWNa3KGjCinTRVYYERkH7A08A3wCDyGd7SigswDGjh3rhxkV8/KZD6No0TqDPrsHqQffyrxut73Bys29sim5MXYTwQz1x//96rKwTTBkqN9B/qpHwURkEHAXcKqqbsZjSGdVna2qU1V1avvwwlFPa4brgUGluLiA1W3W1dWZWRzKbcHs0DiMu96cw11vzvE1PfNXxpic6QbvHDxyOgePLJnx1hAAlptyytMSN0S19M214M4iTVghmx9U1SvybB8H3K+qexU7zuQpk3XevHkV21EOHSu3AjBo1EBoKFG5xvRu6sskPtuSei8rjI0XTvrrN/n9p68BINWVysrMaTCUy5bUe/zhJMsf51u3HJsp1z7l9A+fB2RnTjVk09batqBYKP1CfECa9bd4i1F4GKsqOkdYVCwwIiJYaTc3quqprvJRTtY0ETkN2EdVZxY7VpACI2q5ahZrtTh19rzkEDav3wLABSd/jdN/+iuAgi7NBkNUaWtt5ax/ngPAmdPOoCltZz11PZN0vt3FTuOtGIdmQnD5VCowU0V0vse6ArESmGrGYD4BHAe8JCIL7bKzgWPKDekcJF66w5w6i856KFO20/EfQRKyXd2jBs7MRATY51dHZeKbdaQ2Z0Lzu2cym5aGIQys+Vw/3laQ5yvYtnOrERaDr1TVReYXZbVgUtSma8vdvenz8d9Nrs9k13N3i/mJqPCPd/4OwNd/ezHLLvin7+cw9E+2vNnB4F2Lh0AymBZMPuIx1dXFhy4ub9KW091Vqk7Hmq10rNlaE/Fyp26thbiA1eo6aMRhHDTisEiIS1tra2a9d1P8ffxfuGdpZt2djjrKuO10fx6FyP2t3LHiVu5YcStDxpbnxm8wOMROYAwGg8EQD+LXRRYiO39jClA45EwlnmAGQ71zz2WPAzDjhweEbEltMV1k22NaMGXw9u8WFBQXgBGDdsmsd6c7Gf2daYz+zrSsOk5+mVripVvQYAiaR2d7vY0a6oVYCsxb89fy1vy1WWUdK7cy951HMpFbCyEqNbsBuz1wBiTaWPWbZ1n1m2ez6gQVWsZ5n//1u/9h+ROrWP7EqkDOm48wz10L4jIGEzX2//pHwjahJGlNZ0L/1xsi8nkReVlE0iJSsBUkIoeLyBIRWSYiZ7rKh4nIwyKy1P6/Y6lzxlJgxk4dydipI7PKBo0ZyP7DDy4ZuVVFPbkqxxnnPaoo//jG/7H7fqPZfb/wwsSN2Xen0M5dC7q6Oln8yAoWP7IibFM8c8eKW8M2gSZpDtuEkiQkEZsw/xWwCPhvYG6hCiLSAFwFHAFMxJp2MtHefCbwqKpOAB61Xxcl1lfyQz8zYcDjQLFWW74nRu2L/gPAhw4ex4cOHhe2GZ7xko21ViQaEiQaYn2rqQtUdbGqFg8VD9OAZar6mqr2AnOAGfa2GViT67H/f7bUOSORcOyF51/oaGttK/XG89J2UU27K9qB9bU8QQUYm7wTRbv6r03nlL1H3K7V+yo54AJ4UKzjemGAiLgHs2ar6uxKzluA0cBbrtcrgX3s9ZFOlBZVXSMiJePbREJggCVR9IwQkflRs8vY5J0o2mVs8k4U7aqFTap6uF/HEpFHsJJA5nKOqt7j5RB5yiruUoiKwBgMBoOhSlS12nGDlYA7f8oYYLW9vtaJNSkio4B1pQ5mOkYNBoPB4PAcMEFExotIMzATuNfedi9wgr1+AlCyRRQVgfGzD9FPomiXsck7UbTL2OSdKNoVRZs8ISJHi8hK4OPAX0XkQbt8FxF5AEBVk8ApwIPAYuBOVX3ZPsQlwCEishQ4xH5d/JxRmMlvMBgMhvojKi0Yg8FgMNQZRmAMBoPBUBNCF5hCYQkCtmGsiDwmIovtUArfs8vPF5FVIrLQXgJNYi4iK0TkJfvc8+2yssM1+GzTHq7rsVBENovIqUFfKxG5UUTWicgiV1nBayMiZ9nfsSUicliANl0mIq+KyIsicreIDLXLx4lIl+t6XVsLm4rYVfDzCvFa3eGyZ4WTyDCoa1XkPhDq9yrWqGpoC1b2leXAbkAz8G9gYgh2jAIm2+uDgf9ghUk4H/hBiNdnBdCeU3YpcKa9fibwvyF/fm9jTTAL9FoB+wOTgUWlro39Wf4baAHG29+5hoBsOhRotNf/12XTOHe9EK5V3s8rzGuVs/1y4Nwgr1WR+0Co36s4L2G3YIqFJQgMVV2jqs/b61uwvCfCC95VnLLDNdSQTwHLVfWNoE+sqnOBjTnFha7NDGCOqvao6uvAMqzvXs1tUtWH1PLMAXgaa15BoBS4VoUI7Vo5iIgAXwBu9/u8JWwqdB8I9XsVZ8IWmHxhCUK9sYvIOGBv4Bm76BS7e+PGoLujsGbQPiQiC0Rkll2WFa4BKBmuoYbMJPsmEOa1gsLXJirfs68Af3O9Hi8iL4jI4yLyyRDsyfd5ReFafRJYq6pLXWWBXquc+0DUv1eRJWyB8TUsQbWIyCDgLuBUVd0MXAPsDkwC1mA124PkE6o6GSuy6bdFZP+Az18QsSZhfQb4k10U9rUqRujfMxE5B0gCf7SL1gC7qurewOnAbSIyJECTCn1eoV8r4BiyH1wCvVZ57gMFq+YpM/M+XIQtMMXCEgSKiDRhfan+qKp/BlDVtaqaUtU0cB0BN39VdbX9fx1wt33+tWKFaUA8hmuoEUcAz6vqWtvGUK+VTaFrE+r3TEROAI4EvqR2573drbLBXl+A1X//gaBsKvJ5hX2tGrFCyt/hsjWwa5XvPkBEv1dxIGyBKRaWIDDsPt8bgMWqeoWrfJSr2tFY+RSCsmmgiAx21rEGixdRQbiGGpH1lBnmtXJR6NrcC8wUkRYRGQ9MAJ7Ns7/viMjhwI+Az6hqp6t8uFi5NxCR3WybXgvCJvuchT6v0K6VzcHAq6q60ikI6loVug8Qwe9VbAjbywCYjuWtsRwr4mcYNuyH1bR9EVhoL9OBPwAv2eX3AqMCtGk3LA+VfwMvO9cG2Akr2c9S+/+wEK5XG7AB2MFVFui1whK3NUAf1pPkV4tdG6xg8cuBJcARAdq0DKuf3vleXWvX/Zz9uf4beB44KuBrVfDzCuta2eU3ASfn1A3kWhW5D4T6vYrzYkLFGAwGg6EmhN1FZjAYDIY6xQiMwWAwGGqCERiDwWAw1AQjMAaDwWCoCUZgDAaDwVATjMAYDAaDoSYYgTEYDAZDTfj/pmvFdx5hV8AAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# must give diverging_cmap here; custom map is unknown to Contact Map Explorer\n",
    "diff.residue_contacts.plot(cmap=custom_cmap, diverging_cmap=True);"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}