{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Example - 02 - Atmosphere" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "In this example, you will learn to use add an atmosphere model to planets." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Let us re-use the code from example-01 to create a planet object for Venus." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from AMAT.planet import Planet\n", "planet = Planet(\"VENUS\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "We are now ready to add an atmosphere model to the planet object." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "AMAT stores atmospheric data in the form of look up tables located in ~root/atmdata. Typically, data is stored in the following format." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "|#Z(m)\t|Temp(K)\t|Pres (Nm2)\t|rho(kgm3)\t|a (m/s)|\n", "|-------|-----------|-----------|-----------|-------|\n", "|0\t\t|735.30\t |9.209E+06\t| 6.479E+01\t| 428.03|\n", "|1000\t|727.70\t |8.645E+06\t| 6.156E+01\t| 425.46|\n", "|2000\t|720.20\t |8.109E+06\t| 5.845E+01\t| 422.88|\n", "|3000\t|712.40\t |7.601E+06\t| 5.547E+01\t| 420.27|\n", "|4000\t|704.60\t |7.120E+06\t| 5.262E+01\t| 417.63|\n", "|5000\t|696.80\t |6.666E+06\t| 4.987E+01\t| 415.09|" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "# the atmosphere model provided in atmdata/Venus/venus-gram-avg.dat.\n", "# The columns for height, Temp, pressure, density are 0, 1, 2, 3\n", "planet.loadAtmosphereModel('../atmdata/Venus/venus-gram-avg.dat', 0 , 1 ,2, 3)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcwAAADpCAYAAACpzQe3AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAAgAElEQVR4nOydd5xU1fXAv2fpvUhHikgRsCJFQFEsqNglJsaKjRhjikYTTTQqdlN+GhOjJhKxd+wGbKig0gQLSFMBQUBBOkvd8/vjvHFnhy0zu1Pemz3fz+d+dt67r5y7886cd+899xxRVRzHcRzHKZ+CXAvgOI7jOFHADabjOI7jJIEbTMdxHMdJAjeYjuM4jpMEbjAdx3EcJwncYDqO4zhOErjBdBzHcZwkcIPpOI7jOElQM9kDReQQ4ACgYfx+Vb0l3UI5jrMrroOOk1uSMpgicjfwY+A9oDCuysMEOU4WcB10nNwjyYTGE5Hvgb1V9ZvMi+Q4TiKug46Te5Kdw/wa2JpJQRzHKRfXQcfJMcn2MPsCfwAeB1bG16nqu5kRrWq0aNFCO3funPbrFhUVAVBQEH1/KW9L7pgxY8YqVW2Z7PGug1Unas8IuMyZJFUdhOSdfg4EjgWGsOv8ScdUbpgtOnfuzPTp09N+3cJCa369evXSfu1s423JHSKyOMVT8ksHV62C3XYDkazJE7VnBFzmTFIJHUx6SPYW4ARVbaGqHeJKKBW1XD77DDp3hpNOyrUkjpMK+aODb70FPXvCv/6Va0kcJyWSNZibgFAO+6TMpk2weDEsX55rSRwnFfJHB1essB7mb34D77yTa2kcJ2mSNZh/Au4UkTYiUhBfMilcRti2zf7WqZNbORwnNfJHB884w4zl9u1w4okwc2auJXKcpEhW2cYAFwPLgO1B2RH8jRYbN9rfhg3LP85xwkX+6CDAX/4CP/oRrF8Phx8OkyblWiLHqZBknX72yKgU2WTDBvvrBtOJFvmjgwA1asAjj1gv84UX4Mgj4a67YNSorDoCOU4qJNvDbKCqixMLsFcmhcsIsR5mo0a5lcNxUiN/dDBGnTrwzDNw6aWwdStcfDGceSasXp1ryRynVJI1mC+LSIk3XBE5AXgw7RJlGh+SdaJJ/uhgPDVrwt13w6OPQoMG8PjjsNde8NBDkMQaccfJJskazCuB8SLSFkBETgXuA47PlGAZw4dknWiSPzpYGmecAR99BEOHmgftuefCEUfAhx+64XRCQ1IGU1WfBW4FXheRnwP/AI5R1RmZFC4jbNlifx94wJaXOE4EyCsdLIvu3eHNN2HsWAtq8PbbMHAg9O9v+2KjQ46TI8o0mKW4ro8F/ou5tx8NfBZJl/aTT4YWLeDbb2GffczhwHFCSN7qYHmIwDnnwLx58LvfQfPmMH06jBwJLVvCqafCY48VjxQ5ThYpT9liLuvx5c9Aa2AWUXVpP/BAmD3bFG/DBjjlFPjrX3MtleOURn7qYDLsthvcfjssXQpjxsDgwTY6NG6cOQa1bGnRuh55BNaty7W0TjWhvGUl+eXGHk+rVuadd9tt8Ic/wBVXQO/ecMwxuZbMceLJXx1Mlnr14LzzrCxbBs89B08/bes2X3zRSu3aMGyYres8/ngzto6TAco0mIHLev4iApddZgYToGnT3MrjOAnkvQ6mSvv28MtfWlm+3IznM8/Au+/Cyy9bEYG+feHoo82IHnQQ1KqVa8mdPCG/5j9S5Ykn7G+fPjBgQG5lcRwnedq2hV/8whyDvvnGArkfeaQZx2nT4KabYMgQ81c45RS45x5k5kwLlOA4lSTZSD/5yT//aX9/9SuPLuI4UaV1awt6cPHFllzhnXdg/Hgr8+bB88/D889TF9C6dWH//aFfPyv9+0O3bhDy3I1OOKi+BvPLL837rlEj+MlPci2N4zjpoEEDGD7cCtjSsQkTYOJEiqZOpWDhQlvb+eGHxec0bmzDuDED2q8f7L67v0Q7u1B9DWZsOclxx0HdurmVxXGczNCpE1x0EVx0EVsLC2HNGurNnm3DtrGybJnl6HzrreLzWrcu2Qvt18+diZzKG0wReQuYDPxNVdekT6Qs8fbb9tc9Y52IEnkdzAXNmsFRR1mJ8c03JQ3otGmwcmWxI1GMPfYoaUD79PGIYdWMqvQw3wE6AG8CfdIjTpbYudM86wAOOyynojhOFYiuDoaJdu1sTedJJ9m2KnzxRbHxnDrVwvZ99ZWVp56y4woKoGfPkkZ0331tmYuTl1TaYKrqDekUJKusWGGLnVu1siEbx4kgkdbBMCMCXbta+elPbd+OHTBnTrEBnTYNPv3UgqDMng0PPmjH1a4N++1XbED79YMePSydmRN5kjaYIrIbMBxoq6p3iEg7oEBVl2ZMukwRiye7Y4eVmtV3KteJDnmlg1GjZk3rPe67L1xwge0rLISPPy5pROfNK+6ZxmjYsNipKFY6dXKnogiSlKUQkUOBZ4HpwGDgDqAbcAVwQsakyxQdOtgi6GXLLMrPLbdA/fq5lspxyiTvdDAfqFfPAiMcdFDxvnXrYMaMkkb0669h4kQrMVq2LDaeBx1k68CbNct2C5wUSXbx0Z3AT1T1GCx+JcAUoH+yNxKRS0VkuohsFZEHE+qOEJG5IrJZRN4WkU5xdXVEZIyIrBeRFSJyebL3LJPatS3fnohlee/YEa67zgKyO044yS8dzFeaNIHDD4ff/x6efRaWLLEpoJdegmuvNSfD3XaD776DV1+FG26AY4+1IPM9e1oIwPvug08+MV8LJ1QkazA7q+qbwedYcrptpDYH+g1wEzAmfqeItACeA64FmmNv0E/GHXI99ibdCRgK/E5Equ7aevjh8L//2Rve6tUwerQZzp/9DObPr/LlHSfN5J8OVhdat7YYt6NHw2uvmbH88kuLNHb55TBoENSpA3Pn2lzoxRfbPGjTptQePpya111n3rrffZfrllR7RJNIzioik4HRqjpeRL5X1eYiMgz4g6oeltINRW4CdlfVkcH2KGCkqg4KthsAq4ADVHWuiCwDzlPVCUH9jUA3VT29vPv06dNHJ0+eXLFAqhRMnkzNu+6ixiuv/LB7x6mnsv2ee2xRcxxbgvnPunmwdtPbkjvq168/Q1X7Jnt8Xutglgj1M7JtG/LxxxRMm0bBlCkUTJ1KQSn5eou6daNo4ECKBg+maNAgtEuX0M2Fhvr/HEeqOgjJv53+FnhZRF4B6onIfdi8yUkpylgavYGPYxuquklEvgB6i8hKoF18ffD55NIuFCj+KIAOHTokd3cRig4+mG0HH4zMm0fNO+6g5uOPU/O55ygaMoSdo0ZVqlGOk2byVwcdqF0b7dePnf36sfOSS2zfihXsnDyZmtOmUeujjyiYMYOCBQsoWLDAppQAbd2aokGD2DloEEUDB6L77utOjBkkqf+sqn4oIvsBZ2LDOV8D/dPkndcQSBxrWAc0Cupi24l1pcl5P3A/QN++fbVevXqpSbL//rYGC6BNG2qfeqpN7JdCytcOMd6W8FNtdDALhFGmUtljDwrbtEFHjKBGvXoWOH7mTEtt9t57MGkSsnIlNcaNo8a4cXZOw4YwcCAcfDAccogtb2nQICfiR+b/nAJJv4qo6jLMMy/dbAQaJ+xrDGwI6mLbWxLq0s/WrcUZTO66y7xpHSckVAsddMqmVi0zgP3729ynqvlbxBlQvvgCXn/dClhvs08fM54HH2yJuFu2zG07IkyZBlNEHqbYuaBMVPWcKsowGzg37r4NgD2B2aq6RkSWA/sBwRPAfsE56adOHcujN348nH22PYRXXWVLUBwny1RLHXSSR8SCIvToUbw2dPlymDy52IDOmmXLW6ZOhb/+1Y7Zay9LfTZ0qEU6a9MmZ02IGuV5yS4EvgjKOmzOogawNDjvJGBtsjcSkZoiUje4Rg0RqSsiNYFxwN4iMiKo/xPwiarODU59CLhGRJqJyF7ARcCDKbQxNZ580h6+7dvhH/+Azp0tk/v48VBUlLHbOk4pVE8ddCpP27b2e3XXXbYedO1a++269lozkPXqmTfu/fdbFKO2baFXL8st+uyzsGpVrlsQblS1wgKMBw5J2HcwMD6Z84Pjr8feluPL9UHdkcBcoBCYiLnQx86rg83ZrAdWApcnc78DDzxQq8Snn6qedppqQYGqDX6oduig23/+c93yyiuq27ZV7fohYPPmzbp58+Zci5EWotYWYLomqTtaXXUwzUTtGVHNgMzbtql+8IHqrbeqDhumWr9+8e9brOy7r+pvfqM6frxqYWHuZc4Qqeqgqia9rGQd0EJVt8ftqwWsVtXEuY9Q0LdvX50+fXrVL7Rsma2N+s9/YNGi4v1Nm8Kpp8Lpp9ubWwQ90woLC4H8mJyPWltEJNVlJdVXB9NE1J4RyILM27bZcO3bb1t5/33z5YhRv76tWR8+3AIsdO6ce5nTRKo6CMkHLpgJ3CIi9YIb1QNuBmalJmIEad8e/vhHm0z/4AO2X3EFRT172lDHmDEwbJg9RKNHW5ogx8kM1VcHncxRu7Y5A117reUDXbvW/l59tQVP2LzZgiZccomlN+vVC377W3jnnWoZiSjZHmZn4DGgL7AGaIZFAzlDVRdlTrzKk6m32x/enr76yuY7H3sMFi60ypo14eST7eE67LDQLShOJCpvgskQtbZUoofZGdfBKhG1ZwRCIPOyZRYR7bXXYMIE2BDnHN2ypaVEO+UUOOIIc5okBDInSWV6mEkZzLgbdATaAstVdUmK8mWVjBvM2MOgakMZ99wDzz9f/NbVoweMGgXnnAMtWqRdjnQQlQc7GaLWlsooa3BetdfByhK1ZwRCJvP27TZk+/LLMG6cjbrFaNTIwv+dcgqFhx0GDRuGQ+ZyyJjBFJEyh25VNZSuo1kzmPEsW2ZznfffXzw8W7u29TpPO80CL4coQ3uolLGKRK0tlehhug5Wkag9IxBimVXhs8/McD73nKU5i1XVrUvRscdS4+yzbd4zpCHyMmkwiyhjPZiqhjIzak4MZowdO+CVV8xwvvZacfSgunVtnecxx9hi4n32KTOSUDYIrTJWgqi1pRIG03WwikTtGYEIyfzll2Y8x42zdaAxmjSBESPgjDNsmipEibQzaTA7JexqC1wFvKSqD6Ryw2yRU4MZz5Il8NRT9hb2wQcl62rUsJQ+BxwA3btbUtlOnSxrSvv2Ftkjg0RGGZMgam2phMF0HawiUXtGIJoyb1mwgBrPPEOtp5+2UH4x2raFkSPhoovMgSjHZHwOM+FmTYBpqtq9UhfIMKExmPEsWwYvvGCGc+ZM+PzzsoMhFBSY0ezYsaQh7dzZHrZOnao81BFFZSyLqLWlsnOYCdeoljpYWaL2jEAeyDx3Ljz+ODz6aPGcp4itLrj4YjjxRPutywHZNpgdsGggoUwTHkqDmcjmzfDppxa+6quvYPHi4rJ8efFQblm0bVtsQDt3hq5dbZi3d++khnqjqIxlEbW2pMlgVksdrCxRe0Ygj2RWNYeh++6zEbfYWs+99rJk22eemfERtUQqo4NJrbYvJaZlfWAI8EgqN3MSqF8fBgywksi2bbB0qRnPJUuKDemiRVaWLDGjunz5rkO9BQXQrRvsu68Z0H33tdKpU87e5pyq4TroRBoRC/w+eDDceSeMHWvh++bOhfPOg+uug9/9Di688IflKWEk2TnM6xJ2bQJmqeobGZEqDUSih1kVduywId5Fi6x3+tVXMG+e9VjnzSt9UXGjRmZAAyO6tUcPinr3pl4eBF8OzfeSJJWYw3QdrCJRe0Ygz2Xevt2yQ916q01PgY2W3XijxbnN8Mt9Jp1+BqjqlFL291fVqancMFvkvcEsjy1b7AH85BMrn35qf1euLP34Tp3MiPboYc5H3btbD7Vdu9AHX4gRie8ljkoYTNfBKhK1ZwSqicxFRbaG/dprYc4c27fffnD77baqIENk0mCuLy1epYh8r6rNU7lhtqjWBrMsVq404xkY0KJZs5DPP0fiY0fG06CBGc5u3YqNaMyg7rZbdmWvgKh9L5UwmK6DVSRqzwhUM5l37oSHHoI//cmmo8DWsP/f/yUVwzZV0j6HGSyWFvsoEnyOsSewI2UpndzRurWVI48EYGthIezYQb2lS20R8vz5VhYssL+rV5tD0qxSwpU2a7arEY0Z10aNstyw/MV10Kk21Khh85mnn26pFUePtp7n//5n8byvuCLnQRAqcvrZQbGjQaJiFmHBn50oU7OmrQXt2XPXuu+/LzaeiX/XrIEpU6wk0qZNSSMa+9ylS84f+AjiOuhUL+rVgyuvNM/ZK66wZSnXXmt/H3oIDjwwZ6JVZDD3wN5o38E88mIo8J2qFmZKMCcENG9euhevKqxYUboxXbjQ6lasgHffLXmeiM2XJvZKY0EbIpgiLQu4DjrVk3btLLnFRRfZms05c+y36JprrMeZ5WUoUIHBVNXFwcfEKCNOdUbE1oC2bQtDhpSs27kTvv56V2M6f37xkphFi+D110ueV6uW9UBLG+Zt167aLodxHXSqPUOHWqCXP/zBlqLccAO8+qqt58zA3GZ5lGkwReR+VR0VfH6orONU9ZxMCOZElBo17CHu3BmOOqpk3bZttvylNGO6dKkth5k3b9dr1q9vQRkSjWmPHqFzPkonroOOE1C/vq3fPPlkOPdcmDbN4nE//DAcd1zWxCivh/lV3OcvyjzKcZKldm0zcj167Fq3ebMN58Yb0djn774rXiKTSIsW0KMHtfbcE+3RA/be26KH7LFHToZs0ozroOPEc9hh1ts85xxLcHH88Zbs+qabsjIKVabBVNVb4zbvU9UViceISPRXvDvhoH794ohEiaxdW3qvdN48WLUKVq2iZnyGBLD50K5diw30XnsV/20eylUYu+A66Dil0Lw5vPgi3HGHzWXeeqtFDHrkEfsdySDJelnMB3ZZAwbMAaLx6+NEl6ZNoV8/K/GoWt7RefPY9sknFCxYQM2FC82QLl5sSjR37q7XC3qlJYxojx42hxpexyPXQceJUVAAV11lvwkjRlhasUMPNUPatm3Gbpvsr8Mu4V5EpDHm1p5xRKQ58AAwDFgFXK2qj2Xj3k6IEbGMLu3bs3PgQHYCNWOLpTdvtt7ovHlmNOP/Br1SEnultWrBnnuWbkxz3yt1HXScRI44wmJpH3ccTJ9uXrTjx5e+TC4NVBS44GvMfb2eiCxJqN4NeDwjUu3KP4FtQGtgf+AVEflYVWdn6f5O1Khf38Jr7bdfyf2qFoM35mAUb0yXLCnulb7wQsnzWrTY1YjG5koz2Ct1HXScCujZ09aDn3yyZUQZMgQmTLA8w2mmIk0/C3uzfRU4O26/AitVtRSXxvQiIg2AEcDeqroRmCQiLwbyXJXp+zt5hgjsvruVI44oWRfrlSb2SGO90kmTrMQT65UmGtP09UpdBx2nIlq2tKVqI0ZYZKChQ80paPDgtN6monWY7wCISAtV3ZzWOydPd2Cnqs6P2/cxcGh5JxUVFf0Q0zCdbNmyJe3XzBXelgREipesxBP0SgsWLEDmzbO/8+cj8+dT8PXXZc6VFq5aVWUnBNfB9BHF591lTgEReOIJap93HjXGjUOHDWPbCy9QlEajWd46zNEJ26Uep6p/Sps0pdMQWJewbx2wS8BSERkFjALo0KFDhsVyqg1Br7Ro991h6FBKJE7btAlZuPAHYyoLFlAwb56tOa2isXQddJwUqV2bbWPHUqtBAwqmTqWoW7e0Xr68HmZYnvaN7Ood2BjYkHigqt4P3A+WKSGTEf6jlD2gIrwtVbqhzW8edFDJ/arUq3pqNNfBDBBGmSrCZU6RsWNhzRrqpTmwSXnrMM+r6OQgk0KmmQ/UFJFuqrog2Lcf4M4GTnhJQx5R10HHqSQFBRmJApZUPsxdThLZBzgHOFNV26Vdql3v9wTm5HAh5qH3KjCoPA89EfkOWFxWfRVpgbnW5wPeltzQSVVbVvZk18FKE6VnJIbLnBlS1sGk/eFFpCVwBnAu9nY5Cfh1SuJVnkuAMcC3wGrg5xW5s1flx6giRGR6qolHw4q3JTq4DladKD4jLnN4qGgdZi3gRGAkcDSwEFv31Qk4TVW/zbSAAKr6PXByNu7lOGHCddBxwkNF8x8rgfuAecBBqtpLVW/EFjA7jpN5XAcdJyRUZDA/AZoCA4B+ItIs8yJFgvtzLUAa8baEG9fB9BLFZ8RlDgkVOv2ISCfMueAcoCMwAVuw3FNVl2VcQsep5rgOOk44qNAlXVUXq+qNqtoNOAJYjgV8/lhE7si0gI5T3XEddJxwUNllJXWBU4BzVPXYtEvlOE65uA5mDhE5GLgp2GwHvKKql+VQpHIRkdbAOGA7sBNbarQ8t1Ltiog0AV4HemHz8Z/lWKSUqdSiZ1XdoqqP55uiikgdEXlARBaLyAYRmSkix8bVHyEic0Vks4i8HQyVxZ87RkTWi8gKEbk8N63YFRHpJiJbROSRuH1nBO3cJCLPB+mbYnXNRWRcULdYRM7IjeQlEZHTReTzQK4vROSQYH8kv5eqkK86GAZUdZKqHqaqhwHvA8/nWKSKWAUcrKqHAg8BF+RYnrLYDBwHPJNrQSpLNqKERImawNfY/FAT4FrgKRHpLCItgOeCfc2B6cCTcedeD3TD3P2HAr8TkWOyJ3q5/BOYFtsQkd6Y5+XZWLqmzcA9CcfHUjmdCfwrOCdniMhRwO3AeVgM0yHAlxH/XpwQEyzp6Q+8l2tZykNVd6pqLC9qI0IagUlVt6vqd7mWo0qoqpdyCualOAILKP1+3P4GQCGwV7C9DBgWV38j8EQI5D8deAozHI8E+24BHos7Zk/MQDYK2rUN6B5X/zBwW47b8T5wQSn7I/m9eEnbc3Ep9pK0FXgwoa45NlS5CYs4dEaK1x4O3B0FebHoS1Ow5Uedwvw/Bx7EUsXl/PlJtWQu820eEMwNdMfe2H6OpTQCQFU3icgXQG8RWYnNdXwcd/rH5Hiht4g0BkZjjiLxwzS9MQMEgKp+ISLbsLYWUYlUTplERGoAfYEXRWQhUBcbJrsSa0ukvhcnrXyDzTceDSRG+y4z6bWItKH0ocEfqeqK4PNpwH+jIK+qzgIGiMiPgauBi9MsdzyVakMG5ckabjDLIBiOeRQYq6pzRaQhkDicEEtx1DBuO7Eul9wIPKCqXyekhiovXdPOcupyRWugFvAj4BDMueEF4BqsLVH7Xpw0oarPAYhIX2D32P6Kkl4HRvHgsq4b6H8/0jwfmAl5RaSOqm4NNtdhUywZo7JtyKRM2cINZikEGSAext6ULg12l5fiaGPc9paEupwgIvsDRwIHlFJdXluKyqnLFbEsxHdr4P0nIn/DDOa7ROh7cbJGpZJex3Ek8JYWzw1mmqrI20dEbsdedrcA52dAvmSosA0i8irW8+whIvep6oPZFbFquMFMQKwr9gDWqxmuqtuDqtlY0OvYcQ2wub/ZqrpGRJZjAbFfDw7Jdfqjw4DOwJKgd9kQqCEivYD/YfIBICJdgDpYGqciQpbKKfj/LsWyZSQSte/FyQ5JJ70uDVV9DXgt3UKVQ6XlVdUPMCe4XFNhG1R1eFYlSjPuJbsr/wJ6AieoamHc/nHA3iIyIlgD9yfgE1WdG9Q/BFwjIs1EZC/gImxyO1fcjxmO/YNyL/AKNu/wKHCCiBwSGJjRwHOqukFVN2Fep6NFpIGIDAZOwnrcueS/wC9FpFUQHu43wMtE73txskPSSa9DQtTkLY18aEO5uMGMI1i/9zPMwKwQkY1BOVPNHXoEcDOwBovteXrc6dcBX2CeYe8Af1bV/2W1AXGo6ubAGWBFMP+xEdiiqt8FE/AXY4bzW+wN8JK40y/BJvO/xTJjVJjKKQvciC2NmQ98DswEbo7a9+JkjR+SXsftC/PoQtTkLY18aEO5VCrSj+M4ThgQkZrY1NJ1mAPKRcAOVd0hlUh67fJWTD60obJ4D9NxnChzDeYUdhVwVvD5mqAujCMlUZO3NPKhDZXCe5iO4ziOkwTew3Qcx3GcJHCD6TiO4zhJ4AbTcRzHcZLADabjOI7jJIEbTMdxHMdJAjeYWSQuEMJGESkSkcL44Ai5lq8qBMmZywxmnYX73yYi/4nb7hQkmf5zrmRyHCe/8FiyWURVY9kzEJFFwIWq+kbuJEoOEampqjuicg8R2RN4C3hIVa9NxzUdx3G8hxkiRKSGiFwrIl+KyCoReVREmgZ1e4nIDhG5QESWichqETlfRAaKyGcisjbI4BG71sUi8paI3Cci60VkjogMiatvLiIPBT3Dr0XkuiBLS/y5/xSRNcBVwf0nisj3IvKdiIwVkUbB8U8DrYAJQW/5VyJyjFjuyvj2/dALDXqEj4nIkyKyATi9vPan8D/sgYXAu8+NpeM46cQNZri4EhiG5b3bHcv7+H9x9TWAfYEuwHnA3cAVWPqcfYHzRGRA3PFDsPQ6uwG3Ac+LJZUGiyO7LrhWfyyp8tkJ584CWgB/DfaNBtoA+wA9gD8CqOppWGSPYaraUFX/nmR7RwBjgSbAs0m0vyJixvJvqnpLCuc5TrVDRFREuuZajijhBjNc/AxLFvuNqm4BbgB+EqQcizFaVbeq6ovB9kOqulpVlwDvUzL/5deqeo+qblfVh4ClwNFBkPkhwOVBkPblwN8pGbT8S1X9t6ruVNVCVZ2rqm+p6rYgmPudJJ9bsCzeUdVXVbUoyAyTTPvLYz9smuHZKsrlVENEZFGcX8FKEfmvWOL4akdgTDcF/4tlIvI3EamRa7lyjc9hhoTAKHQAXhWR+HiFBVgPESw56+q4ukJgZcJ2vIIvTbjNYqAd0AmoC3wXZ4sKgPgh1K8T5GsH3AUMwrKbFADLk2lbOfxwjyTavyqJ6z0FbAXeEpEhqrqsivI51Y8TVPUNEWkPjMdipF4Vf0DwrEoWk0vvQjb8CoD9VHVhkBZvIpaN5N4M3zPUeA8zJKgF9V0GHK6qTeNKXVVNxliUxu4J2x2BbzBDtRFoFnefxqraJ16khHP/DGwC9lbVxlg2Ainn+E1A/diGiNQCmicc88M5aWz/pcDbwJsi0jqF8xznB4KXrdeAvQGC+fubRWQysBnoIiJNROQBEVke9MJuivXCRKSriLwjIuuC+fgng/0iInwwmjAAACAASURBVP8nIt8GdZ+ISPw9LozJICIjRWRS3LaKyC9EZAGwINi3l4i8HvgWzBORH1emvSJycODLMLSU/8Vc4L24/8VVYh7oGwLfiFMqc88o4gYzXNwL3CYiHQDEkiWfUIXrdQgceGqKyFmYwZygql8BHwJ3iEgjESkQkW5S/rKQRpiRXS8iHYHLE+pXYvOhMT4HmovIEYGxvIGKn7dy2x84DZ1e5tn8YHhHYbkz3xCRFhXc03F2IXgGh2N5V2OcjT1bjbDRmrHADqArNhUyDHuRBMvfOgFohr243h3sH4ZNh3QHmgI/AeJHjSriZCznay+x5O+vA49hTnc/Be4Rkd4pXA8RORrLLDJCVd8upb4XcAjF/4svgu0mmF4/IiJtU7lnVHGDGS7uAN7AhhQ3YHOSfco/pVzexRT5e8xB5xRVXRfU/RRT2LlB/ZNAeT2yP2HOOOuAcew6T3gzcLOYt+6lQa/w15hz0VJgBRUPq5bZfhGpj/1QTa3gGgRDZSOBOZjnbrOKznGcgOdFZC0wCXMgi3cee1BVZwdDoc2BY4HfqOomVf0Wc1CLvdBtx6Y+2qnqFlWdFLe/EbAXNqz7eeBDkCy3qur3wZz/8cAiVf2vqu5Q1Y8wvfxRCtc7DbgfGK6qibr1kZiX/EvAf4D/Aqjq04GfQZGqPon1dvuncM/ooqpe8rAAFwNv5FqONLbnSOC/uZbDS/4WYBFwZBl1E4GL4rb7A0XA2riyHpgd1LcB/o1NgcwGzo8791fADOA7zFg1jrvHhXHHjQQmxW0r0C1u+3fAtgQZNgL/SrK9io0M/bmMuq5lnHcO5kEfu+cO4IJcf3/ZKO7040QCtQAPoQ/y4OQ18fP0X2MOZi20FOcbNU/yi8DmB7HpgXdVdaHasqu/i0grzFHtSuBaEub9MaNbkQzvqOpRVWjTacADIrJMVe+s6ODAw/7fwBHAB6q6U0RmUdKfIW/xIVnHcZwUURtGnQD8VUQaB34Ae4rIoQAicpqIxJzu1mCGbqeI9BORAcG8/iZgC7AzOG4WcKqI1A/WR15QgRgvA91F5GwRqRWUfiLSM5BhpFhEsfL4BjN+vxKRS5JoeoOgLd8F9ziPwBmoOuAGM09R1XtV9chcy+E4ecw5QG1srnwN8AwQc37pB0wRkY3Ai8Cv1ZztGmM9tDWY49Bq4C/BOf+HDbGuxByKHi3v5qq6AXMiOh0zfCuA24E6wSEdgMkVNUJtDfcRwO/jvXTLOHYOFsjkg0DOfZK5R74gwZi04ziOk0eIyATMUH+ea1nyBTeYjuM4jpMEPiTrOI7jOEngBtNxHMdxksANpuM4juMkgRtMx3Ecx0kCN5iO4ziOkwRuMB3HcRwnCZIOjScih2CBvEskVFXPbO84OcF10nGyS1IGU0TuBn6M5UQrjKvyRZyOkwNcJx0n+yQVuEBEvscSB3+TeZEcx6kI10nHyT7JzmHGIvM7jhMOXCcdJ8sk28PsC/wBy8q9Mr5OVd/NjGiO45SF66TjZJ9knX4OxLKLD2HX+ZKO6RYqHbRo0UI7d+6cazFKpaioCICCgmg4KUdNXgi/zDNmzFilqi2rcInQ62SqOlhUVMSqVato1apV5oTKAt6OcFFWOyqjg8kazFuAE4IkvpGgc+fOTJ8+fdeKDRvg0Udh1CjI0Y9pYaH9vtWrVy8n90+VqMkL4ZdZRBZX8RKh18kydbAMCgsLGTx4cErnhBFvRw4oKrLf9nXrYP16+7tuHQwdSiGU2o7K6GCyBnMTkB/DPJdfDv/5DzzxhP3t2jXXEjlOZcgfnXSqNzt27Gro4ksy+zdsgNKmFxcsgPbt0yZqsgbzT8CdIjIa+Da+QlWL0iZNNjjxRHjxRXjnHdh7b7j2WrjySqhdO9eSOU4q5I9OOtGnsBDWroU1a6yk8nnDhvTI0LAhNGlSsqR5FDFZgzkm+PuzuH2CzZfUSKtEmeaEE2D2bLjiChg7Fq65xoZo//lPGDo019I5TrLkj07Gcf755+dahLQQ2XaomgFbtQpZupQb+ve338lVq3Ytq1cXG7+tVXDYFilp5Bo33tXwVbS/cWOoUcZjX1iYtu8jWS/ZTmXVqWpV52IyQt++fbXCsfe33oKLL7ZuO8Dxx8Ptt0OvXhmVLezza4lETV4Iv8wiMkNV+1bh/NDrZFI6GEfYv7NkCVU7CgtLGrnvvivd+MWX7dtTv0+tWtCsWXFp2rT0z6XVNW6cUX+Ssr6Pyuhgsj3MBqo6J3GniBwNhEI5K8Xhh8Mnn8Bf/mKG8uWX4dVX4cIL4brroF27XEvoOGWRnzrpJMfWrbB8OSxbBt98Y6W0zxs3pn7tBg2gZUuKmjdHd9uNGq1bQ4sWu5bddis2fPXqWU8xz0nWYL4sIkeo6lexHSJyAnA/0DYjkmWLunVtWPaii+D66+Hf/4b774cHH4QLLoDf/x46lfky7zi5In91srqzdSssWQKLFln5+uuSRnDZMhsOTYbataFly9INXllGMOiJbQ1TTzkkJGswrwTGi8ihqrpcRE4F/gEcnznRskzr1vCvf8Gvf22OQM8+a9v//jeccw5cfbV71DphIv91Ml9JNIiJ5Zskoh3WqAFt29ooWPv29jf+c+xvkybVoueXLZIymKr6rIg0Bl4XkX8C1wLHqOonGZUuF+y1Fzz9NMyZA7fcAo8/DmPGWI9z+HD46U/N07Zhwwov5TiZolrpZNTYuBFZutTmCxcvTt0g1qgBHTpA585WOnTY1RC2bFm2k4uTMco0mCKSOAs7FmiOubMPA2aLSEHeurD36gWPPGLDtLfdZp5iL79spX59cxA65hgYNiyt63wcpyyqvU6GgW3bYOlSGyZdssT+xpclS6i3dm351ygoKGkQO3eGPfYo/ty+PdRMOvOik0XK+1Z2sGuqoFjffhZ54MKeFF27WoCDm2+2nufjj8P778NTT1kBW885bBgcfTQccsgPcwCOk2ZcJzOJKnz/vRnCJUusdxj7HCsrVpS+QD7+MnXqoO3bU9CxI3TsWNIYxgxirVrZaJGTZsozmHtkTYoo0Lo1XHqplUWL4KWXYMIEePtt+OwzK3/7mzkRDRlixnPYMOjd2+cQnHThOlkVYr3D0gxhrGzeXP41CgrM4HXsaL3EWInb3tKwIYi4s0weUqbBDMtarlDSuTP88pdWtm61HueECTB+PMycaZ8nTLBj27eHo46CQw+13meXLjkV3YkurpPlEN87LMsgJtE7pHFj84qP9Q4TS7t2FQ+XFhaWX+9EFh8oryp16liEoKFD4dZbYeVKeOMNM54TJpgL+IMPWgFo145agwdTNHgwHHGEzZWGNKOG44SOoiIb4fn8811LRXOHNWoU9w7LKk2aZKUZTjRxg5luWreGM8+0omqBEd58E957z8o331Dz6adtPhTM2/bAA6FvX+jXz/526eLDuI5TVAQzZtgL6KefmlGcOxe2bCn9+PJ6h5062TIMd6ZxqoA/PZlEBPbbz8rll9sPwNy5bHvjDQomT6bmBx+YZ90771iJ0ayZGdGYAe3b1+ZH3Ig6+c7KlfZi+b//2QjNqlW7HtO2rY3M9OxZsrRu7TriZBQ3mNmkoAB69WLnHnuw86KLqFmvnv1ATJ9eXKZNKx7WfSMu1WGrVsXGM2ZI27TJXVscJ5189BG1//AHaowfX3J/p062fGvAADOKe+1l8UcdJwdU2mCKyFvAZOBvqromfSJVM1q3huOOswI2jLtsWbHxjBnSb7+1OLevvlp8bvv2JQ1o374W2sqplkRSJ1evhl/8Ap58khqA1q2LHHaYGcmjj4YePbzX6ISGqvQw3wE6AG8CfdIjjoMI7L67lZNPtn2q5ugQb0BnzDDDumwZvPBC8fl77FGyJ9qnjzsyVB+ipZMffgg//rFNS9Sty/af/Ywdl19OvY4dcy2Z45RKpQ2mqt6QTkGcchAxQ7jHHvYDAzYfumBByaHcmTPhq6+sxJyKALp3L9kTPeAAy0jg5BWR0snp0+Gww2xZVv/+8OST7GjdOtdSOU65JG0wRWQ3YDjQVlXvEJF2QIGqLs2YdE7ZFBTYcFWPHuaRC7Bjh3kRxg/nzpoF8+dbeeyx4nN79iw2oAMHwr77ugdhxIisTq5fD6edZsbynHMswUHt2r5+0Qk9Sf1CisihwLPAdGAwcAfQDbgCOCFj0jmpUbOmhenbe28YOdL2bdtmUYjiHYs+/RRmz7YSWx/aoIE5VgwaZGXgQHeuCDGR1snJk22KAeDee81YOk4ESLZLcSfwE1V9U0RizgRTgP6ZEctJG7Vr2zxmnz4wapTtKyy09aHTp8OUKfDBB7BwIbz1lpUYvXvDoEHU6NePooMOMkPsDhhhIbo6OXCgLQ1ZvtyGYy+7DE4/3Z8tJ/QkG2Kms6q+GXyOxZbaRmpDupeKyHQR2SoiDybUHSEic0Vks4i8LSKd4urqiMgYEVkvIitE5PJk7+mUQb161pv8xS/goYdsLnTlShg3Dq68EgYPNkM7ezb8+9/UHjWKuvvuax69J58Md9wBkyaVvYDcyQZV0smc6mPTppYyr3lzG/244AJo355av/89Mm9eypdznGyRrMGcIyJHJ+w7Evg0hXt9A9wEjInfKSItgOewfH7NsSGmJ+MOuR4bauoEDAV+JyLHpHBfJxlatSppDNevtxi5f/kLO088EW3VyvL7vfAC/P73Fhe3cWPrLfz2t/Dccxar08kWVdXJ3OrjMcdYIPQHH7SXt7VrqXn33dQ94AB7YfvPf+wZdJwQIVpRMGJARA4CXgZeAX4MPITNk5ykqtNSuqHITcDuqjoy2B4FjFTVQcF2A2AVcICqzhWRZcB5qjohqL8R6Kaqp5d3nz59+ujkyZNTES1rbAl6ZnXr1s2xJMmxZcsWUKXe8uUUfPjhD0Vmz0YSnp+izp0pGjiQogEDKDroILR375wkug37/7h+/fozVLVvZc9Pl05mUh9T0UH56CPkvvuoPW4csnEjANqyJVvfew+N2DKTsD97yZLv7aiMDiY1fKOqH4rIfsCZ2Bvp10D/NHnj9QY+jrvXJhH5AugtIiuBdvH1weeTS7tQoOyjADp06JAG0ZwfEEG7dGFnly7sPOMM27duHQXTplHwwQcUTJlCwdSpFCxaRMGiRZY3FNBGjcx4DhxI0eDBFPXtawm4nSqRQZ1Mmz6uWrWKwYMH/7B9/vnnc8EFF5R6U+3Th8K77mLzLbfQ5MQTKZg5E/nuO+TzzyNnMJ3w8sADDzBmzA+DKi1SPT/pOUhVXYZ54qWbhsB3CfvWAY2Cuth2Yt0uqOr9wP0Affv21bDnowu7fInsIm+9enDCCVYAdu40D9z337cyeTKyaBE13niDGrEwf7VqWZzcgw+2Yd3BgzManShq/+NUyJBOpk0fW7VqxfTp01O6ecHrr1Mwc6ZtXHYZdU46KbLOQPny7OVTOy699FIuvfRSAESklEDF5VOmwRSRh9k1u/suqOo5qd40gY1A44R9jYENQV1se0tCnRM2atSA/fe3cskltu+bb2xOdNIkC6r98ccW4eXDD+Evf7FjevWyRexHHGFp0po1y1kTwkyWdDKn+ljw5ZfFG/vvD9u3+7ITJzSU5/SzEPgiKOuwYZcawNLgvJOAChLQJcVsYL/YRjBnsicwO4iHuTy+Pvg8Ow33dbJBu3YWnejvf7dIRGvWwGuvwR//CEOGWD7ROXPgnntgxAjrbfbta45Fr78OmzfnugVhIhs6mVN93HHuueZIBnDuuZas/aabitdtOk4uUdUKCzAeOCRh38HA+GTOD46vCdQFbgUeDj7XBFpiyj8i2Hc78GHcebdhMTKbAXthCntMRfc78MADNaxs3rxZN2/enGsxkiaj8m7Zojppkuro0aqHHqpaq5aqRc+1Uru26lFHqd51l+rCheGQOQ0A0zVJ3SmtVFUns6GPqergD9/Ztm2qY8ao7r13yWdh4EDVv/9ddfHiSv3Ps0XYn71kyfd2VEYHk1XOdUCthH21gPVJ38jc0TWhXB/UHQnMBQqBidgas9h5dTCnhvXASuDyZO7nBjN9ZFXeTZtUx49XvfJK1T59VEVK/mj27Gl1Eyeqbt8eDpkrQRoMZpV0Mhv6WGmDGaOoyJ6F009XrV+/5HPQq5fqZZepvvFGuc9BLgj7s5cs+d6OyuhgsstKJgLTgD+paqGI1ANuAA5S1SEVXiAH9O3bV1N1OMgWhUHMzKhMpudU3lWrLJnwyy/b33Vx/ia77QYnnQQ/+pHNf8bNdYX9fywiVV1WMpGQ62SqOljud7ZpE7z4Ijz1lOWJ3bixuK5FCzjlFItPO3RozmMih/3ZS5Z8b0dldDDZwAUjsXiV6wLX8nXY8E9VHX4cp3xatICzzoInnrDACRMnwhVXWND51astYszw4RZ44dxz7Ue1ekQgGkl10skGDeCnP7VoVKtX23Nw1VWWiWfVKgvgPmyYZfS58UYLu+c4aSapHuYPB4t0BNoCy1V1ScakSgPew0wfoZV3zhx45hkrn8YFuGnYkB3Dh7Pzxz+mzokn2lKWkFHVHmbcdUKrk2ntYZaFqn33zzxja38XLrT9NWta5KqLL4bDD8/q0pTQ6kuK5Hs7KqODyQ7JltkTVdWiVG6YLdxgpo9IyDt/Pjz7rP1wfvRR8f4WLeAnP7Fe6oABoVnTl4Yh2dDrZFYMZjxFRZY84F//shCOO3fa/u7d4Wc/sxGIDK75jREJfUmCfG9HJodkdwDbyyiOk3u6d4err4YZM+CLL9h+ww0U9expw3X//KctVejWDa6/HpaEqiNWWVwnEykogCOPtBenxYvtu27f3l6mfvtby5By/PGWcGBtOlbEOdWNZA3mHkCXuDIYeIkgDJ3jhIouXdhx5ZVsnT7d1n7Gfiy/+AJuuMHmuU48EV591Xol0cR1sjzat4frrrP1m+PGwdFHW4/zlVesp9m6NRx6qK33fe45S7y+bVuupXZCTkpzmCVOFGkCTFPV7ukVKT34kGz6iJq8UIrMO3facN1//2vDttuDjtg++8A115iHZRaHa9M1h5lwzVDpZNaHZCti5Uozjk89Be++u+vLUkEBdOwIXbvaaETXrsWf99gDkgxCHkV9KY18b0fG5jBLPVGkA/CJqoYyjpkbzPQRNXmhApm//da8a//xD1i2zPYddBDcdZclNM4CGTKYodLJ0BnMeFatKg7ROH26DdsuXlz2iIMIdOhQ0pD27GlJ1jt2NGObi3ZkkHxvR2V0MNlks4kxLOsDQ4BHUrmZ44SCVq1sScJll1k+xuuvtx/OAQNsuO7WW20IN8S4TlaRFi1sPvP444v3bdtmQ7gLFpi37cKFxZ8XLbK57yVL4M03S16rQQOLh9y7N+yzDwV7703RfvtZcgInr0h2he/ChO1NwL2q+kaa5XGc7FGnjnlPnnEG3HIL/O1vMHasDdvddJMFkM/xIvhycJ1MN7Vrm/NY91JGtLdvN6MZM6Tz59uyptmzbah32jQrWCgkFbEe6GGH2bKWQw81I+1EmmSXlQxQ1Sml7O+vqlMzIlkV8SHZ9BE1eaGSMn/xhfU6X3rJtg84AO69NyPDtGlYVhJ6nQz1kGw6Wb262HjOmkXR1KnIZ58h2xMclvfd1wznoEFWQp7nM7LfRwK5WIe5XlUTU/4gIt+ravNUbpgt3GCmj6jJC1WU+cUX4Ze/tOE3EbjySosek8Y0U2kwmKHXyWpjMBMoLCyErVup99ln8PbbVt5/H7ZuLXng7rsXG89Bg2C//UKVyiyvvg+ysA5TRApEpIZ9FAm2Y6UbthbMcfKLE0+0HsPvfmcG84477Adt3rxcS+Y6GRXq1LEk6ddea97Za9fa3xtvhGOPhaZNYelS89j9zW9sFKNRI5tH/+Uvba3o3LlRXvaUl1Q0QbODYseCREUsAm5Ou0SOEwYaNIDbb7fg7meeaQER+vQx79qf/CSXkrlORpG6dS0w/NChtl1UZAZx8mTrfb7/vs2LTp1qJUbjxnDggWZQ+/Wzv7vvHpqIVdWNigzmHoBg+e/iMyAo8J2qFmZKMMcJBYMGwaxZ5gD02GNw+ukWu3T06BJLCbKI62Q+UFBgnrW9esFFF9m+tWttiUvMgWjqVFv2FBvWjdG6dUkD2r8/NAvFSqK8p1yDqaqLg4+dsiCL44STJk3gkUdsuOyyy+Dmm23I9rHHkl7Mni5cJ/OYpk0ttN+RRxbvW7682HjGDOnKleaYFnNOA/PsHTCguOy7b6jmQ/OFMg2miNyvqqOCzw+VdZyq5mc6IceJRwR+9StbKvDjH1u4tWOPtSDfjXfxvcmQCK6T1Y62bW1O/cQTbVvVvLljRnTKFEs2MH++lYcftuPq1LEphP79i43oHnv4UG4VKa+H+VXc5y8yLYjjRIKjjoL33rPcixMnWuLq117L1ho718nqjkhxpKGf/tT2bdsGn3xixnPKFDOk8+bBBx9YidGypRnQwYOt9OvnwRVSJNllJW1UdUWy+8OALytJH1GTF7Ig85dfmvH88kubh5o40X6QkiQNy0pCr5PVelkJOW7HmjXWC40Z0SlTLBxgPLVqQd++5s0bM6JxL36haEcayHpoPGA+UNq40xwgFGu+HCerdOkCkybZfNOcOdbjfOutbDpfuE46ZdOsmT2Tw4bZtip89ZWFgJw82Z7dTz8t7oX++c92XK9eMGQIHHKI9UB33z13bQghyRrMXQa+RaQx5sbuONWTtm0truiQIeZJe8wx8Prr2ZrTdJ10kkfEXvK6dLFQkGBeuR9+aMZz0iTrhc6ZY+Xee6kHFHXubNGJDjnEnvOuXav1PGi5BlNEvsbc1euJSGLW3d2AxzMlWIIczYEHgGHAKuBqVX0sG/d2nHJp06bYaE6dCqeeanOatWpl5Hauk07aaNrUXvKOOca2t2619cbvvgvvvYdOmkTBokUWQ3fsWDumTZti4zlkCOy9d66WV+WEinqYZ2Fvsq8CZ8ftV2ClqmYr9Mk/gW1Aa2B/4BUR+VhVZ2fp/o5TNh06mNEcOND+XnIJ3H9/pt7EXSedzFCnTnGYvquuYsvGjcinn1J36lRzdHv3XVixAp5+2gqY0T344GID2qdPxl4Ww0CyTj/1VXVzFuQp7d4NgDXA3qo6P9j3MLBMVa8q67w+ffro5MmTsyRlamzZsgWAullew1dZoiYv5EZmmTqVOsccg2zZwvabb2bHZZeVeWz9+vWr6vQTep1MVQej+JyVRt62QxWZP5+CyZOtvPceBUuXljhH69enaMAAigYPZudZZ6EhCDBf1vdRGR0s02CKyOhkLqCqf0rlhqkiIgcA76tqvbh9VwCHquoJCceOAkYBdOjQ4cB5IYj9WRpRU6ioyQu5k7ngueeoc9ZZaEEBW2fNQrt2LfW4SilrxHSyY8eO2iLO6/L888/nggsuKPO6W7Zs4fHHH+e8887LjOBZojq1Q5YsMeM5aZL9nT+/+PyJE9EsJWQvj/h2PPDAA4wZMwaAmTNnLlbVzqlcqzyD+d9kLqCqGX0qROQQ4GlVbRO37yLgTFU9rKzzfFlJ+oiavJBjmf/6V+jUCX70ozIPqVRqoYjpZGWWlQwePJiPPvooneJmnWrdjpUri52I7rgjFMOzZbUjrctKklE6EcnGbO9GdnWfbwxsyMK9HSd1fvvbjFzWddIJPa1bw4gRVvKQSqWTF5F9gHOAM4F2aZVoV+YDNUWkm6ouCPbtB5TrXDBjxoxVIrK4vGNyTAvMuzAqRE1eCLfMaY0FG0adrKQOthCRsH5nqeDtCBeltSNlHUzK6QdARFoCZwDnYsoxCfiHqj6d6k1TRUSewLwAL8Q88l4FBkXZI09EplfF6SPbRE1eiKbMqeA66TjZpaJ1mLWAE4GRwNHAQmydVyfgNFX9NtMCBlwCjAG+BVYDP3fFdKojrpOOkzsqGpJdiUUOeRC4TlU/AhCRSzIsVwlU9Xvg5Gze03FCiuuk4+SIihwEPgGaAgOAfiLiWUrTx/25FiBFoiYvRFPminCddJwcUa7BDFzE9wQmAFcAK0TkJaABkHt/4QijqpH6MY+avBBNmSsiijopInVE5AERWSwiG0RkpogcG1d/hIjMFZHNIvK2iHRKOHeMiKwXkRUicnluWlGMiDwiIssDmeaLyIVxdVFrSzcR2SIij8TtOyP4rjaJyPNBGMRYXXMRGRfULRaRM3IjeTEiMjFow8agzIurS29bVDXpAhyMvbWvxbwP70jlfC9evKS3REEnMWN+PdAZe0k/HluC0hnzZF4HnAbUBf4MfBh37q3Ae0AzoCewAjgmx+3pDdQJPu8VyHRgRNsyIZDpkbi2bQCGAA2Bx4An4o5/HHgyqDs4aG/vHLdhInBhGd9TWtuStJdsPCJSFzgFOEdVj63oeMdxMkvUdFJEPgFuwALGj1TVQcH+BpjhP0BV54rIMuA8VZ0Q1N8IdFPV03MkeglEpAf2g/1rbKg8Mm0RkdOBU7GUcF1V9SwRuQXorKpnBMfsCXyOfU9FVCJMaaYRkYmYwf9Pwv60t6VSi5xVdYuqPh4FxYwqIlJDRB4NhnXGiEil1sxmExHpLCLfBUMkE4NlD6FFRJqIyNRgGGfvXMtTFaKkkyLSGuiOrdvsDXwcq1PVTcAXQO9gfrZdfH3wuXf2pC0dEblHRDYDc4Hl2LKayLRFLBXcaCAxykZiG77Agux3D8rOmIEJCMX3AdwqIqtEZLKIHBbsS3tbqk9eluhxCvClqg7FlPLUHMuTLO+o6mFB+S7XwlTAZuA44JlcC1JdCJbFPAqMVdW52HDYuoTD1gGNgjoS6mN1OUVVLwnkOAR4DthKtNpyI/CAqn6dsL+iNpRVl0t+D3QB2mPTEy8Fvcm0t8UNZnjZE5gVfP4IU8woMFhE3hORW0TCnWlWVbdHwKjnDWJh+x7G3vIvDXaXF2ZvY9x2Yl3OUdWdqjoJ2B34ORFpi4jsDxwJ/F8p1RW18ChH/QAADMNJREFUIXQhEVV1iqpuUNWtqjoWmAwMJwNtcYOZBkTkUhGZLiJbReTBhLrKepXNAQ4PPh+JOQqkjQzJvBzoik2ytyKNveIMyetkieDl6QEsf+YIVd0eVM3GohTFjmuAvSzOVtU12DO1X9ylKgyLmQNqEshMNNpyGOZwtUREVmDe1iNE5CN2bUMXoA4WDvGHkIhx1wrj96FYztj0tyWX3k35UjDDcDLwL+DBhLoyPbGANlg4s8TSJvjC/wq8BdwN/DXsMidcYzgwOgryYkEA9s71c5TPBbgX+BBomLC/ZfB9jcA8S2+npGfpbcA72AvjXpjRyZlnKfYieHrwrNXAoi1tAk6KSluA+oFexMpfsGmJltgc3npsRKsB8AglPUufCPStATCYHHvJYo5WRwf/75pYLOVNQI9MtCWnSpRvBbgp/sc8+CK2Ad3j9j0M3Jbida8HhoRdZqBx3OdbMY/N0Mobd7wbzAwWLGyfAluwobBYOTOoPxKbpy/EPE47x51bBwvBtx6LcnR5jtvSMjB6awOZPgUuiquPTFvi5LqeYFlJsH0GsCQwPC8AzePqmgPPB3VLgDNC8H1Mw4ZS12IvZUdlqi2VWlbilI6I3ATsrqojg+2kk1+Xcq022BvQTuBNVb0lAjKfgCnfZuAr4HxV3RFWeYNjX8WChy8G7lPVB9Mpr+M4+UPolypEnEp7lanqCmyuIdtUReaXgJcyIVQ5VMlzT1WHp10ix3HyEnf6ySyh9CqrgKjJHDV5HceJKG4wM0tUvMriiZrMUZPXcZyI4gYzDYhIzSA0WQ2ghojUFZGaapE+ngNGi0gDERmMedM9nEt5IXoyR01ex3HyDzeY6eEazCvuKuCs4PM1Qd0lQD0s0e7jhCfRbtRkjpq8juPkGe4l6ziO4zhJ4D1Mx3GcShIkGbiw4iNLPffnIrIyCP6/W7plSxciskhEjsy1HGHADabjOJFDRA4WkfdFZJ2IfB9kqeiXa7mSJQhC/zdgmKo2VNXVpRxzgVgy6g2BYX1FRHId6DztiCV5fizXciSDr8N0HCdSBKmpXsYCnj8F1MbCn23NpVwp0hoL51bqXLuIHArcgoXOmykizYEKA3FElOFYerTQ4z1Mx3GiRncAtfyfO1W1UFUnqOonACIyMuhx3h30QOeKyBGxk4M8qA+IyHIRWSYiN4lIjbj680XkcxFZIyLjRaRTXN1RwfXWicg/sJjPpSIidUTkThH5Jih3Bvu6A/OCw9aKyFulnN4P+EBVZwZt/V5Vx6rqhrg2PCSWf3axiFwjlg0GEbleRB6Jk6OziKgEOXWDYeQbg//RBhGZICIt4o4/O7jmahH5Y9LfCj8M314pIp8EyRAeEJHWIvJacK83xHKDxo4vAI4C/hd4vj8S3HetiEwTy50aGtxgOo4TNeYDO0VkrIgcG/8DHMcA4EugBXAd8FzQSwMYC+zAMuscAAwDLgQQkZOBP2DB/lsC72Ge1wRG5VnMO7sFlhx6cDly/hE4CAu9uB/QH7hGLWlxLFFxU1U9vJRzpwBHi8gNIjJYROok1N8NNMHyQB4KnAOcV44siZwRHN8K66FfEbSxF5bg4Gws4fVuWPqyVBiBGcHuWK/4Nex/2gKzOb+KO7Y/lvd3FXBu0KYOwX0vxrzhQ4MbTMdxIoWqrsey0ijwb+A7EXkxoTfyLXCnWs7TJ7Ee3XHBMccCv1HVTar6LZYX8vTgvJ8Bt6rq50Ec5FuA/YNe5nBgjqo+o5ae7E5gRTminoll7PlWLe/qDZghSqaN72FGuw/wCrBaRP4mIjWC3vBPgKvV8kAuwjIbJXXtgP+q6nxVLcSGtfcP9v8IeFlV31XVrcC1QFEK1wW4W1VXquoy7IVjiqrODK43DntJiXEcxcOx2zFD2TUYOZgRfNehwQ2m4ziRIzBoI1V1d2BvrDd0Z9why7TkmrnFwTGdgFrA8mDYby1wH9bTIqi/K67ue2zYtX1w/tdxMmj8dim0C+6bKEOybXwtSCDQHAvGMRLrCbfAeoWJ126f7LUpaeg3YzGZYzLHt3ETsItDUgWsjPtcWMp2w7jt+PnLh4HxwBPBEPYdgXNUaHCDWY0QkUNEZF7FR1bq2iNFZFImru045aGqcwlStMXtbi8i8fOLHYFvMGOwFWihqk2D0lhVY0OkXwM/i6trqqr1VPV9LH9lh9gFg+t3oGy+wQxwogyptq9IVd/EcuPuDazCemOJ114WfN6E5byM0SaF2yW2sT7W60s7YhmZ2gIfAQSjATeoai9gEHA8NtQcGtxg5oBgYrwwmARfG7jHXxybtM8UqvqeqvZIkCMU66tEpLaIrBKRhhUfvcu5bwfOD+tF5GMROSkTMjrhQET2EpHfisjuwXYH4KdYLsQYrYBfiUgtETkN6Am8qqrLgQnAX0WksYgUiMieYl6pYImurxaR3sG1mwTngw2N9haRUwMHml9RvjF6HLhGRFoG859/wpIYJ9PGk0TkdBFpJkZ/bK7yQ1XdiQ2j3iwijYLh4svjrj0LGCIiHUWkCXB1MvcMeAY4XmzZTm1gNHF2QkQOE5F0RbsZDvwvNhIgIkNFZJ9gyHk99lKwM033SgtuMHPHCaraCHtLvA34PfBAbkXKKUOAWaq6sRLn/hpoq6qNgVHAIyLSNq3SOWFiA+bUM0VENmGG8jPgt3HHTAG6Yb2xm4Efxa11PAcb0pwDrMGMRFsAVR0H3I4NC64PrntsULcKOA3T19X/3965hlhVhWH4+WAsDGxioBAOUnTB6E4/hKAgMBDFShDKSKMoIyqRSOpHBAZClF3oRxZZmKWVCl2wMktCwS5kWljgj7IMJJEKywGTjN5+fOuMe86MM/vo5JnkfWDDOWvvvfba57Lfvdd6v2+V+j8Zop2LgC+B7eRE09tKWR32AXOB70jxWAEslrSyrJ9HPkn+AGwGXiMnp0bSR8CqctytZAhOLUpKyXtKfXtKO3ZXNpkAfFa3vmFoDScZT34X+4Ed5ETdtW4wjht1ZrX2MuKzhO8Crmkpm0QOrl9U3p8MPEHOBL6XvPMdW9ZdTf6I7yfNDXuA2yp1TSMvBr1kN82C6n7l9avleH+SU2Q9QN5Bz2tp13ZgRo1zuhXYXHm/mPwjd5MJ058kL14/AveSho2uyvZPUWagJ2eqXwR8Wtq2luwWWkn+mbZQmcl+kM/xIDCp09+zl84srb9FLyP62b4ITBmBerrK9aC70+fUzuInzFGCpC9IEbyqFD1G2rIvI+3vDbJLp8l4UowawO3AsxV7/UvkOMw4csxjQJyXpDmkGF+rzDTyOGm3n93cJiIuLfXXDiouXVxLgUvILCZ/kHfKU8u5XA7MGGTXaaRgN5lFuv4awDnkXe0y0gCxgwwVqB733Yg4SD5ZbCTv7I0xI4ikOyStH4GqeoCHy/Xhf4MFc3TxM9BTzARzgfuUAcu9pL19VmXbQ6Rl/ZCk98knsYmVdRdExKmS9knaVvP47wDnxeG5JecAqyT9VXP/MeS4TQ8pxAdK+Q3AM5J2S9pHdmn1ERFnA2MkVQ1JyyTtLH+odcBOSRuUVv819LemI2k6MI4U3vWS2rXCG2OOE8pQm+c63Y52sWCOLhqkjf100uW2tWJv/6CUN/mtiEeTqjV8JikcP0XEpoi4os7BlXFSq4HZxYB0E+3NK3kuaX9/pEVk+1nVGWjFr8ZiNWnHmt5s/yFJ68iA7+vaaLc5gZD0sqQrO90Oc+JhwRwlRCaObpDjfr+SonChDlvbuyXVcpBK2iLpetIp+DYpgoNuOkjZcjLgejJwQFI7A/w7yOwh6yJiYqV8D/2zhbRa8Vu7Y4+VLrIb1xhjRgwLZocp1vbpwBvACknflO7EpcDTEXFG2a4REVNq1HdSRNwcEd3KbCT7ObI1ey+ZWquPIpD/kCadfk+XkTkoFw51fEmvk2mwNkREU7RWA/PLOZxGOoKbdY4ljTobhzu3wSghBlMjYmwJIZhNOm43HU19xhhzJCyYnWNtRPSS3ZMPkS7Rai7IB4Hvgc+LvX0Dh8coh2MOsKvsdxcVI08Lj5JxYr9HxIJK+SvAxQy0dE9gaBs9AJKWk/FbH0fEWaT4f0g6br8iu1//JoV8Mplk+mCtMxtIAAtJt/AvZIjJjW2M2xpjTC2iWHyN6SMibgHurI4DlSDxNZJqjYcOU/9U4HlJZ0bEEuBbSUuOtV5jjPkv8ROm6UdJhXU38EK1vDhcj0osS3fptIjoiogGGRLyVln9deW1McaMWvyEafooY6Rvkt2/M1tcuMdS7ynkmOL5pJnpPWC+RtlMBMYYMxQWTGOMMaYG7pI1xhhjamDBNMYYY2pgwTTGGGNqYME0xhhjamDBNMYYY2pgwTTGGGNq8C8ll/7dtd5F6AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "# Let us now use the checkAtmProfiles function to inspect the atmospheric profiles.\n", "planet.checkAtmProfiles()" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "16127.792366356383" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# Compute the scale height at the Venusian surface for illustration.\n", "# planet.density_int is the interpolation function created by planet object \n", "# when atmosphere model is loaded.\n", "planet.scaleHeight(0, planet.density_int)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Congratulations!** Your planet now has an atmosphere model. In the next example, we will compute aerocapture trajectories for a vehicle flying in the Venusian atmosphere." ] } ], "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.3" } }, "nbformat": 4, "nbformat_minor": 2 }