{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Exo-REM\n", "\n", "Example on how to access to the atmospheric grid info to check the parameters and their ranges.\n", "\n", "### Imports" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Open the grid\n", "\n", "You can check the free parameters names, ranges and units together with the wavelength range and spectral resolution." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "Dimensions: (wavelength: 29922, par1: 33, par2: 5, par3: 4, par4: 15)\n", "Coordinates:\n", " * wavelength (wavelength) float64 0.6667 0.6667 0.6667 ... 245.4 248.4 251.6\n", " * par1 (par1) float64 400.0 450.0 500.0 ... 1.9e+03 1.95e+03 2e+03\n", " * par2 (par2) float64 3.0 3.5 4.0 4.5 5.0\n", " * par3 (par3) float64 -0.5 0.0 0.5 1.0\n", " * par4 (par4) float64 0.1 0.15 0.2 0.25 0.3 ... 0.6 0.65 0.7 0.75 0.8\n", "Data variables:\n", " grid (wavelength, par1, par2, par3, par4) float64 ...\n", "Attributes:\n", " key: ['par1', 'par2', 'par3', 'par4']\n", " par: ['teff', 'logg', 'mh', 'co']\n", " title: ['Teff', 'log(g)', '[M/H]', 'C/O']\n", " unit: ['(K)', '(dex)', '', '']\n", " res: [29999.50000004 29998.49999991 29997.50000002 ... 80.5\\n ...\n" ] } ], "source": [ "grid_path = '/EXOREM_native.nc'\n", "ds = xr.open_dataset(grid_path, decode_cf=False)\n", "print(ds)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Get a spectrum for specific values" ] }, { "cell_type": "code", "execution_count": 14, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlAAAAE8CAYAAAAPPI/5AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjkuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy80BEi2AAAACXBIWXMAAA9hAAAPYQGoP6dpAABemElEQVR4nO3dd1gUV9sG8HsBARXBhliw94YEewtgV+y9xBaNvcX2aizYojG22DWaiDUWYjd2QWNJbKjYG3YBFaQJUvZ8f/jthGUX2F12Gcr9u665dM6eOfPMzO7y7MyZMwohhAARERER6cxM7gCIiIiIMhsmUERERER6YgJFREREpCcmUERERER6YgJFREREpCcmUERERER6YgJFREREpCcmUERERER6YgJFREREpCcmUJQlPXr0CC1atICdnR0UCgX2798PALhy5QoaNGiA3LlzQ6FQ4MaNG7LGmZkpFArMmjVL7jCyNTc3N7i5uaXrOtPruHt5eUGhUEjT+/fvTb5OXd24cUMtNm9vb7lDkpQqVQoDBgyQO4xsIdsnUEk/pEmnf/75R9Z4LCwsUKxYMQwYMACvX7/WqO/m5pZs7JUqVdLa7vnz5zXaEUKgePHiUCgUaNu2rdG2Z8eOHfjll19SrTdr1qwUj4Nq0vWPRf/+/eHv748ff/wRW7duRa1atRAXF4du3bohJCQEy5Ytw9atW1GyZEmDt+3Bgwf4/vvv0aBBA1hbW0OhUODZs2ca9Xx9fVPcph9//FGt/sePHzFkyBDY29sjd+7ccHd3x/Xr17XGcPDgQbi4uMDa2holSpSAp6cn4uPj1eqo9m3SP0AvX75E2bJlkT9//mTbz25+/PFHKBQKVKtWTe5Q6P+pPqt58uSRygYMGACFQgFbW1tER0drLPPo0SPp87V48WKt7dasWRMjRoyQ2rOxsUk2BoVCgVGjRknzJUuWxNatW/HDDz8YulkZQkxMDJYtW4a6devCzs4O1tbWqFChAkaNGoWHDx9q1A8NDYWFhQV2794tlcXFxWHFihWoXbs28uTJAxsbG9SuXRsrVqxAXFxcem5OurOQO4CMYs6cOShdurRGebly5WSI5r94YmJi8M8//8DLywvnz5/H7du3YW1trVbX0dERCxYs0GjDzs5Oo8za2ho7duxAo0aN1MrPnj2LV69ewcrKyqjbsWPHDty+fRvjxo1LsV7nzp3V9nVkZCSGDx+OTp06oXPnzlK5g4NDquuMjo7GpUuXMG3aNLUvvfv37+P58+fYsGEDBg8erP/GJHHp0iWsWLECVapUQeXKlZM9m1W5cmVs3bpVo3zr1q04ceIEWrRoIZUplUp4eHjg5s2bmDRpEgoWLIg1a9bAzc0N165dQ/ny5aW6R48eRceOHeHm5oaVK1fC398f8+bNQ3BwMNauXZti7K9fv4a7uztCQkJw6tQpuLi4GLYTspBXr15h/vz5yJ07t9yhUCIdO3ZEqVKlNMotLCzw6dMnHDp0CN27d1d7bfv27bC2tkZMTIzWNt++fQs/Pz/MmTPHoJjy5cuHb775Br6+vpg/f75Bbcjt/fv3aNWqFa5du4a2bduid+/esLGxwYMHD7Bz5078+uuviI2NVVvm+PHjUCgU0ndWVFQUPDw8cPbsWbRt2xYDBgyAmZkZjh07hrFjx2Lv3r04cuRIlv1MMYH6f61bt0atWrXkDkOSOJ7BgwejYMGCWLhwIQ4ePKjxZWFnZ4dvvvlGp3bbtGmDPXv2YMWKFbCw+O/w79ixAzVr1pTtNLmTkxOcnJyk+ffv32P48OFwcnLSedtU3r17BwDImzevWnlwcLDWckO1b98eHz9+RJ48ebB48eJkEygHBwet2zB79myUL18etWvXlsq8vb1x8eJF7NmzB127dgUAdO/eHRUqVICnpyd27Ngh1Z04cSKcnJxw4sQJ6Vja2tpi/vz5GDt2rNoZyMTevHkDd3d3fPjwASdPnkTNmjUN3QVZysSJE1GvXj0kJCRkqMtFWVlUVJTBf1ytrKzQsGFD/PHHHxrfiTt27ICHhwf+/PNPrcsePXoU1tbWaNKkiUHrzgoGDBgAPz8/eHt7o0uXLmqvzZ07F9OmTdNY5q+//kLDhg2l79Dx48fj7NmzWLlypdqP1eHDh2P16tUYNWoUJk6cmOoPuswq21/C05WnpyfMzMxw+vRptfIhQ4bA0tISN2/elMr27NmDmjVrImfOnChYsCC++eYbrZff9NG4cWMAwJMnT9LUTq9evaQ/nCqxsbHw9vZG7969dW7nwIED8PDwQNGiRWFlZYWyZcti7ty5SEhIkOq4ubnhyJEjeP78uXQ6XdsvSX3cv38fXbt2Rf78+WFtbY1atWrh4MGD0uuzZs2SLstNmjRJWueAAQPg6uoKAOjWrZtelwOTkz9/frXLCvq4fPkyHj9+jD59+qiVe3t7w8HBQe2sm729Pbp3744DBw7g8+fPAIC7d+/i7t27GDJkiFoiPGLECAghku2T8fbtW7i7uyM4OBgnTpww+o8GPz8/tG7dGra2trCxsUHTpk21Xga/desWXF1dkTNnTjg6OmLevHnYtGmTxmVQpVKJWbNmoWjRosiVKxfc3d1x9+5do/fzOHfuHLy9vXW63JyakJAQTJw4EdWrV4eNjQ1sbW3RunVrte8I4L9Lu7t378aPP/4IR0dHWFtbo2nTpnj8+LFGu7/++ivKli2LnDlzok6dOvj77791junkyZNo1KgR8ubNCxsbG1SsWFHj8lNwcDAGDRoEBwcHWFtbo0aNGti8eXOqbT9//hwjRoxAxYoVkTNnThQoUADdunXTuJyt6kZw9uxZjBgxAoUKFYKjo6PO26BN7969cfToUXz8+FEqu3LlCh49epTi99mRI0fg7u6OnDlzpmn9utJ3H124cAHjx4+XLuN36tRJ+mGoIoTAvHnz4OjoKH027ty5o1M8//77L44cOYJBgwZpJE/Al+Q06aVPpVKJY8eOwcPDA8CXM7a//fYbmjRpopY8qYwcORLu7u7YuHEjXr16pVNcmQ3PQP2/sLAwjV+dCoUCBQoUAABMnz4dhw4dwqBBg+Dv7488efLg+PHj2LBhA+bOnYsaNWoA+PIBGDhwIGrXro0FCxYgKCgIy5cvx4ULF+Dn52fw2Q/VBy1fvnwaryX3izlnzpwav+5KlSqF+vXr448//kDr1q0BfPk1FhYWhp49e2LFihU6xePl5QUbGxuMHz8eNjY2OHPmDGbOnInw8HAsWrQIADBt2jSEhYXh1atXWLZsGQCk2M8gNXfu3EHDhg1RrFgxTJkyBblz58bu3bvRsWNH/Pnnn9Llvrx58+L7779Hr1690KZNG9jY2MDBwQHFihXD/PnzMWbMGNSuXVu6HPj582dEREToFEPBggUNjj+x7du3A4BGAuXn5wcXFxeYman/tqlTpw5+/fVXPHz4ENWrV4efnx8AaCRARYsWhaOjo/R6YkFBQejatSsCAwNx4sQJtTNfxnDnzh00btwYtra2mDx5MnLkyIH169fDzc0NZ8+eRd26dQH8d/lQoVBg6tSpyJ07NzZu3Kj18vHUqVPx888/o127dmjZsiVu3ryJli1balyaUSqVCAkJ0SlOOzs75MiRQ5pPSEjA6NGjMXjwYFSvXj0Ne+CLp0+fYv/+/ejWrRtKly6NoKAgrF+/Hq6urrh79y6KFi2qVv+nn36CmZkZJk6ciLCwMPz888/o06cP/v33X6nOb7/9hqFDh6JBgwYYN24cnj59ivbt2yN//vwoXrx4ivHcuXMHbdu2hZOTE+bMmQMrKys8fvwYFy5ckOpER0fDzc0Njx8/xqhRo1C6dGns2bMHAwYMwMePHzF27Nhk279y5QouXryInj17wtHREc+ePcPatWvh5uaGu3fvIleuXGr1R4wYAXt7e8ycORNRUVH67FoNnTt3xrBhw7B37158++23AL6cfapUqVKyl6Xj4uJw6tQprZfeTHXmUd99NHr0aOTLlw+enp549uwZfvnlF4waNQq7du2S6sycORPz5s1DmzZt0KZNG1y/fh0tWrTQuOymjepHZ9++ffXahnfv3qFNmzYAvvzdSEhIQL9+/ZJdpl+/fvDx8cGxY8eM0m0iwxHZ3KZNmwQArZOVlZVaXX9/f2FpaSkGDx4sQkNDRbFixUStWrVEXFycEEKI2NhYUahQIVGtWjURHR0tLXf48GEBQMycOVPneE6dOiXevXsnXr58Kby9vYW9vb2wsrISL1++VKvv6uqabPxDhw7VaPfKlSti1apVIk+ePOLTp09CCCG6desm3N3dhRBClCxZUnh4eKQap2rZxIYOHSpy5colYmJipDIPDw9RsmTJVNtL6t27dwKA8PT0lMqaNm0qqlevrta+UqkUDRo0EOXLl5fKAgICBACxaNEitTZ9fHwEALFnzx618pTeA0mn5CxatEgAEAEBAaluW3x8vHBwcBB16tTReC137tzi22+/1Sg/cuSIACCOHTumtr4XL15o1K1du7aoV6+eNO/p6SkAiJIlSwpbW1tx6dKlVGPURdLj07FjR2FpaSmePHkilb1580bkyZNHfP3111LZ6NGjhUKhEH5+flLZhw8fRP78+dX2YWBgoLCwsBAdO3ZUW++sWbMEANG/f3+pTHXMdZl8fHzU2lu1apWws7MTwcHBQogvn6mqVasavF9iYmJEQkKCWllAQICwsrISc+bMkcpU78fKlSuLz58/S+XLly8XAIS/v78Q4r/vFWdnZ7V6v/76qwAgXF1dU4xn2bJlAoB49+5dsnV++eUXAUBs27ZNKouNjRX169cXNjY2Ijw8XCpPety1fRdcunRJABBbtmyRylSfs0aNGon4+PgUY05cX9tnqn///iJ37txCCCG6du0qmjZtKoQQIiEhQRQuXFjMnj072e+B06dPa7Tbv3//VN83I0eO1Igjue+UpPTdR82aNRNKpVIq//7774W5ubn4+PGjEEKI4OBgYWlpKTw8PNTq/fDDDxqfDW06deokAIjQ0NAU6yU2Y8YMte/ycePGCQBqn+Okrl+/LgCI8ePH67yezIRnoP7f6tWrUaFCBbUyc3Nztflq1aph9uzZmDp1Km7duoX379+r9T+5evUqgoODMWvWLLWO3h4eHqhUqRKOHDmC2bNn6xRPs2bN1OZLlSqFbdu2aT3lXapUKWzYsEGjPLnT4927d8e4ceNw+PBhtGrVCocPH9b5zJNK4lPfERER+Pz5Mxo3boz169fj/v370hk5YwkJCcGZM2cwZ84cREREqJ0xatmyJTw9PfH69WsUK1ZM77ZbtmypdknT1E6fPo2goCCtd/BER0drPROjej+p7jhS/Ztc3fDwcI3yoKAg5M+fH0WKFElT/NokJCTgxIkT6NixI8qUKSOVFylSBL1798aGDRsQHh4OW1tbHDt2DPXr14ezs7NUL3/+/OjTpw9WrlwplZ0+fRrx8fHSnVIqo0eP1riNvnDhwjofw8TvzQ8fPmDmzJmYMWMG7O3t9dji5CU+JgkJCfj48aN02Uzb3Y4DBw6EpaWlNK+6XP/06VNUq1ZN+l6ZM2eOWr0BAwZg0qRJqcajOut94MABDBw4UOPsJvClb0vhwoXRq1cvqSxHjhwYM2YMevXqJXUS1ibxd0FcXBzCw8NRrlw55M2bF9evX9c4y/Hdd99pfLemRe/evdGtWzcEBgbi9u3bCAwMTPHy3V9//YUqVapodCewtrbGoUOHtC7TvHnzNMWo7z4aMmQIFAqFNN+4cWMsW7YMz58/h5OTE06dOoXY2FiMHj1ard64ceN06tSu+n7QpwvCX3/9JV2+AyB9B6fUhuo1bd9HWQETqP9Xp04dnfqDTJo0CTt37sTly5cxf/58VKlSRXrt+fPnAICKFStqLFepUiVp+ICEhASN69n58+dX+3JUJXRhYWH4/fffce7cuWTvkMudO7dGwpUSe3t7NGvWDDt27MCnT5+QkJAgdVjW1Z07dzB9+nScOXNG48MRFhamV1u6ePz4MYQQmDFjBmbMmKG1TnBwsEEJVJEiRUySVCRn+/btMDc3R48ePTRey5kzp9TPKTHVJSvVF7Hq3+TqauvbsW3bNnzzzTdo3rw5zp8/j0KFCqVpOxJ79+4dPn36pPW9X7lyZSiVSrx8+RJVq1bF8+fPUb9+fY16Se94VX2ekpbnz59f41K2tbW1Xp8BlenTpyN//vwYPXq03ssmR6lUYvny5VizZg0CAgLU+gWqugQkVqJECbV51baFhoYC+G8/JL4DE/iS4CROVpPTo0cPbNy4EYMHD8aUKVPQtGlTdO7cGV27dpWSqefPn6N8+fIayVXlypXVYtAmOjoaCxYswKZNm/D69WsIIaTXtH0XaLvbOS3atGmDPHnyYNeuXbhx4wZq166NcuXKaR1SBPjS/6ldu3Ya5ebm5ga9h3Sh7z4y9D1hb2+vtZtHUra2tgC+JEG6dCsJDAzE9evX1e5aVCVHKXV/0CXJysyYQOnp6dOnePToEQDA39/foDZevnyp8SXi4+Oj1qk5cULXsWNHNGrUCL1798aDBw/S1I9IpXfv3vjuu+8QGBiI1q1b69U36+PHj3B1dYWtrS3mzJmDsmXLwtraGtevX8f//vc/KJXKNMeXlKrNiRMnomXLllrrGDrkRHR0tM5JX+HChQ1aR+J17du3D82aNdM6JEORIkXw9u1bjXJVmar/jCrhe/v2rUYfmLdv36JOnToabbi6umL37t3o3LkzWrZsCV9fX61DXWRG2n6UJEf1Y+XRo0f49ddf8csvv+DNmzfS6zExMYiLi8OzZ89ga2uL/Pnz6xXL/PnzMWPGDHz77beYO3cu8ufPDzMzM4wbN07rZyO5szGJ/8imRc6cOXHu3Dn4+PjgyJEjOHbsGHbt2oUmTZrgxIkTaT4bNHr0aGzatAnjxo1D/fr1pcFre/bsqXV7jd1x28rKCp07d8bmzZvx9OnTFAf5DAgIwP3799P9rjB995Gp3xOqO3T9/f2lM54pUd216O7uLpWpkutbt26pnU1O7NatWwCgdqIhK2ECpQelUokBAwbA1tZWOlXatWtX6Y4p1d1fDx480Lg99sGDB9Lr2i43pHTJy9zcHAsWLIC7uztWrVqFKVOmpHlbOnXqhKFDh+Kff/5R65ioC19fX3z48AF79+7F119/LZUHBARo1E18ejktVL+0c+TIYfRfibt27cLAgQN1qpvWL7CDBw8iIiJCo/O4irOzM/7++28olUq1swH//vsvcuXKJV1mVn1hXb16VS1ZevPmDV69eoUhQ4Zobb9du3b4/fff0b9/f7Rt2xYnTpwwyh80e3t75MqVCw8ePNB47f79+zAzM5MSvZIlS2q9yyxpmerz8vjxY7UfHB8+fJB+iato+1GSHNWPldevX0OpVGLMmDEYM2aMRr3SpUtj7Nixet+Z5+3tDXd3d/z2229q5R8/fjToJgTVfnj06JHa90pcXBwCAgJ0ulxuZmaGpk2bomnTpli6dCnmz5+PadOmwcfHB82aNUPJkiVx69Ytjffd/fv31WLQxtvbG/3798eSJUukspiYGLU740ytd+/e+P3332FmZoaePXsmW+/IkSOws7PTGAfP1Iy9jxK/JxKfhXz37p3GZ0Obdu3aYcGCBdi2bZtOCZS2uxZbt24Nc3NzbN26NdmO5Fu2bIGFhQVatWqV6joyIyZQeli6dCkuXryIgwcPwsPDA76+vhg+fDi+/vprFCxYELVq1UKhQoWwbt06fPvtt9Ilt6NHj+LevXuYOXMmAMMuN7i5uaFOnTr45ZdfMG7cOI3BNPVlY2ODtWvX4tmzZ1pPZ6dE9esocTIRGxuLNWvWaNTNnTu3US7pFSpUCG5ubli/fj1Gjx6tccnt3bt3BvdhSc8+UDt27ECuXLnQqVMnra937doV3t7e2Lt3r3RZ9f3799izZw/atWsnvaeqVq2KSpUq4ddff8XQoUOlY7J27VooFIoUL8n27dsXoaGhGDt2LLp06YIDBw6o3ZVmCHNzc7Ro0QIHDhzAs2fPpP4lQUFB0sCtqssGLVu2xOrVq3Hjxg0pEQwJCZHuTFRp2rQpLCwssHbtWrU+KKtWrdJYvyF9oKpVq4Z9+/ZpvD59+nRERERg+fLlKFu2rE5tJmZubq6RaO/ZswevX7826CxprVq1YG9vj3Xr1qn1l/Ly8tLpD3BISIjGWTTVflddAm7Tpg1OnDiBXbt2Sf2g4uPjsXLlStjY2EhDgGijbXtXrlypdunS1Nzd3TF37lwUKFAgxbPEf/31F1q0aKE29Ed6MPY+atasGXLkyIGVK1eiRYsW0g9VXZP9+vXro1WrVti4cSNat26Njh07qr0eGxuLH374AYsXL0ZcXBxOnjypMVhz8eLFMXDgQGzcuBFr167F8OHD1V5ft24dzpw5g6FDh6Z5uIqMignU/zt69Kj0ayuxBg0aoEyZMrh37x5mzJiBAQMGSAmHl5cXnJ2dMWLECOzevRs5cuTAwoULMXDgQLi6uqJXr17SMAalSpXC999/n6YYJ02ahG7dusHLywvDhg2TysPCwrBt2zaty6Q0CGX//v0NiqNBgwbIly8f+vfvjzFjxkChUGDr1q1az87UrFkTu3btwvjx41G7dm3Y2NjonbCprF69Go0aNUL16tXx3XffoUyZMggKCsKlS5fw6tUrjXF2dGVoH6iwsDCp07PqlvBVq1Yhb968yJs3r8bYKCEhITh69Ci6dOmS7GXYrl27ol69ehg4cCDu3r0rjUSekJCgcQPCokWL0L59e7Ro0QI9e/bE7du3sWrVKgwePFg6vZ6cMWPGICQkBLNnz0a/fv2wfft2mJmZwdfXF+7u7vD09NT7eWfz5s2TxhsaMWIELCwssH79enz+/Bk///yzVG/y5MnYtm0bmjdvjtGjR0vDGJQoUQIhISHSHwMHBweMHTsWS5YsQfv27dGqVSvcvHkTR48eRcGCBdXObhryo6RgwYIafziA//4IJX1t1qxZmD17tsbl9qTatm2LOXPmYODAgWjQoAH8/f2xfft2nforaZMjRw7MmzcPQ4cORZMmTdCjRw8EBARg06ZNOrU5Z84cnDt3Dh4eHihZsiSCg4OxZs0aODo6SmdihgwZgvXr12PAgAG4du0aSpUqBW9vb1y4cAG//PJLin1Y2rZti61bt8LOzg5VqlTBpUuXcOrUKa39vUzFzMwM06dPT7FOdHQ0fHx8sG7dunSK6j/G3kf29vaYOHEiFixYgLZt26JNmzbw8/OTPhu62LJlC1q0aIHOnTujXbt2aNq0KXLnzo1Hjx5h586dePv2LRYvXozz588jPDxcrQO5yrJly3D//n2MGDECx44dk840HT9+HAcOHICrq6vaWbcsR6a7/zKM1G5h37Rpk4iPjxe1a9cWjo6O0m2kKqpbjnft2iWV7dq1S3z11VfCyspK5M+fX/Tp00e8evVKr3iuXLmi8VpCQoIoW7asKFu2rHQbcErDGCQ+vCm1m5iuwxhcuHBB1KtXT+TMmVMULVpUTJ48WRw/flzjNvHIyEjRu3dvkTdvXulWel1oG8ZACCGePHki+vXrJwoXLixy5MghihUrJtq2bSu8vb2lOvoOY2ColG6d17ad69atEwDEwYMHU2w3JCREDBo0SBQoUEDkypVLuLq6Jnvc9u3bJ5ydnYWVlZVwdHQU06dPF7GxsWp1VMMYaLuNffTo0QKAGDZsmBBCiEOHDgkAYt26daluv7bjc/36ddGyZUthY2MjcuXKJdzd3cXFixc1lvXz8xONGzeW4l6wYIFYsWKFACACAwOlevHx8WLGjBmicOHCImfOnKJJkybi3r17okCBAlLMxpbcMAYTJkwQCoVC3Lt3L8XlY2JixIQJE0SRIkVEzpw5RcOGDcWlS5eEq6ur2pADyb0fVe+rTZs2qZWvWbNGlC5dWlhZWYlatWqJc+fOabSpzenTp0WHDh1E0aJFhaWlpShatKjo1auXePjwoVq9oKAgMXDgQFGwYEFhaWkpqlevrhGDEJrHPTQ0VFrOxsZGtGzZUty/f1+ULFlS7XZ6Xb+DktZPbRiD5CT9Hjh8+LBQKBQiKChI7/aQxmEM0rqPVOtJ/N2akJAgZs+eLb3P3NzcxO3btzXaTMmnT5/E4sWLRe3atYWNjY2wtLQU5cuXF6NHjxaPHz8WQggxceJEUaVKlWTb+Pz5s1i2bJmoWbOmyJ07t8iVK5dwcXERv/zyi8Z3UVaT7RMoIvrPpEmThKOjo9pYW+ll7NixwtraOtUxgkJDQwUAMW/evHSK7IvatWuLrl27pus6szNVMnH9+nXx7t07tfGODDF8+HBRu3Zto8QWHx8v3r17J/bv32/UH2UZUeXKlcWkSZPkDiND4iU8IpL4+PhgxowZRn+odFLR0dFqHVI/fPiArVu3olGjRmp3ICWtB/x3iS2tj+LRR3h4OG7evKnTo03IuFQjir979y5NTwJwdnY2uPtAUv7+/vjqq6+M0lZGFhsbix49emg8a5C+UAhhpPsiiYh05OzsDDc3N1SuXBlBQUH47bff8ObNG5w+fVrtzk4vLy94eXlJj+Q5f/48/vjjD7Ro0QLHjx+XcQvI1N6+fav2bDdXV9c03+xgLJGRkWrPeHRycjLquGqUOTCBIqJ098MPP8Db2xuvXr2CQqGAi4sLPD09NTqCX79+HZMnT8aNGzcQHh4OBwcHdOnSBfPmzTPKeGhERIZiAkVERESkJ82HIhERERFRiphAEREREelJlrvwlEol3rx5gzx58hjtUR9ERERE2gghEBERgaJFi2o8NNtQsiRQb9680XgAKhEREZEpvXz50miPlpElgVI9FuDly5fS87GIiIiITCE8PBzFixdP8bFE+pIlgVJdtrO1tWUCRUREROnCmN2G2ImciIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xAQqGzp+/DgUCgWOHj0qdyhERESZkkIIIdJ7peHh4bCzs0NYWBjHgZJB4nEwZDj8RERE6coUeQfPQBERERHpiQlUNrF//354enryjBMREZERyPIoF0p/nTp1AgDUr19f5kiIiIgyP56BymZat24tdwhERESZHhOobO7MmTNQKpVyh0FERJSpMIHK5po2bYpNmzbJHQYREVGmwgSK4O3tLXcIREREmQoTKCIiIiI9MYEiIiIi0hMTKMKLFy/kDoGIiChTYQJFuHv3rtwhEBERZSpMoIiIiIj0xASKiIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xASKiIiISE9MoIiIiIj0xAQqAzl//jzOnz8vdxhERESUCgu5A6AvPn36hMaNGwMAoqKikCtXLpkjIiIiouTwDFQGERERIf0/KipKxkiIiIgoNUygiIiIiPTEBIqIiIhIT0ygsqCAgAC0adMGPj4+codCRESUJbETeQakVCrTtPw333yDixcv4ujRoxBCGCkqIiIiUuEZqAyoTZs2Bi97+/ZtXLx4Ua0srQkZERERqWMClQFdv37d4GWrV6+uUbZgwYK0hENERERJMIHK4g4dOoTp06fLHQYREVGWwgQqi2vfvr3cIRAREWU5TKCIiIiI9MQEKoNQKBRyh0BEREQ6YgJFREREpCcmUERERER6YgJFREREpCcmUAQAuH//vtwhEBERZRpMoLIApVKJs2fPIjQ01OA2Dh48aMSIiIiIsjYmUFnA5s2b4ebmBhcXF7lDISIiyhaYQGUBu3btAgA8e/ZM3kCIiIiyCSZQGURERISs6+c4VERERLpjApVBPHr0yOBlhRBGjISIiIhSwwQqg5D7DNDDhw9lXT8REVFmwgQqCzDGGaiNGzcaIRIiIqLsgQlUFsBLeEREROmLCRQRERGRnphAZQE8A0VERJS+mEBlEGnpRM4EioiIKH0xgcoCmEARERGlLyZQGQSTICIiosyDCVQWEBUVJXcIRERE2QoTqCzg8uXLcodARESUrTCBygDi4uKwZMkSucMgIiIiHTGBygDWrl2LkydPyh0GERER6YgJVAZw6tQpg5d9+fKlESMhIiIiXTCBktm0adNw6NAhg5Z98OABSpQoYeSIiIiIKDVMoGQ2f/58g5c9cuSIESMhIiIiXTGByuASEhLg7+/PcaL+X3x8PGJjY+UOg4iIsjkmUBmUKmEaNGgQnJyc0nSmKiv4559/sGTJElSsWBGFChViEkVERLJSCBlObYSHh8POzg5hYWGwtbVN79VnKMk9A0+pVEKhUEivW1lZISYmRq3O0qVLMWHCBKPFkpHPciXdT7du3UL16tVlioaIiDITU+QdPAOVib19+1buEIiIiLIlJlCZ2OLFi+UOQTa7d++WOwQiIsrGmEBlUAkJCWrznz9/hr+/P96/fy9TRBnLvHnzEBoaKncYRESUTTGByqDKli2r0efJyckJ9vb2MkWU8YSHh8sdAhERZVNMoDKoFy9e4Pz583KHkaFl5E7vRESUtTGBogxv27ZtWstdXFzw5s2bdI6GiIiICVSGltxYR5cvX8bTp0/TORr59O3bV2t5aGgoJk2alM7REBERARZyB0DJ+/z5s9byunXrpnMkGVdYWJjcIRARUTbEM1CUqR05cgQPHz6UOwwiIspmmEBRplexYkW5QyAiomyGCRRlaFOnTpU7BCIiIg1MoDKwq1evyh2CrIQQ+Omnn+QOg4iISAMTqAzMx8dH7hBkldxdiERERHLjXXiUIYWEhPDsExERZVhMoChDKlCggNwhEBERJYuX8DIwhUIhdwiyuHjxotwhEBERpYgJVAb24cMHuUOQxdChQ+UOgYiIKEVMoGR08uTJFF9/8OBBOkWSsdy+fVvuEIiIiFLEBEpGLVq0kDsEIiIiMgATKCIiIiI9MYEiIiIi0hMTKMpQ4uPj5Q6BiIgoVUygKEPJkSOH3CEQERGligkUERERkZ6YQFGW8OrVK3h7eyMhIUHuUIiIKBtgAkUZRlRUlMHLFi9eHN26dcNvv/1mxIiIiIi0YwJFGUbt2rXT3EZqg5MSEREZAxMomURGRsodQoZz7949uUMgIiLSCRMomfj4+MgdAhERERnIoAQqJiYm2dfevn1rcDBEacVO5ERElB4MSqBcXFxw48YNjfI///wTTk5OaY0pW2jfvr3cIWRJ+/btkzsEIiLKBgxKoNzc3FCvXj0sXLgQwJe7pwYMGIC+ffvihx9+MGqARERERBmNhSELrVmzBh4eHhg8eDAOHz6Mt2/fwsbGBpcvX0a1atWMHSMRERFRhmJQAgUArVu3RufOnbF27VpYWFjg0KFDTJ50xOe9aXr16pXcIRAREenMoEt4T548Qf369XH48GEcP34ckydPRvv27TF58mTExcUZO8YsZ/369XKHkOGMHz/eaG1xOAQiIjI1gxIoZ2dnlC5dGjdv3kTz5s0xb948+Pj4YO/evahTp46xY8xybt++LXcIGcrbt2+xZ88eo7VXpUoVo7VFRESkjUEJ1Jo1a7Bz507kzZtXKmvQoAH8/Pzg4uJirNgomyhatKjcIRAREenFoASqb9++Wsvz5MnDZ5ERERFRlmdQJ/ItW7Yk+5pCoUg2wSJKL1euXDHKs/WIiIi0MSiBGjt2rNp8XFwcPn36BEtLS+TKlYsJVCqEEHKHoFVoaCjy5cuXrusMCAgwSbt16tTB5cuXmUQREZFJGHQJLzQ0VG2KjIzEgwcP0KhRI/zxxx/GjpHSycOHD9N9nWXKlDFZ23zeIBERmYrRHiZcvnx5/PTTTxpnp0iTQqGQO4RsISwsTO4QiIgoizJaAgUAFhYWePPmjTGbzJIy6iW8rGb+/Plyh0BERFmUQX2gDh48qDYvhMDbt2+xatUqNGzY0CiBUdZ36dIluUMgIiIyiEEJVMeOHdXmFQoF7O3t0aRJEyxZssQYcWVpGfXSUnpfWjx37ly6ro+IiMhYDEqglEqlsePIVnbu3Cl3CERERJQGRu0DRaSP69evyx0CERGRQXQ+A6XPw16XLl1qUDCUvezevVvuEIiIiAyicwLl5+enUz3eok9ERERZnc4JlI+PD54+fYpSpUrBzIxX/oiIiCj70isTKl++PN6/fy/N9+jRA0FBQUYPioiIiCgj0yuBSjoA5F9//YWoqCijBkRZ39mzZ9G7d2+5wyAiIjKYQcMYEKWFm5ub3CEQERGliV5noBQKhUYncXYazzp4LDMPIQRCQkIQHh4udyhERNmSXmeghBAYMGAArKysAAAxMTEYNmwYcufOrVZv7969xoswi8nIYx/t27cPtWvXljsMo6pQoQJOnz6N4sWLyx2KURw6dAjnzp1DYGAgtm3bBgDYvHkz+vXrB+DLZ2/JkiXYsWMHSpYsKWeoRERZmkLo8WTbgQMH6lRv06ZNKb4eHh4OOzs7hIWFwdbWVtfVZ3pXrlxBnTp15A4jRenxoGNDz3TNmjULs2bN0nu5nj174o8//jBonRlNcvvOy8sLHTp0QL58+QAALVu2xLFjx9IzNCKiDMsUeYdeCZSxZNcEKjNcIsuoCVTu3LkRGRlp0LKdO3fGn3/+qfdyGYW/vz+cnJz0WqZEiRJ4/vy5iSIiIspcTJF3cEAnyvAsLS3x8uVLucNId35+fqhWrZreyRMAvHjxwgQRERGRCu/Co3RlyGN+WrRoIV2aMkRmfPj1hw8f0K5dO7x+/VruUIiISAuegUonHz58kDuEDGHChAl6L9OoUaM0rfPcuXNpWj69zZ8/HwULFmTyRESUgTGBSie1atWSOwQ1efLkkTuEFLVu3Vr6f1ofHRQdHZ3WcNLVtGnT5A6BiIhSwQQqnTx79kzW9dvY2Ej/79ChAz5+/Ki13urVq02y/jVr1mDIkCE61//mm2+Mtm4Z7pNQc+vWLfTv3z/F90BsbCxmz56NixcvGm298fHxRmuLiIjUMYHKJu7fv682b2ZmhoULF2rUGzVqFMaMGQMASEhIwNKlS3HlypU0r3/kyJHYsGGDzvV79eolnSVr2bJlmtYtdwLl4uKCLVu2oFOnTsnWWb9+PWbNmoWGDRsabb05cuTA8uXLjdYeERH9hwlUNuDh4YFixYpplE+ePFlr/ZUrVwIAtmzZggkTJsgydpVCocDr16/x8OFDg+5CS0zuBCohIQEAcPv27WTrmGqA1XHjxsm+/UREWRETqGwg6R9QXf+g+vv7myKcFHXr1g3BwcEAvvTTKl++fJrbjI2NxdixYzNsIhEXFwcvLy+TtZ/WPmRERKSJ36xGduPGDSgUCixbtkzuUCRyJw763IG4e/du2NvbGz2GFStW4OzZs0Zv11AfPnxA9erVsXDhwnQZTHbhwoUoXbq0XpdRiYgoeUygjOyrr74CAIwfP14qS67DdnoxNIEyVuLVoUMHo7STVjt27EC/fv0QERGRbut89OiR9P/Eo6gvXLgQt2/fxpQpUxATE2PyOKZMmYJnz57p1ZGfiIiSxwQqHTRp0kTW9Rt6Cc9YLly4YLS2Xrx4gX379hm07IYNG7B161b8+OOPRosnNZUrV9Za/v79+3SLIansOKo7EZGxMYEyoW3btiEsLAx+fn6yxiHXGagXL14YPVkrXrw4OnbsmKY2Hj9+bJxgdKDqQJ5Uag/cNqUHDx7Itm4ioqyCCZQJ9e3bN8Vb1zO6tCQ/GzZsQMmSJTFo0CAjRmQc+pz9EUIgJCTEhNGkv/S4ZEhElNUxgTIxHx8fuUNIdwkJCVJfGznPtCQncV+k1IwZMwYFChTAoUOHTBhR+mrXrh0A+W8uICLKzJhAZQO5cuUCAAwYMAAAMHXqVJ2WM/QP7Pnz5w1a7sSJEwYtpy99butftWoVAN33GfBl2ISuXbtqnH2Li4vDoUOH0L59e53bMpW4uDjY2dnB3NyciRQRkQEs5A6ATC937twAgN9//x3Lly83+W3zhj5CpHnz5kaORDt9zkCp6JNkbNmyBX/++afW1zJC8gQAVapUke5GPHjwYIa5U5KIKLPgGahsRKFQ6JU8JU4a6tSpg6tXr+q8HlO7efOmwcsaEp9SqdS5rtzDVugicUf6jh074pdffpEvGCKiTIgJFGkVExOjlkBduXIF7u7uKS6TkJCAq1evpstDbNPyeBddE6hZs2ZJ/0/pDJQQItNfBvv+++/lDoGIKFNhAkVaqfr+JBYZGYnw8HC0atVK6hz+5MkTXL16Fffu3cOUKVNQu3ZtjBw5Uu/1pcdo3Crnz5/HkiVL0KRJE0RGRiZbb/bs2dL/hRD4999/0bRpU9y6dQuRkZFYvnw5nj9/Dnd3dzRo0EA6S5XZkykiIkod+0Clwfv37xEREYHSpUsDALZu3SpzRNoZ8gc9MDBQ63I///wzjh8/juPHj8PR0REtWrTQqGPIOEunT5/WexlDRUdHY+LEiQCA5cuXY9q0aakuI4RAvXr1AAA1atSQyseNGyf9//Xr17hy5UqyD2kmIqKsgwlUGqie2fb27VsULlwY/fr1kzki41myZInW8sT9ezZu3Gi09ZUpU8Zobenj6dOnGmWvX7/WGK1clz5QCoUCXbp0MVps6e3ChQto2LCh3GEQEWUKvIRnBP7+/nKHkG4S9x86evSo0drVp5O2Mf3+++8ayVKHDh2wdu1atbInT56k2lZ6dJ43pUaNGqF79+747bff5A6FiCjDYwJlBJm5z8vw4cP1qp84STDmQ3nlSqAAYPr06Wrz165dM6id9Hhkj5ubm0nb37NnDwYPHmzSdRARZQVMoEgvpjrLYm5ubpJ29XXnzh2Dl1WN8G0qgwYNgo+PD/r27WvS9RARUeqYQJFeTJVAFShQwCTt6iM2NhYuLi5yh4GmTZuibt26GuWq/lWLFi1CkyZNsGvXLgghUKRIkfQOkYgo22MClQXpkwToe/nRFAlUhQoVDFquf//+Ro1j8eLFiI2NNWqbhjh16pTafvbz88PmzZvRqlUrAICDgwNOnz6N7t27A0jboKLJiY+PR0JCgtHbJSLKKphAZTHjxo2Dg4ODWpkx+2jp8xw5XRmaCK1evdpoMXz8+BFHjhwxWnuGWrRoEQCgV69eAIBKlSrB2dkZ/fr1SzZ5tbe3R2BgoFHjyJEjBywsLPDs2TOjtktElFUwgTKCjNSJPOkDbFOjb+dtU5yBMnTcJNUz/owhX758uHjxotHa01eHDh0QGRkpjU81cuRInDhxQueYkibNxpL0DkUiIvqCCVQWU65cOY0kJ6Wk58CBA6YOKVUWFhyObP/+/WoJobm5OZo3b458+fLJGJVxx/oiIspKmEBlMdbW1hplKZ0hCwoK0qv9tNylRqZ15coVo56VU4mOjjZ6m0REmR0TKAPdunVL7hBkcezYMaO2xz/OxlOrVi107tzZ6O3mypUL+fLlw+fPn43eNhFRZsUEykDHjx+XO4RMb/jw4VrPmGU35cuXN1pbpuqP9/HjR1hbW2PhwoUmaZ+IKLNh5xOSjZOTk9whZAjGfBSQqW9omDJlCmrUqIHHjx/D0dERLi4uKFGihEnXSUSUETGB0sOaNWtQpkwZaTyejGb27NkATHdHVkY0YsQIrFmzRu4wDPbdd9/BysrKaO0VLFjQaG0lp3Xr1mrzGekuVCKi9MJLeDq6cuUKRo4cqfHHA0CGGHwRAFxdXQFAr8ssXbt2NVU4qTLGmFL58+c3QiTyUQ2GaSzZKXmmzCU+Pl7uEIiMigmUjl69eqU2v2nTJun/69evT+9wtMqbNy+ALwMr6sqY/W/0Vb9+fdnWnVXxbBBlRLt27ULu3Lmxf/9+AMDVq1dRoUIFaV4bIQQePXoEIQT8/Pxw4sQJvH37Fvfu3VOrExMTgylTpmTbG3tIPkygdJR0LKXEH+KXL1+mdzha1ahRQ+9lTPVsO11w/Cfj7//MmED5+PhwxPMM5t27d5g/fz7evHmj8zIRERHYv38/oqOjERgYqDZIb8+ePREbG4tOnToB+PJcx0ePHqFTp04ICQnRaCsqKgr9+vVDhQoVYGZmBhcXF7Rs2RJFixZFlSpV8Pr1a4SGhqJ06dLImTMnFi5cqPb9t2HDBvj4+KRhDxCljgmUEWT0P1oZPb7s5u7du9L/nZ2djdq2HMf6wIEDuHTpEoQQ8PX1xYMHDwAAmzdvhrOzM54/fy7VfffuHc6ePYvPnz/j9OnTWLt2LZo0aYLSpUune9zZzd27d/Hvv//qVLdDhw6YNm0aihUrhj/++EPj9W3btuH777/Hw4cPpbKePXuiU6dOqFSpEooUKYI+ffpobbtu3bp48eKFNF+gQAFMmzYNFhYWaNKkCaZNmwYbGxts27YtxW1Zv3692nsLAEaNGoU1a9ZgyJAhaNKkCRQKBRQKBZydnRETEwPgy2fkxIkTGlcViPQmZBAWFiYAiLCwMDlWb5BGjRoJAEK1y1T/ByCqVaumUSbHlFji8p49eya7XT/88INs8d67dy/Nx+XHH3+Ufb/rMzk5OQkhhAgPDxeBgYFp3v6k5syZI9u2de7cWfr/2bNnpf+3bdtWCCHEtGnTUn3/KpVK4efnJ2JjY42+b+T2+fNnoVQqRVRUlF7LRUVFiQsXLoiEhASD1qtap2o/BwYGih9++EF07txZJCQkiIiICKnurl27RMuWLTWOzadPn9TaTPxaq1atNMpUU/PmzcWFCxdk/9yppkmTJoljx45pvOcmTpwovLy8DNq/lDmYIu9gAqWjpF/0ieerVasmQkJCZP1iqFChQrLx9urVK9ntSu2PmiknYyRQ4eHhol69emLRokWyfznrMsXHx6d5m1Mye/ZsaV1lypSRfXtVk4eHh071evXqJf3/p59+SnFblUql1vJnz56JDRs2iM+fP5tiFxtk6dKlAoDIkSOHACBev36t87JNmjQRAMTSpUsNWnerVq3U9vGlS5ek/+fKlUv6v7m5ebLH5d27d8LOzk4AECdPntR4fcyYMbK/xwyZvLy8xOnTp6X5zp07i/Xr16vtv8ePH4uNGzeKuLg4g/Y/ZQxMoGSU+EOXdL5KlSri22+/lfWLwNAEavr06bLFbIwEKrltzohTUFCQUbdXm1mzZknr+/Dhg/T/+fPnCx8fH7V4fv/9d9n3SWqTEEJERkaKzZs3iw8fPojPnz+LrVu3ihs3boiSJUsKT09PjX2QM2dOAUDMmTNHCCFEbGys7H/8km7XokWLhBBCfP3116JUqVIa9a9evSoCAgLUli1btqyIj4/XODt34MAB0aZNG7F27VqhVCrFqVOnxFdffSX8/PzE1atXjXIcypYtK/t7wVRT48aNNcrevn0rKlSoIObOnSuV2draipCQEJO/V8g0mEDJKOmXeuL58uXLCxsbG1m/BMqXL59svBk1gbp//75Rj5Gx/liYYoqMjDTqtiYnKChI5M2bVwwePFgIIYS/v79YvHix+Pz5s4iOjhb29vaibt26Ijo6WgiR8ZPO+Ph40bt37xTrREREiBkzZggPDw9RrVo1qbxKlSoiPj5emk98GSokJCTZM1jJiYuLEwsXLhSXL1/WeC0iIkJ8/PhR7NixQ0RERAhfX18xffp0KdlJGvOAAQNEYGCgNL98+XIxd+5ccfToUXHnzh0BQNjY2IiIiAit2+zn5yeEEOLVq1dq5c7OzrIfs6wwjRgxQmt57969pWMeEBAgDh8+rNd7iOTDBEpGiT9ESeczwq+zzJhA6dsXRBePHj2S/Vhom9JTSpcJY2Nj1frSyL1fTD35+fmpza9bt05UqFBBABB9+vTRa7+uWbNG43iqkp3EU58+faT/L168WCQkJGiNLblkx9LSUqdtywjfO9ltKl26tMZn58CBA2LTpk3i+fPner2fKH0xgZJR4g+RUqnU+FDJ/cHObJfw3r17Z5LjFBQUJPux0DZlVKr47O3txcOHD2XfT+k9PXjwQCiVSjFr1iwxduxYsWfPHmnfHD9+XIwYMUJERUWJ2NhY0bp1a2m5+Ph48enTJ1GvXj3Zt4FT+k2JL7cmfc3GxkaOjzDpyBR5BwfiMUDiW3AzO7nGgUqPR45kFL/99pvcIaRq0aJFKF++PA4ePIj27dvLHU66qVixotr88uXLNepoe1QQxzDLnlIarywyMhIfP36EjY0NFAoFzM3N0y8wkgXHgTJA4gHitM1nJnIkUM2bNzdZ23IODJqUn58foqOj8e2338odSrJ27tyJoUOHSmP2tGvXDrt27ZI5KqL/eHp66lx32bJlJozkCz8/Pzx58kTray4uLihfvjycnJw4/l42wATKCDL6ByUjJRXNmzfHxo0bTdZ+jhw5TNa2rhYsWIB///0Xzs7OsLa2ljucFPXo0QPr1q1TO6Ni7OfzEenr5s2b2Lp1K3r16oWZM2dKn6dvvvkGlpaWAIAqVapI9WfOnAlvb2+MGjUKZcqUAQDpX+DLo2QSD/qZFi4uLihXrpzW1wICAvDs2TPcvXsXnz9/Nsr6KAMz2sVAPWT2PlBPnz5Vm3d0dJT92nxKfaAS3zmS1MyZM9MtxitXrpj6MAkhRKp3bZly8vX1TZdtNLXFixeLUaNGZdhO+Zwy93T16lWpE/zGjRs1Xk/Jp0+fxNu3b4UQXzrxr1u3Tu3GiRcvXohly5aJsLAwUa5cOVG0aFHpjsugoCC1mwFMOSUeoJTkx07kMkr64U88X6xYMdm/kAYMGJBsvCklUJ6enukW44sXL0x9mIQQQty9ezdd9/24ceNEuXLlMtTAjcYi9wCxnDLf9PPPP6sN0Lt582ZRr149ce3aNalMCCFCQ0OFj4+PUCqV4vbt2+K7774TdnZ2Yvjw4SZ/X9+5c0ccPnzYpO9v1ZhRSqVS3Lp1K0uOrp+ZMIGSUeIPhpOTk9p83rx5Zf3CmjNnjvj48WOy8Wa3BOrNmzfpuv9NMRxDRqFUKnUeRZyTaaZjx46JuLg4UaJEiWTrNGjQQK82y5cvr7WNqKgoUbVqVY3XEg/NAED873//k/4fHR0tihQpIoAvI84LIcSHDx9E5cqVxaxZs9TeT3fu3BFPnz5N8f2WnpLeUQ1ADBo0yCjHrWHDhsLa2lqt7NWrV9K6z5w5k25n5YkJlKwSfwhS+iKTY0ot3oySQKkGb0wPmzdvTlOsf//9t1i4cGGq9e7evZtu2ySn+/fvy/4+zyrT7NmzxfDhw3Wq26ZNG+kYREVFiUGDBonOnTuLPXv2SJ/dKlWqqI3tpUqOnJ2dxZEjRzTa/N///ieio6NF9erVpbIuXbqoJS/x8fHiw4cPonTp0sLV1VUI8d84WFu3bhVhYWGiT58+4tChQ0IIIZ4+fSpmzpxpsuFJTCkhIUHcuXNH+Pr6infv3onY2FiTHv/t27eLy5cvS/OUPphAyUjuL92UptTiTSmB2rFjh8njCw4OluWLNfGvZF0nb29v6X25du1arXXWrVsnzp8/nyUv2aUmaf+/tE6VKlWS/fNjykn1CJmPHz+Kjx8/qn3nvX//Xly9elUIIcS8efPE1KlTRUJCgnjz5o2YOXOmWLNmTYrHQnVpKOnAqQEBAWLq1KnS8/Z2794tChYsKBYvXizCw8OleiEhIeLQoUMpXlqKj49XS6zS+wyRXHbt2pVu75Hssk/lxgRKRnJ/Eac0pRZvSgnUX3/9ZdLYtm3bZqpDkqrkHoOR0nTq1Clp+XXr1knlN27cEMCXh41md4n31+bNm8X69ev12sdXrlwRo0ePFtu2bRNBQUFi06ZN0kN2M8tkZWUlgoODRadOnUSNGjWkP7j/+9//xNChQwUAMWHCBLkPFaVBer2XYmNj9Xq4NBmGCZSM5P7CTm4qUaJEqvGm9MiKo0ePmjQ+uWm7wyel6dy5c9KyERERoly5cmLkyJHSPH8tfum7UatWLXHt2jUhhFB73lxyU9++fcXDhw9FcHCw1jaTe9xJ4snS0lIsWrRIKJVKMWPGDFGgQIFk6zZt2lTExcUZ/L719/cXmzZtEmfOnBEAxJYtW6TXVqxYobXfm+q9ERsbKy5fvpziI3Uo49u5c6dwdHQUW7ZsEcWLFzfZd6Sq7aFDhwohvjyOKr36i2YnTKBklJ5JkT7T+PHjU403pQTKkLM0uk5bt2411eHQWVRUlHBychLjx48XgwcPVosv8Rkm1ZS4L4kQ2eeSRVoolUphYWGhsS9r1aolOnTooHM7Kb2XSpcuLV0OS26ZxB129+7dK4QQYtiwYQa9d7W5deuW2L9/v0H7iDK/xO+Pffv2mfy7vWbNmmLkyJGiePHiYuvWrfwuSiMmUDLR5Re2MSYXFxe9l0nuMkHiOildwkta15hTRpQ0vsTzfLK64T59+iRCQ0OlfXnv3j2923j69Km4efOm8PLyEsOGDRMLFy4UXbt2FdHR0cn+8VCtz8rKSsTExIgjR46I69evS69HR0eLffv2iYkTJ2q8P1evXi1Wr14tzTdo0EBMnTpVnD592uD9QFnXkCFDBAAxevRoIcSXHw5LlixJl78NqunYsWMy74XMiwmUTA4dOpQuHw59LzcBup2BYgL1n8qVKwsAolWrVkKI/7adp8yNY9u2beLXX39Nt/UtW7ZMp++S2NhYsXfvXvHu3TuxcOFCtb5uN27cEIMGDVK7xZwoqdjYWHHhwgWNTvdxcXFi8uTJAoD4/fffRZUqVaTvlVmzZhn9e9XFxUXkyZNHbNq0SZ4dkUkxgZLJ7t27M0wC9fz5c70TqJQu4SWta6ypdOnSpjgUafb69Wvx888/i/fv3wshvoxX8+TJE5mjIqLMLjIyUgghRGBgoNi8ebOIjo4WkZGRol27dmLJkiXi06dPRv+etba2FrNnzxZPnjwRnz59Eu/evRMHDx7Uerk7uzNF3sFHiusgIz0suESJEnKHkCpfX184OTnJHYZWRYsWxaRJk6T5/PnzI3/+/DJGRERZQe7cuQEADg4O6Nevn1R+8OBBjbqenp6oWLEievfunaZ1xsTEwNPTU+OBy5UrV8aYMWNga2uLZcuW4bfffsuw38mZGRMoHaRXAtWwYcN0WY8pLF68GBMnTgQAuLq6yhwNEVHG888//+DUqVP43//+B3Nzc0RERODAgQNQKpWYOnUqunTpgiJFisDf3z9N67l37x6GDx8uzdeoUQNmZmZo3bo1wsPD4eDggP79+8PZ2RmOjo5QKpVQKpWIj4/X+wHoCQkJuHTpEpydnWFjY5OmuDMbJlA6SK8EqlKlSrh58yYKFy4MBwcHnZYRQpg4qtQtWrQIEyZMQFxcHIoXLy53OEREGVLdunVRt25daX7IkCEYMmSINB8UFAQAGD16NOzt7dGpUydUqlQJSqUSRYoUQVhYmMHrViqVOHLkiDTv7e2dYn0PDw906NABw4YNU/sb6OzsDFtbW3To0AH16tXD2bNn8cMPPwAA5syZg5EjRyIqKgq5cuXC6dOn0bhxY2zduhXffPMNihYtanD8GZLRLgbqIbP1gUqPp3f/+eefautMrl7S18aNG6c15sR1TNkHqnfv3iImJsY4O5qIiFK0atUqUbt2bZP/TTLFpFAoxKRJk6SnRHh4eIjQ0FCxZs0aUaZMGTFlypRUt9/Q4RxMkXcohEj/Uxjh4eGws7NDWFgYbG1t03v1elMoFCZfR9LDkNw6hRBqr40bNw7Lli3TqJe4Tp8+fbBt27Zk152W7ZPh7UNElO19/vwZnz59wvTp0/Hq1Su8evUK169flzssk6levToKFCiAf/75BwUKFEC/fv3QpEkTuLi4wM7ODubm5rhz5w4KFSoEe3t7AMDjx48RHR2N6tWrmyTv4CU8IiKiTMbKygpWVlZYvXq1WrkQAn/99ReOHz+OlStXyhSd8SXuF/b69WssWLAACxYs0FrXxsYGkZGRamX58uUzekxmRm+R9LZixYpkX0vuDaKiyxmg1M4w9ezZM9U2tHn8+LFByxERkWkoFAp4eHhgxYoVEF+GKoJSqYSXlxe8vb3h4+ODpk2bqi0zcuRImaI1jaTJEwCEhoYafT1MoJIRHx+fLusZN24cRo8enezrdnZ2BrV7/vx5nes2aNDAoHWULVvWoOWIiCj9KBQK9O/fH126dIGbmxtOnTolJVdCCKxatQrR0dE4d+4coqOjce3aNaxbtw4AsHLlSvTp00dqa9euXahfvz6AL3cVRkZGYvHixTh69Kgs2yYnXsLT4sGDB3BycsL333+Pn376yaTrUl2rTU6ZMmUMateYQyKcPXsWDx8+xHfffQcAyJUrV7b8sBARZVXW1tZo3LgxAMDFxQUuLi7o378/rK2tMXLkSHTq1Alv3rxB9+7d0b17d7VlJ0yYAAAIDAxE3rx5cfz4cZQtWxZFihSBlZUVrl69ijdv3qB58+bIkycP7t69i1evXqF48eIQQuDTp08ICQnBnTt34O/vj3379iE6OlpqP0+ePIiIiEi/naEjdiL/f+/fv0fjxo3RvXt33LlzB3/++ScAzU7bxqa63TOps2fP4ubNmxg9ejTMzP47UZg0njFjxmD58uVa21bVW7NmjdqYIEm9ePECJUuWTPb1uLg4WFhYQKlUQggBc3PzVLeLiIjIEAkJCbCw+HJ+59OnT8iZMyfevHkDIQSio6ORN29edO/eHb169UKlSpVgZ2cHJycnvHr1CvPnz8fatWsBANOmTUN8fDwWLlyI0qVLIyAgwKh5BxOo/zd16lStZ5sSEhJMmjDosvvLlSuHJ0+eSPV1TaAePXqEc+fOYcCAAaluQ+I2ra2tERkZiRo1aiBPnjy4ePFiutyJSEREBADv3r1DfHw8ihQpotdyQgjcvn0blSpVQo4cOQAAYWFh+Pz5MxwcHHgXnikkd6lu8uTJJlvnyZMndarXp08fzJkzR+trKSVg5cuXR/ny5fWOS3Xq9NatW1AoFEyeiIgoXaXWvSU5CoUC1atXVyuzs7NDeHi4McJSw07kqViyZIle9StVqqRz3WbNmulULz0fjZJ4pFgzMzMmT0RERFowgTKyc+fOGb3NJk2a4NSpU3jx4oXR21ZRnSZt0aKFydZBRESUVfASnhH5+vrC3t4eV69eRa1atbTWsbW1RXh4OEaNGqVX20nH7VAxVhe2f//9F3v37sW3335rlPaIiIiyMiZQMN6YT19//TUAoGbNmnj16hUcHR016oSGhqrdVZdWBQsWNEo7xYsXx9ixY43SFhERUVbHS3gAYmNjjdJO4v5CxYoV0+gP9euvvxoledq7dy8sLCzQunVrTJw4Mc3tERERkX54BgqQBg8ztgYNGuD+/fvSvGogyrTq1KkT4uLijNIWERER6Y9noACDnmBdunRpeHt7p1hn/vz5hoZEREREGRgTKANduHABXbp0wevXrzF58mQEBgZq1HFwcICLi4sM0REREZEp8RKegVS3/RctWhQLFy5Mtl7+/PnTKyQiIiJKJ9n+DJSDg4PW8jlz5qBAgQJQKpUar+nzUMNVq1ahePHiWL16tcExEhERUcaSbRKosLAwreXBwcEaZYsXL8aMGTPw/v17KBQKODk5Sa8dPXoUNjY2Oq+3YsWKePHiBUaMGKF/0ERERJQhZfkEKiEhAQqFAnnz5tV4LElyjykpWbKk2nzNmjWl/7dq1cr4QRIREVGmkuUSKCEEDh8+jN27dwMALCzUu3mpHlUSExOTbBtJl1myZAlGjhyJf/75x8jREhERUWakEMZ6FogewsPDYWdnh7CwMNja2hq1bV0efpuQkABzc3O1Mn9/f+kJzgcPHkS7du2MGhcRERHJwxR5R5Y5AzVnzhydkicAqFOnjtp8fHy82rJJz0ARERERJZbpM4U1a9Zg5MiRqdYTQkhJ0rVr16Ty3r17w9zcHIUKFZLKknsQMBERERGQBRKoJUuWpFpHdZXy48ePyJs3r1ReunRpbN++HQBgb2+P/fv3QwgBe3t7k8RKREREWUOmv4SX9DEs169fx4ULF9C1a1ds2bIFibt42dnZqZ1pmjFjhtqyHTp0QMeOHU0aLxEREWV+Wa4TeWrOnz8vPTz4/fv3KFCgQLqun4iIiNKXKfKOTH8JT1+NGjXC06dPYWNjw+SJiIiIDJLtEijgS98nIiIiIkNl+j5QREREROmNCRQRERGRnphAEREREemJCRQRERGRnphAEREREemJCRQRERGRnphAEREREelJlnGgVIOfh4eHy7F6IiIiykZU+YYxH74iSwL14cMHAEDx4sXlWD0RERFlQx8+fICdnZ1R2pIlgcqfPz8A4MWLF0bbkIwsPDwcxYsXx8uXL9P92X9y4PZmfdltm7m9WVt2214g+21zWFgYSpQoIeUfxiBLAmVm9qXrlZ2dXbY4cCq2trbc3iwsu20vkP22mdubtWW37QWy3zar8g+jtGW0loiIiIiyCSZQRERERHqSJYGysrKCp6cnrKys5Fh9uuP2Zm3ZbXuB7LfN3N6sLbttL5D9ttkU26sQxrynj4iIiCgb4CU8IiIiIj0xgSIiIiLSExMoIiIiIj0xgSIiIiLSk0kSqNWrV6NUqVKwtrZG3bp1cfny5WTrenl5QaFQqE3W1tamCMskzp07h3bt2qFo0aJQKBTYv39/qsv4+vrCxcUFVlZWKFeuHLy8vEwepzHpu82+vr4ax1ihUCAwMDB9Ak6DBQsWoHbt2siTJw8KFSqEjh074sGDB6kut2fPHlSqVAnW1taoXr06/vrrr3SI1jgM2ebM/Dleu3YtnJycpAEF69evj6NHj6a4TGY+vvpub2Y+ttr89NNPUCgUGDduXIr1MvMxTkqXbc7Mx3nWrFkasVeqVCnFZYxxfI2eQO3atQvjx4+Hp6cnrl+/jho1aqBly5YIDg5OdhlbW1u8fftWmp4/f27ssEwmKioKNWrUwOrVq3WqHxAQAA8PD7i7u+PGjRsYN24cBg8ejOPHj5s4UuPRd5tVHjx4oHacCxUqZKIIjefs2bMYOXIk/vnnH5w8eRJxcXFo0aIFoqKikl3m4sWL6NWrFwYNGgQ/Pz907NgRHTt2xO3bt9MxcsMZss1A5v0cOzo64qeffsK1a9dw9epVNGnSBB06dMCdO3e01s/sx1ff7QUy77FN6sqVK1i/fj2cnJxSrJfZj3Fium4zkLmPc9WqVdViP3/+fLJ1jXZ8hZHVqVNHjBw5UppPSEgQRYsWFQsWLNBaf9OmTcLOzs7YYcgCgNi3b1+KdSZPniyqVq2qVtajRw/RsmVLE0ZmOrpss4+PjwAgQkND0yUmUwoODhYAxNmzZ5Ot0717d+Hh4aFWVrduXTF06FBTh2cSumxzVvocCyFEvnz5xMaNG7W+ltWOrxApb29WObYRERGifPny4uTJk8LV1VWMHTs22bpZ5Rjrs82Z+Th7enqKGjVq6FzfWMfXqGegYmNjce3aNTRr1kwqMzMzQ7NmzXDp0qVkl4uMjETJkiVRvHjxVH8JZXaXLl1S2z8A0LJlyxT3T1bh7OyMIkWKoHnz5rhw4YLc4RgkLCwMAFJ8IGVWO8a6bDOQNT7HCQkJ2LlzJ6KiolC/fn2tdbLS8dVle4GscWxHjhwJDw8PjWOnTVY5xvpsM5C5j/OjR49QtGhRlClTBn369MGLFy+SrWus42vUBOr9+/dISEiAg4ODWrmDg0Oy/V0qVqyI33//HQcOHMC2bdugVCrRoEEDvHr1ypihZRiBgYFa9094eDiio6Nlisq0ihQpgnXr1uHPP//En3/+ieLFi8PNzQ3Xr1+XOzS9KJVKjBs3Dg0bNkS1atWSrZfcMc4Mfb6S0nWbM/vn2N/fHzY2NrCyssKwYcOwb98+VKlSRWvdrHB89dnezH5sAWDnzp24fv06FixYoFP9rHCM9d3mzHyc69atCy8vLxw7dgxr165FQEAAGjdujIiICK31jXV8LQyO2Ejq16+v9sunQYMGqFy5MtavX4+5c+fKGBkZS8WKFVGxYkVpvkGDBnjy5AmWLVuGrVu3yhiZfkaOHInbt2+neG09q9F1mzP757hixYq4ceMGwsLC4O3tjf79++Ps2bPJJhWZnT7bm9mP7cuXLzF27FicPHky03SKTitDtjkzH+fWrVtL/3dyckLdunVRsmRJ7N69G4MGDTLZeo2aQBUsWBDm5uYICgpSKw8KCkLhwoV1aiNHjhz46quv8PjxY2OGlmEULlxY6/6xtbVFzpw5ZYoq/dWpUydTJSKjRo3C4cOHce7cOTg6OqZYN7ljrOtnIKPQZ5uTymyfY0tLS5QrVw4AULNmTVy5cgXLly/H+vXrNepmheOrz/YmldmO7bVr1xAcHAwXFxepLCEhAefOncOqVavw+fNnmJubqy2T2Y+xIducVGY7zonlzZsXFSpUSDZ2Yx1fo17Cs7S0RM2aNXH69GmpTKlU4vTp0yleX08sISEB/v7+KFKkiDFDyzDq16+vtn8A4OTJkzrvn6zixo0bmeIYCyEwatQo7Nu3D2fOnEHp0qVTXSazH2NDtjmpzP45ViqV+Pz5s9bXMvvx1Sal7U0qsx3bpk2bwt/fHzdu3JCmWrVqoU+fPrhx44bWRCKzH2NDtjmpzHacE4uMjMSTJ0+Sjd1ox1evLuc62Llzp7CyshJeXl7i7t27YsiQISJv3rwiMDBQCCFE3759xZQpU6T6s2fPFsePHxdPnjwR165dEz179hTW1tbizp07xg7NJCIiIoSfn5/w8/MTAMTSpUuFn5+feP78uRBCiClTpoi+fftK9Z8+fSpy5colJk2aJO7duydWr14tzM3NxbFjx+TaBL3pu83Lli0T+/fvF48ePRL+/v5i7NixwszMTJw6dUquTdDZ8OHDhZ2dnfD19RVv376Vpk+fPkl1kr6nL1y4ICwsLMTixYvFvXv3hKenp8iRI4fw9/eXYxP0Zsg2Z+bP8ZQpU8TZs2dFQECAuHXrlpgyZYpQKBTixIkTQoisd3z13d7MfGyTk/SOtKx2jLVJbZsz83GeMGGC8PX1FQEBAeLChQuiWbNmomDBgiI4OFgIYbrja/QESgghVq5cKUqUKCEsLS1FnTp1xD///CO95urqKvr37y/Njxs3Tqrr4OAg2rRpI65fv26KsExCdYt+0km1jf379xeurq4ayzg7OwtLS0tRpkwZsWnTpnSPOy303eaFCxeKsmXLCmtra5E/f37h5uYmzpw5I0/wetK2nQDUjlnS97QQQuzevVtUqFBBWFpaiqpVq4ojR46kb+BpYMg2Z+bP8bfffitKliwpLC0thb29vWjatKmUTAiR9Y6vvtubmY9tcpImE1ntGGuT2jZn5uPco0cPUaRIEWFpaSmKFSsmevToIR4/fiy9bqrjqxBCCP3OWRERERFlb3wWHhEREZGemEARERER6YkJFBEREZGemEARERER6YkJFBEREZGemEARERER6YkJFBEREZGemEARERER6YkJFBGZzKxZs+Ds7Cx3GBKFQoH9+/frvdyDBw9QuHBhREREGD+oRN6/f49ChQrh1atXJl0PEaUdEyiiTG7dunXIkycP4uPjpbLIyEjkyJEDbm5uanV9fX2hUCjw5MmTdI4yfRk7cZs6dSpGjx6NPHnyGK1NbQoWLIh+/frB09PTpOshorRjAkWUybm7uyMyMhJXr16Vyv7++28ULlwY//77L2JiYqRyHx8flChRAmXLlpUj1EzpxYsXOHz4MAYMGJAu6xs4cCC2b9+OkJCQdFkfERmGCRRRJlexYkUUKVIEvr6+Upmvry86dOiA0qVL459//lErd3d3BwBs3boVtWrVQp48eVC4cGH07t0bwcHBAAClUglHR0esXbtWbV1+fn4wMzPD8+fPAQAfP37E4MGDYW9vD1tbWzRp0gQ3b95MMd6NGzeicuXKsLa2RqVKlbBmzRrptWfPnkGhUGDv3r1wd3dHrly5UKNGDVy6dEmtjQ0bNqB48eLIlSsXOnXqhKVLlyJv3rwAAC8vL8yePRs3b96EQqGAQqGAl5eXtOz79+/RqVMn5MqVC+XLl8fBgwdTjHf37t2oUaMGihUrJpVpO8P1yy+/oFSpUtL8gAED0LFjR8yfPx8ODg7Imzcv5syZg/j4eEyaNAn58+eHo6MjNm3apNZO1apVUbRoUezbty/FuIhIXkygiLIAd3d3+Pj4SPM+Pj5wc3ODq6urVB4dHY1///1XSqDi4uIwd+5c3Lx5E/v378ezZ8+ksyxmZmbo1asXduzYobae7du3o2HDhihZsiQAoFu3bggODsbRo0dx7do1uLi4oGnTpsmePdm+fTtmzpyJH3/8Effu3cP8+fMxY8YMbN68Wa3etGnTMHHiRNy4cQMVKlRAr169pEuUFy5cwLBhwzB27FjcuHEDzZs3x48//igt26NHD0yYMAFVq1bF27dv8fbtW/To0UN6ffbs2ejevTtu3bqFNm3aoE+fPime7fn7779Rq1atFPd/cs6cOYM3b97g3LlzWLp0KTw9PdG2bVvky5cP//77L4YNG4ahQ4dq9HmqU6cO/v77b4PWSUTpRBBRprdhwwaRO3duERcXJ8LDw4WFhYUIDg4WO3bsEF9//bUQQojTp08LAOL58+da27hy5YoAICIiIoQQQvj5+QmFQiHVT0hIEMWKFRNr164VQgjx999/C1tbWxETE6PWTtmyZcX69euFEEJ4enqKGjVqqL22Y8cOtfpz584V9evXF0IIERAQIACIjRs3Sq/fuXNHABD37t0TQgjRo0cP4eHhodZGnz59hJ2dnTSfdL0qAMT06dOl+cjISAFAHD16VOs+EUKIGjVqiDlz5qiVaWt/2bJlomTJktJ8//79RcmSJUVCQoJUVrFiRdG4cWNpPj4+XuTOnVv88ccfam19//33ws3NLdmYiEh+PANFlAW4ubkhKioKV65cwd9//40KFSrA3t4erq6uUj8oX19flClTBiVKlAAAXLt2De3atUOJEiWQJ08euLq6AvjS5wcAnJ2dUblyZeks1NmzZxEcHIxu3boBAG7evInIyEgUKFAANjY20hQQEKC1k3pUVBSePHmCQYMGqdWfN2+eRn0nJyfp/0WKFAEA6fLigwcPUKdOHbX6SedTkrjt3Llzw9bWVmpbm+joaFhbW+vcfmJVq1aFmdl/X7MODg6oXr26NG9ubo4CBQporD9nzpz49OmTQeskovRhIXcARJR25cqVg6OjI3x8fBAaGiolQ0WLFkXx4sVx8eJF+Pj4oEmTJgC+JDMtW7ZEy5YtsX37dtjb2+PFixdo2bIlYmNjpXb79OmDHTt2YMqUKdixYwdatWqFAgUKAPhyp1/Svlcqqv5IiUVGRgL40n+pbt26aq+Zm5urzefIkUP6v0KhAPClX5YxJG5b1X5KbRcsWBChoaGptpuQkKDTunRZf0hICOzt7VNdJxHJhwkUURbh7u4OX19fhIaGYtKkSVL5119/jaNHj+Ly5csYPnw4AOD+/fv48OEDfvrpJxQvXhwA1O7iU+nduzemT5+Oa9euwdvbG+vWrZNec3FxQWBgICwsLNQ6TyfHwcEBRYsWxdOnT9GnTx+Dt7NixYq4cuWKWlnSeUtLS60JjSG++uor3L17V6M8KChIbf7p06dGWR8A3L59W2MICiLKWHgJjyiLcHd3x/nz53Hjxg3pDBQAuLq6Yv369YiNjZU6kJcoUQKWlpZYuXIlnj59ioMHD2Lu3LkabZYqVQoNGjTAoEGDkJCQgPbt20uvNWvWDPXr10fHjh1x4sQJPHv2DBcvXsS0adO0JmPAlw7cCxYswIoVK/Dw4UP4+/tj06ZNWLp0qc7bOXr0aPz1119YunQpHj16hPXr1+Po0aPSmSpV3AEBAbhx4wbev3+Pz58/69x+Ui1btsSlS5c0ErLAwEDMmTMHT58+xZ9//omtW7ciNDQU9+/fN3hdAPDp0ydcu3YNLVq0SFM7RGRaTKCIsgh3d3dER0ejXLlycHBwkMpdXV0REREhDXcAAPb29vDy8sKePXtQpUoV/PTTT1i8eLHWdvv06YObN2+iU6dOyJkzp1SuUCjw119/4euvv8bAgQNRoUIF9OzZE8+fP1dbf2KDBw/Gxo0bsWnTJlSvXh2urq7w8vJC6dKldd7Ohg0bYt26dVi6dClq1KiBY8eO4fvvv1frp9SlSxe0atUK7u7usLe3xx9//KFz+0m1bt0aFhYWOHXqlFp5tWrV8PDhQ1StWhUzZszAxo0bYWlpiYkTJxq8LgA4cOAASpQogcaNG6epHSIyLYUQQsgdBBFRWnz33Xe4f/++yW79X716NQ4ePIjjx48D+DIO1P79+3Hjxg2jr6tevXoYM2YMevfubfS2ich42AeKiDKdxYsXo3nz5sidOzeOHj2KzZs3qw3IaWxDhw7Fx48fERERYdLHubx//x6dO3dGr169TLYOIjIOnoEiokyne/fu8PX1RUREBMqUKYPRo0dj2LBh6bZ+U56BIqLMgQkUERERkZ7YiZyIiIhIT0ygiIiIiPTEBIqIiIhIT0ygiIiIiPTEBIqIiIhIT0ygiIiIiPTEBIqIiIhIT0ygiIiIiPT0f4COOS3d1MJhAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "# Open the grid \n", "grid = ds['grid']\n", "\n", "# get wavelength \n", "wave = grid['wavelength'].values\n", "\n", "# for Exo-REM par1=Teff, par2=logg, par3=[M/H], par4=C/O\n", "# get flux by selecting a value for each parameter\n", "flux = grid.sel(par1=1700,par2=4.0, par3=0., par4=0.55, method=\"nearest\")\n", "\n", "# Plot the spectrum\n", "plt.figure(figsize=(7,3))\n", "plt.plot(wave, flux*wave, 'k')\n", "plt.xlabel('Wavelength (µm)')\n", "plt.ylabel('Flux')\n", "plt.xlim(0.5,5) \n", "plt.yticks([])\n", "plt.title('Exo-REM at Teff=1700K, logg=4, and solar [M/H] and C/O')\n", "plt.show()" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.7.13 ('exo_formosa')", "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.11.9" }, "orig_nbformat": 4, "vscode": { "interpreter": { "hash": "9aec1e17db1a1c759b18bf85d14c335cdcfb4080f740d60b236f2747bc3cb68f" } } }, "nbformat": 4, "nbformat_minor": 2 }