03-WCPS_extension.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# 03-Extension to the WCS core: Web Coverage Processing Service"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "<a href=\"./02-WCS_core_requests.ipynb\"><< 02-WCS core requests</a>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space>&nbsp;<space> <a href=\"./04-WCS_access.ipynb\">04-How to develop and test WCS queries?>></a>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This tutorial is a practical introduction into the Web Coverage Processing Service.  \n",
    "Examples include:\n",
    "- [Basic WCPS request](#basic_wcps)\n",
    "- [Mathematical operations with WCPS](#maths_wcps)\n",
    "- [On-the-fly color-scheme application with WCPS](#colouring_wcps)\n",
    "    "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Web Coverage Processing Service (WCPS) requests"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The [previous tutorial](./02-WCS_core_requests.ipynb) explained the principal setup of WCS requests. A Web Coverage Processing Service (WCPS) is an extension of the WCS core suite. A WCPS has the same setup as a WCS and consists of three components:\n",
    "* Service endpoint,\n",
    "* request type and\n",
    "* query parameter.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    " \n",
    "> **Service endpoint with service description** remained unchanged:  \n",
    "<font color=\"blue\">http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1</font>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "> **Request type:**  \n",
    "The processing extension defines an additonal request type:\n",
    "- [<font color=\"blue\">ProcessCoverages</font>](#wcps_examples)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "> **Query parameter:**  \n",
    "Based on a text-based *SQL-like* query language, data can be accessed, modified and retrieved.   \n",
    "The principal setup of a query is:  \n",
    "\n",
    "\n",
    "   >> <font color=\"blue\">&query=for c in (</font><font color=\"green\"><i><b>coverageId</b></i></font><font color=\"blue\">) return encode (c</font><font color=\"orange\"><i><b>[...]</b></i></font><font color=\"blue\">, \"</font><font color=\"purple\"><i><b>format</b></i></font><font color=\"blue\">\")</font>,</p>\n",
    "\n",
    ">with:\n",
    "* <font color=\"green\"><i><b>coverageId: </b></i></font> <font>  specifies the parameter (coverage) the query shall be applied on.  \n",
    "     \n",
    "     \n",
    "     >> Example: <font color=\"green\"><i><b>temp2m</b></i></font> for ERA-interim 2m air temperature\n",
    "     \n",
    "     \n",
    "     > * <font color=\"orange\"><i><b>[...]</b></i></font>:    specifies the data or subset of data to work with. You can subset any axis of the coverage.\n",
    "\n",
    "\n",
    "     >> Example: the coverage temp2m has three axes: Lat, Long, ansi, and subsets of all three axes can be defined as: <font color=\"orange\"><i><b>[Lat(-10.0:10.0), Long(-10.0:10.0), ansi(\"1999-01-01T00:00\":\"1999-01-31T18:00\")]</b></i></font>.  \n",
    "\n",
    "\n",
    "> * <font color=\"purple\"><i><b>format</b></i></font>: the implied output format can be specified here.\n",
    "\n",
    "\n",
    "     >> Possible format encodings are:\n",
    "* <font color=\"purple\"><b>csv</b></font>\n",
    "* <font color=\"purple\"><b>netcdf</b></font> or\n",
    "* <font color=\"purple\"><b>png</b></font>  \n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## <a id='wcps_examples'></a>Examples: WCPS requests"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The following examples show a variety of WCPS queries and mathematical operations that can be applied to the data of a WCS Server. Examples base on ERA-interim 2m air temperature global fields, 6-hourly values from 1 January 1979 to 31 December 2015"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a id='basic_wcps'></a><font color=\"black\">Basic WCPS requests</font>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* Retrieve **2m air temperature in Kelvin for Italy for 15 January 1999 at 12 UTC** encoded as **csv**     \n",
    "\n",
    "(*Note: csv output is suitable for further data processing and integrating data request into processing routines in Python for example*) \n",
    "\n",
    "<p><font color=\"blue\"><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(c[Lat(35.4:47.1), Long(6.6:18.5),ansi(\"1999-01-15T12:00\")],\"csv\")' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(c[Lat(35.4:47.1), Long(6.6:18.5),ansi(\"1999-01-15T12:00\")],\"csv\")</a></font></p>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "* Retrieve **2m air temperature in Kelvin for Europe for 10 January 1999 at 12 UTC** encoded as **netCDF**  \n",
    "\n",
    "<p><font color=\"blue\"><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(c[Lat(-35.0:65.0), Long(-35.0:80.0),ansi(\"1999-01-12T12:00\")],\"netcdf\")' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(c[Lat(-35.0:65.0), Long(35.0:80.0),ansi(\"1999-01-12T12:00\")],\"netcdf\")</a></font></p>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "  \n",
    "**Note**: Both queries above do not require WCPS per se. Same data can be rerieved with WCS core requests."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a id='maths_wcps'></a><font color=\"black\">Conduct mathematical operations with WCPS</font>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Get the **minimum, maximum and mean temperature in degree celsius** of Berlin for the year 2003\n",
    "\n",
    "There are already predefined <font color=\"purple\"><b>operators</b></font> for retrieving the minimum, maximum and mean: min(), max() and avg(). We additionally convert the 2m air temperature in K into degC via a <font color=\"green\"><b>simple subtraction</b></font>.\n",
    "\n",
    "* **Minimum**  \n",
    "\n",
    ">for c in (temp2m) return encode (<font color=\"purple\"><b>min (</b></font> c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")] <font color=\"purple\"><b>)</b></font><font color=\"green\"> <b>- 273.15</b></font>,\"csv\")\n",
    "<p><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode( min(c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")]) - 273.15,\"csv\")' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(min (c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")]) - 273.15,\"csv\")</a></p>\n",
    "\n",
    "* **Maximum**\n",
    "   \n",
    ">for c in (temp2m) return encode (<font color=\"purple\"><b>max (</b></font> c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")] <font color=\"purple\"><b>)</b></font><font color=\"green\"><b>- 273.15</b></font>,\"csv\")\n",
    "<p><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode( max(c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")]) - 273.15,\"csv\")' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(max (c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")]) - 273.15,\"csv\")</a></p>\n",
    "\n",
    "* **Average**\n",
    "   \n",
    ">for c in (temp2m) return encode (<font color=\"purple\"><b>avg (</b></font> c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")] <font color=\"purple\"><b>)</b></font><font color=\"green\"> <b>- 273.15</b></font>,\"csv\")\n",
    "<p><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode( avg(c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")]) - 273.15,\"csv\")' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode(avg (c[Lat(52), Long(13),ansi(\"2003-01-01T00:00\":\"2003-12-31T18:00\")]) - 273.15,\"csv\")</a></p>\n",
    "\n",
    "\n",
    "With three simple queries, we retrieve the minimum, mamximum and average 2m air temperature of Berlin for the year 2003, which are -16 degC, 35.5 degC and 10.1 degC respectively.\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### <a id='colouring_wcps'></a><font color=\"black\">Apply a colour scheme to a png output</font>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The switch statement allows to classify the png pixel values and to apply RGB colour codes to it. Below, there are two examples how predefined RGB colour codes for 2m air temperature and total precipitation can be applied on-the-fly.\n",
    "\n",
    "* **Colour scale for 2m air temperature on 18 June 2012 at 12 UTC in degrees Celsius**\n",
    "\n",
    "> <font color=\"blue\">for c in (temp2m) return encode (  \n",
    "**switch**  \n",
    "case -48 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128}  \n",
    "case -44 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128}  \n",
    "...   \n",
    "default return {red: 255; green: 0; blue:0},\"png\")</font>   \n",
    "<p><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for%20c%20in%20(temp2m)%20return%20encode%20(switch%20case%20-48%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20128%7D%20case%20-44%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20128%7D%20case%20-40%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20128%7D%20case%20-36%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20128%7D%20case%20-32%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20128%7D%20case%20-28%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20217%7D%20case%20-24%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%2064%3B%20green%3A%200%3B%20blue%3A%20255%7D%20case%20-20%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20128%3B%20green%3A%200%3B%20blue%3A%20255%7D%20case%20-16%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%20128%3B%20blue%3A%20255%7D%20case%20-12%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%20255%3B%20blue%3A%20255%7D%20case%20-8%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%200%3B%20green%3A%20255%3B%20blue%3A%20128%7D%20case%20-4%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20128%3B%20green%3A%20255%3B%20blue%3A%200%7D%20case%20-0%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20218%3B%20green%3A%20255%3B%20blue%3A%200%7D%20case%204%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%20255%3B%20blue%3A%200%7D%20case%208%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%20245%3B%20blue%3A%200%7D%20case%2012%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%20218%3B%20blue%3A%200%7D%20case%2016%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%20176%3B%20blue%3A%200%7D%20case%2020%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%20164%3B%20blue%3A%200%7D%20case%2024%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%2079%3B%20blue%3A%200%7D%20case%2028%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%2037%3B%20blue%3A%200%7D%20case%2032%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%2010%3B%20blue%3A%200%7D%20case%2036%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%200%7D%20case%2040%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%20255%7D%20case%2044%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%20255%7D%20case%2048%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%20255%7D%20case%2052%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%20255%7D%20case%2056%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20-%20273.5%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%20255%7D%20%20default%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A0%7D%2C%22png%22)' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (temp2m) return encode (switch case -48 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128} case -44 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128} case -40 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128} case -36 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128} case -32 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 128} case -28 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 0; blue: 217} case -24 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 64; green: 0; blue: 255} case -20 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 128; green: 0; blue: 255} case -16 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 128; blue: 255} case -12 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 255; blue: 255} case -8 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 0; green: 255; blue: 128} case -4 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 128; green: 255; blue: 0} case -0 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 218; green: 255; blue: 0} case 4 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 255; blue: 0} case 8 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 245; blue: 0} case 12 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 218; blue: 0} case 16 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 176; blue: 0} case 20 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 164; blue: 0} case 24 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 79; blue: 0} case 28 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 37; blue: 0} case 32 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 10; blue: 0} case 36 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 0; blue: 0} case 40 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 0; blue: 255} case 44 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 0; blue: 255} case 48 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 0; blue: 255} case 52 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 0; blue: 255} case 56 > c[ansi(\"2012-06-18T12:00\")] - 273.5 return {red: 255; green: 0; blue: 255}  default return {red: 255; green: 0; blue:0},\"png\")</a></p>\n",
    "\n",
    "\n",
    "* **Total precipitation for 18 June 2016 at 12 UTC in mm of water**  \n",
    "By specifying alpha parameter of the RGB colour code, one can set specific values to transparent. In the case of total precipitation, all values with less than 0.5 mm of water are set to transparent.\n",
    "\n",
    "> <font color=\"blue\">for c in (precipitation) return encode (  \n",
    "**switch**  \n",
    "case 0.5 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 255; green: 255; blue: 255; alpha: 0}  \n",
    "case 2 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 0; green: 255; blue: 255; alpha: 255}  \n",
    "...   \n",
    "default return {red: 255; green: 0; blue:0; alpha: 255},\"png\")</font>\n",
    "<p><a href='http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for%20c%20in%20(precipitation)%20return%20encode%20(%20switch%20case%200.5%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%20255%3B%20green%3A%20255%3B%20blue%3A%20255%3B%20alpha%3A%200%7D%20case%202%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%200%3B%20green%3A%20255%3B%20blue%3A%20255%3B%20alpha%3A%20255%7D%20case%204%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%200%3B%20green%3A%20128%3B%20blue%3A%20255%3B%20alpha%3A%20255%7D%20case%2010%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%200%3B%20green%3A%200%3B%20blue%3A%20255%3B%20alpha%3A%20255%7D%20case%2025%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%20218%3B%20green%3A%200%3B%20blue%3A%20255%3B%20alpha%3A%20255%7D%20case%2050%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%20255%3B%20alpha%3A%20255%7D%20case%20100%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%20255%3B%20green%3A%20164%3B%20blue%3A%200%3B%20alpha%3A%20255%7D%20case%20250%20%3E%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A%200%3B%20alpha%3A%20255%7D%20case%20250%20%3C%3D%20c%5Bansi(%222012-06-18T12%3A00%22)%5D%20*%201000%20return%20%7Bred%3A%2051%3B%20green%3A%2051%3B%20blue%3A%2051%3B%20alpha%3A%20255%7D%20default%20return%20%7Bred%3A%20255%3B%20green%3A%200%3B%20blue%3A0%3B%20alpha%3A%20255%7D%2C%22png%22)' target=\"_blank\"> http://earthserver.ecmwf.int/rasdaman/ows?service=WCS&version=2.0.1&request=ProcessCoverages&query=for c in (precipitation) return encode ( switch case 0.5 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 255; green: 255; blue: 255; alpha: 0} case 2 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 0; green: 255; blue: 255; alpha: 255} case 4 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 0; green: 128; blue: 255; alpha: 255} case 10 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 0; green: 0; blue: 255; alpha: 255} case 25 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 218; green: 0; blue: 255; alpha: 255} case 50 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 255; green: 0; blue: 255; alpha: 255} case 100 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 255; green: 164; blue: 0; alpha: 255} case 250 > c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 255; green: 0; blue: 0; alpha: 255} case 250 <= c[ansi(\"2012-06-18T12:00\")] * 1000 return {red: 51; green: 51; blue: 51; alpha: 255} default return {red: 255; green: 0; blue:0; alpha: 255},\"png\")</a></p>"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The examples above shall give you an idea of what data processing and retrieval with a Web Coverage Processing Service is alike. WCPS can deal with far more complex data processing requests. For more information about the variety of data processing possible with WCPS, have a look to the two WCPS webinars available (see [Ressources](#wcps_ressources))."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The [next tutorial chapter](./04-WCS_access.ipynb) will show you different ways to access and work with a Web Coverage (Processing) Service."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## <a id='wcps_ressources'></a>Ressources"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### WCPS Webinar Part 1: Introduction"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEAAUDBAgICAgICAgICAgGBwgIBwcHBwgICAgICggICAgI\nCAgIChALCAgOCggIDRUNDhERExMTCAsWGBYSGBASExIBBQUFBwYHDwgIDxcUEhQXHR8dHB8cHBce\nHB8bHBweGB4XHBgeGBwcHh0dHRUXFRwfFh4UFxcfHBweHBQcHRwVHP/AABEIAWgB4AMBIgACEQED\nEQH/xAAdAAEAAgIDAQEAAAAAAAAAAAAABwgFBgEDBAIJ/8QAThAAAQQBAwICBgUHCQMLBQAAAgAB\nAwQFBhESEyEHMRQYIlSU1RVBUVNVCCMyYXGV0RYzQkOBkZPT8CShsTVEUmNydZK1wdThFzZiZPH/\nxAAaAQEAAwEBAQAAAAAAAAAAAAAAAQIEBQMG/8QALREBAAICAAQEBQMFAAAAAAAAAAECAxEEEhNR\nBSExkWFxscHwMkHRFFJikqH/2gAMAwEAAhEDEQA/AKZIiICIiAiIgIiICIiAiIgIiICIiAiIgIiI\nCIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiI\ngIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICI\niAiIgIiICIiAiIgIrKP+Rjqj3/AfGZH5euR/Iw1Q/wDz/AfGZH5egrUisr6l+qff8B8Zkfl659S7\nVPv+A+MyPy5BWlFZQvyMNUN/z/AfF5D5enqY6o9/wHxeR+XoK1orJ+plqj3/AAHxeQ+Xrn1MtUe/\n4D4vI/L0Fa0Vlw/It1S/dshp/wCMyPy5fXqVap/ENP8AxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr\n3/T/AMZkflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDGZH5cgrMisz6lWqvf9P8AxmR+XJ6lWqvf\n9P8AxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr3/T/AMZkflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0\n/wDGZH5cgrMisz6lWqvf9P8AxmR+XJ6lWqvf9P8AxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr3/T/\nAMZkflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDGZH5cgrMisz6lWqvf9P8AxmR+XJ6lWqvf9P8A\nxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr3/T/AMZkflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDG\nZH5cgrMisz6lWqvf9P8AxmR+XJ6lWqvf9P8AxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr3/T/AMZk\nflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDGZH5cgrMisz6lWqvf9P8AxmR+XJ6lWqvf9P8AxmR+\nXIKzIrM+pVqr3/T/AMZkflyepVqr3/T/AMZkflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDGZH5c\ngrMisz6lWqvf9P8AxmR+XJ6lWqvf9P8AxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr3/T/AMZkflyC\nsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDGZH5cgrMisz6lWqvf9P8AxmR+XJ6lWqvf9P8AxmR+XIKz\nIrM+pVqr3/T/AMZkflyepVqr3/T/AMZkflyCsyKzPqVaq9/0/wDGZH5cnqVaq9/0/wDGZH5cgrMi\nsz6lWqvf9P8AxmR+XJ6lWqvf9P8AxmR+XIKzIrM+pVqr3/T/AMZkflyepVqr3/T/AMZkflyC97d9\n2XIPtsy4It/Jci396Du8v7f+K4QW/vZN9/JBwbbt+tl0uz7frXo2/wB663bbf/cg+f4d/wBq4YfN\nvtX1/rsvsQJ++3kg+qnbt/x/3r0LpGIm+tvPdd+yAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIi\nICIiAiIgIiICIiAiIgIiIMFqvVtHFvWG5JM0mQlKGpBVo3b9iYwB5JOFehDJLwEG3I+PEd23dt23\nxNrxQwsMk8diezUerHblOW9iMrTrmFSM5bL17VqqEFtxjjM2GEzchAiFiZt17NTYFp8phsgVkImw\npZAngIORWGs1PRvZPm3T4O7E/slvuzdt91FlDwhmE4Xs5fGZBqte7FGd2hkrMl4rVOWm55cpswcV\nkSjskZDCEDmfHiQNuKCRD8U8OIPJI+UibqxwgFjTefglnmkaQwhqwTUBlty8IZjcIRNxCIydmEXd\nYzTnitVlht2bUV8YWy1ypSGrp/OWJ2grR1tyuV4aZzVZHKUn2mCPduzM/F3Ws1fDC/FIF8c5jgnp\nXqVmhBFjcjJixngrZerMc1axmZZpp5Ycs7O8U0e3ocO7F3S74Y3ZTKV83jZbFia3Pdimx+THH9a3\nKBBJUr1M1EcUnCBgd5ZJndhZ24t2cNxv+MWnoWYnt2Jo3ap/tFLD5e9WYrkcEtOJ7NKnJCM8oWqr\njE5c39Ii7bmO/ti8TMScB2ALISDFa9EmgiwGbkuQT9ALLDYx4UntQA8MsRsZxiLtIGzvyZanp3ws\nho4ubEDlAPfKacujKcDMYjhoMBF0ijabuczYIn3Z24+keRcPa69WeG0tvI3MhVzNUGv3vSTo24bs\ntN+eLxGPiKQaGTrHLZD6KMgkcuLDdmHi/wCkg2W34vYCIGkK3ZOP0QbkskGHy1gKtZynDqXzgqE2\nOcSq2WIbHTIehJyZuLrIal1BPBlcBVheMq2YlyA2X48ncIaB2YXjNn2H2xbv33Z1pGB8Lo62LzeO\nfK1y/lDp9sT1QruI1pSbMFJOIHZMpYuWV9mMj5ca/czcnJpEtYgfScSZOzljWnEX225OVR4Sdm+r\nssnG5+hji0/3Uj/a9a/dq4SccZJ543HLb35Z17TqXn1tm5qrxRw7C8guRSOLF5Pswiz9v2/tZR/4\nh+L+SxI0uhhiyfpgWOqcR2A6ZQnEzbjDXkbY2lbbu3eM1tvjoWQiw1i3jI4ZbWOdrLwT1+u01cWd\nrICLOxsTA/Ubg7O7wM3ffZ/J+Tpq98zga9mU4itQzWK9wYQCMAkGYjjZowZmFngkgJv+33d33d+v\njrWMfUmsTETpzb2mb8kTqdIou/lOZKBmefTXQY3dhee5ZhYnbzYXkpNydt28vtU9eGupHzGJo5Io\nWgK/B1SgY3kaN+ZA4sbiPJtx89m81CH5adv0k8Hh6wPPfs2ZJ468e3PY2apXF2f9FpZDkZnft/s5\n7+SmzStGLB4vGY8ic2p1Y6ryizMJHDVkmmmLk7cRLoSl+0mV+Jri6Nb1rqZ+Mz5KYZydS1bTuIbK\niw9fUMMrRlEEsgzkAROLRsxyHXC00bczbYmhJzd32b82Tb77M/RS1ZVmcODS8ZDrCBkIgDtY59E3\nIibiz8Nti2J3MBZncmZYGtn0REHgAP1eX61y7tt9jsm+32s6+HfdB9Rn32XduusAX3/xZA2/+E7e\na4cv9b/w818uTf6b/wBXQfe7f+q7IXffy/aukZB+vd/9fqXYE7v2EfLz/V+tB6EXDOuUBERAREQE\nREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQdU9cJP0x39lx83Z9ncSdmdvLuAv\n/YvP9FV9tum22wN+kXZgbYGbv2Zuzt+tmfz7r2og870ouIhwbiDu4szu22+7l5P3Z933byddTYmu\n39W3mz/pH3dm2bfv32bt3XtRB55aURb8gZ933d93Z9/a2fdn7P7Z/wDiXW+Mg4iPTHiG3Ft32Zm3\nZm8+7d3XsRB4oMXCDs/FyIXZ+Zk7l7P6G+3bYfJu3bv9rrryFgRs0wd/amKdgb7eMLk/+5ZFapqk\n9svgG325TZHdvt2oG6weI4erhiv+dJ9r1n7K3y9KOb5R7zr7trUMan0lldOWXyWk6VWenaOU8tgm\nHpy2JpC3CxDYMn4xR9mGAOAxs8vES6hcZnRdPHlnHPeO3dF6RdFHhL4Z2YLsuodQShbz95ncRDvW\nxkRDw6FZu7PI0b8HNuzDyEXLkZyb/qopBCE44BneOWUnYoTm6btTtPEYgDs/J5ekH19pX8vNswiZ\nMtslt2/PkUpFI1DTijs8yAKFRiOWZgmKkfSL2LLxHJsbPGXVq1xIu7ELxExe0wh8WJJxfnDjoWcw\npywi+OkaRnkuEPCVxJhGSKvxN93biW7+S3RF5rtWPMZBxlca7M8dazNHvTtEMhhweGFmIgNpCZzZ\nx4+bNxcu6zWNmsFLOEzBxhcGjkCKSNpeQCbk3MnbZndx2Z38vNvJe9EGO23XcIf/AMXpYGZ99u64\nKNndB0OX2f6/+V8E/wDr+K7JjYW9lt9vN/qb+xeeRvr+p9kH0Tt9vf7FzxZtn27fWvnZmdvsf61x\n9vfuyDvCLfdtu32rsiiYe7d3+t1xWdtv+K7kBcr5b/W65FByiIgIiICIiAiIgIiICIiAiIgIiICI\niAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgI7IiDzTD3/UTf711CPbi/\n177L1Tju37O7LqcHdvLv57IOoR7bP5t5Lt6T9n/sdc9J+z/Wu5v9fag4Ftv9f8GXP+u65ZlyzIOP\n9d1zsiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIi\nICIiAiIgIiICIiAiIg4dH/19q5RmQcbf2LlmREBERAREQEREBERAREQEREBERAREQEREBERAREQE\nREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREHBOzM7u+zM27u/\nZmb7XUS3Pyg8AEptEOStVIDMJ8rTx5y4+Fx2Z3KZyYzB992IAJibZ2d2cXfMflJZGWrpXMywsTmV\nUK78N+TR2bENWYuzbsIxTG7v9TM/l5rY/D/D1qmIx9OCMGrx0IB4cWcT5RC8hm39MjIiInf9JzJ3\n33QZHGZqpZqBfgsRSU5YeuNpjZoukzO5GRFtwYdi35bceLs+2zqMpvyhtPjI+wZM6QydMsxHjZCx\nwnzYOPU36pPu/bjG++3bzbeJMrfmo6O1vRq8o62P1bbx9Xj+jBSO9AB1x7fzbsLi7f8A7T/ap18S\nMrjcHiYKEuJv38fbiLGjSxFSOy4QdAmfqicocQcW25M7vyJv2oN1ky9Uaj3ysQtSat6W9x5R9H9G\n6fV6/V349Lp+1y3227qMY/yhdPvKLEOSipyn04svNjpAxxnyIdhkd+q36Lu7vGzM2++2z7aPrjPY\n69ovE18RXvVcTkdQUsQcNuRynirBYllkeSYbErtG8kAML9TtzBu3ZlP2b07SuUJcZYrgVGev6MVY\nN4hGFhYQGJ4nZ4XDiLiQOzi4C7OzsyD13Lox15LIsU4RwHOI1+JnMIg8jNDuTCREzdt3Zn3butVr\neJWPk09/KUQsvR9GOx0eMD2/YlKAoOLTdHr9UXDbqcd/rW04fGwU60FStG0VanBHXrwi5OMcMYNH\nHGzk7vxERZu7/UqsWt4sFc0js7cvECvhqwC+7hQntDkK0r7u7iPKA33Zn25Dv9boJ+xPiTQs6fk1\nGEdoaUFe5OVeQIBt/wCyyTRHEwNM8XVM4XYWeRmfmG7tv22PS+XDIUad+OOWGPI1ILUUVgQGYI5o\nxlAZRjMhGRhNt2YnZn+tVezchV8JntKxv0pLevIMVQCN9ujTvTx3Kg+bO4kNaRn37P1OO7st88Zt\naTV8tU09UzdbTNSti2uWsrPBHORu5vBWx8ISOzCXEOo+xC7t35MwOJhKviPq6vgsZZytqKxNBSeB\npIqgxFOXWsw1Q4DNLGD7HMLvuTdmfbd9mfN0bAzRRyizsM8YSCxbcmExYmYtndt9n+p1V7Vuupcr\novVFK1erZOxgr2KibJ1BCMbtOfKUyrWSij9mM3eKcfZZm/NfW7OT7FmspqjT1XG5+7l69nHSyUIs\njgwoRwwUqkwMP5i1yKWaSPdvbfi7k7O/IW4oLCqMNU+OWEo2p6ohfyB0ZHjvy4un6RXpEPJpPSJi\nMRbg4ExcOWzs7P3EmbddeZCSpisnbhEimpYy7YhEG3MpIq0kgCLbPuTkLMzLVPya8bDW0vieiLM9\nqv6VOfmUs8xkckkhP3MvId38hAW8mZkG36R1HSy1OK9QnGxWsM/A2YhdiF+JhIBsxRyC7OziTM7L\nLKIvEQ4tJ4+Opp2tDVvasz8dau8hHJBFdukwnbKOVyYQAQjjEGbgH5puLiHF8fWv6g01lsRBl8yG\nbx2pLT0XkkpQ0pqV4+9fpNGRdQJDIQ2ctvPYRdm5BNqKE8Dr+1i62s6+TslZtaUsy2KUljixSU7c\nPPEwkTMzGTmwi7t5dYWf6t8ZqbWmVxGndN1bGVGvl9SOx3MzkxiJ6Fd+NmzI8TB0ykhGxXgZnHb7\nX3diQT+igDwu8QbXp+Rwsmfr6ij+hJslQzFeCKvJDNG7xS0pAjImMm3GVndydm+vYmEPb+Tyepcl\nj6mbtZtrMJ1chBBipqsYDYnC1Yjis3LwM5jtKBBxjj7BHE+7vyYgnJFVy1rnOYu7ipLGrsdl7V3M\nV6uU0/j69WavWqySmEjhaj2JuLcBHcYzYpA5ObCXLeM1kdQZ/UOWxmIzI4OppeOmMswUIbstu1aj\nKUXIJ3ZmjDpyhtvtvHu4nybgE1oq4YDxQzFDCauvX5gt5DEZyTH1ItmetDYKYa/TgHYSOrHIRmwk\n/IhiZndnfdfOi/EW3Ry2Jgsauoalr52y1K5BFXrwHQtSiI1pIJIdurGU7jF3YBdi7Rs5NxCe9Y5+\nLF0LeRnCWSHHwHPLHXYCmMRbuMbSGIOX/aJm/Wovr/lGYjiEtjG5+lVkYCa9axkb1uJ8eBMVaxIR\nsTEztwEt/q3d2Z9u8ff/ALYzf/ds3/BlDtTxKK3pSHT9PT2oLVybT1bExztjg9BeQqUdT0sLAyk7\nQM79UTIBbsLu4NuTBZWlZjmjjmhMZIp4wliljJiCSMxYgMCbsQuLs7O31Ou1V4zUWex1vRmmqOU9\nBmtYWzWtydILMATQ0ic5RjNm6rwiMjxtuLOUcfLtusnas6jv5NtL0c81c9P4yCbN558fBJZt252Y\n4YQq8mGGPpHGXsmxee5Ft7YToihjSXiJkMfDqijmpIr+Q0XANkLcQDXbIV5qxWaonGDcY5v5sSdm\n2brg3tOJEXg0lW1xNXxWeDMVr4Zb0O1awJUa9WtFj7AtM7Q2nLm0wxGG3ZnYi9opWF2MJ2RV68Sf\nEC1Yz2QxsWp6mlqeBjrCMs9aGzLkrksbTSbdZ24wxNsDsz7M79xPm3Dw5fxjyculad6G1DBYHPNh\ns3mqlULcdSMWcvTY6xNw5SxHWNmJhZ3l4jweSPYLJooT8CdRZGxYtiOo6OpsfHR6sfOA8bl4LfPY\nY3pHHyaobc2eSQ39po+LMzFvGuk/EnOXzCw2rqVTLvaeM9M5jFFRxwbWnj9FPI9M+JcWZuPFpmd3\nDlybk4WNsa6qBn4tOvHY9NnxX0oMzNB6M1fry1+Du83W6vKEn7RuOzt7W/ZtpUW3dQXh13UxnpBj\nRk0t6XJSHi8T2vTLcbzcnFjcuIAPmzbC3ZRpi8zqu9idSZKDUD1IdM5XMGAPRgsWbIVa8Vj0XqyN\nwr144xFw9g3Mp5OfYRQWdRQhqTxDytrEaUgx8sNXMazaIZLbxCQVgihjK/Yghk5D+nIGzFy2Enb9\nLYm8dC3qXF6s07h8jm3ydO9FkZXmGpFUKyIULZtDZiFz7xSxREJsW5c+7vtsgnHL5GCpBNasyjDX\nqxHLPNI+wxxgzkRP+xmUXVPyhMCcoNLHk6tScxCHLW8eUWOkcndm2m5OYD23cjAWFmd3dmZ3bcvF\nrSx5rC5DGRy9GS7CLQyORCLSxyxzxDIQe00JHEIlsz+yRdn8nhbM66u1cY2F1ppm1HjCihqT5XGP\n1K7CBxBBK4w7xxEzjGTcJebuz8I924ILJM/+9FAHilrU4L+O0/Q1BX0/Qr4WK5Nm7IR25LAvtBTq\nRPK7bm4A0rkxC5M+/LYXE8Pa8VchJpTUjR5aGzktNXsbFDm8eEIPapWshWCGw8PF4wMwa1G+w7ex\n23dnJwsuir7qizqfDY2pqK7nDtwtfxFnIYyClFFXq46beO3WA25S2i5S1RE34P7JE7O5Ot0sZ+3d\n1lWx1S0YY7CYg72UCF43jtWLTtHUrTE4uTcYyjnHi7b99927OGy+G+uqmeivS047EY4rK2cXO1lo\nGcrFcITkOPoTSN0nacNuTiXZ9xbtvtKrxg/FLIUdMapytiZ7lrG6muUMc1hhcYgcaAQA4xsPOKJ5\npZNnfd+Ltv37bdpfA6zx9+hLazEObpXCNsxWkp16foO8bcJKkgPym4yP5MIM4g/sM5MQBLKKsOpP\nEu/du5ieHVlHT8eGu2amLw89WvMeQ9GZ2KzNLL3GOeT2RfY2bZ2YWcHI538KdVNm8Nj8psAldr7z\nBE+8YWIzOCyAO7u/EZopW2d3dttn7sg2dERAREQEREBERAREQEREGM1Zg4clRt4+xv0chWlryOG3\nMWMXFjB3Z2aQXdiZ3bzFlD+Ei1/iqkeGgpYnIBVjevRzh2+kMVYGYK72qpm0pyiGzbCJM3Fmd5Nu\nZzmiCMtPeEkEem7eCuzvYmzJTWsnfYWcpMjKQStaBi/SeKSKDjy/S6DO7Nu7LAUpvEGlXjxcdDE3\nSgiaCDPyXeMfSBmCOSzVM2mOxwbd3YXbdu/L+lNiIIgwngzw0g+m7FsfSHI7Ed6ADYa1zrekQFEz\nkxmAEwg5bgRjz24cm44rN1/EDJUyws9TF02sRDXu6hhu8xkhf2ZygrA7TBMYM++8YM/N2bp77hOi\nIPDp/HNTqVqjSzT+h1ooOvZkKWebpxiHVmkN3I5S47u7/W7qHMv4Z3j19XzAV3fEnHDduWPSIWD0\n6GlYpQxNXeXqvIPGufJg47SH33clOKIIMy/hlel19BmBru+JIYLtqx6RCwenQUpqkEbV+r1nkEgg\nPlw47SH7W7kyy3ilo/Kx5qDUWGqU8nL9Gljb+JulFE00TSvPFNXlldoxm5cRdyJtmjBm3Yi2l1EE\nFax0fqO/pTM17NSieTylulJTxWMapXGCrDkas7wyW5SjjmlEBmLcj22Ftndy2XRlNI6szEWOwOVr\nY6LE0Zqct/KxTtJJkK9f9Gt6I5O8VghEee49PkTOJEIuBT4iD4sQhIBxmLEEoEBgXdiAmcSF2+x2\nd2UIYPEay03CWJxVOhmcbHNL9FXLVxoLFKCUjlaO7HIYPOIGRN+bd3dn7cWdgjnJEEHT+EOVLAwA\neSafUNPODqKGxPJLLSHIMTk9WNjHeOq7vz9kBZ5Ny4szuK9uL0/qTOZXGXNQU6eLpacmksw061gL\nUt69twhsOQETQ1w4sYs5MTeTs/LeOZEQQV44eF1/K5/HWKkZljcpFXqak2mrxxNVqXYrcZSxSG0k\n8kjO4i4CfF6wb8Wd+W3eN+irmSHGXsY1Uslp669itBeFnr2YJBEbNUy/oMfTif7H6e24u7GMjogi\nnRWN1HYO/YyePxmKryY2StSxVNq01krJMTFPLdj9kAdn48GJ2fcd2Dg7yefw10Pko9CfQVkPo/JT\nUstXcDmikaE7Ny7JDzlqnIDgUc0bu4O7sxv9bbNLyIKtSeG+pSx+KqNp7FVBwOTx080lGxTbIZUo\nH6J2jkeQYuHTOSQ2kNjIuHEW24rbsgGUx+rM1JpwMfkpcvBQLL4y5YenLQmigYa9tiJm61Nxm5O8\nXUJzmkHYXFlO60bxE8JsFn5gsZKn1LEQNGNiKeaA3jYnJgk6RMMjM7ls5M7tyLZ23dBDnhvombPY\nXWuOlsRPZs6ns7X4x415L9eWOwcgC3Io68p7eW7sFjtu7LetF0NWWL9L6RxmGxdTHu5XrNcK1ibJ\nGwC8ZVRDk9Ruo27uXF29rbfYVJWj9M0cRUjo46uNatE5EMYkZkRk+5HJJITnIb9vaJ3fZmbyZmbM\nINU8YMTYvYHK06kTz2bdGWKvCxxxvJI7eyLHMQgO/wBpEzfrXp8McdPTweGqWY3is0sPjq9mFyA3\niniqQxyxucREBuJiTbi7s+3Z3butiRBGWtdM37Gr9MZOGs8lDFwZQb1ppq4tAU1KxFCzxHI0snIz\nFvzYltv32buvFrnTecx+cl1Dp+vWyD5GjHUymJsTBWOU4Xb0ezXnkcQY2ARB2Im2YX7FyZ45aRBE\nmkPDS5Yr6js504Y8lrSL0ewFL85FQrR1irVRiI+0kwCQk792d4o+77O74PT2C1v0cZgpBqY6hhZa\nsc2do3GeW5QrM8YVoIGd5AMohAXeQBZ34v8AU7PO6IIV1bpDN43OZHL4bG4/M19Qx1fS6N2SGCSp\nZrh0RmiOZ2F4XEjMmZ3InM229kd8pb07qocDXGvcx0eajuhbuVoaleLHXa27seMkd4XdmePp7yDw\ncnBx5iL82lZEEJ+E+g8k2emz17F0NPj9GSUY8ZjDgN55pLPVkuzvX3h5Owt9rltE7szg/LVtYeHe\nrMtCWHyNHD3jOZ+GtZWqhcjq+lBLxevCASjK8Y8HjAGBxbhyfZpHsqiCLZNJ3x1tSyYwGeMq6WHH\nneOeB39KG3aPpFE8nXI3CQC5sHH2vPfdlgtD6Iy1fTWs6E9Io7eavaglxsD2KhPZjtY6GCqTSBM8\ncXOQSHaUgdttyYW7qb0QV51Xpy1QwejDeSrRz+CswRUauQMHhtTygwz4z0qEiiCSTpx7PzYS6ZMx\nM+xj1QSZu1rjS9jMQVak40snK2LqTjYKjW9BtxdW1MxOMk0s8rDtHuLNEGzu/J2nHWmlaGZqFSyV\ncbNcyE+BEYEEg7sMkckZMcZsxE24u3YiZ+zuz4Lw68KMFgJZJ8ZT6ViePpHYlmlmk6XJjeMOqTjE\nLkIu/BmcuAct+LbB6/F7RceoMNdxMhCHpgA8UhhzGOeKUJ4SIfrDnGLEzd+JFsoyztPXuVoHgreP\nxdeOzFHXuZ57gyhLCzj1ZIqoF1GmIR77gLO5Ptw3ZxnlEEJ6v8OsjjLuNyuAqVcp9H4KHB2sVkJI\nwKetX4lXnhmnfg1jcAFyMuzRjtvyLby6z0dqS/pPNVrFSiWTytqlJSxWMapXGvWhyFWfoy2pSjjm\nlEBmLcj22Ftndy2U7Ig1vN6ca/g5sVP7D3MW9Q39knikKv02NvMXIJNiZ+7bg3mtG/Jk0XksXTuW\nc1EcWUyc1cJBlngsSjTp1Y61MJJq8hgZM3V/pO/Fw377qXUQQHp3wnv29Oanw9+L0GfL6juX8dJJ\nLBMJRs1GSrMT1jk6cUh1zjJn2NhIn477b57TMGtMhkMcWWjr4ijiCMrr0rgylmZOmLB+ajImir82\n3cT4v+m2z7s7S8iCA7GjNSYa3locNjcXlKWdyFi/VtXTgjlxc9kfzjSxy7dauDgDCAMXYRd+5EzT\nPpCjYrUasNs4ZbUcI+kyVYRggKZ9yleGIWZhj5E+3Znfzdt3dZVEBERAREQEREBERAREQEREBERA\nREQEREBERAREQEREBERAREQEREBERBE/jRlMz9NaZxOIyr4n6dbNNZsNj6l9t6dWtahd4rQbv2aY\nPZMP57d+XFmXRNkNQ6dt4dstmoM7RzuXixBu+Hix1qtYsxylUkiepIQSx84nYuQtsz/2j4vHqjcs\n6l0XBj730bclbUfo+Q9EiuejuNCqcj+jTO0cvOIZI/aft1eTdxZbJgPDi+9unczuoLGcPEzHPQrt\njqeNqRWCjKIbEkNXd55QAzYXIvZ5Psg2jU+t8Pi5I4sjk6VKWZmKOKzaijkIHdxaTgRcmi3Em5u3\nHdn7r6yussTVYCsZGnCM1QrsUkliNo5KgvGxWAk34lFvLEzEz93kDbfdlHYHBeyeclw2Ex90nvDj\ns7mc3kHaq1mrWiCaCvV6M8kgQxGAFG3o8ZGxd/0iUc+AlSK7JodrYNYGDCajOMJxEx3hy5wwcwdu\nJOAO227Ps4C/mzOgk+Hxioy6nr4uLIY48XYw4HFaCVjKxmbN6vDUpRysXHd65c2Bmdy9IZ9+zbyV\nqDMVqFaS3blGGvBx5yExF7RmMUUYADOUkpyGACAs5ERiIs7uzKLxowN4iiDQxMAaEGUQaMGEZBzj\nCEgjtsxsLMzF5szMpfdt/wDX9rII68G9X3stYzxXIpKoUsnFDToThEM1Wu9OCQRneNt+ufLqkJET\ng8rhvsK3fCZerdjKapPHYiCaWApIi5C00RvHNG7/APSExIXb6nZ1oXg//wAraw/7/i/8uqKQcXja\n1UHjq14K0ZyySlHWhjhApZTeSWVxjZmeQzIiIvN3d3dB6kREBERAREQEREBERAREQEREBERAREQE\nREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEReHN5erSjGa3PHXiOaKAZJS4i8\n0ptHDGz/APSIyEWb63dkHuRF037cVeKSeeWOGGACkmmmMY444xZyM5DN2EAZmd3d32bZBg8/o+vd\nymHy0ktgLGnfpD0WKIomgl9OrjWm9JE4iMuICzjwMNn335N2WxIiDQL3hLjJchPfebJxx3pevfxU\nGSnhxV6fYGea5Tjduu7tGzEDl0zZyYhJnddmj/C3H4qbHzVprj/Q9fJ16sUp13j6WQuvemY2CASf\nhI7iGztsPYub91vax/03U9N+jfSIvT/Q/TvQ+bdb0Tq9D0jh59Lq+xv9qDB6m0HWvZKllms3qV7H\nCMTT4+yMLWqgztYejcCSMxnqEbO7hsz+077s7C7bYiIMFprS0FCzk7UMk5nnLg3LAzFG4RyDBHXY\nYGCMXGPjEz7E5Pu799tmbs0hp8MbAdcLFqy0lqzZeS9MM0olYmKZ4hIQFmhBy4gO3YRZt381k7lu\nKERKaWOITlihEpTEGKWWQYoYhcnZnkOQwAR83cmZu7ry0c3Uns26UNiKS1i+h6dWA2eSt6RG8tfq\nj/R5gzk32syDIIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiI\nCIiAiIgIiICIiAiIgKNPyj/+Sav/AH/hP/MYFJajn8oavJLiqwxRySk2dwpOMUZSEwDkICM3EGd2\nEWZ3d/JmZ3dB5dRawzx5zIYXEVcZIdfF07dezkisxQQHLJOExWiruRzs/GJgijAH/nHeRmZmeOPG\n7Vt/M6JyZywUYHxWZfF6kgaSaUxlqZCq0T409mFyewVQn6vZgeTZ3dhd5RwFeRtZ5qV45GiPBYkQ\nleMmjIxnuOQDI7cSJmdt2Z925N9qivK6cvWtKeINWCpYKxZ1vk7NeH0eVpLEEd7F2Ckrhw5Ti8cM\nji4M/Jx2bdBuNXWOsyyE2FarpyXJ1sdXuzWGPJQYyJpLdoePUdzsSmcA1BEGjHY2sk5OwgBZjJax\n1DdvX6+nqOLlrYCwFW/Zy1ixE9y30BnnqY8a4u0RA0sLPJM/Hc3bbb2l5vDfKDkdUZXIwQ2gq2sF\njAjK3TsVC6kVq8EsRBYAXYxdvL62ISbdiZ3x+GzBaTu6grWcXlbUGazVnN4yzi6E18LJ3IYevTfo\nM/RsjNCewntuxb+WzuGYteKktnGYebFUhPKaltTUqtG9IUcVSxWaf6RO4YDzKCu9aXfgzEbcXZm3\n7apgMjlG1rYlylepBkaOgJxP0KSWWjY4ZYJ47Fd5WGZoS6nBwNmJihkbd24k/RT0lkMTU0xnbFCe\ne3hcnmruWx1H/aLEFfNtO0vSjHvYkrsVbkIb78Ddu27rI4O1NmNZz3ApXqlS1oWenVnvU5oOoz5Y\nNpCCQWeEnI5uMZ7EQRMe3E2dBlNSeJmRr6EDU0cdJ8idLHzvFJBO9LnYuVq8rNCNhpeDBMbs3V3Z\n2bd38nzOr9W5mXKS4bTtXHyWcfUr28lezE1gKcDWJDavVCOozzSWDjhmPl2EWYf0n3Zob1Xlbcmg\nJ9PfQ+XDI4qtSrXmPHTtDH6Pkq3E4JmHjb6jgHEYuT8SI/0B5PJmRtFpzUmYyc9HIWaGqqmMJrWN\npy3Xgu0Ipqr15oYGc4+pFJC4k7cXdnbz32DSPGLPZzJVtLTFQq4+avrKpVerkAt8mz8E1iGlLHIw\nsM+DMOuRSR7m/IGByZid9q05rmTH3NZWMvSxMdjTuN07PkbOFqHFPkbcuJKSaOWzYk52ohnZ4q/U\nYHCM2Ynfu68PiHNmr2HweVyGPljLHa2o5g8fVqTHeq4SGay0HpNeMpCkujHJG5sPFm5dxBxJmxGT\n01NqKfxFjrQzsGdoaUs4s7ME1QbEkOP9LijZ7AjwNpYowMS2cCLY+Ls6DeKWuNRUZqE2oMbjYMdn\nLlalB9G2bMtzGWLLMNaPItNG0cwlK7RucXFhcu++zcvh9Z6muXc5SxFLESPgslFENjISWoYyrPWh\nm9HYYSN57xk8vtu8MYM8W7G7usZqPUs+qPorFV8Rlako5XH3c3JkaE9Wtj4KcwWZomsSCwzzyGIB\nGIb8mJy7MzraPCWvIGV1aRxyAM2diKIjjIRkD6Pqi5xkTbSDuztu27bs7fUg2Dwt1YGcw9DKhG0X\np8LkcTHzaOUJDhnjE3ZnIRljNmd2bszdmWzKLfyUKssOkMRFPFLBKHp/OKeI4pB3yl0h5RyMxDuL\ns7bt3Z2f61KSAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIg\nIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIi\nAiIgIiICIiAiIgIiICIiAiIgIiICIuqxYjj4vIYRsZsAOZiPIy/RAeT9yf6mZB2ouobUbn02kBz2\nJ+mxi57C7Mb8d9+zkLP9nJvtXagLE6h1Lj8f0/TblesVguMEcsotLMX/AEYYv05S7t2Fn8136hyL\nVK0s77bxjsDP5OZOwgztu27cnbf9W6rvq3D4yvWPJZu9duEVo5bMssoxvY6jg0VbaAWOCs3ShjYI\ne2xOziTO3HZwvC9bzn0YuK4uMM8seqeKOs8XNOFYbsQWZt+jVsc6tiV2fZ2ihsiByO3/AOLP9qz6\n/PbxE1q+YsdWZ4Iwr9V42jhlljldiEgcXm2I32bhu/HuxF7LOzNZv8ljxDLJV5cXauw3LWNhjlrz\nDzaaWnv0SadpCIjkjkYWc3fdxni3b6y0cX4bOHHz1n5vLheP61uWybkWs601xSxM1OvZC3LPlGsP\nTr0KU92aX0cYjn4xVxI34jKJeXkxP5C6ymnMwF6BrAQ24BIiFor9OelO3F9nd4LAjIIv9Tu3du7L\nnTjvFYvMeUuhFomdMkiLV9R6+xdCtkbc87vFg5YIcj0YZJDglmKEYo+LD7ZP14nfjvsxd1FKWvOq\nxtM2isbltCLquWQhjOWQuMcQuRls77C3n2Ft3/Yy6MVkQsi5gMos223WhOJyZ23Eh5t7QP8Aa32K\nqXsREQEREBFq38uanUxMLx2Gkz72OiLjFvV9HYWm9M/O/m2aeSCvuHP87Yibyfksvh83XtA8kZcR\n9NuUh6jgLyT1LM9WwMbMT8tpK0369gd9mVprMI3DJIsDjtVQWKlq3DHKfoMuTiOq51o7Eh4+7cx8\nvDrTjCASTUZuBymA7OLk4e1xzI2o3PptIDnsT9NjFz2F2Y34779nIWf7OTfao1MJdqLz+nwbG/Wi\n2iNwkfqhsBs27gb7+yTN32fuvQoBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQER\nEBERAREQEREBRz4xYizcdoY6DWo5sXkYOvHTo2rAyytDxqsWQPo1IJBByI3il5vCLO4OItJIyK1L\ncs7hFo3GkYYLS0gDWsHSZro6gisyWSAHsNC9Vq8srzP7fFwKQH792MvtdSbLIICRETCIM5EROzCI\ns27uTv2Zmb619LgwYmcSZnEmdiF23Z2fs7Oz+bJa3N6kRr0aL4uXRkw809aQJxpzwST+jmMm0fUY\nDIuDvsINJ1H+wY3dVG8TNVZYciMgP06tEP8AZD6AnEYyxx9dpJXZ+TyELg4bt7I7bd93sE+hLGlM\nlZyOHrPbwWVbjlsLELFJXH2n6laN+0sI8z2D+iMhC7cXYg1jU3hZ9JQzSact1bVO1Ge+Ptk4WKzk\nL8Y2c/0hEnZ26nEm4D3Pbd/ouAvhxU1M7ifSfrE9p7ftL5rxCmbJm5q1848pj6TXv8f3hAoZma1O\nBekUajscw8AEwqszEEvPpnKQiJHMZbMzM7x7vu49rA/kmwTSZG9M16pagq0AjMa0W5BJYtSyxAVg\nZzF3EYpuQluW7s7O4v1JtW0v+T9nWkaS2+Kx8Qm7ySzxx2eIDLWLsMUrc3Ia7l3MO0kg7tz3Gy/h\n/pyvjagV6ov0+RSS2CAQktTGTkcpCLNxDd9hb7GFm7Mzl5+IcVjnHy0ne2vgeHyRfdo00nxpw89n\nLYWx9H5e1WxtTJ8psDaCtdisWXpxxsMpWYCjHpwzM/E+7Su222+/XlcRauVcPQ+i7RY6YrhXJdSQ\nwZ3JVZms7Vxmis25Itj5nIM5lO0cYAPEOzNLyLlV4u0VrXX6fn8f5dbpRuZ7q/8AhPoKyVzDtlMX\nMNbGYzLziN8KkkUNmbOCVCp0YvzYSR14BsM0YDGJEDgw8WZvNd8P7kmO1RSfEnHPm9XAcdmvFUjk\nLBvkqk8RQShuPEBhnk4E3Z5HImcidnsSi9J8Ryc3NqP+99/VT+nrrX56aRZS0c1TUuSuVsc8OMsa\nfgr2K1OKGGLJX5L8s81o2GQBktRRPtzNnN+uexb7sWdKjfOi8XGfpRXQeOGR4StSUR/SjLk7xm7O\n7bCXnwdnZ2dmfdkWTLltkmJt217PatYr6NNPAs4Voo47JQ/SgWJRsBAIgIRFy4xRCIhEREzbbN/S\n7d93+87g5escNaPhWykdeOwUfERrtATCTiG2wsUGwszdtxft3Z229F5rNMuYIjtTdUbTByh9CkqD\nVdowERZowKQHOtxJv6LiJNvvv3WV07juNi7Ykh4SSXJGhIxFy6LBGPMC8xE3Z/28W38lnkQRJY05\neEs4/os8gYmeK7hBGT87enfLHqa3HATy7dM5/Q6rCXFmeuY/oOyx2X0XkjCq0scsh2cBUrc6lOhO\ndDMPZsW8jajnvS86JyTz1ZhlBjZixYuXtDEJzYi9ozWhTkhX/Ueir82KytatipILx2Nezy2unUZ8\nlUyT5+OhSCYZOqTzS2sXMwkzMzY4OXH83y3mTTBjIdqOiD2/5TSWgncQCV65QHX59dmc44ijMo3d\nv6Jl2fd95HRROaZTywhDA6RtSWahz4t4YSixAWq/0fj6lQCq2bBlDHXryyyHDD1SZnmlkbYncH2J\n3KbhZmZmZtmZtmZvJm+xlyire83nzK1ioiIqLCIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIi\nICIiAiIgIiICIiAiIgIiICIiAsNkdLY+xJ1ZKsfVd93mj5QyO/2kcTiRP+3dZlFat7UndZ0ralbx\nq0bYungKkTsQxciH9EppJJ3H7HHrEXF/1ssoiJa9rfqnZWla+VY0IiKqwiIgIiICIiAiIgIiICIi\nAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIi\nICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAi\nIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIsNqPU9LHlCFqSQCtDKUIxVbVlyGJ4mkd2rRHwZnmib2\ntt+fbfusDW8WMBLIcUV4pZYf52GKjfklj+r85GFdyD+1mQ23dFqX/wBR8R99a/dOU/8AaryTeLWn\nwmGud9wsS8enWOjfGc+XLjwhKvzLfiW2zd+L/Y6I3DeEWpf/AFHxH31r905T/wBqtZzXj7pupPJX\nlnt9SHjzZsdaH9IBkbtIAk3sm3mzf3bOidpTRRB6x2l/v7f7vn/gnrHaX+/t/u+f+CCX0UQesdpf\n7+3+75/4J6x2l/v7f7vn/ggl9FEHrHaX+/t/u+f+Cesdpf7+3+75/wCCCX0UQesdpf7+3+75/wCC\nesdpf7+3+75/4IJfRRB6x2l/v7f7vn/gnrHaX+/t/u+f+CCX0UQesdpf7+3+75/4J6x2l/v7f7vn\n/ggl9FEHrHaX+/t/u+f+Cesdpf7+3+75/wCCCX0UQesdpf7+3+75/wCCesdpf7+3+75/4IJfRRB6\nx2l/v7f7vn/gnrHaX+/t/u+f+CCX0UQesdpf7+3+75/4J6x2l/v7f7vn/ggl9FEHrHaX+/t/u+f+\nCesdpf7+3+75/wCCCX0UQesdpf7+3+75/wCCesdpf7+3+75/4IJfRUA9dXVX4fp/4PI/MU9dXVX4\nfp/4PI/MUF/0VAPXV1V+H6f+DyPzFPXV1V+H6f8Ag8j8xQX/AEVAPXV1V+H6f+DyPzFPXV1V+H6f\n+DyPzFBf9FQD11dVfh+n/g8j8xT11dVfh+n/AIPI/MUF/wBFQD11dVfh+n/g8j8xT11dVfh+n/g8\nj8xQX/RUA9dXVX4fp/4PI/MU9dXVX4fp/wCDyPzFBf8ARUA9dXVX4fp/4PI/MU9dXVX4fp/4PI/M\nUF/0VAPXV1V+H6f+DyPzFPXV1V+H6f8Ag8j8xQX/AEVAPXV1V+H6f+DyPzFPXV1V+H6f+DyPzFBf\n9FQD11dVfh+n/g8j8xT11dVfh+n/AIPI/MUF/wBFQD11dVfh+n/g8j8xT11dVfh+n/g8j8xQX/RU\nA9dXVX4fp/4PI/MU9dXVX4fp/wCDyPzFBf8ARUA9dXVX4fp/4PI/MU9dXVX4fp/4PI/MUF/0VAPX\nV1V+H6f+DyPzFPXV1V+H6f8Ag8j8xQX/AEVAPXV1V+H6f+DyPzFPXV1V+H6f+DyPzFBZL8p+5DHL\nh4pp2hjtPMFkRn6ViTGtksEWY9GYHaUjHHtbJ3i9sREiHZ2Z1q2swtxYSM71fCVwDGRPpqTDxVI2\n+n2Gd4QxpQ2HP6OImq9Pgz71Hvek8BZVP8afHPMarkxst6KjUkwZ2JKcmLC1XNpJ3rERkctmQmMS\nqROLg4uz8vPttp9zWuRlCUHkgj9KEgsyVqFGpPYAn3kCxZqwBNOBv3JjJ2L690Quy+tMM/dsvinZ\n/LbJ09tvq/rVrcmtMZE4Qen46eH0fUcdun9MYsauWnyliueNntSyzPJHHWiDhI7/AJwGpQNAEg7A\nqV8/1D/cthoa3yMGNlxMMox0rJynKAwx9Quq0LSD1tuoIO0As7M7bsRs+7Pso0iI0unX1hiBABLN\nYsyEBEpHydNnMmZmc3bq9nd+/wDaoQ8T70NjK3Z60sViI3g6c1eUJojcadcSYZI3cS2ISF9n7Ozt\n9Sr1z/UP9yzWK1PYrRDDGELiDk7OYm7+0Tk/kbN9f2JoiNLBPp2D/aeNU3KrKMeJAntk2fhKV4zu\nRjE/UsxRRsNknpbN05Q8v5wvHqfAR1ojeOvKTDDHPNaaSThUsyWI42x0nWEHZ2gkGV4TBrTdeAiY\nACVihH+Wdnv+ard/P83J3/b+c7o2s7P3Vbt2/m5PL/EUrJ2o6YhlrOTxTxSFRrTV5hE5CkOWzjY5\n5JYzIYp4Q9KmdjqEbRRvGNgY5tmP5t6EKIpo5LEgT13cigOmHMo2r4WfpA4WiaTJGWahhjrhyCSS\nPYZtpAJ4L/lnZ+6rf4cn+Yjays/dVv8ADk/X/wBZ+t/73QWBzfh10q422sHWgerRcmv1fR3js2Hq\nVw65POTV4etJMczv7VcRZhCfs6x2p9H+i42C8HpwEQkE0F3GvUkE/SZwGSZisk0LuDRC0Q83dhY3\n4sfaD21lZbyirN+yOT/MQdZ2W8oqzfsjk/zEFgcppPHs5BDZYT9OepzGR7UbRWGbK4y3XjhZ5rLy\nY0Zawxs356zxEPJzLl/D+N44HCye0j1YTnKhOISFY1BNi+pAM00ZNagiGNp6xi3Tkdg5O5tIFfv5\naWvuq31/1cn1+f8AWLj+WVn7qt/hyfZt959nZBOMmiNq0Vx7ErVp4+oMj0g5OxYe7mYwYGtu3W6d\nNgMSIRF7AOBTcZGHVa8PISfaR+P3cXMW/NTSe0XJuP8ANb+T+y0pf0Nijj+Wdnz6Vbf7enJ9u/3n\n2rn+Wtr7uv8A+CT/ADEE06R0/Haq2p2ErM8DyC9OP2SjhGHqNZ6vXBoikldohlJphicSIoJxcnhy\neY0jFCNxir2YjqNY6ZmMzRFUCrkrNO6RyNt6VLJSpxSiO8I/SDAzRzG0daA31pafzirPt9oSf5i4\nbWdlvKKt/hyf5iCcM3hK8dKWetCbsL1T69kZ3eOKWhhJxjjlilaD0w5chMTxHGe8QTELx9IepxS0\nZLJjzskMsUteC7eJnjYwloR46rerzM/UZxjLeeNjAZH6hsJMDRylHCH8tLXn0q2/29OT/MT+Wdn7\nqt57/wA3J5/b/OeaCesdoUXvVq8k0hxyXoYJX9FeMWj9KaEimkaZ/R47Q7jTPYvSjJh2i81htS6Z\nKnVjsc7BvIUQkEtB6zQvLXKcIpjKctrrcJBOAGPpvEXMh5RdWHn1nZ+6rfX/AFcn1+f9Yjazst5R\nVm+rtHJ5f4iCYdZYgKhSsEcgOGQyXTGTnuWNjKoFCyzE3tV5CknEZv0TcXZnfisy3h65GzDZneJ8\nnYxrzfR7ORTww5GU+lDHaIjfljTbY3D2bERM5bSjFBMeurgiYCMIhLx6gC0rCfF9x5i0uxbP3bfy\nXX/LKz91W8tv5uTy+z+c8uzIJ+wfh6duvY4FO1iI6/SlKoz0i5jZYoGsx2CaQgk6ATkAF6OcUos0\nzbkul9CxvMEEdyxIb2pYJG+jI22aN81HzhZr7vMTyYOz7D8PZlidnIuUYwO+srL/ANVW/wAOT/MT\n+Wdn7qt/hyfq/wCs/U39yCas1pQKVmpXs2mbrXRqXpemMUNV+jjLE0kU0kj9aAY8k35yQIv5kncW\nZ17Z9KRvXa3YCxi+PWY6kdKzPKwxUbmRcijyc1c2l6NKUOQEURyOLC8e0jRwSOtLTeUVZv2Ryf5n\n61w2s7LeUVb/AA5P8z9Tf3IJ8yWiYYadqY5yienBNG9izWsBVknj1DbxvWhOv1ZJGaCqwHEMb9Mr\nMDuT8ndtV1PUCvevV4/5urftwRbu7/m4rEkYdyZnf2Rbu7M6i1tZ2W79Kt+3pyfZt959i5bWtr7u\nv/4JP8xBrKIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIg\nIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIi\nAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIi\nICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAi\nIgIiIP/Z\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"700\"\n",
       "            height=\"400\"\n",
       "            src=\"https://www.youtube.com/embed/MnPABAAQnXM\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x108a5dc18>"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from IPython.display import YouTubeVideo\n",
    "YouTubeVideo('MnPABAAQnXM', width=700, height=400)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "  "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### WCPS Webinar Part 2: Advanced Queries"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/jpeg": "/9j/4AAQSkZJRgABAQAAAQABAAD/2wCEAAUDBAgICAgICAgICAgGBwgIBwcHBwgICAgICggICAgI\nCAgIChELCAgOCggIDRUNDhERExMTCAsWGBYSGBASExIBBQUFBwYHDwgIDxcUEhQXHR8dHB8cHBce\nHB8bHBweGB4XHBgeGBwcHh0dHRUXFRwfFh4UFxcfHBweHBQcHRwVHP/AABEIAWgB4AMBIgACEQED\nEQH/xAAdAAEAAgIDAQEAAAAAAAAAAAAABwgFBgEDBAIJ/8QATBAAAQQBAwICBQgGBwQJBQAAAgAB\nAwQFBhESEyEHMRQYIpTVFSMyQVFTVFUIFmFxldEzQoGRk9PwJEOhsURSYmNyksHU4Rc2ZHXx/8QA\nGgEBAAMBAQEAAAAAAAAAAAAAAAECBAUDBv/EAC0RAQACAgAEBAUDBQAAAAAAAAABAgMRBBITUQUh\nMZFhcbHB8DJB0RRSYpKh/9oADAMBAAIRAxEAPwCmSIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiI\ngIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICI\niAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgI\niICIiAiIgIiICKyj/oY6o/H4D3zI/D1yP6GGqH/6fgPfMj8PQVqRWV9S/VP4/Ae+ZH4eufUu1T+P\nwHvmR+HIK0orKF+hhqhv+n4D3vIfD09THVH4/Ae95H4egrWisn6mWqPx+A97yHw9c+plqj8fgPe8\nj8PQVrRWXD9C3VL92yGn/fMj8OX16lWqfzDT/vmR+HIKzIrM+pVqr8fp/wB8yPw5PUq1V+P0/wC+\nZH4cgrMisz6lWqvx+n/fMj8OT1KtVfj9P++ZH4cgrMisz6lWqvx+n/fMj8OT1KtVfj9P++ZH4cgr\nMisz6lWqvx+n/fMj8OT1KtVfj9P++ZH4cgrMisz6lWqvx+n/AHzI/Dk9SrVX4/T/AL5kfhyCsyKz\nPqVaq/H6f98yPw5PUq1V+P0/75kfhyCsyKzPqVaq/H6f98yPw5PUq1V+P0/75kfhyCsyKzPqVaq/\nH6f98yPw5PUq1V+P0/75kfhyCsyKzPqVaq/H6f8AfMj8OT1KtVfj9P8AvmR+HIKzIrM+pVqr8fp/\n3zI/Dk9SrVX4/T/vmR+HIKzIrM+pVqr8fp/3zI/Dk9SrVX4/T/vmR+HIKzIrM+pVqr8fp/3zI/Dk\n9SrVX4/T/vmR+HIKzIrM+pVqr8fp/wB8yPw5PUq1V+P0/wC+ZH4cgrMisz6lWqvx+n/fMj8OT1Kt\nVfj9P++ZH4cgrMisz6lWqvx+n/fMj8OT1KtVfj9P++ZH4cgrMisz6lWqvx+n/fMj8OT1KtVfj9P+\n+ZH4cgrMisz6lWqvx+n/AHzI/Dk9SrVX4/T/AL5kfhyCsyKzPqVaq/H6f98yPw5PUq1V+P0/75kf\nhyCsyKzPqVaq/H6f98yPw5PUq1V+P0/75kfhyCsyKzPqVaq/H6f98yPw5PUq1V+P0/75kfhyCsyK\nzPqVaq/H6f8AfMj8OT1KtVfj9P8AvmR+HIKzIrM+pVqr8fp/3zI/Dk9SrVX4/T/vmR+HIKzIrM+p\nVqr8fp/3zI/Dk9SrVX4/T/vmR+HIKzIrM+pVqr8fp/3zI/Dk9SrVX4/T/vmR+HIKzIrM+pVqr8fp\n/wB8yPw5PUq1V+P0/wC+ZH4cgve3fdlyD7bMuCLfyXIt/eg7vL+3/muEFv72TffyQcG27ftZdLs+\n37V6Nv8Aiut223/4IPn+Xf8AeuGHzb7V9f67L7ECfvt5IPqp27f8/wDivQukYib6289137ICIiAi\nIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgwWq9W0cW9YbkkzSZCUoakFWj\ndv2JjAHkk4V6EMkvAQbcj48R3bd23bfE2vFDCwyTx2J7NR6sduU5b2IytOuYVIzlsvXtWqoQW3GO\nMzYYTNyECIWJm3Xs1NgWnymGyBWQibClkCeAg5FYazU9G9k+bdPg7sT+yW+7N233UWUPCGYThezl\n8ZkGq17sUZ3aGSsyXitU5abnlymzBxWRKOyRkMIQOZ8eJA24oJEPxTw4g8kj5SJurHCAWNN5+CWe\naRpDCGrBNQaW3LwhmNwhE3EIjJ2YRd1jNOeK1WWG3ZtRXxhbLXKlIaun85YnaCtHW3K5XhpnNVkc\npSfaYI927Mz8XdazV8ML8UgXxzmOCelepWaEEWNyMmLGeCtl6sxzVrGZlmmnlhyzs7xTR7ehw7sX\ndLvhjdlMpXzeNlsWJrc92KbH5Mcf1rcoEElSvUzURxScIGB3lkmd2Fnbi3Zw3G/4xaehZie3Ymjd\nqn+0UsPl71ZiuRwS04ns0qckIzyhaquMTlzf0iLtuY7+2LxMxJwHYAshIMVr0SaCLAZuS5BP0Ass\nNjHhSe1ADwyxGxnGIu0gbO/JlqenfCyGji5sQOUA98ppy6MpwMxiOGgwEXSKNpu5zNgifdnbj6R5\nFw9rr1Z4bS28jcyFXM1Qa/e9JOjbhuy0354vEY+IpBoZOsctkPkoyCRy4sN2YeL/AEkGy2/F7ARA\n0hW7Jx+iDclkgw+WsBVrOU4dS+cFQmxziVWyxDY6ZD0JOTNxdZDUuoJ4MrgKsLxlWzEuQGy/Hk7h\nDQOzC8Zs+w+2Ld++7OtIwPhdHWxebxz5WuX6w6fbE9UK7iNaUmzBSTiB2SKWLllfZjI+XGv3M3Jy\naRLWIH0nEmTs5Y1pxF9tuTlUeEnZvq7LJxufoY4tP91I/wBr1r92rhJxxknnjcctvflnXtOpefW2\nbmqvFHDsLyC5FI4sXk+zCLP2/f8AvZR/4h+L+SxI0uhhiyfpgWOqcR2A6ZQnEzbjDXkbY2lbbu3e\nM1tvjoWQiw1i3jI4ZbWOdrLwT1+u01cWdrICLOxsTA/Ubg7O7wM3ffZ/J+jpq98zga9mU4itQzWK\n9wYQCMAkGYjjZowbYWeCSAm/8fd3fd36+OtYx9SaxMROnNvaZvyROp0ii7+k5koGZ59NdBjd2F57\nlmFidvNheSk3J23by+1T14a6kfMYmjkihaAr8HVKBjeRo35kDixuLcm3Hz2bzUIfpp2/STweHrA8\n9+zZknjrx7c9jZqlcXZ/otLIcjM79v8AZz38lNmlaMWDxeMx5E5tTqx1XlFmYSOGrJNNMXJ24iXQ\nlL95Mr8TXF0a3rXUz8ZnyUwzk6lq2ncQ2VFh6+oYZWjKIJZBnIAicWjZjkOuFpo25m2xNCTm7vs3\nzZNvvsz9FLVlWZw4NLxkOsIGQiAO1jn0TciJuLPw22LYncwFmdyZlga2fREQeAA/Z5ftXLu232Oy\nb7fazr4d90H1GffZd266wBff/NkDb/4Tt5rhy/1v/LzXy5N/pv8A1dB97t/6rshd9/L966RkH693\n/wBfsXYE7v2EfLz/AGftQehFwzrlAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREB\nERAREQEREHVPXCT6Y7+y4+bs+zuJOzO3l3AX/sXn+Sq+23TbbYG+kXZgbYGbv2Zuzt+1mfz7r2og\n870ouIhwbiDu4szu22+7l5P3Z933byddTYmu3+7bzZ/pH3dm2bfv32bt3XtRB55aURb8gZ933d93\nZ9/a2fdn7P7Z/wDmXW+Mg4iPTHiG3Ft32Zm3Zm8+7d3XsRB4oMXCDs/FyIXZ+Zk7l7P0N9u2w+Td\nu3f7XXXkLAjZpg7+1MU7A328YXJ/+CyK1TVJ7ZfANvtymyO7fbtQN1g8Rw9XDFf86T7XrP2Vvl6U\nc3yj3nX3bWoY1PpLK6csvktJ0qs9O0cp5bBMPTlsTSFuFiGwZPxij7MMAcBjZ5eIl1C4zOi6ePLO\nOe8du6L0i6KPCXwzswXZdQ6glC3n7zO4iHetjIiHh0KzeTyNG/Bzbsw8hFy5Gcm/6qKQQhOOAZ3j\nllJ2KE5um7U7TxGIA7PyeXpB9faV/LzbMImTLbJbdvz5FKRSNQ04o7PMgChUYjlmYJipH0i9iy8R\nybGzxl1atcSLuxC8RMXtMIfFiScX5w46FnMKcsIvjpGkZ5LhDwlcSYRkir8Tfd24lu/kt0Rea7Vj\nzGQcZXGuzPHWszR707RDIYcHhhZiITaQmc2cePmzcXLus1jZrBSzhMwcYXBo5AikjaXkAm5NzJ22\nZ3cdmd/LzbyXvRBjtt13CH/8XpYGZ99u64KNndB0OX2f6/8AlfBP/r+a7JjYW9lt9vN/qb+xeeRv\nr+p9kH0Tt9vf7FzxZtn27fWvnZmdvsf61x9vfuyDvCLfdtu32rsiiYe7d3+t1xWdtv8Amu5AXK+W\n/wBbrkUHKIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgI\niICIiAiIgIiICIiAjsiIPNMPf9hN/wAV1CPbi/177L1Tju37u7LqcHdvLv57IOoR7bP5t5Lt6T9n\n/sdc9J+z/Wu5v9fag4Ftv9f8mXP+u65ZlyzIOP8AXdc7IiAiIgIiICIiAiIgIiICIiAiIgIiICIi\nAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiIOHR/wDX2rlGZBxt/YuW\nZEQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQERE\nBERAREQEREBERAREQEREBERAREQcE7Mzu77Mzbu79mZvtdRLc/SDwASm0Q5K1UgMwnytPHnLj4XH\nZncpnJjMH33YgAmJtnZ3Zxd8x+klkZaulczLCxOZVQrvw35NHZsQ1Zi7NuwjFMbu/wBTM/l5rY/D\n/D1qmIx9OCMGrx0IB4cWcT5RC8hm39cjIiInf6TmTvvugyOMzVSzUC/BYikpyw9cbTGzRdJmdyMi\nLbgw7Fvy248XZ9tnUZTfpDafGR9gyZ0hk6ZZiPGyFjhPmwcepv1Sfd+3GN99u3m28SZW/NR0drej\nV5R1sfq23j6vH6MFI70AHXHt/RuwuLt/+U/2qdfEjK43B4mChLib9/H24ixo0sRUjsuEHQJn6onK\nHEHFtuTO78ib96DdZMvVGo98rELUmrelvceUfR/Run1ev1d+PS6ftct9tu6jGP8ASF0+8osQ5KKn\nKfTiy82OkDHGfIh2GR36rfRd3d42Zm3322fbR9cZ7HXtF4mviK96ricjqCliDhtyOU8VYLEssjyT\nDYldo3kgBhfqduYN27Mp+zenaVyhLjLFcCoz1/RirBvEIwsLCAxPE7PC4cRcSB2cXAXZ2dmQeu5d\nGOvJZFinCOA5xGvxM5hEHkZodyYSImbtu7M+7d1qtbxKx8mnv1lELL0fRjsdHjA9v2JSgKDi03R6\n/VFw26nHf61tOHxsFOtBUrRtFWpwR168IuTjHDGDRxxs5O78REWbu/1KrFreLBXNI7O3LxAr4asA\nvu4UJ7Q5CtK+7u4jygN92Z9uQ7/W6CfsT4k0LOn5NRhHaGlBXuTlXkCAbf8Assk0RxMDTPF1TOF2\nFnkZn5hu7b9tj0vlwyFGnfjjlhjyNSC1FFYEBmCOaMZQGUYzIRkYTbdmJ2Z/rVXs3IVfCZ7Ssb9K\nS3ryDFUAjfbo0708dyoPmzuJDWkZ9+z9Tju7LfPGbWk1fLVNPVM3W0zUrYtrlrKzwRzkbubwVsfC\nEjswlxDqPsQu7d+TMDiYSr4j6ur4LGWcraisTQUngaSKoMRTl1rMNUOAzSgD7HMLvuTdmfbd9mfN\n0bAzRRyizsM8YSCxbcmExYmYtndt9n+p1V7VuupcrovVFK1erZOxgr2KibJ1BCMbtOfKUyrWSij9\nmM3eKcfZZm+a+t2cn2LNZTVGnquNz93L17OOlkoRZHBhQjhgpVJgYfmLXIpZpI929t+LuTs78hbi\ngsKow1T45YSjanqiF/IHRkeO/Li6fpFekQ8mk9ImIxFuDgTFw5bOzs/cSZt115kJKmKyduESKalj\nLtiEQbcykirSSAIts+5OQszMtU/RrxsNbS+J6Isz2q/pU5+ZSzzGRySSE/cy8h3fyEBbyZmQbfpH\nUdLLU4r1CcbFawz8DZiF2IX4mEgGzFHILs7OJMzsssoi8RDi0nj46mna0NW9qzPx1q7yEckEV26T\nCdso5XJhABCOMQZuAfNNxcQ4vj61/UGmstiIMvmQzeO1Jaei8klKGlNSvH3r9JoyfqBIZCGzlt57\nCLs3IJtRQngdf2sXW1nXydkrNrSlmWxSkscWKSnbh54mEiZmYyc2EXdvLrCz/VvjNTa0yuI07pur\nYyo18vqR2O5mcmMRPQrvxs2ZHiYOmUkI2K8DM47fa+7sSCf0UAeF3iDa9PyOFkz9fUUfyJNkqGYr\nwRV5IZo3eKWlIEZExk24ys7uTs317Ewh7f0eT1LksfUzdrNtZhOrkIIMVNVjAbE4WrEcVm5eBnMd\npQIOMcfYI4n3d+TEE5IquWtc5zF3cVJY1djsvau5ivVymn8fXqzV61WSUwkcLUexNxbgI7jGbFIH\nJzYS5bxmsjqDP6hy2MxGZHB1NLx0xlmChDdlt2rUZSi5BO7M0YdOUNt9t493E+TcAmtFXDAeKGYo\nYTV16/MFvIYjOSY+pFsz1obBTDX6cA7CR1Y5CM2En5EMTM7s77r50X4i26OWxMFjV1DUtfO2WpXI\nIq9eA6FqURGtJBJDt1YyncYu7ALsXaNnJuIT3rHPxYuhbyM4SyQ4+A55Y67AUxiLdxjaQxBy/wDE\nTN+1RfX/AEjMRxCWxjc/SqyMBNetYyN63E+PAmKtYkI2JiZ24CW/1bu7M+3ePv8A9sZv/wDWzf8A\nJlDtTxKK3pSHT9PT2oLVybT1bExztjg9BeQqUdT0sLAyk7QM79UTIBbsLu4NuTBZWlZjmjjmhMZI\np4wliljJiCSMxYgMCbsQuLs7O31Ou1V4zUWex1vRmmqOU9BmtYWzWtydILMATQ0ic5RjNm6rwiMj\nxtuLOUcfLtusnas6jv5NtL0c81c9P4yCbN558fBJZt252Y4YQq8mGGPpHGXsmxee5Ft7YToihjSX\niJkMfDqijmpIr+Q0XANkLcQDXbIV5qxWaonGDcY5v6MSdm2brg3tOJEXg0lW1xNXxWeDMVr4Zb0O\n1awJUa9WtFj7AtM7Q2nLm0wxGG3ZnYi9opWF2MJ2RV68SfEC1Yz2QxsWp6mlqeBjrCMs9aGzLkrk\nsbTSbdZ24wxNsDsz7M79xPm3Dw5fxjyculad6G1DBYHPNhs3mqlULcdSMWcvTY6xNw5SxHWNmJhZ\n3l4jweSPYLJooT8CdRZGxYtiOo6OpsfHR6sfOA8bl4LfPYY3pHHyaobc2eSQ39po+LMzFvGuk/En\nOXzCw2rqVTLvaeM9M5jFFRxwbWnj9FPI9MuJcWZuPFpmd3Dlybk4WNsa6qBn4tOvHY9NnxXyoMzN\nB6M1fry1+Du83W6vKEn7RuOzt7W/ZtpUW3dQXh13UxnpBjRk0t6XJSHi8T2vTLcbzcnFjcuIAPmz\nbC3ZRpi8zqu9idSZKDUD1IdM5XMGAPRgsWbIVa8Vj0XqyNwr144xFw9g3Mp5OfYRQWdRQhqTxDyt\nrEaUgx8sNXMazaIZLbxCQVgihjK/Yghk5D9OQNmLlsJO30tibx0LepcXqzTuHyObfJ070WRleYak\nVQrIhQtm0NmIXLvFLFEQmxblz7u+2yCccvkYKkE1qzKMNerEcs80j7DHGDORE/7mZRdU/SEwJyg0\nseTq1JzEIctbx5RY6Ryd2babk5gPbdyMBYWZ3d2Zndty8WtLHmsLkMZHL0ZLsItDI5EItLHLHPEM\nhB7TQkcQiWzP7JF2fyeFszrq7VxjYXWmmbUeMKKGpPlcY/UrsIHEEErjDvHETOMZNwl5u7Pwj3bg\ngskz/wDFFAHilrU4L+O0/Q1BX0/Qr4WK5Nm7IR25LAvtBTqRPK7bm4A0rkxC5M+/LYXE8Pa8VchJ\npTUjR5aGzktNXsbFDm8eEIPapWshWCGw8PF4wMwa1G+w7ex23dnJwsuir7qizqfDY2pqK7nDtwtf\nxFnIYyClFFXq46beO3WA25S2i5S1RE34P7JE7O5Ot0sZ+3d1lWx1S0YY7CYg72UCF43jtWLTtHUr\nTE4uTcYyjnHi7b99927OGy+G+uqmeivS047EY4rK2cXO1loGcrFcITkOPoTG3Sdpw25OJdn3Fu2+\n0qvGD8UshR0xqnK2JnuWsbqa5QxzWGFxiBxoBADjGw84onmlk2d934u2/ftt2l8DrPH36EtrMQ5u\nlcI2zFaSnXp+g7xtwkqSA/KbjI/kwgziD+wzkxAEsoqw6k8S7927mJ4dWUdPx4a7ZqYvDz1a8x5D\n0ZnYrM0svcY55PZF9jZtnZhZwcjnfwp1U2bw2PymwCV2vvMET7xhYjM4LIA7u78RmilbZ3d222fu\nyDZ0REBERAREQEREBERAREQYzVmDhyVG3j7G/RyFaWvI4bcxYxcWMHdnZpBd2JndvMWUP4SLX+Kq\nR4aClicgFWN69HOHb6QxVgZgrvaqmbSnKIbNsIkzcWZ3k25nOaIIy094SQR6bt4K7O9ibMlNayd9\nhZykyMpBK1oGL6TxSRQceX0ugzuzbuywFKbxBpV48XHQxN0oImggz8l3jH0gZgjks1TNpjscG3d2\nF23bvy/rTYiCIMJ4M8NIPpuxbH0hyOxHegA2Gtc63pEBRM5MZgBMIOW4EY89uHJuOKzdfxAyVMsL\nPUxdNrEQ17uoYbvMZIX9mcoKwO0wTGDPvvGDPzdm6e+4ToiDw6fxzU6lao0s0/odaKDr2ZClnm6c\nYh1ZpDdyOUuO7u/1u6hzL+Gd49fV8wFd3xJxw3blj0iFg9OhpWKUMTV3l6ryDxrnyYOO0h993JTi\niCDMv4ZXpdfQZga7viSGC7asekQsHp0FKapBG1fq9Z5BIID5cOO0h+1u5Mst4paPyseag1FhqlPJ\ny/JpY2/ibpRRNNE0rzxTV5ZXaMZuXEXcibZowZt2ItpdRBBWsdH6jv6UzNezUonk8pbpSU8VjGqV\nxgqw5GrO8MluUgjmlEBmLcj22Ftndy2XRlNI6szEWOwOVrY6LE0Zqct/KxTtJJkK9f6Nb0Ryd4rB\nCI89x6fImcSIRcCnxEHxYhCQDjMWIJQIDAu7EBM4kLt9js7soQweI1lpuEsTiqdDM42OaX5KuWrj\nQWKUEpHK0d2OQwecQMib5t3d2ftxZ2COckQQdP4Q5UsDAB5Jp9Q084OoobE8kstIcgxOT1Y2Md46\nru/P2QFnk3LizO4r24vT+pM5lcZc1BTp4ulpyaSzDTrWAtS3r23CGw5ARNDXDixizkxN5Oz8t45k\nRBBXjh4XX8rn8dYqRmWNykVepqTaavHE1WpditxlLFIbSTySM7iLgJ8XrBvxZ35bd436KuZIcZex\njVSyWnrr2K0F4WevZgkERs1TL+ox9OJ/sfp7bi7sYyOiCKdFY3Udg79jJ4/GYqvJjZK1LFU2rTWS\nskxMU8t2P2QB2fjwYnZ9x3YODvJ5/DXQ+Sj0J8hWQ+T8lNSy1dwOaKRoTs3LskPOWqcgOBRzRu7g\n7uzG/wBbbNLyIKtSeG+pSx+KqNp7FVBwOTx080lGxTbIZUoH6J2jkeQYuHTOSQ2kNjIuHEW24rbs\ngGUx+rM1JpwMfkpcvBQLL4y5YenLQmigYa9tiJm61Nxm5O8XUJzmkHYXFlO60bxE8JsFn5gsZKn1\nLEQNGNiKeaA3jYnJgk6RMMjM7ls5M7tyLZ23dBDnhvombPYXWuOlsRPZs6ns7X4x415L9eWOwcgC\n3Io68p7eW7sFjtu7LetF0NWWL9L5RxmGxdTHu5XrNcK1ibJGwC8ZVRDk9Ruo27uXF29rbfYVJWj9\nM0cRUjo46uNatE5EMYkZkRk+5HJJITnIb9vaJ3fZmbyZmbMINU8YMTYvYHK06kTz2bdGWKvCxxxv\nJI7eyLHMQgO/2kTN+1enwxx09PB4apZjeKzSw+Or2YXIDeKeKpDHLG5xEQG4mJNuLuz7dndu62JE\nEZa10zfsav0xk4azyUMXBlBvWmmri0BTUrEULPEcjSycjMW+bEtt++zd14tc6bzmPzkuodP162Qf\nI0Y6mUxNiYKxynC7ej2a88jiDGwCIOxE2zC/YuTPHLSIIk0h4aXLFfUdnOnDHktaRej2ApfORUK0\ndYq1UYiPtJMAkJO/dneKPu+zu+D09gtb9HGYKQamOoYWWrHNnaNxnluUKzPGFaCBneQDKIQF3kAW\nd+L/AFOzzuiCFdW6QzeNzmRy+GxuPzNfUMdX0ujdkhgkqWa4dEZojmdheFxIzJmdyJzNtvZHfKW9\nO6qHA1xr3MdHmo7oW7laGpXix12tu7HjJHeF3Znj6e8g8HJwceYi/NpWRBCfhPoPJNnps9exdDT4\n/JklGPGYw4DeeaSz1ZLs7194eTsLfa5bRO7M4Py1bWHh3qzLQlh8jRw94zmfhrWVqoXI6vpQS8Xr\nwgEoyvGPB4wBgcW4cn2aR7Kogi2TSd8dbUsmMBnjKulhx53jngd/Sht2j6RRPJ1yNwkAubBx9rz3\n3ZYLQ+iMtX01rOhPSKO3mr2oJcbA9ioT2Y7WOhgqk0gTPHFzkEh2lIHbbcmFu6m9EFedV6ctUMHo\nw3kq0c/grMEVGrkDB4bU8oMM+M9KhIogkk6cez82EumTMTPsY9UEmbta40vYzEFWpONLJyti6k42\nCo1vQbcXVtTMTjJNLPKw7R7izRBs7vydpx1ppWhmahUslXGzXMhPgRGBBIO7DJHJGTHGbMRNuLt2\nImfs7s+C8OvCjBYCWSfGU+lYnj6R2JZpZpOlyY3jDqk4xC5CLvwZnLgHLfi2wevxe0XHqDDXcTIQ\nh6YAPFIYcxjnilCeEiH6w5xixM3fiRbKMs7T17laB4K3j8XXjsxR17mee4MoSws49WSKqBdRpiEe\n+4CzuT7cN2cZ5RBCer/DrI4y7jcrgKlXKfJ+ChwdrFZCSMCnrV+JV54Zp34NY3ABcjLs0Y7b8i28\nus9Hakv6TzVaxUolk8rapSUsVjGqVxr1ochVn6MtqUgjmlEBmLcj22Ftndy2U7Ig1vN6ca/g5sVP\n7D3MW9Q39knikKv02NvMXIJNiZ+7bg3mtG/Rk0XksXTuWc1EcWUyc1cJBlngsSjTp1Y61MJJq8hg\nZM3V/rO/Fw377qXUQQHp3wnv29Oanw9+L0GfL6juX8dJJLBMJRs1GSrMT1jk6cUh1zjJn2NhIn47\n7b57TMGtMhkMcWWjr4ijiCMrr0rgylmZOmLB81GRNFX5tu4nxf6bbPuztLyIIDsaM1JhreWhw2Nx\neUpZ3IWL9W1dOCOXFz2R+caWOXbrVwcAYQBi7CLv3ImaZ9IUbFajVhtnDLajhH0mSrCMEBTPuUrw\nxCzMMfIn27M7+btu7rKogIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiINS\n8ZMvYoYDLXakrw2adCaWvMwRm8cgtuJMEwEBfuIXb9i0ccFrSCkWRDV1e6QUntR4+5pqlDBKTRNM\n0UlirMMoM+3HkP2v2+zav0gv/tfO7tu3yZPu32+ytXxXh3qO3Tjhu6ysHj7lMAlqU8DjaU/QOMWe\nILo8zD2HceTDu/8AezhvOB1zSlwlHN3ZoMfWvU608hWpwiiikmAXeHqybMb83cW+stm2buuyHxAw\nh0ZsmGVonQqGEdm2FmM4oJDMI4wmcX3jMiljZmLZ35j9q1PX1TG1pdPYmtiPlS9jop5sJjjuPWpV\noatUapWL0kjkBALSxALlHKfIuQtuxOoW8YQnhi1u09ehTnkxWl5bNXFSFLXjlbMR9PnMcMZTTtGb\nO5dMN2Me31uE0eKXjBjqFZ2oZLHTXRycNKaLrjMUAjZhjyJEEZdpIYpHd+XYS82fZxUkYfJV7leK\n1VlCevZjaSCeIuUckb/RMCbzF/tUX+PmEp1sLSir1YIghzeEijGOEB2D0+vG49m8nEWZ/tZu6liG\nIQFgARARbYQAWEWb7GZuzMghTKZvV0mPvaljtQY6pjWyNiDTN/CyBLYo05JG6l21OQ2a9iSKvJIA\ngLDtJH32LdpRwOrKVsMfxlGOfM44MjUqSO4zlWKOKQj4OzP7HWjZ/LZyZaPq65LqizYwNAiDD1Ju\nhqXLg7i0pA7PLg8eQ/TsF9GaVn2iF3HuRcVI4YeoMsM7Va7T04Cr1p+hH1oK5cOcEMvHnHE/TDcW\ndmfgP2Mg9yIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIg\nIiICIiAiIgIsVk9SUKtqtSsW4IbeRcmp1pJGaWdx7P0x83WVQEXhmy9ULUVI54xt2YZJ4azl85JD\nE4jLII/WIuYM7/8AaZe5AREQYjWmn4crj7mNsHLHDka515ZK7gMwAbbOUbyAQMf/AIhdv2LI0a4w\nxRxC7uMEYRi5bcnEBYWctmZt9m+pl3Ig1PxB0DUzJVpZZ71K1j3l9FyGKuHTuRBKwNPEMos7dM2A\nWdnbft2dt331rOeBuItQX6/XyUAZijTqXTitRTTSlVv/ACgN2SxdhllmvHJ7BySEbODMzMLsztKK\nIMJrjS9XMUpKFzqtFKUUjSV5XhnililCaGWGUe4SCYC/1t5s7Ozuy9GAxJVKUNN7dy09eHpNeuyR\ny3ZPNmkmlGMRklZnb2nHvxZ33fd37ZsvVC1FSOeMbdmGSeGs5fOSQxOIyyCP1iLmDO//AGmXud9u\n7/Ugh7F+AkFWIYK2p9Y14Y+TjFXzscEbERFJITRxVGBiIyIndm7kRO/d1vmN0eEE+NsPfyU54bGF\njhazbGUbouMAvcyHzbPYvfMb9VnHvLJ279uMX4g4K1a9CrZjGz2+Ti1aG9Achm2+4RiJ/OG2z7iO\n7tt3WzICIiAiLw4TL1bsZTVJ47EQTSwFJEXIWmiN45o3f/rCYkLt9Ts6D3IiICIiAiIgIiICIiAi\nIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiLpv2RhilmPdxgjOQmFnInEBcn4\ni3d32byQV+8QsdPlchqvLVx5S6Nq4yvhi9omK7Qmjz90AFvKR3aKB9u7tJt5bKTtZ6+9HxmNuY6A\nbtnUk9KvhoCk4RSSW4nsBLPKDFwgCAJJCdvPgzbty3bRPC7wjq38VDkcoWTgv547GUvwV8ldpRtJ\ndnksMJ1ozFmkaGSIS5Mzu4vuzeTYrHUMhXxVDhTt25vDnVlqM6wxP6TewzekBHPRDZmskNW5DwEf\npejGP0t0GTrXM6eq8bDk6+PiytbTeXeCeidiXF2Skkr9EwGfhPHxkjcTiJ92bi7H7fbPS+J9s9H1\nM3Xhrllsg1SlDUljmGs+XluBj5oih6nWGIJ2mfhz5cY9t3fu/gxOcky2r8VdioX69CLB5KGGzepz\nVjmleau8u8Ug8oom2Bgc9nN2l2biLEWAw2BufrUODKtM+KxepL2qo7JxO1bjPRiKrXgNh4uUV+5O\nXF33bZ9t9t2CStc5fUnpYUMFQoOQ0xsz5fNFZDHc3l6fokEVPeaSxxFzfd2EWIe5Puy0y94xZFsf\nQmjpVPlBtbhpfLUxeSeMiHrvM1CV5IuE0jDBwOTkIvI+7E3dNfzl+sVuPOUs7exhVKLaeq4YcidK\naVxl9PG9HRkCOay8vBhaxuAxi7vxZ930LF4G1Xr1YPk2Wp0vGelaalXrkUNSkVWIo+k8IdP0KIXa\nNpBZo26e3bbZBMmktYZoMz8j56pjITuY+bI0LOJs2JoWjhmjilqzvajAzmFpRLqMIDs3l37YCPxJ\n1HPWbP0sNUs6ec36NMJLZ6htVGnGH02vXCPo9x5yNA+5u2zb/W+b1bjzn1dit45HgPTuZglmGMnj\nB5JajMJSM3ETduTszvu/F/sWo6V1Zl8Vha+nK2HunqPHM2PqnNQsHhZBCXjHkpL4OIegvBsT7ExM\nW4ebboNwz+r83byVvG6cqY4/kVqr5O7m5LcUDy2IxsDSqx1o3N5mgKMnkL2ReTZ23bZ+ibxVlLEV\nrEOP3zN3MHgBxUk+8UOXimmisjLZjF2epGNeaXqM3tCwM/FyfjH2vdD4iHPZWzqPA5HKVcq9ezjM\njiwyUwxytXigtULENCZihkeWNjjc2diGR25NxXrx2kJ8dhsXlKOAmpPi9SnmpsDDPat3yxslcse5\ncLU0hPkWqjWkeASZvmnbYSd2QZnC3MuescRDm69OK7WwGTd58XJNJQtRyTVuJwNZFpopBKMxKM99\nvYdidi7bL48SHYixGGaQ4otTZmGlfKMiCQscEM1q3CEgOzg8rQBE7s/0ZT81gcTnJMtq/FXYqF+v\nQiweShhs3qc1Y5pXmrvLvFIPKKJtgYHPZzdpdm4ixFuPjDp+3cqVbOOEZMlgMjBlKNeQ2jjtlEMk\nc1M5CbaNpYJpgYvqJwd3Zt3YO7Unhvh7mOPHNQq14xi41JK1eOGWnKLfM2K0kYsUUoEzEzs/fZ2f\ndndn0Dw58U79wdBxP6PIGqKuajys8sUr2HnxVbi0kBjIIRuc8cjnyA2dn9nj5rnVni9lTDjh8Fkz\nstUstcoXsNfgmpz9IihsjcFirzRAYu3RZt5eQ8DZ24lhmp1cfW0DmcY8+XxGmo8lVtWMZVksSn6Z\nSerJb9HH5yMAtQzcwdtwcnF/abZwki7r0quXzle20Q43AYGtlSkjike0/J7ZWGJ+pxkHhXHiIgz7\nu/d9221+HxA1HWjq5fK4zGwYHITUY3irW7B5bGx3ZQgr2LjSxjBOHUmrsQRbELG/YnZxWDkx9rP5\nLVbBUtURzOlKcGOK/AdcyLe30SnAm+YJ5H3eIvbECDkIkTi2kUtB4C1FUow6JyYZ2eSCDJhfmzcG\nPpjyYb1s78dxopK7gMrxvGW58g2bf2XCXsRrHUN/M5OhSq4r0LCZivBbt2ztxSPRkhhkOOuETm0+\nQblOfIulGzNCOzuZEGO034gHVwvWr4yk97I6nyGJxmOx0PoNae092ztPZdnPpv0oZZpZf6xC/wBH\nk22f8Ja8gZXVpHHIAzZ2IoiOMhGQPk+qLnGRNtIO7O27btuzt9Siu3o2W7pyE5cZZyEWK1jlrl7E\nxtNDauY6S3cim9FFiAzmYZY5QYSbmwOzP7TbhK2kNV5mHJw4jUVbHR2MjVs28bdw01g6kzVpIms1\nZYrYtNFYCOxCfLuBM5/RdmZ9b074kakyEQZinicfYwUt56rUYp7ny+EY3ipyWpIii6IuIt1XgdmJ\nmbuWz8mw3hRpHCPna1vA6auUK+MrWZJctlTzNZ2uSxHVirVatyxxsO0c1nm5ATDu2ziXHfA58GnG\neSrgM3htcSzOPPCQ3a+PmtDZfjZmtcmo2seXLkckze1zk33Z93Cz6LrrMbADSOLyMA9RwZ2Fz2bk\n4s/dh33XYgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIi\nAiIgIiICIiAiIgIiIC+IYhBuICINyIuICwtyInMy2b63Iid3+t3d19ogIiICIiAiIgIiICIiAiIg\nIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIuqxYjj4vIYRsZsAOZiPIy\n+iA8n7k/1MyDtRdQ2o3PptIDnsT9NjFz2F2Y34779nIWf7OTfau1AWJ1DqXH4/p+m3K9YrBcYI5Z\nRaWYv+rDF9OUu7dhZ/Nd+oci1StLO+28Y7Az+TmTsIM7b925O2/7N1XfVuHxleseSzd67cIrRy2Z\nZZRjex1HBoq20AscFZulDGwQ9tidnEmduOzheF63nPoxcVxcYZ5Y9U8UdZ4uacKw3Ygszb9GrY51\nbErs+ztFDZETkdv+yz/as+vz28RNavmLHVmeCMK/VeNo4ZZY5XYhIHF5tiN9m4bvx7sReyzszWb/\nAEWPEMslXlxdq7DctY2GOWvMPNppae/RJp2kJyOSORhZzd93GeLdvrLRxfhs4cfPWfm8uF4/rW5b\nJuRazrTXFLEzU69kLcs+Uaw9OvQpT3ZpfRxiOfjFXEjfiMol5eTE/kLrKaczAXoGsBDbgEiIWiv0\n56U7cX2d3gsCMgi/1O7d27sudOO8Vi8x5S6EWiZ0ySItX1Hr7F0K2Rtzzu8WDlghyPRhkkOCWYoR\nij4sPtk/Xid+O+zF3UUpa86rG0zaKxuW0Iuq5ZCGM5ZC4xxC5GWzvsLefYW3f9zLoxWRCyLmAyiz\nbbdaE4nJnbcSHm3tA/2t9iql7EREBERARat+vNTqYmF47DSZ97HRFxi3q+jsLTemfO/Ns08kFfcO\nfztiJvJ+Sy+Hzde0DyRlxH025SHqOAvJPUsz1bAxsxPy2krTft2B32ZWmswjcMkiwOO1VBYqWrcM\ncp+gy5OI6rnWjsSHj7tzHy8OtOMIBJNRm4HKYDs4uTh7XHMjajc+m0gOexP02MXPYXZjfjvv2chZ\n/s5N9qjUwl2ovP6fBsb9aLaI3CR+qGwGzbuBvv7JM3fZ+69CgEREBERAREQEREBERAREQEREBERA\nREQEREBERAREQEREBERAREQEREBERAREQFHPjFiLNx2hjoNajmxeRg68dOjasDLK0PGqxZA+jUgk\nEHIjeKXm8Is7g4i0kjIrUtyzuEWjcaRhgtLSANawdJmujqCKzJZIAew0L1WryyvM/t8XApAfv3Yy\n+11JssggJERMIgzkRE7MIizbu5O/ZmZvrX0uDBiZxJmcSZ2IXbdnZ+zs7P5slrc3qRGvRovi5dGT\nDzT1pAnGnPBJP6OYybR9RgMi4O+wg0nUf7Bjd1UbxM1VlhyIyA/Tq0Q/2Q+gJxGMscfXaSV2fk8h\nC4OG7eyO23fd7BPoSxpTJWcjh6z28FlW45bCxCxSVx9p+pWjftLCPM9g/qjIQu3F2INY1N4WfKUM\n0mnLdW1TtRnvj7ZOFis5C/GNnP6QiTs7dTiTcB7ntu/0XAXw4qamdxPpP1ie09v2l814hTNkzc1a\n+ceUx9Jr3+P7wgUMzNanAvSKNR2OYeACYVWZiCXn0zlIREjmMtmZmd4933ce1gf0TYJpMjema9Ut\nQVaARmNaLcgksWpZYgKwM5i7iMU3IS3Ld2dncX6k2raX/R+zrSNJbfFY+ITd5JZ447PEBlrF2GKV\nubkNdy7mHaSQd257jZfw/wBOV8bUCvVF+nyKSWwQCElqYycjlIRZuIbvsLfYws3ZmcvPxDisc4+W\nk7218Dw+SL7tGmk+NOHns5bC2Pk/L2q2NqZPlNgbQVrsViy9OONhlKzCUY9OGZn4n3aV2223368r\niLVyrh6HyXaLHTFcK5LqSGDO5KrM1nauM0Vm3JFsfM5BnMp2jjAB4h2ZpeRcqvF2ita6/T8/j/Lr\ndKNzPdX/AMJ9BWSuYdspi5hrYzGZecRvhUkihszZwSoVOjF82EkdeAbDNGAxiRA4MPFmbzXfD+5J\njtUUnxJxz5vVwHHZrxVI5Cwb5KpPEUEobjxAYZ5OBN2eRyJnInZ7EovSfEcnNzaj/vff1U/p661+\nemkWUtHNU1LkrlbHPDjLGn4K9itTihhiyV+S/LPNaNhkAZLUUT7czZzfrnsW+7FnSo3zovFxn6UV\n0HjhkeErUlEfpRlyd4zdndthLz4Ozs7OzPuyLJly2yTE27a9ntWsV9GmngWcK0UcdkoflQLEo2Ag\nEQEIi5cYohEQiIiZttm/rdu+7/edwcvWOGtHwrZSOvHYKPiI12gJhJxDbYWKDYWZu24v27s7bei8\n1mmXMER2puqNpg5Q+hSVBqu0YCIs0YFIDnW4k39VxEm3337rK6dx3GxdsSQ8JJLkjQkYi5dFgjHm\nBeYibs/7+Lb+SzyIIksacvCWcf0WeQMTPFdwgjJ87enfLHqa3HATy7dM5/Q6rCXFmeuY/Qdljsvo\nvJGFVpY5ZDs4CpW51KdCc6GYezYt5G1HPel50TknnqzDKDGzFixcvaGITmxF7RmtCnJCv+o9FX5s\nVla1bFSQXjsa9nltdOoz5Kpknz8dCkEwydUnmltYuZhJmZmxwcuPzfLeZNMGMh2o6IPb/WaS0E7i\nASvXKA6/PrsznHEUZlG7t/VMuz7vvI6KJzTKeWEIYHSNqSzUOfFvDCUWIC1X+T8fUqAVWzYMoY69\neWSQ4YeqTM80sjbE7g+xO5TcLMzMzNszNszN5M32MuUVb3m8+ZWsVERFRYREQEREBERAREQEREBE\nRAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQFhsjpbH2JOrJVj6rvu80fKGR3+0jicS\nJ/37rMorVvak7rOlbUreNWjbF08BUidiGLkQ/RKaSSdx+xx6xFxf9rLKIiWva36p2VpWvlWNCIiq\nsIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICI\niAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgI\niICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICIiAiIgIiICLDaj1PSx5QhakkArQylCMVW\n1ZchieJpHdq0R8GZ5om9rbfn237rA1vFjASyHFFeKWWH+lhio35JY/q+cjCu5B/azIbbui1L/wCo\n+I++tfwnKf8AtV5JvFrT4TDXO+4WJePTrHRvjOfLlx4QlX5lvxLbZu/F/sdEbhvCLUv/AKj4j761\n/Ccp/wC1Ws5rx903Unkryz2+pDx5s2OtD9IBkbtIAk3sm3mzf3bOidpTRRB6x2l/v7f8Pn/knrHa\nX+/t/wAPn/kgl9FEHrHaX+/t/wAPn/knrHaX+/t/w+f+SCX0UQesdpf7+3/D5/5J6x2l/v7f8Pn/\nAJIJfRRB6x2l/v7f8Pn/AJJ6x2l/v7f8Pn/kgl9FEHrHaX+/t/w+f+Sesdpf7+3/AA+f+SCX0UQe\nsdpf7+3/AA+f+Sesdpf7+3/D5/5IJfRRB6x2l/v7f8Pn/knrHaX+/t/w+f8Akgl9FEHrHaX+/t/w\n+f8AknrHaX+/t/w+f+SCX0UQesdpf7+3/D5/5J6x2l/v7f8AD5/5IJfRRB6x2l/v7f8AD5/5J6x2\nl/v7f8Pn/kgl9FEHrHaX+/t/w+f+Sesdpf7+3/D5/wCSCX0UQesdpf7+3/D5/wCSesdpf7+3/D5/\n5IJfRUA9dXVX5fp/3PI/EU9dXVX5fp/3PI/EUF/0VAPXV1V+X6f9zyPxFPXV1V+X6f8Ac8j8RQX/\nAEVAPXV1V+X6f9zyPxFPXV1V+X6f9zyPxFBf9FQD11dVfl+n/c8j8RT11dVfl+n/AHPI/EUF/wBF\nQD11dVfl+n/c8j8RT11dVfl+n/c8j8RQX/RUA9dXVX5fp/3PI/EU9dXVX5fp/wBzyPxFBf8ARUA9\ndXVX5fp/3PI/EU9dXVX5fp/3PI/EUF/0VAPXV1V+X6f9zyPxFPXV1V+X6f8Ac8j8RQX/AEVAPXV1\nV+X6f9zyPxFPXV1V+X6f9zyPxFBf9FQD11dVfl+n/c8j8RT11dVfl+n/AHPI/EUF/wBFQD11dVfl\n+n/c8j8RT11dVfl+n/c8j8RQX/RUA9dXVX5fp/3PI/EU9dXVX5fp/wBzyPxFBf8ARUA9dXVX5fp/\n3PI/EU9dXVX5fp/3PI/EUF/0VAPXV1V+X6f9zyPxFPXV1V+X6f8Ac8j8RQX/AEVAPXV1V+X6f9zy\nPxFPXV1V+X6f9zyPxFBZL9J+5DHLh4pp2hjtPMFkRn6ViTGtksEWY9GYHaUjHHtbJ3i9sREiHZ2Z\n1q2swtxYSM71fCVwDGRPpqTDxVI2+X2Gd4QxpQ2HP5OImq9Pgz71Hvek8BZVP8afHPMarkxst6Kj\nUkwZ2JKcmLC1XNpJ3rERkctmQmMSqROLg4uz8vPttp9zWuRlCUHkgj9KEgsyVqFGpPYAn3kCxZqw\nBNOBv3JjJ2L690Quy+tMM/dsvinZ/LbJ09tvq/3q1uTWmMicIPT8dPD6PqOO3T+WMWNXLT5SxXPG\nz2pZZnkjjrRBwkd/nAalA0ASDsCpXz/YP9y2GhrfIwY2XEwyjHSsnKcoDDH1C6rQtIPW26gg7QCz\nsztuxGz7s+yjSIjS6dfWGIEAEs1izIQESkfJ02cyZmZzdur2d37/ANqhDxPvQ2MrdnrSxWIjeDpz\nV5QmiNxp1xJhkjdxLYhIX2fs7O31KvXP9g/3LNYrU9itEMMYQuIOTs5ibv7ROT+Rs31/YmiI0sE+\nnYP9p41Tcqsox4kCe2TZ+EpXjO5GMT9SzFFGw2Sels3TlDy/pC8ep8BHWiN468pMMMc81ppJOFSz\nJYjjbHSdYQdnaCQZXhMGtN14CJgAJWKEf1zs9/mq3fz+bk7/AL/nO6NrOz91W7dv6OTy/wARSsna\njpiGWs5PFPFIVGtNXmETkKQ5bONjnkljMhinhD0qZ2OoRtFG8Y2Bjm2Y/m3oQoimjksSBPXdyKA6\nYcyjavhZ+kDhaJpMkZZqGGOuHIJJI9hm2kAngv8AXOz91W/w5P8AMRtZWfuq3+HJ+3/vP2v/AHug\nsDm/DrpVxttYOtA9Wi5Nfq+jvHZsPUrh1yecmrw9aSY5nf2q4izCE/Z1jtT6P9FxsF4PTgIhIJoL\nuNepIJ+kzgMkzFZJoXcGiFoh5u7CxvxY+0HtrKy3lFWb90cn+Yg6zst5RVm/dHJ/mILA5TSePZyC\nGywn6c9TmMj2o2isM2VxluvHCzzWXkxoy1hjZvnrPEQ8nMuX8P43jgcLJ7SPVhOcqE4hIVjUE2L6\nkAzTRk1qCIY2nrGLdOR2Dk7m0gV+/XS191W+v/dyfX5/7xcfrlZ+6rf4cn2bfefZ2QTjJojatFce\nxK1aePqDI9IOTsWHu5mMGBrbt1unTYDEiERewDgU3GRh1WvDyEn2kfj93FzFvmppPaLk3H+i38n9\nlpS/qbFHH652fPpVt/t6cn27/efauf11tfd1/wDySf5iCadI6fjtVbU7CVmeB5BenH7JRwjD1Gs9\nXrg0RSSu0Qyk0wxOJEUE4uTw5PMaRihG4xV7MR1GsdMzGZoiqBVyVmndI5G29KlkpU4pRHeEflBg\nZo5jaOtAb60tP5xVn2+0JP8AMXDazst5RVv8OT/MQThm8JXjpSz1oTdheqfXsjO7xxS0MJOMccsU\nrQemHLkJieI4z3iCYhePpD1OKWjJZMedkhlilrwXbxM8bGEtCPHVb1eZn6jOMZbzxsYDI/UNhJga\nOUo4Q/XS159Ktv8Ab05P8xP1zs/dVvPf+jk8/t/pPNBPWO0KL3q1eSaQ45L0MEr+ivGLR+lNCRTS\nNM/o8dodxpnsXpRkw7Reaw2pdMlTqx2Odg3kKISCWg9ZoXlrlOEUxlOW11uEgnADH03iLmQ8ourD\nz6zs/dVvr/3cn1+f+8RtZ2W8oqzfV2jk8v8AEQTDrLEBUKVgjkBwyGS6Yyc9yxsZVAoWWYm9qvIU\nk4jN9E3F2Z34rMt4euRsw2Z3ifJ2Ma83yezkU8MORlPpQx2iI35Y022Nw9mxETOW0oxQTHrq4ImA\njCIS8eoAtKwnxfceYtLsWz9238l1/rlZ+6reW39HJ5fZ/SeXZkE/YPw9O3XscCnaxEdfpSlUZ6Rc\nxssUDWY7BNIQSdAJyAC9HOKUWaZtyXS+hY3mCCO5YkN7UsEjfJkbbNG+aj5ws193mJ5MHZ9h+Hsy\nxOzkXKMYHfWVl/8AdVv8OT/MT9c7P3Vb/Dk/Z/3n7G/uQTVmtKBSs1K9m0zda6NS9L0xihqv0cZY\nmkimkkfrQDHkm+ckCL+hJ3Fmde2fSkb12t2AsYvj1mOpHSszysMVG5kXIo8nNXNpejSlDkBFEcji\nwvHtI0cEjrS03lFWb90cn+Z+1cNrOy3lFW/w5P8AM/Y39yCfMlomGGnamOconpwTRvYs1rAVZJ49\nQ28b1oTr9WSRmgqsBxDG/TKzA7k/J3bVdT1Ar3r1eP8Ao6t+3BFu7v8ANxWJIw7kzO/si3d2Z1Fr\nazst36Vb9/Tk+zb7z7Fy2tbX3df/AMkn+Yg1lERAREQEREBERAREQEREBERAREQEREBERAREQERE\nBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERARE\nQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBE\nRAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQEREBERAREQE\nREBERAREQEREBERAREQEREBERAREQEREH//Z\n",
      "text/html": [
       "\n",
       "        <iframe\n",
       "            width=\"700\"\n",
       "            height=\"400\"\n",
       "            src=\"https://www.youtube.com/embed/m1q2AIoQADs\"\n",
       "            frameborder=\"0\"\n",
       "            allowfullscreen\n",
       "        ></iframe>\n",
       "        "
      ],
      "text/plain": [
       "<IPython.lib.display.YouTubeVideo at 0x108a5db70>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "YouTubeVideo('m1q2AIoQADs', width=700, height=400)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "***"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "&copy; 2017 ECMWF\n",
    "\n",
    "This software is licensed under the terms of the Apache Licence Version 2.0 which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. In applying this licence, ECMWF does not waive the privileges and immunities granted to it by virtue of its status as an intergovernmental organisation nor\n",
    "does it submit to any jurisdiction."
   ]
  }
 ],
 "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.6.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}