Skip to content

Instantly share code, notes, and snippets.

@wtbarnes
Created January 9, 2022 21:47
Show Gist options
  • Save wtbarnes/259d5966538af741ebc97d758033c667 to your computer and use it in GitHub Desktop.
Save wtbarnes/259d5966538af741ebc97d758033c667 to your computer and use it in GitHub Desktop.
Notebooks for testing synthesizAR workflows
Display the source blob
Display the rendered blob
Raw
{
"cells": [
{
"cell_type": "code",
"execution_count": 30,
"metadata": {},
"outputs": [],
"source": [
"import numpy as np\n",
"import matplotlib.pyplot as plt\n",
"import astropy.units as u\n",
"import astropy.time\n",
"import astropy.constants as const\n",
"from astropy.coordinates import SkyCoord\n",
"from astropy.visualization import ImageNormalize,SqrtStretch,LogStretch\n",
"from sunpy.coordinates import HeliographicStonyhurst\n",
"import sunpy.map\n",
"import distributed\n",
"import aiapy.response\n",
"\n",
"import synthesizAR\n",
"from synthesizAR.extrapolate import semi_circular_bundle\n",
"from synthesizAR.interfaces import MartensInterface\n",
"from synthesizAR.visualize import plot_fieldlines\n",
"from synthesizAR.instruments import InstrumentSDOAIA"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [],
"source": [
"obs_bundle = SkyCoord(lon=0*u.deg,\n",
" lat=0*u.deg,\n",
" frame=HeliographicStonyhurst,\n",
" obstime=astropy.time.Time('2020-01-01 00:00:00', scale='utc', format='iso'))"
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [],
"source": [
"bundle = semi_circular_bundle(150*u.Mm, 5*u.Mm, 500, observer=obs_bundle, n_points=1000, gamma=0*u.deg)"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [],
"source": [
"obs_view = SkyCoord(lon=0*u.deg,\n",
" lat=0*u.deg,\n",
" radius=const.au,\n",
" frame=HeliographicStonyhurst,\n",
" obstime=obs_bundle.obstime)"
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAOcAAAJOCAYAAABMeBAvAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAABAmklEQVR4nO3dd5wV1f3/8dcbliK9CAsCwgqIUaMiWIgdSywxqDF2JYkRjDWmWVI0muL3l2iixh5jiLFGLMRuRCI2pAioAQVZ6U1AyoLAsp/fH3MWL+v2vTN37vJ5Ph73cWfmzp1zbnnfKffMGZkZzrn0aZLrCjjnKufhdC6lPJzOpZSH07mU8nA6l1IeTudSqiDXFcgnnTt3LisqKlKu65HPiouLbcWKFb5SqAUPZx0UFRVp0qRJua5GXhs8eLD/uNWS/4I5l1KpDqekv0laJun9jGmdJL0saVa47ximS9KtkmZLmi5p34znDA/zz5I0PGP6J5n3zqVJqsMJ/B04tsK0q4BXzKw/8EoYBzgO6B9uI4A7IQozcC1wALA/cG15oJ1Ls1SH08xeA1ZWmDwMGBWGRwEnZUz/h0XeBjpI6g58HXjZzFaa2SrgZb4I/PIK986lRqrDWYVCM1schpcAhWG4BzA/Y74FYVpV0zGz/TLvnUuTfAznVhadUuOn1bhGKR/DuTRsrhLul4XpC4FeGfP1DNOqml4rkh6QVCKppLi4uEEVd64u8jGcY4DyI67Dgaczpp8XjtoeCKwOm78vAsdI6hgOBB0TptWKmZ1rZq3NrHVRUVH2XoVzNUh1IwRJDwOHAztKWkB01PVG4DFJ5wNzgdPC7M8BxwOzgfXAdwHMbKWkG4CJYb7rzaziQSbnUifV4TSzM6t46MhK5jXg4iqW8zfgb1msmnOxy8fNWue2Cx5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lEr1ydbbi7KyMtauXcvq1atZu3Ytq1atYvHixSxbtoxPP/2UVatWsXz5clasWMGaNWsoKSlh/fr1lJSUsHHjRjZv3kxZWRktW7Zk48aNNGnShIKCAlq0aEHr1q1p1aoVrVq1on379nTu3JmuXbvSoUMHdtxxR7p27Uq3bt3o2LEjbdu2pX379rRr144mTfx3O9c8nDH7/PPPWbp0KUuWLGHZsmVbA7do0SIWL17M2rVrWb9+PVu2bGHz5s1I2norn1ZQUEDHjh0pLCxkp512YvHixXz++efstttuDBgwgJ49e9K+fXvMjAkTJvCVr3yFBQsW8OGHHzJz5kyaN29O9+7dMTOWLl3KRx99RGlpKc2bN6dJkyaYGZIoKyujoKCAgoICWrVqRZs2bdhpp53o3r07Xbp0oUuXLhQWFlJYWEi3bt1o2bJlrt/eRs3DmQWbN29mwYIFzJs3j/nz5zN//nwWLVrE559/TvPmzWnTpg1lZWVs2LCB1atXs3HjRjp16sQBBxzAzjvvTNu2bVmyZAkff/wxixYtYuedd2bPPfdkt912o2/fvrRs2ZIPPviAe+65h8LCQn72s58xYMAApC9fE2j06NH88pe/3DpuZsyePZvRo0ezaNEifvzjH7PXXnuxceNGPv74Y2bOnMn777/P3Llz2WmnnSgqKmKnnXZi3bp1zJs3j1mzZjFlyhSaN29O+/bt2WGHHWjSpAnr169n48aNtGjRgh49etCzZ0969erFzjvvTK9evWjWrFmSH0GjpKjrHVcbgwcPtmeeeYbZs2czZ84cPv74YxYvXkzTpk3p2bMnO++8M927d2fjxo0sW7aMWbNmsXTpUrp168aAAQMYMGAAffv2pXPnznz66ae8+uqrjB8/nmbNmnHggQey33770adPn21CV1JSwu23386qVau46KKL6NWrVzU1hJEjR3L33XdX+tiCBQu48847adOmDZdccglt27bd+piZMW/ePN555x3eeustNm3axMEHH8zQoUPp0qULK1euZPbs2cyaNYsPP/yQRYsW0aVLFwYMGEBhYSHNmzdnyZIlW3+gtmzZQrdu3dhll13o27cv/fr1o1u3buy3335MmjSp2iuNSboC+D5Rn8TvEXXW1h14BOgMTAbONbNNktoADwFtgHPMbFFNn2O+8HDWwY477mgjR46kX79+9O3bl1122YVu3boxe/ZsJk+ezJQpU1i7di1f+cpX2Guvvdhzzz3p2rXr1rBt2rSJ8ePH88ILL1BaWsrQoUM55JBD6NChQ6XlzZw5k5tuuokzzjiDoUOHVrqmrKi6cJYbN24c//znP7niiivYY489Kp1n9erVjB8/nrFjxyKJY489lkMPPZQWLVoAUZiXL1/Oe++9x3vvvceMGTNo3bo1AwcOZPDgwey6664sWbJk64/YrFmzWLJkCdOmTas2nJJ6AK8Du5vZBkmP8UXPik+Y2SOS7gKmmdmdki4k6tF/HnC2mV1V1bLzjYezDgYPHmyTJk1i7dq1W9cwxcXF9O/fn8GDB7PvvvvSqVOnLz1v+fLlPPXUU7z99tt87Wtf4/jjj6d79+7VlvXSSy/x9NNP84tf/KLGeTPVJpwAy5Yt4/rrr+f444/n+OOPr3beJUuW8Pzzz/P666+z//77c/LJJ9O1a9cvzffZZ58xefJkJk+ezIcffkifPn048MADOfDAA7eupQcPHlybcL4N7A2sAZ4CbgMeBLqZWamkIcB1ZvZ1SRcDc4jCeZ6ZXVnji88THs466Nevn5100kl8/vnn7L///gwZMoR+/fpVuUabO3cuDz30EAsXLuTkk0/msMMOo6Cg+t18M+Ohhx5ixowZ/OIXv6jzQZfahhNg48aN/P73v6eoqIjzzjuvxjVzaWkpr7/+Ok888QRdu3blrLPOYpdddqnydcyZM4e33nqLt99+m+bNm3PwwQfzu9/9rjabtZcDvwU2AC8BlwNvm1m/8Hgv4Hkz21NSB+BhoCXRpu6CWr34PODhrINOnTrZ+eefT6tWrYDoC5ipfHzVqlVMnz6dTZs2sccee1BYWLjNF7+q55kZ06dPZ/369ey///71+jvjhRde4NhjK141sWpmxsSJE2nevDn77LPP1npWDGrF+i9fvpz333+fgoIC9tprr61bDFU9b8OGDcybN4+XXnqJVatWLQDeAh4xsycqzN8RGA2cDnwG/At4nGhN+aVw1vqF5iE/WlsHvXr14oorrtjmC5j5ZV6+fDmPPPIIpaWlXH/99ey+++6VzltxuNyjjz5KixYtuOyyyyp9vDb7nEuWLOHKK6vesqvsx9jMuOOOO2jXrh3nnHNOtc+pODxz5kweeeQRtmzZwplnnknXrl23+bGpuJyvfvWrAF2AU4BT+XJDmKOAYjNbDiDpCeAgoks6FphZKXW83k3eMjO/1fI2aNAgq8yGDRvsvvvus5EjR9rbb79tZWVllc5XnTFjxti1115rW7ZsqfNzM40YMaJezysrK7MbbrjBRo8eXa/nvvPOO3bhhRfavffeayUlJVXOG97DKt9jooscfwC0AkR0DdZLidagZ4R57gIuqm45jeHmzUAawMx47bXXuOSSS+jUqRO33347BxxwQK3WcJkmTpzI2LFjueaaa3LWMkcSV111FW+88QZvvfVWnZ+73377cfvtt1NYWMill17Kq6++Wh62OjGzCUSbsVOI/kZpAtwDXAn8SNJsor9T7qvzwvNNrn8d8umWueZcunSpXXPNNXbjjTfaqlWrrL7mz59v3//+9xu0jEz1XXOWW7NmjV1wwQX2ySef1HsZq1evtj/84Q925ZVX2uLFi7d5rKY1p9++uOW8Avl0GzRokJWVldm///1vu+CCC2zKlCnWEBs2bLALL7zQZs2a1aDlZGpoOM3M5syZYyNGjKh287Q2pk2bZhdccIE9+eSTWzf1PZy1v/lmbR2UlpZy9dVXU1xczG233cbAgQMbtLxbb72VE088kX79+mWphtlRVFTEqaeeyp/+9KcGLWevvfbitttuY9GiRfzsZz9j+fLlWarh9sHDWQcLFy7ktNNO49JLL93aUqa+xo0bx8aNGznuuOOyVLvsOvroo2nSpAkvv/xyg5bTokULLrroIs455xyuvvrqLNVuO5HrVXc+3QYOHGjZ8Omnn9r5559va9euzcryMmVjs7ZcSUmJff/737elS5dmZXlr1qyxjh07GnAT4T92v1V98zVnHWTjSKqZcfPNN3PhhRfSpk2bLNQqPq1ateKSSy7hpptuwqzhjVXatm1Ljx49AL4JLJTUsP2CRs7DmbCxY8fSsWNHBg8enOuq1Mree+9N9+7deemll7KyvBYtWmBm/YFbgTOystBGysOZoLVr1/LQQw8xcuTIXFelTkaMGMFjjz3G6tWrs7ZMM7vRGlEj9Th4OBN03333cdZZZ21zHmU+aNWqFcOHD691g3qXHR7OhMyaNYt58+YxdOjQXFelXg455BA+/fRTZsyYkeuqbDc8nAkwixqWX3LJJXVu2pcWkrj44ou58847s3JwyNXMw5mA8ePH07NnzyrPfcwXvXv3pn///owdOzbXVdkueDhjVlpayj//+U++973v5boqWXHeeefx8MMPs3nz5lxXpdHzcMbs2Wef5aCDDqJjx465rkpWtG/fniOOOIKnn346luVLGiBpasZtjaQfSuok6WVJs8J9xzB/E0n/kPSmpMo7RMpTHs4Ybdy4kTFjxvDtb38711XJqlNPPZXnn3+eDRs2ZH3ZZvahme1jZvsAg4D1wJPAVcAr4T/SV8I4wDHABOBk4MdZr1AOeThj9NRTT/H1r399a7cmjUWLFi044YQTeOKJJ2qeuWGOBD42s7nAMKITrwn3J4XhpkBZuOXn0bYqeDhj8vnnn/PSSy8xbNiwXFclFieeeCJjx45l/fr1cRZzBlHnXQCFZrY4DC8BCsPwi8BhwBjg5jgrkzTvQ6gOSktLWbZs2TbTKv6tUD7+9NNPc+CBB7Jq1aqs/fVQm+WUlJSwaFHV/Spn668cSRx88MHcf//9fOtb36p02Q0pS1Jzoja4XzqVxcxMkoXhUhppM0APZx0sXbqU2267bZtplX0BS0tLef755znuuOO44447qvyS1vbLW5cv+cyZM7n33ntrnK+2PxhVzWdmbNmyheeee46FCxdW2uVnZT9cxcXFSCrJmNa6iqKPA6aY2dIwvlRSdzNbLKk7sKyK5zUaHs466NGjBzfccEON8z399NPssccenH322QnUaluLFi3i2muvTay8vfbai6ZNm3LqqafWav4XXniBFStWVBXITGfyxSYtRJutw4Ebw308h4tTxPc5s2zLli38+9//brT7mhWdeOKJPP/885SWlmZtmZJaA0cDmUecbgSOljSLqPvMG7NWYEp5OLPstddeY9CgQak/VzNbWrVqxZAhQ7LaasjMSsyss5mtzpi2wsyONLP+ZnaUma3MWoEp5eHMIjNj9OjRWw+QbC9OOeUUnnrqKW9zm2Ueziz64IMP6NatW6UX+WnMOnXqRO/evZk6dWquq9KoeDiz6PHHH6/1gZHG5tRTT+Xxxx/PdTUaFQ9nlixbtoyVK1ey22675boqOdG3b182bNjA4sWLa57Z1YqHM0vGjBmz3RyhrcqwYcNiaxC/PfJwZsHmzZt56623OOyww3JdlZw66KCDmDhxIhs3bsx1VRoFD2cWvPbaawwZMqTGC+M2dgUFBRx66KGMGzcu11VpFDycWfDss8/yjW98I9fVSIUTTjiB5557LtfVaBQ8nA00f/58mjVrRrdu3XJdlVTYcccdadOmDcXFxbmuSt7zcDbQs88+ywknnJDraqTKN77xDZ555plcVyPveTgboLS0lIkTJ/K1r30t11VJlf3335+pU6d6P0MN5OFsgAkTJjB48ODt/kBQRU2bNmXIkCG88cYbua5KXvNwNsBzzz2X2kv45dqxxx7LCy+8kOtq5DUPZz199tlnrF69mj59+uS6KqnUs2dPNmzYwIoVK3Jdlbzl4aynsWPH5u2lFZJy1FFH8Z///CfX1chbHs568nDW7PDDD/cGCQ3g4ayH+fPn065dOzp06JDrqqRa27Zt6dKlC3PmzMl1VfKSh7Me/vOf/3DkkUfmuhp54eijj/ZN23rycNaRmfHmm29y0EEH5boqeeGAAw5gwoQJ3ktCPXg46+ijjz6id+/etGzZMtdVyQvNmzdn11135YMPPsh1VfKOh7OO/EBQ3Q0dOpRXX30119XIOx7OOpo0aRL77bdfrquRVwYOHMjUqVMpKyvLdVXyioezDjZs2ED//v1p1qxZrquSVwoKCth999157733cl2VvOLhrIN169Zx+OGH57oaeemwww7z/zzryMNZByUlJQwaNCjX1chLAwcOZPr06bmuRl7xcNZBixYtfJO2npo2bcqAAQNyXY284uGsg7Zt2+a6Cnnt0EMPzXUV8oqHsw5at67NxbFcVQYNGsScOXOQVO3mh6QOkh6XNFPSDElDJHWS9LKkWeG+Y5i3iaR/SHpT0h7JvJJkeDjrIFsXnt1eNWvWjA0bNgD8oIZZbwFeMLPdgL2BGcBVwCtm1h94JYwDHANMAE4GfhxHvXPFw+kS1blzZ4Dzq3pcUnvgUOA+ADPbZGafAcOAUWG2UcBJYbgpUBZujerX08PpElVYWAiwYzWzFAHLgfslvSvpr+F6nYVmVn6thyVAYRh+ETiM6OK6N8dT69zwzm/qYOnSpdx8c+Wff5wNu+uy7GnTpvHHP/6xxvni3ESvatmSkISZ9ajm6QXAvsClZjZB0i18sQkLgJmZJAvDpcAZ2al5ung466B9+/Ycf/zxlT5W05e9rmGob3gmTZpU52u21PWHpab5q3v8N7/5DZJKMuateJRtAbDAzCaE8ceJwrlUUnczWyypO7CsTpXOQx7OOmjZsmXqryLWvn17+vfvn+tqVKmoqIgVK1ZUedjbzJZImi9pgJl9CBwJ/C/chhNdbn440OivmOThdGl0KfCgpObAHOC7RMdHHpN0PjAXOC2H9UuEh9OljplNBQZX8tB21f2EH611LqU8nM6llIfTuZTycDqXUh5O51LKw+lcSnk4nUupvA2npCskfSDpfUkPS2opqUjSBEmzJT0a/sRGUoswPjs83idMP1zS3yV9R9J1uXw9zlWUl+GU1AO4DBhsZnsSnTZ0BvB/wJ/MrB+wii9OTTofWBWm/ynM51yq5WU4gwJgB0kFQCtgMTCUqKE0bHvOX+a5gI8DRypqWb4JWA1sANYlU23naicvm++Z2UJJfwTmEQXrJWAy8Fk4hQiisxvKT03qAcwPzy2VtBrobGZvAm8mWnnnaikv15yh/5hhRCfm7gS0Bo7NaaWcy7K8DCdwFFBsZsvNbDPwBHAQ0CFs5gL0BBaG4YVAL4DweHugVtdDl/SApBJJJcXFxdl8Dc5VK1/DOQ84UFKrsO9Yfs7fq8CpYZ7Mc/7GhHHC42OtlmcYm9m5ZtbazFoXFRVl7QU4V5O8DGc4S/5xYArwHtHruAe4EviRpNlAZ0InUeG+c5j+Iyp0e+FcGuXlASEAM7sWuLbC5DnA/pXM+znw7STq5Vy25OWa07ntgYfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0qSPpE0nvSZoqaVKY1knSy5JmhfuOYXoTSf+Q9KakPXJb8+zycLq0OsLM9jGzwWH8KuAVM+sPvMIXvfYfA0wATgZ+nHw14+PhdPki8xqrmddebQqUhZuSr1Z8PJwujQx4SdJkSSPCtEIzWxyGlwCFYfhF4DCii1XdnGw145W310rJhU2bNjF37txtptXyYmX1nr+u1qxZQ1yXKowu6Bbf/BkODhdI7gq8LGlm5oNmZpIsDJcCZ9S3oDTzcNbBihUrePDBB7eON+DLF5uPP/6YRx55JNfV2Krij1FxcTGSSjIeb13JcxaG+2WSniS6ONVSSd3NbLGk7sCyeGueex7OOujevTvXXHNNrqtRrU8++YSrr74619Wo0hNPPMGKFSu+FMhykloDTcxsbRg+BrieL66xeiPbXnu10fJwurQpBJ4MWyUFwENm9oKkicBjks4H5gKn5bCOifBwulQxsznA3pVMX0F0BfPthh+tdS6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0LqU8nM6llIfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lYuvxXdKYWsy20sy+E1cdnMtncV6O4SvA96t5XMDtMZbvXF6LM5w/N7P/VjeDpF/HWL5zeS22fU4zeywb8zi3vYr9gJCklyV1yBjvKOnFuMt1Lt8lcbR2RzP7rHzEzFYBXRMo17m8lkQ4yyTtXD4iqTdg1czvnCOZi+f+HHhd0n+JjtAeAoxIoFzn8lrs4QyXDN8XODBM+qGZfRp3uc7luyQOCAk4FtjXzJ4BWknaP+5ynct3Sexz3gEMAc4M42vxxgfO1SiJfc4DzGxfSe9CdLRWUvMEynUuryWx5twsqSnhCK2kLkBZAuU6l9eSCOetwJNAV0m/BV4HfpdAuc7ltSSO1j4oaTJwJNFfKSeZ2Yy4y3Uu3yVxtLYvUGxmtwPvA0dnNudzzlUuic3a0cAWSf2Au4FewEMJlOtcXkuk+Z6ZlQKnAH8xs58C3RMo17m8ltTR2jOB84BnwrRmCZTrXF5LIpzfJWqE8FszK5ZUBDyQQLnO5bUkGiHMBa4wsy1hfB5wWwLlOpfXklhzvgLskDG+A/CfBMp1Lq8lEc6WZraufCQMt0qgXOfyWhLhLAmnjAEgaRCwIYFynctrSexzXg78S9IiohZC3YDTEyjXubwWazhDg/dDgN2AAWHyh2a2Oc5ynWsMYt2sDUdozzSzzWb2frh5MF2NJDWV9K6kZ8J4kaQJkmZLerT8tENJbSSNkTRW0k65rXV2JbHP+Yakv0g6RNK+5bcEynX57XIg8wSJ/wP+ZGb9gFXA+WH6OUTNQi8HLku0hjFLYp9zn3B/fcY0A4YmULbLQ5J6AicAvwV+FLq6GQqcFWYZBVwH3Ak0JTo/uIzomEajkcQpY0fEXYZrdP4M/AxoG8Y7A5+FNtoAC4AeYfhB4GGgJXBugnWMXRJrTiSdAOxB9AYCYGbXV/2MdCopKeHtt98GwKzhXe9mYxkVLVmyhDfffDNry4tWWsktQ9I3gGVmNlnS4TXNHzosP66+dUuz2MMp6S6iRgdHAH8FTgXeibvcOGzYsIEZMxp2nng2vuzVWb16NR999FEsy87Gj0lxcTGSSjKW2brCLAcB35R0PNGPeTvgFqCDpIKw9uwJLGxwZdLOzGK9AdMr3LcBxsddbhy3QYMGWdqNGDEi11WoVngPa/vdORx4Jgz/CzgjDN8FXFTb5eTrLYmjteWtgdaHQ92b8fM5Xd1dSXRwaDbRPuh9Oa5P7JLY53wmdEvyB2AK0ZHaexMo1+U5MxsHjAvDc4DtqjPyJI7W3hAGR4c/lFua2eq4y3Uu38W2WVtZQwMz25gZzIY0RpDUQdLjkmZKmiFpiKRO4Xqgs8J9xzCvJN0aWpdMLy9XUh9J4yQdLunv9a2Lc3GIc5/z/nCh3E5V3WjYfsMtwAtmthuwN1FrkquAV8ysP9F5pFeFeY8D+ofbCKI/r51LtTg3a9sDk6m+1cby+ixYUnvgUOA7AGa2CdgkaRjRET6IWpGMIzqQMAz4h0WH+t4Oa93uwBZgJbAJ8E1tlyqxhdPM+sS1bKCIKNj3S9qb6EfgcqDQzBaHeZYAhWG4BzA/4/kLgB5mNomoV0CA7P1z71wWJPFXShwKgH2BO81sIFDCF5uwAIS1pF9B2+WtfA3nAmCBmU0I448ThXVp2Fwl3C8Ljy8k6sy6XK1bmEh6QFKJpJLi4uKsVN652sjLcJrZEmC+pPITuI8E/geMAYaHacOBp8PwGOC8cNT2QGB1xuZvTWWda2atzax1UVFR9l6EczVIom1tE6KjqTsRtRZ638yWVf+sWrkUeDCcdDuHqH/cJsBjks4n6pLztDDvc8DxwGxgfZjXuVSLLZzhAkZXAkcBs4gO4LQEdpW0nugE2VFmVq9rdZrZVGBwJQ8dWcm8Blxcn3Kcy5U415y/Ifo/cWQIx1aSuhKdOHsu0V8ezrkK4gznSDNbU9kDYbP2zzGW7Vzei/OA0LuSzohx+c41anGGcyhwemjj2i/GcpxrlOJsITQXOFnScUQ98E0k6oSp/PFvxlW2c41B3J1KDwB+AowHbicjnM656sX5V8qNRA3OrzCzF+Iqx7nGKs59zlJgYGYwy3vvds7VLM59zl9UMrlHJdOcc5WItW2tpCaSTsuY9G6c5TnXmMR9IaMyop67y8e/F2d5zjUmSZyV8h9JP5HUq0IXJc65aiTRNWb5hXIzG54bsEsCZTuXt5LoGtNPgnSuHpK6kNGewO5seyGjfyRRtnP5KomTra8l6hFvd6KTno8DXgc8nM5VI4kDQqcSnQC9xMy+S9QrQvsEynUuryVyIaPwl0qppHZEnW71quE5zm33ktjnnBQuZHQvUf+y64C3EijXubyWxNHai8LgXZJeANqZ2fS4y3Uu38V5VkqVFymStK+ZTYmrbOcagzjXnDdV85gR9ZTgnKtCnGelHBHXsp3bHiTxP2cz4AdEVwWD6Mpfd5vZ5rjLdq6isLu1zsw+CuOtgWOAxWb2dk4rV0ESR2vvBJoBd4Txc8O07ydQtnMV3Ua4ZIckETWImQfsKOk5M/ttLiuXKYlw7mdme2eMj5U0LYFynatMezObHYYPAXYws2GSCoCpQGrCmUQjhC3h0gwASNqF6KK1zuVCacbwUKIroGNmpaTse5nEmvOnwKuS5hBd5bo3fiEhlzuzJd1AtJYcAZwBEM4x3r7CaWavSOoPlF+u70Mz2xh3uc5VYQRwHXAecJWZvRamNwdG5qpSlYmzEcJ+wHwzW2JmGyXtA3wLmCvpOjNbGVfZzlXFzFZK+ilhZSGphZltDNd8XZLb2m0rzn3Ou4FNAJIOBW4kOk1sNXBPjOU6VyVJRwAfEh21vRWYJelLl41Mgzg3a5tmrB1PB+4xs9HAaElTYyzXuercBAw1szmw9Tqy/wKqbG6aK3GuOZuGw9MQnc85NuOxRHpgcK4STcuDCWBmHwNNc1ifKsUZkoeB/0r6lOhy8+MBwhXHVsdYrnPVmSjpfuCBMD4cmJjD+lQpzra1v5X0CtAdeCnj6tZNgEvjKte5GvyA6Khs+amMrxG1WEudOI/WtqmsrWJ5m8aMedbFVQfnKgptuv8SbqkW52bt0+HAz9PAZDMrga0thI4ATiPqHeHxGOvg3DYkjSVqDFOpNJ1NFedm7ZGSjifahDhIUkeiplMfAs8Cw8N/S84l6Se5rkBtxXrU1MyeI+oO07lUMLMpkvYC1pjZJ7muT3WSaPjuXGqEI7WjgJclXSipg6Q7anpeLng43fZmCFGDg0HABWb2GTA4pzWqgofTbW/mA13NbA1QIKkJsEOO61SpRMIp6WBJ3w3DXST5xY1crqwGpoXN225E53Om8h+DpK6VMpjoLID7ibos+SdwUNxlO1eJMeEG8DzwPzN7P4f1qVISbVxPBgYCUwDMbJGktgmU69yX5NPV7ZII5yYzM0kGW3s7cy4nMnrk2EYaryObRDgfk3Q30EHSBcD3iFoGOZcLmUdmWxOdzrhjjupSrSS6KfmjpKOBNUT7nb8ys5fjLte5ylTogWMl8EdJk3JVn+okcl5lCKMH0qXV/5PUJFyqMjXiPCtlLdE1USplZu3iKtu5ujCzx3Jdh8rE2fC9LUDohnAx0cmtAs4mOsfTOVeNJBohfNPM7jCztWa2xszuBIYlUK7LQ5JaSnpH0jRJH0j6dZheJGmCpNmSHpXUPExvI2mMpLGSdspt7bMriX3OEklnA48QbeaeCZQkUK7LTxuJOuBaFy6C9bqk54EfAX8ys0ck3QWcT9SDwTlEPT3OAy4Drqpu4aFRTJXM7NdZeA1ZkcSa8yyiE6uXhtu3wzTnvsQi5b1jNAu38uu5ljezGwWcFIabAmXhVuVJ1BkuB9YBa6u4pUYSf6V8gm/GujqQ1BSYDPQDbgc+Bj4L1zMBWAD0CMMPEnUm15LoCnY1WWdm1V3YOTWSaFt7P5UctTWz78VddrbNnTuXESNGbB2PriAX+aL/stwaP378NnXMhYa+L2a2BdhHUgfgSWC3aub9DDiuLtWrc4VyJIl9zmcyhlsStbVdlEC5Wde7d2/uuSfdndWPHDmSu+++O9fVqFLnzp2RtPWYg5lV2ZzTzD6T9CrROZgdJBWEtWdPYGE9q5COX9FaSGKzdnTmuKSHiS5Y6rZDRUVFrFixospASuoCbA7B3AE4Gvg/4FXgVKIDi8OJOo6rj76SvgmsNrP/1nMZichFz+v9ga45KNflh+7AqLDf2QR4zMyekfQ/4BFJvwHeBe6r5/KfBjYDHUNPfH8G7jezkxte9exKYp+zYkuhJcCVcZfr8pOZTSc6xbDi9DnA/lkooreZ7SGpBTDBzK6T1DMLy826JDZr/dxNlyYfStrNzGZKQlJLomMhqRP7/5zhkgw1TnMuIR2Ad8OBpt5E10m5Jac1qkKcDd9bAq2AHUOH0uWHsNvxxX9UziXtuozhz4FZZrYqR3WpVpybtSOBHwI7EbooCdaQB9epcI1TxmXmt5L0ezO7Ohf1qU5sm7Vmdkvo+uEnZlaUcdvbzDycLick/V7SCklrJa2RtAb4aRhPVUDj3KwdamZjgYWSTqn4uJk9EVfZzlVjGFCY0RQQSVPMLHVXto5zs/YwoqtZn1jJYwZ4OF0uTM4MZvBBTmpSgzhPti4/Ned6MyvOfMw7lXa5YmZfahxf2bQ0SOKUsdGVTEtlD9uu8Qv7mVv3N8vHw2OpalYa5z7nbsAeQPsK+5ztSOmfvq7xq67vKjM7OMm61CTOfc4BwDeI/vTN3O9cC1wQY7nOVUlSAbAXkNly7UbgaqDYzObmpGKViHOf82miS88PMbO34irHuTp6kWh3LrPXg92AHwMPAY0/nBnelXQx0Sbu1s3ZfDzZ2jUKncxsm4b14a+Uyv5VyKkkDgg9QHSpta8D/yU6UTZVfbW47crfK5k2KulK1EYSa85+ZvZtScPMbJSkh4DxCZTr3JeY2S2SDgL68MX3/7OcVagaSYRzc7j/TNKeROdz+snWLick/ZMomFOJeuyD6KSM1K09kwjnPeGslF8QXbS0DfDLBMp1rjL7AntYWnpkq0bs+5xm9lczW2Vmr5nZLmbWFfg07nKdq8IMouMeqZeLPoQA/kTlLYeci1tH4H1J7xCdzwlAGo/W5iqcedN3qGt0rst1BWorV+FM/fa+a5zM7DVJ3YADwqQJZrYkl3WqSpxta9+j8hAKKIyrXOeqI2k40drzNaLv562SrjOz+3NasUrEueb8RozLdq6+rgH2Le83SFIn4C1g+wlnmhoQO5dhFdFVxsqtCdNSJ1f7nM7lylTgFUn/CuOnA9PD5i5mlprGCB5Ot73ZRBTQ/mG8vGfIQaSspVAi4QwXpNnZzD5MojznqmJml+W6DrWVRI/vJxL9Ur0QxveRNCbucp2rjKTOkv4paZmkpZIelLRjrutVmSROGbuO6AI0nwGY2VTAO/hyuXI70cqiB9E1Pm8H7splhaqSRDg3m9nqCtO8EYLLla+Y2R/NbDMgM3sT6JXrSlUmiX3ODySdBTSV1B+4DHgzgXKdq0zTzBFJqQwmJLPmvJSoi5KNRH20rCa6hopzufCapL3DcGeiPoWuyGF9qpTEmnM3M/s58PNsLzhc/XgSsNDMvhE6q36E6E2fDJxrZpvChVL/QXS4fAVwupl9Iulw4DvAOKCPmV2X7Tq6dDGzizJGdzOz9TmrTA2SWHPeJGmGpBtCTwjZdDnR+Xnl/g/4k5n1I2r1cX6Yfj6wKkz/U5jPbefSHExI5mTrI4AjgOXA3ZLek/SLhi43XCr8BOCvYVzAUL7oTX4UcFIYHsYXfy4/DhwZ5t9EtJm9gW2bdDmXc4k0Qgin5Nwarib8M+BXwG8auNg/h2WVdw7cGfgs4yI1C/jiIr09gPmhLqWSVgOdw5E6PzjlUimJRghfkXRdOIXsNqIwNKibCEnfAJaZ2eRs1NG5NEpizfk34FHg62a2KEvLPAj4pqTjiTqqbgfcAnSQVBDWnj2J/mQm3PcCFoTu+NsTHRiqkaQHgFMAOnXqlKXqO1ezJPY5h5jZn7MYTMzsajPraWZ9gDOAsWZ2NvAqcGqYbTjwdBgeE8YJj4+tbe9rZnaumbU2s9ZFRd6wySUnzp4QHjOz0yrpEUGAmdleMRR7JfCIpN8A7wL3hen3AQ9Img2sJAq0c6kW52bt5eE+1h4RzGwc0f+UmNkcona8Fef5HPh2nPVwLtti26w1s8Vh8CIzm5t5Ay6q7rnOuWQaIRxdybTjEijXubwW5z7nD4jWkH0lTc94qC3+36JzNYpzn/Mh4Hng98BVGdPXmtnKGMt1rlGIc59ztZl9QvT/48qM/c1SSQdU/2znXBL7nHeybbvVdWGac64aSYRTmX/4m1kZ3uufczVKIpxzJF0mqVm4XQ7MSaBc5/JaEuG8EPgaUfvWBUQXkBmRQLnO5bXYNy/NbBneXM65OkvilLFdJb0i6f0wvlc2TrZ2rrFLYrP2XuBqYDOAmU3H16TO1SiJcLYys3cqTCutdE7n3FZJhPNTSX0Jp41JOhVYXP1TnHNJ/N94MXAPsJukhUAxcE4C5TqX15I4WjsHOEpSa6CJma2Nu0znGoM4z0o5x8z+KelHFaZDtIm7EhhTfvlv59y24tznbB3u21Zya0fU+/rzMZbvXF6Lbc1pZneH+19XNY+k6+Mq37l8l0QjhJ6SngwXK10maXTorR0z+1Xc5TuXr5L4K+V+oq4pdwq3f4dpzrlqJBHOLmZ2v5mVhtvfgS4JlOtcXksinCsknSOpabidQy17W3due5ZEOL8HnAYsIWoZdCrw3QTKdS6vxdoIIVzc9ndm9s04y3GuMYp1zWlmW4DekprHWY5zjVESbWvnAG9IGgOUlE80s5sTKNu5vJXEPufHwDOhrMxWQs59iaRekl6V9D9JH4Q+p5DUSdLLkmaF+45hehNJ/5D0pqQ9clv77Eqi4fuvASS1i0a94burVinwYzObIqktMFnSy8B3gFfM7EZJVxF1VH4lcAwwAfgpUQfm38tNtbMviRZCg8NlAKcD70maJmlQ3OW6/GRmi81sShheC8wAegDDgFFhtlHASWG4KVAWbkq0sjFLYrP2b0RXGusTLnZ7Md5CyNWCpD7AQKI1Y2HGleuWAIVh+EXgMKJWaI3qOEYSB4S2mNn48hEze11SXnZTMnfuXEaOHPml6bW8SHYixo8fz4gRue95NJwauI26vE+S2gCjgR+a2ZrM5ZmZSbIwXEoj7ZMqiXD+V9LdwMNE53GeDoyTtC9A+SZMPujduzd33313rqtRrZEjR6a6jp07d0ZS5lH71hXnkdSMKJgPmtkTYfJSSd3NbLGk7sCyZGqcO0mEc+9wf22F6QOJwjo0gTq4lCgqKmLFihVfCmQ5RavI+4AZFf5uGwMMB24M90/HWtEUSOJo7RFxl+EalYOAc4kOHk4N064hCuVjks4H5hI1CW3UYg+npPZEa81Dw6T/Ateb2eq4y3b5x8xep+qjrkcmWZdcS+po7VqiX7rTgDX40VrnapTEPmdfM/tWxvivMzZXnHNVSGLNuUHSweUjkg4CNiRQrnN5LYk15w+AUWHfU0RdYn4ngXKdy2tJHK2dCuwd2tZiZmviLtO5xiBXnUr7KWPO1SDONWdmp9LOuTrKaafSzrmqxblZe2t1j5vZZXGV7VxjEOdm7eQYl+1coxfnZu2ozHFJrcxsfVzlOdfYJNETwhBJ/wNmhvG9Jd0Rd7nO5bskWgj9Gfg6oZd3M5vGF43gnXNVSCKcmNn8CpO2JFGuc/ksieZ78yV9DbBwhvvlRJ02OeeqkcSa80KiTr16AAuBfcK4c64aSbSt/RQ4O+5ynGts4myEcBtRH0GV8kYIzlUvzjXnpIzhX/PlDr6cc9VIpBGCpB9WbJTgnKteIn+lUM3mrXOuckmF0zlXR3EeEFrLF2vMVpLKe0AQUY/67eIq27nGIM59Tj/J2rkG8M1a51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0LqU8nM6llIfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1J5GU5JvSS9Kul/kj6QdHmY3knSy5JmhfuOYbok3SpptqTpkvYN0/tIGifpcEl/z+FLcu5L8jKcQCnwYzPbHTgQuFjS7sBVwCtm1h94JYwDHAf0D7cRwJ3JV9m5usnLcJrZYjObEobXAjOAHsAwYFSYbRRwUhgeBvzDIm8DHSR1B7YAK4FNwOrkXoFzNSvIdQUaSlIfYCAwASg0s8XhoSVAYRjuAczPeNoCoIeZTQJOCdPejL+2ztVeXq45y0lqA4wGfmhmazIfMzMDLCcVcy4L8jackpoRBfNBM3siTF4aNlcJ98vC9IVAr4yn9wzTalPOA5JKJJUUFxdnp/LO1UJehlOSgPuAGWZ2c8ZDY4DhYXg48HTG9PPCUdsDgdUZm7/VMrNzzay1mbUuKirK0itwrmb5us95EHAu8J6kqWHaNcCNwGOSzgfmAqeFx54DjgdmA+uB7yZaW+fqIS/DaWavA6ri4SMrmd+Ai2OtlHNZlpebtc5tDzyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0qSLpb5KWSXo/Y1pVJ9E3kfQPSW9K2iN3tY6Hh9Olzd+BYytMq+ok+mOIThU8GfhxUhVMiofTpYqZvUZ0Anymqk6ibwqUhVtVzTnzlofT5YOqTqJ/ETiM6Kyjmyt7Yj7Ly4bvubJo0SJ++ctf5roa1Zo4cWLq69gQZmaSLAyXAmfkuEqx8XDWQdeuXbn00ktzXY1qzZ07N9V1vOOOO5BUUj5uZq1r8bSlkrqb2eIKJ9E3ah7OOigoKKBr1665rka1dthhh1TXsaioiBUrVtQmkJnKT6K/kW1Pom/UfJ/TpYqkh4G3gAGSFoQT528EjpY0CzgqjDd6vuZ0qWJmZ1bx0JdOom/sfM3pXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0LqU8nM6llIfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0LqU8nM6llIfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0LqU8nM6llIfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncym13YRT0rGSPpQ0W9JVYdo4SX0kfZLj6rkMVXxWe0h6S9IoSdvF93a7eJGSmgK3A8cBuwNnSto9t7Vylanms/oR8E1gEnBM7mqYnIJcVyAh+wOzzWwOgKRHgGHASmALsLy+C35vwWoeeWcez723mA2bt1BaVpaVCtdXaZvj6HfNszkrv6BJE5o1bcI3996JM/bfma/2bF/XRVT1WTUFDCgDlM06p9X2Es4ewPyM8QXAAWZ2Shjfrz4LfXLKQm5/dTafrFhHaW4zuY1c1qW0rIzPS8t4ZOI8JhSv5OIj+nHyvj3qsohKPyvgFuBZYBZwV9YqnGLbxWZtXK7/9wfMXp6uYKbFFoPZy9dx/b8/yMryzOxdMzvAzM4xsy1ZWWjKbS9rzoVAr4zxnmFajSQ9AJwC0KJFC0aOHLn1sVVtjstiFbNo8+Zc12CrVZs3b/OeFRcXI6kkY5YnzOzcjPF6f1aNjcws13WInaQC4CPgSKIPeiJwlpnV6Wd98ODBNmnSpK3jRVc9S+revc2boVmzXNdiKwHFN56wdXzw4MFMmjSpyn3GbH1WjcF2sVlrZqXAJcCLwAzgsWx82F/fo9v2cWSiAb6+R7c6zR/XZ5WPtpfNWszsOeC5bC7zrnMHsevPn2PTltStP1OheVNx17mD6vy8OD6rfLTdhDMuH/32eG79zyxenrGUBSvX57o6rN28gbY7tMppHXp2asXRXynksqP657Qe+c7DmQWXHdU/NV/EkSNHcvdNd+e6Gi4Ltot9TufykYfTuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKeXhdC6lPJzOpZSH07mU8nA6l1IeTudSysPpXEp5OJ1LKQ+ncynl4XQupTyczqWUh9O5lPJwOpdSHk7nUsrD6VxKeTidSykPp3Mp5eF0LqW2iwsZZcuOO+5offr0qddzi4uLKSoqym6FclROQ8ooLi62FStW+EqhFjycCZFUYmatG0M5Sb2W7Z3/gjmXUh5O51LKw5mcJxpROUm9lu2a73M6l1K+5nQupTycWSTpcEmrJU0Nt19lPHaspA8lzZZ0Vcb0cZL6SPqknmV+abkNXWZYxieS3guvY1KY1knSy5JmhfuOYfp3JF0Xbt+pb5luWx7O7BtvZvuE2/UAkpoCtwPHAbsDZ0ravaEFxbXcDEeE1zE4jF8FvGJm/YFXwriLiYczGfsDs81sjpltAh4BhoXHVgJbgOVZXG5DllmdYcCoMDwKOCkMbwDWhduGLJe53fLLzmffEEnTgEXAT8zsA6AHMD9jngXAAQBmdkqYtl89yqp0uQ1cZjkDXpJkwN1mdg9QaGaLw+NLgEIAM3u0AeW4Kng4s2sK0NvM1kk6HngK6J/bKtXbwWa2UFJX4GVJMzMfNDMLwXUx8c3aBpJ0cfkBIKCNma0DMLPngGaSdgQWAr0yntYzTGuouJaLmS0M98uAJ4k2oZdK6g4Q7pdloyxXOQ9nA5nZ7eUHgIAySQKQtD/R+7sCmAj0l1QkqTlwBjAmC8XHslxJrSW1LR8GjgHeD8seHmYbDjzd0LJc1XyzNrtOBX4gqZTowMgZFrXyKJV0CfAi0BT4W9gXbRAzi2W5RPuST4bfmQLgITN7QdJE4DFJ5wNzgdOyUJargrcQci6lfLPWuZTycDqXUh5O51LKw+lcSnk4nUspD6dzKVVjOCWtqzD+HUl/qeE5W+eRdKGk8xpWzRrrWK8yJHWQdFHG+E6SHs9CffpIer+hy6mhjDczyjqrHs+v8XOs5DkDJd0XhgslPSNpmqT/SXquhuc26D2RdGuFU/B+Lun2MPxHSUOreN64cErdN+tbdjZJOiS8XzW/F2ZW7Q1YV2H8O8BfanhOjfPU9QYUZHN5YZl9gPfzZblVlHU48Ew9nlfnzwj4F7B3GL4buDzjsb2y+Z5U/LyBdsAcYJdwKwY6hMd6Ay9VsZxxwOA6vs6sf9fq8140aLNWUhdJoyVNDLeDKpnnOkk/CcP7SHpb0nRJT2acrDtO0i2hjer7oelb+XMfkPQG8ED49R0bnv+KpJ0rKaOvpBckTZY0XtJuYXphKHNauH0NuBHoG8r9Q+ave6jnHhmvY5ykwaFp298kvSPpXUnDqCVJR4bnvBeW0SJM/0TSryVNCY+V17mLopOaP5D0V0lzFbXVzdyiuRE4JLyGKyquEcPa7fAw/F1JH0l6BzgoY57afI5tiQI4LUzqTnQWDABmNj3Mp/Bevh9ey+mVLKtP+GymhNvXwvTDw/QxwP8yn2Nma4CfA38Jt1+Z2WfhsblAZ0ndavEZXBBe47TwmluF6X+XdJekCcD/k9RP0n/CfFPC96q7pNcyvqeHhOceI+mtMN+/JLUJ0/eT9GZYxjvhPay9WqR8CzA14zaP8IsLPER09gLAzsCMir/KwHVEp04BTAcOC8PXA3/O+HW7NwwfSvhVCc+dDOwQxv8NDA/D3wOeqqSMV4D+YfgAYGwYfhT4YRhuCrSnwi9Y5jhwBfDrMNwd+DAM/w44Jwx3AD4CWtf0ywi0JDq9a9cw/o+M+nwCXBqGLwL+Gob/Alwdho8lOo1rxzC+LtwfTsaakwprROCZME/38Nl1AZoDb9T0OVao/xHA6IzxrwOfAa8ShWanMP1bwMvhPS4MZXav8N62AlqG4f7ApIzXUgIUVfN9fAt4vZLp9wLfqmnNCXTOGP5Nxvv+9/BeNQ3jE4CTMz67VsCPgZ9nfIfaAjsCr5V/B4ArgV+F93gOsF+Y3o6wRq7s+1HZrTZtazdY1KgbiPZVgPIz448CdlfUBhOgXfmvRkWS2hNthvw3TBpFtJlU7mEAM3tNUjtJHcL0MWZWfgLvEKD8XMUHgP9XoYw2wNeAf2XUqUW4HwqcF8rYAqxWWHNX4THgJeBaojak5fuixwDfLF9TE31wOwMzqlkWwACg2Mw+CuOjgIuBP4fx8h7tJme8xoOBk0OdX5C0qoYyqnMAMM7MlgNIehTYNTxW6edo4QyboDsZJ2+b2YuSdiH60TgOeFfSnqHOD4f3eKmk/xKdVzo9Y1nNgL9I2ofox3/XjMfeMbPiyl6ApJ6hHmWV1G8ZsFMt3oc9Jf2G6Ie1DVG75HL/MrMtYQ3Xw8yeDK/181D+ROBvkpoRrRimSjqMqBeKN8L715zoB2QAsNjMJoZlrKlF3bbR0IbvTYADyytfLuNDrouKjXzLx0vqWJ/PMn9M6suicxlXSNoLOB24MDwkol/oDxtaRgUbw/0WGva5lLLtgb6WtXhOpZ9jBRsqLsvMVhKtdR+S9AzRVk9tXAEsBfYOZWeWWwIgqRfRlhLAXWZ2F3AL0Y/lV8L9TzOe15La9cLwd+AkM5sWVjSHVyy7KmHFcShwAvB3STcDq4CXzezMzHklfbUWdalWQ/9KeQm4tHwk/BJWysxWA6vKt9OBc4H/ZsxyeljGwcDqMH9FbxKdFgVwNjC+QhlrgGJJ3w7LkqS9w8OvAD8I05uGNflaok2TqjwK/Axob2GfiuiX9lJp66lhA6t5fqYPgT6S+oXxiq+/Mm8QzvyQdAxQ2Zq+4mv4BNhHUpPwBd8/TJ8AHCapc/jl/3bGc2rzOc4A+mXMMzRjf60t0JdoE3Y8cHp4j7sQBfadCstqT7RWKSN6H5pWLMzM5tsXfTHdJek4oCvR7sANwCnatr+kXYlOa6tJW2BxeA/OrmwGM1sLLJB0Unh9LSS1ktQbWGpm9wJ/BfYF3gYOKv9cFR2T2JXo8+4uab/y90hSnX50GxrOy4DBig7Q/I8v1i5VGQ78QdJ0YB+i/c5yn0t6F7gLOL+K518KfDc8/1zg8ozHyte0ZwPnK+oq5AO+6KvncuAISe8RbTrubmYriDZH3pf0h0rKe5zox+CxjGk3EG2WTZf0QRivzABJC8pvwInAd4k2ud8DysJrrc6vgWMUHaT6NlHXIGsrzDMd2BIOOlxBFOhiogMqtxL1zoBF3YtcR7TJ9QbbbobX+Dma2UygfcZBjUHApPBZvEW0nzyR6MTs6cA0YCzwMzNbUmFxdwDDw2e0GzWssSS1JNr8v8giJURrzfK/65oR/XBMqm45wS+JfqjeAGZWM9+5wGXh9b0JdCNay04L39PTgVvCbsJ3gIcz3ovdLOrT6XTgtvA6X6Z2WzFbpeKUMUnjiA7o1ObNrez5twFTzOz+rFYsxxQdzd1i0XmbQ4A7s7HJ3oD6XAGsNbO/5qoOlZF0MrCvmf2yksfG0YDvVhwk9SE6iLdndfPlfQshSTcQHezIRs8CabMzMDH88t4KXJDj+tzJF/vGaVIA3FTFYyuJ9g9T0wiBaF/60xrnTcOa0zn3ZXm/5nSusfJwOpdSHk7nUsrD6VxKeTidSykPp3Mp9f8BH5EDsWLWKBYAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 720x720 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig,ax = plot_fieldlines(*bundle, \n",
" observer=obs_view,\n",
" figsize=(10,10),\n",
" plot_kwargs={'color': 'C0','alpha':0.2,},\n",
" imshow_kwargs = {'extent':[-35, 35, 950, 1020]},\n",
" draw_grid=True,\n",
" check_visible=False,)"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {},
"outputs": [],
"source": [
"skeleton = synthesizAR.Skeleton([synthesizAR.Loop(f'loop{i}', b, np.nan*np.ones(b.shape)*u.G)\n",
" for i,b in enumerate(bundle)])"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [],
"source": [
"martens = MartensInterface(1e-3*u.erg/u.s/(u.cm**3))"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [],
"source": [
"client = distributed.Client()"
]
},
{
"cell_type": "code",
"execution_count": 9,
"metadata": {},
"outputs": [],
"source": [
"status = skeleton.load_loop_simulations(martens, 'martens.zarr')"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [],
"source": [
"aia = InstrumentSDOAIA([0,1]*u.s, obs_view, )"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [],
"source": [
"aia.observe(skeleton, 'aia', channels=[aiapy.response.Channel(171*u.angstrom)])"
]
},
{
"cell_type": "code",
"execution_count": 26,
"metadata": {},
"outputs": [],
"source": [
"m_aia = sunpy.map.Map('aia/m_171_t0.fits')"
]
},
{
"cell_type": "code",
"execution_count": 38,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAtAAAADcCAYAAABZCpbOAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAA+iElEQVR4nO3debxd09kH8N8vk8gkZkErqbGoqKk1e/G2qFar5il5UWpI06pWqSGiRb2qgiBBGqJVlJbXq0UpRfuWCKFaOkjUFGQQSUhkeN4/9t5nP/fctc7e5+aee09uft/P53zOOmvvvfY6073rrP2stWhmEBERERGRcrp1dgVERERERFYkakCLiIiIiNRBDWgRERERkTqoAS0iIiIiUgc1oEVERERE6qAGtIiIiIhIHdSAFhERERGpgxrQIiIiIiJ16BHbQPLeEsfPNrPh7VcdEREREZHmFm1AA/gkgBNrbCeAse1bHRERERGR5larAf19M3us1sEkL2zn+oiIiIiINDWaWWfXQURERERkhVE4iJDkQyQHuserk3ygobUSEREREWlSZWbhWMvM3ssemNkcAOs0rEYiIiIiIk2sTAN6GcmPZw9IbgRAcR8iIiIislKqNYgw830AT5B8DMnMG7sDOKmhtRIRERERaVKlBhGSXAvAZ9OH/2dmMxtaKxERERGRJlVmECEB7AdgOzO7D0Afkjs1vGYiIiIiIk2osAea5HUAlgHY28w+SXJ1AA+a2Y4dUUERERERkWZSJgb6M2a2HclngWQWDpK9GlwvEREREZGmVGYWjsUkuyOdeYPk2kh6pEVEREREVjplGtBXAfgVgHVI/hDAEwAubmitRERERESaVNlZOLYAsA+SaeweNrO/NbpiIiIiIiLNqMwgwo0BvG5mi0juBWAbALf41QlFRERERFYWZUI47gKwlOQmAMYB+BiAnze0ViIiIiIiTarUUt5mtgTAwQCuMbPvABjU2GqJiIiIiDSnsrNwHAngOAD3pXk9G1clEREREZHmVaYB/V8AdgbwQzObRnIIgEmNrZaIiIiISHMqM4iwL4CFZrY0fdwdwCpm9kEH1E9EREREpKmU6YF+GMCq7vGqAH7XmOqIiIiIiDS3Mg3o3mY2P3uQpvs0rkoiIiIiIs2rTAN6AcntsgcktwfwYeOqJCIiIiLSvMrEQO8A4HYAbyJZiXA9AIeb2TONr56IiIiISHPpUWtjOmBwdwBbANg8zX7ZzBY3umIiIiIiIs2oTA/0U2a2UwfVR0RERESkqZVpQP8EycIptwNYkOWb2ZTGVk1EREREpPmUaUD/PpBtZrZ3Y6okIiIiItK8ChvQIiIiIiKSqzmIMEPyCwC2AtA7yzOz0Y2qlIiIiIhIsyqcB5rk9QAOBzACyTR2hwLYqMH1EhERERFpSmVioJ83s23cfT8AvzGz3TumiiIiIiIizaPMSoTZqoMfkFwfwGIAgxpXJRERERGR5lUmBvo+kgMB/DeAKQAMwA2NrJSIiIiISLOqaxYOkqsA6G1mcxtXJRERERGR5hUN4SC5XXWemS3yjefQPiIiIiIiXVm0B5rkVAB7IZl5I+ZhM/t0A+olIiIiItKUasVArwbgGdRuQL/bvtUREREREWluWolQRERERKQOZaaxExERERGRlBrQIiIiIiJ1UANaRERERKQOhQupkOwGYCiA9ZGsSvgXM3un0RUTEREREWlGtaax2xjAWQD2BfAPJDNu9AawGYAPAIwDcLOZLeuYqoqIiIiIdL5aDejbAFwH4HGr2onkOgCOAjDHzG5ueC1FRERERJpErQb0ADN7v4PrIyIiIiLS1GoNInyW5BEdVhMRERERkRVArQb03gAOJ/kQyU06qkIiIiIiIs2scCVCkvsDmAjgaQCVAYNm9qWG1kxEREREpAnVnMaO5OYAzgTwOICxcA1oERFJkOwNYFH1gGsREemaog1okpcCOAjAt8zstx1XJRGRFQfJVQF8G8BMANd3cnVERKQD1IqBXgLg077xTPK+xldJRGSFsgzApwD06+yKiIhIxyiMgW6xM/msmX26gfUREREREWlqtXqgQbIbycNc1rMNro+IrCRIPkpyDslVqvKnk9y3zL6BMk8nOZnkIpITA9vnV92Wkry67PFln0PVPquQvInkqyTnkXwuHZzt91mD5K9ILkj3O6rs8bWOjdSn5v5teA3arbwSZZV+ru1ZVtH+9ZYlIiu+mg3odJnu77rHxze8RiLS5ZEcDGB3AAag5ow+9ewL4E0APwAwIbTRzPplNwDrAfgQwJ1lj29jvXoAeA3AngBWA3AugDvS4zNjAXwEYF0ARwO4juRWJY+vdWxI0f6lX4MGlFdUVj3PtT3LKtq/3rJEZAVXZhq7S5EMjrkdwIIs38xmN7ZqItJVkTwfwOcB/BnAZmZ2oNs2HcCJZva7on1rlP8DABua2fAa+wwDcAGAjatnzyh5fN31csc+D+BCM7uLZF8AcwBsbWZ/T7dPAvCGmX2v1vEAflvPsfWcq+Rr0G7lFZVV57narayi8gBcVE9ZItI11OyBTh0O4DQAfwDwTHqb3MhKiUiXdxyAn6W3z5Nct532rccwALcsx9RzbapXut9mAF5MszYDsCRrfKWmAgj2YFYdX9exbdi/SHuWV1RWze0kryV5bXuUVWd57f2aisgKoOY80ABgZkM6oiIisnIguRuAjQDcYWYzSf4LwFEAfrI8+9ZZh42QhESc0Mbj21Qvkj2RNLhvNrOX0ux+AN6v2nUugP5Fx5Pcveyx9Z6rpPYsr6ismtvN7NT2KqvO8tr7NRWRFUCZHmiQ3JrkYSSPy26NrpiIdFnDADxoZjPTxz9P85Z333ocC+AJM5vWxuPrrhfJbgAmIYmVPd1tmg9gQNXuAwDMK3F8zWNJHs18wORvyp6rxnNo1/KqFJVVz7nas6yi/dvzNRCRFURhDzTJCwDsBWBLAPcD2B/AEwBuaWjNRKTLYbLoyGEAupOckWavAmAgyaFmNrUt+7bBcQAubcuBbakXSQK4CckgswPMbLHb/HcAPUhuamb/SPOGIg/xqHV8zWPNLAsxycrpW3SuWtq7vCpFr0Ph69Sgsor2r7csEekCyvRAHwJgHwAzzOy/kPxhWK2htRKRrurLAJYi+UG+bXr7JIDHkTRq27ovAIBkDybLandH0sDtTbJH1T67ANgALWffKH18W+oF4Lp0ny+a2Yd+g5ktAHA3gNEk+5LcFckqsJOKji95bF3nKvkatHt5RWXV81zbs6yi/estS0S6CDOreQPwVHr/DJLLUgTwUtFxuummm27VNySzRvw4kH8YgBlIropNB7BvmX0D20YhmVbO30ZV7TMOScMnVL8yx9dVLySx0gZgIZLL/dntaLfPGgB+jWSmo38DOKrs8bWOjTzHmvuXeQ0aVV6Jsmq9TtcDuL49yqq3vHrfA910023Fv5WZxu5aAOcAOALAt9M/3M9Z0hstIiIiIrJSqXcp78EABpjZ8w2rkYiIiIhIE4s2oEluV+tAM5vSkBqJiIiIiDSxWg3o39c4zsxs78ZUSVZ0JCeZ2bGdXQ8R0fdRpNnoO9k11BXCIVIGyQVm1rez6yEi+j6KNBt9J7uGMvNA9wRwCoA90qxHAYyzlvOYioiIiIisFMrMwnEjgJ4Abk6zjgWw1MxObHDdYvVpSJf59lutmT/ouRYA4JV/vFzJ+sTGHw8fmO775it5SPiSpfnmZa62y5YF8iLpzNKlrfOq982S/q2MpgN5IiIisuIjW6cZ2x7I99u7Bcry+aFzAUB3v2+31mX16J6n1/+EG263OF3Y1ZbkWfPerKR7DtjI1ThpxDzz/L/RaGussYbNmjWr1bophT3QAHY0s6Hu8SMkl2cFsKY0+c4v5A/WPwkAcOgBu1Wy7rzz++EDBx0PABh1VM9K1sz3883zF+bpBWl6oeu7/2BRnv7wo9bFz10QPu0Cd9zStGG+yJXrG/E+vXhp67zseEANaxERkRWJb5D6dK+0hecbt7171t7X5/VZJU/3DOT7fVdx5fbrnadX7ZXc9181z1uzf56+8LY/5w/empDcfzSjkjXjsQsq6fU+d16+ryWNGG5wMhptyJAhDOWXaUAvJbmxmf0LAEh+AskqXCske/PGSvrgg/NO9K33y1cmf29Bkn79+XHFBS6dDwAY9WO3bw+3UOM6h7c6ZOLI/L2YPT/Pn+fWJ1uYNqZ9o9qnFy9x6fTd8I3xUKPZH7c40sD+aEnrfL/vMtfY9g3vbB9/vLcssK9vrPvj/Jd9SWDfZZHe9aw+S1bYT6eIiDRCrLc06xn1jcFYz2vWSPT/Y/xxvpHZPbDOczeX5xuy2Tn8Mb5B6uuT7eO3++N6uVZddj5/Lr+9Z4/W+VmDt3rf3i4/a0D3dQ3stQbk6eFjAj1xb45vnRfTLW+Br7ffrZX0mHOPqaRveDC5P+gz+WG//lXevuP6jQ+SKNOA/g6A35N8BUmP/0YAtIiKiIiIiKyUChvQZvYwyU0BbJ5mvWxmi2odIyIiIiLSVUUb0CR3BPCamc0ws0UktwXwVQCvkhxlZrPbckKS0wHMQxIGssTMdiA5EcBEAKMADDez6W0pu5QeAyvJu/+UX2ZYcEd+jeScScn9i3fnsTVbfelHlfTxx59VSd94erJPtwPequS9MnFQJf2J4z9fSV80bHUAwHmXja3kfXfEaZX03A/yambxzD5sY1EkdjoLXfBhED7Ew5eRXSpa6PJi4RNZvi8rFLYB5CEaSyMhHj7sIhSWISIi0iixgfXZv6kyoX8LFrbOYzA6NtcjEtYRCtHwZcXCQbKwi16R1puPW86eZ18Xk+zDQULhHD6Ew9fBh3BkMc7+XGeen8cnv3NLfpJ1jkjbRpNdrPJ//iNPPz6wkjx+dDLoa8L4syt5iyfn7a1/5OMJMTYtbs9z3Js5+wF0pECUTsU4AB8BAMk9AFwK4BYAcwHUEcwS9B9mtq2Z7bCc5YiIiIiIdKhaIRzdXS/z4QDGm9ldAO4i+Vw712Muksb6bDRogGJlur63bqrknf6F/FfSW3Pyfe+6Kw1EH3RCsKwJP3dTqSx8Nbl/545K1rS38823HLt6JT3q4vQX2iobVPIuuzv/9TTr1rw+tz6W3PtfvH5GD98bnfUUh/Kq06FBhB9FepiXBGbs8Pv6X/FZebEZPUJT+8XOGxpc6H+Bx+qQlRsrKzTgcFkgD2j5Kz07n+8d8IMiRUSk/flBd6Fe2tjUaFm+//se67H1xxWVm527xUC9wAA/v6/PKxpw6J9vUQ90z8jMGz4/tK/v+V4lMLjQ90D77X4WjY+vndwffIl7gWfcXEmus98vKunzjk2uyK/j5lYYscMj+YMNTq0kJ9yxT5IYuGf+HD5xcSV9zc55++3QQ5JBgrc/kb8B196UDyL0UzSz6DJBG9Xqge5OMvuY7APAPeNSgw9jDMCDJJ8heRIAmNlIM/ujmR1sZq8tR9kiIiIiIg1VqyF8G4DHSM4E8CGAxwGA5CZIeozbajcze4PkOgAeIvmSmf2h1gEkJwE4eDnOKSIiIiJSl2nTpoFktirH3WZ2LFCwEiHJzwIYBOBBM1uQ5m0GoJ+ZTYkeWBLJUQDmm9nldRxTc+iZn+d578/n8wC++m5y7+danuPmYA4NILA38rmd994vD4B/1y2U8sITY5JEd7es/bpH5+m3f5anl6Uj/wbukef13bqSvGlEfpnh7feSez9Y0Kf9YizZgMCFkYVUfGhHFqbgBxbGQjQyft9QiIdPx0IiQgu3xBZzCa3euDQSMhEb1Fh9vIiIrByyK/a+ecOCUItY2IYXCgNhZFBedo5YWb4OWdhErCxftyysIjY40R+X7RsL4QgNIvRzO8dCOLJwjBOvcS/wotfz9Mx783QWyLDecXnejHzdDb/q4Kd2TSZVWNvNKf3Ib/N2WGjRlBaLvbjns97APL1umn7sobbNE7399ttj8uTJreJAas3C0c/M/q8638z+XrXP/Op9apTZF0A3M5uXpj8HYHTZ40VEREREOlutGOh7SP6Y5B5pYxdAshIhyRNIPgBgvzrPty6AJ9KlwJ8C8L9m9tv6qy0iIiIi0jmiPdBmtg/JAwCcDGBXkqsDWALgZQD/C2CYmc2IHR8p8xUAQ5ejvuFy/+CGd87Jxzo+MsVNrZHN/+xGil557vBK+ps3uDiHd24HALz3cL5s5CPPu0sVH/4zT6+6SXLvlqkcPyK/1nHS6Ptc3R5O7uc+medtmYd4nHCqu8SRXu648bx9K1l+2W8fzpGFXSyMhFqEwjViS3mHwjliczuHQjB8CEgoFMPvGys3lI7N0RkK/fD7xmbOyJ6/vzQWCyPJLgPGnk92Dv/cu0UG/WZhNv7SWez1zy7RaWlyEelMoWWwY3/jvNjMF5mi5a5DS2qHjqnWMxA2EQulyPJjy3f782Xp2PaimTVCZfl0KASkuu6hfWOzcGT5PhQjFs6RzfPsZ+FYrU+ePvmCfEltfPgvtPLWhDw9f2qe/vh3AQDjTs/bRSefm4dl+Fk4Xnh++yQxwK3P/Zv8xfYhutksaWO+lm8fedHEfPt6w1rX8ak8ZNamHlBJc+j9rfctoebH28zuB9C2kkVEREREuqASv+VERERERCRTcxaOZhSahaNFt37/7fN0v23z9NI0/qF7vzwvDdUAACye1fpk7tKCX4DFL7AydUxy+WDooe6SxPontT4vALyXro6y5hfyvA9eytN9tsjTy1qvGfrrc/NLIH5RldAsHEWLqhSFbfh0aGGT6uOy8IfFsZk1AmEZft/QIifV58gUHVdmxo5QWEZsAZZadfFlFB3v9w0tABM7bmkgXKS6Ptlxse3+El/2HvpLlssi9c3OHQt1qXVM2XOULTc2ot4LvUf+UmjRa+2F6hu7bF1Urq9vaHGf0CVyIBwqVKbcItm+/nNRFKIUO1coP/b+hI6LXTqPfY7K7hu71O+fZ2hxDl/HUBmxz3EwHCHy+Q/NAlHmuFB9fVn+dcgu1YeOry4jVK4XO1+o3NBx3SNhGVndfFhB0WwZMaEyWny+/WsTmOEi9nkpWkjFf158SETo+KJwD79vLJwjq3toZo7qMnoGQjh8uT5cIwvh6OeW/f7y91/IH7gZwyo+clG8vdbL0+//2eWvm55gcJ7nQl5/My6fWWP/C1sveNdiQTuf/1G6Yl0/FxXs21Z+VpBVNkzuF07P87L2GAAOGo5aYrNwqAdaRERERKQOpVYUJLkbgE3N7Kck10YyD/S0xlatmC2ZlyR8r7L/1fHaFXm655rJvf+VNP2iPL3dE3k6G3D4+MA8b/f3gnUYOuzRJOF/2bje49vO6l9JH3l+Oivgvy/L93XB+Ndckf8qO310GrDv5pT2vyJnurmos95oP1jQ9xSHeqPL9ECH5nb2Zfn80LLfsTmas1/ssYGDoSXCY8uCe0WDE72svNDc0dXnCy37HTpvbNnwogGJRQMZfZ7viQjOXR6pox90WiTU4+0/F74HI5Qfmks8JtbDGuqpbOvFsvYchNnWacWL6h777BSVsbwXEMu8Ntk52rtele9C8a51ve7Zvu3xvod6PWO970U9yf57FVp22v8t6hm5MlAZmFZwXiD/zocG1LWqZ6Bcf95QfWID4vzrXjSIMHsd/P+uWH1DA/h8fRcFeq5j+/q/Udk+sfP6MrL6+v+DsXKzOsRep9BnoGek17lofunQeYG8h9nX18+b7D9z2fs60C1p0aLXedb/VpLXnHMgAOD0k3fJt6/91Ty9xucqydvOHQIAOPJSN6lD1isNYP9T3eDEjO91dgP/sNNf8nTW5vKDFxe9kad7uMklsjZZ1rYDWgwyNMvT9Sz7XdgDTfICAGcBODvN6gkg8IxFRERERLq+MiEcXwHwJQALAMDM3gTQv+YRIiIiIiJdVOEgQpJPmdlOJKeY2Xbpoip/MrNtOqaKrepTqbC9NbHV9uOPGV5JT7jVbc9CN1bbtZI16pj8d8CoSy7M982C0tc6KFwJP8AvGyQ48578vJfkS0T6pcOzS2pnXzwmz3T1Ad01mWxAYTcX0f+3/DLDhOvypTAXBEI4FhUMKPSXfPylxVAIR2zf0CA2P5CxaGBgbN7l0L5FA/iq61nZXjC3c9Hczz4/FmYS+gr5S2NFA/BiISehK0mx55Md5y/fxcIuss9J7DJxPQMZGxFWINJV1DPwslsdI5JiYQGhcK+iQZG+PrHQEF9GFjoQu9JdNMgwNKg3NlgwNviwui6xOvjj/XOIhUqEzuVDJ0ODHotCNPzzCc3R7I8rszx3to/f3jsyiDBL+8GCviwfzpGlTxjmJmLYYXKeXvJenv7g5eQ+C40FgNkPVpKXnH1aq3L9UuDHn+0mfhi4p3sig5N7RmJ/3s7XzcCSuQCAUWfm5xo1yYWJ+Lmo03CP4UcPr2RN/NnE4ClCAwqXZxDhHSTHARhI8msAfgfghhLHiYiIiIh0OYWDCM3scpL/CeB9AJsDON/MHmp4zUREREREmtAKPQ90JYQjtGRjtXQJ78u+N7ySNcdN0XzJ1fklhbNHJCEYs+bl28ff4sZNupkxgvxchdn8gwDQd6vk3o8E9TOIBDw0Or9q4Gfe8LMqZJfkY2EbPpwgyy8zC8figpk1QqEdyyLbvdCy36GZN4D8UlpsifDQ7Bxl5voNzb4RCxMJhXt4oRCPotk/vKJ9Y+ETRfNax85bT4hGqG5FfzJiYTH+8mVs9pNa5yi6BC7SHrpFLvtn6pkzusy+2T6x2XaK5sAuM9d0xl/e93/vsjCFMnXIzuFDAWLHhUJDQttDISLVxwXrFbmGHprb2Zfr65s9j1h4S9FS3rH6ZuX69yQWahGaXzpW9+zvaGwpb5/O3tfY8t0+vWpgKe813Wi3z48KfBkC61UAaDkrWRpqgblP5nkfOyN8XMbN+HHKMQdW0j4M5LKxSZvtrNPzkNk1XHPqrB9NzB+sc3hy70NiPRcawvWOabU5FsIR7YEmOQ9A9F+TmQ2IbRMRERER6aqiDWgz6w8AJC8C8BaASQAI4GgAgzqkdiIiIiIiTabMLBxTzWxoUV5H2WHrtW3yXQcnD7JQiS1uih+QmXlvcu9HZm50diU58qT8MsCY8cmlgW+enOf5cA8/s0YWWnDfVRvlmTtPz9PTzs/TQ0Yn9/6yh7+k8Ma1eXq945J7PxrV7Xvfea1DO/wsHD692IU/ZLNk+BAPfzkrlB9bPGVxINyjxeIpkdk7ihY8Cc3CEZulIzQrRSwkIhRSEpsppGhZaV9Wt8Al2FgoxbJAubEZOxq17Het42PniF0aDr1msWV9Q0uL+zx/OTb0OomsKLLPd2zxlNBiI7F9Q6Efse9j6O9Ht0i4Qeh8ZRZSyY7zIQSxhaGKXofQMbFZOEKLyLQoI/A6xGbACIXAxGb0CIVoxMIrQvv6c/nwCC8URuJnCvF1z/b1z9eHhoSO87Nt+LANn5+FP3zxB/6fomv4eDPSWcA2ODW83Xv9quR+w2/kef/MQzj2P/wnlXS2jHhv9zr5hV2uviEPtc3aZ1eOy/MWPJ632fru5GZW6+3aZwEfPTm8ku61fXIch1xQyVueWTgWkDyaZHeS3UgejXRO6LYguR/Jl0n+k+T30rxHSQ4mOb2t5YqIiIiIdIQyDeijABwG4O30dmiaVzeS3QGMBbA/gC0BHElyy7aUJSIiIiLSGcpMYzcdQGRFkbrtBOCfZvYKAJD8RVr2bABLAbxbWIItBj6akaT3ikzzELLmFwAAf79x7UrWZlvnE4aPueMf+b6rbgIAuPJetx57OosHAEy8ZHglPfyscQCAcT84uZJ38qDR+XHZzBtAHrqRXf4AAAs/h+tOT4a/nnLNvOD20KXz0CIosfwyoRah8AlfbvBSf4kwhlrHV6uM0nbltrjMGDgmtnBJLPQgdBwK9g1N0B+7POrToVHjvSL1rezqL2+WmGWjVh2B+AIIoX2z97tH5JjQzBqx1z90iTRWl9Al3TIhJ0VlFYXmFL1OHT3jRzPMNlL0+rdX+dXnKHr9Y5+dooVJQiEARZ9TIP8u+EvgRX+3YmWFwglii4mEFgBZFqlvKGyrTGhCz0CoRWyWjV4FrYbQ+WIzVYTew9gCT6G8ehaG8Vp8jgJlxP6mFP3tDIXmeLH/R9n5/N/3pZH/eaG/YT4MMFRH3wbw510cmd0q5NrT8ik5Tj1vbJLwYac+3DQLQQXy0FO377WX52Ebv/mfcZX0hEuTdtTx33WLzfnQD+fKSelic73WreT1dSEpLeqQLdAS0eurbiY3/5wKFDagSf4Ugdk4zOz40mfJbQDgNff4dQCfMbM0qBk7tqFMEREREZEOU9iABnCfS/cG8BUAbzamOmEkJwE4uPJ4m2RAYPYrcvg++U8uvzzjsKOGV9I3/y4Z+bfZKZGfWVmgOwB0T6LWr74gD0gfccoelfTwMX5UV9KrfPK5edaNF+e90SeePTbfkA1gHJiXhTmPVJLXjM6XpDz9kKwu4XmiQ0uC+l+8sXl/QwPTYoPUQmLHZb96e0YGloR6eFr8wo4MsMmSsQE4LQacsPV2r8Wv+MDxsd6k7Bd7tEegYPBiUe9BPT17sTpUejBcnu91CD3PWG+ef/1Cz937KPBexAYyhnpBSg1kDCxTHvsMVI7xx5fohWrL9o6wMtQhVv7yvj+xgag+u1vgcxj7jmVif19afA4LrgiFeiJb9KxGemlD5459N0Pfx9h8zNl3M3YVLdRbGiurnr/1Ie0xH3ZWRs/AMtxAfEBh0XmzdLfI1YLQVYTYuULHxZbyDl1lCF1BqD4uaxPE5owuWtLctz9OPSGfN2LMqKStMnKUa9+sunGezqIEgMpV9vEX5+2bU8/Pe52xzmGV5PFXnZQkpudX8cd8LX+hRo52S4AveiO5H+wmaviim+HBTbqQtRF9+7DFst4P5x+O34zP65mZNm0aSPqxf3eb2bFlQjju8o9J3gbgiaLjIt4A8DH3eMM0r6gOxwI4Nj1/E/xLEREREZGubsiQIZg1a1bf6vwygwirbQpgnTbW42kAm5IcQrIXgCMA3NvGskREREREOlyZeaCrVyScAeDs6p7p0ickDwBwJYDuACaY2Q/rOb7FPNCbj6u579Qxedf/0JHJU7julDzvPdch//4HebqyrPfq++SZi2fl6XnP5OlBSSj4uNPya0Wz3dSJZ5+0Zv5gtV2T+61uz/PemhCufGhw4aA87HzciDygf25a99Dy3kDLgQLZpbrYgEM/GCHL9/NAe6Flv/1l0xYDFQPLgS+LbPeyfWODHbpFQgTKlhu7dOtl+8TKCs0DHQsxCM2rHAqZiNVncWTfxQWXYEPhGs0QHiAi7Sf0tyY2v3ooXMD/fekZCbUoWp67aEBm0ZLnZZbUrlW+3zc2iDA0wC82iLPoOcQGVWb5sbmo/fuS5cdCOPzA1dC81X57aNluv6+fi9qnB/RJ7ltMXOBDW91gPSxzjY2Ma5/gDy70dO2vAgAu+dEdlSw/t/Mp17oGRtYe8ufy7bB38rbTOd9IQmz7uaU0VnenPeW6/B/cbWcmL9qRl5f4p/fmeAAAN8hDceteyjuTrUjYXszsfgD3t2eZIiIiIiIdpTCEg+TDZfJERERERFYG0R5okr0B9AGwFsnVkc9GOwDJdHSd4pkXZ4JbJF3sZoEQDjdf89Bj3AQi/74MAPD2e3mWD+Hwy3NXzHG/Ewa5OaF9eMXfhrUq6+xLxwX3/dE5yejO2fNXzfPudCeeslvrOuwwOU8ver2S9EtdZstzx2Y/cNEpwRHHRXPx+u0+nCM0MrjFaOvInNBZuT4yZJXI3JTZTzwrmJvV18fXIbQELeBmZiiY19OfIzY/aSgUokcd+8bmW83ELm/6srLzLQuEaohI11LP/OCx+e19GdnfnVj4hFc0g05bQjRis2mUCe0IbQ+FcMT+HhbNWx2akcPPXhELZcnSLY6PhH5k/1tiS3n7c4TCMkJLffu0D/Hw6d6hcA7f7tnonDx9j6vEliOS+40vq2SddWjerlnDhVKcdfGeAICzL81DMX50dh4egWc+68q9Nbnvs0WeN/OePG35P/zs8/fO3HyzD13Fk+tVkkde/FaScO1DrOfmfnZ86EaRWiEcJwP4JoD1AUxx+e8DuKb0GUREREREupBoA9rMxgAYQ3KEmV3dgXUSEREREWla0Vk4SO5tZo+QPDi03czubmjNIvw80DbtQgDAwskXVLb33m1iJX3kYcMr6dtuT2fW6JV367cY3fn2z/J0/3SJb7+ISbq8NwBc+/X8Usap1yWBCGNOyq+L+Bk9zpvkQjSyEIyntq5kXfbTfDTrgoX5rhf+LBkJe+XJ+RjOme/n2xe6+IcsrMKHV/gZKXx+Fo7gR1v7ff3HIctfEphBozo/Ky82k8XiwIwdsQVcQsuFx2b38EKzcPhyi8qILb+dXVaLvU4iItI+ugXCILwys3SEFjwpmpHD/6/oVbDASChUA2gZJpKFSliJcosWUgmFdsRmAgmdw+f5EA6fn4Vz+Jk5/KwWZ4zPZ+e44OikXdLHLZ19/L55eu0j5+QP0gVNLjo2D/HonyfxzRvyBso1X08qd/o496L5Wc96ulnNslCTtQ7K8xa8mKc//GclefhhyYwdt985sZL38A3DK+l9js8XhOGGrRdSacssHHsCeATAFwPbDECnNKBFRERERDpTrRCOrFt3tJlN89tIDmlorUREREREmlSZhVSmmNl2VXnPmNn2Da1ZvD6tKmxL3MTfc5/M02t8vnV+tpgJgDvPynvkZ7siTj4vCfe44twTK3ln3JSf9uqT3HHpoikXXB5Z1KVgQZQWC6lscGqeTi9PTBqdXxeZ60JDFrrRpgvSKJDQgimxdGwhFR+6kO3jQx986Ij/6GRhGb6sUDiIF8qrzi9a8CQUBhKqV7Wsnv6S25LIvqHFT0L1iS3wEgoNiW1fFvg6xkJdQouj+M+AvxwYWthFYSgiXUvRQir+O9/W2TBC231ZywLnKCq3zMwbWd1jZYVm04gtYuKF6hhboCW0aEpscZTQLBxF+xaFePi6xcI2egRm5AjlAS3DNbLQDp/nFzw59nu35g/WPTq5f+um1nlAeIG4xTPzdJ/NK8lRI4+opLOQkZE35B+iHx+fvzHf/uGNlfR1o5P22ZpupZLDLnMfvgV/ydN907DZyCwcLJiyqu4QDpJbANgKwGpVcdADAPQOHyUiIiIi0rXVioHeHMCBAAaiZRz0PABfa2Cd6sYe4cUS7X/dg82+ldx/8FIl69BTLsy398871EcMOxAAMMv1SuPdfOXyEePdr5z0F82V5w6vZH3zh+6Xmg9wz/gBi643+sbT8x84J45JTu6X544t1Z31RscGDvpe2MWBHuhQj6/Pb7Ekd8G+/odcaOAgEF4a2/eMhHqCQ0uBVx9Xz8DAUM927Lhsn9g826HtRct+Fy313R6KetRFpGsJfbeXRQZdx/4+tEWsA6/sIMBYr3JogF50bmj3Py/r6Y2tYVA0gM//T/T7Zv/TYvv2DAwC9HXwvb+h3nH/noR6nX0ZsQH9ft9snxZrGEQ+D6Grny163F1bJmurnHiVmyRhxi152l9lX/Jecu8GA17x3QMr6VFX5L3KlfU23nuskvWcCyA+/cQ8KuCaG91xGd8jPvvBPJ1e0efus9zOw1sfX6daMdD3ALiH5M5m9qflPpOIiIiISBdQqwc68yzJ05CEc1RCN8zs+PghIiIiIiJdU5lBhHcCeAnAUQBGAzgawN/MbGTjqxesT+kL0P65nXpAcsnh2vvzvCcuza+h3O362K+4Pr000GNgnvnBy3n6xe/n6U2TJS3PODNfa2a+m895/C2BwHsf3P7aT/J0t3xSxcsvuS7ZNRLC4QcRhpby9vuGQjB8iMGiyMDAjwL7xgYGZpeQYvM5h8InYoMIlwUG1cXK8rL8ou0+XWYwX1afRoVaiIhIbdnf39i8y7H5ljOx4zKxOZi9ULlFczfHBjqGjutZsN2n/XzOsTmhszL8vv4coTmhfZ6f59mnv3vGAekJBuaZm/zYFebW20gncDjx0N0qWQPcPNBXXHxY/mD2A8m9DwFZ8wt5euH0SvKMryfhHOuvkW8+86et23xA3u4rGiwYExtEWLC6PABgEzM7D8ACM7sZwBcAfKZNtRARERERWcGVaUBn/ZPvkdwawGoA1qn3RCT3IjmX5HPp7Xy3bTrJwSQfrbdcEREREZGOVCYGejzJ1QGcC+BeAP0AnNfG8z1uZgcW79Y+fHe9vZmGZbi5CHf7Xt7dv9v85yrpK0Z+GgDw7ty8rEuuyud5PvOGPH/BwiR0w4dt+PAKrP3VPJ1dnvDhIFu4UaPOgD5JCEdsKeoQP0dzb3cZJrQEdexKht83u/wTm3vYXyrKLlO1CINw6W6BeTX98b7cRe64rA6xkAn/PLJ6xpZs9T8Xs30WR2bDgJ/TOT2u6PUXEZH2E5pxI/a/KzR3c4vlvQuWAA8tBV6dX8+y3lkZPSNzO/tQiuw5xcI2QrNw+O1F4RyxsA1fblZPH6rh54RerU+extbprGTd3IzGr12Rpwfumaf7DQUAfNnFLdzplus46Yw7KukB6Tkuv2rDSt4/bs/Xwrj3qfy4KybclyR8iIdz7U35LB1tDd0oUtgDbWY3mtkcM/uDmX3CzNYBMLPouDq9C2ApgNntXK6IiIiISLsqE8IR8pPiXYJ2JjmV5G9IbpVlmtmOZvaamR1c62ARERERkc5WOAtH8CDyNTP7WJ3HDACwzMzmkzwAwBgz27TksZMAZI3rPrX2jak8z6XzK3kf3p0vwLLqofnrcMlxSXe/D+G44p58+4+G5ZcDzvrhmCThFmJB743y9Mx783S6VPekb+XHv+PO4cNAPggsz+0XTwktmhJaMKU6P5vNwpflZ7gIzaLhtxfNyBFbRjs0M4af/cNfMgvNzhGb3aNoSe3YoiuaRUNEpGvrEQiTAMJLh4dCKoDw7B1FYR3+3EVhGz4/Vt/QEt8twi8CoRixcn1ZodAPH8Lhy+3twjmyffq6fddeLU8Pu9L9g80Wjus92FVi3Tw955FK8uKzTgYAnDMpP/6Mg/IXYq0B+WHn3JI0IJ76cf4kdjoz/I99eUM41lhjDcyePfuD9OHdZnYs0PYe6FLND5KnZYMGAfQzs/kAYGb3A+hJcq1SJzM71sz6mlnf4r1FRERERJbfkCFDkLVBs8YzUGMQIckXEG4oE8C6gfxWzGwsgLFpeeuRpJkZyZ2QNN5n1SygPaVLPH7xoHwpyP+558ZW2wHg5TeS+4mT8oGDeHN8Jfm6q/UxxyTTYd/6e9ed+oD7WTfkpEpy5BeTX0Fjrh9byRtx8mmV9Dy3KmbWO9ti7me/PPeS1vmx+ZxDy3ovicwD7YXmgfZ8HYqW8hYREelI9SxXvrB4l7pk/wv9/+KiuZ1jgx5Dvdi+9zjUmw24QYSBXunqMrJ9fR39IEK/b5Y/0HVpnvkD1645MK/Q1RftkSQ+ejvfead8LYyjjjy5kl6jX5pw7a0rrsvbYcOOyffF61clRR33i0rWgTvm573v3sBS3+2s1iwc7T1bxiEATiG5BMCHAI6wtsSPiIiIiIh0omgD2sxebc8Tmdk1AK5pzzJFRERERDpamwYRdqZ6lvKuh72RXyY49NDkMsGdd7oQjjU+l6f9pYj5U1uVNerb+WWGOfmYxUqIhh8s6C8x+YF9WdqHV3wUCJnwZfjjY0tuh8I1fJhIaM7n2EA8ERERaX8MDHQE8gGQsfmcuwUGH/oQjt6R0I/QUt4+hMOHkWSDCP2ARD/IsF/v1mkf7nHhFa5t5WWDC/080m9NqCRf+NW3K+lPHZKUwQ1cWEeDLM9S3iIiIiIikirVgCa5KsnNG10ZEREREZFmVxjCQfKLAC4H0MvMhpDcFsBoM/tSB9SvlbXWWssGDx7cGaeWkqZNm4YhQ4Z0djVEBPo+ijQbfSdXLNOmTbNZs2a16nAu04B+BsDeAB41s0+neS+Y2acaUlNZ4ZFcoDm7RZqDvo8izUXfya6hTAjHYjObW5WnoWQiIiIislKqNQ905kWSRwHoTnJTAN8A8MfGVktERNrJ3Z1dARFpQd/JLqBMCEcfAN8HkM3j9gCAH5hZey/cI10EyUl+uUsRERGRrqRMA3o7M5vSQfUREREREWlqZRrQvwewHoBfArjdzP5S8wARERERkS6s1EqEJNcDcBiAwwEMQNKQ/kGD6yYiIiIi0nRKLaRiZjPM7CoAXwfwHIDzG1kpaX4k9yI5l+Rz6e18t20/ki+T/CfJ77n8R0kOJjm9UyotspIIfQf1/RPpOCSnk3wh/f84Oc1bg+RDJP+R3q+e5g8nOSq9De/UiktphQ1okp9M39QXAFyNZAaODRteM1kRPG5m26a30QBAsjuAsQD2B7AlgCNJbtmZlRRZmeg7KNI0/iP9/7hD+vh7AB42s00BPJw+lhVUmWnsJgC4HcDnzezNBtdHVnw7Afinmb0CACR/AeAgAH8FMBvAUgDvdl71RLq82HdQ3z+RznUQgL3S9M0AHgVwFoAPAcxP8z/s8FpJmxQ2oM1s546oiKyQdiY5FcCbAM40sxcBbADgNbfP6wA+AwBmdnCat2OH1lJk5RL8Dur7J9KhDMCDJA3AODMbD2BdM3sr3T4DwLoAYGa3d1IdZTlEG9Ak7zCzw9LQDT/SkADMzLZpeO2kmU0BsJGZzSd5AIBfA9i0c6skIiLSFHYzszdIrgPgIZIv+Y1mZmnjWlZQtXqgR6b3B3ZERaT5kTwNwNfShwdkIT1mdj/Ja0muBeANAB9zh22Y5olIx9B3UKSTmdkb6f07JH+FJLTqbZKDzOwtkoMAvNOplZTlEh1E6C4znGpmr/obgFM7pnrSTMxsbDZoEMAykgQAkjsh+SzNAvA0gE1JDiHZC8ARAO7trDqLrIT0HRTpRCT7kuyfpZGs5PwXJN/DYeluwwDc0zk1lPZQZhDhfyIJcvf2D+TJyuUQAKeQXIJk0MMRlkwqvoTk6UiWfO8OYEIaGy0iHcDM9B0U6VzrAvhV2sfUA8DPzey3JJ8GcAfJEwC8imR9DVlBRRdSIXkKkp7mjQH8023qD+CPZnZ046snIiIiItJcajWgVwOwOoBL0HKuwnlmNrsD6iYiIiIi0nQKl/Im+VkAL5rZvPTxAACfNLM/d0D9RERERESaSpkG9LMAtkvjW0GyG4DJZrZdB9RPRERERKSpFC7ljaSRXWllm9kylBt8KCIiIiLS5ZRpQL9C8hske6a3kQBeaXTFRERERESaUZkG9NcB7IJkIv5sWeaTGlkpEREREZFmVdiANrN3zOwIM1vHzNY1s6PMTKvniKzkSM6vejyc5DUFx1T2Ifl1ksc1uI5tOgfJgSRPdY/XJ/nLdqjPYJJ/Wd5yCs7xR3euo9pwfOH7GDjm0yRvStPrkryP5FSSfyV5f8Gxy/WakLyK5Pnu8fdJjk3Tl5PcO3LcoyRfJvmltp67PZHcPX29Gvr5EJH2UdiAJrkZyYezLzXJbUie2/iqiUhXZmbXm9kty1sOyeiYjOU4x0C4FVfN7E0zO6QN5XQ4M9slTQ4GUHcDuo3OAXBVmh4N4CEzG2pmW6LlNKjLLfB+nwtgOMlPkPwEgBMBfD/ddnXB+Y82s9KrNNb6rC0vM3scwAGNKl9E2leZEI4bAJwNYDEAmNnzSJaGFREJIrk2ybtIPp3edg3sM4rkmWl6W5L/R/J5kr8iuXqa/yjJMSSfI/mXdNn47NhJJJ8EMCntxXwkPf5hkh8PnGNjkr8l+QzJx0lukeavm55zanrbBcClADZOz/vfvpc0redW7nk8SnKHdPneCSSfIvksyYPqeL32SY95IS1jlTR/OskLSU5Jt2V1XpvkQyRfJHkjyVdJrpVuy64MXApg9/Q5fKu6ZzntJd4rTf8Xyb+TfArArm6fMu9jfwDbmNnUNGsQknA/AJX/GWDiv9P38QWShwfKGpy+N1PS2y5p/l5p/r0A/uqPMbP3kTSYr0lv55vZe+m2VwGsSXK9Eu/B19LnODV9zn3S/Ikkryf5ZwCXkdyE5O/S/aakn6tBJP/gPqe7p8d+juSf0v3uJNkvzd+R5B/TMp5KX0MRWYGUaUD3MbOnqvKWNKIyIrJCWTVtMDxH8jkkPY+ZMQB+YmY7AvgqgBsLyroFwFlmtg2AFwBc4Lb1MbNtkfQIT3D5WwLY18yORNLTeHN6/M+Q94Z64wGMMLPtAZwJ4No0/yoAj5nZUADbAXgRSa/lv8xsWzP7TlU5tyNdgpfkIACDzGwykkbcI2a2E4D/APDfJPsWPG+Q7A1gIoDDzexTSGY5OsXtMjOdNvS6tN5IX59HzGwrAL8E8PFA0d8D8Hj6HH5S4/yDAFyIpOG8G5LXNVPmfdwBgA87GAvgJpK/ZxJOsX6afzCAbQEMBbAvktdnUFVZ7wD4z/T5Ho6W7+N2AEaa2WbVFTCz25As/DXAzCZVbZ4C96OghrvNbMf0c/A3ACe4bRsC2MXMzkDy+Rqb7rcLgLeQ9PQ/kH5OhwJ4Lv1Bcy6Sz+h2ACYDOINkLySfoZFpGfsC+LBE/USkiZS5HDWT5MYAsnmgD0HyB0NEVm4fpg0GAEnsLJLGFJA0CrYkmW0ekPW+VWOy6ulAM3sszboZwJ1ul9sAwMz+QHIAyYFp/r1mljU8dkbSQAOASQAuqzpHPySNnTtdnVZJ7/cGcFx6jqUA5jLtAY+4A8CDSBqxhyFpwALA5wB8iWmPN4DeSBq2f6tRFgBsDmCamf09fXwzgNMAXJk+vju9f8Y9x90AfCWt829Jzik4Ry2fAfComb0LACRvB5A1UoPvo5n5+PdBAN7NHpjZA0xCKfYDsD+AZ0lundb5tvQ1fpvkYwB2BPC8K6sngGtIbgtgqasHADxlZtNCT4Dkhmk9lgXq9w6A9UPHVdma5A+QhO/0A/CA23anmS1Ne4o3MLNfpc91YXr+pwFMINkTwK/N7DmSeyL5MfJk+vr1AvAnJO/3W2b2dFrG+yXqJiJNpkwD+jQkPTdbkHwDwDQAxzS0ViKyousG4LNZAyPjGmL1qF7tKXu8oM76vOcb/G1lZm+QnEVyGyS9pF9PNxHAV83s5eU9R5VF6f1SLN8c/EvQ8qpj7xLHBN/HKh9Wl2VmswH8HMDPSd4HYI+SdfwWgLeR9OJ2A+DPuwAASH4MwP+kedeb2fVIesovAPDJ9N5fNeiNcj28EwF82cympj8G96o+d0z6424PAF8AMJHkFQDmIIkFP9LvS/JTJeoiIk2uzCwcr5jZvgDWBrCFme1mZtMbXjMRWZE9CGBE9iDtUQwys7kA5mRxowCOBfCY2+XwtIzdAMxN96/2R+RjM44G8HjVOd4HMI3koWlZJDk03fww0pAJkt3THvF5AGrFpd4O4LsAVstifJH0WI5g+iuB5KdrHO+9DGAwyU3Sx9XPP+RJ5GEkn0MSvlCt+jlMB7AtyW5pI3SnNP/PAPYkuWbag3qoO6bM+/g3AJu4ffZ28cP9AWwM4N9I3pPD09d4bSSN6urwwNWQ9M4uQ/I6dK8+mZm9loalbGtm15PcH8A6SMKALgJwMEkfhrIZWoaYxPQH8Fb6Ghwd2sHM5gF4neSX0+e3Csk+JDcC8LaZ3YAkzGU7AP8HYNfsfWUSI78Zkvd7EMkds9eIDRycKCKNEf3SkjzGzG4leUZVPpD0AM1Gcgl1eS4dikjX9A0AY0k+j+TvzB+Q99SGDANwfdrwegXAf7ltC0k+i+Ty/vGR40cA+CnJ7yAJJ/DHZz3WRwO4jsksQj0B/ALAVAAjAYwneQKSXt5TzOxPJJ9kMnDwN0jier1fIun1vMjlXYQk7OJ5kt2QXK07MFDXzUm+7h5/K63vnWlD6mkA10eeZ+ZCALeRPBZJWMAMJA1m73kAS0lORdK7emVap78iafROAQAze4vkqLSc9wA858oofB/N7CWSq5HsnzYwt0cShpH1eN9oZk+TnIwk1GYqkvfku2Y2g+RgV9y1AO5iMvXgb1HQ85vGj18J4JB0xdwF6WfgGgB7p43hTZDEHxc5D8mPiXfT+9gPqGMBjCM5Gsng+kMB7A7gOyQXA5gP4Dgzezftyb6N6aBQAOea2d+ZDKC8muSqSHrH902PE5EVBN0q3S03kCeb2TiSFwR3ANYEsJOZfbZhtRORlRrJRwGcmQ7Sa8vxVwOYYmY/bdeKdbK0QbbUzJaQ3BnAde0RnrIc9fkWgHlmVjRYtEOR/AqA7czsvMC2R7Ecn61GSH9M3GdmW3d2XUSktlrzp45L7y+M7ZP+AhcRaTokL0IyQG5UJ1elET4O4I60p/sjAF/r5Ppch5ahH82iB4AfR7bNRhKvfE49c0E3ShrCdC2AmZ1dFxEpFu2BruyQjG6+Gvk0QI8jmX7n9fhRIiIiIiJdU5l5oH8K4F4k0wCtj2T0c5e6HCoiIiIiUlaZHujnqmPrQnkiIiIiIiuDMj3Qs0gek0491J3kMQBmNbpiIiIiIiLNqEwP9EZIYqB3RjL10B8BfMPM/t346omIiIiINJeaDWiS3QHcYmbBSeVFRERERFY2NUM4zGwpgI1I9uqg+oiIiIiINLUyy4e+AuBJkvfCrQplZlc0rFYiIiIiIk2qTAP6X+mtG+JLm4qIiIiIrBQKBxFWdiQHADAzm9fYKomIiIiINK/CaexI7kDyBQDPA3iB5FSS2ze+aiIiIiIizafMNHbPAzjNzB5PH+8G4Foz26YD6iciIiIi0lTKLKSyNGs8A4CZPQFgSeOqJCIiIiLSvMr0QF8JYFUAtyFZSOVwAAsB3AoAZjalsVUUEREREWkeZRrQv6+x2cxs7/atkoiIiIhI8yo9C4eIiIiIiJSbhWM1kleQnJzefkxytY6onIiIiIhIsykziHACgHkADktv7wP4aSMrJSIiIiLSrMrEQD9nZtsW5YmIiIiIrAzK9EB/mM79DAAguSuADxtXJRERERGR5lWmB3pbADcDWA0AAcwGMNzMpja8diIiIiIiTab0LBwkBwCAmb3f0BqJiIiIiDSxaAOa5DFmdivJM0LbzeyKhtZMRERERKQJ9aixrW96378jKiIiIiIisiLQQioiIiIiInWI9kCTvKrWgWb2jfavjoiIiIhIc6sVwvFMh9VCRERERGQFUc8sHH3M7IMG10dEREREpKkVLqRCcmeSfwXwUvp4KMlrG14zEREREZEmVGYlwisBfB7ALABIF1DZo4F1EhERERFpWmUa0DCz16qyljagLiIiIiIiTa/WIMLMayR3AWAkewIYCeBvja2WiIiIiEhzKhxESHItAGMA7AuAAB4EMNLMZjW+eiIiIiIizUULqYiIiIiI1KHWQipXA4i2rrWQioiIiIisjGrFQE926QsBXNDguoiIiIiINL1SIRwknzWzT3dAfUREREREmlqpaexQI5RDRERERGRlUrYBLSIiIiIiqBHCQXIe8p7nPgA+yDYBMDMb0PjqiYiIiIg0F01jJyIiIiJSB4VwiIiIiIjUQQ1oEREREZE6qAEtIiIiIlIHNaBFREREROqgBrSIiIiISB3+H/hr+NyaBjcYAAAAAElFTkSuQmCC\n",
"text/plain": [
"<Figure size 864x432 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"fig = plt.figure(figsize=(12,6))\n",
"ax = fig.gca(projection=m_aia)\n",
"m_aia.plot(axes=ax,norm=ImageNormalize(vmin=1,vmax=1e3,stretch=LogStretch()))\n",
"ax.grid(alpha=0)"
]
},
{
"cell_type": "code",
"execution_count": 34,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"(23, 244)"
]
},
"execution_count": 34,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"m_aia.data.shape"
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAXoAAAD9CAYAAACyYrxEAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAArRUlEQVR4nO3dd3hc5ZX48e+ZGfVuNavZki33hgu2Mc20UI2BZAmGJARYTEhI2WyyIZv8kmyWTa8kpJgSSgKEUIJJHHq1AYNtbGxjW5a7ZKtZ0qiXmXl/f8youMiWVe7V3Dmf5+FBmpFG5/VIZ86ct1wxxqCUUsq5XHYHoJRSanhpoldKKYfTRK+UUg6niV4ppRxOE71SSjmcJnqllHI4TfRKKeVwmuiVUsrhhjzRi8gUEfmDiDwpIrcP9eMrpZQ6Nf1K9CLygIhUiciWo26/RER2iEipiNwJYIzZZoz5HHAtcObQh6yUUupUSH+OQBCRc4Am4GFjzPTQbW6gBLgIKAPeB5YZYz4SkSuB24FHjDGPnuzxMzIyTGFh4YAHoZRSkWj9+vU1xpjMk32dpz8PZox5U0QKj7p5PlBqjNkNICKPA0uBj4wxK4GVIvJP4LiJXkSWA8sBxowZw7p16/oTilJKqRAR2defr+tXou9DHnCg1+dlwAIRWQxcA8QAq/r6ZmPMCmAFwLx58/RkNaWUGiaDSfTHZYx5HXh9qB9XKaXUwAxm1U05UNDr8/zQbf0mIktEZIXX6x1EGEoppU5kMIn+fWCCiBSJSDRwHbDyVB7AGPOcMWZ5SkrKIMJQSil1Iv1dXvkY8A4wSUTKROQWY4wPuAN4AdgGPGGM2Tp8oSqllBqI/q66WdbH7as4wYSrUkop+9l6BIL26JVSavjZmujt6tFvLvOyfl+dpT9TKeVMW8q9BAIje4V4RB5qdtc/P+IbT31odxhKqTC37VADV/xmNS9vq7Q7lBOKyNZNZUMbu6ubaO3wW/pzlVLO0tUZ2He4xeZITiwiWzdVje0EDGyvaLD05yqlnOXDsnoADnnb7A3kJCKuddPc7qMlVMlvPaiJXik1cB+WBbsRFQ2tNkdyYhGX6Ksa27s/7ivRbz3opT+neiqlnO+Qt5Xd1U3H3N7S4aOksjH0NVrR98mOHn1VQ/AJiXa7+OjQsYl+/b46Lr97Nf/48BAHalsorTr2CVZKOV9Xrvjy4xu55aFjT9fderCBgIH0hGgqNNH3bTh79KVVjRyobTmmMu+q6BeMG8X2Qw00tfuOuP/d3YcBeHpDGZ/903t89YmNQx6bUmpk21PTzPwfvMJ9b+3m/b217Klp5pD3yPbM2lCuuGhqNlWN7fj8ATtC7RdHtm42Hajnwl+8ydk/eY0/rdl7xH1dif7ms4po9wX47aulR9y/bm8tAK/tqGZXdTNldSO796aUGnpdrZofrNpGV624dndt9/3GGP6+8SCnF6YxPS8Ff8BQ09RxxGN88o/vsPinr/HQ23utCrtPjkz0vVfT7KhoPOK+6sZ2otzC4omZfGJuPvev3s2B2uDSqEDAsH5fHbPye95h1DZ30O7TZZhKRZLKhmBBGDAwLiOBpFgPa/f0JPqtBxsorWriqtl55KTEAlDR0NO+6fQHeH9vLXsPt/CLl0psn/NzZKLfe7iFKLcwNSeZQw1H9s6qGtvITIxBRPjcuePp9BveCb0F21nVREObj0+fUciy+QV8bGp28Hsa2o/5GUop56poaMMlMCs/hU8tHMv8wlGs3XO4+/5nN5YT5RYun5HD6K5E36u1U+FtI2BgRl4K3tZODtrcw3dmoq9ppmBUPPlpcUf840Owos9MDj4xRRkJxHhclISq/nX7gq/Y88am8cNrZnL9gjFAcIMVwHef3cJrO6qsGoZSykL3vFba3bqt9LaRkRjDs3ecxc1nFTG/aBS7q5upbQ62Z1aXHmZBUTqp8dHkpMQBR668OVgfzDsXhYrFj2xeyu3IVTd7D7dQmJ5ATkrsMcueqhrayUqKAcDtEoqzEikJrazZXOYlNT6KsenxAGQn97wla2jr5KF39vHi1oohjVUpZT9jDL94qYQ/vxu8BGtFQ1t3pQ4wPS/Yzt12qIHGtk52VDQwd2waAGnxUUR7XEcm+lCBef7kLEQiPNEPx6obYwz7DjdTmJ7A6JQ4Gtt8R6ysqW5qJzOU6AEmZSd1V/RbDnqZnpuCiAAwOpToKxva2RlaL3v4qAkXpVT4a2jz4Q8YdoaKvsqGNrKSehL95NFJQDDRbzrgJWDoTvQiwriMhCPmAw/WB5P++MxEijIS+OiQvSf0Oq51U93YTkuHn8KM+J5JktArbWVDG7XNHRSGKnaAiaOTqGhoo6apnR0VjUzLS+6+LzX0Sl3Z0EZJZfAX4HCzJnqlnKa+Jfh3XVrVhD9gQhV9T0GYnhhDZlIM2w41sn5fHSJw2pjU7vtn5qewubxno2V5fSujEqKJi3YzNSf5uHt2rOS4RL+nphmAsekJvSZJgon+jZJqAM4qzuz++onZiQCs2nyITr9hWm7PuwsRITs5JpTog6/WtZrolXKcrr/rdl+A0qom6ls6u9/Rd5k8OontFQ2s31/HpOwkkmOjuu+bkZ9KbXNH93Lsg/Wt5KYGv39qbjIHalupaWrn1ofX8cg7e60ZVC/9usLUSOcPGK6/910qGtqIcgdfu4rSEzAEX127Njq8WVJNZlIMU3KSur93Ynbw46fWlwEwPTe590MzOjmWCm9bd8vmcJOuwFHKaepbOrs/Xl1aA/TM0XWZmpPMn9bsZXd1Mx+fm3fEfTNDPfzN5V4KRsVTXtdKUUYCABdOyeYnz+/ga3/bxOs7qqlubOfTZxQO42iO5YiKfku5l7V7ahmdHEtLu4/MpBhyU2N7JlO9bfgDhrd21nDOhMzuHjxAXmocealxbCrzkhDtpjA94YjHzk6OPaKib2jz0eEbuTvglFKnrq6l55366p3Bd/69J2MBJuck0eEPIAK3nTP+mPui3MKHZcH2zcH6VvLSgqtxJmYncWZxOq/vCD7ulnKv5Uek21rRi8gSYElxcfGgHuet0BPzuxvmkBYfTcAYPG4XHjeMSojm5W2VPL+1Am9rJ4snZR7xvSLC48sX8pW/biQ3NQ6XS464Pzs5ln9tqcAfMIxNj2ff4RbqWjqOebVXSoWvulBFnxjj6d5Xc/Tf+GkFabgE/vuyKRSMij/ivhiPm8mjk3lvz2HqWzpp7vCTlxrXff9nFxWxpvQws/JT2FTmZVNZPQvHpQ/zqHo4YtXNmyU1TM9LJj0xBpdL8Lh7hjU6OZZNZV7aOv1854qpXD4j55jvLxgVz1O3L+I3y2Yfc9+ErET8ocuEXTA5uCZWV94o5Sx1zR24BE4vTKOtM0CUW8jtlaghuO9m/bcv4lMLxx73Ma6clcuG/fVcdvdbAEzJ6WkDXzgli98sm83vPzUX6DlqxSph36NvbOtkw/46lp8z7rj3XzxtNJNGJ3HXVdNJiDn14X7y9ALOLM4gLtrNrqomHlizRydklXKYupYOUuOj+eE1M9lS7qUoM4HE4+SLtIToPh/j388uor61g4ff3sfP/20WZxZndN8nIiyZlQsEF4Css/ia1WGf6D862IAvYDi9aNRx7//yhRMG9fgi0v02rWvC5nCzTsgq5ST1LZ2kxkcxOiX2mN58f4kIX794Ml+9aBLuo1rAvU0enczGA/UDjHRgwn4yti00MZocO/yvWemhV3Nt3SjlLLXNHYyK77taPxUnSvIAsVEuyxd0hH2i7zoD2uMa/qGkxEXhdom2bpRymK7WjRWiPS46LT67PuwTfdc/WJR7+Ificglp8VHaulHKYepbOkmLjzr5Fw6BKLdW9Kes0x9cERPlPvHbpaGSnhCjrRulHMQYQ11LxwknWodStMdFRyRV9ENxeqWVFT0E1+Vr60Yp52jt9NPuC5BmUesmxh1M9FZejCTs19H7QhW9x6KKPiHGQ7PFu9qUUsOna7OUla0bY8AXiJBEPxS63gJFW1TRe1xCwMInSCk1vLpOrky1KNFHe4K5ysoJ2bBP9N2rbixK9G634AvoWTdKOUXXxGiMx23Jz+tqM1s5IRv2id7qyVi3SPeRCEqp8Nf193yy9e9Dpauit3JCNuwTfYfFk7Eel1jaW1NKDa+uRO+xOtFrRd9/vu6K3qLWjfbolXIUyyt6bd2cus7Q+dBWPUket1b0SjlJ19+zVSv3eiZjddVNv3UGApZV8xB8QdEevVLO0VPRW5NHdDJ2ADp9hiiLqnkITsZqRa+Uc3T9PbtFJ2NHLF8gQJTHyorepRW9Ug7iDy2Xtqr927VCMGIq+qE6AsHK1o3Hra0bpZykq7C2qkcfE2kV/VAcgdDpt7h1oz16pRzFZ3FFH+0ObszqjJSKfih0+i1u3YjujFXKSaxeRx/lCbVuIqWiHwo+v7HsCYLQOnqDpSfPKaWGT9dkrMuqyVi3nnVzyjqs7tGHXlS0faOUM/gtXkffla/atXXTfz6LE7079MugSyyVcgard8bG6BEIp67Tbyw70Ax61tpqRa+UM/T06K1Jh3pM8QB0+gOWHVEMPa/6WtEr5Qw+iyt63Rk7AJ3+gGUXHYGeHr0ebKaUM1i9YUor+gHwBYxlkygA7tCLilb0SjmDz+LllV0/Ryv6U9Dhs3gyVnv0SjlKwOLWjYgQ7XHRrhV9//kC1k7Gerp79LppSiknsPpQM4AYt4tOnx5T3G9Wn3Xj7u7RW/YjlVLDyB8wuARcFm68jPK46PD7Lft5YZ/ogztjrT3UDLSiV8opfAFrcwgEd8dqRX8KOvwBoj3WvRK7tEevlKP4AwaL8zxRHtGzbk6Fzx+wtqLXdfRKOYrfporeylU3nuF4UBG5CrgcSAbuN8a8OBw/B7p2xlrfo9eKXiln8AeMZStuukR73COzoheRB0SkSkS2HHX7JSKyQ0RKReROAGPM340xtwKfAz45tCEfKTgZa+GqG7cmeqWcxBcIWHoCLkC0W0bsOvoHgUt63yAibuAe4FJgKrBMRKb2+pJvh+4fNtavutENU0o5iT0VvWtk7ow1xrwJ1B5183yg1Biz2xjTATwOLJWgHwP/MsZsON7jichyEVknIuuqq6sHFLw/YAgY644XBd0wpZTT+PzWJ/ooi3v0gy2F84ADvT4vC932ReBC4BMi8rnjfaMxZoUxZp4xZl5mZuaAfnjXK6IdPXpdXqmUM/iNPRW9lT36YZmMNcbcDdw9HI/dW1f7xI4eveZ5pZwhuOrG6h59eFX05UBBr8/zQ7f1i4gsEZEVXq93QD+86+K6WtErpQbKZ0OPPsriin6wGfJ9YIKIFIlINHAdsLK/32yMec4YszwlJWVAP7wzlGwtPY9ee/RKOYrf4t31EDrrZiQmehF5DHgHmCQiZSJyizHGB9wBvABsA54wxmwdnlCP1ekPJttoKydjdcOUUo5iS0U/UjdMGWOW9XH7KmDVQH64iCwBlhQXFw/k2/GFXhHtOOtGLzyilDME7JqMDaMe/aAMunXTterGo0cgKKUGxo6KPriOXg8165euf6goC58kPdRMKWfx27AzNtzW0dvKjnX0Ht0Zq5Sj2LFhqmsdvTHW5BFbE/2gl1d2VfQWtm7c2qNXylH8Fl93GnoWkFjVvnFGj97CV2Pt0SvlLH5juluyVokOFadWraUP69aNz4aKvqdHrxumlHICO3bGdrWbOy3q04d1ou/sXl6pFb1SamCCPXqLLzwSSRX94Hv0NhyBoOfRK+UodlT0XTvsA5EwGTv4Hn3XoWbWr6PXRK+UM/gCge4CzqnCunXTdbCYladX6hEISjlLwFjb/rVDWCf6DjtOr9QNU0o5ii8Q6P67dqqwTvQ959HbcUyxJnqlnMBvw4YpqzliMtbKzQ4igtslumFKKYfw2bBhymo6GTsAbpdoRa+UQ9hxcXCrhXXrpmd5pfVLo3TDlFLO4AtYf+ERq4X16Hw2rKOH4Ay9VvRKOUMgYP0RCFYL60TfEWrdWL7Zwa09eqWcIhJ69P2+wtRIdNGUbHJTYhGLX421olfKOSKhR29roh/spQRn5KcwI39gE7mD4RLRdfRKOYTPhguPWC2sV93YRSt6pZzBGEPA4PiKPqx79HbRHr1SztD1zlx3xqpjeFwureiVcoCuv2M91EwdwyV61o1STtD1d6w9enWMYEWvG6aUCnfdFb1umFJHc7t01Y1STqAVveqTx62JXikn6J6M1UQ/fAZ7eqVd9FAzpZxBE70FwnUdvVs3TCnlCF1zbZro1TG0olfKGbRHr/rk0Q1TSjmCT1s3qi9u3TCllCMEuit6Z6dCZ49umLh1w5RSjtBT0dscyDBz+PCGh1b0SjmDXzdMqb549OLgSjmCTydjVV/cbtEjEJRyAL8ur1R90XX0SjlD6LLTWtGrY+mFR5RyBt0wpfrk1h69Uo6gRyBYIFzPuvG4taJXygl0w5QFwvWsG704uFLO4PfrhinVB+3RK+UMfqMVveqD2+XSil4pB+g+1EyvGauOphceUcoZut6Zu0QTvTqK9uiVcoauDVO6jl4dI9ij152xSoU7n1979KoPbpcQMGCMVvVKhbOA0R696kPX2zxt3ygV3nQdveqTO/Tqr0sslQpv3TtjdTJWHa3rl0IreqXCm083TKm+JMR4AKhv7bQ5EqXUYHSGjq90a49eHW1KTjIAW8vD64wepdSRdlc3kxYfRUK02+5QhpUm+gGYmpOMS2DLwQa7Q1FKDcKmsnpm5qci2qNXR4uLdlOclagVvVJhrKXDR0llI7MKUu0OZdgNeaIXkXEicr+IPDnUjz2STM9NYctBTfRKhast5Q0EDMzKD6/TcweiX4leRB4QkSoR2XLU7ZeIyA4RKRWROwGMMbuNMbcMR7AjybS8FCob2qlqbLM7FKXUAGw6UA/AzPxUW+OwQn8r+geBS3rfICJu4B7gUmAqsExEpg5pdCPY9NzghOymA1rVKxWONh6oJy81jsykGLtDGXb9SvTGmDeB2qNung+Uhir4DuBxYOkQxzdizSpIJT7azRslVXaHopQ6RZ3+AG/trGbhuHS7Q7HEYHr0ecCBXp+XAXkiki4ifwBmi8g3+/pmEVkuIutEZF11dfUgwrBHbJSbs4ozeG17tZ55o1SYWbe3joY2HxdNzbI7FEsM+WSsMeawMeZzxpjxxpgfnuDrVhhj5hlj5mVmZg51GJY4f3IW5fWt7KhstDsUpdQpeHV7JdFuF2dPCM/cc6oGk+jLgYJen+eHbosY500OVgOvbNP2jVLhwhjDK9uqWDg+vXuXu9MNJtG/D0wQkSIRiQauA1aeygOIyBIRWeH1hueEZnZyLHPHpvHsxnJt3ygVJrYebGB3TTMXT8u2OxTL9Hd55WPAO8AkESkTkVuMMT7gDuAFYBvwhDFm66n8cGPMc8aY5Skp4buO9RNz8ympbGKzbp5SKiw8vaGcaLeLK2bk2h2KZfq76maZMSbHGBNljMk3xtwfun2VMWZiqB//f8Mb6sh0+cwcYjwu/rauzO5QlFIn4fMHWLnpIOdPziIlPsrucCxj6xEI4d66AUiOjeKS6aNZuekgbZ1+u8NRSp3AW6U11DS1c/WcPLtDsZStid4JrRsItm+8rZ06KavUCPf3D8pJjY/ivEmRsayyix5qNgQWjc8gJyWWJ9cfOPkXK6Vs0dTu44WtFVw+I4doT2Slvsga7TBxu4Rr5uTxRkk1B+tb7Q5HKXUcqzYfoq0zwDUR1rYB7dEPmWXzxwDwyLv7bI5EKXU0YwwPrN7DxOxE5oxJszscy2mPfojkp8Vz8bTRPLp2P60dOimr1Ejy+o5qtlc0cts54x1/kZHj0dbNELrpzCK8rZ08/YEutVRqJPn9G7vITYnlytMiZ+18b5roh9DphWlMz0vmT2v26k5ZpUaIDfvreG9PLbecPY4od2SmPO3RDyER4aZFRZRWNfHWzhq7w1FKAX94fRcpcVFcd3rByb/YobRHP8SumJVDZlIM97xWqlW9UjYrrWripW2V3HjG2Ig5wOx4IvN9zDCK8bj5wuLxrN1Ty5rSw3aHo1REW/HmLmI8Lm5cVGh3KLbSRD8Mli0YQ25KLD97cYdW9UrZ5JC3lWc+KOfaeQWkJzr/coEnool+GMR43HzpgglsPFCvxyIoZZOfvVCCINx69ji7Q7GdTsYOk4/PzWdsejw/f6mEQECreqWstKXcy9MflHHTWYUUjIq3Oxzb6WTsMIlyu/jqRRPZdqiBpz+IqAtvKWUrYwzf/8dHpMVH84Xziu0OZ0TQ1s0wWjIzl1kFqfz0he00t/vsDkepiLBqcwXv7anlPz82keTYyDlz/kQ00Q8jl0v4zhVTqWxo549v7LI7HKUcr63Tzw9WbWPy6CSuO32M3eGMGJroh9ncsWlcOSuXP765m3I92VKpYXVv6O/su0um4XZF3pk2fdFEb4FvXDoZgJ88v93mSJRyrgpvG797fReXTh/NGePT7Q5nRNFVNxbIS43jtnPG8ezGg2zYX2d3OEo50o+f347fGP77sil2hzLi6Kobi9x27niykmL4n5Vb8etyS6WG1Pp9dTzzQTm3nl2kyymPQ1s3FkmI8fCty6ewqczLw+/stTscpRwjEDB8/7mtZCXF8PnFupzyeDTRW+jKWbmcNymTn76wg7K6FrvDUcoRnlh3gE1lXu68dHJEH1x2IproLSQi3HX1DAC+9cwWPQdHqUEqrWrkf577iIXjRnHVaZF3Ldj+0kRvsbzUOP7r4km8UVLN39brlaiUGiifP8BX/rqR+Gg3v75uNi5dTtknTfQ2+PQZhSwoGsX3Vm5lb02z3eEoFZYefmcfW8ob+J+l08hOjrU7nBFNE70N3C7hl588jSi3iy8//gGd/oDdISkVViq8bfz8xR2cOzGTy2fk2B3OiKfr6G2SmxrHj66ZwaYyL798qcTucJQKK9//x1Z8AcP/Lp2OiLZsTkbX0dvo0hk5fHJeAb9/Yxfv7NKrUSl1MsYYfvvqTlZtruBLF0xgTLqume8Pbd3Y7DtLplKUnsB//HUjFd42u8NRakS7+5VSfvZiCVedlsvyc/SCIv2lid5mCTEefnP9bBrbOvnU/Wupa+6wOySlRqT1+2r59SslXD07r3uOS/WP/kuNANNyU3jgs6ez73Az3//HR3aHo9SIU9vcwZce20huahzfXzpN+/KnSBP9CLFgXDq3nzueZz4o5/Udep1Zpbr4/AG++NgGqpva+d0Nc0jSi4mcMk30I8jnzytmQlYit/95A2+UVNsdjlIjwk9f3MGa0sPcddV0Zuan2h1OWNJEP4LERrl59NaFFGUkcNsj69h6MPKWnSrV21/W7uOPb+zmhgVjuHZegd3hhC1N9CNMZlIMD908n9S4aG57ZD2VDboSR0WmR97dx7ee2cIFk7P47pJpdocT1jTRj0CZSTGs+Mxc6po7WLbiXaoaNdmryPLo2v38v79v4cIpWfzuU3OI9miqGgz91xuhZuan8tDN86loaOPmB9+nud1nd0hKWWJLuZfvrtzC4kmZ/O6GucR43HaHFPb0CIQRbF7hKH57/Ww+OtjATX96n/oWXWOvnG1zmZc7Ht3AqIRofnntaVrJDxE9AmGEO39yNr+6bjYbD9Tz8d+/zYFavWCJcqYn1h1g6T2raWr385tlc0hLiLY7JMfQl8swcOWsXB65ZT7Vje1c/bu32VnZaHdISg2pJ9eX8Y2nPuTM4gxe/dq5zC8aZXdIjqKJPkwsGJfO059fhEtg2b1rKa1qsjskpYbEY+/t5+tPbuLM8Rnc+5l5JOuGqCGniT6MFGcl8eitCwG4/t53NdmrsGaM4cE1e/jm05tZPDGT+26cR2yUTrwOB030YaY4K5HHbl2AP2C48rer+ev7++0OSalT5m3p5PY/b+B7z33ERVOz+cOn52qSH0aa6MPQhOwknvviWcwZk8Y3ntrMj/61nUBALzSuwkOFt41r//gOr2yv5BuXTOb3N8zRJZTDzGN3AGpgclPjeOjm+Xzn2S384Y1dVDa08cNrZmhVpEa00qombnzgPbytnTx003wWFWfYHVJE0EQfxtwu4a6rppOTEsvPXixh7e7DfH/pdC6cmm13aEodIRAwrHhrN796uYTEGA+PL1/I9DxdVm0Vbd2EORHhjvMn8Oi/LyAlPprb/ryeZz4owxht5aiRoandx1f+upEf/Ws7507MZOUdZ2mSt5gmeodYVJzB3z53BnPGpPIff93EJ/7wDiW63l7ZyBjDo2v3c85PXmPlpoN845LJ/OFTc8lNjbM7tIijid5BEmM8/OXfF/J/V09nb00zV/xmNfe9tVsnapXl9h9u4Yb71vLfz2xmQlYiz37hTG5fPF6vDGUT7dE7TLTHxQ0LxnLxtNHc+dRm7vrnNl7ZVsXPrp1FnlZSygIb9tfxmfvfA+AHV89g2fwCTfA204reoTISY7j3M3P58cdn8GFZPRf94g1++VIJ7T6/3aEphyqra+FXL5dw4wPvkZ4YzfNfOZvrF4zRJD8CaEXvYCLCJ08fw6LxGfzoX9v59Ss7eWFrBXcvm83E7CS7w1MO4fMH+MGq7Tz8zl78xnDGuHR+8omZ5KfF2x2aCtFEHwEKRsVzzw1zuGZbJf/15Idc8ZvV3HRmIdfOK2B8ZqLd4akw1dDWyUNr9rJmVw3v7q7l+gVjuOO8Yp1sHYGGvHUjIgki8pCI3CsiNwz146uBu2BKNs9/5RwunJLFfW/t4dJfvcUj7+6jSS9qok5RIGD4j8c38ouXS9hd3cz/Lp3GD66eoUl+hOpXRS8iDwBXAFXGmOm9br8E+DXgBu4zxvwIuAZ40hjznIj8FfjL0IetBiozKYbf3TCXqsY2/vOJTfy/v2/heyu38rWPTeL2xePtDk+NcIe8rfxg1XY+2F9HWV0r31sylc+eWWR3WOok+tu6eRD4LfBw1w0i4gbuAS4CyoD3RWQlkA9sDn2ZzvyNUFlJsTx403zWlNbw6Nr9/Pj57azfV8ul03O4anYebpdOoKkedc0d3PvWbh58ey8BY7hgSjafXVTIjYsK7Q5N9UO/Er0x5k0RKTzq5vlAqTFmN4CIPA4sJZj084GNnKA1JCLLgeUAY8aMOdW41RBwu4RzJmZyZnEGP39xB3//oJyXt1Xxp7f3cPOZRVwxM1cv5RbhvC2d3L96Nw+s2Utzh4/LZ+Tw9YsnMTY9we7Q1CkYzGRsHnCg1+dlwALgbuC3InI58Fxf32yMWQGsAJg3b57u6LGR2yX81yWT+frFk1i56SC/eKmErz6xid++Vsqdl0zmwinZuLTCjxid/gDv7amlurGdu/65jZqmdi6bMZovXzCRSaN1tVY4GvJVN8aYZuCmoX5cNfxEhKWn5bFkZi6v7ajif//xEcsfWc+4zARuPXsc18zJ0+NkHa6lw8cdj37Aq9urAJiYnciDN52uZ9OEucEk+nKgoNfn+aHb+k1ElgBLiouLBxGGGmoul3DBlGzOmZjJqs2HWPHmbr759GZ+/fJOzp+SxWXTczhrgh4v6xT+gOGJdQdYXVrDmzuqaerw8e3LpzAjL4VZBal69LUDDCbRvw9MEJEiggn+OuD6U3kAY8xzwHPz5s27dRBxqGES5Xax9LQ8rpyVy+rSGu57aw/PbTwYPKhqYibXnV7A+ZOzNBGEobZOP//48BDbDzWw5aCXd3fXkpsSy8XTR3PtvAK9OLfD9Hd55WPAYiBDRMqA7xpj7heRO4AXCC6vfMAYs3XYIlW2ERHOnpDJ2RMyaff5eWD1Xv60Zg+f/8sGkmI8FGcncv6kLD63eDxRbp28HamMMeypaWbtnlp++VIJVY3txHhcRLld/OQTM7l2XsHJH0SFpf6uulnWx+2rgFVDGpEa0WI8bm5fPJ7l54zj7V01rNp8iJ2VTfz8pRKe2lDGouIMrp6dx7yxaXrGyQhxuKmdd3fX8si7e3l3dy0As/JT+NUnT+OM8ekA+lw5nK1HIGiPPny5XT1VPsDzWyr4y9p9rAy1djISo/nYtGAbYPLoJG3vWGxPTTMvbq3gkLeNx9/fT1tngIzEaL512RQWjktnWm6yrqSKILYmeu3RO8cl00dzyfTRtHT4eH5LBa/vqOap9WU8unY/Hpdw3uQszp2YyewxqUzKTsKjLZ4h1dTu47XtVVQ2tPHWzhre3FmNMeASuHxmLrecVcS03GRtrUUoPdRMDan4aA/XzMnnmjn51DZ38PauGj7YX89zmw7y0keVAMRFuZlVkMJNZxYxZ0waafFRmvhPUVVDG3tqmjnkbWNNaQ3/3HyIlo7gRvS81Di+eF4xNywcS0ZijO5yVtq6UcNnVEI0V8zM5YqZuXz78ikcqG3lgwN1bDxQz6vbq7jtkfVA8Pyda2bnMSE7icL0eKbkJJMQozVIb+X1rby6vYpKbxuVDW08uaGMrssCJ0S7uXJWLp+Ym09xViIpcVHac1dH0NaNsoSIMCY9njHp8Sw9LY9vXTaFl7dVUtnQzqvbq7hv9R78oUseelzC9LwUFowbxcKidOYWpuH3GzxuISk2yuaRDK/mdh9lda0crG+lpLKRHZWNVDe2s7q0BmNABNwifHZRIedPziI3NY681DidA1EnpGWTsoXH7eKS6TkA3LiokE5/gLK6VnZXN7F+Xx3v7anlgdV7+OMbu4/4vsmjk7hgShZTcpIpTE9gbHo8iTEe2n2BsEl2h7ytbC7zsqg4gz3VzeysauSQt41NB+p5ZXtV9wsewOjkWJLjPHx+8Xg+PiefsekJdPgCxEWHx1jVyKCJXo0IUW4XRRkJFGUkcMGUbABaO/x8cKCODfvqiI1y09rhZ82uGn7/+i56X+88yi10+g0fm5rNhOxEjIHEWA+pcdGcOymTvNQ4Siob2VPTzNkTMoiPHtiv/YHaFnJT447oedc1d5Aa39Mq6fAF2He4mV3VzbT7/KTFR+M3hn9+eIjWDj81Te1s2F9Hp9/gEo4YR3ZyDDefWcjM/FRyU+MoTI8nPTHmmDg0yatTpYlejVhx0W4Wjc9g0fie4xa+eMEEWjp87K1pYd/hZvYebqG+tQO/3/C39WW8EjqjpXdVHB/t7p6ojItyMyM/hfSEaPbUNFPb3MFnzhjL5nIv2ysaiXK7yE+Lo7yuFb8xnDcpiy+eX8zKTQf5zrNbyU2J5d/mFRAT5eKp9WXsqm4mPy2O/LQ4qhvb2Xe4BV/g2DP6UuKiyEiMJjU+mk8vLOTsiRm8XVrD5NHJnDYmlZyU2AG/ACl1MmKMfQdH9pqMvXXnzp22xaGcIRAwdM1BtvsClNe38uLWSg43tZOfFkdxVhIvb6vkw7J6Gtt8ZCbF0OkP8P7eOjKTYlg4Lp2Wdh8HvW0UpMXhDxhe3VGFEKy8zyxOxyXS3S9fUDSKs4oz2FTmpaG1k7SEKIqzEoP/ZSYRH+OmqqGd5nYfZ0/M0APhVLfH39vPnU9v5p1vnk9OysCvyiUi640x8072dToZqxyj9wag2Cg34zMTuX3xkdfEPfowNmMMu6qbGZsef9w15tsONfCvLRX4/AG+dMEEYqPclNe30ukLUJhx8jPZ9Zq8aiTQ94oqookIxVl9J+MpOclMyUk+4rY8vS6qCjO6S0UppRxOE71SSjmcrYleRJaIyAqv12tnGEop5Wi2JnpjzHPGmOUpKXqZMqWUGi7aulFKKYfTRK+UUg6niV4ppRxOE71SStnEqoMJbD0CoTsIkWpg3wC/PQOoGcJwwomOPfJE6rghcsd+onGPNcZknuwBRkSiHwwRWdefsx6cSMceeWOP1HFD5I59KMatrRullHI4TfRKKeVwTkj0K+wOwEY69sgTqeOGyB37oMcd9j16pZRSJ+aEil4ppdQJaKJXSimHC+tELyKXiMgOESkVkTvtjmc4icheEdksIhtFZF3otlEi8pKI7Az9P83uOIeCiDwgIlUisqXXbccdqwTdHfod+FBE5tgX+eD1MfbviUh56LnfKCKX9brvm6Gx7xCRi+2JevBEpEBEXhORj0Rkq4h8OXS745/3E4x96J53Y0xY/ge4gV3AOCAa2ARMtTuuYRzvXiDjqNt+AtwZ+vhO4Md2xzlEYz0HmANsOdlYgcuAfwECLATW2h3/MIz9e8DXjvO1U0O/9zFAUejvwW33GAY47hxgTujjJKAkND7HP+8nGPuQPe/hXNHPB0qNMbuNMR3A48BSm2Oy2lLgodDHDwFX2RfK0DHGvAnUHnVzX2NdCjxsgt4FUkUkx5JAh0EfY+/LUuBxY0y7MWYPUErw7yLsGGMOGWM2hD5uBLYBeUTA836CsffllJ/3cE70ecCBXp+XceJ/nHBngBdFZL2ILA/dlm2MORT6uALItic0S/Q11kj5Pbgj1KJ4oFeLzpFjF5FCYDawlgh73o8aOwzR8x7OiT7SnGWMmQNcCnxBRM7pfacJvqeLiLWykTTWkN8D44HTgEPAz22NZhiJSCLwFPAVY0xD7/uc/rwfZ+xD9ryHc6IvBwp6fZ4fus2RjDHlof9XAc8QfKtW2fV2NfT/KvsiHHZ9jdXxvwfGmEpjjN8YEwDupedtuqPGLiJRBBPdX4wxT4dujojn/XhjH8rnPZwT/fvABBEpEpFo4Dpgpc0xDQsRSRCRpK6PgY8BWwiO98bQl90IPGtPhJboa6wrgc+EVmEsBLy93uo7wlG956sJPvcQHPt1IhIjIkXABOA9q+MbCiIiwP3ANmPML3rd5fjnva+xD+nzbveM8yBnqy8jOEO9C/iW3fEM4zjHEZxl3wRs7RorkA68AuwEXgZG2R3rEI33MYJvVTsJ9h9v6WusBFdd3BP6HdgMzLM7/mEY+yOhsX0Y+iPP6fX13wqNfQdwqd3xD2LcZxFsy3wIbAz9d1kkPO8nGPuQPe96BIJSSjlcOLdulFJK9YMmeqWUcjhN9Eop5XCa6JVSyuE00SullMNpoldKKYfTRK+UUg73/wGwun0hKB1OBAAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 1 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"plt.plot(m_aia.data.mean(axis=0))\n",
"plt.yscale('log')"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python [conda env:million-degree-loops]",
"language": "python",
"name": "conda-env-million-degree-loops-py"
},
"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.8.5"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Display the source blob
Display the rendered blob
Raw
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment