viewer.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# LOFAR single station imaging"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Based on version 1.4 of V.A. Moss's tutorial, tutorial from M.A. Brentjens 8/8/2017.\n",
    "\n",
    "Updated February 2020 by T.J. Dijkema with aerial imagery backdrops."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "from lofarimaging import make_xst_plots, read_acm_cube, get_station_type\n",
    "import os\n",
    "import datetime\n",
    "from glob import glob\n",
    "from tqdm import tqdm_notebook as tqdm"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "station_name = \"DE603\"\n",
    "station_type = get_station_type(station_name)\n",
    "xst_filename = 'test/20170720_095816_mode_3_xst_sb297.dat'\n",
    "xst_data = read_acm_cube(xst_filename, station_type)[0]\n",
    "obsdatestr, obstimestr, _, rcu_mode, _, subbandname = os.path.basename(xst_filename).rstrip(\".dat\").split(\"_\")\n",
    "obstime = datetime.datetime.strptime(obsdatestr + \":\" + obstimestr, '%Y%m%d:%H%M%S')\n",
    "subband = int(subbandname[2:])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Maximum at -6m east, 70m north of station center (lat/long 50.97998, 11.71118)\n"
     ]
    }
   ],
   "source": [
    "sky_fig, ground_fig, leaflet_map = make_xst_plots(xst_data, station_name, obstime, subband, rcu_mode, caltable_dir=\"test/CalTables\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "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+CjxoZWFkPiAgICAKICAgIDxtZXRhIGh0dHAtZXF1aXY9ImNvbnRlbnQtdHlwZSIgY29udGVudD0idGV4dC9odG1sOyBjaGFyc2V0PVVURi04IiAvPgogICAgPHNjcmlwdD5MX1BSRUZFUl9DQU5WQVM9ZmFsc2U7IExfTk9fVE9VQ0g9ZmFsc2U7IExfRElTQUJMRV8zRD1mYWxzZTs8L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS40LjAvZGlzdC9sZWFmbGV0LmpzIj48L3NjcmlwdD4KICAgIDxzY3JpcHQgc3JjPSJodHRwczovL2NvZGUuanF1ZXJ5LmNvbS9qcXVlcnktMS4xMi40Lm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvanMvYm9vdHN0cmFwLm1pbi5qcyI+PC9zY3JpcHQ+CiAgICA8c2NyaXB0IHNyYz0iaHR0cHM6Ly9jZG5qcy5jbG91ZGZsYXJlLmNvbS9hamF4L2xpYnMvTGVhZmxldC5hd2Vzb21lLW1hcmtlcnMvMi4wLjIvbGVhZmxldC5hd2Vzb21lLW1hcmtlcnMuanMiPjwvc2NyaXB0PgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2Nkbi5qc2RlbGl2ci5uZXQvbnBtL2xlYWZsZXRAMS40LjAvZGlzdC9sZWFmbGV0LmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL21heGNkbi5ib290c3RyYXBjZG4uY29tL2Jvb3RzdHJhcC8zLjIuMC9jc3MvYm9vdHN0cmFwLm1pbi5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9tYXhjZG4uYm9vdHN0cmFwY2RuLmNvbS9ib290c3RyYXAvMy4yLjAvY3NzL2Jvb3RzdHJhcC10aGVtZS5taW4uY3NzIi8+CiAgICA8bGluayByZWw9InN0eWxlc2hlZXQiIGhyZWY9Imh0dHBzOi8vbWF4Y2RuLmJvb3RzdHJhcGNkbi5jb20vZm9udC1hd2Vzb21lLzQuNi4zL2Nzcy9mb250LWF3ZXNvbWUubWluLmNzcyIvPgogICAgPGxpbmsgcmVsPSJzdHlsZXNoZWV0IiBocmVmPSJodHRwczovL2NkbmpzLmNsb3VkZmxhcmUuY29tL2FqYXgvbGlicy9MZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy8yLjAuMi9sZWFmbGV0LmF3ZXNvbWUtbWFya2Vycy5jc3MiLz4KICAgIDxsaW5rIHJlbD0ic3R5bGVzaGVldCIgaHJlZj0iaHR0cHM6Ly9yYXdjZG4uZ2l0aGFjay5jb20vcHl0aG9uLXZpc3VhbGl6YXRpb24vZm9saXVtL21hc3Rlci9mb2xpdW0vdGVtcGxhdGVzL2xlYWZsZXQuYXdlc29tZS5yb3RhdGUuY3NzIi8+CiAgICA8c3R5bGU+aHRtbCwgYm9keSB7d2lkdGg6IDEwMCU7aGVpZ2h0OiAxMDAlO21hcmdpbjogMDtwYWRkaW5nOiAwO308L3N0eWxlPgogICAgPHN0eWxlPiNtYXAge3Bvc2l0aW9uOmFic29sdXRlO3RvcDowO2JvdHRvbTowO3JpZ2h0OjA7bGVmdDowO308L3N0eWxlPgogICAgCiAgICA8bWV0YSBuYW1lPSJ2aWV3cG9ydCIgY29udGVudD0id2lkdGg9ZGV2aWNlLXdpZHRoLAogICAgICAgIGluaXRpYWwtc2NhbGU9MS4wLCBtYXhpbXVtLXNjYWxlPTEuMCwgdXNlci1zY2FsYWJsZT1ubyIgLz4KICAgIDxzdHlsZT4jbWFwX2QwNWFlY2IzZDZkZjQxMTM4NTgyNzA2NGQwOGRlYzJjIHsKICAgICAgICBwb3NpdGlvbjogcmVsYXRpdmU7CiAgICAgICAgd2lkdGg6IDEwMC4wJTsKICAgICAgICBoZWlnaHQ6IDEwMC4wJTsKICAgICAgICBsZWZ0OiAwLjAlOwogICAgICAgIHRvcDogMC4wJTsKICAgICAgICB9CiAgICA8L3N0eWxlPgogICAgPHN0eWxlPgogICAgICAgIC5sZWFmbGV0LWltYWdlLWxheWVyIHsKICAgICAgICBpbWFnZS1yZW5kZXJpbmc6IC13ZWJraXQtb3B0aW1pemUtY29udHJhc3Q7IC8qIG9sZCBhbmRyb2lkL3NhZmFyaSovCiAgICAgICAgaW1hZ2UtcmVuZGVyaW5nOiBjcmlzcC1lZGdlczsgLyogc2FmYXJpICovCiAgICAgICAgaW1hZ2UtcmVuZGVyaW5nOiBwaXhlbGF0ZWQ7IC8qIGNocm9tZSAqLwogICAgICAgIGltYWdlLXJlbmRlcmluZzogLW1vei1jcmlzcC1lZGdlczsgLyogZmlyZWZveCAqLwogICAgICAgIGltYWdlLXJlbmRlcmluZzogLW8tY3Jpc3AtZWRnZXM7IC8qIG9wZXJhICovCiAgICAgICAgLW1zLWludGVycG9sYXRpb24tbW9kZTogbmVhcmVzdC1uZWlnaGJvcjsgLyogaWUgKi8KICAgICAgICB9CiAgICAgICAgPC9zdHlsZT4KPC9oZWFkPgo8Ym9keT4gICAgCiAgICAKICAgIDxkaXYgY2xhc3M9ImZvbGl1bS1tYXAiIGlkPSJtYXBfZDA1YWVjYjNkNmRmNDExMzg1ODI3MDY0ZDA4ZGVjMmMiID48L2Rpdj4KPC9ib2R5Pgo8c2NyaXB0PiAgICAKICAgIAogICAgCiAgICAgICAgdmFyIGJvdW5kcyA9IG51bGw7CiAgICAKCiAgICB2YXIgbWFwX2QwNWFlY2IzZDZkZjQxMTM4NTgyNzA2NGQwOGRlYzJjID0gTC5tYXAoCiAgICAgICAgJ21hcF9kMDVhZWNiM2Q2ZGY0MTEzODU4MjcwNjRkMDhkZWMyYycsIHsKICAgICAgICBjZW50ZXI6IFs1MC45NzkzNDc4NzAxOTk2MDQsIDExLjcxMTI2NzQ0NjM1NzY4M10sCiAgICAgICAgem9vbTogMTksCiAgICAgICAgbWF4Qm91bmRzOiBib3VuZHMsCiAgICAgICAgbGF5ZXJzOiBbXSwKICAgICAgICB3b3JsZENvcHlKdW1wOiBmYWxzZSwKICAgICAgICBjcnM6IEwuQ1JTLkVQU0czODU3LAogICAgICAgIHpvb21Db250cm9sOiB0cnVlLAogICAgICAgIH0pOwoKCiAgICAKICAgIHZhciB0aWxlX2xheWVyXzA2MjA2Y2U0OGNjZDQzMzhiNTcwYzYxYmE3MjcwNDk4ID0gTC50aWxlTGF5ZXIoCiAgICAgICAgJ2h0dHA6Ly9zZXJ2aWNlcy5hcmNnaXNvbmxpbmUuY29tL2FyY2dpcy9yZXN0L3NlcnZpY2VzL1dvcmxkX0ltYWdlcnkvTWFwU2VydmVyL01hcFNlcnZlci90aWxlL3t6fS97eX0ve3h9JywKICAgICAgICB7CiAgICAgICAgImF0dHJpYnV0aW9uIjogIkVTUkkiLAogICAgICAgICJkZXRlY3RSZXRpbmEiOiBmYWxzZSwKICAgICAgICAibWF4TmF0aXZlWm9vbSI6IDE4LAogICAgICAgICJtYXhab29tIjogMTgsCiAgICAgICAgIm1pblpvb20iOiAwLAogICAgICAgICJub1dyYXAiOiBmYWxzZSwKICAgICAgICAib3BhY2l0eSI6IDEsCiAgICAgICAgInN1YmRvbWFpbnMiOiAiYWJjIiwKICAgICAgICAidG1zIjogZmFsc2UKfSkuYWRkVG8obWFwX2QwNWFlY2IzZDZkZjQxMTM4NTgyNzA2NGQwOGRlYzJjKTsKICAgIHZhciB0aWxlX2xheWVyX2NiYzA1ZGE4NzY2ZDRkMDdiYTZiMWY4NDAyYTkyNjU1ID0gTC50aWxlTGF5ZXIoCiAgICAgICAgJ2h0dHBzOi8ve3N9LnRpbGUub3BlbnN0cmVldG1hcC5vcmcve3p9L3t4fS97eX0ucG5nJywKICAgICAgICB7CiAgICAgICAgImF0dHJpYnV0aW9uIjogbnVsbCwKICAgICAgICAiZGV0ZWN0UmV0aW5hIjogZmFsc2UsCiAgICAgICAgIm1heE5hdGl2ZVpvb20iOiAxOCwKICAgICAgICAibWF4Wm9vbSI6IDE4LAogICAgICAgICJtaW5ab29tIjogMCwKICAgICAgICAibm9XcmFwIjogZmFsc2UsCiAgICAgICAgIm9wYWNpdHkiOiAxLAogICAgICAgICJzdWJkb21haW5zIjogImFiYyIsCiAgICAgICAgInRtcyI6IGZhbHNlCn0pLmFkZFRvKG1hcF9kMDVhZWNiM2Q2ZGY0MTEzODU4MjcwNjRkMDhkZWMyYyk7CiAgICAKICAgICAgICAgICAgICAgIHZhciBpbWFnZV9vdmVybGF5XzgxNDY5OGJlZDlhNDQwYTk5NWFiMzcxZDJkOTJhMDU5ID0gTC5pbWFnZU92ZXJsYXkoCiAgICAgICAgICAgICAgICAgICAgJ2RhdGE6aW1hZ2UvcG5nO2Jhc2U2NCxpVkJPUncwS0dnb0FBQUFOU1VoRVVnQUFBSllBQUFDV0NBWUFBQUE4QVhIaUFBQXFza2xFUVZSNDJ0V2RlYkJsVjFYR0QrVlVXbGhhVHFXRmltVnBZU2xsb1NDRGdJQk1Ja05FakFRU1FrS0hJY3FrekFVVmhTQ0lsRlJBUkNPakNBZ3lSTUFBaGlsQUlCQ1NFTHBKNkhTL3pweE9KNTEwcHA2NzM3dmU3OWhmMTVlVnRmWmUrNXh6Myt2OGNTcnA5OTQ5OTl4emZuZXRiNjlwZDA4OTVySWZQdnJwVi96SVU0Njc4a2VmZlB6VlAvWW5KMXh6OXlldHUvYkhuL0RzNjM3aThjL2Q5cE9QZmY3Mm4zN1VpM2Y4M0IvKzljMC8vd2V2dU8wZUQzNzFybDkrMENtNzcvbkF2OTM3Sy9kLzNiNWZ2ZC9yRC95YUhyLzdob08vamtOL2hyL0QzK04xRDMzVnpsOTgrTXR1L1FXYzg5RXZ2T2xuSDNmeURUK0ZBKzlsRDF5RFBmZzd2QWF2eDNsNGZUZ3ZyaEh2d1d2RXRmemVtdy84eG9QZnR1ZmVqemo5OXQ5NS9BZHVldURUL3VlYVI3N3c2MXVlL09hTGxrNzgySmJOTHpuditnMXZ2SHJuK2Y5eTg3NHJQcnozNEcxbnJpd2YrTXJLN3UzbnJkeDZ4VVVyTjI5WjN4KzNYZjJkbFgyM2YrUEE4dDZ6OEhlWDMzYkIyOCs2K3Z1dmVmdjNscDZEY3gxOXhuVVBmK1E3Yjd2dkEvNXgvMi9pOCtJKzRWcHdYZnljdUg3Y1c5eG4zRy9jZDl6LzBvRy80Zk9KRHUrNTRmM3c3SGpZKzh6N2lkZmd3T3Q1OEdmZStYalBlYjk1ejNIdzMvaGRweDlnYXJnVU1zS2xOeHpuNUlmMndDckJWZnVnZUIrOEg5Nzd2bTg2ZUs4SG5iYjN0eDcyamwzM2VkejdkendBRUR6dmExdWU5SVlMbDQ3LzBPYk5MLzc2dHUrZFNyRDJMZS84N01yQmZWOWIyYm50L0I2b216WmV2TEpqODRZZXJMMjNuRXV3TnQxeTRWcy9lOVhHVjc5bC9kSTZuT3ZKSDczK29UZy8zZ3Z2NjMyQmNPMjR0N2pIQkdZcXNFcFFXYkFzVkFwV0ZpcDd2eFVxL0UxM3pORmJmbkJxdUN4WUpiaHdzU1c0TWxiTHdxVldDKytKOTRjbGVlamJkLy8ybzk5OTYvMmU4dkZ0RDF2MzVjc2UvOXJ6bDQ3OXdLYk5ML3JhZFJlLzdzcmJMM2pIanIyWGZlZ3dXTGR2dmFBSDZzYUxMK2tQV0s5RFlPSHZOdDc4bmRNK2ZjV2xyM3pUZDVaT3dMbWUrS0ViZngvbngzdFp5NHpQaCt2bEExeVV0Um9MbFlJMUJpcjhiUStXQnhjT2hZdHZnQmVxVlJnQ2wzZmpNMkR4V2pJZkd1ZFhxMFYzQ0hkMTFFZTJQK1RFTDE3K3VGZWZ0M1RNZXpZdVBmOUwxMTd5TjNCdGg4SGF2L3VjM2tJQnFPc3Z1blJsKy9xTnZmWGFzK05iQU91bXZVc2Z2SGpIUlcvNXhHV2JYbmJxQlV2UGVPYm5yL2dqV0VKWVJieVg1d1lKMWlMY1lNWmE2ZjJjQ2lvOXJBYzZETlphd2RWaXRld0hMN2xFYTdYb0RxR3pZRjJlOGI5WFBlYmwzMXc2K3ZSTGxrNyt3alhmUHdXdURjRDBZTTIxVkcraEFOUzJDemIzY01GNnpjSGF2N3o3Yy9pN0RUZGQ5SS9RWjZkOGUrbHB4MzdtNmtmREVnTGVyQnZNUW1YQjBuL1RXZzExZ1MxUVdibFJnZ3F2dndOWVU4R2xnQkVvUEZnY0NsZjB6YzZBMVdxMTFCM0N1a0RBdi9nYlMzOEs4UTJ0Qk5kMkI3QmdvUURWMW5Ndlc3bnV2S1ZlYTgwRlBjQzZjYyttLzdodysvcC9nRDU3NWJlV252clVUMjE5QklEdDRaM2ZCNnNocDNLRDltZmVjNWtLcXBLR3RWQ3BZZUI3ZEU5L3l1WWZLTUhGRDFLQ2l3L1FnOHVDcFhCUmI5VmNvbldIR2JnaUVROTNDT3NDQWYrWDUydzVDdUliV2dtdURjRDBZTTIxVkE4U2dMcjI2NWYzY01FdEhnSnIrNTZONzhkS0VtNzBKZWN1L1JrMEc0Q2xHOVFiUHNRTmV0Yksvc3hhS3d0VjFnVzJlQUVMVkFRVnp0dURsWVZML2JrdTk0ZkFoWU53cVV2TVdxMnMyYllpSHU0UU9ndXJ1T2Q4ZGNzVElMNmhsYjU3NDNmZkRHQjZzT1l1cndjSlVGMzE1YXRXcnY3cWxiMWJQQVRXRFhzdWZ0KzUyemI4SGR3bzRJUnJ4WG54SHZ3cytrVVo2d1k5MER6QnpvV01aNjJtaG9wdW5oWlpvY0xuUEF4V0NiQ3A0QUpNMENFOGFMblVmZFNzbGdXcjFXckJIU0lzQUJnZzRDRytQN0swK2EvZzJnRE1ZYkFBRXFDNjRxeHIrLzlDWiszY2RqNSt2MjMzZDkrTmxTVGNLT0NFYThYbjBjOWhiM3JXRFNwRTN1RTlBd3VWM3NQSUJiYmN3MWFvWExEV0FpN3FyYXhMckZrdGUyT3NpS2ZPZ29ESHloQWhCN2kyNjNkdmVHOGZJRVZ3RlBycXlpOWV2YngwNXJZZUxyakZPVmo0L1hXN3Z2TXVyQ1FSWU1XS0VLNjFkN1B5R2Fad2d4RlkxZ1h5TTl0RjBGcEJGWUsxS0xnVUxGZ09ISVFyNHhLbnNGcDBoNDk5M3kzM2g0REh5aEJhQ1VGU1dLTERZTTFCV3I3c2M5Y3RiL3prOWg0dTZLeERZRzNkOWUzVElmZ1JZTVdLRU1JZFh4eHZsYXZSOWxZM1dMSldyUzV3U3FqMC9EeXZRbFVFeXdQTUM2UmF1T3lLeklZaWFuQUJnc2hxZVdCRk1aeklhdEVkUW1kQmRDTWQ4NDZMbDU1Nzl0YUxYd3RMMUlPRnFQc2NwT1ZObjc1aCtYc2Yzd0c0ZXIwMS8vbWVnN2Q4R2xGNkNINkVHckFJd0lJQTErNVpxNm5kSUsyVmhjcXpWa09qNmhtb1BDdkZCY1VkeFBzaTRJcmlYTXpmNFFIRGVoQXU2aTI4ZGlxclpjTWlkSWZRV1FpVUloM3oxZzFMSnlHV0JVc0VjUHFvKzF5dzkxQmQ5T0ZibHk4NTQ2WmV3TS9CMm4xZ3g2Y1FUSVhnUjZnQml3QmNQODd0eGVSYXdnd3QxcW9WcW13QWRDeFVQSDhLckJwY2RvV1NnWXRXQzJBUkx2eU1laXNTOGhtcjVkMHdCa3pwRG1GbGtKQ0dnTWZLOE13ck43NEthUjJBMDRNMUYrdzlWT2U5Zi9meSt2KzZwUmZ3YzdCMkhianhrd2ltSW9hRk9Cak8wZXUyK2JtdHBXMXhnNjNXU3FHeW1rZDF6NWlvdXJlNnpVRFZCMGl6WUkyRlM0T1ZxcmRvdGZDZ0NaZDFpYldZVm1TMTlNYnhPdFFkUW1kQndDTm4rTitYYjNyNTBxMFgvaFBBNmNHYUMvYmw4eit3Yy9rYjc5NjNmT0dIYm9lUUIxZzdEOXh3Qm1KZTBHV3dkaER1T0IrdHJBVXJHMmJJV2lzTFJJc0x6RUxGZncrRnFrL3B0SUJsQVpzQ0xsb3RnRVc0bU1qMVhHSUVWdFpxOGYzeHZ0QlppSnJEcFNIa2dPZzd3RUdlc0JmdTMzenZudVd6MzdFTXE5V3ZET2RnM2JaLzY4Y1E4NEl1ZzdXRFM4WDVQQXViZFlNMWErVzVRRUpsOWR5UkFCVmVOd2dzejNxVjRGS0xZWlBEdEZxQUNtRUF3c1VvTmg5WTFoMlc0RklSRDNqeGZ0QkljR24vZnVuU0M1RC91MzMvdGsvMFlFRzRmL1gwZzh0bm5UWmJQdWVkK3dFYXdMcDEzelVmUldnQ0VYdXNDUHN2eFB3NmJXMlpyYjBhWTYyc0M2UjI5SmI5WTZIS2hoUDRIbEc0WnpCWVU4SkZxNFVIVGJqd00xZzF1c1JXcXhYZFNCWHhlRThFU3VIU0VFVS8vNGIxZncrTDFJT0ZNTU1YM2paYlB2Tk5zOTVxSWVTd2UvdDVxTVZDY0JTaEJxd0ljWjA0cDVkQ3liakJWbXRsb1ZLTE1oUXF0VkplWURkcnBUU0dPQXFzSVhDcGtOYllFc3RhNEZwdzRQK1oxS1ZMeklLVnNWcDBoOUJaY0dsWUdTSk5jOHUrcXo3U2c0VVY0ZWZlTWxzKzQ5UlpEOWpjZ2dFc0pLcXhna1NvQVZEaTJ2R1p2THljZFlNcUgxb0Z1MExsaFRROHNWNkRTbzhwb05Mckd3MldwN3N5Y0dtc1MvTjRzRmlFQy8vdVYzREdKZHJDdndnc0R5NjFXamczZ3B0d2FiQkFpR1doSmd2aUhTdkI1VSsvY1hid1AwK1pMWC9temJNK2xyVm54N2VRVDhRS0VvRlZRQWxBQ2IyOXByRmdxV0N2dWNBcG9NcUdFeUlydFJDd2hzSmw0MHUwV25qZytDOXJ5UEYzOWx1YXNWcmV6U1hVZmRKNERqSmNHaXdRb3Vtb2NPakJtcThFRDM3c3RiTUQ3M3ZWRElEMXNheTl0NXlMdEE5aVdFZys0eHB4N1RpM0IzdEpYMFV3Mlh0bnEzZHRqZGRxUWxXelVnc0RLd3VYM2lndkRJQ0hqWWVHQS8vUEFqcjhyV2UxSXJBaXVHaTFjRTY4SDF3YXlsOFFjZ0E0RU9rSU5SejQ0R3RtKy83NUpUTUFCdGVJT3EydHV5NzROOFN3NEQ3eEpjQTE4WUh3dmZYaHRGZ3IvUjFlcXk1UXkxUXNCSjZ1bWdvcXIwUTlzbEozQU91NEoyKzZtd2NJZnM1alVYQjVNU2E2UkZncmhBTnc0QUhTSmFxUWI0RkxYU0xmRStlRGxZUkxnNEFITUVqclFFc2h4TEQvM2ErWTdYM0xDMmNBckErUzd0OTlEb0tvNy9yKzBsL0FmUUo0bk10Ny94SllOYWh3djNpditFV3doWVBXV3JWQXhZZGZDeWQ0TlhpMk5NbWU4dzVnUlFEcDcxb0JLOEVWZlhpS2FycEVRSVVBSkN3WDY1M29FcWV5V3JDU09EK3FGQUFNd0lHV1FtQVUxbXJQNjU4M0EyQTlXTXNIdm9Lb080USs4b3g0TGM3SGgyTEJpdlRWV0d0Vmdzckd1a29pM1NiSk0zcXFCR29JbGdXbzlMdEZ3Z1dBWUExZ3RRQVc0ZXFqOUhPcnBxdXdtb2lQNEtMVnd2a0FNUUtsQUFiZ3dPVWhkZ1ZydGV0Vko4NzIvK3RMWjBqdnJNeFd2b1RnS0lRK1VqbDRMYzVwSHpDLytWbDk1VmtyWEtjK1JGcnFtcTZxV1pYc3lpOXp2Z2lxSHF6amo3cTBQenlJU3NjWXVLTE9FbG9Td05VbnBBODFQOEJWNFdEZkhuN3ZDZm1TMWJJdVVkOExFQ05RaXFJL2dJUDJMNEFGYTdYelJjZk85cDcyb2g2c2d5djdQNC9nS0dxNFlFMWg4WEJlRFJpMnVrRUxIQVc3WmkyODJ2bEZRZFdxMGFwZ0hXbHdNUnpBR25VVTV1Ry8rRGR6Y3pVaG4ybDJwVHZFK1pIYVFTeHJ0ckp5TnFMdXUxK3pibmJMYzQrZDdYbkR5VE9zRWxHV2pKQUVTbTJ3cU1CMTZzUFJCOTRLbGxvclhKc0dMMjJuVDdRQ3pFS1YxVk9SbFNwQmRUanlybUFOQVd4cXVOU2F3Tld3eEFVV0N3K2ZiZXo0blNma2EyQlp1UGcrc0l3SUlTRDRDY3NFc0hhOTlQalpEY2NmUDl0OXlrazlXQ2lwUVF3TERhcDljSFgrZXJvOEMxWkdYMFhXaW02YTFtb29WRjR3dFFTVkZ4TEtBR1VyTFVLd3BnQU1ON01WTG51amNHTVpFZ0JRMERXQWk2MVdkSW0xVHA1U2ZBdXY3L09IYzdjTEFZK1FBK3JhQWRiT0Z6eHRkdTJmbnpqYjlZb1RlckNRb0VZdklTcFBZZVZ3RHR2anAyQzFhQ3NLZG5XQkdoUzJZajBEVlZRYmx0RlRHYWc4b0E2RFJUQVdCWmNGTElMTDNqRHFETGlxdmtacjdnSmh0UkJ6QWx4d1JmaTU1eEk5cTFVcXI4RTU0SFlSS0VYOU8wcG5vTEhnQmk5LzRrbTl6Z0pZU0VDalhBWjZESzRhNTdPZndRcjNqTFhpb2diWG8wbGhmbWs4WGRVQ1ZSUk84RnlmQjVTRnFnUVUzN2ZUMWQ4VTFtc3F1UGp3Q1JjZVBGd2lyQmJnQW1SOWtXREZKV1pxdC9CNldFYWNHK1V3QUFqaGhodWY5WXpacHNjOGQzYnJ5ZjhQRnZLRVdEbml2ZkhnY2IwUldGazNTR3ZGejVweGdVT2d5cm8rcjN5bUJTZ2VIVC84SWdFYkFoY1AzbkJZQ09nYUNIaUFoUVB1RVM1UlkxdTJoQ1NUOXNHTmdEdUVpMFVvQVdrZDFHSkJYMjE0eU1tem05WWQxNjhLRVpYSHlyRVBnOHpmVDF2ZWEyQ1ZVamVlWUdlYUtJTEtFOWVabForZXgxcXBDS29zVE5wNjF1a055QUkyQmk0RnJBU1h2Wm40Z0xCT3NGSndoYWhYaDRXaFMvU1MxRmJBUjNEUkJlSGNTQzRqK282VXpuVlBQMkgycmZzOFozYjljYy9zQTZTb2lVZnFCNERqTmJ4M0piQksra3F0bFFwMmRZR1pGVnRVOFJEcEtTK3RabXV5U2tCRk1PblJlU3NZRFpDT3RWNUQ0YktBOGVZRElqYWNRdXR3eW91WDdxbUJwWm9PTnd5Q0hBMm9LRkVHU0JEdVo5OXJYUy9na1N0RThCUy94elhnZFY3THV3cjNXclNkMW9vdWtNMnUrQXllV0M5QlZSUHBVY05MMWtwNU1GbWd0Qml6SzVWNExNSjYyWDdGR2x3S0dTNllRcHZqaUhCWWwraEY1YjBtQTV0ZkE1eFk4YUVyR21VeTF4ejlyTmxuZittRTJWVlBlVlpmTm9PZjkxV2o4Ny9ENjIxR29RVXNLOWc1WjBwZFlHUmxiS1ZEU2FSSHJpK0NxaFdtcU5tbDg0ckZGZ1ZZQkZkVUxHZ0JJMXg0Q0d5R2dOWGkwTFBJSlVaZ1diaHdrM0V1QkVGUjR3Nmd6dmlaWi9RclExU1FvbklVSU9NOWNEMGVXQmw5WmNNTEt0aVpJcG9LS205d1hzMzF0UUxsM2Q4dTZ2SW9BWloxajk3UGg4QmxCNDNSdXJDRml5NlIvWDFSQlVTdDN3NzNBZ0llc1N3MHJBS3NEOS85Mk5uU0h6OTdCdEJRcjRYRkErNFhyUkloc1dEVnJKVUtkbldCRVZTNmNzdlVwZHM2K1ZyUFlGWTdsWG9QOUI1MFVUQXRHc1NhQVd4S3VMeFJQcmdtWEN0Q0JBQkJYU0xjWk9RU280cEwxWEdBRTdFcVRKdUJLM3ovRHoydER6bWdyeENkUEQzTTg3K3pBbHh2YXNrTjhzdEJ6YWd1VU1NQ1hsZE9xWVpLZ2N5NHZoSlFGcVlXb0E2MzJFY21zUVpZaTN1TWZ0NEtsejRnZkJCY0swQ0NGVkdYeUtKQXpiRmx3TUx2NEU0eDdBTlZwRnVQT1dIMm5yc2RNN3YwVVhPdzVoWU1aVFhRY2d5TVpvUzcvWElRR2dabWVaMU1aa2ZORTVuMGpBZFZiYkJISzB3MW9QZzV1MmlnVm1SdUZiQ3NleXhCMXdLWFhhcmpRK0ZhKytYL29SQUVkSmUydk92MWU4ME85a0E4Qzhsb0ZQVWh6SEFZck8zck53STRRTXZBcUUxTDhlY2xsMjd6Z2JoR3E2dHFVRm5YcDE4TXI1OVROWndkZFZTQ0tRdVVKMW02MGdqQVdpUlh2NTFqQU12QXBZRHB2M0U5ZlZQR0liMEZxd1gzMlBjTkhtcDdWK3RybTBsdFdBT2ZHU0dGZnRib3V1TjZzSHBYZVBPVzlhaUx4emxwbWV3UU5LdXZMSHhXc0JQOERGU1JTQys1UGtMbFBjZFdtTEpBSFFhck5tTXl1ckFwQVd1Rnl3WWVjVE1RUEdYS0I5WUxvSGxOR0RXcmhiOUZ5QUdOcTBoQ1EyUDE0bjNudHZNUkhHVXFKd0xMdWtHRmlnVjh2TmQ0TDlWVlVadVhqYVJiOTIxZFgyU2xTakRWckZNV3FNTWFLenZFTkFPWTFWOUQ0Q29CWnFzRitITjhJRndQUU9MdUV4RHluTEtuaVdvOXZKZ1cvZzR1RldrZDFHTmhWWGpGVWV0bWNJMm9pOGU1OUF1bFZnUC8xbXUxK2t1ckZ4aGF5RUlWRmZoRnJzKzJ5YmZBRkFFVlFXVmpqajFZUTJkNlIvMzk2ZzRXWWIwOHVQQmZ2Qyt1bmNGVFZweHFYMkxKYXZGQjlmcGovbHJNd0VKcDhwbjNlT1lNSW41NTVlQVgwWm1EejY4enJ6eDlaVmVMQ2dIdkxVTUxHbnJJclB3aWdhN2JqbmlyL05LbURDVjNGd0hsd2FUWDJKVmFmR3FBbGVZblpRQ0xvRk80MURKRmdQSGZLdWJaNGFORFoxVnYyWWVtQ1c5WVBjeHhRSGZPT2ZjK3FVOUNZeUFiSnRQZzRXbThTY0hTYTdIMS9WcHV6TkNDUnQremZZTXRWcXEwdTBmR09sbWdTaURabzh0czlKUGRtTWNEVFBWWEtmUlFneXNER0VER2V6SXliL3NTVlcveDBKdk1nQ1VzSGFMc2FMRy81T0hQNnd2OW9MbVl6dEdZazlWWGFxMjhMV01ZV3RCZ3IwSmxKeTNYWEovVlVxVjBWZ3RRMGVGQjVPVmdPMi9VYzYybnJBVXdkUnVlL29xQ3FTMXc2ZS94UHJpeHNGQnNmbVdyUGp0OHJKRGxEV2MxS2NCQ0NUSXFISkNBUmpNRld1OWhzV0Q5OFBrVUxINCthNjI4TDZqcUtnMDcxQm9tb3RpVWJXSXRnZVF0Vm9iQTVJRjBKN0JLT3lKRVZxd0VtSVdNdmw0QnMrN1I1aFF0WUMxdzBWb3c3Y05KTmh4RnFTN1JqaDFpd1I5Y0tLTHNHTGJXVjQrZWNlb010VmpRV0FpUTRuTlppMkxqYTE3TG1lb3F2bCtwQkxtVWxva2FXS09qRmFoV2tHd2xTWmZkZGlQVHhHZ0I4M3cvSDRqVlh6VzQ4RE5QVzVXbTNlRDkyTzFzQitqeVc4NWdJYTBIL2diaGlqNnRzMlB6QmpTdG91Z1BZaDZyUWtES3RucWJ4ckZqaDJ5OGlOSkFTMlcwcU0rdVVxTzBqS2ZERmcxVDVGcWozczJ1RkR2SzdGQlJHdlpoRDl0NHlRZml0ZlJIdWNZSUxyVnFlbzJzUHVWUVhUdmpsRjhPNml1NFQzUTVZMTRXT3ArWHYvVDJQaytJV2l3TWFJT2JwSURYWkQzZlY0SHdvR0tlazVCb0o3V0Z5ck5TSGxCZUYzWUdxTWpGWlNDcVRmV3BnbFdibUJ4cHNjaE42cEpZYjdnRkxDbzJqT0RTR0pnVno5Uk9IRkhKbWZLMG9qcGRHWUZWYUNtQTFZK05QRFI3Rkt0RU5xcXFnTmNxRUZ1Um9HNldEeE0vczFCNWFSbjlra1liUHBXT1ZxM1VDcEtGeVRiTFZyYzhhWVdzWnNYd1FIRkVnSG0xOFI1ZytGa05MbDRQNGRJSnp0eEY0azdUQnVkaUh5NFBDV2RzS0lEUlJXZ0hRNUVmaHVBaTRRMHdLZUExbm1XaHN0a0pGaW9TcWxwYVpncWdNbFlwMnQzV08wb3dhVmRQTjhWVW1WS3ZvR2ZGRkREdTdjY1ZrOVZmdHRxMEJTNjZKbzBqNmZ0ekRDVU9Eb0dEcTBQcUJtMWdzRkxvMkVFVUhxMzFCQXR1Rlo5RmMzWTYwZGpiU2c3L1Q0QXRWRmFnUjBDVmx2d3RNTFdBRk1GVW11ZlFnMVdxQUYzazRMVVNZS3EvcklXMDhUREM1VFdGNEwrMEpKRVZwVzdoTUJMRXZqaEZHZnZsWUpJeTlCVzJSSUVyQkZqUVlVenQ4S0hyamFjRjFnZE5xTFY3Mm82QnRLR2EwbEJaQzFUTnhWbElhb25vR2t5MUhzUkJqUktaYVRSZUNLQ1VqcURGNEg0MEdwMjJLMVhxc1JwY0dwRzMrL29wWkFRTHEwZUFBN0RROFl6SmZZaG5JVmdLeUNEZThYdk9iZERJdHRZMnFhYmk2cFNmaVZCRkEydnROblNaS0hmTktrVXdsWTRTVUxXNURkVVcrekVEUTd4Nkt3dVpkUUVLbUMzWHRiMTRCRTVqWVY3MWc5VmJmRC9lZE41QVJ1dVJ3SVlyeEJBMmhCd2c0bEdxRE1EUWdvODhJZ1M4Z3NCejhwdXZibzVDSFA4bFZIYTB0Z2VVVFdCSCtxcG1sV3JGZkRXZ1BMQ2lEaDYxZmwxci9kUVkwTlNhMmJBQUg3SUNSaDJrN3RHV291RGZ0RjRLVjlURHA5WkZId1pERGFoczRNWk42SHBHY1I5bXdPTkFHSUwxWGh5QnBPRUsyd3RJWGFjekh1eTRSUnM4amhwM00xdnRldHYzbHVxdlNxK3I5UlBXbWlxS0taV3hvTlVxUnozSWRMbk5JYlM2U1NiZG8xY0h4WFBUU3RsNmVYMFB2VG1NZHlHSXlybnYwRk1vN0FOWTJQQVNIZEFJbkdLK1F6L2QrZEFNTDE0cmI2Z05HN0Qva0tFR1czVVFBZVhWUlVXN2I5UzZhbXFpdkhTZTFuTGx3eW1kbXBVWkE5dVFFbVhOdFdtUWtPRUp0VjRxZ0FrY3RSL2hpa284VkhQeHBzRmlRVVBCNWFHS0ZJQmhCd29jaUcwQkxDU2lPU01lQVZkZUk2OURReHlzaldkUU5Dcnp0aE4zTWh0K0RvVXBxNzJ5ZFZ0aHJuQ0lTMXN0eUtKNW5OUXNkc2xPQzhDVm93Y1hQN2kramxxT1E5NWd0V0NaVUVrS1M0VjVEbXlReFVoSkFBN0lHSUhYUmcxYmwwYnh6dXZWcExYS2dLZ0RwZ1pVeVIxbGdweTFXcTBzU0pPQXRScXdXY2pzM0hOZGtxc1Y0cytvdmJRUVVLc0sxQTFSUTJCVkNIaFlQdys0c0JKa2poQlE5ZFVTcCt5K0o4QmliSXF1VkJQU0xPR3g2UnpyS3FPZHUwcEFlU0MxSm9scnh4Q1k5RDI2c1ZQN1ZnTkdxNTNVbFZtZHdvZkptaWRiWmNxRExwWG5vWTREV0hCekFBZ3VFUlVObkRrUDhkN25IbUhaRG9Vb0dOVlhLMFZMcGU1T0xaTU5wRVpBUlFJNXlnbG1hNlZhQUl2R1F0V09MZ1BIME0wRXBvWk5yMGUxbUhZZ0szaTJGb3lXajB0K0hmREdGTS9oYmU3bS80YUZBbWpRVXVoYlpJNlI1YzRFU2t1VGJkV2xSdVM5ZEk5Tkl0czJ0ZEtXYmkwVkNrTmdHbk5NdmpQRm9vQXJoUzVzYUVGalYvck5wbnVrN21MQUVwRGdRUUlxd01VOUU3bDVGUGYzWWZrTnB3eHFQRXhYcXdxVTdSRHk2c285RUtLdW1CSk1ZMU0zazRMbEpaaW4zQkZzekxsYnhpSkZVMnhzbTVZMlBQQ2hNRzZHbTh1SnpOeDhpV1BCdWZNcmY4OWNuc0tnOFN2YmVlME42YStsYVZxcUVtb2hpT3lLYnpLd1NoV1pZNkJvS2NmSm5qdUtnVVdKYUhXYkdyVW5jQm80WlNpRExoRVdpZE5yZEI4ZlRhWlRnTnVpUDlXQVZnaDc0Njh6dTVoRmVUd2JIYTlGeEJjRmtndFdxU3B6Q3RCYUlhdWQyMHNaMlZTUFYzNURrUGd6aFV1MUQwdHJtSmpXcWdUbU1yazZ0WUZSdFJqZXpBaGQxVnFMRWtXOGE5VUV0Zyt4MXFYVEt1YUh1Tkp1eUVPZnduV09mUThWKzNiVmFHZEw2SG4wZDFhUHdaTGdBV0hWUjQzRWxCSTNSMmY2aDd0U1VNQnpoY2ZWcGxvazdWbTBBQ2cwMFpUaTJzYVZObVJSbTJBNEZKd1d1THJWQWlxNzJod0NtTmQ4WWN0czlEV0VTeFBGdU9rTUgrQjEydExPTGVqd045dzFnNER4ZFp5UnluaVcxeFRSRXZpMDFpdXlSSm1ZVlV2RXZYYWt3VnJyMWQ5VTlXRGUrZXpFWm11MU9PZUswWEVHVnBuL0kwZ0tDZitPVm9OaVhlTlZtaWVrKzdPbE1OazRVNm1XUFZNcVUwci9aUE9EZzhDYUt2Sys2RWg5WnBhOE4vblpxMjVsQVNCK2pnZkQvK0puWEEweTNJQUhTR3RGVjRuZnd6cXhPVlp6bDVyOFZ2MWxLeGFHZE0rVVZvSVJSTFh4ajltU21sYVh1R1pnVFFWYVZGam9WWlRhdG4wV0FmTEdBUkw4bStNbThSQVJYZ0FrZElHTWUya2poRDVvT3lCRnF6QXlNeEphNGxOalhOcFljSTVJc0taeXBWRlp0RGNKa0ZwTFJUNGgwTlNLYmd6ZzdYYktua0xxTWhzZHo4SXp0a1ZyaW56Z1ZBTGVBM1ZOSXU5VGhTUmFoTDdHck9qNjlNYndvZEFDYWR5SmNQRjNmWFdEaWFacmxOeFdmVnBRbzFuMjJUYXVJZk1WU2xCT0Jkd2R3RnFOME1KWXVFcmQycDY3OHdLcUdsYlFBanlLZWMwMXdvSkJRNmsxZ3l1MHRWTjg4TFkrekJibEthRGVPQ0J2aDlic2JLck1lS0dXWVI1RE9wOWRWN2dJUzdISUl4b0RFQVY1VmJUcmNGNDJaUEQvT1lhSUR3S3ZaZXlKQ1djN0hWa3JVRmxOb1dFTS9SWkhlY0sxT0laV1JiUkEyTldzd1pGK2VEUFV2YjF0K0crQ3Bla2U3ZjNqS2hGZ3FRYlR5VFNxdlFpZlFxTWdhMmRTYmJCdTVxR3VsdHNjY29SampLWSt4cGp1MG1IbnZudnZxMzlyVjRnTUJRQXN4cTQwSHFVVkVEcHVpTUZURFNkb1hsRGh0VldyWHBWcEJxcE1BZDZVY0V4MWRGT0RNcld2YngyeWFsMmo1eGJ0ckZER2VYU0doRGFhNnNaSjNzWUJUUDlvNjVjTm1HckxWN1JWU1FheUlaWnN6Y0JxQWFiVjRyVDY4akZtM3JOT0ZQSWFlaUE4REpSaWxjZmZNVVpHOThYclloeUw4SG1mazBGS3I2WktXL0h0ekt5U2k4d0k2cWtCRzFxV2ZBZFhPQ1ZFUTMxM3lkM1YzSjRYREkxbVNoQW9oUW9QRnJrKy9GNWpaTGFEaGlVMWdJR2dVckRyNTlBVWpqWldXUDFsQjZoa0FjczJRNHlGYll4dzc4WDdsTnBuakY2S0FLcVY5V1RtU09ncWtPQm8xQjM2U1FPb1RGVGJCOG9oSXRSUk9vWmJPNWh0eU1GdVZHQzNQeWtCRmtIV01vOGhPelp5eXBCRXQxWUExU0RLQWxScjdWZW9tQTVpQXByNkNkZkgxNm1HMHBZdGR1UndGY2pQWkR1RkxGeWFrRllCYjF2cG8wYUtGaXVXYWEwZkExc0xaTjFxUVdSQmlpQnFUZlZFZzNEMWZiUnVpd0ZRN2F4bTJFRlhqWXpPTTJlb0hkcDRuVzN4MHB1dm93QnNVNjBObkhxekdpd3NYcWdpeWlsR1UvWmEydVBINnF2QllHWDBXTVlhdGRSc2xjcHR2RkpscjhxQnVraUhkUER2VktkcFphbXU1UEJ6VnBQeU0xcHI0cms4WFVGR2N4bjB2bnFseXkyZE81RW04eXphMEJiNlNjVjdCcUNzTldwSk1IdFFsWFlQODFJNnRDNnMxbVFnbEwrenU2THFBMlpLaHk2TDllNmFvTllINmdWS2F3UDRvOTBmU3JOQnZkVm5xVUtpdFRvaXU1dkZLTEJLcTdHTVN4dGE5VkFxQ3JRaTNYdHZEWWl5UG9tQlRsc3ZyK2ZSQkxMV05mSG4rdjkyTEtSdFZvMFdSalkvYUd1NnZGbXZXY0NpOXh0YWh0UGFpdCsxTE9kTEs3U2h2WU5EdEpRbnpxTUdDdHVGVE91aVZhWGFhTXJQb3dsbWFpU3NDaGxJVlcybEQ4Yk9mN2NQdVRSSDM2dUdpRGJjckFHVzZVdk16aVROMUczZFNid1BnV2hLa0dxdU1nTlVGaXJOOVZtaGIzZVhzR0tjKy9UQURlb0FOWDA0WHNJNWs0cXkwMmE4SFcxcm16dTA3b0RxUVRiRjhGdEMxdFZLVVVxcnRiRWdsWFJSYmV1VDBrSkFDL25zdzdidFlXcWR0T3JUOWdEUzVla2NVUy9SYk4vUG5zTmFMOTJxeGJOZWVsMDIrVjRETEF0WnRuSzFCYlN1MUU4NEJVelorcTdhNXVRdFVLbW04a1l3cXI3aWplWERwUVdKd09Id2Y1eFhvZkJTTlJyanNvVi9IbHgySUxEbkhoY0IyQlJ6SUlyMVdGT1VDR2VyUWJOQVJWSDBxSjdkUW1XRGZ1b0d0YUpVVy9EMW02d0JUTUxDRm51ZWk3bERIVTJrMDJRMDNlUEI1Vmt2L1h0UDNHY0F5N2pKc1pDRjlWaXIwUXVZQmFvRkttK3V1ejVjYndkVlR2elRwRFRQcFJOanZQSWJyYjFpTzc0dUFxakxkR0F1VjVMYWdwK0Z5MW92QzVoWHdURUVzQ0hGaHhuSXVpbkxpR3NyeHRvbVRLWFpERjdDV2V2WlMvTVM2STQwR0dwSEcrbndYQVhMM2tDZDE4NnEwMGpYZVhCNUxzNWFTRHNodWhhYThFYVFSeXZKYk4vaVdNaTZSY0dVQmFvR2xlMjJzY0tiVU5sVmxkNDgzUXhLdFpTZEEyOHRnYmU4MXprT3JNL1NlNkZXMDRNcld2MTV4WWtseUcwNG93Ull4b29OYVU4clZiTjBVelEzREFYS2c4cCs0NnliOE9KUDl0dWlPVGFGeXM1dTBKQ0V0d3U5amZld2twUzdXZGlwZ2Rxd29VTFdqdXkydWI5b2hWWUtsaXBnSmJocW0ydFp5RElGbGhuSXVxbGhtaElvNzRiWmFMbTkwYlFTZEQwV0ttOThKRytpL2s0SGhhaDQxOGt1Q3E3ZHZkN202UlF1RzNBc3daV0p4bWVzMTJvRDF0WEVkMjNFVVMxSlBHVFd1NzJoWGllenphZHA2UzhqNGxyWVo1dFh2ZFdYdFdSMnRpbkxiTlQ2NkRBMzc3eTJIVDRhVk9zbG9qVnhYYk5lWHJTL0ZiQk1nVisyc3JjckFaWHBMUndDVkdTbExGVGV5cyt6YWpvaWlOcEhvYks2eW5PajNyUS9MWDlSNFc1MzV2TDIrNGtxRkR6dFpjTVNVWTZ2QmJBYVhHTUJxMW14Ym1pemFzdjhoVmFndkpWZkpFYnQzQ2xOMmRpWldlcm1WQXpiQWJqcUJuVUtNdDZEVmFRMjhleUZDcUtncGJmMWNSWXVEekN2NVcydEFPTjdka01IbmswSmxMMEJYalRkZ2tXbzFEMUZVUEd3TzEvWWljczJ2YUoxV1FvV1Y0V2FjRmFyNVgwK1Q4Y1FMdXEyQ0M0UEdpLzNXQUp4S0dEZTlXY0E2NGJDTkRRbXBTa1QrNkdqWkxMOTBCWXFQQlJ1aEdSbmo5ckV1cDJtWEhJdmRtZDZ6bmJYamFNOE9LT0g3QVVUZFJjdzYySWp1S0lDd1Zick5RU3dMR1RkRkJVTG1aMitQRE5iY24ybGZRKzFYdDJyQnJWanVsWC8yRzFQb21vREc2amswRFh1U3NiNmQ4OGk2SGxxMy9xeGNMVUNsclZlWXdIckZqR0hQZVAyU2xZcWdzcHUzS1I3NnRqV0xtL3ZRcTAvajZiQWFKakRXa251YTRnMk1BNjRMWUZWZXhoMm8zYjlYRVBoaWlwUlZ4dXdiaEZBamJWU0hpQ3RVSGxKV2gyU0Zya3VHNUJWNjhOaHR4elpyUnRGMVRSUGFhck1WSmJMczE2bCtpdnZDMUVETFB1bDZhWjBlYVhWWGlrMlZZTEtiakduTjExalQ3V0V0UVl0YllqQUs1V3hPVXE2WUFERkhXQjFxbCtweVNTQ3krc1ZMTUZWRXZXMkJUOExXTFFpSHdOWUw5Nm5tTW9YeFpjeVlZUXhVTEhpc3dhVmJ0M3I1UVNqSVNMMlM4S21DdTVkclhzbnRvSlYyb1pYd3hCMjk5WU1YR01BOHlvdHNvRHB0WFJUQUZXeVVsR3Vyd1VxdXlUM2NuU2xBVzAxYTVVRmkzc2Jxb0NQd2dPbEIxaURTMU5BRmk1YnRGaUN5MnNUYXdFczR4NGp3THFoUUUxdHBXcFE2ZjdKTmo5WGdzcGFLNTJSWmFQa2RrSGhnWVZySVZnNmhTWURWblpFcElXTGgzWkd0OEJsQWN0K0FWcXNsd1dzRzFJWkd2bllFdTB0VUhrM1ZxSFMxRXpOQlZwckZZbjJMRmk0Rm03cVpBT2x0WWRWaWpIWnlUUzZ0YkFkZmVSMVJVOEptUDB5MU1aRGVjK2hhNm5Ec2c4c3VuSFJLS0d4VUttMXlicEFPNXlqVkVpbjBIbUpjVzV0d2szUVdVbHEzV0h0L3RUZzhuWmhYUlJjMmM3MlZzQzZWcUM4bXhhWnpDeTAyZ0xsUWFXN3FPcjVheU1rdGNMQTdvRHFKWEd0TmF5QnBYbERMK0JhY2k4WnE2VXJ4YW5naXRyMHA0YXJxeFg0ZVVDVjNxaGx1SWNIbFowcm9GQlphMUk2YkwyNDF3Z2FWVzU2ZzNHMTNvdWJOSEc3WCtzT28vTzNXSzBJTHEva3BnYVhWOEt6YU1DNkVsQlpxR3JEUHNaQTVWVVBaQWJlMnBLWXJMaTIxYUFXTEZhUndtb3hJYTFseHkxV3E1VFRzekV2cjU0ckduVlVlaDl2ZU5zWTkrZzFPaDhHYXloUVE2MlVkVlZUUStVOXNLeFF0V0RaOCtnT0ZRUkx5MmhLSlMydFlGbTkxUUpYeGpWNkhjelJGN0dGaDE2OGV3SzRkSkxJL0MwS3FxaVByZ1pXYlpKTERTeHIrZlFCYzc0V3dhSTd0RzU3Q3FzVnpjeXlqYkhSTHZkRDRSb0NtSDdtTG1xeGF2R25pNENxMVAzYllxMHk4NzIwV3JRRkxFQkZzQmhubXRvZFJvUFZzcU1oUzdxcnhYclZ4cURmQ2F3c1VGTkRwUzFTM2pDTkVsVFI4Qkp2bFpjRnkxc2NlQTlkbXlvSWxuV0hHUkcvMW5CbHJWZkdQWHBIRjBGVkEyb3NWTjU0SE52cEc4RlRnOG8rdkZKalF0Umk1bFdhcXBEV1RjYzlkOWhxdFVxMVU2VkpmYTF3dGJyR21uWkxnZFV5YkxZbDhOa0tsV2V0YXFPV292cjVXa05DdE1JZEFsYk5IWTYxV3JVR2pkSWt2aGE0TW9EVjNHTTNCS2dNV0ZFcmxCMzFNd1FxNy9kUjBEYmJpRkNibmVEcHJJdzdYSVRWeW96a0xzSFZxcnV5THJZSVZnYW96SXoxTEZRMm9lclZxR2V0bFMzVzh4bzdXOER5bWtJMTFXTGpXUmwzT0pYVml0eFhiYlJqclp5bUJGZHRnclBleDY3VlNxMDFWSkUxSyswRjZFMXFzYS96M3R1RFVnVjg1QTY5cXM4c1dCbXJWWnFNM0xKclJSYXUwdno1eUhwMXJVQU5oU3J5MTFGRFEzWVdhaWxvNTRFVnBYVmF3Ykk2eXd1V2pyRmFyUzV4RVhDVlltaVI5VHJjVjlnQ1ZLYnBvUWFWcFR3N1hIZXN0WW9nTGdsM095eE44MjVNbG50Z01YZTRHbFpyREZ3WlVWOXpqUkZnM1ZwQ2xiRldualdKckZWcFpFK3BoN0JVTlJydDVHVUZ2S2V6TXExWlk2eFdOSlpvU3JneXVzdURxNXNhS3EzV0hBTlZLZXJ2V1N0dmxsUnAxR0xKR25ybjhNQ3lBcDdWRGd3N1JKV2V0VUxBMHV5RzdDRGFsZ25IclhDVjRtaE5lMEtYb0xJV1lqV2dLZ243MHZ4MHJ4NXJMRmpVV1FUTDZxeldob2F4TG5FS3VGcDFsMmU5dXFGUWxZWmZUQVZWelZxVndMSWRLOUZnZmd1V3RudEZZT204ckVqQTF6cHJwclJhcFkwRjFncXViaTJnS2lVeXgxaXJhSHpRa05WbkJxeG9aV2pkNFJpdEZYV1JaL1hXVUxoYWRaZDFqZDJSQ3RVWWE2VXVlUWhZSmZIdnJReHJPbXNxcTFXRGE0cmRKWWJBNVRIUTNaV2d5b0RsV2F2c0NqUzdxaXl0REV2dXNHUzFXdURLREprZHM2dEVLMXdlQzkwVVVHWDIwMnVCYWlyUjNpTGNoNExGbGFFbjRMWEp0UFFsbXhLc0l3V3U3cTRNbGRmZkdIVUN0K2lyVnJCYWROWnF3bFhhUFdMUmNIVlRRRlVMdHJWQWxTbndpM2JFaWxyTXg3cEJ6OXJVVWpzS1ZvdlY4aW8wRm1HMXNvQ1Y5dUlwWlR5Nk1WQkZqWkExWGRYU1VsUXI4cXRkM3hUNkt0cGV0eWJnSTNkNHBMakVMR0F0K3lIeUdydlZobXFJQzh4VWp0WUdZcXcyV0txemRCZjcxYlJhaTRETGU5YmVkWFpyQ1ZWVVF4MkIxYnFsMnBRcndnZ3NYUmw2QWo1eWg3WGM1VkM0U21PVE1uQ1ZBR3VCcTZzTlRxdGQwQmlvU2tNbldydHlTak1LdkJyKzB2Vmx3V29SOExYV3FscUhkZ3RZdFVtQ3F3RlhkeVJBMWRybWxSSHRKZUZlQWl2YUNId29XT29PVjlOcXJTWmMzb0t0V3l1b3ByUldMVzZ3RkE5YkxiQXlicm9rR2JKZ3JTVmMzVkNvaG9RVnBvS3F4UTFHWUhsd1Jic3lsTUNLQkx3MldGaDNPSVhWcW9uMldtM1hvdUhxcG9JcTg2MWJMV3MxQkt3by85VUtWcVN6ZEV1NW9WYXJ0RVBZR0t1MUNMaTZ1eHBVcldCbEU5cXRZRVVoaHhKWW1TYUVMRnd0WUswRlhGMlVwVjQwVkVNRWU0dG9qeDVZQ2F6YVBQWUlyRmFkTmRacVpZOVd1R3BndGNEVlRRM1ZrV1N0b2hWaEMxaWxXRTFKd0d0TG1PcXNWaEUvQnE1U2VjNVF2WldGcTdzclFYVlhBTXQyUjZ1QVgyMnIxVEp6ZmdxWHFKK3BLK21xS2FGYVRXdlZDbFpVWDk0Q2xxZXpiQ1ByVUxCV0E2NW9BNm1oY0hVWnFNYkVxdGJTV3EwV1dKR0F0KzZ3TmlodHFHNXRCV3MxNE9wYW9ScnlvYWVDYW9nYkhBTld0RnRXQ2F5YXpySldLeXFwbWNwcVpTWVpMZ0t1L3dPSG9SbzF6SElDN0FBQUFBQkpSVTVFcmtKZ2dnPT0nLAogICAgICAgICAgICAgICAgICAgIFtbNTAuOTc3OTk5MDU1MzA0MDQsIDExLjcwOTEzMjYxMDU1MjI4XSwgWzUwLjk4MDY5NjY0NTY2NTM4LCAxMS43MTM0MDI0MDU3MDYyM11dLAogICAgICAgICAgICAgICAgICAgIHsKICAiYWx0IjogIiIsCiAgImNsYXNzTmFtZSI6ICIiLAogICJjcm9zc09yaWdpbiI6IGZhbHNlLAogICJlcnJvck92ZXJsYXlVcmwiOiAiIiwKICAiaW50ZXJhY3RpdmUiOiB0cnVlLAogICJvcGFjaXR5IjogMC43LAogICJ6SW5kZXgiOiAxCn0KICAgICAgICAgICAgICAgICAgICApLmFkZFRvKG1hcF9kMDVhZWNiM2Q2ZGY0MTEzODU4MjcwNjRkMDhkZWMyYyk7CiAgICAgICAgICAgIAogICAgCiAgICAgICAgICAgIHZhciBsYXllcl9jb250cm9sX2VmMThhZjljMjA0YTQ5ZGY4Y2YxMGRkYzJkY2UyNTI1ID0gewogICAgICAgICAgICAgICAgYmFzZV9sYXllcnMgOiB7ICJodHRwOi8vc2VydmljZXMuYXJjZ2lzb25saW5lLmNvbS9hcmNnaXMvcmVzdC9zZXJ2aWNlcy93b3JsZF9pbWFnZXJ5L21hcHNlcnZlci9tYXBzZXJ2ZXIvdGlsZS97en0ve3l9L3t4fSIgOiB0aWxlX2xheWVyXzA2MjA2Y2U0OGNjZDQzMzhiNTcwYzYxYmE3MjcwNDk4LCJvcGVuc3RyZWV0bWFwIiA6IHRpbGVfbGF5ZXJfY2JjMDVkYTg3NjZkNGQwN2JhNmIxZjg0MDJhOTI2NTUsIH0sCiAgICAgICAgICAgICAgICBvdmVybGF5cyA6IHsgIk5lYXIgZmllbGQgaW1hZ2UiIDogaW1hZ2Vfb3ZlcmxheV84MTQ2OThiZWQ5YTQ0MGE5OTVhYjM3MWQyZDkyYTA1OSwgfQogICAgICAgICAgICAgICAgfTsKICAgICAgICAgICAgTC5jb250cm9sLmxheWVycygKICAgICAgICAgICAgICAgIGxheWVyX2NvbnRyb2xfZWYxOGFmOWMyMDRhNDlkZjhjZjEwZGRjMmRjZTI1MjUuYmFzZV9sYXllcnMsCiAgICAgICAgICAgICAgICBsYXllcl9jb250cm9sX2VmMThhZjljMjA0YTQ5ZGY4Y2YxMGRkYzJkY2UyNTI1Lm92ZXJsYXlzLAogICAgICAgICAgICAgICAge3Bvc2l0aW9uOiAndG9wcmlnaHQnLAogICAgICAgICAgICAgICAgIGNvbGxhcHNlZDogdHJ1ZSwKICAgICAgICAgICAgICAgICBhdXRvWkluZGV4OiB0cnVlCiAgICAgICAgICAgICAgICB9KS5hZGRUbyhtYXBfZDA1YWVjYjNkNmRmNDExMzg1ODI3MDY0ZDA4ZGVjMmMpOwogICAgICAgICAgICAKICAgICAgICAgICAgICAgIHRpbGVfbGF5ZXJfY2JjMDVkYTg3NjZkNGQwN2JhNmIxZjg0MDJhOTI2NTUucmVtb3ZlKCk7CiAgICAgICAgCjwvc2NyaXB0Pg==\" 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 0xb336fd438>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "leaflet_map"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "cd97f46622934492bb28b5bd62183f3c",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "HBox(children=(IntProgress(value=0, max=1), HTML(value='')))"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Maximum at -6m east, 70m north of station center (lat/long 50.97998, 11.71118)\n",
      "\n"
     ]
    }
   ],
   "source": [
    "for xst_filename in tqdm(glob(\"test/*.dat\")):\n",
    "    xst_data = read_acm_cube(xst_filename, station_type)[0]\n",
    "    obsdatestr, obstimestr, _, rcu_mode, _, subbandname = os.path.basename(xst_filename).rstrip(\".dat\").split(\"_\")\n",
    "    obstime = datetime.datetime.strptime(obsdatestr + \":\" + obstimestr, '%Y%m%d:%H%M%S')\n",
    "    subband = int(subbandname[2:])\n",
    "    sky_fig, ground_fig, leaflet_map = make_xst_plots(xst_data, station_name, obstime, subband, rcu_mode, caltable_dir=\"test/CalTables\")"
   ]
  }
 ],
 "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.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}