127 lines
20 KiB
Plaintext
127 lines
20 KiB
Plaintext
{
|
|
"cells": [
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": 1,
|
|
"metadata": {},
|
|
"outputs": [
|
|
{
|
|
"data": {
|
|
"text/plain": [
|
|
"Text(0.5, 1.0, 'vel_init')"
|
|
]
|
|
},
|
|
"execution_count": 1,
|
|
"metadata": {},
|
|
"output_type": "execute_result"
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAAEICAYAAABhxi57AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAYnElEQVR4nO3dfYxd9X3n8fdnxoYEmmDADfHaLHE3FhG72gTW4qGJWILTLLgIuipNnDYpICJCRBLINipQaZPdqCs1u9mmRKxwR1ACLcQlLgR2hQiUgJquNtYaQ5KCyYbwaGPA5hlCAjPz2T/OGXx9mblzx3PuOXeOPy/paO6558zvnOtrvv7yPb8H2SYiIuo30vQNRETsqxKAIyIakgAcEdGQBOCIiIYkAEdENCQBOCKiIQnAMZQknS3pH5u+j4hBSgCOiGhIAnBEREMSgGOgJF0saWPXe5dJ+qakgyRdJWmHpO2S/lTSaFP3GlG3BOAYtA3AWknvACgD7MeA64FvAePAe4GjgY8Cn27mNiPqlwAcA2X7MWAL8O/Lt04GfgE8AqwFLrL9qu1ngG8A6xq50YgGLGr6BmKfcD3wCeBa4PfL/SOAxcAOSVPnjQBPNHGDEU1IAI46fAf475JWUGTCJwAvAL8Cltoeb/DeIhqTEkQMnO2dwN3A1cAjtrfa3gHcThGY3ylpRNK/kPRvm7zXiDolAEddrgc+Uv6c8ofAfsADwPPARmBZ/bcW0QxlQvaIiGYkA46IaEgCcEREQxKAIyIakgAcEdGQWvsB76f9/TYOrPOSEbFAvczzu2z/+t7+/r/78IF+9rmJvs6958e/+p7tU/b2Wnur1gD8Ng7kOK2p85IRsUD9vTc+Np/f3/XcBJu+t6Kvcxcv+/nS+Vxrb9U7Ek5C++9f6yUjYoH65XwbMBOerOJOBqbWACwJjWa2wYgYPAOTDPc4h2TAETGcXp1/E5MkA95NgkWZ/yciBs+YN1KC6CDQopQgImLwDEykBNFBgv0W13rJiNh3pQbcSYI8hIuIGhiYGPLJxmoPwN4/GXBE1KPKCrCkL1KsWWjgJ8A5tn/Zcfxs4L8B28u3Lrd9Za82638iNpLRzxExeMaV1YAlLQe+ABxl+zVJN1CsX/itrlP/1vbn+m231gDsETG5X3pBRMTg2fBGtRWIRcDbJb0BHAA8WUWD9RrV7OdERMybmKDveLNU0uaO/THbY1M7trdL+jrwOPAacLvt26dp53clnQj8P+CLtnsuMlt7DXhycR7CRcTgGZjsPwPeZXv1TAclHQycAaykWFD2O5I+aftvOk77n8C3bf9K0meAa4CTe1209n7ATgYcETWZQwY8m49QLCi7E0DSjcBvAm8GYNvPdpx/JfBfZ2u0/l4Qi/MQLiIGrxiIUVkAfhw4XtIBFCWINUBnyQJJy8rVvgFOB7bO1mi9D+GAyWTAEVEDA2+4moTP9iZJG4EtwDhwLzAm6avAZtu3AF+QdHp5/Dng7Nnarb0EMZkMOCJqYMREhYv+2P4K8JWut7/ccfxS4NK5tJkacES01qSHO97UW4IQON2AI6IGFdeAB6L2h3ATi4f7DyQi2kJMVFQDHpTa81GPJABHxOAVK2IkAL/JgsnMxRMRNbDF6x7ugV/1P4RLBhwRNZlMDXg3CyaSAUdEDYqHcClB7Gm4/zwiojXyEG5PIr0gIqIWeQjXTSQDjojaTGQgxm4GJjMQIyJqYMQbQz7yq4GhyLVeMSL2UXkI103JgCOiHkYpQXRLBhwRdclDuA5OCSIiamKTbmh7EEwurnaZ0oiI6RQP4YY742tgMp66rxgR+6o8hOuQyXgioi5GmZB9D0oGHBH1SQbcSTC5X2rAETF4BiYrzPgkfRH4dNn0T4BzbP+y4/j+wLXAvwGeBT5u+9FebTYwGU8CcETUQZUtSSRpOfAF4Cjbr0m6AVgHfKvjtHOB522/V9I64GvAx3u1W3MGbJxeEBFRg2JZ+kp7QSwC3i7pDeAA4Mmu42cA/6l8vRG4XJJszxj0+grA06XewDJgA3AocA/wKduv926IZMARUQtbcylBLJW0uWN/zPbY7ra8XdLXgceB14Dbbd/e1cZy4Iny/HFJL1LEx10zXXTWANwj9V4LfMP2BknrKdLvK3o3BiyenO2SERGVmMNAjF22V890UNLBFBnuSuAF4DuSPmn7b+Zzf/2WILpT7x3AycDvl8evoUi9ZwnAZmRRAnBEDF4xH3Bl3dA+AjxieyeApBuB3wQ6A/B24HBgm6RFwEEUD+NmNGsAni71pig5vGB7vDxtG0X6/RaSzgPOA9j/Xe9k+btemO2SERE8Ou8WKl0R43HgeEkHUMTBNcDmrnNuAc4C/g9wJvD9XvVf6K8E8ZbUGzil37su6yhjAMv/5RKfsuyBfn81IvZh/3uev190Q6smA7a9SdJGYAswDtwLjEn6KrDZ9i3AVcBfS3oIeI6iVNtTPyWI6VLvDwJLJC0qs+AVFOl3T5JZPDI+22kREfNW9VwQtr8CfKXr7S93HP8l8HtzabOf/PzN1FuSKFLvB4C7KNJsKNLum+dy4YiIQZtkpK+tKbNe2fYmij5tWyi6oI1QlBQuBv5DmW4fSpF+R0QMhWI6SvW1NaWvXhAzpN4PA8fO5WICRkk/4IioRybjiYhoQDEbWibjiYioXTEUOQE4IqIByYC7mBFlJFxE1KPCkXADkQw4IlppqhfEMEsAjojWSgmiQ7qhRURdsiZcRERDDIwnA46IaEZKEBERTXBKEG+RbmgRUYeKJ2QfiGTAEdFayYAjIhpQ5YTsg5IAHBGtZMT4ZB7CvUmYUVIDjoh6DHsNeLj/eYiI2FsuShD9bLORdKSk+zq2lyRd1HXOSZJe7DjnyzM096aUICKilSpelPOnwAcAJI1SrIF50zSn/sD2af22mwAcEa01oIdwa4Cf235svg010A84c0FExOAZMdH/Q7ilkjZ37I/ZHpvh3HXAt2c4doKkHwFPAl+yfX+viyYDjojWmsNDuF22V892kqT9gNOBS6c5vAU4wvYrktYC3wVW9WovD+EiopVc4UO4DqcCW2w//dbr+SXbr5SvbwUWS1raq7EGpqNMN7SIqIerrwF/ghnKD5LeDTxt25KOpUhwn+3VWEoQEdFS1U7GI+lA4LeAz3S8dz6A7fXAmcBnJY0DrwHrbPd86JUAHBGtVWUGbPtV4NCu99Z3vL4cuHwubSYAR0Qr2TAxOdwj4WpfFXk001FGRE2GfShyMuCIaCUzkIdwlUoAjoiWyooYERGN6d0HoXm19wMeST/giKhJShAREQ0oekEM92DfBOCIaK2UIDoJRjMbWkTUJCWIiIgGGA19AO6rQCJpiaSNkh6UtFXSCZIOkXSHpJ+VPw8e9M1GRMyF+9ya0m+F+jLgNtvvA94PbAUuAe60vQq4s9yPiBgOBk+qr60ps5YgJB0EnAicDWD7deB1SWcAJ5WnXQPcDVzcs62sihwRNWpDCWIlsBO4WtK9kq4sp2U7zPaO8pyngMOm+2VJ50naLGnzy8+NV3PXERF9sPvbmtJPAF4EHANcYfto4FW6yg3lnJfTfgzbY7ZX2179jkPyzC8i6jE1F0Q/W1P6CcDbgG22N5X7GykC8tOSlgGUP58ZzC1GROwFA1Z/W0NmTUltPyXpCUlH2v4pxZLMD5TbWcCflT9v7ueCI5mOMiJq0paBGJ8HritXBH0YOIcie75B0rnAY8DHBnOLERF7o9keDv3oKwDbvg+YbsnmNZXeTURElYY8Ax7umSoiIvaWq3sIJ+lISfd1bC9JuqjrHEn6pqSHJP1Y0jGztdvAsvRD/k9SRLRHReGmfP71AQBJo8B24Kau004FVpXbccAV5c8ZJQOOiBZTn9ucrAF+bvuxrvfPAK514YfAkqmeYjNJx9yIaK/+O10tlbS5Y3/M9tgM564Dvj3N+8uBJzr2t5Xv7ZjmXKD2EkRWRY6Imkz1A+7PLtvTdTTYQ9kT7HTg0nnc2ZuSAUdEaw2gH/CpwBbbT09zbDtweMf+ivK9GaUGHBHtVf18lJ9g+vIDwC3AH5a9IY4HXuyYL2dayYAjor0qHGZcTkL2W8BnOt47H8D2euBWYC3wEPALigFrPdUegLMqckTUpcoV0Gy/Chza9d76jtcGLphLm8mAI6KdLGjDUOSIiAVpyMd9JQBHRHslAO+WfsARUasE4IiIBsxtIEYjEoAjorWq7AUxCAnAEdFeCcB7ynSUEVGXZMAREU1JDTgiogFzn+ehdrWviJGhyBFRmwTgiIhmDPuwgwTgiGivZMAREfWT0wuiS4YiR0SN0gsiIqIhyYAjIpqREkRERBOcXhB7EBmKHBE1GvJwk1WRI6K9KlwVWdISSRslPShpq6QTuo6fJOlFSfeV25dnazMliIhorYprwJcBt9k+U9J+wAHTnPMD26f122ACcETELCQdBJwInA1g+3Xg9fm2W28AFowMe1U8Itqj/wx4qaTNHftjtsc69lcCO4GrJb0fuAe4sFyqvtMJkn4EPAl8yfb9vS6aDDgi2mluvSB22V7d4/gi4Bjg87Y3SboMuAT4jx3nbAGOsP2KpLXAd4FVvS6ah3AR0V7VPYTbBmyzvanc30gRkHdfyn7J9ivl61uBxZKW9mo0ATgiWknsng9itm02tp8CnpB0ZPnWGuCBPa4nvVuSytfHUsTXZ3u123cJQtIosBnYbvs0SSuBDcChFPWQT5WF6ZnbwIxmPuCIqEu1vSA+D1xX9oB4GDhH0vkAttcDZwKflTQOvAass93zDuZSA74Q2Aq8s9z/GvAN2xskrQfOBa6Yy6eJiBiYimdDs30f0F0nXt9x/HLg8rm02VcJQtIK4LeBK8t9ASdT1EEArgF+Zy4XjogYuMk+t4b0WwP+C+CP2X2rhwIv2B4v97cBy6f7RUnnSdosafMLz6b8EBH1qaoGPCizliAknQY8Y/seSSfN9QJlX7oxgPf96/2duSAiojZDHm76qQF/EDi97Nf2Nooa8GXAEkmLyix4BbB9cLcZETFHC2BV5FlLELYvtb3C9nuAdcD3bf8BcBfFUz+As4CbB3aXERF7YcGXIHq4GNgg6U+Be4GrZvsFkaHIEVGjIc+A5xSAbd8N3F2+fhg4tvpbioioxrDne5kLIiLaaQHUgBOAI6KVVG7DrPYAnG5oEVGbIQ83yYAjorWyKnJERFMSgCMiGpBl6fdUTEc55P8kRUR7DHm4SQYcEa2VGnBERFMSgCMimpEMuMvIsP+JREQ7mEYnW+9HFuWMiFaqclFOAElLJG2U9KCkrZJO6DouSd+U9JCkH0s6Zqa2pqQEERHtVe3/cF8G3Gb7zHJhzgO6jp8KrCq34yjWyDyuV4M1d0PLUOSIqI96L0rcfzvSQcCJwNkA5Qrw3avAnwFcW66E/MMyY15me8dM7aYEERHt5DlssHRq7cpyO6+rtZXATuBqSfdKulLSgV3nLAee6Nifca3MKQnAEdFac6gB77K9umMb62pqEXAMcIXto4FXgUvme38JwBHRWprsb+vDNmCb7U3l/kaKgNxpO3B4x/6sa2XW3w0tNeCIqEtF4cb2U5KekHSk7Z8Ca4AHuk67BficpA0UD99e7FX/hfSCiIi2qn7Bzc8D15U9IB4GzpF0PoDt9cCtwFrgIeAXwDmzNZgAHBHtVWEAtn0fsLrr7fUdxw1cMJc2E4AjopWmBmIMs/r7AQ/7n0hEtIYmhzveJAOOiHbKqsgREc3JihgREU1JBrxbliSKiDoN+yOnZMAR0U4GKpqMZ1ASgCOitVID7pLJJyKiDukHHBHRFDsliIiIpiQDjohoSgLwbsVQ5DqvGBH7smTAERFNMDAx3BF41k4Jkg6XdJekByTdL+nC8v1DJN0h6Wflz4MHf7sREf2rcln6QeinV9g48Ee2jwKOBy6QdBTFekh32l4F3EkF6yNFRFRqqifEbFtDZi1BlEtq7ChfvyxpK8VKn2cAJ5WnXQPcDVw8W3vpBxwRdRn2GvCc4qGk9wBHA5uAwzrWO3oKOGyG3zlvaqnn554b8mEpEdEec1uWvhF9B2BJvwb8HXCR7Zc6j5VLcUz7MWyPTS31fMghyX8joh4CNOG+tqb01QtC0mKK4Hud7RvLt5+WtMz2DknLgGcGdZMREXtDFdZ3JT0KvAxMAOO2V3cdPwm4GXikfOtG21/t1easAViSgKuArbb/vOPQLcBZwJ+VP2+etS1gdLaTIiKqMJjywodt7+px/Ae2T+u3sX4y4A8CnwJ+Ium+8r0/oQi8N0g6F3gM+Fi/F42IGLwWzAVh+x8pktfprKn2diIiqjOHXhBLJW3u2B+zPdZ1joHbJRn4y2mOA5wg6UfAk8CXbN/f66L1joSTGFXGIkdETfrPgHd113Sn8SHb2yW9C7hD0oO2/6Hj+BbgCNuvSFoLfBdY1avBdEuIiHZytb0gbG8vfz4D3AQc23X8JduvlK9vBRZLWtqrzQTgiGivivoBSzpQ0jumXgMfBf6p65x3l50WkHQsRXx9tle7mYwnIlqrwm5ohwE3lfF1EXC97dsknQ9gez1wJvBZSePAa8C6cozEjGqfjjIpd0TUpqIAbPth4P3TvL++4/XlwOVzaTcZcES0k4Ehn/0gATgiWkm40pFwg5AAHBHtNTncKXDtAXh0xjEdEREVSgkiIqI5KUFERDQlAXg3ASMZihwRtWjBZDwREQvSAlgVOQE4IlorNeCIiKYkAO8p3dAiohYGJhOAIyIakIdwERHNSQCOiGiAgYnhHgpXcz9gMZIJKSOiFgYnAEdENCMliIiIBqQXREREg5IB7ynL0kdEbRKAIyIaYMPERGXNSXoUeBmYAMZtr+46LuAyYC3wC+Bs21t6tZkAHBHtVX0G/GHbu2Y4diqwqtyOA64of86ogVWR0w0tImpSbwniDODacin6H0paImmZ7R0z/UKiYUS0lIteEP1ssFTS5o7tvOkb5HZJ98xwfDnwRMf+tvK9GaUEERHtZHD/AzF2ddd0p/Eh29slvQu4Q9KDtv9hPreYDDgi2mtisr+tD7a3lz+fAW4Cju06ZTtweMf+ivK9GdWeAY9kOsqIqINd2bL0kg4ERmy/XL7+KPDVrtNuAT4naQPFw7cXe9V/ISWIiGiz6h7CHQbcVPQ0YxFwve3bJJ1fXMbrgVspuqA9RNEN7ZzZGk0AjojWckUZsO2HgfdP8/76jtcGLphLuwnAEdFSmZB9D0KMKs/9IqIGC2AynnlFQ0mnSPqppIckXVLVTUVEzJcBT0z0tTVlrwOwpFHgf1AMvzsK+ISko6q6sYiIeXE5IXs/W0PmkwEfCzxk+2HbrwMbKIbiRUQMBU+6r60p8wnAfQ27k3Te1PC+nc82l+pHxD5oyDPggT+Esz0GjAFI2jm67KFXgZlmE1rIltK+z9XGzwTt/Fxt/ExHzOeXX+b57/29Ny7t8/RG/uzmE4DnPOzO9q9L2tzHmOsFp42fq42fCdr5udr4mebL9ilN38Ns5lOC+L/AKkkrJe0HrKMYihcREX3Y6wzY9rikzwHfA0aBv7J9f2V3FhHRcvOqAdu+lWL881yMzeeaQ6yNn6uNnwna+bna+JlaTx7yoXoREW2VccEREQ1JAI6IaEitAbgNc0dIOlzSXZIekHS/pAvL9w+RdIekn5U/D276XudK0qikeyX9r3J/paRN5ff1t2VvlwWlXBhxo6QHJW2VdEJLvqsvln///knStyW9rQ3f176mtgDcorkjxoE/sn0UcDxwQfk5LgHutL0KuLPcX2guBLZ27H8N+Ibt9wLPA+c2clfzcxlwm+33UcznupUF/l1JWg58AVht+19R9EJaRzu+r31KnRlwK+aOsL3D9pby9csU/0Evp/gs15SnXQP8TiM3uJckrQB+G7iy3BdwMrCxPGUhfqaDgBOBqwBsv277BRb4d1VaBLxd0iLgAGAHC/z72hfVGYDnvGTzsJP0HuBoYBNwWMf6T09RLGGykPwF8MfA1MD4Q4EXbI+X+wvx+1oJ7ASuLksrV5breS3o76pcHPLrwOMUgfdF4B4W/ve1z8lDuL0k6deAvwMusv1S57FyaZIF079P0mnAM7bvafpeKrYIOAa4wvbRwKt0lRsW2ncFUNasz6D4B+afAQcCQz/sNt6qzgA857kjhpWkxRTB9zrbN5ZvPy1pWXl8GfBMU/e3Fz4InC7pUYrS0MkUtdMl5f/iwsL8vrYB22xvKvc3UgTkhfxdAXwEeMT2TttvADdSfIcL/fva59QZgFsxd0RZG70K2Gr7zzsO3QKcVb4+C7i57nvbW7Yvtb3C9nsovpfv2/4D4C7gzPK0BfWZAGw/BTwh6cjyrTXAAyzg76r0OHC8pAPKv49Tn2tBf1/7olpHwklaS1FrnJo74r/UdvGKSPoQ8APgJ+yul/4JRR34BuCfA48BH7P9XCM3OQ+STgK+ZPs0Sb9BkREfAtwLfNL2rxq8vTmT9AGKB4v7AQ9TLBU+wgL/riT9Z+DjFL1y7gU+TVHzXdDf174mQ5EjIhqSh3AREQ1JAI6IaEgCcEREQxKAIyIakgAcEdGQBOCIiIYkAEdENOT/AzTtGehVMKdgAAAAAElFTkSuQmCC",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 2 Axes>"
|
|
]
|
|
},
|
|
"metadata": {
|
|
"needs_background": "light"
|
|
},
|
|
"output_type": "display_data"
|
|
},
|
|
{
|
|
"data": {
|
|
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAWAAAAEICAYAAABhxi57AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAZJElEQVR4nO3dfbAd9X3f8ffnSgIbbPOkGKsSxSTW4KGd2lAND7FLMXJcUCikU2rLsR1g8GA82AY3HgOZ1m48yUycOnHwkEG9A+EhAROsgCEtxRAME6dTqxUC20HCtcyjhACJ5ycb7r2f/rErdHS4D+fcu2f33NXnNbNzz57d+9s9OuKrL9/9Pcg2ERFRv5GmbyAiYk+VABwR0ZAE4IiIhiQAR0Q0JAE4IqIhCcAREQ1JAI6IaEgCcNRC0pmS/qGH816S9Ks9ttnzuRHDaGHTNxDRyfbbZnOupKuALbb/0yDuK2IQkgFHRDQkATj6IulCSWu73rtE0rck7SfpCknbJG2V9AeSFvTZviW9p3x9laQ/l/Q/JL0oaZ2kX+s+V9I5wCeAL5dlib+t4rNGDFoCcPTremCVpLcDlAH2o8B1wFXAGPAe4EjgI8Cn53i91cDvAwcAm4E/7D7B9ihwLfDHtt9m+9/O8ZoRtUgAjr7YfgTYAPy78q0TgVeAh4BVwAW2X7b9FPBNigA6FzfZ/j+2xyiC7Pvn2F7E0MhDuJiN64CPA9cAv13uHwosArZJ2nneCPDYHK/1RMfrV4CeH9JFDLsE4JiN7wB/ImkZRSZ8HPAc8EtgcZmt1i3zqsa8kxJE9M32duBu4ErgIdubbG8DbqcIzO+QNCLp1yT965pu60kgfYJjXkkAjtm6Dvhw+XOn3wH2AjYCzwJrgSU13c8VwBGSnpP03ZquGTEnyooYERHNSAYcEdGQPISLWkn6V8D/nOxYP8OQI9ogJYiIiIbUmgHvpb39Fvat85IRMU+9yLM7bP/KbH//33xoXz/9zHhP597z419+z/ZJs73WbNUagN/CvhyjlXVeMiLmqb/z2kfm8vs7nhln3feW9XTuoiU/XzyXa81WvTVgCe29d62XjIh56hdzbcCMe6KKOxmYWgOwJLSgr8mxIiJmxcDEkA+QTAYcEcPp5bk3MUEy4F0kWJiebxExeMa8nhJEB4EWpgQREYNnYDwliA4S7LWo1ktGxJ4rNeBOEuQhXETUwMD4kA80qz0Ae+9kwBFRjyorwJK+SLHEloGfAGfZ/kXH8TOB/wpsLd+61Pbl07VZ/xOxkcz/ExGDZ1xZDVjSUuALwBG2X5V0A8VyW1d1nfrXtj/Xa7u1BmCPiIm90gsiIgbPhterrUAsBN4q6XVgH+DxKhqs1wLNfE5ExJyJcXqON4slre/YHy1X2wbA9lZJ3wAeBV4Fbrd9+yTt/HtJxwP/D/ii7WnXRKy9BjyxKA/hImLwDEz0ngHvsL1iqoOSDgBOAw6jWP/wO5I+afuvOk77W+Dbtn8p6TPA1RSrhk+p9n7ATgYcETXpIwOeyYcp1j/cDiDpRuDXgTcCsO2nO86/HPjjmRqtvxfEojyEi4jBKwZiVBaAHwWOlbQPRQliJdBZskDSknJxWoBTgU0zNVrvQzhgIhlwRNTAwOuuJuGzvU7SWmADMAbcC4xK+hqw3vYtwBcknVoefwY4c6Z2ay9BTCQDjogaGDFe4bKXtr8KfLXr7a90HL8YuLifNlMDjojWmvBwx5t6SxACpxtwRNSg4hrwQNT+EG580XD/gUREW4jximrAg1J7PuqRBOCIGLxiRYwE4DdYMJG5eCKiBrZ4zcM98Kv+h3DJgCOiJhOpAe9iwXgy4IioQfEQLiWI3Q33n0dEtEYewu1OpBdERNQiD+G6iWTAEVGb8QzE2MXARAZiREQNjHh9yEd+NTAUudYrRsQeKg/huikZcETUwygliG7JgCOiLnkI18EpQURETWzSDW03golF1S5TGhExmeIh3HBnfA1MxlP3FSNiT5WHcB0yGU9E1MUoE7LvRsmAI6I+yYA7CSb2Sg04IgbPwESFGZ+kLwKfLpv+CXCW7V90HN8buAb4l8DTwMdsPzxdmw1MxpMAHBF1UGVLEklaCnwBOML2q5JuAFYDV3WcdjbwrO33SFoNfB342HTt1pwBG6cXRETUoFiWvtJeEAuBt0p6HdgHeLzr+GnAfylfrwUulSTbUwa9ngLwZKk3sAS4HjgIuAf4lO3Xpm+IZMARUQtb/ZQgFkta37E/ant0V1veKukbwKPAq8Dttm/vamMp8Fh5/pik5yni446pLjpjAJ4m9V4FfNP29ZLWUKTfl03fGLBoYqZLRkRUoo+BGDtsr5jqoKQDKDLcw4DngO9I+qTtv5rL/fVaguhOvbcBJwK/XR6/miL1niEAm5GFCcARMXjFfMCVdUP7MPCQ7e0Akm4Efh3oDMBbgUOALZIWAvtRPIyb0owBeLLUm6Lk8JztsfK0LRTp95tIOgc4B2Dvd76Dpe98bqZLRkTw8JxbqHRFjEeBYyXtQxEHVwLru865BTgD+N/A6cD3p6v/Qm8liDel3sBJvd51WUcZBVj6z/b3SUs29vqrEbEH+19z/P2iG1o1GbDtdZLWAhuAMeBeYFTS14D1tm8BrgD+UtJm4BmKUu20eilBTJZ6fwDYX9LCMgteRpF+T0syi0bGZjotImLOqp4LwvZXga92vf2VjuO/AP5DP232kp+/kXpLEkXqvRG4iyLNhiLtvrmfC0dEDNoEIz1tTZnxyrbXUfRp20DRBW2EoqRwIfAfy3T7IIr0OyJiKBTTUaqnrSk99YKYIvV+EDi6n4sJWED6AUdEPTIZT0REA4rZ0DIZT0RE7YqhyAnAERENSAbcxYwoI+Eioh4VjoQbiGTAEdFKO3tBDLME4IhorZQgOqQbWkTUJWvCRUQ0xMBYMuCIiGakBBER0QSnBPEm6YYWEXWoeEL2gUgGHBGtlQw4IqIBVU7IPigJwBHRSkaMTeQh3BuEWUBqwBFRj2GvAQ/3Pw8REbPlogTRyzYTSYdLuq9je0HSBV3nnCDp+Y5zvjJFc29ICSIiWqniRTl/CrwfQNICijUwb5rk1B/YPqXXdhOAI6K1BvQQbiXwc9uPzLWhBvoBZy6IiBg8I8Z7fwi3WNL6jv1R26NTnLsa+PYUx46T9CPgceBLtu+f7qLJgCOitfp4CLfD9oqZTpK0F3AqcPEkhzcAh9p+SdIq4LvA8unay0O4iGglV/gQrsPJwAbbT775en7B9kvl61uBRZIWT9dYA9NRphtaRNTD1deAP84U5QdJ7wKetG1JR1MkuE9P11hKEBHRUtVOxiNpX+A3gM90vHcugO01wOnAZyWNAa8Cq21P+9ArATgiWqvKDNj2y8BBXe+t6Xh9KXBpP20mAEdEK9kwPjHcI+FqXxV5QaajjIiaDPtQ5GTAEdFKZiAP4SqVABwRLZUVMSIiGjN9H4Tm1d4PeCT9gCOiJilBREQ0oOgFMdyDfROAI6K1UoLoJFiQ2dAioiYpQURENMBo6ANwTwUSSftLWivpAUmbJB0n6UBJd0j6WfnzgEHfbEREP9zj1pReK9SXALfZfi/wPmATcBFwp+3lwJ3lfkTEcDB4Qj1tTZmxBCFpP+B44EwA268Br0k6DTihPO1q4G7gwmnbyqrIEVGjNpQgDgO2A1dKulfS5eW0bAfb3lae8wRw8GS/LOkcSeslrX/xmbFq7joiogd2b1tTegnAC4GjgMtsHwm8TFe5oZzzctKPYXvU9grbK95+YJ75RUQ9ds4F0cvWlF4C8BZgi+115f5aioD8pKQlAOXPpwZzixERs2DA6m1ryIwpqe0nJD0m6XDbP6VYknljuZ0B/FH58+ZeLjiS6SgjoiZtGYjxeeDackXQB4GzKLLnGySdDTwCfHQwtxgRMRvN9nDoRU8B2PZ9wGRLNq+s9G4iIqo05BnwcM9UERExW67uIZykwyXd17G9IOmCrnMk6VuSNkv6saSjZmq3gWXph/yfpIhoj4rCTfn86/0AkhYAW4Gbuk47GVhebscAl5U/p5QMOCJaTD1ufVkJ/Nz2I13vnwZc48IPgf139hSbSjrmRkR79d7parGk9R37o7ZHpzh3NfDtSd5fCjzWsb+lfG/bJOcCtZcgsipyRNRkZz/g3uywPVlHg92UPcFOBS6ew529IRlwRLTWAPoBnwxssP3kJMe2Aod07C8r35tSasAR0V7Vz0f5cSYvPwDcAvxO2RviWOD5jvlyJpUMOCLaq8JhxuUkZL8BfKbjvXMBbK8BbgVWAZuBVygGrE2r9gCcVZEjoi5VroBm+2XgoK731nS8NnBeP20mA46IdrKgDUORIyLmpSEf95UAHBHtlQC8S/oBR0StEoAjIhrQ30CMRiQAR0RrVdkLYhASgCOivRKAd5fpKCOiLsmAIyKakhpwREQD+p/noXa1r4iRocgRUZsE4IiIZgz7sIME4Ihor2TAERH1k9MLokuGIkdEjdILIiKiIcmAIyKakRJEREQTnF4QuxEZihwRNRrycJNVkSOivSpcFVnS/pLWSnpA0iZJx3UdP0HS85LuK7evzNRmShAR0VoV14AvAW6zfbqkvYB9JjnnB7ZP6bXBBOCIiBlI2g84HjgTwPZrwGtzbbfeACwYGfaqeES0R+8Z8GJJ6zv2R22PduwfBmwHrpT0PuAe4PxyqfpOx0n6EfA48CXb90930WTAEdFO/fWC2GF7xTTHFwJHAZ+3vU7SJcBFwH/uOGcDcKjtlyStAr4LLJ/uonkIFxHtVd1DuC3AFtvryv21FAF516XsF2y/VL6+FVgkafF0jSYAR0QriV3zQcy0zcT2E8Bjkg4v31oJbNztetK7JKl8fTRFfH16unZ7LkFIWgCsB7baPkXSYcD1wEEU9ZBPlYXpqdvALMh8wBFRl2p7QXweuLbsAfEgcJakcwFsrwFOBz4raQx4FVhte9o76KcGfD6wCXhHuf914Ju2r5e0BjgbuKyfTxMRMTAVz4Zm+z6gu068puP4pcCl/bTZUwlC0jLgN4HLy30BJ1LUQQCuBn6rnwtHRAzcRI9bQ3qtAf8Z8GV23epBwHO2x8r9LcDSyX5R0jmS1kta/9zTKT9ERH2qqgEPyowlCEmnAE/ZvkfSCf1eoOxLNwrw3n+xtzMXRETUZsjDTS814A8Ap5b92t5CUQO+BNhf0sIyC14GbB3cbUZE9GkerIo8YwnC9sW2l9l+N7Aa+L7tTwB3UTz1AzgDuHlgdxkRMQvzvgQxjQuB6yX9AXAvcMVMvyAyFDkiajTkGXBfAdj23cDd5esHgaOrv6WIiGoMe76XuSAiop3mQQ04ATgiWknlNsxqD8DphhYRtRnycJMMOCJaK6siR0Q0JQE4IqIBWZZ+d8V0lEP+T1JEtMeQh5tkwBHRWqkBR0Q0JQE4IqIZyYC7jAz7n0hEtINpdLL1XmRRzohopSoX5QSQtL+ktZIekLRJ0nFdxyXpW5I2S/qxpKOmamunlCAior2q/R/uS4DbbJ9eLsy5T9fxk4Hl5XYMxRqZx0zXYM3d0DIUOSLqo+kXJe69HWk/4HjgTIByBfjuVeBPA64pV0L+YZkxL7G9bap2U4KIiHZyHxss3rl2Zbmd09XaYcB24EpJ90q6XNK+XecsBR7r2J9yrcydEoAjorX6qAHvsL2iYxvtamohcBRwme0jgZeBi+Z6fwnAEdFamuht68EWYIvtdeX+WoqA3GkrcEjH/oxrZdbfDS014IioS0XhxvYTkh6TdLjtnwIrgY1dp90CfE7S9RQP356frv4L6QUREW1V/YKbnweuLXtAPAicJelcANtrgFuBVcBm4BXgrJkaTACOiPaqMADbvg9Y0fX2mo7jBs7rp80E4IhopZ0DMYZZ/f2Ah/1PJCJaQxPDHW+SAUdEO2VV5IiI5mRFjIiIpiQD3iVLEkVEnYb9kVMy4IhoJwMVTcYzKAnAEdFaqQF3yeQTEVGH9AOOiGiKnRJERERTkgFHRDQlAXiXYihynVeMiD1ZMuCIiCYYGB/uCDxjpwRJh0i6S9JGSfdLOr98/0BJd0j6WfnzgMHfbkRE76pcln4QeukVNgb8ru0jgGOB8yQdQbEe0p22lwN3UsH6SBERldrZE2KmrSEzliDKJTW2la9flLSJYqXP04ATytOuBu4GLpypvfQDjoi6DHsNuK94KOndwJHAOuDgjvWOngAOnuJ3ztm51PMzzwz5sJSIaI/+lqVvRM8BWNLbgL8BLrD9QuexcimOST+G7dGdSz0feGDy34iohwCNu6etKT31gpC0iCL4Xmv7xvLtJyUtsb1N0hLgqUHdZETEbKjC+q6kh4EXgXFgzPaKruMnADcDD5Vv3Wj7a9O1OWMAliTgCmCT7T/tOHQLcAbwR+XPm2dsC1gw00kREVUYTHnhQ7Z3THP8B7ZP6bWxXjLgDwCfAn4i6b7yvd+jCLw3SDobeAT4aK8XjYgYvBbMBWH7HyiS18msrPZ2IiKq00cviMWS1nfsj9oe7TrHwO2SDPy3SY4DHCfpR8DjwJds3z/dResdCSexQBmLHBE16T0D3tFd053EB21vlfRO4A5JD9j++47jG4BDbb8kaRXwXWD5dA2mW0JEtJOr7QVhe2v58yngJuDoruMv2H6pfH0rsEjS4unaTACOiPaqqB+wpH0lvX3na+AjwD92nfOustMCko6miK9PT9duJuOJiNaqsBvawcBNZXxdCFxn+zZJ5wLYXgOcDnxW0hjwKrC6HCMxpdqno0zKHRG1qSgA234QeN8k76/peH0pcGk/7SYDjoh2MjDksx8kAEdEKwlXOhJuEBKAI6K9JoY7Ba49AC+YckxHRESFUoKIiGhOShAREU1JAN5FwEiGIkdELVowGU9ExLw0D1ZFTgCOiNZKDTgioikJwLtLN7SIqIWBiQTgiIgG5CFcRERzEoAjIhpgYHy4h8LV3A9YjGRCyoiohcEJwBERzUgJIiKiAekFERHRoGTAu8uy9BFRmwTgiIgG2DA+Xllzkh4GXgTGgTHbK7qOC7gEWAW8Apxpe8N0bSYAR0R7VZ8Bf8j2jimOnQwsL7djgMvKn1NqYFXkdEOLiJrUW4I4DbimXIr+h5L2l7TE9rapfiHRMCJaykUviF42WCxpfcd2zuQNcruke6Y4vhR4rGN/S/nelFKCiIh2Mrj3gRg7umu6k/ig7a2S3gncIekB238/l1tMBhwR7TU+0dvWA9tby59PATcBR3edshU4pGN/WfnelGrPgEcyHWVE1MGubFl6SfsCI7ZfLF9/BPha12m3AJ+TdD3Fw7fnp6v/QkoQEdFm1T2EOxi4qehpxkLgOtu3STq3uIzXALdSdEHbTNEN7ayZGk0AjojWckUZsO0HgfdN8v6ajtcGzuun3QTgiGipTMi+GyEWKM/9IqIG82AynjlFQ0knSfqppM2SLqrqpiIi5sqAx8d72poy6wAsaQHw5xTD744APi7piKpuLCJiTlxOyN7L1pC5ZMBHA5ttP2j7NeB6iqF4ERFDwRPuaWvKXAJwT8PuJJ2zc3jf9qebS/UjYg805BnwwB/C2R4FRgEkbV+wZPPLwFSzCc1ni2nf52rjZ4J2fq42fqZD5/LLL/Ls9/7Oaxf3eHojf3ZzCcB9D7uz/SuS1vcw5nreaePnauNngnZ+rjZ+prmyfVLT9zCTuZQg/i+wXNJhkvYCVlMMxYuIiB7MOgO2PSbpc8D3gAXAX9i+v7I7i4houTnVgG3fSjH+uR+jc7nmEGvj52rjZ4J2fq42fqbWk4d8qF5ERFtlXHBEREMSgCMiGlJrAG7D3BGSDpF0l6SNku6XdH75/oGS7pD0s/LnAU3fa78kLZB0r6T/Xu4fJmld+X39ddnbZV4pF0ZcK+kBSZskHdeS7+qL5d+/f5T0bUlvacP3taepLQC3aO6IMeB3bR8BHAucV36Oi4A7bS8H7iz355vzgU0d+18Hvmn7PcCzwNmN3NXcXALcZvu9FPO5bmKef1eSlgJfAFbY/ucUvZBW047va49SZwbcirkjbG+zvaF8/SLFf9BLKT7L1eVpVwO/1cgNzpKkZcBvApeX+wJOBNaWp8zHz7QfcDxwBYDt12w/xzz/rkoLgbdKWgjsA2xjnn9fe6I6A3DfSzYPO0nvBo4E1gEHd6z/9ATFEibzyZ8BXwZ2Dow/CHjO9li5Px+/r8OA7cCVZWnl8nI9r3n9XZWLQ34DeJQi8D4P3MP8/772OHkIN0uS3gb8DXCB7Rc6j5VLk8yb/n2STgGesn1P0/dSsYXAUcBlto8EXqar3DDfviuAsmZ9GsU/MP8E2BcY+mG38WZ1BuC+544YVpIWUQTfa23fWL79pKQl5fElwFNN3d8sfAA4VdLDFKWhEylqp/uX/4sL8/P72gJssb2u3F9LEZDn83cF8GHgIdvbbb8O3EjxHc7372uPU2cAbsXcEWVt9Apgk+0/7Th0C3BG+foM4Oa67222bF9se5ntd1N8L9+3/QngLuD08rR59ZkAbD8BPCbp8PKtlcBG5vF3VXoUOFbSPuXfx52fa15/X3uiWkfCSVpFUWvcOXfEH9Z28YpI+iDwA+An7KqX/h5FHfgG4J8CjwAftf1MIzc5B5JOAL5k+xRJv0qRER8I3At80vYvG7y9vkl6P8WDxb2ABymWCh9hnn9Xkn4f+BhFr5x7gU9T1Hzn9fe1p8lQ5IiIhuQhXEREQxKAIyIakgAcEdGQBOCIiIYkAEdENCQBOCKiIQnAEREN+f92MlRP6De+wAAAAABJRU5ErkJggg==",
|
|
"text/plain": [
|
|
"<Figure size 432x288 with 2 Axes>"
|
|
]
|
|
},
|
|
"metadata": {
|
|
"needs_background": "light"
|
|
},
|
|
"output_type": "display_data"
|
|
}
|
|
],
|
|
"source": [
|
|
"# final model\n",
|
|
"fpath_final_model = './OUTPUT_FILES/final_model.h5'\n",
|
|
"fpath_true_model = \"./test_model_init_100-100-100.h5\"\n",
|
|
"\n",
|
|
"import h5py\n",
|
|
"\n",
|
|
"# read vel, xi, eta dataset \n",
|
|
"with h5py.File(fpath_final_model, 'r') as f:\n",
|
|
" vel = f['vel'][:]\n",
|
|
" xi = f['xi'][:]\n",
|
|
" eta = f['eta'][:]\n",
|
|
" #zeta = f['zeta'][:]\n",
|
|
"\n",
|
|
"# read true model\n",
|
|
"with h5py.File(fpath_true_model, 'r') as f:\n",
|
|
" vel_true = f['vel'][:]\n",
|
|
" xi_true = f['xi'][:]\n",
|
|
" eta_true = f['eta'][:]\n",
|
|
" #zeta_true = f['zeta'][:]\n",
|
|
"# plot\n",
|
|
"import matplotlib.pyplot as plt\n",
|
|
"import numpy as np\n",
|
|
"\n",
|
|
"# plot vel\n",
|
|
"plt.figure()\n",
|
|
"plt.imshow(vel[:,:,5], origin='lower', aspect='auto')\n",
|
|
"plt.colorbar()\n",
|
|
"plt.title('vel')\n",
|
|
"\n",
|
|
"# plot xi\n",
|
|
"#plt.figure()\n",
|
|
"#plt.imshow(xi[:,:,10], origin='lower', aspect='auto')\n",
|
|
"#plt.colorbar()\n",
|
|
"#plt.title('xi')\n",
|
|
"#\n",
|
|
"## plot eta\n",
|
|
"#plt.figure()\n",
|
|
"#plt.imshow(eta[:,:,10], origin='lower', aspect='auto')\n",
|
|
"#plt.colorbar()\n",
|
|
"#plt.title('eta')\n",
|
|
"\n",
|
|
"# plot vel_true\n",
|
|
"plt.figure()\n",
|
|
"plt.imshow(vel_true[:,:,5], origin='lower', aspect='auto')\n",
|
|
"plt.colorbar()\n",
|
|
"plt.title('vel_init')\n"
|
|
]
|
|
},
|
|
{
|
|
"cell_type": "code",
|
|
"execution_count": null,
|
|
"metadata": {},
|
|
"outputs": [],
|
|
"source": []
|
|
}
|
|
],
|
|
"metadata": {
|
|
"kernelspec": {
|
|
"display_name": "Python 3.9.1 64-bit ('3.9.1')",
|
|
"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.9.1"
|
|
},
|
|
"vscode": {
|
|
"interpreter": {
|
|
"hash": "02f83e1f4cd9619657a6845405e2dd67c4de23753ba48bca5dce2ebf57b3813a"
|
|
}
|
|
}
|
|
},
|
|
"nbformat": 4,
|
|
"nbformat_minor": 2
|
|
}
|