UGCA442.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "bb1552ac-8155-47f7-a3ce-4372b5d41629",
   "metadata": {},
   "outputs": [],
   "source": [
    "from src.tools import *\n",
    "from src.main import *"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "12b4dc46-0b54-4f6f-9578-b6d250ef2582",
   "metadata": {},
   "outputs": [],
   "source": [
    "massmodels = import_massmodel_data('./sparc/MassModels_Lelli2016c.mrt')\n",
    "mm_dict = {mm['ID']: massmodels[massmodels['ID'] == mm['ID']] for mm in massmodels}\n",
    "metadata, refs = import_galaxy_sample_data('./sparc/SPARC_Lelli2016c.mrt')\n",
    "btfr_data = import_btfr_data('./sparc/BTFR_Lelli2019.mrt')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "57ee2c5b-3e19-4cb5-988e-024e7ca6fa5a",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "UGCA442\n",
      "N = 26\n"
     ]
    }
   ],
   "source": [
    "# Specify Model Galaxy\n",
    "gal_id = 'UGCA442'\n",
    "n_rows,_ = get_single_galaxy_N(gal_id,datapath='./sparc/sfb_LTG/')\n",
    "print(str(gal_id)+'\\nN = '+str(n_rows))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "10304822-566e-474c-95ee-2764083df5b4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "ID\tT\tD\te_D\tf_D\tInc\te_Inc\tL\te_L\tReff\tSBeff\tRdisk\tSBdisk\tMHI\tRHI\tVflat\te_Vflat\tQ\tRefs\n",
      "UGCA442\tSm\t4.35\t0.22\t2\t64.0\t7.0\t0.14\t0.005\t1.71\t7.6\t1.18\t116.1\t0.263\t4.37\t56.4\t2.1\t1\t['Co00']\n",
      "(1) Cote, S. et al. 2000, AJ, 120, 3027\n"
     ]
    }
   ],
   "source": [
    "# Print the overview for this galaxy\n",
    "print('\\t'.join(['ID'] + [d_name for d_name in metadata[gal_id].dtype.names]))\n",
    "print('\\t'.join([gal_id] + [str(item) if name != 'T' else str(hubble_types[item]) for item,name in zip(metadata[gal_id],metadata[gal_id].dtype.names)]))\n",
    "print('\\n'.join(['({}) '.format(cit_idx+1)+refs[cit_ref] for cit_ref,cit_idx in zip(metadata[gal_id]['Refs'],range(len(metadata[gal_id]['Refs'])))]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "6938a7a4-f740-4f78-8588-57fb5ccbda9a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX8AAAEWCAYAAACOv5f1AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAeeElEQVR4nO3de5gcVZ3/8fcnJAGEJECIMdxMVkHlwRjYEY0gDAzwWyAKqy7gBSPLiroqqKCie5FF3Z+3FXDxFgkaERAMBsJlkezIgOgIJIAhEBSEEIIJCZgLESSEfPePc4Z0mrn0TKamp6c+r+fpZ7pO3b7VM/Ot06dOnVJEYGZm5TKs3gGYmdnAc/I3MyshJ38zsxJy8jczKyEnfzOzEnLyNzMrISd/M7MScvK3wkjaS9J6SdvUO5ZGI+lHkr5U47JLJB3Rxbw2SX+VdGv/Rji45c/vWUnL6h3LYOXk3yDyP/hKSTtUlP2TpLZ+3k+bpH/qj21FxNKI2DEiXuiP7VXKcYakN1SVz8nlzf29zy7i+EDe33lV5cfl8h8NRBw9+FhEHNIxkeN6deUCks6R9JOK6VGSvpn/7v4iaamk2ZLeVLGMJJ0uaVFeZpmkn0l6fSfbjsp1q0m6uLO48ry98wmsMr5jJd0maY2kFZIukjSqY35EfAA4uvaPqHyc/BvLNsAZ9Q5iEPkD8P6OCUljganAqgGO44/ACZKGV5RNJ8XXcCRtC/wSeD0wDRgNvA74KVsm1AtIf4+nA7sA+wBXA8dWbEuk39GfqfhdVe3vYOBV3YT0beDOqrIxwJeA3XJsuwNfr+HwLHPybyxfB86StFNnMyW9VtI8SX+W9HtJJ+TySbmGNCxP/0DSyor1LpH0CUlfBt4KXJibay7M898i6U5Ja/PPt1Ss2ybpi5J+LelpSTdJ2jXPm5hrc8Pz9CmSFuflHpb0oYrtNOea45n5G85ySaf08HlcCpxY0az0bmAOsKFiuwdKas/Hv1zShZJG5nmSdF7e3zpJ90raL887RtL9OdbHJZ3VTRwrgHuB/5fX3QV4CzC36vfzdkn35VjaJL2uYt7+ku7K+7sC2K5q3WmS7snr/kbS5B4+m61xMrAHcHxELIqIFyLiLxExOyLOyfHsDXwUeHdE/DIinouIZyLi0oj4SsW23gpMIJ0gTur47CuOazjw38DHOwtE0knAGqC1sjwiLouIG/M+VwM/AA7a+kMvDyf/xjIfaANekoiUmoPmAZcBLwdOAr4jad+IeARYB+yfFz8EWF+RfA4FbomIfwF+RWom2DEiPpYT2fXAt4CxwDeB63Mtu8N7gFPyfkd2Fl+2ks01yVOA8yQdUDH/FaQa3e7AqcC3Je3czefxJ+B+4Kg8/X7gx1XLvAB8EtiV9K2gBfjnPO+o/Fnsk/d7AvBUnjcT+FBEjAL2I9WEu/NjNtdsTwKuAZ7rmClpH+By4BPAOOAG4FpJI3NCvBq4hFSD/hnwzop19wcuBj5E+h18H5iba+hFOAL4RUT8pZtlWoBlEXFHD9uaDlwLXJmn31Y1/5PArRGxsHpFSaOBc4FP1RDzIcB9NSxnmZN/4/l34OOSxlWVTwOWRMQPI2JjRNwNXAX8Q55/C3CopFfk6dl5ehIpGf+ui/0dCzwYEZfk7V4OPMCW/8Q/jIg/RMSzpH/yKZ1tKCKuj4g/RnILcBOpZtjheeDciHg+Im4A1gOv6eHz+DHwfkmvBXaKiPaqfS6IiN/m2JeQEuehFfsbBbwWUEQsjojlFfP2lTQ6IlZHxF09xDEHaJY0hs5PQicC10fEvIh4HvgGsD3pG8KbgRHA+fnYZ7NlM8dpwPcj4vZcC59FOrG8uYeY+mpX0rcZACRNyd841kn6fS4eCyzvdO3N672M9Pd3WT7m2WzZTLcn6YT2711s4ovAzIjo9qKtpCNJJ5mutmOdcPJvMBGxCLgOOLtq1iuBN+V/0jWS1gDvJdWmISX/ZlIN6VbSN4hD8+tXEbGpi13uBjxaVfYoqXbeYUXF+2eAHTvbkKSjJf02N0utAY4hJZoOT0XExlq2VeHnwOHAx0g15+p97iPpOqWLguuA/+zYZ0T8EriQ1Ka8UtKMXNuEVPM+BnhU0i2SpnYXRD7xXQ/8KzA2In5dtcgWn2P+vB8jfY67AY/HlkPsVn7mrwTOrPrd7pnX64sXSCebSiNIJzxI334mVMR6T0TsBLwD2LazZbrw98BG0rccSM10R1dUXM4nnezXVq8oaQrpG8h51fOqlnsz6dvuuyKiIa+x1IuTf2P6AvBBtkzAj5GabnaqeO0YER/J828h1bKb8/vbSG2kh+bpDtVjfP+JlHwq7QU83puAcxPFVaQa7/icTG4A1JvtVIuIZ4D/AT5CJ8kf+C7pm8reETEa+HzlPiPiWxHxt8C+pOafT+fyOyPiOFJT1tVsbrbozo+BM4GfdDJvi89RkkgJ/HFSDXr3XNZhr4r3jwFfrvrdvix/C+uLpcDEqrJJbD7htAJHqaJnWSdagT0kNXWzzHTSyXuppBWk5qwRpGZCSE1HX88n5o4KRLuk95D+TidWrHsW8E5JL34Dy81hc4F/jIgtrglYz5z8G1BEPARcQbqI1uE6YB9JJ0sakV9v7GjXj4gHgWeB95FOEuuAJ0g13Mrk/wTwNxXTN+TtvkfScEknkhLldb0MeySp1rgK2CjpaDa31W+tzwOH5madaqNI1zvW56ahjpMh+fN5k6QRwF+AvwKbcjv8eyWNyc0V64CuvhlVugU4knQBs9qVwLGSWvL+ziQ13fwGaCfVkE/Pv7d3AAdWrPsD4MM5VknaQamr46jqndToCuBfJe0haZjSPQJvIzXLQDqJLQfmSNpP0jaStgNeTPT57+k7wOVKF+tHStpO0kmSzpa0Oym5TyM1A04B3gB8lc1NP/vkso755DjmADNIPYA65n2P9M2q46L6fsCNwMcj4to+fg6l5uTfuM4FXqyZRcTTpGR6EqmWuYL0j1Z5UfAWUtPKYxXTAirbsy8A3iVptaRvRcRTpH/gM0lf9T8DTIuIJ3sTbI7vdFISXE2q/c3tdqXat/2niLiti9ln5X09TUqiV1TMG53LVpNqvU+xubvgycCS3FT0YVITWk9xRES0RsSfO5n3e9KJ97+BJ0lJ7m0RsSEiNpCaVD5A6hJ5Iqk5q2Pd+aRvehfmWB/Ky/bVuaSTzm15e18D3pubFImIvwKHkS6mX086+f0eeCPponiH09ncbLaG1OX170kXeE8G7omImyJiRceL1HFgsqT9ImJl1TyAJyPi2dyLp3LeeuCvEdHRjfdM0oXzmUo909ZL8gXfXlD4SV5mQ5akm0i9nOZHxGH1jmegSJpJuti8MiJecuOYOfmbmZWSm33MzErIyd/MrISG97zI4LDrrrvGxIkT6x2GmVlDWbBgwZMRUX1TaOMk/4kTJzJ//vx6h2Fm1lAkVd+kCbjZx8yslJz8zcxKyMnfzKyEnPzNzErIyd/MrISc/M3MSsjJ38yshJz8zcxKqPDkL2knSbMlPaD08O6pknZRetD4g/lnd89pNTPrP83N6VVyA1HzvwC4MSJeS3pww2LSIwhbI2Jv0hOBqh9JaGZmBSo0+eeHWR8CzATID65YAxwHzMqLzQKOLzIOMzPbUtE1/0mkx/b9UNLdki7KzwUdHxHL8zIrgPGdrSzpNEnzJc1ftWpVZ4uYmVkfFJ38hwMHAN+NiP1Jz0ndookn0tNkOn2iTETMiIimiGgaN+4lg9KZmVkfFZ38lwHLIuL2PD2bdDJ4QtIEgPxzZcFxmFmRfBG14RSa/PODlx+T9Jpc1EJ6KPRcYHoumw5cU2QcZma2pYEYz//jwKWSRgIPA6eQTjpXSjoVeBQ4YQDiMDOzrPDkHxH3AE2dzGopet9mZtY53+FrZlZCTv5mVi5r18LSpdDeXu9I6srJ38zKo70dFi6ERx6BlpZSnwCc/M2sPNraYNOm9H7DhjRdUk7+ZlYezc0wLKe9kSNLfW+Ck7+ZlcfUqTB5MkyaBK2taXowK/DmOSd/MyuXMWNgr70Gf+KHQi9OO/mbDUaNNlyCe9D0v4IvTjv5m9nWcQ+aYhR8cdrJ38y2jnvQFKPgi9NO/ma2ddyDphgFX5weiIHdzGwo60hSa9fCpZc2xoXURnH33YVt2snfzLbemDHp5cTfMNzsYzYYufeMFcw1f7PBpqP3zKZNqfdMI9yM1Eh8QRpwzd9s8HHvGRsATv5mg417z9gAcPI3G2wabfwZa0hu8zcbjNx7xgrmmr+ZWQk5+ZuZlZCbfcxs67lHUsNxzd/MrIRc8zcbjFyTtoK55m9mVkJO/mZmJeTkb2ZWQk7+ZmYl5ORvZlZCTv5mZiXk5G9mVkJO/lYezc0eHtksK/wmL0lLgKeBF4CNEdEkaRfgCmAisAQ4ISJWFx2LmZklA1XzPywipkREU54+G2iNiL2B1jxtZmYDpF7NPscBs/L7WcDxdYrDzKyUBiL5B3CTpAWSTstl4yNieX6/Ahjf2YqSTpM0X9L8VatWDUCoZmblMBADux0cEY9LejkwT9IDlTMjIiRFZytGxAxgBkBTU1Ony5iZWe8VXvOPiMfzz5XAHOBA4AlJEwDyz5VFx2HG2rWwdCm0t9c7ErO6KzT5S9pB0qiO98BRwCJgLjA9LzYduKbIOMxob4eFC+GRR6ClxScAK72im33GA3Mkdezrsoi4UdKdwJWSTgUeBU4oOA4ru7Y22LQpvd+wIU374ehWYoUm/4h4GHhDJ+VPAS1F7ttsC83NMGxYOgGMHOmbvaz0fIevlcPUqTB5MkyaBK2trvVb6fkxjlYeY8aklxO/mWv+ZmZl5ORvZlZCTv5mZiXk5G9mVkK+4Gvl0dZW7wjMBg3X/M3MSsjJ38yshJz8zcxKyMnfzKyEnPzNzErIyd/MrISc/M3MSsjJ38yshJz8zcxKyMnfzKyEnPzNzErIyd/MrIRqGthN0suBg4DdgGeBRcD8iNhUYGxmZlaQbpO/pMOAs4FdgLuBlcB2wPHAqyTNBv4rItYVHKeZmfWjnmr+xwAfjIil1TMkDQemAUcCVxUQm5mZFaTb5B8Rn+5m3kbg6v4OyMzMilfTBV9JZ0garWSmpLskHVV0cDbINTenl5k1nFp7+/xjbtc/CtgZOBn4SmFRmZlZoWpN/so/jwEuiYj7KsrMzKzB1Jr8F0i6iZT8fyFpFOBunmZmDaqnrp4jIuJ54FRgCvBwRDwjaSxwygDEZ2ZmBeipq2e7pGXAjcCNEbEGICKeAp4qODYzMytIt80+EdEEfCJPni/pTknnSTpK0raFR2eD29q1sHQptLfXOxIz66Ue2/wjYklEfC8ijgfeAlwLHAH8StL1Bcdng1V7OyxcCI88Ai0tPgGYNZheDewWEc9HxC8j4jMRcSBwWi3rSdpG0t2SrsvTkyTdLukhSVdIGtmH2K2e2tpgU77mv2FDmjazhlHrTV7TcvJeLWmdpKclrYuIx2vczxnA4orprwLnRcSrgdWkC8rWSJqbYVj+8xk50jd7mTWYWmv+5wPTgV0iYnREjIqI0bWsKGkP4Fjgojwt4HBgdl5kFmmgOGskU6fC5MkwaRK0tqZpM2sYNQ3pDDwGLIqI6MM+zgc+A4zK02OBNXlsIIBlwO592K7V25gx6eXEb9Zwak3+nwFukHQL8FxHYUR8s7uVJE0DVkbEAknNvQ1O0mnk6wp77bVXb1c3M7Mu1Jr8vwysJ43l35uLswcBb5d0TF53NHABsJOk4bn2vwfQ6bWDiJgBzABoamrqy7cOMzPrRK3Jf7eI2K+3G4+IzwGfA8g1/7Mi4r2Sfga8C/gp6VrCNb3dtpmZ9V2tF3xv6OchnD8LfErSQ6RrADP7cds2UNra3MXTrEHVWvP/CHCWpOeA50kjekatPX5IC7cBbfn9w8CBvYrUzMz6TU3JPyJGVZflLptmZtaAar3J69yq6WHATwqJyMzMCldrm/+ekjou3G4LzAEeLCwqMzMrVM2PcQRen08A1wI3R8Q5hUVlZmaF6ulhLgdUTF4AfB/4NXCrpAMi4q4igzMzs2L0dMH3v6qmVwP75vIgjdFjZmYNptvkHxGHDVQgZmY2cLpt85f0vu66dEp6laSD+z8sMzMrUk/NPmOBeyQtABYAq0hj9LwaOBR4Eji70AjNzKzf9dTsc4GkC0lt+wcBk4FnSQ9mOTkilhYfopmZ9bce7/CNiBeAefllZmZDQK+e4WtmZkODk7+ZWQk5+ZuZlVCtA7udIWm0kpmS7urn8f3NzGwA1Ty2T0SsA44CdgZOBr5SWFRmZlaoWpN/x41exwCXRMR9FWVmZtZgak3+CyTdREr+v5A0CthUXFhmZlakWh/jeCowBXg4Ip6RNBY4pbCozMysULU+xnGTpInA+yQFcFtEzCk0MjMzK0ytvX2+A3wYuBdYBHxI0reLDMzMzIpTa7PP4cDrIiIAJM0C7i8sKjMzK1StF3wfAvaqmN4TP8PXzKxh9fQYx2tJT+waBSyWdEeedSBwR5crmpnZoNZTs883BiQKMzMbUD2N539Lx3tJ44E35sk7ImJlkYGZmVlxau3tcwKpmecfgBOA2yW9q8jAzMysOLX29vkX4I0dtX1J44D/BWYXFZiZmRWn1t4+w6qaeZ7qxbpmZjbI1Frzv1HSL4DL8/SJwA3FhGRmZkWrdXiHT0t6B3BwLprh4R3MzBpXrTV/IuLnwM87myepPSKm9ltUZmZWqP5qt9+us0JJ20m6Q9LvJN0n6T9y+SRJt0t6SNIVkkb2UxxmZlaD/kr+0UX5c8DhEfEG0pDQfyfpzcBXgfMi4tXAatKQ0WZmNkAK7bETyfo8OSK/gjRQXEc30VnA8UXG0TCam9PLzKxg/ZX8u3yko6RtJN0DrATmAX8E1kTExrzIMmD3LtY9TdJ8SfNXrVrVT6GamVm3yV/StyUdVMN2Tu5qRkS8EBFTgD1IA8K9ttbgImJGRDRFRNO4ceNqXc3MzHrQU83/D8A3JC2R9DVJ+3e2UEQs6mlHEbEGuBmYCuwkqaOn0R7A47WHbGZmW6vb5B8RF+QunIeS7uq9WNIDkr4gaZ+eNi5pnKSd8vvtgSOBxaSTQMfYQNOBa/p+CGZm1ls1tflHxKMR8dWI2B94N+kC7eIaVp0A3CxpIXAnMC8irgM+C3xK0kPAWGBmX4IfctauhaVLob293pGY2RBX001euYnmaOAkoAVoA87pab2IWAi8pKkoIh4mtf9bh/Z2WLgQNm2ClhZobYWpvm/OzIrR0wXfIyVdTOqR80HgeuBVEXFSRLippj+1taXED7BhQ5o2MytITzX/zwGXAWdGxOoBiKe8mpth2LB0Ahg50v39zaxQPT3J6/CBCqT0pk6FyZNTu/+ll7rJx8wKVfPAbjYAxoxJLyd+MyuYH8hiZlZCTv5mZiXk5G9mVkJO/mZmJeQLvoOJ+/ab2QBxzd/MrISc/M3MSsjJ38yshJz8zcxKyMnfzKyEnPzNzErIyd/MrISc/M3MSsjJ38yshJz8zcxKyMnfzKyEnPzNzErIyd/MrISc/M3MSsjJ38yshJz8zcxKyMnfzKyEnPzNzErIyd/MrISc/M3MSsjJ38yshJz8zcxKyMnfzKyECk3+kvaUdLOk+yXdJ+mMXL6LpHmSHsw/dy4yDjMz21LRNf+NwJkRsS/wZuCjkvYFzgZaI2JvoDVPm5nZACk0+UfE8oi4K79/GlgM7A4cB8zKi80Cji8yDjMz29KAtflLmgjsD9wOjI+I5XnWCmB8F+ucJmm+pPmrVq0amEDNzEpgQJK/pB2Bq4BPRMS6ynkREUB0tl5EzIiIpohoGjdu3ABEamZWDoUnf0kjSIn/0oj4eS5+QtKEPH8CsLLoOMzMbLOie/sImAksjohvVsyaC0zP76cD1xQZh5mZbWl4wds/CDgZuFfSPbns88BXgCslnQo8CpxQcBxmZlah0OQfEbcB6mJ2S5H7NjOzrvkOXzOzEnLyNzMrISd/M7MScvI3MyshJ38zsxJy8jczKyEnfzOzEnLyNzMrISd/M7MScvI3MyuhoZ/8m5vTy8zMXjT0k7+Zmb2Ek7+ZWQk5+ZuZlZCTv5lZCTn5m5mV0NBP/mvXwtKl0N5e70jMzAaNoZ3829th4UJ45BFoafEJwMwsG9rJv60NNm1K7zdsSNNmZjbEk39zMwzLhzhypG/2MjPLhnbynzoVJk+GSZOgtTVNm5kZw+sdQOHGjEkvJ34zsxcN7Zq/mZl1ysnfzKyEnPzNzErIyd/MrISG/gVf9+03M3sJ1/zNzErIyd/MrISc/M3MSsjJ38yshJz8zcxKqNDkL+liSSslLaoo20XSPEkP5p87FxmDmZm9VNE1/x8Bf1dVdjbQGhF7A6152szMBlChyT8ibgX+XFV8HDArv58FHF9kDGZm9lL1aPMfHxHL8/sVwPiuFpR0mqT5kuavWrVqYKIzMyuBut7hGxEhKbqZPwOYASBplaRHByy4/rEr8GS9g+gnQ+lYYGgdz1A6FhhaxzMYjuWVnRXWI/k/IWlCRCyXNAFYWctKETGu4Lj6naT5EdFU7zj6w1A6FhhaxzOUjgWG1vEM5mOpR7PPXGB6fj8duKYOMZiZlVrRXT0vB9qB10haJulU4CvAkZIeBI7I02ZmNoAKbfaJiHd3MaulyP0OIjPqHUA/GkrHAkPreIbSscDQOp5BeyyK6PJ6q5mZDVEe3sHMrISc/M3MSsjJvwCdjWnUqCTtKelmSfdLuk/SGfWOqa8kbSfpDkm/y8fyH/WOqT9I2kbS3ZKuq3csW0PSEkn3SrpH0vx6x7O1JO0kabakByQtljS13jFVcpt/ASQdAqwHfhwR+9U7nq2R78WYEBF3SRoFLACOj4j76xxar0kSsENErJc0ArgNOCMiflvn0LaKpE8BTcDoiJhW73j6StISoCki6n1TVL+QNAv4VURcJGkk8LKIWFPnsF7kmn8BuhjTqCFFxPKIuCu/fxpYDOxe36j6JpL1eXJEfjV07UfSHsCxwEX1jsU2kzQGOASYCRARGwZT4gcnf+sFSROB/YHb6xxKn+UmkntId5bPi4iGPZbsfOAzwKY6x9EfArhJ0gJJp9U7mK00CVgF/DA3yV0kaYd6B1XJyd9qImlH4CrgExGxrt7x9FVEvBARU4A9gAMlNWyznKRpwMqIWFDvWPrJwRFxAHA08NHcfNqohgMHAN+NiP2BvzDIhq938rce5fbxq4BLI+Ln9Y6nP+Sv4Dfz0udNNJKDgLfntvKfAodL+kl9Q+q7iHg8/1wJzAEOrG9EW2UZsKzim+Vs0slg0HDyt27li6QzgcUR8c16x7M1JI2TtFN+vz1wJPBAXYPaChHxuYjYIyImAicBv4yI99U5rD6RtEPuUEBuHjkKaNjechGxAnhM0mtyUQswqDpJ1HVI56Eqj2nUDOwqaRnwhYiYWd+o+uwg4GTg3txWDvD5iLihfiH12QRglqRtSBWfKyOiobtHDiHjgTmprsFw4LKIuLG+IW21jwOX5p4+DwOn1DmeLbirp5lZCbnZx8yshJz8zcxKyMnfzKyEnPzNzErIyd/MrISc/G3IkPRCHhFykaRrO/r092L9NklN+f0NvV2/i23uL2lmfn+OpLP6YZv/K2nnrd2OlZuTvw0lz0bElDyS6p+Bj/Z1QxFxTD8NxPV54Fv9sJ1KlwD/3M/btJJx8rehqp08+qikAyW15wG2ftNx16Wk7SX9NI+1PgfYvmPlPLb8rpImVj6XQdJZks7J70/PzzlYKOmn1QHkO1YnR8TvOpn3QUn/k2Nok3RBxbeWA/MyO0r6YR7jfqGkd+bV5wJdPR/brCa+w9eGnHwHbwt5OF3SEA5vjYiNko4A/hN4J/AR4JmIeJ2kycBdvdzV2cCkiHiuiyaiJjoZokDSx0hDSxyf14U01vuUPJjZxcB+wL8BayPi9Xm9nQEiYrWkbSWNjYinehmzGeDkb0PL9nkIit1Jzx2Yl8vHkIZ12Js0bPCIXH4IuUkmIhZKWtjL/S0k3b5/NXB1J/MnkIb1rfR+4DFS4n++ovzyHMetkkbnk8kRpDF7yPNWVyy/EtgNcPK3PnGzjw0lz+bhml8JiM1t/l8Ebs7XAt4GbNeLbW5ky/+TynWPBb5NGq3xTknVlalnO9nXvcBE0pDSlarHWelp3JXt8vbN+sTJ34aciHgGOB04MyfkMcDjefYHKha9FXgPQB7Xf3Inm3sCeLmksZK2Babl5YcBe0bEzcBn8z52rFp3MfDqqrK7gQ8BcyXtVlF+Yt7uwaSmnrWkby4vXrTuaPbJI62+AljS3edg1h0nfxuSIuJuUrPMu4GvAf9f0t1s2dT5XWBHSYuBc0nPJ67ezvN53h2kZNwxBPQ2wE8k3UtK6N+q7h0UEQ8AYzqGKq4ovw04C7he0q65+K85vu8Bp+ayLwE754vAvwMOy+V/C/w2Ijb24iMx24JH9TQrkKRPAk9HRJfP2JXUBpwVEfNr3OYFwNyIaO2fKK2MXPM3K9Z3gef6eZuLnPhta7nmb2ZWQq75m5mVkJO/mVkJOfmbmZWQk7+ZWQk5+ZuZldD/Ae4lYBqFnlZDAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "D\tR\tVobs\te_Vobs\tVgas\tVdisk\tVbul\tSBdisk\tSBbul\n",
      " 4.35\t 0.42\t14.20\t 1.91\t 4.87\t 4.78\t 0.00\t11.05\t 0.00\n",
      " 4.35\t 1.26\t28.60\t 1.82\t13.14\t10.76\t 0.00\t 5.81\t 0.00\n",
      " 4.35\t 2.11\t41.00\t 1.74\t19.65\t13.60\t 0.00\t 2.75\t 0.00\n",
      " 4.35\t 2.96\t49.00\t 1.91\t22.42\t13.29\t 0.00\t 1.00\t 0.00\n",
      " 4.35\t 3.79\t54.80\t 2.05\t22.82\t12.56\t 0.00\t 0.72\t 0.00\n",
      " 4.35\t 4.65\t56.40\t 3.12\t21.37\t12.33\t 0.00\t 0.43\t 0.00\n",
      " 4.35\t 5.48\t57.80\t 2.83\t18.73\t12.04\t 0.00\t 0.18\t 0.00\n",
      " 4.35\t 6.33\t56.50\t 0.65\t16.75\t10.62\t 0.00\t 0.00\t 0.00\n"
     ]
    }
   ],
   "source": [
    "plot_single_galaxy_mm(massmodels,gal_id)\n",
    "\n",
    "# Print the raw data for this galaxy as well\n",
    "labels = mm_dict[gal_id].dtype.names[1:]\n",
    "print('\\t'.join(labels))\n",
    "for mm_line in mm_dict[gal_id]:\n",
    "    print_arr = ['{:>5.2f}'.format(arr_item) for arr_item in [item for item in mm_line][1:]]\n",
    "    print('\\t'.join(print_arr))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "7d5df59d-2001-4308-8bff-10fbf283b2d4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAhDUlEQVR4nO3de5wcVZn/8c83V7KGHRQiNxMG5KJRENhZJYI6Gi+Aiygs3lbCzR+gqLAGVMBdUFbRnysKXlAQidG4iAsiCyhiZFyRER1CCIRwT0iEAAFl5JqQzLN/1Omk0/TM1Eympme6vu/Xq1/pPlV16qlKz9OnT586pYjAzMzKY0yjAzAzs+HlxG9mVjJO/GZmJePEb2ZWMk78ZmYl48RvZlYyTvxmZiXjxF8ykjokfbjRcRRF0lOSdmpwDP8h6TFJD0ualmIam5YN6vxLapfUk+raf+ijHpkk7ZqOeV0zv2+HmxN/E5K0TNKz6Q/mEUlzJE0eonrfOhQxDmLfuRJmREyOiPsHWHd1Un1S0l2SjhpknNOA2cD0iNgmIpanmNYNpr4aD6W6fpn2daakH9WJISTtXPX6bZKuT8f2uKSFkj4tabOqdXaV9NP0gdUtaZGkT1Y+sNI6k9M5+kUfx7+LpOfqxZWWf786PkkTJV0k6YEU30JJB1TWj4i7I2Iy8LsBnivrgxN/8zoo/cHsDbQBn21wPIWSNG4Tq3gona+/Bz4NXChp+iD2Mw14PCIe3cR4hoSkw4D/Bn4M7BARWwLvA14GTE3rvBy4CVgB7B4RLcBhZO+bzauqOxRYDbxN0ja97PJbwJ96iWU/4OU1xePSft8EtJC9Ty+V1DqgA7UBceJvchHxIPAL4NVVxTtI+n1qYf1K0laVBZLeJWmxpCdSK/uVqfyHZEntf1Kr71N9rZ+WLZN0Smo9Pp1adltL+kXa968lvbhq/X0k3ZjqulVSeyr/AvAG4Jtp399M5SHpBEn3APdUlVVak5MkfTW1Jrsl3SBpUj/nKyLiCuCvwHRJR6Zz9TVJjwNnSmqRNFfSqlT3ZyWNSd+GrgO2S3HOkdSaYqr7gSHpaElLJP1V0rWSdujv/zQvSQLOAT4fERdGxF/SMd4VER+PiHvSqp8DboyIT0bEyqp1PhgRT1RVeQTwHWAR8KE6+3s/8AQwv86yccA3gI9Xl0fE0xFxZkQsi4ieiLgKWAr8w6Ycu/UjIvxosgewDHhrej4VWAyclV53APcBuwKT0usvpWW7Ak8DbwPGA58C7gUm1NY7gPX/AGwNbA88CiwA9gI2A34DnJHW3R54HDiQrEHytvR6SlXcH645ziBLtC8BJlWV7Zyefytttz0wFng9MLHO+WoH/pyejwHeAzwP7AYcCawlS1jj0jmbC/ycrDXcCtwNHFNbV3rdmmIaV3scwMHpfL0y1f1ZsgRc7/90o3pT2ZnAj+qsG8DOwCvS89Z+3i8PA0f1s84OQA8wnawra1HN8r9P5+Fl9eICTgHOrf0/qrOfrYHngFfUlL/g/9+PwT/c4m9eV0h6ArgB+C3wxaplF0fWd/oscCmwZyp/H3B1RFwXEc8D/0mW6F7fyz7yrP+NiHgksm8evwNuiohbIuI54GdkHwKQtSCviYhrImv5XQd0kX0Q9OXsiPhLOpb1JI0BjgZOjIgHI2JdRNwYEat7qWe7dL4eA84ADo+Iu9KyhyLiGxGxFlgDvB84NSKejIhlwFeBw/uJs57jU/xLUt1fBPYcwlZ/5Zvcw5UCSZekb1TPSKrEvCWwsp+6DidL9ncAlwCvkrRX1fKzgIsi4s+1G0qaChwH/HtfO5A0HpgH/CAi7uwnHtsEm9ovaiPXuyPi170se7jq+TNA5Yff7YAHKgsiokfSCrIWcz151n+k6vmzdV5X9r0DcJikg6qWjweu72XfFSt6Kd+K7FvFff1sX/FQRLwsxz62SnE9UFX2AL2fo77sAJwr6atVZUp1PVB/k42sTbFs2DhLnpB9Y3k8Pd+WrPuEiHh/Wu8Gsm9BpPW27Wdfs4ALUx0PSvotWdfPLZL2BN7Khg/xWl8n627q7q3y9EH9Q7IP1o/1E4ttIrf4rdpDZMkIWN9HPBV4MBXVzuHd3/oDsQL4YURsUfV4UUR8qZd900/5Y2RdBrU/Jg5G9T4eI0uq1a3yaQz+mI+rOeZJEXFjzu2Xk3UlVduR7APhQeCu9O8h/dTza7IfbuuS9HpgF+BUZUNUHwZeB3ww9d23pziWp2UnA4dKWpCqmAl8pWpbgE5JH0z1C7iIrJvn0PTt0QrkxG/VLgXeKWlmajnOJhvFUUlEjwA7DWD9gfgRcJCkd0gaK2kzZcMsK63w2n33KSJ6gO8D50jaLtU5Q9LEQcRWXe86suP+gqTNU7fMJ1P8A/UdsmT6KoD0o/FhA9j+l8ArJB0uabykl5B1F10WEWvTOZgNnCHp/0l6sTK7kCXZijOA10v6itJoHUk7S/qRpC3IWvbXkfXv75kerybr1jsAuIDsA7ay7DvA1cA7Uv27Aq+pWg5wEFlXH8D5ZL9zHFTbZWfFcOK39VKf9ofIRl88RvbHeVBErEmrnA18NvURn5xj/YHsewXZj52nAavIWsOnsOE9ei7wz2n0y3k5qz0ZuI1seOFfgC8zNO/5j5P9qH0/2W8oPyb7kBmQiPhZiukSSX8DbidLpHm3fzStfxzZD+e3k42q+UjVOj8B3kv2/7SC7P/pUrJk/dO0zn3ADLJW+2JJ3cBlZL+xPJ+2/0ZEPFz1WErWNXNERDxTvQx4CnguIlZV4qxZDvBYRDybPjiPI/tAeDiNhnpK0r/kP5M2UIrwHbjMRjpJbwSuJftG9b6IuLbBIQ2L9O3kT8AE4KMRMaexETUHJ34zs5JxV4+ZWck48ZuZlcyoGMe/1VZbRWtra6PDMDMbVW6++ebHImJKbfmoSPytra10dXU1Ogwzs1FFUt0LAd3VY2ZWMk78ZmYl48RvZlYyTvxmZiXjxG9mVjJO/GZmJePEb2ZWMk78ZmYl48RvZjZCtc9pp31O+5DX68RvZlYyTvxmZiXjxG9mVjJO/GZmJePEb2ZWMk78ZmYl48RvZlYyTvxmZiXjxG9mVjJO/GZmJVNY4pc0VdL1ku6QtFjSiTXLZ0sKSVsVFYOZmb1QkTdbXwvMjogFkjYHbpZ0XUTcIWkq8HZgeYH7NzOzOgpr8UfEyohYkJ4/CSwBtk+LvwZ8Coii9m9mZvUNSx+/pFZgL+AmSQcDD0bErf1sc6ykLkldq1atGo4wzcxKofDEL2kycBlwEln3z2nAv/e3XURcEBFtEdE2ZcqUYoM0MyuRQhO/pPFkSX9eRFwOvBzYEbhV0jLgZcACSdsUGYeZmW1Q2I+7kgRcBCyJiHMAIuI24KVV6ywD2iLisaLiMDOzjRXZ4t8XOBx4i6SF6XFggfszM7McCmvxR8QNgPpZp7Wo/ZuZWX2+ctfMrGSc+M3MSsaJ38ysZJz4zcxKxonfzGyE6l7dzfLu5XSu6BzSep34zcxGoM4VnSx6ZBFLn1jKzLkzhzT5O/GbmY1AHcs66IkeANasW0PHso4hq9uJ38xsBGpvbWeMshQ9YewE2lvbh6zuIufjNzOzQZoxdQZ7bL0H3c91M++QecyYOmPI6nbiNzMboVomttAysWVIkz64q8fMrHSc+M3MSsaJ38ysZJz4zcxKxonfzJpG+5x22ue0NzqMEc+J38ysZJz4zcxKxonfzKxknPjNzErGid/MrGSc+M3MSsaJ38ysZJz4zcxKxonfzKxknPjNzErGid/MrGSc+M3MSsaJ38ysZJz4zcxKxonfzJpG9+pulncvp3NFZ6NDGdGc+M2sKXSu6GTRI4tY+sRSZs6d6eTfByd+M2sKHcs66IkeANasW0PHso7GBjSCjetroaTNgH8C3gBsBzwL3A5cHRGLiw/PzCyf9tZ2xmgMPdHDhLETaG9tb3RIm6zjyI5C6u018Uv6HFnS7wBuAh4FNgN2Bb6UPhRmR8SiQiIzMxuAGVNnsMfWe9D9XDfzDpnHjKkzGh3SiNVXi/+PEXFGL8vOkfRSYFoBMZmZDUrLxBZaJrY46fej18QfEVf3tWFEPEr2LcDMzEaRXn/clTRW0nGSzpK0b82yz/ZXsaSpkq6XdIekxZJOTOVnSVokaaGkX0nabtMPw8yK1D6nnfY57Y0Ow4ZIX6N6vgu8CXgcOE/SOVXLDslR91qy3wCmA/sAJ0iaDnwlIvaIiD2Bq4B/H1TkZmY2KH0l/tdGxAcj4uvA64DJki6XNBFQfxVHxMqIWJCePwksAbaPiL9VrfYiIAYdvZmZDVhfiX9C5UlErI2IY4GFwG+AyQPZiaRWYC+y0UFI+oKkFcC/0EuLX9Kxkrokda1atWoguzMzsz70lfi7JO1fXRARnwcuBlrz7kDSZOAy4KRKaz8iTo+IqcA84GP1touICyKiLSLapkyZknd3ZmbWj14Tf0R8KCJ+Waf8exExPk/lksaTJf15EXF5nVXmAYfmDdbMzDZdv1M2SBo7mIolCbgIWBIR51SV71K12sHAnYOp38zMBqe/KRs2B/6L7AregdoXOBy4TdLCVHYacIyk3YAe4AHg+EHUbWZmg9TXlA3bAlcAXxhMxRFxA/VH/1wzmPrMzGxo9NXi/x1wSkRcOVzBmJlZ8frq4/8rsP1wBWJmI5dvcNJc+kr87cABkk4YpljMbATyDU6aT1+TtD0t6V1kUzeYWUnVu8HJSJ39sqj565tNn8M5I2JdRHx4uIIxs5GncoMToGlucFJ2ecbx75huulJ5PSlNwWBmJVC5wcmOW+zI/FnzR2xr3/LLc8/dn5KNua9Yl8rMrCRaJrYwrWWak36TyJP4x0XEmsqL9HxCH+ub2Sby/PdWpDyJf1X6kRcASQcDjxUXkpmZFanPKRuS44F5kr5FNnf+n4FZhUZlZmaF6TfxR8R9wD5pemUi4qnCozIzs8LkGdWztaSLgJ9GxFOSpks6ZhhiMzOzAuTp458DXAtUbop+N3BSQfGYmVnB8iT+rSLiUtKQzohYSzak08zMRqE8if9pSVuSboouaR+gu9CozMysMHkS/yeBK4GXS/o9MBf4eKFRmZWcZ8O0Iiki+l9JGgfsRnZjlbsi4vmiA6vW1tYWXV1dw7lLs4bpXNHJfhfvR0/0MGncJE+TYIMm6eaIaKstzzOq5zBgUkQsBt4N/ETS3kMfoplB/dkwzYZSnq6ef4uIJyXtB8wku4H6+cWGZVZeng3TipYn8VdG8LwTuDAirsZz9ZgVxrNhWtHyTNnwoKTvAm8DvixpIvk+MMxskFomttAyscVJ3wqRJ4G/l+wCrndExBPAS4BTigzKzMyKk6fFvwdwXUQ8mV4/jcfxm5mNWnla/OcD1ROzPYV/3DUzG7XyJH5F1WD/iOgh3zcFMzMbgfIk/vslfULS+PQ4Ebi/6MDMzKwYeRL/8cDrgQfJbsLyOuDYIoMyM7Pi9NllI2ks8LWIeP8wxWNmQMeRHY0OwZpYny3+iFgH7CDJF2yZmTWJPD/S3g/8XtKVZEM5AYiIcwqLyszMCpMn8d+XHmOAzYsNx8zMipbnZuufG45AzMxsePSb+CVNAT4FvArYrFIeEW8pMC4zMytInuGc84A7gR2BzwHLgD8VGJOZmRUoT+LfMiIuAp6PiN9GxNGAW/tmZqNUnsRfuc3iSknvlLQX2QydfZI0VdL1ku6QtDhd8Yukr0i6U9IiST+TtMXgwzczs4HKk/j/Q1ILMBs4Gfge8K85tlsLzI6I6cA+wAmSpgPXAa+OiD2Au4FTBxW5mZkNSq8/7kr6APCriLgqFXUDb85bcUSsBFam509KWgJsHxG/qlrtD8A/DzhqMzMbtL5G9UwDfippPDAf+AXwx+qZOvOS1ArsBdxUs+ho4CcDrc/MzAav166eiPhyGrJ5IHArWZJeIOnHkmZJ2jrPDiRNBi4DToqIv1WVn07WHTSvl+2OldQlqWvVqlX5j8jMzPqkgTbgUz/9AcDbI+Id/aw7HrgKuLZ6igdJRwLHATMj4pn+9tnW1hZdXV0DitPMrOwk3RwRbbXleS7g2rtO8RXAuf1sJ+AiYElN0t+f7IKwN+VJ+maN1D6nHfBsmdZc8szV821gb2ARIODVwGKgRdJHan6srbYvcDhwm6SFqew04DxgInBd9tnAHyLi+EEfgZmZDUiexP8QcExELIb1XT2fJ2u1Xw7UTfwRcQPZB0WtawYXqpmZDYU84/h3rSR9gIi4A3hFRPj2i2Zmo1CexL9Y0vmS3pQe3wbukDSRDVf1mjWl7tXdLO9eTueKzkaHYjZk8iT+I4F7gZPS4/5U9jwDuKDLbLTpXNHJokcWsfSJpcycO9PJ35pGnvn4nwW+mh61nhryiMxGiI5lHfREDwBr1q2hY1kHM6bOaHBUZpsuz3DOXYCzgelsPB//TgXGZdZw7a3tjNEYeqKHCWMn0N7a3uiQzIZEnq6ei4Hzya6yfTMwF/hRkUGZjQQzps5gj633YMctdmT+rPlu7VvTyJP4J0XEfLKrfB+IiDOBdxYbltnI0DKxhWkt05z0rankGce/WtIY4B5JHwMeBCYXG5aZmRUlT4v/RODvgE8A/0B2Ne4RRQZlZmbFyTOqp3J/3aeAo4oNx2xk8Rw91ozyjOppA04HdqheP91By8zMRpk8ffzzgFOA24CeYsMxM7Oi5Un8qyLiysIjMTOzYZEn8Z8h6Xtkt19cXSmMiMsLi8rMzAqTJ/EfBbwCGM+Grp4gm5LZzMxGmTyJ/x8jYrfCIzEzs2GRZxz/jenmK2Zm1gTytPj3ARZKWkrWxy8gPJzTzGx0ypP49y88CrNB8s3QzQau18QvaXJEPBURD/S3TjGhmZlZEfrq4/+5pK9KeqOkF1UKJe0k6RhJ1+JvA9ZgvjWi2cD1mvgjYibZ2P3jyO672y3pcbK5+LcBjoiI/x6eMM1eyLdGNBucPvv4I+Ia4JphisVsQHxrRLPByTOc02xEqtwaEfCtEc0GIM+oHrMRqXJrxO7nupl3yDy39s1ycuK3Ua1lYgstE1uc9M0GIM98/NPqlUfE8qEPx2xgPH7fbODytPivJpuUTcBmwI7AXcCrCozLzMwKkufWi7tXv5a0N/DRwiIyM7NCDXhUT0QsAF5XQCxmZjYM8vTxf7Lq5Rhgb+ChwiIyM7NC5enj37zq+VqyPv/LignHzMyKlqeP/3PDEYiZmQ2PPF09bcDpwA7V63s+fjOz0SlPV8884BTgNjbcc9fMzEapPIl/VURcWXgkZmY2LPIM5zxD0vckfUDSIZVHfxtJmirpekl3SFos6cRUflh63ZO6kawk2ue0r79jlpk1Tp4W/1HAK4DxbOjqCeDyfrZbC8yOiAWSNgdulnQdcDtwCPDdwYVso1X36m66n+umc0Wn59Yxa6A8if8fI2K3gVYcESuBlen5k5KWANtHxHUAkgZapY1ilZum9EQPM+fOZP6s+U7+Zg2Sp6vnRknTN2UnklqBvYCbBrDNsZK6JHWtWrVqU3ZvI0C9m6aYWWPkafHvAyyUtBRYTTZZW+QdzilpMtkFXydFxN/yBhYRFwAXALS1tUXe7Wxkqtw0pSd6fNMUswbLk/gHfUN1SePJkv68iOjvNwFrYr5pitnIkefK3QcGU7GyTvyLgCURcc5g6rDmcstxtzQ6BDOj2Dtw7QscDtwmaWEqOw2YCHwDmAJcLWlhRLyjwDjMzKxKYYk/Im4g+z2gnp8VtV8zM+vbgOfjt3LyxVdmzcOJ33LpXt3N8u7ldK7obHQoZraJnPitX5WLr5Y+sZSZc2c6+ZuNck781i9ffGXWXJz4rV+Vi68AX3xl1gSKHM5pTWLG1BnccNQNdCzroL213RdfmY1yTvyWy4ypM5zwzZqEu3pKZK/v7sVO5+7kH2fNSs6JvyQ8MsfMKpz4S8Ijc8yswom/JDwyx8wq/ONuSXhkjplVOPGXiEfmmBm4q8fMrHSc+Ee5zhWdnP27sz1Kx8xyc1fPKNa5opP9Lt6Pnuhh0rhJzJ813105ZtYvt/hHMQ/RNLPBcOIfxTxE08wGw109o5iHaJrZYDjxj3IeomlmA+WuHjOzknHiNzMrGSf+EcpTKJtZUZz4RyBPoWxmRXLib4D+WvMen29mRXLiH2Z5WvMen29mRXLiL0hvrfo8rfnK+PwvvuWLnobBzIacx/EXoNKq74keZs6duVHybm9tZ9K4SaxZt6bP1rzH55tZUZz4C1CvVV9J4jOmzmD+rPm+2tbMGqYUib9zRSeHXnooa9at4bDphzHrNbMKTbj9terdmjezRlJENDqGfrW1tUVXV9egtu1c0cm+39+XYMNxThw7keuPuH6j5Nu5onNIW+FDXZ+Z2UBJujki2mrLm77F37GsY6OkDy/sfqmd1/7r+3+dx595fJOStlv1ZjZSNX3ib29tZ8LYCaxZt2Z9WRBs+Xdbrn9d3Se/eu1qjr/qeILwzU3MrCmVYjjn0Xsezbt3ezdvnPbG9WUn/fKk9UMtK33yYzWWMWPGrP+GkPfiKU+vYGajSVMn/koXzndu/g7X3nct06dMX7+sOqlXRtqc9eaz+NaB3+rz4qnaJO/pFcxstCmsq0fSVGAusDUQwAURca6klwA/AVqBZcB7I+KvRcRQO6wS2Gi0zcULL+bCBRcy75B5G/XJ7/7S3df/MPvRaz5K93PdzDtkHsALxuf3NXTTzGwkKrLFvxaYHRHTgX2AEyRNBz4DzI+IXYD56XUhqrtwJoydwKzXzFrfst9y0pbc85d76rbUZ0ydwalvOBVgo9b83FvnviDJe3oFMxttCmvxR8RKYGV6/qSkJcD2wMFAe1rtB0AH8OkiYujrYqnTf3P6+uer162u21Kv941hjMbQEz3rk7xvf2hmo82wjOqR1ArsBdwEbJ0+FAAeJusKKky9YZW1QzzHamzdlnrthVizXjOLWa+Z9YIk76GbZjaaFJ74JU0GLgNOioi/SVq/LCJCUt0ryCQdCxwLMG3atCGNqdI90xM9jBszjm8e+M26ibu3bwxO8mY2mhV65a6k8cBVwLURcU4quwtoj4iVkrYFOiJit77q2ZQrd3vjK2vNrNkN+5W7ypr2FwFLKkk/uRI4AvhS+vfnRcXQV3J394yZlVWRXT37AocDt0lamMpOI0v4l0o6BngAeG8RO6+dhsFX4JqZZYoc1XMDoF4WzyxqvxUeX29mVl/TXrlbO4bf4+vNzDJNO0mbb3hiZlZf0yZ+8A+4Zmb1NG1Xj5mZ1efEb2ZWMk78ZmYl48RvZlYyTvxmZiXjxG9mVjJO/GZmJVPo7JxDRdIqsnl98toKeKygcEYbn4sNfC428LnYWLOejx0iYkpt4ahI/AMlqaveVKRl5HOxgc/FBj4XGyvb+XBXj5lZyTjxm5mVTLMm/gsaHcAI4nOxgc/FBj4XGyvV+WjKPn4zM+tds7b4zcysF078ZmYl01SJX9L+ku6SdK+kzzQ6nkaQtEzSbZIWSupKZS+RdJ2ke9K/L250nEWQ9H1Jj0q6vaqs7rErc156ryyStHfjIh96vZyLMyU9mN4bCyUdWLXs1HQu7pL0jsZEXQxJUyVdL+kOSYslnZjKS/negCZK/JLGAt8CDgCmAx+QNL2xUTXMmyNiz6pxyZ8B5kfELsD89LoZzQH2rynr7dgPAHZJj2OB84cpxuEyhxeeC4CvpffGnhFxDUD6O3k/8Kq0zbfT31OzWAvMjojpwD7ACemYy/reaJ7ED7wWuDci7o+INcAlwMENjmmkOBj4QXr+A+DdjQulOBHxv8Bfaop7O/aDgbmR+QOwhaRthyXQYdDLuejNwcAlEbE6IpYC95L9PTWFiFgZEQvS8yeBJcD2lPS9Ac2V+LcHVlS9/nMqK5sAfiXpZknHprKtI2Jlev4wsHVjQmuI3o69rO+Xj6Xui+9XdfmV5lxIagX2Am6ixO+NZkr8ltkvIvYm+7p6gqQ3Vi+MbPxuKcfwlvnYk/OBlwN7AiuBrzY0mmEmaTJwGXBSRPytelnZ3hvNlPgfBKZWvX5ZKiuViHgw/fso8DOyr+yPVL6qpn8fbVyEw663Yy/d+yUiHomIdRHRA1zIhu6cpj8XksaTJf15EXF5Ki7te6OZEv+fgF0k7ShpAtmPVVc2OKZhJelFkjavPAfeDtxOdh6OSKsdAfy8MRE2RG/HfiUwK43g2Aforvra35Rq+qnfQ/begOxcvF/SREk7kv2o+cfhjq8okgRcBCyJiHOqFpX3vRERTfMADgTuBu4DTm90PA04/p2AW9NjceUcAFuSjVq4B/g18JJGx1rQ8f8XWRfG82T9ssf0duyAyEaB3QfcBrQ1Ov5hOBc/TMe6iCy5bVu1/unpXNwFHNDo+If4XOxH1o2zCFiYHgeW9b0REZ6ywcysbJqpq8fMzHJw4jczKxknfjOzknHiNzMrGSd+M7OSceK3piFpXZp18nZJ/yNpiwFu3yGpLT2/ZqDb91LnXpIuSs+PT7ND3i3pzFS2X5peY7Gkn6ex9C+V9Os0y2qXpJ3Tuv8p6S2bGpOZE781k2cjm3Xy1WQTlJ0w2Ioi4sCIeGIIYjoNOC89v5dsnpjdgSMkvQx4jmzc/KuAZ4DDgHHAyRGxO9kVtpVZI79B886sasPIid+aVSdpYi1Jr5XUKekWSTdK2i2VT5J0iaQlkn4GTKpsrOy+BltJaq2Z0/7kqtb6J9Ic74skXVIbQLqKeo+IuBUgIn4d2cyxIkvuayKiK7LpNQAmAs9FxEMRsbC6LG3/ALClpG2G7CxZKY1rdABmQy3NJT+T7DJ9gDuBN0TEWklvBb4IHAp8BHgmIl4paQ9gwQB39Rlgx4hY3Uu3UBsbpkWodgHZNMjr50ySdAywDVXTaUjaEzgJqO7eWQDsSzbvjNmgOPFbM5kkaSFZS38JcF0qbwF+IGkXskv3x6fyN5K6YSJikaRFA9zfImCepCuAK+os3xZYVV0g6V2p/MiqsinAGcDeEfF81erfB46MiGVVZY8C2w0wTrONuKvHmsmzEbEnsANZd0qlj/8s4PrU938QsNkA6lzLxn8n1du+k2xOl72BP0mqbUg9W2dfewC/imyGzIrdgNsi4rGadXeO7IYq1TZL9ZoNmhO/NZ2IeAb4BDA7JeMWNkyre2TVqv8LfBBA0qvJknKtR4CXStpS0kTgn9L6Y4CpEXE98Om0j8k12y4Bdq4pu4IXzhp7N/ClOvs+qk7ZrtTvPjLLzYnfmlJE3ELWFfMB4P8DZ0u6hY27N88HJktaAnweuLlOPc+nZX8k6zq6My0aC/xI0m3ALcB5taOAIuJOoKUyVXayH/C6mt1MIxvNU2t29Ys0p/zOQFf9ozbLx7NzmhVI0r8CT0bE94agrveQ/Q7wb5semZWZW/xmxTofWD1EdY2jZLdLtGK4xW9mVjJu8ZuZlYwTv5lZyTjxm5mVjBO/mVnJOPGbmZXM/wFp4/WSNvMP6QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Photometric Data\n",
      "Radius    mu        Kill      Error   \n",
      "    2.40    19.906     False     0.025\n",
      "    2.88    20.077      True     0.029\n",
      "    3.46    20.173      True     0.031\n",
      "    4.15    20.175      True     0.031\n",
      "    4.98    20.160      True     0.031\n",
      "    5.97    20.109      True     0.036\n",
      "    7.17    20.015      True     0.040\n",
      "    8.60    20.040      True     0.036\n",
      "   10.32    20.146      True     0.031\n",
      "   12.38    20.117      True     0.037\n",
      "   14.86    20.047      True     0.040\n",
      "   17.83    20.057      True     0.049\n",
      "   21.40    20.077      True     0.049\n",
      "   25.68    20.246      True     0.039\n",
      "   30.81    20.293      True     0.037\n",
      "   36.98    20.325      True     0.038\n",
      "   44.37    20.497      True     0.041\n",
      "   53.25    20.671      True     0.048\n",
      "   63.90    20.809      True     0.054\n",
      "   76.68    20.994      True     0.065\n",
      "   92.01    21.380      True     0.094\n",
      "  110.41    21.810      True     0.132\n",
      "  132.49    22.486      True     0.222\n",
      "  158.99    22.791      True     0.283\n",
      "  190.79    23.171      True     0.381\n",
      "  228.95    23.740      True     0.578\n"
     ]
    }
   ],
   "source": [
    "photo_data = import_photometric_profile(datapath='./sparc/sfb_LTG/',gal_ID=gal_id)\n",
    "\n",
    "plot_single_galaxy_photo_profile(photo_data, gal_id)\n",
    "\n",
    "# Print the raw photometric data for this galaxy\n",
    "photo_labels = photo_data.dtype.names\n",
    "print('Photometric Data')\n",
    "print('  '.join(['{:<8}'.format(item) for item in photo_labels]))\n",
    "for photo_line in photo_data:\n",
    "    photo_line_zip = zip([item for item in photo_line],photo_labels)\n",
    "    photo_print_arr = ['{:>8}'.format(str(item)) if name == 'Kill' else '{:>8.2f}'.format(item) if name == 'Radius' else '{:>8.3f}'.format(item) for item,name in photo_line_zip]\n",
    "    print('  '.join(photo_print_arr))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "ddc79bc7-bd3c-462f-ba6a-8702337a6be7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEWCAYAAABhffzLAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAh1ElEQVR4nO3deZxcVZn/8c+XkBBIWEMGwhKCLCKigrZAEDFj0B+bwG9QNoHojzHoCKjAsMzoT3AZxB0UwSCrIBEE2UEQaRBtkASQxaAiSwgQQtgkbCHhmT/OKah0uqurO3W7um59369Xv7r61q17n6rufu6p55w6RxGBmZm1j+WaHYCZmQ0uJ34zszbjxG9m1mac+M3M2owTv5lZm3HiNzNrM078LUxSSNq4Sef+pKTr69jvHEnfGIyYWllfr6ekD0r6awHnnSTpDUkLJO3U6OMPVZI2zc95saR/b3Y8g82Jv4kkPSLplfwH+JykqyWt38R4/pr/Ic6RtFDSi/nrPkknSlq1sm9EXBARH23w+TslvZrP+U9JMyUdK2mFRp5nKOr+ena/qEfE7yPi7QWd/omIGB0R1+VzHy/p/O47dY9J0kck3ZR/X89IulvSMZJGVu2zqaSLJc2X9IKkeyQdIWlY1T6j8//Atb0FKGmT/LexVFz5/rOq45O0gqQzJT2a47tb0s6V/SPibxExGvh9P1+rUnDib76P5T/AccBTwI+aEYSkjYBhEfG3vOnbEbEyMBb4NLAt8AdJowoO5dB83nHAkcC+wDWSVPB5rR8kfQL4FfALYIOIGAPsA6wHrJ/32Qi4HXgMeFdErAp8AugAVq463F7Aa8BHJK3dyylPBe7oJZbtgY26bV4+n/dDwKrAl4GLJE3o1xMtKSf+ISIiXiX9I21e2ZZbwP9e9fOnJN3a0+MljZF0ZW4p3yHpG9X7StpM0g2Sns0t+727HWJX4Jqe4oqIO4DdgTGki8ASsSj5gaR5+fz3StqihxhXzi3EU/pK5BHxUkR05vNOzPEhabn8LuAfuZV5kaQ1qs6xvaQ/Snpe0mOSPpW3ryrpPElP51bglyUtV/Vc/pCfw/OSHpK0Xd7+WH5eU6rOcY6k0/Pr+aKkmyVtUHX/dvl38EL+vl3VfZ/Kx39R0sOSPtnD63lL3v3PuSW8j1JJZk7Vcd6R/z6el3S/pN27xXeq0jvIFyXdnpNwQ+Tf3feBr0XEGRHxLEBE/DUiDouIv+ddTwD+GBFHRMSTVfvsHxHPVx1yCnA6cA9wQA/n2xd4Hrixh/uWJzWWDqvenv9+jo+IRyLijYi4CngYeN+yPPeycOIfIiStRGox3TbAQ5wKvASsTfpHqk5Uo4AbSK2zfyG1on8iafOqx+8CXN3bwSPixXyMD/Zw90eBHYBNSa2rvYFnqneQNIb0j/uHiDg86pwrJCJmAzOqznsYsCepJbcO8BzpuZOT77WkRDAW2BK4Oz/uRzm2t+XHHkS+iGXbkBLPGNLrNB14P7AxKRn9WNLoqv0/CXwdWDOf44Icwxqk1/GUfKzvA1fnC/OovH3n/K5mu6r4qp/zDvnme3IJ5pfV90saDlwJXE/6fR4GXCCpuhS0Lynxrg48CHyz+3mWwdtJLftL+thvR1Jjplf5dzaJ9PpdQPq9VN+/CvA14IheDvEl4JaIuKeP86xF+vu8v4+Y24ITf/NdJul54AXgI8B3+nsApXrpXsBXI+LliPgLcG7VLrsBj0TE2RGxKCLuIv3TfiI/fiVSkuvs41RPAGv0sP110lv3zQBFxKxKCy9bB7gZuDgivtzf59ftvJ8F/jsi5kTEa8DxwMdzy29/4LcRcWFEvB4Rz0TE3fn12Rc4LiJejIhHgO8BB1ad4+H8+iwGfkkqV3wtIl6LiOuBhaSLQMXVEXFLjuG/gYlK/TO7An+PiJ/n1/pC4AHgY/lxbwBbSFoxIp6MiIEkom2B0cC3ImJhRPwOuArYr2qfX0fEnyJiESmhbjmA8/Rmzfx9bmWDpOn53cfLkiqv6xjgyaUevaQDgXvy3+x04J2Stqq6/+vAmRExp/sD8+t9CPD/a50gXygvAM6NiAf6iKctOPE3354RsRowEjgUuFm91zl7M5a3apoV1bc3ALbJ/5jP5wvNJ0nvDgAmk96Sv9bHedYFnu2+MSeeH5Na3vMkTcsttYpdgRVJb+cHovq8GwC/rnoes4DFwFqkZP2PHh6/JjAceLRq26P5uBVPVd1+BSAium+rbvG/+fpGxIIc3zr5q/o8b54rIl4ivav7LPBkLsVs1vNTrmkd4LGIeKPG85lbdfvlbrH3ZRHp9XpTTp6QLvKVd3PjKvdHxL757/hOoNJx+0z1Pr04iPxuKSIeJzUQpuRzbkl61/CDXh77Q9LF+YXeDp7LeT8nXbgP7SOWtuHEP0RExOKIuJSUxLbPm18CVqrarbcLwtOkf9b1qrZVjw56DLg5Ilar+hodEZ/L9+9CD/X9arnMsSO9jIKIiFMi4n2kPopNgf+suvsM4DpSJ22/Oodzq+59Ved9jFQqqX4uI3PSeIylO/kA5pMS1gZV28YDj/cnlm7efH3za7MG6Z3JE93Os8S5IuI3EfERUkJ8gPTa9NcTwPqVPoru52iA2cCEbts2JP2NPQ78NX//tz6O81vSO9Ee5b6PTYDjJM2VNJdUcts/v4OblOOYne87CthL0p35EJOB71Q9FqBL0v75+ALOJDUK9oqI1/uIt2048Q8RSvYg1WRn5c13A/8maSWlYWoH9/TYXJ64FDg+77sZS9ZKrwI2lXSgpOH56/2S3pHv35le6vtKw+LeB1xGqqef3cM+75e0TW4VvgS8SippVDuUlDCulLRizRcjHXMlSR8CLgf+xFsXptOBb1Y6UyWNza8bpJbjjpL2lrR8rqtvmV+fi/LjVs6PPQLocWhgnXZR6kgeQSpH3BYRj+U4N5W0f45hH9LF8CpJa0naI1/8XgMWsPTrVPEUqT+iJ7eTWvFH59/lJFIpafoyPJ9q1wGbVf29rAH8D3BJLl+9QRpx9VVJn5G0ev773YSUZCu+Cmwn6TuVd7GSNpZ0vqTVSC37G0ivz5b5awvSu8OdgWmkC3nlvtNJf6f/Jx9/U+A9VfeTX4df59unAe8gjZx7pSGvTEk48TfflZIWAP8kdcBNqar7/oD0FvUpUs3+ghrHOZTUeTmX9Nb2QlJyqXTMfpRU534i73MSsILS6JsFuRO12tGSXiS9XT8PmAlsl8sV3a1Cark+Ryo5PEO3vorcmTsVmANcrqqx3t38OJ/3KdJb+UuAnarKGicDVwDX5/1uI7USKx3Bu5CS0rOkC+d78uMOI12UHgJuJXXgntVLDPX4BSmxPUt6R3JAjuEZUp/KkaTX4Whgt4iYT/p/O4L0O3iW1Mn8uaWOnBwPnJtLWkuMwIqIhaQEtzPp3cxPgIMaVb+OiHn52IcA84D7SKNqPle1zy9JnfgHkN5pzSddXKcBF+d9/kEakTUBuF/SC6Tf5wzSO7C9gR9FxNyqr4dJf79Tcn/Vm/eRLpSvRsTTlTi73Q8wPyJeyRf3Q0gXhLlKo6MWKI+ianeqc3CFtRhJJwFrR8SUPvY7GlgzIo4enMhan6RzgDkD7KgeUiTtAPyG1EjYJyJ+0+SQBkV+d3IHMAL4j4g4p7kRDa7lmx2ANUYu74wA7iWN0DkYqOej6I+QhgZaG4qIW0illbYS6bMGqzU7jmZx4i+PlUnlnXVIZZLvkerjNUXERQXHZWZDjEs9ZmZtxp27ZmZtpiVKPWuuuWZMmDCh2WGYmbWUmTNnzo+Isd23t0TinzBhAjNmzGh2GGZmLUVS90+RAy71mJm1HSd+M7M248RvZtZmnPjNzNqME7+ZWZtx4jczazNtl/i7uuDEE9N3M7N21BLj+BulqwsmT4aFC2HECLjxRpg4sdlRmZkNrrZq8Xd2pqS/eHH63tnZ7IjMzAZfWyX+SZNSS3/YsPR90qRmR2RmNvjaqtQzcWIq73R2pqTvMo+ZtaO2avFDSvbHHZduu5PXzNpRoS1+SV8irQIVpJWhPg2MIy0KPYa0juuBeQ3RQdNXJ29Xl98VmFl5Fdbil7QucDjQERFbAMNIi32fBPwgIjYmLc59cFEx9KZWJ2/lovCVr6TvfkdgZmVTdKlneWBFScsDKwFPAh8GfpXvPxfYs+AYllKrk7f7ReG881wSMrNyKazUExGPS/ouMBt4BbieVNp5PiIW5d3mAOv29HhJU4GpAOPHj29obLU6eSsXhYUL04Xh7LNh0SKP+zez8igs8UtaHdgD2BB4HrgY2Knex0fENGAaQEdHR8MXBp44seckXn1RmD0bzjhjyZKQE7+ZtboiSz07Ag9HxNMR8TpwKfABYLVc+gFYD3i8qAAGOj1DZeTPQQctXRLylA9m1uqKHNUzG9hW0kqkUs9kYAZwE/Bx0sieKcDlRZy8EdMzdC8Jgad8MLPWV1iLPyJuJ3Xi3kkayrkcqXRzDHCEpAdJQzrPLOL8jZqeodL6nzjRUz6YWTkUOo4/Ir4KfLXb5oeArYs8LyzZSduo6RmKOKaZ2WAr7ZQNRUzP4CkfzKwMFNHwATMN19HRETNmzGh2GGZmLUXSzIjo6L69rebq8YgcM7MSl3q68/w8ZmZJ2yT+nkbkVBK8V+Yys3bSNqWe/szP42GaZlZmbdPir3d+Hg/TNLOya5vED/XNz+Mav5mVXVsl/p5Ud+pWVuYyMyuztk787tQ1s3bUFp27vY3fd6eumbWj0rf4a7Xq3alrZu2o9In/vPPg1VchYunx++7UNbN2VOrE39UFZ52Vkj7Acsst3arvbaSPmVlZlbrG39mZ1suteOONpoViZjZklDrxT5qUWvkVEe7ANTMrdeKfOBFOPRWGD08XgBVWcAeumVmpa/wAU6fCu97lDlwzs4rSJ35wB66ZWbVSl3rMzGxpTvxmZm3Gid/MrM048ZuZtRknfjOzNuPE3wC9zf5pZjYUtcVwziJ5Tn8zazVu8S8jz+lvZq3GiX8ZVeb0HzbMc/qbWWtwqWcZeU5/M2s1TvwN4CkhzKyVuNRTII/2MbOhyC3+gni0j5kNVW7xF8SjfcxsqHLiL4hH+5jZUFWz1CNpPWBf4IPAOsArwH3A1cC1EeFVbHvh0T5mNlT1mvglnQ2sC1wFnATMA0YCmwI7Af8t6diIuGUwAm1FtUb7dHX5omBmzVGrxf+9iLivh+33AZdKGgGMLyascnPHr5k1U681/l6SfvX9CyPiwcaHVH7u+DWzZuo18UtaRdK3JP1c0v7d7vtJ8aGVlzt+zayZao3qOTt/vwTYV9IlklbI27at5+CSVpP0K0kPSJolaaKkNSTdIOnv+fvqy/IEWlGl4/frX3eZx8wGX63Ev1FEHBsRl0XE7sCdwO8kjenH8U8GrouIzYD3ALOAY4EbI2IT4Mb8c9uZOBGOO85J38wGX63O3RUkLVcZshkR35T0OHALMLqvA0taFdgB+FR+/EJgoaQ9gEl5t3OBTuCYAcZvZmb9VKvFfyXw4eoNEXEOcCSwsI5jbwg8DZwt6S5JP5M0ClgrIp7M+8wF1up31GZmNmC1RvUcHRG/7WH7dblM05flgfcCp0XEVsBLdCvrREQA0dODJU2VNEPSjKeffrqO05mZWT36nLJB0soDPPYcYE5E3J5//hXpQvCUpHH52ONIHwxbSkRMi4iOiOgYO3bsAEMwM7PuaiZ+SZVP7vZbRMwFHpP09rxpMvAX4ApgSt42Bbh8IMc3M7OBqTVlwzuB6cBnluH4hwEX5E/5PgR8mnSxuUjSwcCjwN7LcHwzM+unWqN6bgL2jIjbBnrwiLgb6OjhrskDPaaZmS2bWqWeO4D/O1iBmJnZ4KiV+HcHVpX07cEKxszMildrOOfiiJgKLBjEeMzMrGB9DueMiK8NRiBmZjY46hnHv231WP48a+c2xYZlZmZFqWfN3dNYstyzIG8zM7MWVE/iV55aAYA8aVvNtXrNzGzoqifxPyTpcEnD89cXSB/GMjOzFlRP4v8ssB3weP7aBphaZFBmZlacPks2ETEP2HcQYrFuurrSeryTJnnBFjNrnD4Tv6S3kVbS2pY0hXIX8KWIcLmnQF1dMHlyWox9xAgv0WhmjVNPqecXwEXAOGAd4GLgwiKDstTSX7gQFi9O3zs7mx2RmZVFPYl/pYj4eUQsyl/nAyOLDqzdTZqUWvrDhqXvkyY1OyIzK4t6hmVeK+lY0hTNAewDXCNpDYCIeLbA+NrWxImpvOMav5k1mqqG6Pe8g/RwjbsjIt7W2JCW1tHRETNmzCj6NGZmpSJpZkQsNTV+rYVY1omIJyJiw2JDMzOzwVSr1POzXM7pBK4Dbo2IRYMSlZmZFabXxB8Ru0gaCUwiLcjyXUmzSReB6yJi9uCEaGZmjVSzczciXiUnegBJGwI7Az+WtHZEbF18iGZm1kh1T7YmaRXgBdLonul4gRYzs5ZUzyd3DwFOAF4lDeeEQRrNY2ZmjVdPi/8oYIuImF90MGZmVrx6Prn7D+DlogMxM7PBUU+L/zjgj5JuB16rbIyIwwuLyszMClNP4v8p8DvgXuCNYsMxM7Oi1ZP4h0fEEYVHYmZmg6KeGv+1kqZKGidpjcpX4ZGZmVkh6mnx75e/H1e1LQAP5zQza0H1LL3oSdqGOC/RaGb9Uc8HuA7qaXtEnNf4cKy/vESjmfVXPaWe91fdHglMBu4EnPiHgJ6WaHTiN7Na6in1HFb9s6TVSHP12BBQWaKx0uL3Eo1m1pe6J2mr8hLguv8Q4SUazay/6qnxX8lbk7MtB2wOXFxkUNY/Eyc64ZtZ/epp8X+36vYi4NGImFNQPGZmVrB6avw3d98maXZEjC8mJDMzK1I9n9ztiRoahZmZDZqBJv7oexczMxuKei31SOptYjYBo4sJx8zMilarxr9yjftObnQgZmY2OHpN/BFxQiNOIGkYMAN4PCJ2k7Qh6QNgY4CZwIERsbAR5zIzs77VKvWcUuuB/ViB6wvALGCV/PNJwA8iYrqk04GDgdPqPJY1gCd1M2tvtUo9M5f14JLWA3YFvgkcIUnAh4H98y7nAsfjxD9oPKmbmdUq9Zxb/bOk0Xn7gn4c/4fA0bzVXzAGeD4iFuWf5wDr9vRASVOBqQDjx/sjA43iSd3MrM/hnJK2kHQXcD/wF0kzJb2zjsftBsyLiAG9c4iIaRHREREdY8eOHcghrAeVSd2GDfOkbmbtqp4pG6YBR0TETQCSJgFnANv18bgPALtL2oU0nfMqpNFAq0laPrf61wMeH1joNhCe1M3M6kn8oypJHyAiOiWN6utBEXEcebnGfLE4KiI+Keli4OOkkT1TgMsHELctA0/qZtbe6vnk7kOSviJpQv76MvDQMpzzGFJH74Okmv+Zy3AsMzPrp3pa/P8POAG4NP98S95Wt4joBDrz7YeArfvzeDMza5x6Zud8Djgc3vww1qiI+GfRgZmZWTHqGdXzC0mr5Lr+vaSRPf9ZfGhmZlaEemr8m+cW/p7AtaRlFw8sMigzMytOPYl/uKThpMR/RUS8jqdlNjNrWfUk/p8CjwCjgFskbQC4xm9m1qL6TPwRcUpErBsRu0REALOBfy0+NDMzK0KviV/SAZKWuj+SRZI2krR9seGZmVmj1RrOOQa4S9JM0kydT5OmXtgY+BAwHzi28AhtSPGUzmatr9bsnCdL+jFpGuUPAO8GXiHNrX9gRMwenBBtqPCUzmblUPMDXBGxGLghf1mb85TOZuVQz6geM8BTOpuVRT1z9ZgBntLZrCz6TPySVoiI17ptWyMini0uLBuqPKWzWeurp9Rzaf7kLgCSxuGav5lZy6on8V8GXCRpmKQJwG/IC6yYmVnrqWda5jMkjSBdACYAh0TEHwuOy8zMCtJr4pd0RPWPwHjgbmBbSdtGxPcLjs3MzApQq8W/crefL+1lu5mZtZBan9w9ofu2PHfPaK/AZWbWuvq7Atd9eAUuM7OW5hW4zMzajFfgMjNrM16By8yszdQzjv8U4JSqTY9K8gpcZmYtqtY4/gMi4vxu4/mreRy/mVkLqtXiH5W/e9y+mVmJ1BrH/9P8fanx/GZm1rpqlXpO6e0+gIg4vPHhmJlZ0WqN6qkssj4SeC/w9/y1JTCi8MjMzKwQtUo95wJI+hywfUQsyj+fDvx+cMIzM7NGq2cc/+rAKlU/j87bzMysBdWz5u63gLsk3USannkH4Pgig7Ly6+ry2r1mzVLPB7jOlnQtsE3edExEzC02LCuzri6YPBkWLoQRI9IC7k7+ZoOnZqlH0vKSlBP9ncBwYNygRGal1dmZkv7ixel7Z2ezIzJrL70mfkmfAeaRpmj4DHAj8HFguqRjBik+K6FJk1JLf9iw9H3SpGZHZNZeapV6vghsRPrk7ixgg4iYL2kl4A7gpOLDszKaODGVd1zjN2uOWol/YUQ8Bzwn6cGImA8QES9LWjg44VlZTZzohG/WLLUS/4qStiKVg0bk28pfIwcjODMza7xaiX8ub83AWX278rOZmbWgWp/cnbQsB5a0PnAesBZpxa5pEXGypDWAXwITSAu87J1LSmZmNghqjep5v6S1q34+SNLlkk7Jybsvi4AjI2JzYFvg85I2B44FboyITUgjhY5dtqdgZmb9UWsc/0+BhQCSdiB9gvc84AVgWl8HjognI+LOfPtF0sigdYE9gHPzbueS1vI1M7NBUqvGPywins239yGVai4BLpF0d39OImkCsBVwO7BWRDyZ75pLKgX19JipwFSA8ePH9+d0ZmZWQ60W/zBJlQvDZOB3VffVM8cPAJJGA5cAX4yIJRZpj4gg1f+XEhHTIqIjIjrGjh1b7+nMzKwPtRL4hcDNkuYDr5CnYpa0Manc0ydJw0lJ/4KIuDRvfkrSuIh4UtI40qeDzcxskPTa4o+IbwJHAueQ5uOPqscc1teBJQk4E5gVEdVDQa8ApuTbU4DL+x+2mZkNVM2STUTc1sO2v9V57A8ABwL3VvUJ/Bepk/giSQcDjwJ71x2tmZkts7pr9f0VEbeSPuXbk8lFndfMzGqrZwUuMzMrESd+M7M248RvpdHVBSeemL6bWe8Kq/GbDSYv52hWP7f4rRS8nKNZ/Zz4rRS8nKNZ/VzqsVLwco5m9XPit9Lwco5m9XGpx8yszTjxm5m1GSd+M7M248RvZtZmnPjNzNqME7+ZWZtx4jczazNO/GZmbcaJ38yszTjxm5m1GSd+M7M248RvZtZmnPjN+sGrfFkZeHZOszp5lS8rC7f4zerkVb6sLJz4zerkVb6sLFzqMauTV/mysnDiN+uHZq7y1dXli441hhO/WQtwx7I1kmv8Zi3AHcvWSE78Zi3AHcvWSC71mLWAodix7D6H1uXEb9Yimtmx3J37HFqbSz1m1m/uc2htTvxm1m/uc2htLvWYWb8NxT4Hq58Tv5kNyFDqc6iHO6Pf4sRvZqXnzuglucZvZqXnzuglOfGbWem5M3pJLvWYWem1Ymd0kX0STUn8knYCTgaGAT+LiG81Iw4zax+t1BlddJ/EoJd6JA0DTgV2BjYH9pO0+WDHYWY2VBXdJ9GMGv/WwIMR8VBELASmA3s0IQ4zsyGp6D6JZpR61gUeq/p5DrBN950kTQWmAowfP35wIjMzGwKK7pMYsp27ETENmAbQ0dERTQ7HzGxQFdkn0YxSz+PA+lU/r5e3mZnZIGhG4r8D2ETShpJGAPsCVzQhDjOztjTopZ6IWCTpUOA3pOGcZ0XE/YMdh5lZu2pKjT8irgGuaca5zczanadsMDNrM078ZmZtRhFDf6SkpKeBR/vxkDWB+QWF00xlfV7g59aKyvq8oDzPbYOIGNt9Y0sk/v6SNCMiOpodR6OV9XmBn1srKuvzgnI/N3Cpx8ys7Tjxm5m1mbIm/mnNDqAgZX1e4OfWisr6vKDcz62cNX4zM+tdWVv8ZmbWCyd+M7M2U6rEL2knSX+V9KCkY5sdT6NIOkvSPEn3NTuWRpO0vqSbJP1F0v2SvtDsmBpB0khJf5L05/y8Tmh2TI0maZikuyRd1exYGkXSI5LulXS3pBnNjqcopanx5yUd/wZ8hLS4yx3AfhHxl6YG1gCSdgAWAOdFxBbNjqeRJI0DxkXEnZJWBmYCe7b6702SgFERsUDScOBW4AsRcVuTQ2sYSUcAHcAqEbFbs+NpBEmPAB0RUYYPb/WqTC3+0i7pGBG3AM82O44iRMSTEXFnvv0iMIu0SltLi2RB/nF4/ipHKwuQtB6wK/CzZsdi/VemxN/Tko4tn0DaiaQJwFbA7U0OpSFyKeRuYB5wQ0SU4nllPwSOBt5ochyNFsD1kmbm5V9LqUyJ31qYpNHAJcAXI+KfzY6nESJicURsSVplbmtJpSjTSdoNmBcRM5sdSwG2j4j3AjsDn89l1tIpU+L3ko4tKtfALwEuiIhLmx1Po0XE88BNwE5NDqVRPgDsnuvh04EPSzq/uSE1RkQ8nr/PA35NKiGXTpkSv5d0bEG5E/RMYFZEfL/Z8TSKpLGSVsu3VyQNOnigqUE1SEQcFxHrRcQE0v/Z7yLigCaHtcwkjcoDDJA0CvgoULqRdFCixB8Ri4DKko6zgIvKsqSjpAuBLuDtkuZIOrjZMTXQB4ADSa3Gu/PXLs0OqgHGATdJuofUKLkhIkoz7LGk1gJulfRn4E/A1RFxXZNjKkRphnOamVl9StPiNzOz+jjxm5m1GSd+M7M248RvZtZmnPjNzNqME7+VgqTFeSjofZKurIyh78fjOyV15NvX9PfxvRxzK0ln5tvHSzqqAcf8raTVl/U41t6c+K0sXomILfPspc8Cnx/ogSJil/xp22X1X8ApDThOtZ8D/9HgY1qbceK3MuoiT9AnaWtJXXne+D9KenvevqKk6ZJmSfo1sGLlwXlO9jUlTaheA0HSUZKOz7cPz2sI3CNpevcA8idA3x0Rf+7hvs9IujbH0Cnp5Kp3K1vnfUZLOjvPDX+PpL3yw68A9mvUC2XtaflmB2DWSHldhsmkaSAgTZPwwYhYJGlH4H+AvYDPAS9HxDskvRu4s5+nOhbYMCJe66Us1EEPH/eXdChp+oY982MBVoqILfOEYGcBWwBfAV6IiHflx60OEBHPSVpB0piIeKafMZsBTvxWHivmKZDXJU3ZcUPevipwrqRNSFPuDs/bdyCXYSLinjy1Qn/cA1wg6TLgsh7uHwc83W3bQaSpw/eMiNertl+Y47hF0ir5QrIjaR4c8n3PVe0/D1gHcOK3AXGpx8rilTwF8gaAeKvG/3Xgplz7/xgwsh/HXMSS/yPVj90VOBV4L3CHpO6NqFd6ONe9wATSzLHVus+b0tc8KiPz8c0GxInfSiUiXgYOB47MyXhV3pqe+1NVu94C7A+Q58l/dw+Hewr4F0ljJK0A7Jb3Xw5YPyJuAo7J5xjd7bGzgI27bbsLOAS4QtI6Vdv3ycfdnlTeeYH0juXNDupKqSfPZro28Eit18GsFid+K52IuItUitkP+DZwoqS7WLK0eRowWtIs4GuktX67H+f1fN+fSIm4Mq3yMOB8SfeSkvkp3UcBRcQDwKqVaX6rtt8KHAVcLWnNvPnVHN/pQGXm1W8Aq+cO3z8D/5q3vw+4Lc9GazYgnp3TrCCSvgS8GBG9rksrqRM4KiJm1HnMk4ErIuLGxkRp7cgtfrPinAa81uBj3uekb8vKLX4zszbjFr+ZWZtx4jczazNO/GZmbcaJ38yszTjxm5m1mf8Fn1cCB5Umo5AAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "Bulge/Disk Decompositions\n",
      "Radius     SBdisk     SBbulge  \n",
      "  0.06074   78.34838    0.00000\n",
      "  0.07297   71.71831    0.00000\n",
      "  0.08752   71.58632    0.00000\n",
      "  0.10503   72.58218    0.00000\n",
      "  0.12590   76.07291    0.00000\n",
      "  0.15121   82.95260    0.00000\n",
      "  0.18137   81.06437    0.00000\n",
      "  0.21764   73.52415    0.00000\n",
      "  0.26109   75.51444    0.00000\n",
      "  0.31339   80.54341    0.00000\n",
      "  0.37602   79.80499    0.00000\n",
      "  0.45131   78.34838    0.00000\n",
      "  0.54158   67.05482    0.00000\n",
      "  0.64976   64.21404    0.00000\n",
      "  0.77989   62.34908    0.00000\n",
      "  0.93574   53.21451    0.00000\n",
      "  1.12301   45.33462    0.00000\n",
      "  1.34761   39.92363    0.00000\n",
      "  1.61713   33.66899    0.00000\n",
      "  1.94044   23.59555    0.00000\n",
      "  2.32848   15.87925    0.00000\n",
      "  2.79413    8.51981    0.00000\n",
      "  3.35300    6.43324    0.00000\n",
      "  4.02365    4.53346    0.00000\n",
      "  4.82842    2.68429    0.00000\n",
      "  5.63319    0.99147    0.00000\n"
     ]
    }
   ],
   "source": [
    "decomp_data = import_bulge_disk_decomps(datapath='./sparc/BulgeDiskDec_LTG/',gal_ID=gal_id)\n",
    "\n",
    "plot_single_galaxy_decomp(decomp_data, gal_id)\n",
    "\n",
    "# Print the bulge/disk decompositions for this galaxy\n",
    "dec_labels = decomp_data.dtype.names\n",
    "print('\\nBulge/Disk Decompositions')\n",
    "print('  '.join(['{:<9}'.format(item) for item in dec_labels]))\n",
    "for dec_line in decomp_data:\n",
    "    dec_line_zip = zip([item for item in dec_line],dec_labels)\n",
    "    dec_print_arr = ['{:>9.5f}'.format(item) for item,name in dec_line_zip]\n",
    "    print('  '.join(dec_print_arr))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "0af373ec-8dd5-4c89-944b-decebdc37f29",
   "metadata": {},
   "outputs": [],
   "source": [
    "def plot_btfr_data(btfr_data, sample_data, include_irregulars=False):\n",
    "    sp_types = [1, 2, 3, 4, 5, 6, 7]\n",
    "    if include_irregulars:\n",
    "        sp_types += [9, 10, 11]\n",
    "    spirals = np.array([[btfr_data[gal_id]['Vf'],btfr_data[gal_id]['log_Mb'],btfr_data[gal_id]['e_Vf'],btfr_data[gal_id]['e_log_Mb']] for gal_id in btfr_data if sample_data[gal_id]['T'] in sp_types])\n",
    "    lentics = np.array([[btfr_data[gal_id]['Vf'],btfr_data[gal_id]['log_Mb'],btfr_data[gal_id]['e_Vf'],btfr_data[gal_id]['e_log_Mb']] for gal_id in btfr_data if sample_data[gal_id]['T'] == 0])\n",
    "\n",
    "    plt.scatter(spirals[:,2],spirals[:,1],color='k', alpha=0.3)\n",
    "    plt.scatter(lentics[:,2],lentics[:,1],color='b', alpha=0.4)\n",
    "    plt.xlabel('log(M_b)  (solMass)')\n",
    "    plt.ylabel('Flat Rotation Velocity (km/s)')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "3c06332e-92f7-46a4-9c3e-7bf29cbed736",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYgAAAEHCAYAAAC0pdErAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/MnkTPAAAACXBIWXMAAAsTAAALEwEAmpwYAAAs4ElEQVR4nO3de3xcdZ3/8denLb2kaZq2iaHtQAsYUawUNeEiQhG5VERrWW+4KnVZu7ruY3F/u/4W97e/XS/Lqquu911/RRB1hV0VY3DlIkUFlpWSttK1UCUtFUlo05Q2TUPSFsjn98c5E6bpZHIymTNzZub9fDzmMTNnzpnzOen0fM75Xs3dERERGW1KqQMQEZFkUoIQEZGslCBERCQrJQgREclKCUJERLKaVuoACqmhocGXLl1a6jBERMrGpk2b9rp7Y7bPKipBLF26lI0bN5Y6DBGRsmFmT4z1WWxFTGZ2o5ntMbOtGcveZmaPmNmwmbXk2PZ3ZvZrM3vYzHTGFxEpgTjrIG4CVo5athW4Argvwvavc/cz3H3MRCIiIvGJrYjJ3e8zs6Wjlm0DMLO4disiIgWS1FZMDvzUzDaZ2dpcK5rZWjPbaGYbe3t7ixSeiEjlS2qCeK27vwp4A/AhMzt/rBXdfZ27t7h7S2Nj1op4ERHJQyJbMbl7d/i8x8zagDOJVm8hIlWmsxPWr4fubli8GC66CJqbSx1VZUjcHYSZzTazOenXwCUEldsiIkfp7IQbboCDByGVCp5vuCFYLpMX2x2Emd0CXAA0mFkX8PfAPuArQCPwEzN72N0vNbNFwDfc/TKgCWgLK7KnATe7+51xxSml0dXVRUdHB729vTQ2NtLa2koqlSp1WFJm1q+H+fOhvj54n35ev153EYUQZyumK8f4qC3Luk8Bl4WvHweWxxWXlF5XVxft7e3U19fT1NTEwMAA7e3trFq1SklCJqS7O7hzyFRXFyyXyUtcEZNUvo6ODurr66mrq2PKlCnU1dVRX19PR0dHqUOTMrN4MfT3H72svx8WLSpNPJVGCUKKrre3l9ra2qOW1dbWombKMlEXXQT79kFfHwwPB8/79gXLZfIS2YpJKsvo+oYpU6YwMDBAXV3dyDoDAwOombJMVHMzXH31C62YFi2C1atV/1AoShASq2z1Dbt378bMWLp0KbW1tQwMDNDX18eKFStKHa6UoeZmJYS4KEFIrDLrGwDq6uo46aSTGBwcpKamhp6eHhobG1mxYoUqqEUSRgmizCW9uWhvby9NTU1HLautreWZZ55h9erVJYpKRKKo+gSR9BNsLuXQXLSxsVH1DSJlqqpbMaVPsIODgzQ1NTE4OEh7eztdXV2lDi2Scmgu2traSl9fH/39/QwPD9Pf309fXx+tra2lDk1ExlHVdxDZysfTy5NyBQ5H3+WYGWbG8PAwmzdvprW19air89raWnp6emLb/0TvslKpFKtWraKjo6Ng9Q3lfNcnUk6qOkGMVT5e6BPsZGQWI02dOpX7778fd+f8889n5syZ3HfffVxwwQU0NDQAhS++KUQxViqVKtgJvByK1UQqRVUXMaXLxzMlrXw88y5nx44dLFiwgIaGBnbs2MHy5csxM7Zs2RJb8U3SirGSFo9IJavqBFEO5eOZvY4PHDjArFmzmDlzJgcOHKChoYHzzjuPw4cP09PTQ01NTcGvpJPW6zlp8YhUsqouYoqjfLzQMlsBzZ07l6GhIQDmzp0LwMyZM7nkkktiazKatFZISYtHpJJVdYKAwpaPx6G1tZX29nYATjnllKPqINJ3PHH2QM7cfxJ6PSctHpFKZu5e6hgKpqWlxTdu3FjqMApurFZMxWrBk7RWQ0mLR6Scmdkmd2/J+pkShJQ7JQyR/OVKEFVdSS3lr9w7O4okWdXXQSRRtV4R53Pc5dLZUaQc6Q4iYQp5RdzV1UVbWxvr1q2jra0t0VfV+R63mr2KxEcJImEK1RGs3Ipe8j3ucujsKFKulCASplBXxOXW4zjf4y6Hzo4i5UoJImEKdUVcbkUv+R53urNjevKhOHqTi1QrVVInTKE6go3X4zhpFeGTOe6kd3YUKVfqBxGTyZyAC3Hyzhz1NPOEu2rVKoAxP0tSB7hUKkVXV1dikphIJVJHuSLLdXIu5glurETT1tbG4ODgUXcX/f391NTUJGYa0KT8DUUqXa4EoSKmGCSlbf5YRS/lMA9GUv6GItVMCSIGSToBZ7uLyFY/8cQTT9Dd3c26detiKc6ZaLFZMf6GSauHEUma2FoxmdmNZrbHzLZmLHubmT1iZsNmlvWWJlxvpZn91sy2m9m1ccUYl6S0zR+rL0QqlTqqaejOnTt54IEHWLx4cSx9JvLpkxH337Dc+omIlEKczVxvAlaOWrYVuAK4b6yNzGwq8DXgDcBpwJVmdlpMMcYiKW3zx+oL0dXVdVTT0O7ubs4991xOOumkWPpM5NMnI+6/Ybn1ExEphdgShLvfB+wbtWybu/92nE3PBLa7++PufgT4d2BVTGHGIilt83P1hUilUqxevZq1a9dyyimnsGTJkqzrxR3HWOL+G5ZbPxGRUshZB2FmM4HLgfOARcAQwV3AT9z9kZhiWgw8mfG+Czgrpn3FJglt86POvhb3LG35fH/c9QOamU5kfGMmCDP7OEFy+AWwAdgDzAReAnw6TB5/6e7/U4Q4x2Rma4G1ACeeeGIpQ0mcqJ3PCjlLW7YT+0S/P7OJa1NTEwMDA7S3txf0DiJKTKrElqSL+zeaq4jpIXd/tbv/pbvf7O7r3f0/3f2f3f1NwB8C0wsWyQu6gRMy3qfCZVm5+zp3b3H3Fl39HS1qMU2hinPGqvgFJvT9xagfGO+YVYktSVeM3+iYdxDu/pPRy8xsClDr7v3uvofgrqLQOoBmMzuJIDG8E3hXDPupClGLugpRJJar78Lq1asjf3+xmgnnOmb1w5CkK8ZvdNxKajO72czqzGw2Qf3Do2b2kQjb3QL8EjjVzLrM7GozW21mXcA5wE/M7K5w3UVmdjuAuz8H/BlwF7AN+F6M9R1SQIWq+E1CM2FVYkvSFeM3GqUV02nu3g+8BbgDOAl4z3gbufuV7r7Q3Y9z95S73+DubeHrGe7e5O6Xhus+5e6XZWx7u7u/xN1Pcffr8js0KbZCndiT0Ew4CUlKJJdi/EajJIjjzOw4ggRxm7s/C1TOAE4JVE4zwWUq1Ik9Cc2Ek5CkJH7l+n8NivMbHXewPjP7c+CvgS3AG4ETgX9z9/MKFkWBJGWwvsmYyCB1SWxlk8SY8lVJxyLHqoQBIQvxG81rNFczOwd40EetYGYGTA3rChKlEhJE1JFWK+HHHZVO1BKHchjVuBjyHc31vcDXzOwx4E7gTnffHSaMxCWHcjHeyS5qC55CtWBI+sm3GH0ipDolaVDNpBqzDsLdP+jurwI+BswDbjKzX5rZP5rZ+eGYSTIBUdotR614KkQLhnJo668xkyQuaogwvnErqd39N+7+BXdfCVwI/BfwNoLe1TIBUU52USueCvHjLoeTr5qbSlzUEGF8kQbrM7N5ZnY68DJgN/DNscqsZGxRTnZRW/CM/nHv3LmTe++9l8cee2ykNcZ4LTTK4eSrqzyJSxJayyXduBMGmdkngTXA48BwuNgJ7iZkAqIOEBelV3P6x93R0cG2bdt4/PHHWbZsGUuWLGFgYICbbroJM2Pp0qVjlt2Xw4B1hRwnSmS0JAyqmWRRZpR7O3BKOPS2TEKhT3bpH3dbWxsnnHDCURXWe/fuBeD0008fWQZHV2IX4+Q72UrwzETY09NDY2MjK1as0H9qkSKIkiC2AvXEM+5SVYnrZJetNcaRI8fm89EtNOI++RaqBZKu8kRKI0qC+BTwq3Dq0MPphe7+5tiiqmCFONmNvio3s2OKiqZPP3ag3XyLs/KlAe9EyluUBPEt4DPAr3mhDkJKJNtV+Z49e3B3TjrppJGiooaGBsyM/v7+kpXdZ97Z7N27l87OTvr6+hgeHk5cfwsROVaUBDHo7l+OPRKJJNtV+dKlSxkaGhppjdHY2MiaNWtG1i9V2X26EvzIkSNs2LCB2bNnM2PGDNxdnd1EykCUBHG/mX0KuI2ji5g2xxaVjGms3p/PPPNM1uEBxjsBx9mTOl0J3tnZSU1NDQBDQ0OcddZZTJ8+XUVNIgkXpR/EK4GzgX8EPh8+PhdnUDK2QvYLiLsndboS/PDhwxw6dIgZM2Zw1lln0dDQkLj+FiJyrCh3EG9w90OZC8xsQUzxVK2oV/KFbJpajErkVCrFJZdccsygaGMltaSPDSVSTaLcQdxqZiOJxMyOB34aX0jVZyJX8oXs/VmsntRRhzQoh7GhRKpJlDuIHwHfN7O3AicQ1EX8VZxBVZuJXMlHvcKOsl6xelJH7W+hZrEiyTJugnD3681sOkGiWAr8ibv/d8xxVZWoww53dXVx0003sXfvXo4cOcL06dN55JFHWLNmzVEn0GxNYa+//hfMmPFGBgfnsXgxXHTR5IqrJloUFKW/hYZfFkmWMYuYzOx/pR/ATIKZ5B4Gzg6XSYFErXi+44472L59O1OnTmX+/PlMnTqV7du3c8cddxy13uhRWgcGjmfz5lfym990kUrBwYNwww0wNHR0cdXQ0BDTp0/n9ttvzzn9YlxFQRqYTyRZctVBzMl41AI/BLZnLJMCiVpGv3nzZubNm0dNTQ1mRk1NDfPmzWPz5qNbHI+uW3j44Ubmz3eGh/czZQrU18P8+bB+fXBlv3r1ai677DIOHz5MTU3NuCf9qMOET3S+Xw2/LJIsYxYxufvHixlINYtaRu/uBDO+vsDMGD1t7Oi6hX37ZjJ79n5qal7I63V10N39wjYTKf+PUhSUzzhMGphPJFnGTBBmdj3wJXffmuWz2cA7gMPu/t0Y46samWX06fL922+/faR8H2DKlCn88pe/pLGxkSVLljB9+nSefvppzjvvvKO+a3TdQk3Nfvbte45zznnJyDr9/bBo0QvbTKT8P0rldr4VzhqYTyQ5clVSfw34OzN7BcGIrr0EdRHNQB1wI6DkUGDZrrzTczu89KUvpb+/n97eXrZs2UJzczMveclLWLly5VHfMfpKvKVlFo8+ejZmdQwPB8lh3z7I7Hg9kRZNUSq3H3vssZHiorlz59Lc3Mz8+fNV4SxSRnIVMT0MvN3MaoEWYCEwBGxz998WJ7zqk+3KO3Nuhzlz5tDZ2cnu3btpaGjgqquuynrFPfpKvLMzqHPo7g7uHFavhubmF9afSIum8YqCurq62LFjB9OmTWPevHkcOnSIDRs2cNppp3HCCScU7G8lIvGK0sx1APhF/KEIjD+3Q0NDAw0NDQwPD9PT0xO5OKa5+eiEMNpEy/9zFQV1dHSwbNkyHn30UQ4dOsSsWbMYGhpi69atXHHFFZHiFZHSi9JRTgpovP4D2Yp6cs3tUMihKQpV/t/b28uSJUtG7nb2798/0tJJ9Qsi5SPKUBtSIFH6D2Rr6tnQ0EBjY+MxzT9TqVQih6ZIJ7mGhgbOOeccVq5cySte8QpOPfXUksYlIhMzboIIK6knzMxuNLM94Ux06WXzzexuM+sMn+eNse3zZvZw+Lgtn/0nUZT+A9nGWlqzZg1XXXXVMeMvdXV1ReqPUGzqzyBSGaIUMf2Lmc0AbgK+6+4HIn73TcBXgW9nLLsWuMfdP21m14bv/zrLtkPufkbE/ZSNqE1JxyrqGb3s9ttvT+TQFOrPIFIZolRSn2dmzcAfAZvM7CHgm+5+9zjb3WdmS0ctXgVcEL7+FkHld7YEUZEaGxt54okn2L17NwcOHGDu3Lkcf/zxebfsKfRge0mszxCR0olUB+HuncDfEpzMVwBfNrPfmNlEm6Q0ufuu8PVuoGmM9Waa2UYze9DM3pLrC81sbbjuxqRPQJNKpXjggQfo6+tj7ty59PX18cADD+R9Ii1kUY6G2haR0aLUQZxuZl8AtgEXAm9y95eFr7+Q7449GB/Cx/h4ibu3AO8Cvmhmp+T4nnXu3uLuLUkf1K2rq4tzzz2X+vp6Dhw4QH19Peeee27eJ+FCzg0RdXwlEakeUeogvgJ8A/gbdx9KL3T3p8zsbye4vx4zW+juu8xsIbAn20ru3h0+P25mvyCY9nTHBPeVOOnmnyeddNLIsnR/hnwVsmlqEuszRKR0ohQxtbn7dzKTg5ldA+Du35ng/m4DrgpfXwW0j17BzOaFleKYWQNwLvDoBPeTSEkezjrJsYlIaURJEO/NsmzNeBuZ2S3AL4FTzazLzK4GPg1cbGadwEXhe8ysxcy+EW76MmCjmW0Bfg582t0rIkEkuflnkmMTkdKw0UNFj3xgdiVBHcBrgfszPpoDDLv76+MPb2JaWlp848aNpQ4jp0K2FCq0JMcmIvEws01hne8xctVB/DewC2gAPp+x/CDwP4ULr7oUuvlntTRNzRxsMD1laq6xpURk8sa8gyhH5XAHUUiZQ4NnjsCab0umYoua3Do7gylS588PJjpKD1d+9dVKEiKTlesOItec1P8VPh80s/6Mx0Ez648rWImunJumTqTfxfr1QXKor+eYKVNFJD5jJgh3f234PMfd6zIec9y9bqztpHhGzz0NQdPUpHcYhIklt+7u4M4hU10dPPVUkYIVqVJROsqdbWZzMt7PMbOz4g1LoijnpqkTSW6LFwfFSplGT5kqIoUXpZnrvwKZZ6FnwmVSYuXcNHUiye2ii4I6h74+GB4OnvftC5aLSHyiJAjzjJpsdx9GEw0lQiGH2ii2iSS35uagQnrOnKC4ac4cVVCLFMO4rZjM7IcEo66m7xr+FHidu78l1sjyUG2tmMqd+l2IlF6+/SDSPgB8mWA0V4D1wNoCxSZVLMn9LkQk2nwQe4B3FiEWERFJkHEThJmlCEZ0PTdcdD9wjbtrooACUVGLiCRRlErqbxKMwroofPw4XCYFoIl6RCSpotRBNLp7ZkK4ycw+HFM8VSezwxgw8tzR0RHLXYTuVkQkqih3EE+b2bvNbGr4eDfwdNyBVYti9obW3YqITESUBPFHwNsJ5pDeBbwVeF+cQVWTYvaGLuexm0Sk+MZNEO7+hLu/2d0b3f1F7v4Wd/99MYKrBsXsDV3OYzeJSPGNWQdhZl8BxuxF5+5/HktEVSbdG7qjo4Oenh4aGxtZsWJFLPUC6buVuoyR78pl7CYRKb5cldTqklwkxeow1traSnt7MA145vwRK1asiH3fIlJ+xkwQ7v6tzPdmVuPug/GHVB1K0ZqomHcrIlL+onSUOwe4AagFTjSz5cCfuPufxh1cpcqcCa6pqYmBgQHa29uLMtCehrcQkaii9IP4InApQWc53H2LmZ0fZ1CVLkrfB/VXEJFSizRst7s/aWaZi56PJ5zq0NvbS1NTEwB79+5l06ZNbN++ncHBQZ5++mmWL1/OQw89VJI7DBGRtCgJ4kkzew3gZnYccA2wLd6wKlu6NdGRI0e455572LVrF8cddxxz587l/vvv52c/+xlnn3120XpXi4hkM2Y/iDAZQDDc94eAxUA3cEb4XvKU7vuwZcsWDhw4wPTp0wE4+eSTaWhoYPfu3Tw1asJl9VcQkWLL1VGu28y+ASwH3u3uTWFHuXe7u4bamIR0a6JDhw7R39/P7NmzWbx4MU899RTbtm1j37597Ny586ht1F9BRIotV4J4GdBBMFHQk2b2JTM7qzhhVb5UKsWll17KK1/5SubPn093dzfPPvvsSFFTV1cXO3fuLLu5pkWkcoyZINz9aXf/f+7+OuBM4HHgi2a2w8yuK1qEFay1tZWGhgY6OzsBcHcOHjzIiSeeyCWXXEJ3d3fZzTUtIpUjymB9uPtTBH0h/hU4CPxxlO3M7EYz22NmWzOWzTezu82sM3yeN8a2V4XrdJrZVVH2V25SqRRr1qwZOfEfOXKE0047jQsvvJAzzjiDU045hbVr17J69WolBxEpOnMfc7glzGwm8CbgSuA1wJ3AvwN3u/u4TV3D/hIDwLfdfVm47J+Afe7+aTO7Fpjn7n89arv5BEN9tBCMB7UJeLW778+1v5aWFt+4sfxGCGlra2NwcPCoMZL6+/sZHBykoaFBfSFEJDZmtsndW7J9lqsV083A7wmG+v4usNTd17j7nVGSA4C73wfsG7V4FZAexuNbwFuybHopQRLaFyaFu4GVUfZZjrKN6Lpz50727NmjuRtEpGRy9YO4k2BIjYMF3meTu+8KX+8GmrKssxh4MuN9V7isImWOkbRt2zb279/PU089xYIFC1i0aNHI3A2gvhAiUjy5Kqm/HUNyGL0PJ8eQ4lGY2Voz22hmG8u5n0AqlaK1tZU5c+Zw+umnM2fOHKZMmcKGDRvYu3cvoL4QIlJckYbaKLAeM1vo7rvMbCGwJ8s63cAFGe9TwC+yfZm7rwPWQVAHUdhQJyaf8ZMyt9mxYweLFy8ement8OHDzJ49m87OThoaGiq2L4TGnRJJpkitmArsNiDdKukqoD3LOncBl5jZvLCV0yXhsoLr6uqira2NdevW0dbWlncZfz7zPY/epre3l0ceeYS9e/fS3NzMM888w/Dw8Ej9RCX2hdA82SLJFSlBmNlrzOxdZvbe9CPidrcAvwRONbMuM7sa+DRwsZl1AheF7zGzlrDnNu6+D/gkQUe9DuAT4bKCKuTJKZ/5nkdvc/zxxzN16tSRO4azzjoLd2d4eLhi+0JonmyR5IoyH8R3gFOAh3lhFFcHvj3etu5+5RgfvT7LuhvJ6F/h7jcCN463j8mIMux2VJkjtKbV1tbS09MTeZvm5mYefPBBenp6GB4eZvr06TQ3N+edGMqh6Cafv5uIFEeUOogW4DTP1WGiTBXy5JTPfM+jt2loaODlL3/5SA/qycz4VspJiSZC82SLJFeUIqatwPFxB1IK6ZNTpnxPTtn6MoxXZ5C5zZ49e7jnnnt48MEHefGLX8xll102qR7U5VJ0k8/fTUSKI0qCaAAeNbO7zOy29CPuwIqhkCendF+GmpqayOMnpbcZGhpi/fr1ALz+9a+npqZm0hW1vb291NbWHrUsic1k8/m7iUhxRCli+ljcQZRKZge1yRbppL9votumUikWLFjA5ZdfflQxC0yuU1w5Fd1onmyRZBo3Qbj7vWbWBKQvqx9y92x9F8pSEk5OcVTUtra20t7ePvJdAwMD9PX1sWLFiknFKiLVI0orprcDnyXoqGbAV8zsI+7+g5hjqxr5XO2P10Kp0HdHIlJ9ohQx/R+gNX3XYGaNwHpACaJAJnq1H7WFUhLujkSkfEWppJ4yqkjp6YjbSUQTragtlxZKIlLeotxB3GlmdwG3hO/fAdweX0jVaSJX+1HqLMqhk5yIJNu4dwLu/hGCwfBODx/rRk/wI8U1Xv8NjW8kIoUQaTRXd78VuDXmWCSi8eosCjmEiIhUr1wzyv1X+HzQzPozHgfNrL94Icpo49VZlEsnORFJtjHvINz9teHznOKFI1HlqrMop05yEp3qlaTYxq2DCEdzHXeZJIfGN6o8qleSUojSXPXlmW/MbBrw6njCkULQ+EaVR02bpRTGLGIys48CfwPMyqhzMOAI4RSfklzqJFdZNG+GlEKuOohPAZ8ys0+5+0eLGFNRqVxXyoHqlaQUovSD+Gg4N/SZZnZ++lGM4OKmcl0pF6pXklKIUkn9x8B9wF3Ax8Pnj8UbVnGoXFfKheqVpBSidJS7hmCo7wfd/XVm9lLgH+MNqzgqrVxXxWWVTfVKUmxRWjEdcvdDAGY2w91/A5wab1jFUcgpR0tNxWUiUmhREkSXmdUDPwLuNrN24Ik4gyqWSirXVXGZiBRalBnlVocvP2ZmPwfmAnfEGlWRVNKkOpVWXCYipRdlRrnvuPt7IJh+NL0MeE/MsRVFXOW6xa4PUDNIESm0fHpST0U9qXNK1wc8+eST7Ny5kx//+Md84hOfiLW4p5KKy0QkGXKN5vpRMzsInJ4xiutBYA/QXrQIy1BHRwfPP/88jz76KEeOHGHhwoVMmzaN66+/PrZKYzWDFJFCq/qe1HHo7e1l165dzJ49m5qaGgDmz5/Prl27Yp2TQc0gRaSQovakfrOZfS58XF6MwMpZY2MjPT09zJw5c2TZ0NAQjY2NmpNBRMpGlJ7UnyLoLPdo+LjGzCbVUc7MrjGzrWb2iJl9OMvnF5jZATN7OHz83WT2V2ytra1MmzaN/fv34+4MDg7yzDPPsGjRIlUai0jZiFJJ/UbgYne/0d1vBFYCed9FmNky4P3AmcBy4HIze3GWVe939zPCxyfy3V8ppFIp3v/+9/Pcc8+xa9cujjvuOE477TSmTp2qSmMRKRuR5qQG6oF94eu5k9zny4AN7j4IYGb3AlcA/zTJ702U1tZWFi5cqKEvRKRsRUkQnwJ+FXaSM+B84NpJ7HMrcJ2ZLQCGgMuAjVnWO8fMtgBPAX/l7o9k+zIzWwusBTjxxBMnEVbhqdJYRMqZufv4K5ktJBiwD+AhYLq7/z7vnZpdDfwp8AzwCHDY3T+c8XkdMOzuA2Z2GfAld28e73tbWlp848ZsuabyaaA+EcmHmW1y95Zsn+WsgzCzc8zsrcDz7n4b8Dvgn4EHJhOQu9/g7q929/OB/cBjoz7vd/eB8PXtwHFm1jCZfVYyDdSXXVdXF21tbaxbt462traq/3uITFSujnKfBW4E/gD4iZn9A/BTYAMw7tV8Lmb2ovD5RIL6h5tHfX68mVn4+swwzqcns89KpoH6jqWkKTJ5ueog3gi80t0Pmdk84Elgmbv/rgD7vTWsg3gW+JC795nZBwDc/evAW4EPmtlzBPUU7/QoZWFVSgP1HSszaQIjz3F2VMxFRYBSjnIliJF5INx9v5l1Fig54O7nZVn29YzXXwW+Woh9VQMN1HesJCXN9N1MfX09TU1NDAwM0N7erqFQJPFy1UGcbGa3pR/ASaPeS0JooL5jJWkyKBUBSrnKdQexatT7z8cZiIxtvOKJSprXolBaW1tpbw/GlKytrWVgYIC+vj5WrFhR9FiSdDcjMhG5Buu7t5iBSHZRiyfU5+JoSUqaKgKUchW1J7WUSNIqW8tJUpJmku5mRCZCCSLhMosn9u7dS2dnJ319fQwPD6slTJlI0t2MyEREmXL0be7+/fGWSTzSxRNHjhxhw4YNzJ49mxkzZuDuaglTRpJyNyMyEVFGc802WZAmEBpHoXrxplsobdmyZWTyoaGhIZYvX66WMCISqzHvIMzsDQQD6S02sy9nfFQHPBd3YOWskO3e08UT1113HWZGfX09y5Yto6GhgeHh4UgtYdRJS0TykauI6SmCUVbfDGzKWH4Q+Is4gyp3ha5YTqVSXHLJJQwODk64JYw6aYlIvnI1c90CbDGzm9392SLGVPbSFcvpSuUDBw6MdI7KV74tYdQKSkTyFaUOYqmZ/cDMHjWzx9OP2CMrY42NjTzxxBNs2LCBw4cPU19fz4EDB9ixY0fedRHpoqaamhp6enqoqamJdBfQ29tLbW3tUctqa2s1N7aIjCtKM9dvAn8PfAF4HfA+oiWWshBH+Xxrayt33HEH06ZNY9asWQwNDTE8PMyyZcsmdeWeT0sYddISkXxFOdHPcvd7CCYXesLdP0Yw0mvZi2tI6FQqxcknn0xdXR379+9nxowZnHXWWSxZsqToV+4ap0lE8hXlDuKwmU0BOs3sz4BuoHacbcpCnOXzp5566jGVyv39/UW/clcnLRHJV5QEcQ1QA/w58EngQuCqOIMqljgHUUvS8ArqpCUi+Rg3Qbh7uifWAEH9Q8WYbPl8rvoLXbmLSLnL1VHux8CYs7i5+5tjiaiIJnOVH6V/ga7cRaSc5bqD+FzRoiiRyVzlq3+BiFS6XAlip7v/vmiRlEi+V/maBEZEKl2uZq4/Sr8ws1vjD6W8JGlKSxGROORKEJbx+uS4Ayk36l8gIpUuV4LwMV4L+Q99ISJSLnLVQSw3s36CO4lZ4WvC9+7udWNvWh3USklEKlmu0VynFjMQERFJlooZdE9ERApLCUJERLJSghARkayiDNZXcGZ2DfB+ggrv6939i6M+N+BLBHNiDwJr3H1zseNMMs0zLSJxK/odhJktI0gOZwLLgcvN7MWjVnsD0Bw+1gL/WtQgEy6ueSxERDKVoojpZcAGdx909+eAe4ErRq2zCvi2Bx4E6s1sYbEDTarMcaCmTJkyMt91R0fH+BuLiERUigSxFTjPzBaYWQ1BMdIJo9ZZDDyZ8b4rXCZonmkRKY6iJwh33wZ8BvgpcCfwMPB8vt9nZmvNbKOZbayWE6TGgRKRYihJKyZ3v8HdX+3u5wP7gcdGrdLN0XcVqXBZtu9a5+4t7t5SLSdIjQMlIsVQkgRhZi8Kn08kqH+4edQqtwHvtcDZwAF331XkMBNL40CJSDGUpJkrcKuZLQCeBT7k7n1m9gEAd/86cDtB3cR2gmauFTXVaSFoHCgRiVtJEoS7n5dl2dczXjvwoaIGJSIiR1FPahERyUoJQkREslKCEBGRrJQgREQkKyUIERHJSglCRESyUoIQEZGslCBERCQrJQgREclKCUJERLJSghARkaxKNVhfRdI80SJSSZQg8pAtEQC0t7dTX19PU1MTAwMDtLe3axhuESlbShAT1NXVlTURTJ8+fWSeaGDkuaOjQwlCRMqS6iAmqKOjYyQRTJkyhbq6Ourr69m8ebPmiRaRiqIEMUG9vb1ZE4G7a55oEakoShAT1NjYmDURtLS0aJ5oEakoShAT1NramjURrFy5UvNEi0hFUSX1BKVSKVatWkVHRwc9PT00NjayYsWKkUSghCAilUIJIg+pVEqJQEQqnoqYREQkKyUIERHJSglCRESyUoIQEZGslCBERCQrc/dSx1AwZtYLPJHn5g3A3gKGk2Q61spUTccK1XW8cR7rEnfPOuRDRSWIyTCzje7eUuo4ikHHWpmq6Vihuo63VMeqIiYREclKCUJERLJSgnjBulIHUEQ61spUTccK1XW8JTlW1UGIiEhWuoMQEZGslCBERCSrqk8QZrbSzH5rZtvN7NpSx1NoZnajme0xs60Zy+ab2d1m1hk+zytljIViZieY2c/N7FEze8TMrgmXV9zxmtlMM3vIzLaEx/rxcPlJZrYh/D3/h5lNL3WshWJmU83sV2b2n+H7ijxWM/udmf3azB42s43hspL8hqs6QZjZVOBrwBuA04Arzey00kZVcDcBK0ctuxa4x92bgXvC95XgOeAv3f004GzgQ+G/ZyUe72HgQndfDpwBrDSzs4HPAF9w9xcD+4GrSxdiwV0DbMt4X8nH+jp3PyOj70NJfsNVnSCAM4Ht7v64ux8B/h1YVeKYCsrd7wP2jVq8CvhW+PpbwFuKGVNc3H2Xu28OXx8kOJkspgKP1wPpuW+PCx8OXAj8IFxeEccKYGYp4I3AN8L3RoUe6xhK8huu9gSxGHgy431XuKzSNbn7rvD1bqCplMHEwcyWAq8ENlChxxsWuTwM7AHuBnYAfe7+XLhKJf2evwj8b2A4fL+Ayj1WB35qZpvMbG24rCS/Yc0oV+Xc3c2soto6m1ktcCvwYXfvDy42A5V0vO7+PHCGmdUDbcBLSxtRPMzscmCPu28yswtKHE4xvNbdu83sRcDdZvabzA+L+Ruu9juIbuCEjPepcFml6zGzhQDh854Sx1MwZnYcQXL4rrv/MFxcsccL4O59wM+Bc4B6M0tf+FXK7/lc4M1m9juCYuALgS9RmceKu3eHz3sIEv+ZlOg3XO0JogNoDltDTAfeCdxW4piK4TbgqvD1VUB7CWMpmLBc+gZgm7v/c8ZHFXe8ZtYY3jlgZrOAiwnqXH4OvDVcrSKO1d0/6u4pd19K8H/0Z+7+h1TgsZrZbDObk34NXAJspUS/4arvSW1mlxGUb04FbnT360obUWGZ2S3ABQTDBfcAfw/8CPgecCLB8Ohvd/fRFdllx8xeC9wP/JoXyqr/hqAeoqKO18xOJ6isnEpwofc9d/+EmZ1McJU9H/gV8G53P1y6SAsrLGL6K3e/vBKPNTymtvDtNOBmd7/OzBZQgt9w1ScIERHJrtqLmEREZAxKECIikpUShIiIZKUEISIiWSlBiIhIVkoQIiKSlRKElA0zGxh/rZzb/yBsZ54eUvn+UZ8/nDksepbt15jZV7Msn25m92X06s0Vg5nZz8ysLo/4f2FmLfnGn8f+PmdmFxbq+6T8KEFIVTCzlwNT3f3xjMVzzOyE8POX5fvd4UjA9wDviLD6ZcAWd+/Pd38ZChJ/Dl+hMoZGlzwpQUjZCa/CP2tmW8OJVd4RLp9iZv9iZr8JJ1W53czSQzH8IccOT/A9XjipXwncEmH3J4RX8p1m9vcZy38U7mM8I3GEwyr8xIJJf7ZmHMfrw4lxfm3BhE8zxviuMeM3s6Vmdr+ZbQ4frwmXLwzvdh4O93leOCrsTRl/z78AcPcngAVmdnyE45IKpAQh5egKgklylgMXAZ8NBzC7AlhKMPnTewgGr0s7F9g06ntuDbcBeBPw4wj7PhP4A+B04G3pIh+C8XJaI2yfGcdK4Cl3X+7uy4A7zWwmwSRP73D3VxAMt/DBMb4rV/x7gIvd/VUESeTL4fJ3AXe5+xkEf7+HCf6Wi919WbjPb2Z8z+YwZqlCShBSjl4L3OLuz7t7D3Avwcn5tcD33X3Y3XcTDOaWthDoHfU9TwP7zeydBAPdDUbY993u/rS7DwE/DPeZHnr7SHqgtRzmh5MZQTBm1MVm9hkzO8/dDwCnAjvd/bFwnW8B54/xXbniPw643sx+DXyfIGlCMEDl+8zsY8ArwlgeB042s6+Y2Uogs/hrD7BonGOSCqUEIdViCJiZZfl/EEw7G6V4CYLJXMZ6PwM4NM72z5nZFIAwCbyKIFH8g5n9XcQYMo0V/18QDM64HGgBpof7vI8g4XQDN5nZe919f7jeL4APEM7aFppJ8LeTKqQEIeXofuAdYdl5I8EJ7yHgAeAPwrqIJoJRbNO2AS/O8l1twD8Bd0Xc98UWTCA/i2DaxwcAwtE297r7s+Ns/1sg3ZJqETDo7v8GfJYgWfwWWGpm6VjfQ3CHNJax4p8L7HL34fA7pob7XAL0uPv1BIngVWbWAExx91uBvw3jSHsJQfGZVCHNKCflqI2gfmELwRX8/3b33WZ2K/B64FGCqWQ3AwfCbX5CkDDWZ35RWMTyGQDLmHkuh4cIyv5TwL+5+8Zw+evCfYwnHcd24BUE9SfDwLPAB939kJm9D/h+2Gy2A/j6WF+WI/5/AW41s/cCdwLPhMsvAD5iZs8CA8B7Cabq/Gb6zgb4aPh9xxEk1Y1IVdJw31JRzKzW3QfCK/qHgHPD5DGLoE7i3LC+oND7/SFwbUbdwVjrLQS+7e4XFzqGQjOz1cCr3P3/ljoWKQ3dQUil+U8LZlqbDnwyrKzG3YfCZqmLgd8XcocWzEb4o/GSQxjHLjO73szqCtQXIk7TgM+XOggpHd1BiIxiZpcSFttk2Onuq0sRj0ipKEGIiEhWasUkIiJZKUGIiEhWShAiIpKVEoSIiGT1/wFfQgpyN7SmBAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot_btfr_data(btfr_data, metadata) # blue = lenticular galaxy"
   ]
  }
 ],
 "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.8.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}