test_io.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "dfab84bb-b622-4641-85ec-57cae732c600",
   "metadata": {},
   "outputs": [],
   "source": [
    "%load_ext autoreload\n",
    "%autoreload 2\n",
    "\n",
    "import numpy as np\n",
    "\n",
    "from matplotlib import pyplot as plt\n",
    "from astropy    import units  as u\n",
    "\n",
    "from common import hallmark as hm\n",
    "from common import io_ipole as io\n",
    "from common import viz"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "0e1cc6f5-7965-4268-bb9a-2266c9a7311f",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Rhigh [  1  10  40 160]\n",
      "aspin [0.]\n",
      "mag ['M']\n",
      "inc [10 30 50 70]\n",
      "snapshot [1000]\n",
      "window [1]\n"
     ]
    }
   ],
   "source": [
    "pf = hm.ParaFrame(\n",
    "    'test/ipole_86GHz/{mag}a{aspin:g}_w{window:d}/img_s{snapshot:d}_Rh{Rhigh:d}_i{inc:d}.h5'\n",
    ")\n",
    "\n",
    "for k in set(pf.keys()) - {'path'}:\n",
    "    globals()[k] = np.unique(pf[k])\n",
    "    print(k, globals()[k][:16])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "616a52df-ac34-42d5-9e85-23957df93f8a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>path</th>\n",
       "      <th>mag</th>\n",
       "      <th>aspin</th>\n",
       "      <th>window</th>\n",
       "      <th>snapshot</th>\n",
       "      <th>Rhigh</th>\n",
       "      <th>inc</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>test/ipole_86GHz/Ma0_w1/img_s01000_Rh10_i70.h5</td>\n",
       "      <td>M</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "      <td>1000</td>\n",
       "      <td>10</td>\n",
       "      <td>70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>test/ipole_86GHz/Ma0_w1/img_s01000_Rh160_i70.h5</td>\n",
       "      <td>M</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "      <td>1000</td>\n",
       "      <td>160</td>\n",
       "      <td>70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>test/ipole_86GHz/Ma0_w1/img_s01000_Rh1_i70.h5</td>\n",
       "      <td>M</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "      <td>1000</td>\n",
       "      <td>1</td>\n",
       "      <td>70</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>test/ipole_86GHz/Ma0_w1/img_s01000_Rh40_i70.h5</td>\n",
       "      <td>M</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1</td>\n",
       "      <td>1000</td>\n",
       "      <td>40</td>\n",
       "      <td>70</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                               path mag  aspin  window  \\\n",
       "3    test/ipole_86GHz/Ma0_w1/img_s01000_Rh10_i70.h5   M    0.0       1   \n",
       "7   test/ipole_86GHz/Ma0_w1/img_s01000_Rh160_i70.h5   M    0.0       1   \n",
       "11    test/ipole_86GHz/Ma0_w1/img_s01000_Rh1_i70.h5   M    0.0       1   \n",
       "15   test/ipole_86GHz/Ma0_w1/img_s01000_Rh40_i70.h5   M    0.0       1   \n",
       "\n",
       "    snapshot  Rhigh  inc  \n",
       "3       1000     10   70  \n",
       "7       1000    160   70  \n",
       "11      1000      1   70  \n",
       "15      1000     40   70  "
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "sel = pf(inc=70)\n",
    "display(sel)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "138402de-2a3d-4c29-8836-91fc0235ba41",
   "metadata": {},
   "outputs": [],
   "source": [
    "imgs = io.load_mov(sel.path)\n",
    "#imgs = io.load_mov('model/Illinois_thermal/86GHz/Ma+0.5_w5/img_s05999_Rh160_i30_hres.h5')\n",
    "meta = imgs.meta"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "103ae26e-5b35-4f03-bea1-d26a254de965",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAOUAAAEKCAYAAADtrJL1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAABA2klEQVR4nO29e7As213f9/mt7pnZ+7zuQ7oSegASWJhHwCIoPCqBwrYIiklFKCkcVaoCKRPLuExMVZyHMGWbBFwBAqFMJSYGFxUcx1ZIYiOVTZCRyg5FFQJhISEJSUggCV100ZV0de8995w9e6Z7/fLHb63u1T09r71n7+mZ09+qOWempx+rZ69v/57r9xNVZcCAAf2B2/cABgwY0MRAygEDeoaBlAMG9AwDKQcM6BkGUg4Y0DPk+x7AdWMsEz3h5r6HMWAAU+4x03Npb3/gSHnCTb5O/uzVXkQWfuflGEJSDyx+Q9/euX1QXwcM6BkGUu4b20jVAQ8EBlL2AQMxByQYSNkXDMQcEDCQsk8YiDmAgZT9w0DMBx4PXEjkIDCEVB5oDKQ8dLQJPJD04DGor8cGkUEFPnAMkvJYkRJzkJ4HhYGUDwIGgh4UBlI+aBgI2nsMpLwuSMt8V3+54y9yjoVzLrE9B7LuFQMprwNdhIrb2sTq2neT816WoI3zDmTdJwZS7hvbkHCT8+ySnAvXOLDwi0j/x9iBISRyFdjnRNgVyTe6lvQvBNMeU5/GtiEGSXmMuCq1duU1w+Tf6wNpCQEPTGIOpDx2XNbBtPX19kTOA5SIyzCQcg8Qt/0EUr+jSX5dUvS6yHlEZIwYbMprxkUIGY+Lrx0OZnfnWnqNK7Q5tznvAZF3kJTXiF0RKp5nJ9JT3PXYnV1JC9dNlAOxLQdSXhVU10+6SyYEpCS/FEGvi5jV9Q5Hau0DvVJfReTnRORJEXlfsu0HReSPROTd4fXnku++X0Q+IiIfEpFv3c+oLwBxy1XH+F362uiUl5zo1xlKGbASfftL/G/Aazq2/6SqvjK8fglARL4ceD3wFeGYvysi2bWN9DqxIUEHYh4HevVXUNVfBZ7acPfXAm9S1XNV/SjwEeBrr2xwfcFAzKPHofwFvldEfieot4+EbS8BPpHs83jY1h90ORXSSe9k9WsZ1kjNS3tpt1CbB+weh/DL/zTwxcArgSeAnwjbu2Zdp7dDRN4gIr8lIr815/xKBrk1NiHNOoJeNXEGcu4Fvf/FVfVTqlqqqgd+llpFfRz4/GTXlwKfXHKOn1HVV6nqq0ZMrnbAm+AiUmwZOddIzJ1gIOa1ove/toi8KPn4OiB6Zt8CvF5EJiLycuAVwG9e9/iuHQMxV2NLr3Uf0as4pYj8Y+CbgeeLyOPA3wK+WUReiammHwP+EoCqvl9EfgH4XaAA/oqqlnsY9ubomCiyJmanXXapE2jHJa871thHtH/fA/1NpPOPfsS4I4/qlbfCSyFSTRbJQsQmSK91hEzR+XfqShhYMgl3kv1z6UoHV5gcv0wytq/Ro/n+G/p2ntWnFibB4cr4BwzbELjz+F2osRdVCZepk7tQMdepqgeoxvZKfT1mtEnRIJlbMXF8/aQXkabE3FKNFSe7W22yCa7DO3yEOM67OiSsImT8PtlnQWJet+NnUyIMhLwwjvfODhUXWeq0D2Kuy93d5ly7uvaurrFnDOprX9AmYnuRsHOVKrugxl7kcrtQZXdZ9Gud0+cC6YXXqqrvEAMp94l1qis01wCuIuahh0kuQfCdLvzuAQ5Lrh8JFuzCBTuxHW+7mkl36JN5XY7vod7fQMrrRnuitL2wkZAbSNGNnD7rznGAE3clGQ/MfuzC4d/BoaJNui4SbrptFY5gkqboJGM7tS55f4gPncGmvCaoV/q6BHuX8ctNSLDptdae68geOBEDKfeJqH52xCF3mv54hQ6fi0iiK8kuSs95oF7XiIGU1wH1y5/qSxIDKu9q4nG9SmwqLfeqDnb9hl3jiZ7oQ/I+JxhIuQ902IVrc1svWx5xgwm6jJhbE/G61ErXfIhF2MOsI0R0IBhIuS8sI2GcaF6XS8srlJ69JWCKVXnEtGK44oB+r+hrYyDlVWNV/dcgMatJlU62zmSAHhUT3lSVvMohrNAuRAR12MPNCT1fadvAQMo9o5OQre97s+a1XfRr1a4XWby96TC6zp2aBF1axAHZlwMp+4B0gqfNX4O0tKf+9Th8OrGEjJdZ47ns2GVkXXqtZbHcff1WO8BAyuuCelABWsHKOMk3LWVxURX2IlKig4yLWUQXtCmXkGYrol/02j3HQMrrhtcFXlZox9rEgfPdXsSrlAbLyLhLEsRzXeQeusbRJ3v7kujVo2ZJL5FHReRXROTD4f9Hku8Oo5fIuskibtFGazt94OolQ3utohNExGoLOVev9bzsq3GfbmEh91K09+tqpd5OyDjANLtekZLuXiJvBN6uqq8A3h4+H3YvkVXSoTGxhHUhh43Vva6mQcsaCYUas5V0zDJ7bUOg1YPu9kin5+96rTt+4TLJPgfi5IGekXJJL5HXAj8f3v888O3J9uPqJZJOoo7Y226vtSzDKEjHPIfRCLLMPodXvd8aAq0iVX1T2y9LWzWGI8Eh3MkLVfUJgPD/C8L2jXuJ9KJtgerqNLZKBUsl2RIn0CpcVF1zgmTOVNUsQ7IsfHaN6uwL5NzqGkvuYxNytvdZtYJm2bkOxOY8ZEdP1y/f+aur6s8APwNW9/UqB7UW6ln6LEzJF72vMYnACeKpQyPJBNuoCkEX2k6cyn5sqc1R9UvOuykxF0IcqxxUm5J9G6l4gKl2h0DKT4nIi1T1idDC4MmwfeNeIr1G28nT8MA6wBuHuzJSVk3wDSTmAhkjIVNELzDUD4lN0CJv86FxQc/rqqp+8RrJbxKveWi1eg5BfX0L8F3h/XcBb062H08vkYScC1k+scq6S1THlpdxG7Wy8qiORkieIycTGI8gcx3e0VZiw7pXetyOEg3sfKsJuRYHorpCzyTlkl4iPwL8goh8N/CHwHfAgfYSUQ9kqCoCiwRLtlUSJlVj491dIEG9ulbwokoWvKpJW4VqjA2vpW4nHVOpuit0EbJDonemJKo/KM8rMPQSuVYEAsgoNwmVZTDKbXvwcrYJoao22X1p/6slE9h2X+0HNIipqskEdsF50yLiKrKlKt+6Sd2eQ/HYeFz4vDDXNlFfNyFkct3qd1FF5wWoR8uyl5JyWS+RXknKBxbiakLGIH2YXAK22kFl0baMEjLu3ypBWZExzyHPa/V3E+9klkqdEP5t22aRdPF8ayb+okNqjYRfRsgub/SqxkY9JOQqDKTsE9KsGTCpVpZGzLhPSL0Tv+jYqM4B5rwZjyFzRsosyatokzKdtK3eJel3mnJhE1W17UFmCTFXYNG+XpKp4+13odTqgXZoamvEQMp9of30Tu1LJxDnUyUNgyc20jOGSDqIWXlSxyMkBtbb54cOyadNkqRSuI1VYZf0u5SYcXwdXEmJ2lk6MyVjp6RPSB/Op1HdPzAMpNwHQlJ6dPhUaATHo0oY1FDngzRQ+18TNRWLX0qwTQm2Y4OQkYyNa8TxJGPYVNXrImRK4JT4C0WuFokinWFnFiWjtAgaxxF/lwOrMtCFgZTXCVWMAYkq6ZWQFVBvS4mTOSg91Up6DcpsXEFSHRJijK7t0EnIuMyWdHEcCanaTqSFY1rSEBY9twBZl6TNNg/qLyFj/TAK+3XZ3AcoJWEg5f7gOyZwLB2SqpreGzE1Q8rSQiQuEifxoK4iY1BJNcZBG57VJe8jEqJ2FqdqoPWwqc6R3mMgSraZc6iRatgRp412r4WOAiN9stztwJw8MJByP0hT7cLi5y5VVp0gmkqKaA9mdlyWhFIiGbNEXQ3naBBdFQ3Sq7Ltyg6Jsqa8hsgaUmVtwka7sS0lN+y2lXqOu+zexjU4WCkJAyn3BlVFKlUxURlbK0U0cybZsiCtAEQBV8U31xKx7dypiNQRLuhScaPndlup0yJzdeZ2uGWVKtsO46ROq/g5vY73/alpdEEMpLxuBK+gSLAlvVJVF0gTaxKnjGZSEzhFbtk5CySEJhEbYQcfnLgJISs1dInNWdmOK0IpXciyjsQC37hStxocL5c+VFoOq4jSh/srk4ebHmR8MmIg5T5QqawS4pBBHVXXJAvU9mCmDfVNnSxKjeqYKDGdET21C1Ph1eX8aV8fljhrAFaEReK5FlL2XEMjaJCzGlaLjMnv0HBaxfOr1ul0B0rEFAMp94RKWoZUOik9Kr7OffXB9mtITIKHNRWpiyS24xya5tTGgLoELyWuTkhICeFZzNBZVbt2lcOmHWrpyDyqrts+f8uh0yBj+p1q8EjHax5u0kDEQMp9IIRGtMRsJqicNOoFiSqZEzTLjEQjc+SoJKps4WuFUwTNwv5Z07FjKrMgzqEuEJTgbGoRTiUl+RIixm+7kg/S/1MJ2yZoO0WwceKW2i1NQmrmqvOIl/rBFHNf4aCJOZByn4jJ5ZQgghRNIgGQOzS8AFMXS7MLJROiz1ZzU2U1ryWklN6C8iKARzFvrpZYYoG0SNS2HZdJxzh8l5CpRUhphyQqskvtbIo2Z9u72yZjIiVTW7vxO5WlJZ4f2NrJLgyk3BeChNJoU5alZeVQQFmaE0cEHWX4kQvhEUC0+l9DGEQzk6BhPViYrJikhVC636RiZWcmZJAyIU5EVJlXEFPakjGuaCEQNiGNpI6XSNBIoNQh1PYYtyVkak+GeyOSMWw7tEXNbQyk3CdSNRYQVwYSeZvwkXBOLIzgMUEjzkgnTS9tPKeEFHZNJrjxWMEZqdV7JLh79YJkbIQeEjW1QdaoViYkrQiaqsopkRpZPB2ETHODg01eL9k6XLU1YiBlH6At+zJ98leSJzpopGHyVapqkI5SgpqiWvs/ItEwlVbDfu26Eyp0hicadqxXO09bPY2ruJx9XxPSti0QtJLS8cDFB8NCzLUVDqnIXZamaejhJqGnGEi5b1ReSW8TC5CigKKAokS8STbvjGTakSkjCpRhkkbi+rCvtEgVnDuaZbVJ2SJvV5EYk+D1mCsOxYdBfIDE7YlqrG2CJg6eiuCpStqVALFQqiSe0wcbe0kl+QPEQMo+IKqqXhFKtPRGyHlRrQ/UTPCZoHmqqoYJ7SFWvKu+ykICe9hWrUd2wkLNP0mIGci/oBJrsD0DAWvzUKEMslmDShslse1QE1J88PVkiNQxRU2dQLAY9uhIFYz/a1BZNcYqjwAHQ0oR+RhwF1sHUKjqq0TkUeD/BF4GfAz486r6uX2N8VKowhTm9KEoYF4gZYn4rJKSflSrpaLg5gBad2F20fET7c14/jRm2drkTJpqRqdK7Irw3HBi8VMfzhFfQRJrWO4pCoqzbWD5u6LVNvMEu4ZK3PCktr2uYVtzXGoSsgwPsSPCwZAy4E+r6meSz7GlwY+IyBvD5/92P0PbAF0xuRbUK+LDZJvPkXmJFKFkYgY+l2q5knlNpVp5VQXyA8nMEbTiWoG0PrP31XmDVBQvuFJtPzMUTRV2IYE+qM2RnIKCF7NpBZOeLtwTAi5au7ZgW31sjYBpBKmqGn+vJNWwQhnU1qI4OtUVDo+UbbwWq34H1tLgX9FnUsJqYoqr2ptblk9p0nKe4wqlxAjkgwpruyqoNNcRB2lXETJqfYlaapLXSFuOogQCVwRV0wvqFI8YN1ShDOT0RlhFIROzGzGCm80YJXRiY3b9Dk7DQBJCpt8n0rHxu5XeHFZFWTl5FnDA6XaHREoF/oWIKPD3QtXzRksDEXlB14Ei8gbgDQAn3Liu8S5HOx+0c5dAlKKEWYGbeaTMgrSs9qr/08bBduqY+BIlZkJWdeDHYbuCK4EiSEwxO0+8lRxxRSAPat5dZwQRmg8YCbZjpYp23Vr6fXwvrRtI7cfGw0Zq54734MvFDJ4jsCsPiZT/tqp+MhDvV0Tkg5se2Ku2Besgrq5ULiENrihx8xIp7c8VCeYxgngPEu08okSsCRbVWBXwo2ib2jEE6agOdGT7GRFrj20kLkhQa8Or1HBdDQJbNpeU7XhmhGtJSKhJ6kMmU1HCbA7zogqFAAefNBBxMKRU1U+G/58UkX+Kddha1tKgn4hSoaFr1mU3JFYVd85U2cxVThApPK5MFiYLaB7ms5q6CcHJEomY1VLS8mKhnMSxBOmoUIYiXFJCFgnmMNXVKa4QKBUXyVhJz4RYlYocbNmw3EwDiWOebnTSVM6aNLMHmkvY4gqX5KeSsoSiRIuYVufrpVrp73zAWF3frycQkZsicju+B/5d4H0sb2nQXyxz3QeyWksB14jLxSC5lIqba+Wc8VkgXm4SsHqNjXzFiVCcCsUNYX4b5rdM9fW5HVOObV+TntTnHdnx5TjYr5pIYSf4XCgnjvLEUZxklGNXpQJq8P5WtxU9tF0SsIptNr2vC9lEQW0V72FeoPO5OXliJk+quh44IeFwJOULgX8a7Kwc+Eeq+ssi8k46Whr0HjH80dl/wzW3h4luxwWVU2pvZDUFo8CqYprgJzXJIlJJGuOJ2XltY5qKCtm0Jv/sjlTeWYmqchHU3NLsTFdCNvO4meJmHhfil1L64BhqhTzigFMVNiVjlJJVLq9WsVstwitNQD8CWzLiIEipqn8A/KmO7Z8F9tSDYAdQ31JlqeNzTpoFlDHCpYuUK4ePUG1XpJKY5Qn4TNHMyCde7LMDHLhzI0F5QsWRbIaFSk8khEogP6sfBjH84apQiFQqbXaekU+V0X0bvxQh6QAf9q1LoKgkRbBSssZ0wlYCgRTBuVOUlS0Z2zhUqusRSEk4EFIeJVqSQVqV2qoqA41MFoKDxd5rTiORAIIzJ6ikxYltdGVUc5VyErNzBD8J9mNGpaLOJYRGZsFO9VDcjNeu99M5lR1KYWz1Y5jlwvxGxvieI5t5eBbcXHCUtmwshlfCfZLYzdW29J7jbxWl5Gxe25KRkEeitkYMpNwnUjW2C62SH1FSusLUxnkenCrp/M2hnGiwG01XVcDnaiQeqa2nFNCxorEgVylQSqWexrBHJKF4I5cU1AkDPkrgJBc2OJhmtwQpzd4c3SvJ7wvZWYGIooVJTxtHkIweW+TdLnES1GCCg8cSBiyL5xiSz7swkLIPCC3yGuha1OupQxGl5ZqmNqPZhWqvPEgxFO+gPPVWEc+LVagcK4y8qaRFXNpl6q86YOwrr647r72zUgaSVmSNYRHqYgYxJupgLkI5znC3Mm48SWVvKi6oteHeMqkXamdSh0yCg6eSkkVhaXUH2OJuUwyk3DeSZHQLyC9+F9cPRtXVlTVBdAx+ZHajhTwU8iAB5xKcNQoTj2QenWWQKW5U4gtnWTpeTFJ6sWOd2rrLuasqjlcJCHkwhUcm4FUEV2CeYTAiF/UtmLdXYAKuGJGfecafUxxlonsHqRzrC4VsHw1pfhQlnM9gPoP5vGlHHiEGUvYBUY2tFv8m6mxqX/ng2SwcrgyqZLAJyxM11fSkZpE6NQKOSyQzwvuorgpGwjIpQ5IrOva1izZ4af1EbXhjxc2MhIh5aSXH1FEv6AyymZLNTIJj0ZQqLHJ+x1FOxMj5nKngMZndhpzYk9EpFG3J8xk6L7rV1iOyJ2Eg5cWwQZrchaC+Woe4UKcGkNLj5j5IScWVYlKzAH/iYeIZnc4pS4eWIfCfKZOTeYgqCGXp7HUe1OWoao48MinJclvt4ac5elrCTUVnJlFl5oyQKrUaW9YOID8yqVicQD41qZ7NQzw1txUtPhPOnpdzIjBSS1ZQ7+tF141kAUvMZ16gs5m9XxbjPSJiriVlWB61Dl5Vn778cA4QO54QjfYFzpIIYtEsI4HHzZVsHlLeCsHNTe3MJiWnpzPmRUZZOpxT8rzk1sk5RZnhlfC/UIwyZuc5WgqMPNnIk+Ulo1GJiCK3p6gK3gtnf3SL/J5j/LRUJJTSHE7ZeR2miSEa0ZCgMBJcUTuByok9QLIZ+JGjuJlbca8yqKujDHUSHjyWxSTTGXo+O8olWsuwiaT8ZHjJin0y4At2MqK+o51tkm67JDnrWrBxqVaoZOfC2kS1ILwr4itIyqAeiig3JjOKPKP0QuaUcV5we3zOtBhReEeZl5TeUWQl87lJy5u3zrgxnnN7ck4mntx5cvG89w9fjHviBH2kgPuOR36vJD/zFWns3uHpL54wu2OEHd1TxnfVpGYrrc+PBBFzQhWnguYZo+fMlrVEexdqESkyx9LpplP0bFo5d47V45piE1J+QFW/etUOIvLbOxrPYWMXUlOT8hbxnKGUpCWYWyDdzSUscA5xRm9SbVZknIwKxlnJaT5n7EpO8jnTrKBQx3mZV8QTUYoy46HTKY+e3OeFp8/yhSdP8fHpo/yL934FNz885s7HPKdPerLpfbL7c8thjSU4Ql3a5z99DrnDjzPuvXjCcy/NOPmMkk+VfArlRCjHkN+3xIXZHWHyNORnyvzOGDfz5PfmxDS8KoMprJnUouj6pY4Wm5DyG3a0z+GjU0rG9VG+uU8jpYzmtjXQooAit1X+YZsfZSEG6EKaW8iDDeseye3696cmmjJRXnh6l7GzCf3o+D4Ac59xrxjz7PyElz/8FGNX8vD4jMfGd3kkv8fbPv1lfPCTL+T5vzbi5qfmTD4zJf/0szA9R8/OaBQ8DnChBV92espDnxpz5/dOeOqrHmZ2Wzj5nMcVSn4G8xtWDNrnwvwGqHNMnrFzlGWOHwc13StyPoenn0Vjb04tl3tcj8iehA1IqapTABH5DuCXVfWuiPwN4KuBH1bVd8V9jh5di3EjYhvx6vMF1Fz1qHcwmyF5HhbwKsQ45MhRjqxyuiiMzpTZuXlE8ebq9F44Ox9TlBnlHWHiCm5l54xcyUhKbmVTnpg9zEfvPY8vufUkzx/d5dHsOb54/CQvzu7zd37tW7j14REPfWzK+Ml7yFPPWAJ4Ys+1+1Rq6e0BUj6HFCdI6XneuzzlzTH3X2rrV82WVMvqyUJe7SiUN8kd85sZbq64ucedzc2WnM1hPq8eAlXJk/bf5Miwjff1b6jq/yUi/w62SuPHgZ8Gvu5KRnYoSPNX21JzGZapufF4r2hRIOcz5NYNk5ZxBUbuwpIsyznNZpBNQaYOn2W4kwJVoSgdz80nTFzJnXzKI/k9brspr5j8MR90L+ap2U2+YPJZXjb6NK+cPM27zh/lZ595JQ+/d8TDH54x/uSzyPksOJuy+qETyRmqu1cdlWPmzXQK8zniPdm84GbpmT1ywuyh3PJlJTxEQiaSz01ilmMhm5m66569D2fTuhBYLLP5ABAStlu6FWsufBvw06r6ZmC8+yH1HJtMhLjKY+Vrhd8slps8n1k5kMJX9qRGCTOW4KmE0XPK6K5D7mUU8yzE4ZVCM87KEc8WJzxXnjDXjMfcfV4yeorHxnf5vPxpXpA9x/PcKb87fQlv+ehXcvOPSyafndYlNvIMGY2Q8cikd56HZWXJYmxq6amlt8yb+1Pk3hnu/sy8saHkiCtgdN/S9/wYyrFQnAiz20bKydNzOJuadE7itcecLNDGNpLyj0Tk7wGvBn5URCYcyHrMnWOFGhtr7MDiRGp+FxM+O09SxejkfGYB9ptj8387i/XF5VnqwrKp+/blPJtQ3p7jQp+PQh1PzW4ykhKH8ml/g3t+wsQV3HFTbrs5npwP3HsRZx+/zWN3C2ReWnZNTl102cekA0vtUYmpbiGDXX29UkUV5jNUPTIdkd+bMz5xlBMXYpZS5dNOH3VhhYkyfrYgf3pqxFe16gILv/1xe15hO1L+eeA1wI+r6tNhpf9/fTXDOgCkxKyLqjZ2SUnY+gJxlu+aErdBWlWIKqxzuFmJ5kKprqqxEwtfqViQnqngn3MUkjP1wlO3bnA+ypnkBafZHCfKH8xewFPFLT47v8mz/oSn/ZS7s4I/fO4R8nsW+0QERjl4b/YiUDFOxB4WmUNVAllLIKvJGZv4lLYoOXvmjJNSmT52AieO4kRsEfXE/s/vwcnnPNn98EDIMwsBtb2uR5zFk2JjUqrqfRH5l8ArROSbwuYHw8GzDB0SU30o67Fs5UdE+F4ymqROl3BFUhQl7v4MzQQ3crgTR1nVzQmHlpYxk00FJKMshScntzk5nfG8m/fJpaRQx8ez53PucwrN+MT8eTxd3uBTxUM8cfc2bh5qxo5Ctk8hwfOpdUuFqqJcacsrHVTJ9BIkp1fb5AKBpzNTqR47CQulFZ85yrHZw6N7yuTpAld4WyWSZ5ZUMBrZ50jOaNceMSFhC1KKyH8OfB/wUuDdwNcDvw78mSsZ2aEgDX0Ep48R01ckWwWb9Ildli5dUgVfVi0M3LQgzxx+4lBxFowPqzXcHNDwBxUAx/xzE4rznEyUaZHzzOyUuc+4M5rycH6fJ+d3+Fj5fN779Iu5+8e3eegpu6wfZ5bdl4kRE6D0VS+SdHziLVfWGuy4kC/r6wdWUuZjftNyX8WbJxaE8XPK6L6vm+SOMmRmxETV7NYsA59IzSNLq2tjG/X1+4B/C3iHqv5pEflS4L+7mmEdIKpJGLJxfFBRfUv6RSTdmKtvAiklb/5ZtCjNuaKKm5W484xMhHISMr4V3Cws5/J2RkssyCgL4Wm9ST4pyUcld88nvODmc5zenvF0cYMnp7f48B+9gNNP5Nz5RIGbW/FlHWdI6XAukNB7ZC6os8SGSmqWskhOMqo2CqNRI49XwwqQ/EwZ3bdYqyuUcpIhc09WeCNkWDtZ/RZliFNGaXnExNyGlFNVnYqtip+o6gdF5E9e2cg2hIi8Bvg7mML091X1R/Y2mJbUNGdOXHEfqtNBUxpGOBckpUCeh+JZWXU+nReWNZM7W10x84yfhfktxyw4TiiVzGNlP3KQUqy41XxEcdscLfNZzt37Jzxx9zYA5/MR/t4IV1pQ3819VZEdLEdVBEQzG3fI6NG4kiOQFG/OKfFBSibFs3Qywt8Yc/LZOcWNjLNHM/JzI6OUFt4pTyCbuirsQ+aM+JlD1FlDIh4MYm5DysdF5GHgF7G6q5/DcmL3BhHJgP8F+BbgceCdIvIWVf3dfY6rnii11KwiSk5s4kZiBilSEVJqqcJ4VLUSl9KHHh6KD81+wNTAsVrlOh+q2rl5SFQvzMbMzoT5TChPMspTT+lypm7C5M453gv50xlSWsbNZOTICo8rFm03rXJTM8Q5I2fo5yFV81ZXq6zOoacT/I0x5S2LnolXXKjKl52HaghOKHOrkudGGUwt5Y7J2JxF4nCTia0UKYpFYjZ+88PHNo6e14W3PxgcPg8BH7iSUW2OrwU+EgprISJvwloZ7JeUEW1ylgRHpdlLlUpbqXehtGQsxTE9R04mITHdPJJuWpBNLD5YjIVyZGpsNleyOejUYn+aW9kNl1lh5Wxm2TPzWxKIALP7mekXtzzlvYxyBMXNDM3FPKGlGjljHm70R6la23cVyJ1VbtfgbY2EHOXoKKO8NbF7nXn8SYaKkE+NmBHiTcL7kaNUyO4Scn4z8/J6qX43gUViwlFJzW0cPX+ztSnD1MaX73RE2+ElwCeSz4/TkWG097YFy8gZy/9nWegX6ZOV/sEWnc2NvHlmztbcITNvamoo72jqqlR5paJAAc6HOOZcq8ro2blUJUSsPbulwJUn8OwroLiZM7qr3Pi0S6rUhT6Zcx8kotbk84SeIEbUSNz582/gxw4VyKcl7ryspFoVzwxV9iL8xBLcxy5xkE3GiAg+xixjoWrXUX3gSIi5jfp6L3k/Ar4J+IXdDmdrdLk2F/4qvWlb0KHWRjvJVEFTaSV6HHOLFUbbUuYFUpTIvMRPRvhHT8JSrlCNXELburBwuRwLsXpANg2hiKQsZYxz+lzI71MRobgp3J1kuLmdL59mlcppi5s11NCxe7L6PGYjWvkRCRXxPO68rLo1u5mFPOanjtEZZOdWByj2OClHAmOhuD0hEyE7r4koeW5xy3lHQsGRYRv19SfSzyLyY8A7dz6i7fA48PnJ55eyZzt3KyTOoDRWqa5AsgwpS5QJ4hw6yoOX09REmZdkZ4WVBpk5zh/JLbk7KZicnRtJNNTSAcjPQx6pJF25wnNMlGpbLJ+jztY+Wn1Xa1fgSvvfCiyLLXQGyhOzJ6VUI6OaHaq5A2eqcTkJtvC5J79fUo5diI0GCSrK/LatkHHPOWuc6705vzQ8uLoWOx+BhIy4TDmQR4A/3tVALoh3YskMLwf+CHg98J/sd0gboCMbSEuCkyd4a2NamxNLV8szdDxCT8ehL0hW5cMCVQXz2Gagjbj637uaXBqdux2d5NqlK+M56rbqdpwUWl0/9jtxsyrfrvoeCA8NIT9XsqnHTUtLGQz34cK+xanDFTnZzRNkNk+82tE5JuZQomvgh49tbMr3Uv+5Beue/Nm4XVW/avfDWw1VLUTke4G3Yjbuz6nq+697HDtDGuOMFSdjulpIt9NRbt6bUI7R54Ifu6oIVdr2Dqib/aSXyWrbU0rI55FEyU4pISMnlCAdzU6NNWHdLNiaYvuqE1NZ1VRdKVtPCSX0RSkRn9d9MbGxl2OhOHVktyeM7k2ty1ZR2O8TQkfqCA+v1pK5I8A2kvLfv7JRXAKq+kvAL+17HFtjWVJ7bByb1ntVb0u55oUF1kdZI1Mmti2QUo3MgZxVIeU2JyqVFqvTKtVpwg7U1dCjrRi2V+csfFWGJAt2Y+y4hQju3Fa3SOEtJHIjZ/qQZfxMntWqGRBQlc608pKhorsTips5+Y2JEXw6tQeUX0LAI3HywHY25cevciAPJNJJlBI0y8zZk2T9ABYPnM0tfulctdKiWmwSSZGSKlQ5TyHhn9g2XbOQCdTat+pdkjh1qqFnoVeIV1tXWSrig4dVBHdm/SNlOkNvjpndHlkDochvZ/amqFXPk7IeU3ZuN1ROBH86siT1Z1fEI9sLzA8ca5deici7drHPgA0QpGQjJS8kFEieG0mtJJ1NVG+1Ua16XEgeCEkFUWVs2pyx4jqV1IteWp/0E6na1yXGikrzHLHlncT6rLE6QOGRWRHWgVrPDz/Jmd1xZOdUqXWi1CGSeM1g62Yh20ed4CeZEXM0qrzQlmebkPOICAmbScovE5HfWfG9YIkEAy6CVYudoU4oiNk/WZ3O5mYF/iSrClNFYsauy34U1lpOpWFjRqL6UR0SycJ6n2yWNIWFirjxfewbYjHSELssLXYpPqToFaUREtA7N5k9PGL6sGN81yoluLnW6zTTW1WgBDf3eFwI4ThkkpOdTEKl9PNL/NiHgU1I+aUb7HOcbrDrRLQlo9oq0ky9S4szxxQ9TyWxypGEZq9Sqa6aRyGiuJlVJI8JBT495jR4RjO4+aTHnScZN7F5j4IrfSBforZGbmWCTItGVo9Oxtx/2cNMH84sUSCouZVjCOoW8L6udyveavWAJR/oyNL1ZF6g07oqQpU8cGTq6yaFswZb8qqRELJqgRezWpK1lakzQ3NHeWdMcZpRhiVcrqAhgcooETOqPh9oiFHGtgI5zB7WOmYpGZNnlBufmgXiUaumSiierKailnW/ScAkZCiPWT72ELOHxtz9gtwWYc+06hammVgR6NTJpEAR1PHg4c18WanK/nSEm42R6TisGDnespND24J9oa22htbqUVLWvSpdvWIkz9BRjr8xwZ/mzG6PmN9yzG8I85uW0+pzKnXTnceAvFCcWpuDbJp4XgOfyoni7xSMb8146qEJo8/lPO89Y06eKhg/MzP1NHbAConq0aZteENV8bdOKG9NeOZP3GB+0yRhfmZLtfKpdXrOzn14IEglkWMhMKkeAFo/QCQsvJ6MkdPTqmI6xfFIxxQDKa8LK0pTViEQl6iplSobpOUot+SBSY6/MaI8yWzZ1k2huCGUpwnZAvI5Vfvz2K0r2orZPMQYS7j5Ccf0sREzgFyZ3yl55k/kzD854ubYcfLpqdmO89II6n3taBGx1gpZRvG8U84eG3Pv87JK8sY2B40sIKg8txpCKOISiZtyLUr0zOHHOe7GCfLcPZOWR6a2RmyTPPA24K+p6nuucDzHi1U1Y1uoJWNeeV01c9Xk9yNbG2m1bswujJLPhR4fUlgF8phyFyVPDKOYEyWk3p0pp58Sxs+Muf+SEh0pZy8pmN3JeO4Lckb3bjF+Rnnoo3OLSYbEdB+qon/mqybM7sDoOaraQW4eHwShfV+hZKHlgZuXVd9NI6lapmFm6XhVPi3m9KliBLlDJzluMgb1yGz2wCcP/DfAT4rIx4G/rqpPXNGYHgwsq+HjPYxH4DIj5HiE5kF1DesZy0moDJdTTdiqrpU3QriZks3rfpYVtCZqfuZtdUlWJxucfjKjuKHMXjSnfFgpH4LzQsjOHNNHx1WSuniCYwmKUw2OnLoTtJTWeevkKW/t8c4VF4s2L3s4qVqvTI3jVHva+BDWERDn4MSWg3HvjIaP8UgSCLZJHngX8GdE5D8CfllE/gnwY6p6dmWjewBQVbDzilJa0gBU6qzmWVBdc/w4w09yihsZxYmEoszJyaIUrJZ01VKx2kUs8O9D+Ujx4GbxYCE/AzcXstmY8kQpJ0p5y1Pe8Ewfg7jqhIp8tpjajjEJ7QqYPB2JCNnUk98rao9rjKWmsdAUMf5aalVWRENVeAAd5Ug5RsYjmGF5w0ckLbeyKcW8Dx/CKqP/MPAXReT7VfV/v4rBPZCIPRldJKd13dKgKpanWRW+gFodhSAli/oVCRlDJBCEiagl3nishMhcLd90AsWNQLSSqs2elBLUUiM7GvJmC8jOxcItRaymp4zuwejMWvaNnistBFJqvUYyKagVl3XFbtVI8PCq/RZWZzb8NlEKhgeVRImp3oh5JNjGpvw14IuA9wPvAP4z4IPA94nIN6rqG65khMeETe3KdmA9pKT5cUY5DgWNQ4pclFhRSlopkERKVvmqWicChIRuH/JkwSSoeULtu3JCRfzsPEjk2OZdqTo6Z1MjYnwIZOdKfq7k9zzZzJOdxfKQSfijSs3zNSm9mM1YJSpo8/8qDdB+Q80ccnIS+lYGT/CRSMttJOX3AO/Xdssl+C9EZN9lQQ4PaQ+SZd+nHzOH5rFdAbV3UxVVKwNZOXGS/NKYUxdD81W6XAblSCrHTJVKJ8EuLCArwVyyyYMknD87Nw9qNjPb0RWmNmczSw7I79eLm22gUsUfKzJGBxSEKgRalZCtpGlIfqhDREkcd5TXzrCyPBppuY1N+b4VX3/bDsbyYKC1ltJKUUprF61LOIaJqbmEgllJFkyABCFYJweEhJ8s8r5FyFhhPdTqiV5YCGGLGYzvUnlrq+VfyQJqa+8evKrTkKkTsnXceep8Cf97rYtxhZS8RvU/L3axlrCrSJsmw4uEMFLICR7l1qFLlGNg5k7ilLFw1YCLoSre7K0lQUpIKayMf7UsCqr/YxI3QQ42VNawn5J8Jkzy4M2UWG1jpg2niyuU0VmQamFZVXVdrc8jIU/VxRo+hTenTGO1itbJBykZI/lSiU/tzKn2SdXYlLASnD+xbIqTo0n2HJIH9oF1tmUMzpdlcPzUszEuNBYJkzCqrEV09GhFzLgca9HLqU0bL857aZE6IJv50KQ25Kx2LAdDsKybGA8tEmdNQkq7XiL1vFjhLRH7Pr5vrT7p7IWUFK+WzApAqz98Zg6k7BvipC9LW9is42SC2i5SgsNsyUisdHVHmiwQ919AS+rRIm5VZSAmnQvWaTlxGlUFtOJ+RVmTNo7BN23HzgeStwsKgaQk+1bjiQ+DdJASVs2EQmNlyTHEKgdS7gtrpKXGcIEP3sbKDk3UVtXqFFVBq7k2iJOW8mijSnkrtZNsduLwf/xcSd7EaROIWKfgtUjYde9xDNGbWgJO64T8tIRlvIE0LzYipCZWzp6lv+jhoPekFJEfBP4i8Omw6a+HEiCIyPcD340pcn9VVd+6l0FeBbw3D2wW11Pa5mg3ooCryRZXiMSVHQuSLmKZpKtItDgU0eQ8kXSJbdhQMzeRUmnIJ5YEAUu1E0yFjbHXytETj2VRYrogMWPHafxBS8vekzLgJ1X1x9MNIvLlWPW6rwBeDLxNRL5E9YDcb6qNaEOFpNZrTB6onCbB8ZL0ErK44Tyuc6SSZlLoovRrSzpNHDHQtPsSabWgPvo1JFxlM6fH2LLJ+jBoOmyiRpGOoWrxnownyyrb0jJ8DmcatHEopOzCa4E3qeo58FER+QjWxuDX9zusLbEqXpllIRE9lt2gKpMR090qOzIuPi5rx061FKqsidiQnlVWjS6+T1sQpOpmsjh5AWkscVNJlcQiq3YNHftUV0uWi6VxUIkJBZklqKserm15KKT8XhH5TuC3sJUqn8NaFrwj2efxsG0Be29bcAFUrv48M0kJwQkkdTnHUCs1xgjxTSdPJFijUkDL/mssxUpI2CjXkcYT4/i8dhMz7t8lKdeRpOy4ZnXBhGSqiRe4lpSU3lrvlf6gY5a9IGVYFvZ5HV/9AJZn+0OYtfNDwE8Af4Fuxa/zr96btgXLoB7ImlLTidV5FTH1rkoup6r8hktU1kjIJC7YIGPMoCnLWhUMha0ak7wLCcEWiLhqxccm25LtXXV72tddOk6RsOQtvNQfbNZdL0ipqq/eZD8R+Vngn4WPh92yoIWqLTvU3YvHo3rSR5IVHhFz/NjiZa3VVw0VyeMqi1glIPWKVg4arQkZVdU2WoTbmJBhvI1d10jexv4dYxFdsm97PFWNI3ew4ZFekHIVRORFydrN1wEx3e8twD8Skf8Jc/S8AvjNPQxxd4jt2J2zVSJVfZ6giuJx4hD1FuiPdmPRlIixM5b4hJSV06appq5EVEMrSZZ8tcT+W5B2a8jW2N52IqWIDXdXETIWFIsxywNF70kJ/JiIvBJTTT8G/CUAVX2/iPwC1ouyAP7KQXleIzrsr0bdV++RssQV1g9StSmxKhImIYtGWtsaAnaqjKvIE2u1tlsRbHKO9rb259Sj25am7aY+y+zZpPiYHmh4pPekVNX/dMV3fxv429c4nKuBmhS0VfbOEgfUm81XlJbpkmWVCtfI0PGJUyYJezTS1KKEpP5uwVmyToqlnzcsa7KUdKv21dq7WqFdKV6kWvzc3m5eWFP/xRUH6evpPSkfGKjaukBRc8akIQnMFlQc0k4CrUIKhoqIXU6OdYRcWMe5hROnC6tU2a5tkYwLZPaLxIQQ44wpTRLS9aQ2AcQdpBd2IGWf0XKEmI3YIkqq5nVJtTTmWJ2nRchlEqwj82YltrUnW+Up7b8OO7W9b7Qd4zXTscVC1pkLJUTk4GKWAyl7hqpKetaceFKGcoxJ1Kex2iKNITZUwSYxVqqsq5woq1TPruMr6dshstNxNyRlxzWcUK0vrfZLJGdbla0cPVZKROXw+lgOpOwLpG6rXpWYjBMt2FmNrJb24Wk19U3QtjPT/8N4GvuuO1eKDglYf1zm3FkSVPRGPm2n4nURM0UWauiqmhp7QCrsQMo+QCTYP67K5JFURYNmvmdE6nTRRFKk38f/21KykYS+Qqptg00J2NhpiR1ZH2S/TTy29aBqEJP4NkjILDOV3x2WCjuQskeoJGS7N2WLPA1va9LqrvOcHV7XRqwSljtYNkSXHbiSgM2Du49J7yfNdPLW90Q2kOTW8Xltt8feYSBlT1DVnJFmf8qqgjh0q5kdgf2VoYiu0MMKJ0snVoY2VkjbTckbt8cEhcaSGFeNtaHGts+ddsOOXtgDiVkOpOwLYlqYa0mAtn22LszQFefrQpezZSNnznbn7sQmUlRamkJSbAzvlqux7XOE7mWicjBm5UDKPiHG1iDkqjpWSsm4X7K90eV4XRnLNnZBuDbWEX3VNVOvajvDJ1Vjlz2IYldskdBH9zAcPgMp+wBxDQnZUM0iOuy+hXBC22nig8oGm8UZL4tNJC1sRv72ORvkrJ0/DWJ2IU0kOJDQyEDKfSPmkjbyXcNkS5Kql8bzOonI4sSPp1qVtbONVL0INpWKcd90PMu+93XD2/gbLtjGwSw4lESCgZQ9QOrkqZBMtvh5YVIvI+ISidXpyInlP0SAMgnC74ig66Sibz1sNiFm6vyJ4ZKQZLBw3mBXHlIiwUDKfUPcIgEqcrV6L24SeI/lOpZJgyVOoGrv8CQQaamNcaybYhsyrjvPOmJClWTQBYk2ZbvBUE8xkHLPkNR1H5FOuFVpaJuScU2qW+PaYV+bxM6Wf1V1LDvO0/CS7oCIXas/ljmsGnm/zZDJwlicC3Vh+99kdiDlPpGSoe3SX2dDdny3mEGzZPJtHI9cQVAIEmtHErHrmBUJBI3v2yGTNmLX6Cgxe46BlPuGtGxJ2DwbZhUhV15TWAgxrEMkaPrw6CLpqjo7K8a48Bt0nqBDlY3X3SQxwVkzJWvJ3l/bciDlvhAzTaLLPiKdeGuIWB2yNBzglkvLbZLXG9ePIZY6s8ZO103ITR8Wqtr9cOpSZVOkebGdJ07GG9XXnjt7BlL2ER224tJd1036qjzlju2o1vrGjcjXGeCvH0idxFyH9LfaxhHVY2dPL7J1ReQ7ROT9IuJF5FWt775fRD4iIh8SkW9Ntn+NiLw3fPdTsvVfc/8QJ93SIX2twFYqa+i5cSXwfrPXqmMDNlpZsgzqF18d52j3A+0bekFKrELdfwj8arqx1ZrgNcDfFZHY6/ensQLLrwiv11zbaHeBjqf6MpKpaufrQkjJGd+nr10gXZGy7LXy8A5ibvpah54TEnpCSlX9gKp+qOOr1xJaE6jqR4GPAF8rIi8C7qjqr4d27/8A+PbrG/El0SXUoyNlV+Rbh6skYIplUrK9b2ufTe678zdqk3MZUa86e+kS6LtNuaw1wTy8b2/vxCG2LbhyXMTGvKjd2N6+4QOhy8ZcpU1A4sVtk9G3bc/+9rK8NlKuak2gqm9edljHNl2xvRO9bVvQflqvIsqqibxs6dKyfTfFlvHMjZGON40xtjy79vWGmUld3uB15+gpro2Um7YmaGFZa4LHw/v29oNBw9kQJ+m6luur0DGhtzo+4qqI2HV8FzHb32163SWhmmUQF9ZX9lBa9lexNrwFeL2ITETk5YTWBKGNwV0R+frgdf1OYJm07R+W2TNbOkQ6sY3ns+u6m55z5ek2dEytKLDVuE77uqt+p66xxf4pbfTUruzFqETkdSLyOPANwD8XkbcCqOr7gdia4Jdptib4y8Dfx5w/vw/8v9c+8B1grWp1UXJucs5l574kAVft23mt9ri2+b7ru20eSD2EHJq+fVnckUf16+TP7m8AIZPHSkK6enV8G5vYiBdtH9DGiol7VfPjwve8bKzpsct6Y3pPVZlBPdq1MPwa8Rv6dp7VpxYG23fv6wOB7hSzDbyV206ky5DvIonlsDQuuPKe0/vdVNK1bdRN0cNVIwMpe4I2Kdqdtxq4pKd1JQEvSr5NzucWwxudErNj3MvGvPA7beogaufA9sjhM5DyOtE5AZOVDgm6XPv1Mds/2ZeHFjaYiJtIkk2cJh33etkQxmJ8sqVhrPitbMVIP4iYYiDlHtCZe7lEolyJTbcrIq7afxVJO1Z/XEh6t4uNrTIB1l2jRxhIed1or9TvmrxdS5Yui4sWUG7sVp9jbVL3OpIu0RAWvl+F1jlWrTLZiJA9UWEHUu4B6nX9pN62DV3XcasHscEuK8Ib2668aJfriLjIeLvOkRATtlSHe+bsGUi5b2yyHnAXds8FpOF2pw9k2EaCrrNDl425i+At7WJTG7qPduVAyj5hmTS5zLnW7rbbCXkpFXf7i62uDXugGEh5Xdiqd+QFyXlJlXTX6LrWzhcYb0PMVQXIeoSBlHvCRnblOlVvzcS6FAG3mbRbPDzW2aKrxrz099qEmOt+i2hX9sDZM5Byj7gKVe/aiLjsuC2l+zbj3dhuhd0nQVwjBlL2BFsRdM3xK3ba+rwXGEjz83WtxNimw1iPVVcYSHk92LKm10aqbbLvii+3uu6VYNkYLkHWpb/PJsTcKDtpvyrsQMrrxoZP9HXEXErGPhBxE3SNc0vbdGtidlyzb+EQGEi5H2xBTKjV2d5LxctiS6/zSmIeMAZS7gtbOEeOnoxtbGEfbqPqdx3bR/Si8sADj7Rw8LbHHSu2uLe+kuuiGCRln7Cp+nYhAq+L07Va8S377jqxpcSE/lc/3wS9kJTL2haIyMtE5ExE3h1e/2vy3cG3LbgQtiHkNgW4Vu17FXWCNsWWD6BjkJq9ICVL2hYE/L6qvjK8vifZfthtC1ahS53dRsW9KhLti5wXIOY6cvaZvL0g5Yq2BZ04+LYFm6LdqGblvltIxUuPa0/EvCA52yTtMyGhJ6Rcg5eLyG+LyP8nIt8Ytr2ELdoWHDX2qVruA5dwbm0iQfuAvrcteAL4AlX9rIh8DfCLIvIVbNm24Ch7iTxIRGzjErm2G2OPWT29blugqufAeXj/r0Xk94EvYcu2BXvvJdIuy7+L8w0wXAdBrxm9vgsReSz2oxSRL8IcOn9wkG0LdkGkB01V3RYXjff2DL0g5bK2BcA3Ab8jIu8B/m/ge1T1qfDd4bUtuIwzpm9k7Nt4UmzjIOshhrYF+8Y6tbbvf59DCQ9fJCHjin/7ZW0LeiEpByxB3wkJh6NSr5KaPZOoQ5rdvtHlBDqESd5GHHOfJWfPyLcMAyn7gEMk4TIcAjl7joGUA64G6YNmIOhWGEg54OrRl1UnB4KBlAOuH9uq6w8YiQdSDug/NiXxkZB3CIkMOB7sMjyzR+fbQMoBx4dDiZ0uwUDKAceLAyXnYFMOOH5sE57pAYkHUg54sLAsPNMDMkYMpBzwYKNHZIwYbMoBA3qGgZQDBvQMAykHDOgZBlIOGNAzDKQcMKBnGEg5YEDPMJBywICeYSDlgAE9w0DKAQN6hgeuxKSIfBr4+DVd7vnAZ67pWteJY70vuN57+0JVfay98YEj5XVCRH5LVV+1fs/DwrHeF/Tj3gb1dcCAnmEg5YABPcNAyqvFz+x7AFeEY70v6MG9DTblgAE9wyApBwzoGQZSDhjQMwykvCBE5OdE5EkReV+y7VER+RUR+XD4/5Hku+8XkY+IyIdE5Fv3M+rLQ0ReE+7hIyLyxn2P56IQke8QkfeLiBeRV7W+6/xbicjXiMh7w3c/FRoW7x6qOrwu8MIa2v6bwPuSbT8GvDG8fyPwo+H9lwPvASbAy7Emt9m+7+EC95yFsX8RMA739OX7HtcF7+XLgD8J/CvgVcn2pX8r4DexxsaCNSn+965ibIOkvCBU9VeBp1qbXwv8fHj/88C3J9vfpKrnqvpRrPv0117HOHeMrwU+oqp/oKoz4E3YvR0cVPUDqvqhjq86/1Yi8iLgjqr+uhpD/wH133enGEi5W7xQVZ8ACP+/IGx/CfCJZL/Hw7ZDw7Hcxyosu8eXhPft7TvHUM3uetBlexxiLOqg7kNE3gZ8XsdXP6Cqb152WMc2XbF95xhIuVt8SkRepKpPBHXnybD9ceDzk/1eCnzy2kd3eRzUfajqqy9w2LJ7fDy8b2/fOQb1dbd4C/Bd4f13AW9Otr9eRCYi8nLgFZjT4NDwTuAVIvJyERkDr8fu7ZjQ+bcK5shdEfn64HX9Tuq/726xby/Yob6Afww8Acyxp+h3A88D3g58OPz/aLL/D2CevA9xRV67a7rvPwf8XriXH9j3eC5xH68Lf7dz4FPAW9f9rYBXAe8L3/3PhIy4Xb+GNLsBA3qGQX0dMKBnGEg5YEDPMJBywICeYSDlgAE9w0DKAQN6hoGUAwb0DAMpBwzoGQZSDliAiLxMRM5E5N07POepiLxbRGYi8vxdnfcYMZBywDL8vqq+clcnU9WzcL7e5sr2BQMpH0CIyL8UkW8J739YRH5qzf4va1VY+K9E5AfD+18UkX8dVvG/IWy7KSL/XETeIyLvE5H/+Apv5+gwrBJ5MPG3gP9eRF4AfDXwH1ziXH9BVZ8SkVPgnSLy/wDfDHxSVb8NQEQeuuyAHyQMkvIBhFrVBAH+S+D1qlpe4nR/VUTeA7wDW/L0CuC9wKtF5EdF5BtV9ZlLD/oBwkDKBxAi8pXAi4BzVb276WHJ+1E4zzcDrwa+QVX/FPDbwImq/h7wNRg5/wcR+Zs7GvoDgYGUDxjC4uv/A6tFc2+LynpfKCKPiYjDioZlwEPA51T1voh8KfD14RovBu6r6j8EfhwrMDZgQww25QMEEbkB/BPgr6nqB0Tkh4AfBd66weGfxYpFfR7wNmyR7z8EchH5HWzt4TvCvl8J/I8i4rH1pn95pzdy5BjWUw5YgIi8DPhnqvpvdH2+5Lk/hpV0PNb+lpfGoL4O6EIJPHQVyQOYPep3dd5jxCApBwzoGQZJOWBAzzCQcsCAnmEg5YABPcNAygEDeoaBlAMG9AwDKQcM6BkGUg4Y0DP8/9LvpspThG9mAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "ax = viz.show(imgs, s=0)\n",
    "#ax.set_xlim( 99,-99)\n",
    "#ax.set_ylim(-99, 99)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "51d32862-7905-4e55-bacf-399bfbd767c7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$4141166.6 \\; \\mathrm{M_{\\odot}}$"
      ],
      "text/plain": [
       "<Quantity 4141166.64059526 solMass>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/latex": [
       "$8127.0011 \\; \\mathrm{pc}$"
      ],
      "text/plain": [
       "<Quantity 8127.00110226 pc>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/latex": [
       "$86 \\; \\mathrm{GHz}$"
      ],
      "text/plain": [
       "<Quantity 86. GHz>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/latex": [
       "$[5000.0034,~5000.0034,~5000.0034,~5000.0034] \\; \\mathrm{M}$"
      ],
      "text/plain": [
       "<Quantity [5000.00336439, 5000.00336439, 5000.00336439, 5000.00336439] M>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/latex": [
       "$59.646318 \\; \\mathrm{M}$"
      ],
      "text/plain": [
       "<Quantity 59.64631839 M>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/latex": [
       "$79.528425 \\; \\mathrm{M}$"
      ],
      "text/plain": [
       "<Quantity 79.52842452 M>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(meta.mass, meta.dist, meta.freq)\n",
    "display(meta.time, meta.width, meta.height)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "80e448ea-6613-4a22-a9b8-cd7a80a61a4d",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/latex": [
       "$299.99972 \\; \\mathrm{\\mu as}$"
      ],
      "text/plain": [
       "<Quantity 299.9997182 uas>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "text/latex": [
       "$399.99962 \\; \\mathrm{\\mu as}$"
      ],
      "text/plain": [
       "<Quantity 399.99962426 uas>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "display(meta.width.to(u.uas, meta.geom),\n",
    "        meta.height.to(u.uas, meta.geom))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0158cdb3-3ee7-4def-af32-1a142dcc60d0",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.9.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}