{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# StarKit Python3 Test Notebook" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "The first step is to load a spectral grid to work with. Using the open access grid from StarKit Collaboration (https://zenodo.org/record/2557923#.XFoc0KfMzYo)" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "from starkit.gridkit import load_grid\n", "import pylab as plt\n", "import matplotlib" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "[\u001b[1mstarkit.gridkit.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Reading index (\u001b[1mbase.py\u001b[0m:266)\n", "[\u001b[1mstarkit.gridkit.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Discovered columns teff, logg, mh (\u001b[1mbase.py\u001b[0m:269)\n", "[\u001b[1mstarkit.gridkit.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Reading Fluxes (\u001b[1mbase.py\u001b[0m:272)\n", "[\u001b[1mstarkit.gridkit.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Fluxes shape (688, 13184) (\u001b[1mbase.py\u001b[0m:275)\n", "[\u001b[1mpy.warnings \u001b[0m][\u001b[1;33mWARNING\u001b[0m] /Users/devinchu/anaconda2/envs/starkit_python3_test/lib/python3.7/site-packages/astropy/units/format/generic.py:474: UnitsWarning: 'erg/s/angstrom' contains multiple slashes, which is discouraged by the FITS standard\n", " core.UnitsWarning)\n", " (\u001b[1mwarnings.py\u001b[0m:110)\n", "[\u001b[1mstarkit.gridkit.base\u001b[0m][\u001b[1;37mINFO\u001b[0m ] Initializing spec grid (\u001b[1mbase.py\u001b[0m:309)\n", "[\u001b[1mstarkit.gridkit.base\u001b[0m][\u001b[1;33mWARNING\u001b[0m] **** NO WAVELENGTH TYPE SET DEFAULTING TO GRID (vacuum) ****\n", "\n", " (\u001b[1mbase.py\u001b[0m:75)\n" ] }, { "data": { "text/plain": [ "" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "grid = load_grid('phoenix_t4000_10000_w3000_9000_r3000.h5')\n", "grid" ] }, { "cell_type": "code", "execution_count": 9, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Flux')" ] }, "execution_count": 9, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAERCAYAAACU1LsdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXeYFeX1+D9nO2WXtoBIWzpSBVaaotgRE40Go8QeS9TYY/yhJsYaazQx0Shi9KtRjDUWUFFRRASkdxEEhKUjVYqw8P7+mLl35947t+3euWXv+TzPPjvzzjszZ3bvfc+cc973HDHGoCiKoigAOakWQFEURUkfVCkoiqIoflQpKIqiKH5UKSiKoih+VCkoiqIoflQpKIqiKH4yUimIyL9FZJOILIyh7+MiMtf++VZEtidDRkVRlExEMnGdgogcC/wIvGiM6RHHedcBfYwxv/FMOEVRlAwmIy0FY8wXwFZnm4h0EJEPRWSWiEwWka4up44ExiZFSEVRlAwkL9UCJJDRwFXGmGUiMgB4CjjBd1BE2gLtgIkpkk9RFCXtqRVKQUTqA4OB10XE11wY1O084A1jzMFkyqYoipJJ1AqlgOUG226MOTJCn/OA3yVJHkVRlIwkI2MKwRhjdgIrReQcALHo7TsuIl2ARsDUFImoKIqSEWSkUhCRsVgDfBcRqRCRy4DzgctEZB6wCDjTccpI4FWTiVOtFEVRkkhGTklVFEVRvCEjLQVFURTFGzIu0FxaWmrKyspSLYaiKEpGMWvWrC3GmKbR+mWcUigrK2PmzJmpFkNRFCWjEJHvY+mn7iNFURTFjyoFRVEUxY8qBUVRFMWPKgVFURTFjyoFRVEUxY8qBUVRFMWPKgVFURTFjyoFJSzvz1/H9j37Uy2GoihJxDOlEK2Osp3J9AkRWS4i80Wkr1eyKPGzfsdern1lDlf/Z3aqRVEUJYl4aSm8AAyLcPw0oJP9cyXwLw9lUeKk8qCVKHHW99tSLImiKMnEM6XgVkc5iDOBF43FNKChiLTwSh4lPnJzrAp2+w8eCtvn4CHNsKsotY1UxhRaAmsc+xV2WwgicqWIzBSRmZs3b06KcNlO8IA/c9XWgLbHPv6WDreP58OFG5ItmqIoHpJKpSAuba6vnsaY0caYcmNMedOmUZP8KQnAqQC+XrmVEU9P5V+fLwdg5ZbdPPHpMgA+WLg+4Lwdew5w1lNT+P6H3ckTVlGUhJFKpVABtHbstwLWpUgWJYjKQ1Vuo/U79gKwdOOPABxwuJQqHcrjo0UbOO/ZacxZvZ1/TFyeJEkVRUkkqUyd/S5wrYi8CgwAdhhj1kc5R0kSzsH+kF2dL9e27XKkysibu3q7bRnsYevuqumr2/ccSI6giqIkFC+npIbUURaRq0TkKrvLeGAFsBx4FrjGK1mU+PHNPnLyv7mWIecLQgOs3b6XOau3BygEgE+WbGTxup0BbYc0MK0oaY9nloIxZmSU4wb4nVf3V6rPxp37eHbyCtdjCyp2sGHnvpiu882GnXQ7vASATbv20f/+T3ng7J6M7N8mYbIqipJYMq7ymuI9V/9nFrNXb3c99vN/fhnzdZwuqDVb9wBw93uLqNi2h1tO6YKI21wDRVFSiaa5UEIIjgfc9N951brOuPnreWX6agC/Ath34BBPfvYd7W4bzx//t6BmgiqKknBUKSieMenbzdz+tjXw57hYBf+ZtjrZIimKEgVVCkoIXoSDc2LwFL03bx279umsJUVJJaoUlBBWbkn8wjNxXasI89ZYsYuZq7Zy3dg5/P616rmqFEVJDKoUspBXpq/m3GemJvWe4WLKb82uAGDE05Y8ExZvTJZIiqK4oEqhFrDvwEFufWNeyFqBcNz+9gKmrwzNVfjajDWUjRqXaPGA8ErhkLEC0k6M0fUMipIqVCnUAt6es5bXZlbw8Iff1Og6z325MkESBXLwkCHcOJ+bIyxctyOgrd1t45m8TBMfKkoqUKVQC/AFcYPfuOPFq2UD3e780J8qI5i8HKHSJT33B5p9VVFSgiqFWoBvuueunyrjOi/YTeM2bTQR/FR5iHAZLnJzhQMuKTVemb6a68bO8UQeRVHCo0qhFlDdwTy4ZoKXnvwde92nmuaKBGRddfLePE2aqyjJRpVCLSA3lkUALjiVwI69B1iyfmfYvjXl4n9/7dr+1OffhVUKkdi7/yD3vr/Yv65h8bqdLN/0Y41kVBRFlUJGM2X5Fj6OMIXz0CHDXe8uClvwxunnv+HVxLhqnru43L99+/CuMZ3z2syKmK+/Y+8Bduw5wOuz1vDclysZ/YWVuG/4E5M56bFJ8QmrKEoIqhQymPPHTOeKF2cGvGn/ZfwSTn9iMgBLNuzkha9Wcc3Ls13Pd4YUPl8a32yfX5W3cm0/umOpf7tPm0ZxXTMW+t37Mb3vmcD+SuuZ/zFxOZ99synh91GUbEWVQi2gIK/q3zj6ixUsCqpjsGjdTlcXzQ+791M2ahxffBv/9M8jWpS4tuc5XFl51XRrOfnPtO+ZvXqbf9+XedWZYfUhx1RcX5W4bzbsjGhFKYrijiqFWkDzkqKofaat+CGkbb6dYiJc7YRI9A1jBeTmCJP+MJR3fnc0+bk1/3j98X8LOfupr9iwY1/AwrpvHPEPZ4ruQQ9MBGDY3yZzxYsza3x/Rck2VClkKM6ZQ763Yydlo8bx476qKaoXPhca6D1oQt+6Y6G8bSN6t27oekxEaNukHr1bNyQvN3FTXAc+8GnA/uuzquIQwbOonPzomKZrjGGnJtxTlIioUshQ/jPte/92uHoHHy6KvACsutUxi/JzY+qXl5Ocj1dwAj9fkj2Aa16e7VcaT33+Hb3umlCt2U6Kki2oUshAZq/exp/fXRS13979ByMe9y1ei9f1H6thkZ9ASyEenDGVL77dTIfbxzNn9TYe+WgpANv2xJYjSlGyEVUKGcjU70LjA268OmNNxONz7JKbvpk8wQzt0jQuuR4Z0Stgv7rrJ2pK5aHQ5znrqa/82zVNB6IotRlVChlIdbNRBC9Oe+GrVQB8FUbJNCsu9G8/eHZPx/3dBTinvHXAvldpM6Lx7cZdEY/f/d5iAMZ+vZpNu/ZFjEkoSrahSiEDqe5ge9rfJ7u2t2jgPnvJWRgnx/HW73b3Pm1CA8+pshSilfls26QuFdv2cNtbC+h//6fc+N+5SZJMUdIfVQoZSKLH2vU79rm2O3XPz3q1oE3juiHtAG9dM5gXLu0fcn6qLIVo9GzZgK8d9SScOZb2HTjI7W8vYN320BldipINqFLIQFIx2Bbl5fKnn3UDrCR2Tvq2aUSDOvkh56TKUojG+/PXc3NQ2c/b317A7W8v4KNFG3hl+moGPziRzbt+ck3rrSi1mbxUC6BkBiJVs5ViXdeQSRXUXpluuZya1q+Koxx1/yec3bclj/3qyFSJpShJRy2FDGPH3gPcN25JUu4VEEcQ8a9riNUAyMT47d8/XRaw/9bstUxf8QPvzlunAWklK1ClkGGMqUZKiuoSPPY3K7HeoruGyXsUTGn9Ajo2q59gqZLPuaOncf3YOfS+e4I/Vbei1FZUKWQYG8IEhb0gOHbRt00j3rx6EDec2Cmm80WEO4Yf4YVoKeHHnyr5x8TlqRZDUTxFlUKGUZlEF0aPlqEWQb+2jeMKIKfpBKRqM/qLFfz62WmpFkNRPEOVQoZR0zH2hx9/irnvr+zFaE3qFbgev/6EjlEtgXSdlloTvvruB/43Zy1fLtuSalEUJeHo7KMso999n8TcV0SYetsJ1C1w/5jcfEqXqNeojUoB8C94W/Xg6SmWRFESi1oKiiu+UpotGtRxXYMQK+mmEzo3r17g+xdHHu7aXjZqHGWjxvGwo9CPomQyqhQUV7of3iAh1/FCKTSqG7+SumBgGwAa1g11hRXYxYBaNqwT9n4Pj+jNSUc0C3v9pz7/Lm6ZFCUdUaWQYSRrlXBw8rzqUlP30Yh+obWg86pR0W1Y9xbWuTkScs1D9iK7O3/ezfXcNk3qUZCXw+ZdkeMxW3fv5+Ahw7vz1rF1t6bnVjITT5WCiAwTkaUislxERrkcbyMin4nIHBGZLyLDvZSnNnB81/Bvq7ESyzi9KcoAGCs1VQpuJT2rU/t5QPvGnNOvFQ/9shdn920ZcMw3n6u4KHKILdpK7r73fszbc9Zy/dg5HP/o53y4UFN0K5mHZ0pBRHKBJ4HTgG7ASBEJfhX7I/CaMaYPcB7wlFfypDvLN+3i+x92R+2XLEvBl/yupoQTN9YArZtR4FbmM1pBn/zcHB45pzetG9cNUTS+dBwlRe5uqSEdSwP6ReLxj78FrJXnV/1nNjv26mI3JbPw0lLoDyw3xqwwxuwHXgXODOpjAN9k+AbAOrKUkx77guMe+Txqv2SlEzq1+2EJuY7z7fqXfUNdQVHPd5mE6yvz+bvjO3DhwLYA/LyXeyD4l31bMfH3xwWdH3hN39KP+oWhlsJntwzl5pM7B/SLxNqg7Krb1I2kZBheKoWWgLP0V4Xd5uQu4AIRqQDGA9e5XUhErhSRmSIyc/PmzV7ImkEkRytUx0XjhvMytw6LPoX1/37Tn/evO4bfn9w5xJ3jq/vgky0W11T7pvVo3zRwxlG42tH1XdxH7Urr+XNAhct91NelloSPDTuTtwJdURKBl0rB7Rsb/K0aCbxgjGkFDAdeEpEQmYwxo40x5caY8qZN4ysRWdtIxILmWIb73ATVV3ZaCs4MpOE4rnNTerRswHUndmLBXacGHPMpAZ8LTUT88ZFwf5ZoKb3/c9kAv/JxWgof3jiEN64aFHDeoTBmWrvS8NNczxs9jXfmrg17XFHSDS+VQgXgrM/YilD30GXAawDGmKlAEVDqoUwZT7iBKdEE10yoLk5LwXfJc4PKdsZLrt9SqGoL9vePvrAfD/2yJyP7twmVyf7Ud25en2M6lfLSZQO4ZmgHivJz/X26HlZCeVnjgPMePac3xwfVre7YrD73nNk9orw3vDrXP5tr3prtrN+hBXyU9MVLpTAD6CQi7USkACuQ/G5Qn9XAiQAicgSWUsh2/1BEEmEpxHKNRAW0nS4eEWHR3afyF0e95+rgtxRcIw5Vfc49qk3E5/CdfWTrhtw6rGvU+/Zo2YDnHRXm3rpmMJ/cfBz1CvP4c5jprD5WbLYmEZz55BQGPTCRh3Sxm5KmeKYUjDGVwLXAR8ASrFlGi0TkHhE5w+72e+AKEZkHjAUuMZlUmSUFHPIwIV7Xw4r924mKKQQbHPUK80IG6q/vODGua/lcUgGWQpT7OknUJ8z5N7r06Ha8efXgsH2XbtjJxG82+vf/pYvdlDTF09xHxpjxWAFkZ9udju3FwNFeylDb8NJ91KR+1WpfLyyFYObdeQqSE34qKAQO7r5tn2g5ORJ27UB4G6JKKYQTrXGYBIDBBP+N+rVtxPOXHkWX5sUUF+Wxa18lgx+cCMATLim391ceoiBP148q6YUmxMswvMyc7RxIYy25GY1ISqFBnOkqfPL54h3OSxcEL2iIIL5vnUPzkqKQY+OvH+IvJhQNN8V5fJeqxYXRKrVd8/Isxlx8VEz3UpRkoa8pGYaX7iMv8GKtXY7LlNQjWpRw/1k96N0qes6mzs2L+ctZPfnbuaG1l7sdXkJpDLOkILqLLVo6jk+WbGL3T5UAzFi1lacnqUtJST1qKWQYyZp9lChiNTgm3HRs1DQTZ/dtyd8+WcZh9ht+8Jh8/oC2fLRoo8uZofx6QOispHiJtk4ilrhM9z9/FHLOOeWta5SZVlFqgloKGcaotxZ4dm0vMprG6obq3LyYFg3cs5T6uHpoB765d5g/9hFpUE5GMpBwi+CqjscvxX3jltD77gl8syExCQkVJV5UKSieksgiO4JQlJ/rek2f/ZTMyWtRdEKNgvXD/jaZslHjqDx4qNrXUJTqoEpB8ZRExhRygmYfHTImrHWTqEB5JKJZComQoeMdHzB+wXpd8KYkDVUKSgBDOiV2QXlNLYWZq7b5t6vWJ1i/nUZBKpa3xGIJ/PH0IzjHpSZEPFzz8mwGPTARY4xaDornqFKopbRtEn/q67vO6M6/LzmKhXefGr1zjNT0ZXmxo9iPBG0cMpHXI3hNLErh8iHtOaJFSUDb7D+dzLBqZKF9dMJSOt7xAXv2V8Z9rqLEiiqFWsrpPVvE3b9D0/rk5+a4ppCuLol041QtXrM23GZi+dYGJCp3UyRijRkEy9m4XgGPu0yHjcaTn1lTVk94dBIzV22N+3xFiQVVCrWUEyPUE3bFozE0kTGF4PQWhw6FxhR8SiFaEDiZuK1arlNQlXyvcb2CuCy7DTv3MeLpqepKUjwhjb46ihv7Kw9x03/nxlSVzUl+bk7UamTJIJGzj3zk+i2Fqjbfy/hVx3UAoFuQy8YLimO0qM47qg3XndAx7PHLjmnHk7/uC1huvwk3HcvntwyNet23ZmtKbiXxqFJIc2as2srbc9Yy6s341yfE47rxSn14ufbhkDEheYyO79qMVQ+eTsO6seUvqgk5MZpBBXk5/uptbpTWL/Cn3ijIzaFz82LKSutFve5d7y0C4MtlWzjl8Um8NmNNlDMUJTq6ojkNmL16W9hjvmHHxFhxrW5BLnv2HwQgP0dIdTFILyyFqtlHmbO6201BL71vGK/PrGBEv9YsWLsDgML82N/T9uw/yNw127nguekA3PrmfGas2srVQzuEVJtTlFhRSyFF7Nh7gOMf/ZzPl27i7Ke+8rf/dcJS9lc6fMW+ymIGfqo8GPW6zkE4Wu6dZOCNUrB+H8wgpeDDOeW3MC+XCwa2JTdH6NK8mJ4tG3DvmT38x8/ofXhYN9iIfq1oUq+AXzw5JaD99VkVnPDXSVz7yuyMy5OlpAepHzWylKnf/cDKLbv5y/glAe3/mLic12dVuQF8Uy4N0OWPH/rbv77dvQZBoR3UzBGJK82CV4u9vEiIl6nM+uNJjLm43PVYnYJc3rvuGPq0aeRve2JkH8bfMIRx1x8T0r+4KI+d+w6Evdf789cz5bstNRdayTpUKaQYt5fdfQeqLIUcv6UQ1DHMYPuPX/fhmqEd6H54id9PnUq8XEcgCOcPbEPDuvkMj3MKbipoUr+Qwrzc6B2D6H54aObXeWu2c+BgZEvg3bnr+HbjLjbu3Bf3PZXsRWMKaYhTAfiCmdFy8/to2bCOv7RktDQMScFjvdShaX3m3nmKtzdJA84tb82Xy7ewdvteRCz3YzRen1XB67MqAPj7eUdy5pEtvRZTqQWkwaihBONUAD5LoTJIKYR7A3f68OPRCV6N3V64j7LRU/7QiF78+xKrIE+uCH8/r09c59/w6lymLFd3khIdVQppSGAJSmvH6SroX9Y4pqmescRhzzuqdbzixYWXiemSsGjZlTEXlfPCpcmvmOZzB+aIhJQMvfGkTlHPP3/MdKZ+9wMzdDW0EgF1H6UhTisgx68UAmckhRsPnQPlb49tz5/eWRTbPT0aYGt62ZO7NefjxbEVzkkWJ3VrnpL75tumX05O4CrpO3/WjYsHl7Fyy24uHlxG60Z1eX7KSp76PLSS28hnpwFwbOemjLmoXGtEKyHoJyLFxOoKcaY0EMK/gTvbLxxUxj9/bbkZwiVg8/ptu6ZTUof3jD9xXG0l12Ep5DumGzeok09ujuVS6tumEU2LC/2WRO9WDXj6gr4h1/ri2810/uMHlI0axztzdWW0UoUqhQTz9cqtrNj8Y0DbJ4s3smFH8AyQ2NSBL5maM6YgEr4sZ7ghOHx8wVutkCoXT20k3w7Q5IpQ4FAKbon5fJ+P/u0aM6xHC777y3C+vuNE19ToN7w6l7lrtnsktZJpqFJIML96Zion/HVSQNvlL85k4AOfcsKjn/vbfGN6tECsbybSgQBLQQLdSQ7ifTP3dU/XsTsD16d5hm/wz82VgLxWbuk2fB8P37HcHKFZcREvXTbA9dq/eHIK142dowveFFUKyWTFlqqkdle/PDtsP2dKC9931BloFoG6+e7hoLCVyMIM+14Pul6saM5WReH7LOTlSIB14LZI0WcpuKUQb1LPPS/Ue/PW0f728cxctTWjUogoiUWVgsdE+3JFW9zlO/1ApXNBm9Cgbj69W4Uuagp3tWi5k7yaJVTTy3ZuXhzSFmseqNpGaf0CLhlcxkuXDQj4f7kpXn9dCReF4YxHfHzTsSHHRzw9lee+XJkIkZUMRJWCx0SzxqMNmuFiCgBNiwtdLhiPdN5TU3F6tAxVfFWutzR7WI8REe46o3tIJTe3gd+nFNyU/bMXWak2/nj6EXRqXsxhJUUhff464VsAXpq6ik27dEV0NqFKwWOCLYVZ34fPiOrGeaOtKYRuyd/cjJBwlkeqylZ6MXD7FGW2KYVwxOs+6tmqAasePJ3Lh7QH4PWrBvHRjcdy+THt/H32HjjI23Mq+NM7izjxr5N4bMJSdSllCaoUPCbYUvjlv75y7xiFgCmpEQbDcIHrqO6jakkVnUSM25NvPZ43rx7k3y9rYtUaaNc0es2BbMAt0OxTCrGk4m7duC5dDivmjz/rFtB+03/nAbBrXyVPTFzOqzPWhJ3goNQedPGaxyTK/+1ULr4hwG3ADVYYqSxsD4mJVbRuXJfWjavKVZ7dtyVlpXXp68goms24WQNXHtuBrbv3c9Ggtgm7z21vLaBeYR6ndGtOUX78if2UzECVgsdEs7irM2j6TnF3HwXdP2qAOe7bpxwRoV/bxqkWI21wW4PSoE4+D5zdK+H3un7sHAD+MbIP9QvzOL5rnLXAlbRH3Uce47bIzLmCtDpjcqRzMnGQV2pGIrPh3n1Gdw4rKeLnvQ8HrDxPblw3dg6XvjAjYfdV0ge1FDzG7W3+hlfn1uiavgBrItxHfvlUmWQsiSywd/HgMi4eXAbAg2f3pF5hHh/eOIRhf5vs2n/Ljz+xc+8Bdu6r5MjWDRMniJIyVCl4TLh0FDUhovsoaHCPNaaR6tiDUn28moVVr9AaHroeVsJLl/Xn758sY2bQ7Lny+z7xb4+5qJxerRvQrDh0iquSOXjqPhKRYSKyVESWi8ioMH1+JSKLRWSRiLzipTxec+rjX4S0RRuSq/d9Dn9SvJdTd1Pmk4xiSkM6NeWNqwdH7HP5izPpf/+nnsuieItnnyYRyQWeBE4DugEjRaRbUJ9OwG3A0caY7sCNXslTEz77ZhPnj5kWdZ720o27/NufLN7I1t37MVFm8FVnUPbnK0qk+0jJWJJZYC9cbXAnM1Zt5b73FwdMo1YyBy8/Tv2B5caYFcaY/cCrwJlBfa4AnjTGbAMwxmzyUJ5q89uXZjFl+Q/8VOn+Id++Zz+3vbUgoO3yF2dy6QszYlgfEH7QXlCxI8w5Fr4B/aQjmoUc8xGz+0gthozFbUWzVzQrKeKUKPUkznl6KmO+XEnHOz7g5enfJ0kyJVF4qRRaAmsc+xV2m5POQGcRmSIi00RkmIfyVJ8o37nHPv6WsV+vDmlftWV31DQXC9a6D/wA2/fudxcnSJ7zjmrj347fv2ynQ4jzrHjo0ryYe8/s7uEdshu3Fc1e8peze4a0TfrD0IB03j7ueHshZaPGsWT9zmSIpiQAL5WC2yc1eIjMAzoBQ4GRwBgRCZnCICJXishMEZm5efPmhAm4bOMu9od5+3fi6/Pm7ArX4+GCyTv2HmDM5BVxy/WfaaspGzWObXvci7NHsi6CdUKsAWQvLYWPbjqWCweVeXeDLCfZ6T5K6xcy3eFG6tSsPm2b1GP2nSeHPef5KZpgL1PwcvZRBeAsANwKWOfSZ5ox5gCwUkSWYimJgAnQxpjRwGiA8vLyhHjBN+7cx8mPf8HI/m14wOXNx42p3/3A+QNCV4hGGnjdSiJGY/XWPQBs2LHX9XikMSD+2UfqN8p0kuk+8tG8pIjHz+1NSVE+A9s3AaBeQfhVzq/NrOD9+evZs/8gvVo14NZTu3KMS8EfJfXEZCkEB4jttqFRTpsBdBKRdiJSAJwHvBvU53/A8fb1SrHcSfG/WleDHXutt/B0LmL+1mz3MomR3gx1amn2karEgGf1acWJRzT3T111TnJo07guI/u3ZkC7qpXne/YfBGB+xQ4ueG46i9epSykdidVSeE1EXgIeBors3+XAoHAnGGMqReRa4CMgF/i3MWaRiNwDzDTGvGsfO0VEFgMHgT8YY36o/uMklr37D/Lxkqqi8eHeub36Tn6zYZf7gQj3C35pjNl9pMokY8nLTb//3Re3Hu/fLhs1zrXP8CesBXG/7NuKW4d1oblLCm8l+cSqFAYADwFfAcXAy8DR0U4yxowHxge13enYNsDN9k9KiPR1um/cYl6eXhVADhd/SPZXMnKai9RUWFNSh1tCvEzizdkVTFvxA1NGnZBqURRiVwoHgL1AHSxLYaUx0WbgZz4bdwYWF9l34KBrP6+qloXDd7+2dgrphnXzq44F9fXFFMKLqNoi03FLnZ0qZv7xpGrVeV67fS9lo8bRrUUJY68YyE+VB2mmlkNKiHX20QwspXAUcAzWQrQ3PJMqTQgO4LVxpG/2UbFtD3NWx1c4p6b4pPp/p3Xh35eUU15W5beNt0Zzu1JLsfRwKe2pZAbJnpIaidL6hVEH84d+2ZMTw2RXXbx+J73vmUD/v3zKLa/P80JEJQqxWgqXGWNm2tsbgDNF5EKPZEobgtMHlJeF5u8/5qHPkiWOH9/AX5iXywldmwcdi2+AGNS+lI9vOpaOzeonSjwlyaSTpeDGS5f1J0eEeRXbOemI5nRuXkxxUT6ffmOtVX3sV725+bVQBfDGrAoWrdvJvWd2D3jxUbwlVkthk4i0cf4Ak7wULB0IthTSxS/vNgQ8fm5vuh4WWuQ+6rUEOjUvTroLTEkc6R5TGNKpKUd3LOWaoR3p3Nz6jJ7W4zCevaiciwa15aw+LVn5wHDXc5es38mIp6fyz4nL2Lv/oJYETQKxWgrjsJzPghVTaAcsBWr1MtW0VQoug8BZfVpxVp9WcV9L6xxnPqlYp1BTRISTuzXnZEfKjKcv6Mv8ih2ua3senfAtj074luO7NOX8AW05KUqqDaX6xGQpGGN6GmN62b87YeU1+tJb0ZJDpHE+eMD8dtMuNu3aF6Z38kjkON5e6xxnPJmoFNwY1qMFI/u3idjns6W3S08WAAAdV0lEQVSbufzFmYyZvKJaAW0lOtVKc2GMmY0VdK7VBAfwnpm0Ii1SA8ezpiCSdTOs+2Faa7cWkO7uo3goLrKcF6X1CyL2u2/cEtrfPp5XpofmHFNqRqwrmm92/Nxi1z1IXBKiFDDiX18BVf75NVv3sGFHoBUQLoD3hxTPiqjWGCCw/P7TuO6Ejv6mTs01uFwbSPdAczw0rFvAmIvKeeWKgf62/145MGz/299ewE3/ncud7yzkoFoOCSFWS6HY8VOIFWMIToOdUezcVwlUuY+GPPwZAx8ItALClTl8fZZ7YrxkUd0hIC83x+9qaN24Djec2ClxQilKgjipW3NaN6qa/j2gfRPm3XlK2P5vz1nLi1O/p8Pt49m22z2zsBI7MQWajTF3ey1IujB52WbqFebRt02jiBWttu3eT597P06iZFXEYynk25qtKM9yE/lqNGza+RN5iSzuqygJpCg/8LPZoG4+OQKHDEy97QTWbtvL3e8tDkk93+fejxl//RC6HV6STHFrFRGVgoi8R4RYrDHmjIRLlGSWb/qRpydVzXa48LmvAVj14OkRZ+ZUbHPPYJoM4okpnNytOdef0JHLjmkPwCI7CVm4gkGKkg6ICF/84fiAF6DTex3Oe/PW0ahuAS0a1OG24V359bPTyRGY+Puh3Pjfucxds53hT0zmZ71a8M2GXdx8cmeG92yRugfJQKJZCo8mRYoU8+AH37i2R3qRrjyUukE1HkshN0e4+ZQu1TpXUVJJmyaBGQQe/mUvbjixk39yRI+WDejQtB73ntmDstJ6jL6on38iyPvz1wNwzcuz6dGyhOtO6MSp3Q9L7gNkKNGUwkpjTNaG9/cdCD/wpzKoVZOFZuEKAilKulOnIDdg5X1JUT6f/n6of79p/cKA/v3aNmLW99tYuHYnv31pFr89rj1DOjbVOg5RiOZU/p9vQ0Te9FiWtOPwhnXCHkutUqj+uVpLXamtiAhjHbOWHhnRK+D4M5NWcMFz09O6hko6EM1ScA4/7b0UJB0pzAuvM88dPS2JkgRSEw+QWgpKbaZv26pqviV18l37nPP0VAAa1yvghUuPolerkArAWU00S8GE2c4KpizfkmoRXKmJpVCppoJSiynMy+Xa4621OCVF+bx59WCeu7jcte/W3fs5459TOPPJKeza514PPRuJphR6i8hOEdkF9LK3d4rILhGp9bX0fFkc042aVEnTBT61h4m/P45XIyzsylauP7ETX99xIgV5OfRr24gTj2jOlceGOjpeuqw/APPWbKfnXRPYs78y2aKmJRHdR8YYzYGQhtTEUjhwUJVCbaF90/q0b6qr0oMpyMuhWXFgTYfLh7Rjw4593H9WD56dvJKfKg8yuEMpfdo0ZM7q7QCc8OgkNuzcx5BOpTz2qyNpWlzodvlaj65eykBqElP4y9k9EyaHomQKzYqLeGJkH4qL8rn55M7cdtoR5OYIb19zNAPaWbUaNtiVFicv28JR93/CJc9/nZVuJVUKYVi1ZXeqRQhLTaakDrGn4+WnYbF3RUkF9QvdHSafL91Mz7sm8MOPPyVZotSiSiEMl74wI9UihKUm7iPfuToJSVEs7jojtCzMa78d5N/ud98nvDh1Fc9PWZkV6bpjLbKTdew7cDDVIoSlJoHmmpyrKLWR1o7a6yP6teLPP+9GcVE+vz2uPc9MWgHAne8sAuDu9xbTtLiQ5y85ih4ta2ddc1UKYVi/I/XFdMJRE0vBl2W59r/vKErszL3zZESEBo61DVcMqVIKTjbv+omf/eNLrj+xE02LC7lwYNtkiuo5Wec++seny8LmOsoUavKur7WYFSWUhnULAhQCQGn9Qp6/tKqWWPvSwCqFT3y6jD/9byFlo8bx/vx1SZEzGWSdUvjrx98GZEXNRBJiKWhQQVGiMqRjKcN7HsarVw7k/37TP2y/a1+Zw8mPTWLzrswPSqv7KAOpydu+71y1GBQlOnm5OTx1fj+gKs54Ru/D6dWqAfeNWxLQd9mmHxkzeQVrt+9l7prtvPbbQRHzp6UrqhQykJq5j2p+DUXJRoryc/n8lqG0bFSH/NwcLjumHe1uGx/Q55kvqmIQgx+cyJBOpfRp3TAgfX26k3Xuo9pATd7yc/yWQqKkUZTsoay0nr+aoYjw5K/7hsQanExetoUnJi73VzzMhCmtqhQykBqtU/D/Vq2gKDXl9F4tKAtSCo3q5vP/hnUNaLv8xRkcOHiII++ZQNmocbw4dRWQnkpClUIGcGr35gH7NRnOc9R/pCgJJTjJZKdmxVw8uC13n9Gd2X86GYCNO3+i0x0fsHOflXTvzncWsXLLbtrfPp7y+z7m86Wb0mbyhyqFDKCkKHCqXCJWNKtOUJTE4KtR4qu/kpcr1C3I4+LBZTSuV8BFg9zXMRz/6OcAbPlxP5c8P4N2t41Pi9T2qhRi4D+XDUjp/fOC8hTVxPXjexnRmIKiJIZKO/Owr6hPbk7gl6t1o7oh54Sj4x0f8PSk7ygbNY7XZqxhw459SbcgVCnEQKpruuYEjeA1GdAL7LeZYJ+noijV46A9aPsWv/kC0T6K8gP3n7mwX8Tr+RbX3vrmfAY+8CntbhvPuu17EyVuVLJKKXy8eGOqRagWeTnBlkL1yc0RVj14Opce3a5mQimKAlTFFEqK3Gf45wUpiVO7H8bgDk0AWHLPMNdzghXJ4AcnMnfN9pqKGhOeKgURGSYiS0VkuYiMitBvhIgYEXGvm5cgrnhxppeX94ycYKWgvh9FSRt8MYV6dgru4G9nrv19rVeQy4w7TgLgX+f3Y8JNx1KnoKqO2VvXDPZvf3PvaSH3+cWTU/jqO+9LBHu2eE1EcoEngZOBCmCGiLxrjFkc1K8YuB6Y7pUsmU5uAt1HiqIkFt+00noF7sOpL8bQrKTIX82tQd18GtQNnEByZKuGAftPjOzDm7Mq+Hnvw1m2cReTvt3sj194iZcrmvsDy40xKwBE5FXgTGBxUL97gYeBWzyUhW837vLy8p4SHLjSNQaKkj74Ygo+l0/wS5tvoki0b22wR+CM3odzRu/D/fu3DT+iZoLGiJfuo5bAGsd+hd3mR0T6AK2NMe9HupCIXCkiM0Vk5ubNm6slzNokBmqiMf76IXH1D1EKqhMUJW3wvb37JnEED//R1gZ1alafsibWDKUxF5Xzzu+O9kLMmPHSUnD7E/htHxHJAR4HLol2IWPMaGA0QHl5ebXsp3QaR9s2iX2KGrhZCoqipAu+mEKVUgjEN1Ek3Pf2oxuP9V/jpG7Nw/RKHl5aChVAa8d+K8CZdLwY6AF8LiKrgIHAu14Hm9OBeN/0EzklVVGUxHL/WT3p26YhnZoVA6HfT99LXbgJIjk5EjJDKZV4KckMoJOItBORAuA84F3fQWPMDmNMqTGmzBhTBkwDzjDGZOYUoTiINyYQ6j5SraAo6cJRZY1565qjqZNvzSQKmX0UxVJINzxTCsaYSuBa4CNgCfCaMWaRiNwjImd4dd9UM6Bd46h94h3TOzStX01pFEVJFlUWQWzt6YqnNosxZrwxprMxpoMx5n677U5jzLsufYems5Xw+Lm9aRQ0hcyNYzs3jdpHBC4ZXBbTfUf2b0N+boZ8mhQliwlOR+Nvz7FnJWWIrZA+jiyPqanL5aw+rTiqLLoVEJMsCHed0T2kfdWDp7v0NiExBUVR0g9/rZLg2Uf2KJspX+OsUQrpRDwfDmMy58OkKNlMuAJWfkshQ77IqhRSgNtHY+wVA8P2V0tBUTKHcIvX0qVeQjSyRilUd1h9YmQf//aQGOIFMcniMsiX1HFfMlKUnxsivOoIRUk/DO6Dfn5OZg2zmSVtCujXtpF/+4IBbRJyTbcxPVwQ6pZTu2RIeEpRsht/rZKgb2zwlPJ0J2uUQnXfrp3/z0T5BN0uE+7S9Qvz1H2kKBlAVbqGwPYMMxSyRylUFy8G5HiVi+oERclcfJZDhoQUVClEI1lv6ZFuo5aCoqQ/vkBy8LfV9/U9lCFaQZVCFBLtDvzFkYe7tkda2KIqQVHSH9+YH/wSV9cupFNWWi/ZIlULVQpRCA4S/a+GaW1H9Gvt2h7RGAiefaRqQlHSjkq72E7wyuZWjeoy5qJy/vqr3qkQK268TJ2dVlR3IE30gpNwl4t0F3UfKUr6061FCQAndg1Nf50OKbFjJWuUQnVJtPsorFKIcB9VCYqS/nQ7vITF95xK3TBlOTOFrHEfVfdlO9EFbga2axL3OcFl+hRFSU8yXSFAFimF6hLsuoll/kDzkkK6HlYc0i5SvQFeVYKiKMki89Wax8Trz88RYfrtJwHQ6Y7xHDhY82lomZJIS1GUzEcthShEerF3swbyHCfEMy05Ut9gnaA6QlEUr1ClAJwfIadRSH1kx3Zwoe6LB7XlwkFtw16rumO56gBFUZJF1iiFiMsAIs38CTrmfKEPfru/+8weVlZTm0fPiX1eciSjQqekKoqSLLJGKUQi4mriCANytGXrv+jTMu5rFeaF/ktC3EcR76ooilJ9VCnUAKdOqOnLvE8xjb6o3N9WXGjNA1BLQVGUZKFKgfDFMaJxWo/DEiwJFNmWQrPiQr6+46SEX19RFCUSqhSILwWGs+fvju/I5FuPD2kPR2n9grDHfIrJTp9CWWk96tiJtEKC3Wo4KIriEdmjFGIcSP99SXnE406bIidHKKmTb10+ykj98IhevHHV4Kj3b9mwDgBDu1SV/lQloChKstDFazXEN2BHG7d/Ve6eHdWHLz7Rpkldvr7jRJrWLwy5R3BfRVGURJM9lkIEnDGFWF1JvVs1sPvb51Xjbf6fv+7j2t6suCjA8tBAs6IoySJrlEKkwT7gzVugY7P6sV/XHrCrk5q7naPoRsQVzUH7aigoiuIVWaMUYiVHhA9uGBL2uAkavSVkI/Fo7iNFUZKFxhSCyBHIz83hllM6U+rw64fgsxBijCnUBNUJiqIki6xWCq0a1aFi297AGUX2CHztCZ3iulZ1Bm6n0VG/MPy/IsR9pP4jRVE8ImvcR26Ddic7dhCwMjnO6/rOrU5M4eChqhu3aVI3bD8NNCuKkiyyRim44RuSnWNuvP57t2s4iRS09hX6PrJ1w4j3CE3Kp6aCoijekNXuI7dx3K1GgpPg4bgg19KrZwUlv/Px0Y3HhgSnfdQrtFYst3fMQnIjpPqb6gRFUTwiq5VCMEe0KKFRvfCpKJz4humCvBzm33UK9cLUZrVqPLubEV0PK+HZi8o5pmNpNaRVFEVJPKoUEkBJUX61zz25W/OofapT11lRFKU6eBpTEJFhIrJURJaLyCiX4zeLyGIRmS8in4pI+LJlNZXFXT4g/d0xqhIURUkWnikFEckFngROA7oBI0WkW1C3OUC5MaYX8AbwsFfyJIpUKBCdfKQoSrLw0lLoDyw3xqwwxuwHXgXOdHYwxnxmjNlj704DWnkoT1hKiiwvWpfm8aS38EqaUEIDzWlu2iiKkrF4GVNoCaxx7FcAAyL0vwz4wO2AiFwJXAnQpk2bagnjNoz6BtfWjesy9oqBUaeGpgo1FBRFSRZeWgpuY5nrK66IXACUA4+4HTfGjDbGlBtjyps2berWpcYM6tDEX9QmEk3s2Um9WyVPgWjuI0VRkoWXlkIF4Cwi0ApYF9xJRE4C7gCOM8b85JUwifK4lJXW4/3rjqFz88jrGRKJ1lNQFCVZeGkpzAA6iUg7ESkAzgPedXYQkT7AM8AZxphNHsoScRVwvC/iPVo2oCAveYvB1U5QFCVZeDayGWMqgWuBj4AlwGvGmEUico+InGF3ewSoD7wuInNF5N0wl8tqQgLNKZJDUZTaj6eL14wx44HxQW13OrZP8vL+gcIk7U4JR0MKiqIki6xJiHfQxRFflG8FlvNz0vvPEBxo1gXOiqJ4Rdakuah0pKnu1aoB8yt28Oefd6ddaT3O6uuezC5dCLYUItVeUBRFqQlZM7ocPFilFG4ffgRdmhfTqF4Btw7rmkKpYsMXU8jPFf7fsK5cMNCzbCCKomQ5WaMUnJZCbo7EnA01ncjPzeHyIe1TLYaiKLWY9HamJxBnlbNMm+fvs3LyNJigKIrHZJGlcMi/nWm5g4qL8hjW/TAuPbos1aIoilLLyRqlULFtb6pFqDY5OcLTF/ZLtRiKomQBWeM+WlCxw7+dWXaCoihK8sgapfD0hf3oHEdqbEVRlGwka5QCQMO61oyjDAspKIqiJI2sUgq+uTuRkuMpiqJkM9mlFHRGp6IoSkSySikoiqIokckqpSA+B5J6jxRFUVzJLqWgOkFRFCUiWaUUfOjsI0VRFHeySin4lIEGnBVFUdzJKqXgq6usieUURVHcyZrcRwCPjOjFi1O/56iyxqkWRVEUJS3JKqXQrKSIW07tkmoxFEVR0pasch8piqIokVGloCiKovhRpaAoiqL4UaWgKIqi+FGloCiKovhRpaAoiqL4UaWgKIqi+FGloCiKovgRk2HZ4URkM/B9NU8vBbYkUJxUos+SftSW5wB9lnSlJs/S1hjTNFqnjFMKNUFEZhpjylMtRyLQZ0k/astzgD5LupKMZ1H3kaIoiuJHlYKiKIriJ9uUwuhUC5BA9FnSj9ryHKDPkq54/ixZFVNQFEVRIpNtloKiKIoSAVUKiqIoip+MVwoiUiQiX4vIPBFZJCJ32+3tRGS6iCwTkf+KSIHdXmjvL7ePlzmudZvdvlRETk3R8+SKyBwReT/Dn2OViCwQkbkiMtNuaywiH9vP8rGINLLbRUSesGWeLyJ9Hde52O6/TEQuTtGzNBSRN0TkGxFZIiKDMu1ZRKSL/b/w/ewUkRsz7TkcMtxkf98XishYexzI1O/KDfZzLBKRG+221P1fjDEZ/QMIUN/ezgemAwOB14Dz7Pangavt7WuAp+3t84D/2tvdgHlAIdAO+A7ITcHz3Ay8Arxv72fqc6wCSoPaHgZG2dujgIfs7eHAB/b/ciAw3W5vDKywfzeytxul4Fn+D7jc3i4AGmbqs9iy5AIbgLaZ+BxAS2AlUMfefw24JBO/K0APYCFQF6sS5idAp1T+X5L+gfT4D1wXmA0MwFr1l2e3DwI+src/AgbZ23l2PwFuA25zXMvfL4nytwI+BU4A3rflyrjnsO+7ilClsBRoYW+3AJba288AI4P7ASOBZxztAf2S9Bwl9gAkmf4sjnufAkzJ1OfAUgpr7AEwz/6unJqJ3xXgHGCMY/9PwK2p/L9kvPsI/C6XucAm4GMsjb/dGFNpd6nA+iBB1QcK+/gOoImz3eWcZPE3rA/EIXu/CZn5HAAGmCAis0TkSrutuTFmPYD9u5ndHk7mdHiW9sBm4HnbrTdGROqRmc/i4zxgrL2dcc9hjFkLPAqsBtZjffZnkZnflYXAsSLSRETqYlkCrUnh/6VWKAVjzEFjzJFYb9r9gSPcutm/JcyxcO1JQUR+BmwyxsxyNrt0TevncHC0MaYvcBrwOxE5NkLfdH6WPKAv8C9jTB9gN5Y5H450fhZsP/sZwOvRurq0pcVz2P71M7FcPocD9bA+Z+HkSttnMcYsAR7Cepn9EMudVRnhFM+fpVYoBR/GmO3A51i+toYikmcfagWss7crsDQx9vEGwFZnu8s5yeBo4AwRWQW8iuVC+huZ9xwAGGPW2b83AW9jKeuNItICwP69ye4eTuZ0eJYKoMIYM93efwNLSWTis4A1eM42xmy09zPxOU4CVhpjNhtjDgBvAYPJ3O/Kc8aYvsaYY225lpHC/0vGKwURaSoiDe3tOlgfmCXAZ8AIu9vFwDv29rv2PvbxicZywr0LnGfPVGiHFez5OjlPAcaY24wxrYwxZVjm/URjzPlk2HMAiEg9ESn2bWP5sBcGyRz8LBfZMysGAjtsk/kj4BQRaWS/HZ5ityUNY8wGYI2IdLGbTgQWk4HPYjOSKtcRZOZzrAYGikhdERGq/icZ910BEJFm9u82wNlY/5/U/V+SGVTxKFDTC5gDzMcaeO6029tj/YOXY5nKhXZ7kb2/3D7e3nGtO7DiEUuB01L4TEOpmn2Ucc9hyzzP/lkE3GG3N8EKpC+zfze22wV40pZ5AVDuuNZv7GdcDlyaov/HkcBM+zP2P6zZHRn3LFgTMX4AGjjaMu45bBnuBr6xv/MvYc0gyrjvii3DZCylNg84MdX/F01zoSiKovjJePeRoiiKkjhUKSiKoih+VCkoiqIoflQpKIqiKH5UKSiKoih+VCkoCUdEHvdle7T3PxKRMY79v4rIzQm+54+JvJ59zSNFZLhj/y4RuSWG80REJopIiaPtLBExItLVAznLROTXib5ulHv2FJEXknlPJTmoUlC84CusFaaISA5QCnR3HB8MTEmBXPFyJFYumngZDswzxux0tI0EvsRamJhoygBXpeBY4ZtQjDELgFb2giulFqFKQfGCKdhKAUsZLAR22astC7FyU80Rkfoi8qmIzBar9sKZACLykIhc47uY/Yb+e3v7DyIyw84lf7fbzd362G/TS0TkWbHy1k+wV8AjIkfZfaeKyCNi5bYvAO4BzhWr/sC59uW7icjnIrJCRK4P8/znU7UCFRGpj5XG5DIcSkFEhtrX8tVqeNleoYuIDLfbvhQrf76vvsZxUlUTYY69cvxBYIjddpOIXCIir4vIe1hJCcXxXAt8z2Lff5KIvCYi34rIgyJyvlj1SRaISAe73zn2ufNE5AvHc76HN0pOSSWpWMGnP7X/Byt1dhvgt8BVwL1Yb9BHA1/YffKAEnu7FGslpgB9gEmOay22r3UKVuFywXqheR841u7zo/3btQ/W23QlcKTd7zXgAnt7ITDY3n4QWGhvXwL80yHHXVhWUKEt7w9Avsuzfw8UO/YvAJ6zt78C+trbQ7EydrayZZ0KHIO1AncN0M7uN5aqFe7vYSUbBKhv/w2H+o475K6gahXsL7ESruUCzbHSRLSwz9tubxcCa4G77XNuAP5mby8AWtrbDR33ORp4L9WfNf1J7I9aCopX+KyFwViD3VTH/ld2HwH+IiLzsYqLtMRKGTwHaCYih4tIb2CbMWY11oB/ClZak9lAV6x8NU4i9VlpjJlrb88CysTKm1VsjPHJ9EqU5xpnjPnJGLMFK0lZc5c+jY0xuxz7I7GSHGL/Huk49rUxpsIYcwiYi6W8ugIrjDEr7T7OXEVTgMdsK6WhqUoVHczHxpit9vYxwFhjZRPeCEwCjrKPzTDGrDfG/ISVOmGC3b7AlsV3zxdE5AosxeJjE1aWUqUW4Ym/UVGoiiv0xHoTXwP8HtgJ/Nvucz7QFOhnjDkgVobYIvvYG1jJyw6jakAV4AFjzDMR7uvaR6wSjD85mg4CdXBPORyJ4Gu4fYcqRSTHGHNIRJpgZbztISIGa1A1InJrhOuFlckY86CIjMOyuqaJyElhuu52bEd6Ruf9Dzn2D9myYIy5SkQGAKcDc0XkSGPMD1j/q70Rrq1kIGopKF4xBfgZsNV+Q92KVcZyEJbVAFYK4022Qjgeqzykj1ex/NUjsBQEWFkff2P76BGRlmJnmHQQSx8/xphtWPGOgXaT00e+CyiO56FtlmIlZ8OW/0VjTFtjTJkxpjVWJbdjIpz/DdBeqmoJ++IZiEgHY8wCY8xDWEn6usYg5xdYsZFcEWmK5U6LORuofc/pxpg7saqW+VI0d8ZS+EotQpWC4hULsPzu04LadtiuF4CXgXIRmYllNXzj62iMWYQ10K01VRWoJmC5d6aKyAIsZREwGMbSx4XLgNEiMhXrrXqH3f4ZVmDZGWiOhXFY/nqwXEVvBx1/kzCzhexn2ItVV/hDEfkS2OiQ6UZf0BfrLf0DrOytlXYg+CaXS75t95kHTARuNVZK8Fh5xA48L8RSMPPs9uPtZ1VqEZolVcl6RKS+MeZHe3sUVm3cG2pwvRZY1sHJNZXJno30JLDMGPN4da+XaOxZZJOAYyLENZQMRC0FRYHTbWtgITAEuK8mF7Mtm2fFsXitGlwhVt3xRVhutkhxlFTQBhilCqH2oZaCoiiK4kctBUVRFMWPKgVFURTFjyoFRVEUxY8qBUVRFMWPKgVFURTFz/8H7XcsRh8oRPwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "grid.teff = 5780.\n", "grid.logg = 4.4\n", "grid.mh = 0.0\n", "wave, flux = grid()\n", "plt.plot(wave, flux)\n", "plt.xlabel('Wavelength (Angstroms)')\n", "plt.ylabel('Flux')" ] }, { "cell_type": "code", "execution_count": 11, "metadata": {}, "outputs": [], "source": [ "from starkit.base.operations.stellar import DopplerShift, RotationalBroadening" ] }, { "cell_type": "code", "execution_count": 15, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'Flux')" ] }, "execution_count": 15, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAERCAYAAACU1LsdAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4yLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvOIA7rQAAIABJREFUeJztnXeYFeX1+D9nO2WXtoBIWzpSBVaaotgRE40Go8QeS9TYY/yhJsYaazQx0Shi9KtRjDUWUFFRRASkdxEEhKUjVYqw8P7+mLl35947t+3euWXv+TzPPjvzzjszZ3bvfc+cc973HDHGoCiKoigAOakWQFEURUkfVCkoiqIoflQpKIqiKH5UKSiKoih+VCkoiqIoflQpKIqiKH4yUimIyL9FZJOILIyh7+MiMtf++VZEtidDRkVRlExEMnGdgogcC/wIvGiM6RHHedcBfYwxv/FMOEVRlAwmIy0FY8wXwFZnm4h0EJEPRWSWiEwWka4up44ExiZFSEVRlAwkL9UCJJDRwFXGmGUiMgB4CjjBd1BE2gLtgIkpkk9RFCXtqRVKQUTqA4OB10XE11wY1O084A1jzMFkyqYoipJJ1AqlgOUG226MOTJCn/OA3yVJHkVRlIwkI2MKwRhjdgIrReQcALHo7TsuIl2ARsDUFImoKIqSEWSkUhCRsVgDfBcRqRCRy4DzgctEZB6wCDjTccpI4FWTiVOtFEVRkkhGTklVFEVRvCEjLQVFURTFGzIu0FxaWmrKyspSLYaiKEpGMWvWrC3GmKbR+mWcUigrK2PmzJmpFkNRFCWjEJHvY+mn7iNFURTFjyoFRVEUxY8qBUVRFMWPKgVFURTFjyoFRVEUxY8qBUVRFMWPKgVFURTFjyoFJSzvz1/H9j37Uy2GoihJxDOlEK2Osp3J9AkRWS4i80Wkr1eyKPGzfsdern1lDlf/Z3aqRVEUJYl4aSm8AAyLcPw0oJP9cyXwLw9lUeKk8qCVKHHW99tSLImiKMnEM6XgVkc5iDOBF43FNKChiLTwSh4lPnJzrAp2+w8eCtvn4CHNsKsotY1UxhRaAmsc+xV2WwgicqWIzBSRmZs3b06KcNlO8IA/c9XWgLbHPv6WDreP58OFG5ItmqIoHpJKpSAuba6vnsaY0caYcmNMedOmUZP8KQnAqQC+XrmVEU9P5V+fLwdg5ZbdPPHpMgA+WLg+4Lwdew5w1lNT+P6H3ckTVlGUhJFKpVABtHbstwLWpUgWJYjKQ1Vuo/U79gKwdOOPABxwuJQqHcrjo0UbOO/ZacxZvZ1/TFyeJEkVRUkkqUyd/S5wrYi8CgwAdhhj1kc5R0kSzsH+kF2dL9e27XKkysibu3q7bRnsYevuqumr2/ccSI6giqIkFC+npIbUURaRq0TkKrvLeGAFsBx4FrjGK1mU+PHNPnLyv7mWIecLQgOs3b6XOau3BygEgE+WbGTxup0BbYc0MK0oaY9nloIxZmSU4wb4nVf3V6rPxp37eHbyCtdjCyp2sGHnvpiu882GnXQ7vASATbv20f/+T3ng7J6M7N8mYbIqipJYMq7ymuI9V/9nFrNXb3c99vN/fhnzdZwuqDVb9wBw93uLqNi2h1tO6YKI21wDRVFSiaa5UEIIjgfc9N951brOuPnreWX6agC/Ath34BBPfvYd7W4bzx//t6BmgiqKknBUKSieMenbzdz+tjXw57hYBf+ZtjrZIimKEgVVCkoIXoSDc2LwFL03bx279umsJUVJJaoUlBBWbkn8wjNxXasI89ZYsYuZq7Zy3dg5/P616rmqFEVJDKoUspBXpq/m3GemJvWe4WLKb82uAGDE05Y8ExZvTJZIiqK4oEqhFrDvwEFufWNeyFqBcNz+9gKmrwzNVfjajDWUjRqXaPGA8ErhkLEC0k6M0fUMipIqVCnUAt6es5bXZlbw8Iff1Og6z325MkESBXLwkCHcOJ+bIyxctyOgrd1t45m8TBMfKkoqUKVQC/AFcYPfuOPFq2UD3e780J8qI5i8HKHSJT33B5p9VVFSgiqFWoBvuueunyrjOi/YTeM2bTQR/FR5iHAZLnJzhQMuKTVemb6a68bO8UQeRVHCo0qhFlDdwTy4ZoKXnvwde92nmuaKBGRddfLePE2aqyjJRpVCLSA3lkUALjiVwI69B1iyfmfYvjXl4n9/7dr+1OffhVUKkdi7/yD3vr/Yv65h8bqdLN/0Y41kVBRFlUJGM2X5Fj6OMIXz0CHDXe8uClvwxunnv+HVxLhqnru43L99+/CuMZ3z2syKmK+/Y+8Bduw5wOuz1vDclysZ/YWVuG/4E5M56bFJ8QmrKEoIqhQymPPHTOeKF2cGvGn/ZfwSTn9iMgBLNuzkha9Wcc3Ls13Pd4YUPl8a32yfX5W3cm0/umOpf7tPm0ZxXTMW+t37Mb3vmcD+SuuZ/zFxOZ99synh91GUbEWVQi2gIK/q3zj6ixUsCqpjsGjdTlcXzQ+791M2ahxffBv/9M8jWpS4tuc5XFl51XRrOfnPtO+ZvXqbf9+XedWZYfUhx1RcX5W4bzbsjGhFKYrijiqFWkDzkqKofaat+CGkbb6dYiJc7YRI9A1jBeTmCJP+MJR3fnc0+bk1/3j98X8LOfupr9iwY1/AwrpvHPEPZ4ruQQ9MBGDY3yZzxYsza3x/Rck2VClkKM6ZQ763Yydlo8bx476qKaoXPhca6D1oQt+6Y6G8bSN6t27oekxEaNukHr1bNyQvN3FTXAc+8GnA/uuzquIQwbOonPzomKZrjGGnJtxTlIioUshQ/jPte/92uHoHHy6KvACsutUxi/JzY+qXl5Ocj1dwAj9fkj2Aa16e7VcaT33+Hb3umlCt2U6Kki2oUshAZq/exp/fXRS13979ByMe9y1ei9f1H6thkZ9ASyEenDGVL77dTIfbxzNn9TYe+WgpANv2xJYjSlGyEVUKGcjU70LjA268OmNNxONz7JKbvpk8wQzt0jQuuR4Z0Stgv7rrJ2pK5aHQ5znrqa/82zVNB6IotRlVChlIdbNRBC9Oe+GrVQB8FUbJNCsu9G8/eHZPx/3dBTinvHXAvldpM6Lx7cZdEY/f/d5iAMZ+vZpNu/ZFjEkoSrahSiEDqe5ge9rfJ7u2t2jgPnvJWRgnx/HW73b3Pm1CA8+pshSilfls26QuFdv2cNtbC+h//6fc+N+5SZJMUdIfVQoZSKLH2vU79rm2O3XPz3q1oE3juiHtAG9dM5gXLu0fcn6qLIVo9GzZgK8d9SScOZb2HTjI7W8vYN320BldipINqFLIQFIx2Bbl5fKnn3UDrCR2Tvq2aUSDOvkh56TKUojG+/PXc3NQ2c/b317A7W8v4KNFG3hl+moGPziRzbt+ck3rrSi1mbxUC6BkBiJVs5ViXdeQSRXUXpluuZya1q+Koxx1/yec3bclj/3qyFSJpShJRy2FDGPH3gPcN25JUu4VEEcQ8a9riNUAyMT47d8/XRaw/9bstUxf8QPvzlunAWklK1ClkGGMqUZKiuoSPPY3K7HeoruGyXsUTGn9Ajo2q59gqZLPuaOncf3YOfS+e4I/Vbei1FZUKWQYG8IEhb0gOHbRt00j3rx6EDec2Cmm80WEO4Yf4YVoKeHHnyr5x8TlqRZDUTxFlUKGUZlEF0aPlqEWQb+2jeMKIKfpBKRqM/qLFfz62WmpFkNRPEOVQoZR0zH2hx9/irnvr+zFaE3qFbgev/6EjlEtgXSdlloTvvruB/43Zy1fLtuSalEUJeHo7KMso999n8TcV0SYetsJ1C1w/5jcfEqXqNeojUoB8C94W/Xg6SmWRFESi1oKiiu+UpotGtRxXYMQK+mmEzo3r17g+xdHHu7aXjZqHGWjxvGwo9CPomQyqhQUV7of3iAh1/FCKTSqG7+SumBgGwAa1g11hRXYxYBaNqwT9n4Pj+jNSUc0C3v9pz7/Lm6ZFCUdUaWQYSRrlXBw8rzqUlP30Yh+obWg86pR0W1Y9xbWuTkScs1D9iK7O3/ezfXcNk3qUZCXw+ZdkeMxW3fv5+Ahw7vz1rF1t6bnVjITT5WCiAwTkaUislxERrkcbyMin4nIHBGZLyLDvZSnNnB81/Bvq7ESyzi9KcoAGCs1VQpuJT2rU/t5QPvGnNOvFQ/9shdn920ZcMw3n6u4KHKILdpK7r73fszbc9Zy/dg5HP/o53y4UFN0K5mHZ0pBRHKBJ4HTgG7ASBEJfhX7I/CaMaYPcB7wlFfypDvLN+3i+x92R+2XLEvBl/yupoQTN9YArZtR4FbmM1pBn/zcHB45pzetG9cNUTS+dBwlRe5uqSEdSwP6ReLxj78FrJXnV/1nNjv26mI3JbPw0lLoDyw3xqwwxuwHXgXODOpjAN9k+AbAOrKUkx77guMe+Txqv2SlEzq1+2EJuY7z7fqXfUNdQVHPd5mE6yvz+bvjO3DhwLYA/LyXeyD4l31bMfH3xwWdH3hN39KP+oWhlsJntwzl5pM7B/SLxNqg7Krb1I2kZBheKoWWgLP0V4Xd5uQu4AIRqQDGA9e5XUhErhSRmSIyc/PmzV7ImkEkRytUx0XjhvMytw6LPoX1/37Tn/evO4bfn9w5xJ3jq/vgky0W11T7pvVo3zRwxlG42tH1XdxH7Urr+XNAhct91NelloSPDTuTtwJdURKBl0rB7Rsb/K0aCbxgjGkFDAdeEpEQmYwxo40x5caY8qZN4ysRWdtIxILmWIb73ATVV3ZaCs4MpOE4rnNTerRswHUndmLBXacGHPMpAZ8LTUT88ZFwf5ZoKb3/c9kAv/JxWgof3jiEN64aFHDeoTBmWrvS8NNczxs9jXfmrg17XFHSDS+VQgXgrM/YilD30GXAawDGmKlAEVDqoUwZT7iBKdEE10yoLk5LwXfJc4PKdsZLrt9SqGoL9vePvrAfD/2yJyP7twmVyf7Ud25en2M6lfLSZQO4ZmgHivJz/X26HlZCeVnjgPMePac3xwfVre7YrD73nNk9orw3vDrXP5tr3prtrN+hBXyU9MVLpTAD6CQi7USkACuQ/G5Qn9XAiQAicgSWUsh2/1BEEmEpxHKNRAW0nS4eEWHR3afyF0e95+rgtxRcIw5Vfc49qk3E5/CdfWTrhtw6rGvU+/Zo2YDnHRXm3rpmMJ/cfBz1CvP4c5jprD5WbLYmEZz55BQGPTCRh3Sxm5KmeKYUjDGVwLXAR8ASrFlGi0TkHhE5w+72e+AKEZkHjAUuMZlUmSUFHPIwIV7Xw4r924mKKQQbHPUK80IG6q/vODGua/lcUgGWQpT7OknUJ8z5N7r06Ha8efXgsH2XbtjJxG82+vf/pYvdlDTF09xHxpjxWAFkZ9udju3FwNFeylDb8NJ91KR+1WpfLyyFYObdeQqSE34qKAQO7r5tn2g5ORJ27UB4G6JKKYQTrXGYBIDBBP+N+rVtxPOXHkWX5sUUF+Wxa18lgx+cCMATLim391ceoiBP148q6YUmxMswvMyc7RxIYy25GY1ISqFBnOkqfPL54h3OSxcEL2iIIL5vnUPzkqKQY+OvH+IvJhQNN8V5fJeqxYXRKrVd8/Isxlx8VEz3UpRkoa8pGYaX7iMv8GKtXY7LlNQjWpRw/1k96N0qes6mzs2L+ctZPfnbuaG1l7sdXkJpDLOkILqLLVo6jk+WbGL3T5UAzFi1lacnqUtJST1qKWQYyZp9lChiNTgm3HRs1DQTZ/dtyd8+WcZh9ht+8Jh8/oC2fLRoo8uZofx6QOispHiJtk4ilrhM9z9/FHLOOeWta5SZVlFqgloKGcaotxZ4dm0vMprG6obq3LyYFg3cs5T6uHpoB765d5g/9hFpUE5GMpBwi+CqjscvxX3jltD77gl8syExCQkVJV5UKSieksgiO4JQlJ/rek2f/ZTMyWtRdEKNgvXD/jaZslHjqDx4qNrXUJTqoEpB8ZRExhRygmYfHTImrHWTqEB5JKJZComQoeMdHzB+wXpd8KYkDVUKSgBDOiV2QXlNLYWZq7b5t6vWJ1i/nUZBKpa3xGIJ/PH0IzjHpSZEPFzz8mwGPTARY4xaDornqFKopbRtEn/q67vO6M6/LzmKhXefGr1zjNT0ZXmxo9iPBG0cMpHXI3hNLErh8iHtOaJFSUDb7D+dzLBqZKF9dMJSOt7xAXv2V8Z9rqLEiiqFWsrpPVvE3b9D0/rk5+a4ppCuLol041QtXrM23GZi+dYGJCp3UyRijRkEy9m4XgGPu0yHjcaTn1lTVk94dBIzV22N+3xFiQVVCrWUEyPUE3bFozE0kTGF4PQWhw6FxhR8SiFaEDiZuK1arlNQlXyvcb2CuCy7DTv3MeLpqepKUjwhjb46ihv7Kw9x03/nxlSVzUl+bk7UamTJIJGzj3zk+i2Fqjbfy/hVx3UAoFuQy8YLimO0qM47qg3XndAx7PHLjmnHk7/uC1huvwk3HcvntwyNet23ZmtKbiXxqFJIc2as2srbc9Yy6s341yfE47rxSn14ufbhkDEheYyO79qMVQ+eTsO6seUvqgk5MZpBBXk5/uptbpTWL/Cn3ijIzaFz82LKSutFve5d7y0C4MtlWzjl8Um8NmNNlDMUJTq6ojkNmL16W9hjvmHHxFhxrW5BLnv2HwQgP0dIdTFILyyFqtlHmbO6201BL71vGK/PrGBEv9YsWLsDgML82N/T9uw/yNw127nguekA3PrmfGas2srVQzuEVJtTlFhRSyFF7Nh7gOMf/ZzPl27i7Ke+8rf/dcJS9lc6fMW+ymIGfqo8GPW6zkE4Wu6dZOCNUrB+H8wgpeDDOeW3MC+XCwa2JTdH6NK8mJ4tG3DvmT38x8/ofXhYN9iIfq1oUq+AXzw5JaD99VkVnPDXSVz7yuyMy5OlpAepHzWylKnf/cDKLbv5y/glAe3/mLic12dVuQF8Uy4N0OWPH/rbv77dvQZBoR3UzBGJK82CV4u9vEiIl6nM+uNJjLm43PVYnYJc3rvuGPq0aeRve2JkH8bfMIRx1x8T0r+4KI+d+w6Evdf789cz5bstNRdayTpUKaQYt5fdfQeqLIUcv6UQ1DHMYPuPX/fhmqEd6H54id9PnUq8XEcgCOcPbEPDuvkMj3MKbipoUr+Qwrzc6B2D6H54aObXeWu2c+BgZEvg3bnr+HbjLjbu3Bf3PZXsRWMKaYhTAfiCmdFy8/to2bCOv7RktDQMScFjvdShaX3m3nmKtzdJA84tb82Xy7ewdvteRCz3YzRen1XB67MqAPj7eUdy5pEtvRZTqQWkwaihBONUAD5LoTJIKYR7A3f68OPRCV6N3V64j7LRU/7QiF78+xKrIE+uCH8/r09c59/w6lymLFd3khIdVQppSGAJSmvH6SroX9Y4pqmescRhzzuqdbzixYWXiemSsGjZlTEXlfPCpcmvmOZzB+aIhJQMvfGkTlHPP3/MdKZ+9wMzdDW0EgF1H6UhTisgx68UAmckhRsPnQPlb49tz5/eWRTbPT0aYGt62ZO7NefjxbEVzkkWJ3VrnpL75tumX05O4CrpO3/WjYsHl7Fyy24uHlxG60Z1eX7KSp76PLSS28hnpwFwbOemjLmoXGtEKyHoJyLFxOoKcaY0EMK/gTvbLxxUxj9/bbkZwiVg8/ptu6ZTUof3jD9xXG0l12Ep5DumGzeok09ujuVS6tumEU2LC/2WRO9WDXj6gr4h1/ri2810/uMHlI0axztzdWW0UoUqhQTz9cqtrNj8Y0DbJ4s3smFH8AyQ2NSBL5maM6YgEr4sZ7ghOHx8wVutkCoXT20k3w7Q5IpQ4FAKbon5fJ+P/u0aM6xHC777y3C+vuNE19ToN7w6l7lrtnsktZJpqFJIML96Zion/HVSQNvlL85k4AOfcsKjn/vbfGN6tECsbybSgQBLQQLdSQ7ifTP3dU/XsTsD16d5hm/wz82VgLxWbuk2fB8P37HcHKFZcREvXTbA9dq/eHIK142dowveFFUKyWTFlqqkdle/PDtsP2dKC9931BloFoG6+e7hoLCVyMIM+14Pul6saM5WReH7LOTlSIB14LZI0WcpuKUQb1LPPS/Ue/PW0f728cxctTWjUogoiUWVgsdE+3JFW9zlO/1ApXNBm9Cgbj69W4Uuagp3tWi5k7yaJVTTy3ZuXhzSFmseqNpGaf0CLhlcxkuXDQj4f7kpXn9dCReF4YxHfHzTsSHHRzw9lee+XJkIkZUMRJWCx0SzxqMNmuFiCgBNiwtdLhiPdN5TU3F6tAxVfFWutzR7WI8REe46o3tIJTe3gd+nFNyU/bMXWak2/nj6EXRqXsxhJUUhff464VsAXpq6ik27dEV0NqFKwWOCLYVZ34fPiOrGeaOtKYRuyd/cjJBwlkeqylZ6MXD7FGW2KYVwxOs+6tmqAasePJ3Lh7QH4PWrBvHRjcdy+THt/H32HjjI23Mq+NM7izjxr5N4bMJSdSllCaoUPCbYUvjlv75y7xiFgCmpEQbDcIHrqO6jakkVnUSM25NvPZ43rx7k3y9rYtUaaNc0es2BbMAt0OxTCrGk4m7duC5dDivmjz/rFtB+03/nAbBrXyVPTFzOqzPWhJ3goNQedPGaxyTK/+1ULr4hwG3ADVYYqSxsD4mJVbRuXJfWjavKVZ7dtyVlpXXp68goms24WQNXHtuBrbv3c9Ggtgm7z21vLaBeYR6ndGtOUX78if2UzECVgsdEs7irM2j6TnF3HwXdP2qAOe7bpxwRoV/bxqkWI21wW4PSoE4+D5zdK+H3un7sHAD+MbIP9QvzOL5rnLXAlbRH3Uce47bIzLmCtDpjcqRzMnGQV2pGIrPh3n1Gdw4rKeLnvQ8HrDxPblw3dg6XvjAjYfdV0ge1FDzG7W3+hlfn1uiavgBrItxHfvlUmWQsiSywd/HgMi4eXAbAg2f3pF5hHh/eOIRhf5vs2n/Ljz+xc+8Bdu6r5MjWDRMniJIyVCl4TLh0FDUhovsoaHCPNaaR6tiDUn28moVVr9AaHroeVsJLl/Xn758sY2bQ7Lny+z7xb4+5qJxerRvQrDh0iquSOXjqPhKRYSKyVESWi8ioMH1+JSKLRWSRiLzipTxec+rjX4S0RRuSq/d9Dn9SvJdTd1Pmk4xiSkM6NeWNqwdH7HP5izPpf/+nnsuieItnnyYRyQWeBE4DugEjRaRbUJ9OwG3A0caY7sCNXslTEz77ZhPnj5kWdZ720o27/NufLN7I1t37MVFm8FVnUPbnK0qk+0jJWJJZYC9cbXAnM1Zt5b73FwdMo1YyBy8/Tv2B5caYFcaY/cCrwJlBfa4AnjTGbAMwxmzyUJ5q89uXZjFl+Q/8VOn+Id++Zz+3vbUgoO3yF2dy6QszYlgfEH7QXlCxI8w5Fr4B/aQjmoUc8xGz+0gthozFbUWzVzQrKeKUKPUkznl6KmO+XEnHOz7g5enfJ0kyJVF4qRRaAmsc+xV2m5POQGcRmSIi00RkmIfyVJ8o37nHPv6WsV+vDmlftWV31DQXC9a6D/wA2/fudxcnSJ7zjmrj347fv2ynQ4jzrHjo0ryYe8/s7uEdshu3Fc1e8peze4a0TfrD0IB03j7ueHshZaPGsWT9zmSIpiQAL5WC2yc1eIjMAzoBQ4GRwBgRCZnCICJXishMEZm5efPmhAm4bOMu9od5+3fi6/Pm7ArX4+GCyTv2HmDM5BVxy/WfaaspGzWObXvci7NHsi6CdUKsAWQvLYWPbjqWCweVeXeDLCfZ6T5K6xcy3eFG6tSsPm2b1GP2nSeHPef5KZpgL1PwcvZRBeAsANwKWOfSZ5ox5gCwUkSWYimJgAnQxpjRwGiA8vLyhHjBN+7cx8mPf8HI/m14wOXNx42p3/3A+QNCV4hGGnjdSiJGY/XWPQBs2LHX9XikMSD+2UfqN8p0kuk+8tG8pIjHz+1NSVE+A9s3AaBeQfhVzq/NrOD9+evZs/8gvVo14NZTu3KMS8EfJfXEZCkEB4jttqFRTpsBdBKRdiJSAJwHvBvU53/A8fb1SrHcSfG/WleDHXutt/B0LmL+1mz3MomR3gx1amn2karEgGf1acWJRzT3T111TnJo07guI/u3ZkC7qpXne/YfBGB+xQ4ueG46i9epSykdidVSeE1EXgIeBors3+XAoHAnGGMqReRa4CMgF/i3MWaRiNwDzDTGvGsfO0VEFgMHgT8YY36o/uMklr37D/Lxkqqi8eHeub36Tn6zYZf7gQj3C35pjNl9pMokY8nLTb//3Re3Hu/fLhs1zrXP8CesBXG/7NuKW4d1oblLCm8l+cSqFAYADwFfAcXAy8DR0U4yxowHxge13enYNsDN9k9KiPR1um/cYl6eXhVADhd/SPZXMnKai9RUWFNSh1tCvEzizdkVTFvxA1NGnZBqURRiVwoHgL1AHSxLYaUx0WbgZz4bdwYWF9l34KBrP6+qloXDd7+2dgrphnXzq44F9fXFFMKLqNoi03FLnZ0qZv7xpGrVeV67fS9lo8bRrUUJY68YyE+VB2mmlkNKiHX20QwspXAUcAzWQrQ3PJMqTQgO4LVxpG/2UbFtD3NWx1c4p6b4pPp/p3Xh35eUU15W5beNt0Zzu1JLsfRwKe2pZAbJnpIaidL6hVEH84d+2ZMTw2RXXbx+J73vmUD/v3zKLa/P80JEJQqxWgqXGWNm2tsbgDNF5EKPZEobgtMHlJeF5u8/5qHPkiWOH9/AX5iXywldmwcdi2+AGNS+lI9vOpaOzeonSjwlyaSTpeDGS5f1J0eEeRXbOemI5nRuXkxxUT6ffmOtVX3sV725+bVQBfDGrAoWrdvJvWd2D3jxUbwlVkthk4i0cf4Ak7wULB0IthTSxS/vNgQ8fm5vuh4WWuQ+6rUEOjUvTroLTEkc6R5TGNKpKUd3LOWaoR3p3Nz6jJ7W4zCevaiciwa15aw+LVn5wHDXc5es38mIp6fyz4nL2Lv/oJYETQKxWgrjsJzPghVTaAcsBWr1MtW0VQoug8BZfVpxVp9WcV9L6xxnPqlYp1BTRISTuzXnZEfKjKcv6Mv8ih2ua3senfAtj074luO7NOX8AW05KUqqDaX6xGQpGGN6GmN62b87YeU1+tJb0ZJDpHE+eMD8dtMuNu3aF6Z38kjkON5e6xxnPJmoFNwY1qMFI/u3idjns6W3S08WAAAdV0lEQVSbufzFmYyZvKJaAW0lOtVKc2GMmY0VdK7VBAfwnpm0Ii1SA8ezpiCSdTOs+2Faa7cWkO7uo3goLrKcF6X1CyL2u2/cEtrfPp5XpofmHFNqRqwrmm92/Nxi1z1IXBKiFDDiX18BVf75NVv3sGFHoBUQLoD3hxTPiqjWGCCw/P7TuO6Ejv6mTs01uFwbSPdAczw0rFvAmIvKeeWKgf62/145MGz/299ewE3/ncud7yzkoFoOCSFWS6HY8VOIFWMIToOdUezcVwlUuY+GPPwZAx8ItALClTl8fZZ7YrxkUd0hIC83x+9qaN24Djec2ClxQilKgjipW3NaN6qa/j2gfRPm3XlK2P5vz1nLi1O/p8Pt49m22z2zsBI7MQWajTF3ey1IujB52WbqFebRt02jiBWttu3eT597P06iZFXEYynk25qtKM9yE/lqNGza+RN5iSzuqygJpCg/8LPZoG4+OQKHDEy97QTWbtvL3e8tDkk93+fejxl//RC6HV6STHFrFRGVgoi8R4RYrDHmjIRLlGSWb/qRpydVzXa48LmvAVj14OkRZ+ZUbHPPYJoM4okpnNytOdef0JHLjmkPwCI7CVm4gkGKkg6ICF/84fiAF6DTex3Oe/PW0ahuAS0a1OG24V359bPTyRGY+Puh3Pjfucxds53hT0zmZ71a8M2GXdx8cmeG92yRugfJQKJZCo8mRYoU8+AH37i2R3qRrjyUukE1HkshN0e4+ZQu1TpXUVJJmyaBGQQe/mUvbjixk39yRI+WDejQtB73ntmDstJ6jL6on38iyPvz1wNwzcuz6dGyhOtO6MSp3Q9L7gNkKNGUwkpjTNaG9/cdCD/wpzKoVZOFZuEKAilKulOnIDdg5X1JUT6f/n6of79p/cKA/v3aNmLW99tYuHYnv31pFr89rj1DOjbVOg5RiOZU/p9vQ0Te9FiWtOPwhnXCHkutUqj+uVpLXamtiAhjHbOWHhnRK+D4M5NWcMFz09O6hko6EM1ScA4/7b0UJB0pzAuvM88dPS2JkgRSEw+QWgpKbaZv26pqviV18l37nPP0VAAa1yvghUuPolerkArAWU00S8GE2c4KpizfkmoRXKmJpVCppoJSiynMy+Xa4621OCVF+bx59WCeu7jcte/W3fs5459TOPPJKeza514PPRuJphR6i8hOEdkF9LK3d4rILhGp9bX0fFkc042aVEnTBT61h4m/P45XIyzsylauP7ETX99xIgV5OfRr24gTj2jOlceGOjpeuqw/APPWbKfnXRPYs78y2aKmJRHdR8YYzYGQhtTEUjhwUJVCbaF90/q0b6qr0oMpyMuhWXFgTYfLh7Rjw4593H9WD56dvJKfKg8yuEMpfdo0ZM7q7QCc8OgkNuzcx5BOpTz2qyNpWlzodvlaj65eykBqElP4y9k9EyaHomQKzYqLeGJkH4qL8rn55M7cdtoR5OYIb19zNAPaWbUaNtiVFicv28JR93/CJc9/nZVuJVUKYVi1ZXeqRQhLTaakDrGn4+WnYbF3RUkF9QvdHSafL91Mz7sm8MOPPyVZotSiSiEMl74wI9UihKUm7iPfuToJSVEs7jojtCzMa78d5N/ud98nvDh1Fc9PWZkV6bpjLbKTdew7cDDVIoSlJoHmmpyrKLWR1o7a6yP6teLPP+9GcVE+vz2uPc9MWgHAne8sAuDu9xbTtLiQ5y85ih4ta2ddc1UKYVi/I/XFdMJRE0vBl2W59r/vKErszL3zZESEBo61DVcMqVIKTjbv+omf/eNLrj+xE02LC7lwYNtkiuo5Wec++seny8LmOsoUavKur7WYFSWUhnULAhQCQGn9Qp6/tKqWWPvSwCqFT3y6jD/9byFlo8bx/vx1SZEzGWSdUvjrx98GZEXNRBJiKWhQQVGiMqRjKcN7HsarVw7k/37TP2y/a1+Zw8mPTWLzrswPSqv7KAOpydu+71y1GBQlOnm5OTx1fj+gKs54Ru/D6dWqAfeNWxLQd9mmHxkzeQVrt+9l7prtvPbbQRHzp6UrqhQykJq5j2p+DUXJRoryc/n8lqG0bFSH/NwcLjumHe1uGx/Q55kvqmIQgx+cyJBOpfRp3TAgfX26k3Xuo9pATd7yc/yWQqKkUZTsoay0nr+aoYjw5K/7hsQanExetoUnJi73VzzMhCmtqhQykBqtU/D/Vq2gKDXl9F4tKAtSCo3q5vP/hnUNaLv8xRkcOHiII++ZQNmocbw4dRWQnkpClUIGcGr35gH7NRnOc9R/pCgJJTjJZKdmxVw8uC13n9Gd2X86GYCNO3+i0x0fsHOflXTvzncWsXLLbtrfPp7y+z7m86Wb0mbyhyqFDKCkKHCqXCJWNKtOUJTE4KtR4qu/kpcr1C3I4+LBZTSuV8BFg9zXMRz/6OcAbPlxP5c8P4N2t41Pi9T2qhRi4D+XDUjp/fOC8hTVxPXjexnRmIKiJIZKO/Owr6hPbk7gl6t1o7oh54Sj4x0f8PSk7ygbNY7XZqxhw459SbcgVCnEQKpruuYEjeA1GdAL7LeZYJ+noijV46A9aPsWv/kC0T6K8gP3n7mwX8Tr+RbX3vrmfAY+8CntbhvPuu17EyVuVLJKKXy8eGOqRagWeTnBlkL1yc0RVj14Opce3a5mQimKAlTFFEqK3Gf45wUpiVO7H8bgDk0AWHLPMNdzghXJ4AcnMnfN9pqKGhOeKgURGSYiS0VkuYiMitBvhIgYEXGvm5cgrnhxppeX94ycYKWgvh9FSRt8MYV6dgru4G9nrv19rVeQy4w7TgLgX+f3Y8JNx1KnoKqO2VvXDPZvf3PvaSH3+cWTU/jqO+9LBHu2eE1EcoEngZOBCmCGiLxrjFkc1K8YuB6Y7pUsmU5uAt1HiqIkFt+00noF7sOpL8bQrKTIX82tQd18GtQNnEByZKuGAftPjOzDm7Mq+Hnvw1m2cReTvt3sj194iZcrmvsDy40xKwBE5FXgTGBxUL97gYeBWzyUhW837vLy8p4SHLjSNQaKkj74Ygo+l0/wS5tvoki0b22wR+CM3odzRu/D/fu3DT+iZoLGiJfuo5bAGsd+hd3mR0T6AK2NMe9HupCIXCkiM0Vk5ubNm6slzNokBmqiMf76IXH1D1EKqhMUJW3wvb37JnEED//R1gZ1alafsibWDKUxF5Xzzu+O9kLMmPHSUnD7E/htHxHJAR4HLol2IWPMaGA0QHl5ebXsp3QaR9s2iX2KGrhZCoqipAu+mEKVUgjEN1Ek3Pf2oxuP9V/jpG7Nw/RKHl5aChVAa8d+K8CZdLwY6AF8LiKrgIHAu14Hm9OBeN/0EzklVVGUxHL/WT3p26YhnZoVA6HfT99LXbgJIjk5EjJDKZV4KckMoJOItBORAuA84F3fQWPMDmNMqTGmzBhTBkwDzjDGZOYUoTiINyYQ6j5SraAo6cJRZY1565qjqZNvzSQKmX0UxVJINzxTCsaYSuBa4CNgCfCaMWaRiNwjImd4dd9UM6Bd46h94h3TOzStX01pFEVJFlUWQWzt6YqnNosxZrwxprMxpoMx5n677U5jzLsufYems5Xw+Lm9aRQ0hcyNYzs3jdpHBC4ZXBbTfUf2b0N+boZ8mhQliwlOR+Nvz7FnJWWIrZA+jiyPqanL5aw+rTiqLLoVEJMsCHed0T2kfdWDp7v0NiExBUVR0g9/rZLg2Uf2KJspX+OsUQrpRDwfDmMy58OkKNlMuAJWfkshQ77IqhRSgNtHY+wVA8P2V0tBUTKHcIvX0qVeQjSyRilUd1h9YmQf//aQGOIFMcniMsiX1HFfMlKUnxsivOoIRUk/DO6Dfn5OZg2zmSVtCujXtpF/+4IBbRJyTbcxPVwQ6pZTu2RIeEpRsht/rZKgb2zwlPJ0J2uUQnXfrp3/z0T5BN0uE+7S9Qvz1H2kKBlAVbqGwPYMMxSyRylUFy8G5HiVi+oERclcfJZDhoQUVClEI1lv6ZFuo5aCoqQ/vkBy8LfV9/U9lCFaQZVCFBLtDvzFkYe7tkda2KIqQVHSH9+YH/wSV9cupFNWWi/ZIlULVQpRCA4S/a+GaW1H9Gvt2h7RGAiefaRqQlHSjkq72E7wyuZWjeoy5qJy/vqr3qkQK268TJ2dVlR3IE30gpNwl4t0F3UfKUr6061FCQAndg1Nf50OKbFjJWuUQnVJtPsorFKIcB9VCYqS/nQ7vITF95xK3TBlOTOFrHEfVfdlO9EFbga2axL3OcFl+hRFSU8yXSFAFimF6hLsuoll/kDzkkK6HlYc0i5SvQFeVYKiKMki89Wax8Trz88RYfrtJwHQ6Y7xHDhY82lomZJIS1GUzEcthShEerF3swbyHCfEMy05Ut9gnaA6QlEUr1ClAJwfIadRSH1kx3Zwoe6LB7XlwkFtw16rumO56gBFUZJF1iiFiMsAIs38CTrmfKEPfru/+8weVlZTm0fPiX1eciSjQqekKoqSLLJGKUQi4mriCANytGXrv+jTMu5rFeaF/ktC3EcR76ooilJ9VCnUAKdOqOnLvE8xjb6o3N9WXGjNA1BLQVGUZKFKgfDFMaJxWo/DEiwJFNmWQrPiQr6+46SEX19RFCUSqhSILwWGs+fvju/I5FuPD2kPR2n9grDHfIrJTp9CWWk96tiJtEKC3Wo4KIriEdmjFGIcSP99SXnE406bIidHKKmTb10+ykj98IhevHHV4Kj3b9mwDgBDu1SV/lQloChKstDFazXEN2BHG7d/Ve6eHdWHLz7Rpkldvr7jRJrWLwy5R3BfRVGURJM9lkIEnDGFWF1JvVs1sPvb51Xjbf6fv+7j2t6suCjA8tBAs6IoySJrlEKkwT7gzVugY7P6sV/XHrCrk5q7naPoRsQVzUH7aigoiuIVWaMUYiVHhA9uGBL2uAkavSVkI/Fo7iNFUZKFxhSCyBHIz83hllM6U+rw64fgsxBijCnUBNUJiqIki6xWCq0a1aFi297AGUX2CHztCZ3iulZ1Bm6n0VG/MPy/IsR9pP4jRVE8ImvcR26Ddic7dhCwMjnO6/rOrU5M4eChqhu3aVI3bD8NNCuKkiyyRim44RuSnWNuvP57t2s4iRS09hX6PrJ1w4j3CE3Kp6aCoijekNXuI7dx3K1GgpPg4bgg19KrZwUlv/Px0Y3HhgSnfdQrtFYst3fMQnIjpPqb6gRFUTwiq5VCMEe0KKFRvfCpKJz4humCvBzm33UK9cLUZrVqPLubEV0PK+HZi8o5pmNpNaRVFEVJPKoUEkBJUX61zz25W/OofapT11lRFKU6eBpTEJFhIrJURJaLyCiX4zeLyGIRmS8in4pI+LJlNZXFXT4g/d0xqhIURUkWnikFEckFngROA7oBI0WkW1C3OUC5MaYX8AbwsFfyJIpUKBCdfKQoSrLw0lLoDyw3xqwwxuwHXgXOdHYwxnxmjNlj704DWnkoT1hKiiwvWpfm8aS38EqaUEIDzWlu2iiKkrF4GVNoCaxx7FcAAyL0vwz4wO2AiFwJXAnQpk2bagnjNoz6BtfWjesy9oqBUaeGpgo1FBRFSRZeWgpuY5nrK66IXACUA4+4HTfGjDbGlBtjyps2berWpcYM6tDEX9QmEk3s2Um9WyVPgWjuI0VRkoWXlkIF4Cwi0ApYF9xJRE4C7gCOM8b85JUwifK4lJXW4/3rjqFz88jrGRKJ1lNQFCVZeGkpzAA6iUg7ESkAzgPedXYQkT7AM8AZxphNHsoScRVwvC/iPVo2oCAveYvB1U5QFCVZeDayGWMqgWuBj4AlwGvGmEUico+InGF3ewSoD7wuInNF5N0wl8tqQgLNKZJDUZTaj6eL14wx44HxQW13OrZP8vL+gcIk7U4JR0MKiqIki6xJiHfQxRFflG8FlvNz0vvPEBxo1gXOiqJ4Rdakuah0pKnu1aoB8yt28Oefd6ddaT3O6uuezC5dCLYUItVeUBRFqQlZM7ocPFilFG4ffgRdmhfTqF4Btw7rmkKpYsMXU8jPFf7fsK5cMNCzbCCKomQ5WaMUnJZCbo7EnA01ncjPzeHyIe1TLYaiKLWY9HamJxBnlbNMm+fvs3LyNJigKIrHZJGlcMi/nWm5g4qL8hjW/TAuPbos1aIoilLLyRqlULFtb6pFqDY5OcLTF/ZLtRiKomQBWeM+WlCxw7+dWXaCoihK8sgapfD0hf3oHEdqbEVRlGwka5QCQMO61oyjDAspKIqiJI2sUgq+uTuRkuMpiqJkM9mlFHRGp6IoSkSySikoiqIokckqpSA+B5J6jxRFUVzJLqWgOkFRFCUiWaUUfOjsI0VRFHeySin4lIEGnBVFUdzJKqXgq6usieUURVHcyZrcRwCPjOjFi1O/56iyxqkWRVEUJS3JKqXQrKSIW07tkmoxFEVR0pasch8piqIokVGloCiKovhRpaAoiqL4UaWgKIqi+FGloCiKovhRpaAoiqL4UaWgKIqi+FGloCiKovgRk2HZ4URkM/B9NU8vBbYkUJxUos+SftSW5wB9lnSlJs/S1hjTNFqnjFMKNUFEZhpjylMtRyLQZ0k/astzgD5LupKMZ1H3kaIoiuJHlYKiKIriJ9uUwuhUC5BA9FnSj9ryHKDPkq54/ixZFVNQFEVRIpNtloKiKIoSAVUKiqIoip+MVwoiUiQiX4vIPBFZJCJ32+3tRGS6iCwTkf+KSIHdXmjvL7ePlzmudZvdvlRETk3R8+SKyBwReT/Dn2OViCwQkbkiMtNuaywiH9vP8rGINLLbRUSesGWeLyJ9Hde52O6/TEQuTtGzNBSRN0TkGxFZIiKDMu1ZRKSL/b/w/ewUkRsz7TkcMtxkf98XishYexzI1O/KDfZzLBKRG+221P1fjDEZ/QMIUN/ezgemAwOB14Dz7Pangavt7WuAp+3t84D/2tvdgHlAIdAO+A7ITcHz3Ay8Arxv72fqc6wCSoPaHgZG2dujgIfs7eHAB/b/ciAw3W5vDKywfzeytxul4Fn+D7jc3i4AGmbqs9iy5AIbgLaZ+BxAS2AlUMfefw24JBO/K0APYCFQF6sS5idAp1T+X5L+gfT4D1wXmA0MwFr1l2e3DwI+src/AgbZ23l2PwFuA25zXMvfL4nytwI+BU4A3rflyrjnsO+7ilClsBRoYW+3AJba288AI4P7ASOBZxztAf2S9Bwl9gAkmf4sjnufAkzJ1OfAUgpr7AEwz/6unJqJ3xXgHGCMY/9PwK2p/L9kvPsI/C6XucAm4GMsjb/dGFNpd6nA+iBB1QcK+/gOoImz3eWcZPE3rA/EIXu/CZn5HAAGmCAis0TkSrutuTFmPYD9u5ndHk7mdHiW9sBm4HnbrTdGROqRmc/i4zxgrL2dcc9hjFkLPAqsBtZjffZnkZnflYXAsSLSRETqYlkCrUnh/6VWKAVjzEFjzJFYb9r9gSPcutm/JcyxcO1JQUR+BmwyxsxyNrt0TevncHC0MaYvcBrwOxE5NkLfdH6WPKAv8C9jTB9gN5Y5H450fhZsP/sZwOvRurq0pcVz2P71M7FcPocD9bA+Z+HkSttnMcYsAR7Cepn9EMudVRnhFM+fpVYoBR/GmO3A51i+toYikmcfagWss7crsDQx9vEGwFZnu8s5yeBo4AwRWQW8iuVC+huZ9xwAGGPW2b83AW9jKeuNItICwP69ye4eTuZ0eJYKoMIYM93efwNLSWTis4A1eM42xmy09zPxOU4CVhpjNhtjDgBvAYPJ3O/Kc8aYvsaYY225lpHC/0vGKwURaSoiDe3tOlgfmCXAZ8AIu9vFwDv29rv2PvbxicZywr0LnGfPVGiHFez5OjlPAcaY24wxrYwxZVjm/URjzPlk2HMAiEg9ESn2bWP5sBcGyRz8LBfZMysGAjtsk/kj4BQRaWS/HZ5ityUNY8wGYI2IdLGbTgQWk4HPYjOSKtcRZOZzrAYGikhdERGq/icZ910BEJFm9u82wNlY/5/U/V+SGVTxKFDTC5gDzMcaeO6029tj/YOXY5nKhXZ7kb2/3D7e3nGtO7DiEUuB01L4TEOpmn2Ucc9hyzzP/lkE3GG3N8EKpC+zfze22wV40pZ5AVDuuNZv7GdcDlyaov/HkcBM+zP2P6zZHRn3LFgTMX4AGjjaMu45bBnuBr6xv/MvYc0gyrjvii3DZCylNg84MdX/F01zoSiKovjJePeRoiiKkjhUKSiKoih+VCkoiqIoflQpKIqiKH5UKSiKoih+VCkoCUdEHvdle7T3PxKRMY79v4rIzQm+54+JvJ59zSNFZLhj/y4RuSWG80REJopIiaPtLBExItLVAznLROTXib5ulHv2FJEXknlPJTmoUlC84CusFaaISA5QCnR3HB8MTEmBXPFyJFYumngZDswzxux0tI0EvsRamJhoygBXpeBY4ZtQjDELgFb2giulFqFKQfGCKdhKAUsZLAR22astC7FyU80Rkfoi8qmIzBar9sKZACLykIhc47uY/Yb+e3v7DyIyw84lf7fbzd362G/TS0TkWbHy1k+wV8AjIkfZfaeKyCNi5bYvAO4BzhWr/sC59uW7icjnIrJCRK4P8/znU7UCFRGpj5XG5DIcSkFEhtrX8tVqeNleoYuIDLfbvhQrf76vvsZxUlUTYY69cvxBYIjddpOIXCIir4vIe1hJCcXxXAt8z2Lff5KIvCYi34rIgyJyvlj1SRaISAe73zn2ufNE5AvHc76HN0pOSSWpWMGnP7X/Byt1dhvgt8BVwL1Yb9BHA1/YffKAEnu7FGslpgB9gEmOay22r3UKVuFywXqheR841u7zo/3btQ/W23QlcKTd7zXgAnt7ITDY3n4QWGhvXwL80yHHXVhWUKEt7w9Avsuzfw8UO/YvAJ6zt78C+trbQ7EydrayZZ0KHIO1AncN0M7uN5aqFe7vYSUbBKhv/w2H+o475K6gahXsL7ESruUCzbHSRLSwz9tubxcCa4G77XNuAP5mby8AWtrbDR33ORp4L9WfNf1J7I9aCopX+KyFwViD3VTH/ld2HwH+IiLzsYqLtMRKGTwHaCYih4tIb2CbMWY11oB/ClZak9lAV6x8NU4i9VlpjJlrb88CysTKm1VsjPHJ9EqU5xpnjPnJGLMFK0lZc5c+jY0xuxz7I7GSHGL/Huk49rUxpsIYcwiYi6W8ugIrjDEr7T7OXEVTgMdsK6WhqUoVHczHxpit9vYxwFhjZRPeCEwCjrKPzTDGrDfG/ISVOmGC3b7AlsV3zxdE5AosxeJjE1aWUqUW4Ym/UVGoiiv0xHoTXwP8HtgJ/Nvucz7QFOhnjDkgVobYIvvYG1jJyw6jakAV4AFjzDMR7uvaR6wSjD85mg4CdXBPORyJ4Gu4fYcqRSTHGHNIRJpgZbztISIGa1A1InJrhOuFlckY86CIjMOyuqaJyElhuu52bEd6Ruf9Dzn2D9myYIy5SkQGAKcDc0XkSGPMD1j/q70Rrq1kIGopKF4xBfgZsNV+Q92KVcZyEJbVAFYK4022Qjgeqzykj1ex/NUjsBQEWFkff2P76BGRlmJnmHQQSx8/xphtWPGOgXaT00e+CyiO56FtlmIlZ8OW/0VjTFtjTJkxpjVWJbdjIpz/DdBeqmoJ++IZiEgHY8wCY8xDWEn6usYg5xdYsZFcEWmK5U6LORuofc/pxpg7saqW+VI0d8ZS+EotQpWC4hULsPzu04LadtiuF4CXgXIRmYllNXzj62iMWYQ10K01VRWoJmC5d6aKyAIsZREwGMbSx4XLgNEiMhXrrXqH3f4ZVmDZGWiOhXFY/nqwXEVvBx1/kzCzhexn2ItVV/hDEfkS2OiQ6UZf0BfrLf0DrOytlXYg+CaXS75t95kHTARuNVZK8Fh5xA48L8RSMPPs9uPtZ1VqEZolVcl6RKS+MeZHe3sUVm3cG2pwvRZY1sHJNZXJno30JLDMGPN4da+XaOxZZJOAYyLENZQMRC0FRYHTbWtgITAEuK8mF7Mtm2fFsXitGlwhVt3xRVhutkhxlFTQBhilCqH2oZaCoiiK4kctBUVRFMWPKgVFURTFjyoFRVEUxY8qBUVRFMWPKgVFURTFz/8H7XcsRh8oRPwAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "doppler = DopplerShift(vrad=0)\n", "doppler.vrad = 20\n", "new_wave, new_flux = doppler(wave, flux)\n", "plt.plot(new_wave,new_flux)\n", "plt.xlabel('Wavelength (Angstroms)')\n", "plt.ylabel('Flux')" ] }, { "cell_type": "code", "execution_count": 18, "metadata": {}, "outputs": [], "source": [ "from starkit.base.operations.stellar import RotationalBroadening, CCM89Extinction\n", "rotation = RotationalBroadening.from_grid(grid, vrot=200) # it is imperative to load from_grid here as it needs to\n", "extinction = CCM89Extinction(a_v=1.0, r_v=3.1)\n", "my_model = grid | rotation | doppler | extinction\n", "my_model\n", "my_model.param_names\n", "my_model.vrot_1 = 300" ] }, { "cell_type": "code", "execution_count": 20, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "\n", "Name: CompoundModel11\n", "Inputs: ()\n", "Outputs: ('wavelength', 'flux')\n", "Fittable parameters: ('teff_0', 'logg_0', 'mh_0', 'vrot_1', 'limb_darkening_1', 'vrad_2', 'a_v_3', 'r_v_3')\n", "Expression: [0] | [1] | [2] | [3]\n", "Components: \n", " [0]: \n", "\n", " [1]: \n", "\n", " [2]: \n", "\n", " [3]: " ] }, "execution_count": 20, "metadata": {}, "output_type": "execute_result" } ], "source": [ "from starkit import assemble_model\n", "my_model = assemble_model(grid, vrad=20, a_v=0.3, vrot=300)\n", "my_model.__class__" ] }, { "cell_type": "code", "execution_count": 22, "metadata": {}, "outputs": [ { "ename": "AssertionError", "evalue": "", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mAssertionError\u001b[0m Traceback (most recent call last)", "\u001b[0;32m\u001b[0m in \u001b[0;36m\u001b[0;34m\u001b[0m\n\u001b[0;32m----> 1\u001b[0;31m \u001b[0;32massert\u001b[0m \u001b[0;36m1\u001b[0m\u001b[0;34m==\u001b[0m\u001b[0;36m0\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m", "\u001b[0;31mAssertionError\u001b[0m: " ] } ], "source": [ "assert 1==0" ] }, { "cell_type": "code", "execution_count": null, "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.7.2" } }, "nbformat": 4, "nbformat_minor": 2 }