notebooks/Roskilde sources.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Sources in Roskilde Fjord\n",
    "\n",
    "1. Extract coordinates of all sources from pfs file\n",
    "2. Convert from local projected coordinate system to geographical coordinates\n",
    "3. Plot the location of each source on a map"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read the PFS file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from pfs2yaml import pfs2dict, pfs2yaml\n",
    "\n",
    "pfsfilename = \"RoskildeFjord_c9_2011_coarse.m3fm\"\n",
    "\n",
    "with open(pfsfilename) as f:\n",
    "        pfs = f.read()\n",
    "\n",
    "y = pfs2yaml(pfs)\n",
    "d = pfs2dict(pfs)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get projection"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 18,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'PROJCS[\"ETRS_1989_Kp2000_Zealand\",GEOGCS[\"GCS_ETRS_1989\",DATUM[\"D_ETRS_1989\",SPHEROID[\"GRS_1980\",6378137.0,298.257222101]],PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",500000.0],PARAMETER[\"False_Northing\",0.0],PARAMETER[\"Central_Meridian\",12.0],PARAMETER[\"Scale_Factor\",0.99995],PARAMETER[\"Latitude_Of_Origin\",0.0],UNIT[\"Meter\",1.0]]'"
      ]
     },
     "execution_count": 18,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "wkt = d[\"FemEngineHD\"][\"DOMAIN\"]['coordinate_type']\n",
    "wkt"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<CRS: PROJCS[\"ETRS_1989_Kp2000_Zealand\",GEOGCS[\"GCS_ETRS ...>\n",
       "Name: ETRS89 / Kp2000 Zealand\n",
       "Ellipsoid:\n",
       "- semi_major_metre: 6378137.00\n",
       "- semi_minor_metre: 6356752.31\n",
       "- inverse_flattening: 298.26\n",
       "Area of Use:\n",
       "- UNDEFINED\n",
       "Prime Meridian:\n",
       "- longitude: 0.0000\n",
       "- unit_name: Degree\n",
       "- unit_conversion_factor: 0.01745329\n",
       "Axis Info:\n",
       "- Easting[E] (east) EPSG:9001 (metre)\n",
       "- Northing[N] (north) EPSG:9001 (metre)"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pyproj\n",
    "\n",
    "inProj = pyproj.CRS.from_string(wkt)\n",
    "inProj"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Read source data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "metadata": {},
   "outputs": [],
   "source": [
    "from collections import namedtuple\n",
    "\n",
    "Source = namedtuple('Source',['name','lon','lat'])\n",
    "\n",
    "outProj = pyproj.Proj(init='epsg:4326')\n",
    "sources = d[\"FemEngineHD\"][\"HYDRODYNAMIC_MODULE\"][\"SOURCES\"]\n",
    "\n",
    "n_sources = sources['number_of_sources']\n",
    "\n",
    "out = []\n",
    "\n",
    "for i in range(n_sources):\n",
    "\n",
    "    key = f\"SOURCE_{i+1}\"\n",
    "    \n",
    "    source = sources[key]\n",
    "    name = source['name']\n",
    "    x,y,z = source['coordinates']\n",
    "    \n",
    "    lon,lat = pyproj.transform(inProj,outProj,x,y)\n",
    "    s = Source(name=name,lon=lon,lat=lat)\n",
    "    out.append(s)\n",
    "    \n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 45,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[Source(name='Arresø Kanal', lon=12.00640530070797, lat=55.96360022598667),\n",
       " Source(name='Havelse Å', lon=12.057519866061158, lat=55.877360558793754)]"
      ]
     },
     "execution_count": 45,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "out[0:2]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Create an interactive map\n",
    "\n",
    "If the map does not render below, take a look at the [prebaked version](map.html)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div style=\"width:100%;\"><div style=\"position:relative;width:100%;height:0;padding-bottom:60%;\"><iframe src=\"data:text/html;charset=utf-8;base64,PCFET0NUWVBFIGh0bWw+CjxoZWFkPiAgICAKICAgIDxtZXRhIGh0dHAtZXF1aXY9ImNvbnRlbnQtdHlwZSIgY29udGVudD0idGV4dC9odG1sOyBjaGFyc2V0PVVURi04IiAvPgogICAgPHNjcmlwdD5MX1BSRUZFUl9DQU5WQVM9ZmFsc2U7IExfTk9fVE9VQ0g9ZmFsc2U7IExfRElTQUJMRV8zRD1mYWxzZTs8L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS4zLjQvZGlzdC9sZWFmbGV0LmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2FqYXguZ29vZ2xlYXBpcy5jb20vYWpheC9saWJzL2pxdWVyeS8xLjExLjEvanF1ZXJ5Lm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvanMvYm9vdHN0cmFwLm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuanMiPjwvc2NyaXB0PgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS4zLjQvZGlzdC9sZWFmbGV0LmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9jc3MvYm9vdHN0cmFwLm1pbi5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvY3NzL2Jvb3RzdHJhcC10aGVtZS5taW4uY3NzIi8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vZm9udC1hd2Vzb21lLzQuNi4zL2Nzcy9mb250LWF3ZXNvbWUubWluLmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9MZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy8yLjAuMi9sZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9yYXdjZG4uZ2l0aGFjay5jb20vcHl0aG9uLXZpc3VhbGl6YXRpb24vZm9saXVtL21hc3Rlci9mb2xpdW0vdGVtcGxhdGVzL2xlYWZsZXQuYXdlc29tZS5yb3RhdGUuY3NzIi8+CiAgICA8c3R5bGU+aHRtbCwgYm9keSB7d2lkdGg6IDEwMCU7aGVpZ2h0OiAxMDAlO21hcmdpbjogMDtwYWRkaW5nOiAwO308L3N0eWxlPgogICAgPHN0eWxlPiNtYXAge3Bvc2l0aW9uOmFic29sdXRlO3RvcDowO2JvdHRvbTowO3JpZ2h0OjA7bGVmdDowO308L3N0eWxlPgogICAgCiAgICA8bWV0YSBuYW1lPSJ2aWV3cG9ydCIgY29udGVudD0id2lkdGg9ZGV2aWNlLXdpZHRoLAogICAgICAgIGluaXRpYWwtc2NhbGU9MS4wLCBtYXhpbXVtLXNjYWxlPTEuMCwgdXNlci1zY2FsYWJsZT1ubyIgLz4KICAgIDxzdHlsZT4jbWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1IHsKICAgICAgICBwb3NpdGlvbjogcmVsYXRpdmU7CiAgICAgICAgd2lkdGg6IDEwMC4wJTsKICAgICAgICBoZWlnaHQ6IDEwMC4wJTsKICAgICAgICBsZWZ0OiAwLjAlOwogICAgICAgIHRvcDogMC4wJTsKICAgICAgICB9CiAgICA8L3N0eWxlPgo8L2hlYWQ+Cjxib2R5PiAgICAKICAgIAogICAgPGRpdiBjbGFzcz0iZm9saXVtLW1hcCIgaWQ9Im1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSIgPjwvZGl2Pgo8L2JvZHk+CjxzY3JpcHQ+ICAgIAogICAgCiAgICAKICAgICAgICB2YXIgYm91bmRzID0gbnVsbDsKICAgIAoKICAgIHZhciBtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUgPSBMLm1hcCgKICAgICAgICAnbWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1JywgewogICAgICAgIGNlbnRlcjogWzU1LjgsIDEyXSwKICAgICAgICB6b29tOiAxMCwKICAgICAgICBtYXhCb3VuZHM6IGJvdW5kcywKICAgICAgICBsYXllcnM6IFtdLAogICAgICAgIHdvcmxkQ29weUp1bXA6IGZhbHNlLAogICAgICAgIGNyczogTC5DUlMuRVBTRzM4NTcsCiAgICAgICAgem9vbUNvbnRyb2w6IHRydWUsCiAgICAgICAgfSk7CgogICAgCiAgICAKICAgIHZhciB0aWxlX2xheWVyXzM1NTQ1NmFmZmI4NzQ3ZDI4NTAzYTQ1ZmIwZTJmMTY4ID0gTC50aWxlTGF5ZXIoCiAgICAgICAgJ2h0dHBzOi8ve3N9LnRpbGUub3BlbnN0cmVldG1hcC5vcmcve3p9L3t4fS97eX0ucG5nJywKICAgICAgICB7CiAgICAgICAgImF0dHJpYnV0aW9uIjogbnVsbCwKICAgICAgICAiZGV0ZWN0UmV0aW5hIjogZmFsc2UsCiAgICAgICAgIm1heE5hdGl2ZVpvb20iOiAxOCwKICAgICAgICAibWF4Wm9vbSI6IDE4LAogICAgICAgICJtaW5ab29tIjogMCwKICAgICAgICAibm9XcmFwIjogZmFsc2UsCiAgICAgICAgIm9wYWNpdHkiOiAxLAogICAgICAgICJzdWJkb21haW5zIjogImFiYyIsCiAgICAgICAgInRtcyI6IGZhbHNlCn0pLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAKICAgICAgICB2YXIgbWFya2VyX2MwY2MyNjQ3NWM5MjQ3MjQ5N2MxOTJhMmM5MzlhZmJhID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS45NjM2MDAyMjU5ODY2NywgMTIuMDA2NDA1MzAwNzA3OTddLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfYzBjYzI2NDc1YzkyNDcyNDk3YzE5MmEyYzkzOWFmYmEuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgQXJyZXPDuCBLYW5hbGAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyXzM0MzRlNTI1MTVhNTQ2Yzk5ZThjYzVlMGRiMzhkYjdlID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS44NzczNjA1NTg3OTM3NTQsIDEyLjA1NzUxOTg2NjA2MTE1OF0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl8zNDM0ZTUyNTE1YTU0NmM5OWU4Y2M1ZTBkYjM4ZGI3ZS5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBIYXZlbHNlIMOFYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfNjI3YmI3OGYxOTk5NGU5NmIwZTAwNmUwMDBjOWEyMTggPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1Ljg2Mzg4NzU4MTIxMzIxLCAxMi4wNTc0OTk5NTY0NjA0NzldLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfNjI3YmI3OGYxOTk5NGU5NmIwZTAwNmUwMDBjOWEyMTguYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgR3LDpnNlIMOFYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfYTRhMjY2ZDViNDVlNGUxYWIxZjQ5NTZhOTUxNjMwNWEgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1Ljc4MDc5MDA4MjAxMDMyLCAxMi4wODEyODQ4NDU2MjQzMV0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl9hNGEyNjZkNWI0NWU0ZTFhYjFmNDk1NmE5NTE2MzA1YS5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBWw6ZyZWJybyDDhWAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2NhODhjNDEzZDU1ZTQ0OGU5YzBiNDQyNjkzOGRiNDg0ID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS43MDM1NzgzODE1MTE4NiwgMTUuMTAzNzcyNzQ0NjY5Nzk3XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyX2NhODhjNDEzZDU1ZTQ0OGU5YzBiNDQyNjkzOGRiNDg0LmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYEhvdmUgw4VgICsgYDwvZGl2PmAsCiAgICAgICAgICAgIHsic3RpY2t5IjogdHJ1ZX0KICAgICAgICApOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIG1hcmtlcl85NmQ4MDYxMWY4Mjk0NmQ0YjZiOWJmZGQ3NTM2MzE5NSA9IEwubWFya2VyKAogICAgICAgICAgICBbNTUuNjg5MDIxNDQ0ODYxNTQ1LCAxNS4xMDQ0MDAzMTM1MDc4NjldLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfOTZkODA2MTFmODI5NDZkNGI2YjliZmRkNzUzNjMxOTUuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgTWFnbGVtb3NlIMOFYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfODVlNWQwZjU3MzA1NDQzOTk3OWQyMTYzYjdkZmU4OTEgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjY3Mzk4NzUxMzUyMTM5LCAxMi4wMTU1MDg3ODUzMDMxMDVdLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfODVlNWQwZjU3MzA1NDQzOTk3OWQyMTYzYjdkZmU4OTEuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgTGFuZ3ZhZCDDhWAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2E5MTQ2ZGVkZjcxNDQxOGE5ZDRjZTg4ZGY3MmVlOTgwID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS44MjcwNjA1MzU1NTM5MjYsIDEyLjA1OTA0MTMzNDEyMTY1OV0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl9hOTE0NmRlZGY3MTQ0MThhOWQ0Y2U4OGRmNzJlZTk4MC5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBTaWxsZWJybyDDhWAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2FlODU5OTk1MTc4ZjQ2YmI5NTg3MGRmMzQ5M2EyYjMxID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS44MDU2ODY3ODkyMzkzNiwgMTIuMDUxMDM0Nzk1NDU4NjUxXSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyX2FlODU5OTk1MTc4ZjQ2YmI5NTg3MGRmMzQ5M2EyYjMxLmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYE1hZGVtb3NlIMOFYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfZDk2ZThhYzljMTBhNDcwN2JmMmQyNjlmNDdjOGE1MDkgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1Ljk1NzMwMzczMzE5MDE2NiwgMTEuOTUxOTY4MDQwNzgxNTg5XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyX2Q5NmU4YWM5YzEwYTQ3MDdiZjJkMjY5ZjQ3YzhhNTA5LmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYE1lbGJ5IGNhdGNobWVudGAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2FmODliNDAyMjAxNzQ0Y2Y5MTYwNjNmODY1ZTlmZjc2ID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS42NjUzODU3MTI3Njg0OSwgMTEuOTUwNzM3MDU2MzI4ODIzXSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyX2FmODliNDAyMjAxNzQ0Y2Y5MTYwNjNmODY1ZTlmZjc2LmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYExlanJlIMOFYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfYzVlZDhkYTNjMzE0NDYyYzg0NDg3ZDQ5Y2NjOTkzYjMgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjY0OTIwOTkxNzM2OSwgMTEuOTM0Mzk2MDk5MTczNDc4XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyX2M1ZWQ4ZGEzYzMxNDQ2MmM4NDQ4N2Q0OWNjYzk5M2IzLmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYExlanJlcmVuZGVuYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfNWU4MDFiYTExZDQyNGNkZmFhODAxOWU0MjI4NWZlYjQgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjY1OTA4OTA2NDMzMjUsIDExLjkzMTgzNzM4Mjc5MDIxNF0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl81ZTgwMWJhMTFkNDI0Y2RmYWE4MDE5ZTQyMjg1ZmViNC5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBMdW5kYnkgY2F0Y2htZW50YCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfOTM2NDM3NGY0YWJkNDVmYzg3YjhhYWQ3MDc1YTlhMTQgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjY5Mjc4NTA4OTc2MTczLCAxMS45NjAyNDQwNDEyNTY1NF0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl85MzY0Mzc0ZjRhYmQ0NWZjODdiOGFhZDcwNzVhOWExNC5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGDDmHJiw6ZrYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfNjMxMzA3OGYyNzE4NGMxZThjYWY2ZDEwMzg3YzkxYTggPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjcyODI3MDQ3ODczNDcsIDExLjk4NTY3NDg3Nzc1NTk1OF0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl82MzEzMDc4ZjI3MTg0YzFlOGNhZjZkMTAzODdjOTFhOC5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBIw7huZXBpbHNncsO4ZnRgICsgYDwvZGl2PmAsCiAgICAgICAgICAgIHsic3RpY2t5IjogdHJ1ZX0KICAgICAgICApOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIG1hcmtlcl8wNzNlMzUxNDY0MmQ0ZDk0YmVmMzQ5ODFhYThjZTJiMSA9IEwubWFya2VyKAogICAgICAgICAgICBbNTUuNzM5MDQ5NTcxNjk2MzEsIDExLjk5MDQ0NzI4NjY0NzQ3N10sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl8wNzNlMzUxNDY0MmQ0ZDk0YmVmMzQ5ODFhYThjZTJiMS5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBTZWxzw7hgICsgYDwvZGl2PmAsCiAgICAgICAgICAgIHsic3RpY2t5IjogdHJ1ZX0KICAgICAgICApOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIG1hcmtlcl8wY2VjMTkyNDZjN2Y0YWY1YTAzOWUyMDc5YzNkYmYzMiA9IEwubWFya2VyKAogICAgICAgICAgICBbNTUuNzMwMDU4NDczNjQxNjIsIDEyLjA0Nzc1MjU5MTUyNTc4NF0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl8wY2VjMTkyNDZjN2Y0YWY1YTAzOWUyMDc5YzNkYmYzMi5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBTw7huZGVyYnkgY2F0Y2htZW50YCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfZmUwNjBlZDRiNzdjNDlkYWE4OTYzYzEyZDk2ODBlNzkgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1Ljg4MDk2NjUyNjg2MTg3LCAxMi4wMDc5ODk2MTAyOTYyOTRdLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfZmUwNjBlZDRiNzdjNDlkYWE4OTYzYzEyZDk2ODBlNzkuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgSG9ybnNoZXJyZWQgTiBjYXRjaG1lbnRgICsgYDwvZGl2PmAsCiAgICAgICAgICAgIHsic3RpY2t5IjogdHJ1ZX0KICAgICAgICApOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIG1hcmtlcl9kOTQ2OTRmMTVjZDQ0YjE4YmQxZDljMWU0MDRlZjNjNyA9IEwubWFya2VyKAogICAgICAgICAgICBbNTUuOTE1MDkwNjI0ODUwODI2LCAxMi4wNDMxODE3ODI4MjA5NjhdLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfZDk0Njk0ZjE1Y2Q0NGIxOGJkMWQ5YzFlNDA0ZWYzYzcuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgw5hsc3RlZCBXV1RQYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfMjZhODVkOGE2MTZmNDFlNGI1MGEyZTJkNjViNzE2YzIgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjgyNzA2MDUzNTU1MzkyNiwgMTIuMDU5MDQxMzM0MTIxNjU5XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyXzI2YTg1ZDhhNjE2ZjQxZTRiNTBhMmUyZDY1YjcxNmMyLmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYEZyZWRlcmlrc3N1bmQgV1dUUGAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyX2FkNjI1OGRkYWNiNzQwNWFiNDA0NGFjMjI0ZjQzYjEwID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS43MzQ1MTc5MTk0MTk2NDUsIDEyLjEwMDI5MTg4ODUwNzE1NF0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl9hZDYyNThkZGFjYjc0MDVhYjQwNDRhYzIyNGY0M2IxMC5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBKeWxsaW5nZSBXV1RQYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfMzdjNzZiZTQwY2U5NDAzODkwMmU1MWVlYTJjNDE4ZDQgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjY5OTQ5MTEzMzA0NTQwNiwgMTIuMDk1NDMwNjQ5MzE5MDc0XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyXzM3Yzc2YmU0MGNlOTQwMzg5MDJlNTFlZWEyYzQxOGQ0LmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYFJpc8O4IFdXVFBgICsgYDwvZGl2PmAsCiAgICAgICAgICAgIHsic3RpY2t5IjogdHJ1ZX0KICAgICAgICApOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIG1hcmtlcl9kNWVhYWRlMjc3MjU0YzA5OWQxOGI1YTYxYjA5NzA0MSA9IEwubWFya2VyKAogICAgICAgICAgICBbNTUuNjU5MDkxNTM2NjA1Mzk1LCAxMi4wNjM1NTQ4OTExMzYwNzddLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfZDVlYWFkZTI3NzI1NGMwOTlkMThiNWE2MWIwOTcwNDEuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgQmplcmdtYXJrZW4gV1dUUGAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyXzQzYjBhYjFiY2Y4NDRkODJhM2Y2OTY4MjM1ZGUyY2UwID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS42NTYwODAzMjc2NTA2MSwgMTQuOTQxMjEzOTczMjY0NzQ1XSwKICAgICAgICAgICAgewogICAgICAgICAgICAgICAgaWNvbjogbmV3IEwuSWNvbi5EZWZhdWx0KCkKICAgICAgICAgICAgICAgIH0KICAgICAgICAgICAgKS5hZGRUbyhtYXBfZTQ3MDI4MWVlMzc5NGY2MWFiODIzY2FjNWQzZDAyYjUpOwogICAgICAgIAogICAgCiAgICAgICAgbWFya2VyXzQzYjBhYjFiY2Y4NDRkODJhM2Y2OTY4MjM1ZGUyY2UwLmJpbmRUb29sdGlwKAogICAgICAgICAgICBgPGRpdj5gCiAgICAgICAgICAgICsgYEx1bmRieSBXV1RQYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKICAgIAogICAgICAgIHZhciBtYXJrZXJfMjNmZTM2NDIxYzU5NDY5ZWE2ZGU3N2IzYjZlY2NkMjIgPSBMLm1hcmtlcigKICAgICAgICAgICAgWzU1LjY4Mjk2NDc3ODk3Nzk1NCwgMTQuOTY5MzM0ODY2NTk4NTZdLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfMjNmZTM2NDIxYzU5NDY5ZWE2ZGU3N2IzYjZlY2NkMjIuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgU8OmYnkgV1dURWAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyXzM3NTI0YTg1ZjE1ODRmYWVhZWUxNzU0YmEyM2E0ODIwID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS43NjA2MDAzNDQ3MjAwMjUsIDE1LjA1NzU4ODQwNTQ3NTQyNV0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl8zNzUyNGE4NWYxNTg0ZmFlYWVlMTc1NGJhMjNhNDgyMC5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBIeWxsaW5nZXJpaXMgV1dUUGAgKyBgPC9kaXY+YCwKICAgICAgICAgICAgeyJzdGlja3kiOiB0cnVlfQogICAgICAgICk7CiAgICAgICAgCiAgICAKICAgICAgICB2YXIgbWFya2VyXzY1MDRmMWI0NzFjNTQ1MTQ5NzhhZGQyMDhkY2M1ZDdjID0gTC5tYXJrZXIoCiAgICAgICAgICAgIFs1NS44NDAyNjMzODg3MzU5LCAxNS4wMDY0MTIwMzU3OTQyNjFdLAogICAgICAgICAgICB7CiAgICAgICAgICAgICAgICBpY29uOiBuZXcgTC5JY29uLkRlZmF1bHQoKQogICAgICAgICAgICAgICAgfQogICAgICAgICAgICApLmFkZFRvKG1hcF9lNDcwMjgxZWUzNzk0ZjYxYWI4MjNjYWM1ZDNkMDJiNSk7CiAgICAgICAgCiAgICAKICAgICAgICBtYXJrZXJfNjUwNGYxYjQ3MWM1NDUxNDk3OGFkZDIwOGRjYzVkN2MuYmluZFRvb2x0aXAoCiAgICAgICAgICAgIGA8ZGl2PmAKICAgICAgICAgICAgKyBgTnIuRHLDpWJ5IFdXVFBgICsgYDwvZGl2PmAsCiAgICAgICAgICAgIHsic3RpY2t5IjogdHJ1ZX0KICAgICAgICApOwogICAgICAgIAogICAgCiAgICAgICAgdmFyIG1hcmtlcl81MWJmYzRiZDdmMmU0MmI0OTViNzk0YmEzODVkYWRmNSA9IEwubWFya2VyKAogICAgICAgICAgICBbNTUuNzk5ODkyNzEzMDYyNzYsIDE1LjAzOTAzNTM4MTE2NTE5Ml0sCiAgICAgICAgICAgIHsKICAgICAgICAgICAgICAgIGljb246IG5ldyBMLkljb24uRGVmYXVsdCgpCiAgICAgICAgICAgICAgICB9CiAgICAgICAgICAgICkuYWRkVG8obWFwX2U0NzAyODFlZTM3OTRmNjFhYjgyM2NhYzVkM2QwMmI1KTsKICAgICAgICAKICAgIAogICAgICAgIG1hcmtlcl81MWJmYzRiZDdmMmU0MmI0OTViNzk0YmEzODVkYWRmNS5iaW5kVG9vbHRpcCgKICAgICAgICAgICAgYDxkaXY+YAogICAgICAgICAgICArIGBUw7hyc2xldiBXV1RQYCArIGA8L2Rpdj5gLAogICAgICAgICAgICB7InN0aWNreSI6IHRydWV9CiAgICAgICAgKTsKICAgICAgICAKPC9zY3JpcHQ+\" style=\"position:absolute;width:100%;height:100%;left:0;top:0;border:none !important;\" allowfullscreen webkitallowfullscreen mozallowfullscreen></iframe></div></div>"
      ],
      "text/plain": [
       "<folium.folium.Map at 0x256bfc549e8>"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import folium\n",
    "m = folium.Map(location=[55.8, 12])\n",
    "\n",
    "for s in out:\n",
    "    folium.Marker([s.lat, s.lon], tooltip=s.name).add_to(m)\n",
    "\n",
    "m"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [],
   "source": [
    "m.save(\"map.html\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3.6",
   "language": "python",
   "name": "py36"
  },
  "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.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}