{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "## Model setup XBeach 1D\n", "\n", "Setup a 1D XBeach model with the Python toolbox." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "**no xbTools installation found in environment, adding parent path of notebook to see if it works\n" ] } ], "source": [ "# import default modules\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "from matplotlib import cm\n", "import sys\n", "import os\n", "\n", "# method to import xbtools with try routine\n", "try:\n", " import xbTools\n", "except ImportError:\n", " print('**no xbTools installation found in environment, adding parent path of notebook to see if it works')\n", " sys.path.append(os.path.abspath(os.path.join('..', '..', 'xbeach-toolbox')))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Import the toolbox." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "## import xbeach tools\n", "#sys.path.append(os.path.abspath(os.path.join('..' )))\n", "\n", "from xbTools.grid.creation import xgrid, ygrid\n", "from xbTools.grid.extension import seaward_extend\n", "from xbTools.xbeachtools import XBeachModelSetup\n", "from xbTools.general.wave_functions import offshore_depth" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Data\n", "Load and plot the bathymetry data" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0.5, 1.0, 'bathy')" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEWCAYAAABv+EDhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAkfElEQVR4nO3deXyU9bn38c+VQAIkYU/CTkD2RRQQQaSCWsHWrVas+9Jaejyny/O0drH2aNvT9pz21LbnnJ7ntLUKVtDU1gWlVosVFFlllx0kKxAgQAghe3I9f2Q4RgyQbbgzM9/36zUvZu6Ze+a6SDLfmd/M73ebuyMiIrEnLugCREQkGAoAEZEYpQAQEYlRCgARkRilABARiVEKABGRGKUAEAHMLNvMrm6l+3IzG9Ia9yUSTgoAkRYws6Vm9kDQdYg0hwJARCRGKQBEPnSJmW0zs2NmNtfMOphZNzNbZGaHQ9sXmVk/ADP7MTAN+LWZlZjZr+vd19Vmtju0z39bnUQzO2pmY0/dyMzSzKzMzFLPc68iCgCReu4EZgIXAMOA71H3NzIXGAgMAMqAXwO4+yPAMuDL7p7s7l+ud1/XAZcA44BbgZnuXgFkAnfVu93twJvufjiMfYk0SAEg8qFfu3ueux8Ffgzc7u5H3P0Fdy919xOh7Vc04r7+zd2L3D0XWAJcFNr+NHCHmZ3627sbeKZ12xBpnHZBFyDShuTVO58D9DGzTsAvgVlAt9B1KWYW7+41Z7mvgnrnS4FkAHdfbWYngSvM7AAwBHiltRoQaQoFgMiH+tc7PwDYD3wDGA5c6u4FZnYRsAGw0O2as5zu09QNAxUAf3b38mZXLNICCgCRD/2TmS2i7hX7d4E/AinUjfsXmVl34LHT9jkIDG7i4zwDbAZOUDcEJBIIfQYg8qFngb8Be0OnHwG/AjoChcAq4PXT9vkP4JbQt33+szEP4u75wHrq3j0sa5XKRZrBdEAYkfPPzJ4C9rv794KuRWKXhoBEzjMzywBuBi4OuBSJcRoCEjmPzOxfgC3Av7t7VtD1SGzTEJCISIzSOwARkRgVUZ8B9OzZ0zMyMpq178mTJ0lKSmrdggIQDX1EQw8QHX1EQw8QHX2Es4d169YVuvvH1puKqADIyMhg7dq1zdp36dKlTJ8+vXULCkA09BENPUB09BENPUB09BHOHswsp6HtGgISEYlRCgARkRilABARiVEKABGRGKUAEBGJUQoAEZEYpQAQEYlRCgARkTDJPVLKgtU5FJZUBF1KgyJqIpiISCT52Rs7WLT5AI8t3MpVI9O4e3IGU4f0wMzOvfN5oAAQEQmDyupa3t55mKtHpjM4NYkX1uXzxtaDDEtP5vNTB3HTxX3p0D4+0Bo1BCQiEgars45woqKa2y7pz3c/NZIVD1/Jz2ePIz4uju+8+D6X/dtbPP63nRwqDu6Q0HoHICISBm9uO0iH9nFMHdITgMR28dwyoR+fHd+X1VlHeerdLH69ZA+/efsDrruwD+M61Jz3GhUAIiKtzN15c/shLh+SSseEjw7zmBmTB/dg8uAe5Bw5ybwV2Tz/Xh4vVdbwlwMr+PzUQVwzuhfxceH/nEBDQCIirWz7gRPsKyrjmlHpZ73dwB5JPHb9aFZ+9ypuH5FAQXE5Dy5Yz4yfL+UPK7MpqwzvuwIFgIhIK1u87SBmMGNEWqNu37lDe2ZmtGfpQzP4zV3j6ZGcwKMLtzLtZ0uYtzyLiurwBIECQESklb25/SDjB3QjNSWxSfvFxxmzxvTmxQcv4/kvTWFIWhLff3UbV/78bdZkHW31OhUAIiKt6HhpFe/vO870YR87AFejmRmTBnXnuS9O5pkvTKJXlw707tKhFausow+BRURa0Y6CYgDG9OvS4vsyM6YNTWXa0OaHydnoHYCISCvaUXACgFG9OwdcybkpAEREWtGOgmK6dWpPWhPH/4OgABARaUXbDpxgRK/ObWa9n7NRAIiItJKaWmdXwQlG9E4JupRGUQCIiLSS3KOllFXVMDICxv9BASAi0mp2HKj7BtDIXgoAEZGYsr3gBHEGQ9OTgy6lURQAIiKtZPuBYgb1TAp8nf/GUgCIiLSSHQXFETP+DwoAEZFWcaK8iryjZQoAEZFYs+tg3QzgEb0i4yugoAAQEWkV2w6EAkDvAEREYsvqvUdITUmkTxhW7QwXBYCISAvV1DrLdhdyxbDUiFgC4hQFgIhIC23KL+J4WRVXtOAYAEFQAIiItNDbOw8TZ3D5kJ5Bl9IkCgARkRZauusw4/p3pVtSQtClNIkCQESkBY6erGRzfhHThzXuAPBtSWABYGb9zWyJmW03s61m9rWgahERaa5luw/jDlcMj6zxfwj2mMDVwDfcfb2ZpQDrzGyxu28LsCYRkSZ5e+dhunVqz9i+LT8G8PkW2DsAdz/g7utD508A24G+QdUjItJUtbXOO7sPM21oKvFxkfP1z1PM3YOuATPLAN4Bxrh78WnXzQHmAKSnp0/IzMxs1mOUlJSQnBwZS7SeTTT0EQ09QHT0EQ09QHB95BTX8NiKcr44NoGpfdu36L7C2cOMGTPWufvEj13h7oGegGRgHXDzuW47YcIEb64lS5Y0e9+2JBr6iIYe3KOjj2jowT24Pv5n6R4f+O1FXnC8rMX3Fc4egLXewHNqoN8CMrP2wAvAAnd/MchaRESa6t3dhQxPTyG9c+Qs/1BfkN8CMuBJYLu7/yKoOkREmqOssoY12UeZNjSyJn/VF+Q7gKnA3cCVZrYxdPpUgPWIiDTamuyjVFbXMi3Cln+oL7Cvgbr7u0DkfWwuIgIs23WYhPg4JmV0D7qUZtNMYBGRZnh3TyGXDOpGx4TIOP5vQxQAIiJNdKi4nB0FJ5g2NHKHf0ABICLSZG/vOgwQ0R8AgwJARKTJXtm0n75dOzKyV+Qc/rEhCgARkSYoOF7O8j2FfHZ8X+IicPmH+hQAIiJNsHDjPmodPjO+X9CltJgCQESkkdydF9bnM35AVwb1TAq6nBZTAIiINNLW/cXsOljCzVHw6h8UACIijfbC+nwS4uO4/sI+QZfSKhQAIiKNUF5Vw8KN+7l6VBpdOrVs6ee2QgEgItIIL2/Yx9GTldw9OSPoUlqNAkBE5Bzcnd+/m8XoPp2ZPDhy1/45nQJAROQc3t51mD2HSvjC5YOoW8k+OigARETO4cl3s0hLSeS6KPnw9xQFgIjIWewoKGbZ7kLuvSyDhHbR9ZQZXd2IiLSyecuz6dA+jjsmDQi6lFanABAROYNjJyt5acM+PnNxP7olJQRdTqtTAIiInEHme3lUVNdy32UZQZcSFgoAEZEGVNfU8szKbC67oAfDe6UEXU5YKABERBqweNtB9h8vj9pX/6AAEBFp0NwV2fTr1pGrRqYHXUrYKABERE6z6+AJ1mQd5a7JA4mP8IO+nI0CQETkNM+uziUhPo7ZE6Jj2eczUQCIiNRTXlXDi+vzmTmmFz2SE4MuJ6wUACIi9SzafIDi8uqonPh1OgWAiEg9z67OYXBqUlSt+nkmCgARkZBNeUWszy3ijkkDomrVzzNRAIiIAPnHSpnzzFrSUhK5Jco//D1FASAiMe9ISQX3PLmGssoa/vCFSXTtFH3r/jSkXdAFiIgEqabWmfPMOvYVlbHggUsZ0atz0CWdNwoAEYlpTyzby7qcY/zqcxcxMSP6P/itL9AhIDObZWY7zWyPmX0nyFpEJPbsOXSCXyzexczR6dx4UXQd7asxAgsAM4sH/hu4FhgF3G5mo4KqR0RiS3VNLd/402aSEuL50U1jY+JbP6cL8h3AJGCPu+9190ogE7gxwHpEJEYcL6viC0+vZVNeET+8cQypKdE94/dMgvwMoC+QV+9yPnBpQLWISIzYc6iEOX9YS96xUv715rFcPy72hn5OMXcP5oHNZgMz3f2B0OW7gUnu/pXTbjcHmAOQnp4+ITMzs1mPV1JSQnJycsuKbgOioY9o6AGio49o6AEa38fyfVX8YVslCfHwlYs7MKxb/HmornHC+bOYMWPGOnef+LEr3D2QEzAFeKPe5YeBh8+2z4QJE7y5lixZ0ux925Jo6CMaenCPjj6ioQf3c/dRWlHtX//jRh/47UU++zcrfH9R6fkprAnC+bMA1noDz6lBDgG9Bww1s0HAPuA24I4A6xGRKHSouJwH/rCW9/cd52tXDeWrVw2N6jX+myKwAHD3ajP7MvAGEA885e5bg6pHRKLP1v3H+eLTaykqq+KJuydy9ajoPbpXcwQ6EczdXwNeC7IGEYk+pZXV/Mffd/Pksix6Jify/JemMKZvl6DLanM0E1hEosr+ojJm/2Yl+4rKuHViP75z7Ui6J8XG2j5NpQAQkajh7nz3pfc5erKS5780hUmDYmtph6bSaqAiEjVe2rCPpTsP861Zw/Xk3wgKABGJCodPVPDDRdsYP6Ar90zJCLqciKAAEJGo8MNF2yitqOFnt1yor3k2kgJARCLe9iM1vLppPw9Ov4AhaSlBlxMxFAAiEtGqamqZv72C/t078uD0C4IuJ6LoW0AiEtH+sDKHfSXOE/eMpkP7trO2TyTQOwARiVj5x0r51eJdXNgznqtHpgVdTsRRAIhIRCqvquEfF6wH4K5RCTF5QJeWUgCISET6wavb2Jx/nMdvHUdaJz2VNYf+10QkohSVVvKz13fw3JpcHpx+AdeM7hV0SRHrrB8Cm1nxOfY34IC7D2u9kkREPi7nyEmefDeLP63Np6yqhuvH9eEbn9RTT0uc61tAH7j7xWe7gZltaMV6REQ+Yl3OUX73zl7+tu0g7eKMGy/qywPTBjGiV+egS4t45wqAzzbiPhpzGxGRRqupdd7YWsATy/ayIbeILh3b84/TL+DeKRmkde4QdHlR46wB4O576182s87193H3o6ffRkSkuU5WVPP82jyeWp5F3tEyBvboxA9vHM0tE/rRKUHTllpbo/5HzexLwA+BMuDUUeQdGBymukQkhhwsLmfeimwWrMqhuLyaCQO78cinRvLJUb20rk8YNTZSHwJGu3thOIsRkdiy/UAxTyzby6ub9lNT68wa04svXD6YCQO7BV1aTGhsAHwAlIazEBGJDe7O27sO8/tlWby7p5BOCfHceelAPj91EAN6dAq6vJjS2AB4GFhhZquBilMb3f2rYalKRKJORXUNCzfu58llWew8eIL0zol8a9Zw7pw0kC6d2gddXkxqbAD8FngLeB+oDV85IhJtjp2sZMHqHJ5emcPhExWM6JXC47PHcf24PiS001zUIDU2AKrd/ethrUREokp24UmeWv7hxK1PDEvll7cOZuqQHlq3p41obAAsMbM5wKt8dAjoaFiqEpGI5O6syznGE8s0cSsSNDYA7gj9+3C9bfoaqIgAUF1TyxtbD/LEsr1szNPErUjRqABw90HhLkREIs/Jimr+tDaPJ0MTtwZ078QPbhjN7ImauBUJzrUY3Hh3X9/S24hIdNHErehwroiea2bTqVv180yeBM66YJyIRIfTJ27NHN2LB6Zp4lakOlcAdAHWcfYAONx65YhIW+PuvLO7kN8v28uy3Zq4FU3OtRhcxnmqQ0TamNMnbqWlJPLNmcO561JN3IoW+pRGRD6iqLSSBatzmbci+38nbv189jhu0MStqKMAEBHg40fcmja0J7+4dRyXD+mpiVtRSgEgEuPW5RzliXeyeGNbgSZuxZjGHg/g78Dj7v5avW2/c/c5zXlQM/t34HqgkrqVRu9396Lm3JeINN2pI249vrKMD15fqYlbMaqx7wAGAd82s0vc/QehbRNb8LiLgYfdvdrMfkrdDONvt+D+RKQRTj/iVmpH0xG3Ylhjf+JFwFXAf5rZq8BdLXlQd/9bvYurgFtacn8icnZnmriVcHgHV07JCLo8CYi5+7lvZLbB3S8Onb8P+AbQzd37tbiAukD5o7vPP8P1c4A5AOnp6RMyMzOb9TglJSUkJyc3u862Ihr6iIYeIDL6yDtRy+tZVaw6UE2tw4T0eGYNas+QrvFAZPTQGNHQRzh7mDFjxjp3//iojbuf8wR86bTLE4CnzrHPm8CWBk431rvNI8BLhILoXKcJEyZ4cy1ZsqTZ+7Yl0dBHNPTg3nb7qK2t9aU7D/ldv1/lA7+9yEf+81/9sYVbPKfw5Mdu21Z7aKpo6COcPQBrvYHn1MYuBvfb0y6vAz5/jn2uPtv1ZnYvcB1wVahAEWmBssoaXtyQz7zl2ew+VEJaio64JWcXyKc+ZjaLug99r3B3HWtYpAXyj5XyzKocMtfkcbysitF9OmviljRKUB/7/xpIBBaHJpiscvd/CKgWkYjj7ryXfYy5y7N4Y2sBADNH9+K+yzKYNKi7Jm5JowQSAO4+JIjHFYl05VU1vLppP/NWZLN1fzFdOrbni58YzD1TMujbtWPQ5UmE0Rd/RSJAwfFy5q/K4bk1uRw5Wcnw9BR+8pmxfObivnRMiA+6PIlQCgCRNsrdWZ9bxLwV2fz1/QPUuHPViHQ+PzWDKRfowOrScgoAkTamorqG194/wNzl2WzOP05Kh3bcd1kG90zJ0Pr70qoUACJtxKET5Ty7Opf5q3IpLKlgcGoS/3LjaG4e34+kRP2pSuvTb5VIwDbnFzF3eTaLNu+nqsaZPjyV+6cOYtqQnsTp+LoSRgoAkQBU1dTy+pYC5i7PYn1uEUmhwyzeM2Ugg1Mje0kDiRwKAJHz6EhJBZnv5fHMyhwKissZ2KMTj143itkT+5HSQbN15fxSAIicB9v2FzNvRRYvb9xPZXUt04b25Cc3j2H6sDQN80hgFAAiYVJdU8vibQeZuyKbNVlH6dg+nlsn9uO+yzIYkpYSdHkiCgCR1lZUWvm/wzz7isro160jj3xqJLdO7K9F2aRNUQCItJKdBSeYtyKblzbkU15Vy+TB3Xn0+lFcPTKdeA3zSBukABBpgZpa560dh5i7PIsVHxwhsV0cn7m4L/delsHI3jqourRtCgCRZigur+KN7CoeXbOU3KOl9O7SgW/OHM7tkwbQPSkh6PJEGkUBINIEHxwu4ekV2fx5XT6llTVckpHEt2eNYObodNrFa+19iSwKAJFzqK113t59mHnLs3l712ES4uO4flwfxiYWct+NlwVdnkizKQBEzqCkopoX1uXz9Ips9haeJC0lka9/chh3XDqAnsmJLF26NOgSRVpEASBymuzCkzy9Mps/rc2npKKai/p35T9uu4hrx/TWIRYlqigARKhbe//dPYXMW57NWzsPEW/Gp8b25v6pGVw8oFvQ5YmEhQJAYlppZTUvrt/HvBXZ7DlUQo+kBL4yYwh3Th5IeucOQZcnElYKAIlJeUdLeWZVDplrcikur2ZM3848Pnsc143rTWI7HWJRYoMCQGKGu7M66yhzl2exeNtBzIxZo3tx/9QMJgzspkMsSsxRAEjUK6+qYeHGfcxdns2OghN069Sef7jiAu6eMpDeXToGXZ5IYBQAErUOHC/jmZU5PLcml2OlVYzolcJPPzuWGy/qS4f2GuYRUQBIVHF31uUcY+7ybF7fWoC7c/XIdO6fOojJg7trmEekHgWARIWK6hpe3XSAeSuy2LKvmM4d2vGFywdx9+SB9O/eKejyRNokBYBEtEPF5cxflcOza3IpLKlkSFoyP7ppDDeP70unBP16i5yN/kIkIm3MK2Lu8iz+svkANe5cOTyN+6ZmcPmQnhrmEWkkBYBEjMrqWv665QBzl2ezMa+IlMR23DMlg3umDCSjZ1LQ5YlEHAWAtHmFJRU8tzqXZ1blcOhEBYN7JvH960dxy8T+JCfqV1ikufTXI23Wln3Hmbs8m1c37aeyppYrhqXy01syuGJoKnE6xKJIiykApE2prqnlja0Hmbcii/eyj9EpIZ7bJvXn3ssyuCA1OejyRKKKAkDahGMnK3nuvVyeWZnDgePlDOjeie99eiSzJ/anS8f2QZcnEpUCDQAzewj4dyDV3QuDrEWCsf1AMU+vyOalDfuoqK5l6pAe/PDGMVw5Io14DfOIhFVgAWBm/YFPArlB1SDBqHXn9S0FzFuRxaq9R+nQPo6bx/flvssGMbxXStDlicSMIN8B/BL4FrAwwBrkPDpeWsXza/P47TtlFJato2/Xjjx87Qg+d0l/unZKCLo8kZhj7n7+H9TsBuAqd/+amWUDE880BGRmc4A5AOnp6RMyMzOb9ZglJSUkJ0f+h4iR2Mf+kloW51SxfH81lTUwpLMza3AHLk6Lj+hhnkj8WZwuGnqA6OgjnD3MmDFjnbtPPH172ALAzN4EejVw1SPAd4Fr3P34uQKgvokTJ/ratWubVc/SpUuZPn16s/ZtSyKpjzVZR/mvt3azbHchCe3iuOmiPtx7WQaHd22ImB7OJpJ+FmcSDT1AdPQRzh7MrMEACNsQkLtffYZCxgKDgE2hKfv9gPVmNsndC8JVj5w/B4vL+dfXtvPyxv2kpSTyzZnDue2S/vRITgRg6a6ACxQRIIDPANz9fSDt1OWmvAOQtm/hxn1876UtVFTX8pUrh/CP04fQMUFr74u0RZoHIK2ipKKaRxdu4cX1+5gwsBuPzx6n9XlE2rjAA8DdM4KuQVpmU14RX83cQN7RUr521VC+cuUQ2sXHBV2WiJxD4AEgkau21vntO3t5/G87Se/cgT9+aQqXZHQPuiwRaSQFgDTLweJy/u8fN7LigyN8emxvfvKZsXTppCUbRCKJAkCa7M1tB/nmnzdRXlXLTz87llsn9tdBWEQikAJAGq28qoYf/2U7z6zKYXSfzvzn7RdrhU6RCKYAkEbZUVDMV5/bwK6DJXxx2iAemjmcxHb6eqdIJFMAyFm5O39YmcOPX9tO5w7tefrzk7hiWGrQZYlIK1AAyBkdPVnJt/68iTe3H2L68FR+PnscPUOzeUUk8ikApEHv7i7k689vpKi0ikevG8X9UzP0Qa9IlFEAyEdUVtfy+OKd/O6dvQzumcTc+y9hdJ8uQZclImGgAJD/lV14kq9mbmBz/nFunzSAf75uJJ0S9CsiEq301y24Oy+s38djC7cQH2f8z53juXZs76DLEpEwUwDEuCMlFTz2ylYWbT7ApYO688vPXUSfrh2DLktEzgMFQIxyd17dfIDvv7KVE+VVPHTNMB6cPiSij9AlIk2jAIgxx8uqeHF9PgtW57LnUAnj+nXhZ7dM1sHYRWKQAiBGvJ9/nPmrcli4aR/lVbWM69+Vn88ex00X9dHSzSIxSgEQxcoqa3h1034WrM5hU/5xOraP56aL+nLnpQMZ209f7RSJdQqAKLTnUAkLVufwwrp8isurGZqWzA9uGM1nxvelcwct2SwidRQAUaKyupbF2w4yf1UOK/ceoX28MWtMb+68dACXDuquWbwi8jEKgAi3r6iM51bnkvleHoUlFfTr1pFvzhzOrRP7k5qidXtE5MwUABGottZ5e9dh5q/KYcnOQzhw5fA07po8kE8MS9VXOUWkURQAEaSwpIJFeyv53uol5B8ro2dyAg9Ov4DbJw2gX7dOQZcnIhFGAdDGuTvvZR9j/qoc/rrlAFU1zuTBKXzn2hFcM6oXCe30FU4RaR4FQBtVXF7Fyxv2MX9VDrsOlpDSoR13XjqQYXEHueO6KUGXJyJRQAHQxmzZd5wFq3NYuHE/pZU1XNivCz/97FiuH9eHTgntWLr0cNAlikiUUAC0AeVVdRO25q/OZVNeER3ax3HDuD7cNXkgF/brGnR5IhKlFAAB2nu4hAWrc/nzunyOl1VxQWoSj143is9O6EeXjpqwJSLhpQA4z6pqanlz20Hmr85h+Z4jtIszZo7pxV2XDmTyYE3YEpHzRwFwnuwvKiPzvTwy1+Ry6EQFfbt25KFrhnHrJf1JS+kQdHkiEoMUAGFUW+u8s/swC1bn8vftB3Fg+rBUfnLpQGaMSNOELREJlAIgDI6erOT5tXk8uzqX3KOl9EhK4EtXXMAdkwbQv7smbIlI26AAaCXuzrqcuglbr71fQGVNLZMGdecb1wxj1pheJLaLD7pEEZGPUAC0UElFNS9t2MeCVTnsKDhBSmI7bp/UnzsnD2RYuo6yJSJtV2ABYGZfAb4MVAN/cfdvBVVLc2zbX8yC1Tm8vGEfJytrGN2nM/9681huGNeHpETlqoi0fYE8U5nZDOBG4EJ3rzCztCDqaKryqhpee/8A81flsD63iMR2cVwfmrA1rl8XfYVTRCJKUC9VHwT+zd0rANz9UEB1NEp24UkWrM7hT+vyKSqtYnDPJL736ZHcMqEfXTslBF2eiEizmLuf/wc12wgsBGYB5cBD7v7eGW47B5gDkJ6ePiEzM7NZj1lSUkJycnKjb19T62w8XMNbuVVsPVJLvMHFafFcOaA9I7vHBfZqv6l9tEXR0ANERx/R0ANERx/h7GHGjBnr3H3i6dvD9g7AzN4EejVw1SOhx+0GTAYuAZ43s8HeQBq5+++A3wFMnDjRp0+f3qx6li5dSmP2LTheTuZ7uWSuyaOguILeXTrw9U8O4LZL+pPWOfgJW43toy2Lhh4gOvqIhh4gOvoIooewBYC7X32m68zsQeDF0BP+GjOrBXoCgSx1WVvrLP+gkAWrclm8/SC17nxiaCr/ctMYZgxPpV281twXkegT1GcALwNXAkvNbBiQABSG68H+6++7eW5lKUnr327w+uLyKg4WV9A9KYEHpg3izkkDGdBDE7ZEJLoFFQBPAU+Z2RagEri3oeGf1pKakkifpDjS0hoeX4uPi+PqkWmasCUiMSWQAHD3SuCu8/V4t00aQK/SvUyfPuF8PaSISJunwW0RkRilABARiVEKABGRGKUAEBGJUQoAEZEYpQAQEYlRCgARkRilABARiVGBrAbaXGZ2GMhp5u49CeNyE+dRNPQRDT1AdPQRDT1AdPQRzh4Gunvq6RsjKgBawszWNrQcaqSJhj6ioQeIjj6ioQeIjj6C6EFDQCIiMUoBICISo2IpAH4XdAGtJBr6iIYeIDr6iIYeIDr6OO89xMxnACIi8lGx9A5ARETqUQCIiMSomAgAM5tlZjvNbI+ZfSfoes7EzJ4ys0OhI6Wd2tbdzBab2e7Qv93qXfdwqKedZjYzmKo/ysz6m9kSM9tuZlvN7Guh7ZHWRwczW2Nmm0J9/CC0PaL6ADCzeDPbYGaLQpcjsYdsM3vfzDaa2drQtojqw8y6mtmfzWxH6O9jSuA9uHtUn4B44ANgMHXHHt4EjAq6rjPU+glgPLCl3rafAd8Jnf8O8NPQ+VGhXhKBQaEe49tAD72B8aHzKcCuUK2R1ocByaHz7YHVwORI6yNU29eBZ4FFkfg7FaotG+h52raI6gN4GnggdD4B6Bp0D7HwDmASsMfd93rdoSgzgRsDrqlB7v4OcPS0zTdS94tD6N+b6m3PdPcKd88C9lDXa6Dc/YC7rw+dPwFsB/oSeX24u5eELrYPnZwI68PM+gGfBn5fb3NE9XAWEdOHmXWm7gXek1B3WFx3LyLgHmIhAPoCefUu54e2RYp0dz8AdU+uQFpoe5vvy8wygIupe/UccX2Ehk42AoeAxe4eiX38CvgWUFtvW6T1AHXh+zczW2dmc0LbIqmPwcBhYG5oOO73ZpZEwD3EQgBYA9ui4buvbbovM0sGXgD+j7sXn+2mDWxrE324e427XwT0AyaZ2Ziz3LzN9WFm1wGH3H1dY3dpYFub+FkAU919PHAt8E9m9omz3LYt9tGOuuHd/3H3i4GT1A35nMl56SEWAiAf6F/vcj9gf0C1NMdBM+sNEPr3UGh7m+3LzNpT9+S/wN1fDG2OuD5OCb1VXwrMIrL6mArcYGbZ1A19Xmlm84msHgBw9/2hfw8BL1E3HBJJfeQD+aF3kQB/pi4QAu0hFgLgPWComQ0yswTgNuCVgGtqileAe0Pn7wUW1tt+m5klmtkgYCiwJoD6PsLMjLpxzu3u/ot6V0VaH6lm1jV0viNwNbCDCOrD3R92937unkHd7/1b7n4XEdQDgJklmVnKqfPANcAWIqgPdy8A8sxseGjTVcA2gu4h6E/Gz8cJ+BR130b5AHgk6HrOUudzwAGgirpXAF8AegB/B3aH/u1e7/aPhHraCVwbdP2hmi6n7q3qZmBj6PSpCOzjQmBDqI8twKOh7RHVR73apvPht4Aiqgfqxs83hU5bT/0NR2AfFwFrQ79TLwPdgu5BS0GIiMSoWBgCEhGRBigARERilAJARCRGKQBERGKUAkBEJEYpAEREYpQCQKQZzCzDzMpCawU1Zb/PhZb4XRSm0kQaTQEg0nwfeN1aQY3m7n8EHghPOSJNowAQOY2ZXWJmm0MHhUkKHRDmbAvBnXpHsCO0yuMWM1tgZleb2fLQwT7a8rLKEqPaBV2ASFvj7u+Z2SvAj4COwHx333KO3QCGALOBOdStQXUHdUtj3AB8lw/XehdpExQAIg37IXVP4uXAVxu5T5a7vw9gZluBv7u7m9n7QEZYqhRpAQ0BiTSsO5BM3WEtOzRyn4p652vrXa5FL7akDVIAiDTsd8A/AwuAnwZci0hY6FWJyGnM7B6g2t2fNbN4YIWZXenubwVdm0hr0nLQIs0QOt7xInc/67eDzrDvdOAhd7+ulcsSaRINAYk0Tw3QpTkTwYD/BxwLR1EiTaF3ACIiMUrvAEREYpQCQEQkRikARERilAJARCRG/X8ALY3Xyn9oIgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "## load data\n", "bathy = np.loadtxt('..//examples//data//input//bathy.dep')\n", "bathy = bathy[0,:]\n", "\n", "## set bathy grid\n", "nx = 124\n", "ny = 72\n", "dx = 5\n", "dy = 20\n", "\n", "x = np.linspace(0,(nx-1)*dx,nx)\n", "\n", "## plot\n", "plt.figure()\n", "plt.plot(x,bathy)\n", "plt.xlabel('x [m]')\n", "plt.ylabel('z [m]')\n", "plt.grid('on')\n", "plt.title('bathy')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Create x-grid\n", "Create spatially varying x-grid resolution. In this example a minimum grid resolution of 2 m is applied.\n", "\n", "Note that the following settings can be specified:\n", "\n", "`xgrid(x,z,\n", " ppwl=20,\n", " dxmin=5,\n", " dxmax=np.inf,\n", " vardx=1,\n", " wl = 0,\n", " eps = 0.01,\n", " Tm = 8,\n", " xdry=None,\n", " zdry=None,\n", " dxdry = None,\n", " depthfac = 2,\n", " maxfac = 1.15,\n", " nonh = False)`" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEGCAYAAABsLkJ6AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAtiUlEQVR4nO3deXzU1b3/8dfn+50JQUBQtFiWFmwvEFDZIuJWQdxQFIsiXlDRXi9X29vNVquo7bUbtjzotb96bcutVVsXFMVIrS3iEsUVgYCAARSLSrCIegGjhMxyfn/MwmSykJCE70zm/Xw8eJB8M/Odc1jO53zP8jnmnENERAqPF3QBREQkGAoAIiIFSgFARKRAKQCIiBQoBQARkQIVCroALXHYYYe5/v37t/o+n376KV26dGl9gQKkOuQG1SE3qA5NW7FixYfOucOzr+dVAOjfvz/Lly9v9X3Ky8sZO3Zs6wsUINUhN6gOuUF1aJqZvdPQdQ0BiYgUKAUAEZECpQAgIlKg8moOoCGRSIQtW7ZQU1PT7Pd0796dysrKdixV+0vVobi4mL59+xIOh4MukojkmbwPAFu2bKFbt270798fM2vWez755BO6devWziVrX5988gldu3blo48+YsuWLQwYMCDoIolInsn7AFBTU9Oixr8jMTN69uzJ9u3bgy6KiGQpq6hizuINbN2xm949OnPtmYM4f0SfoItVR94HAKAgG/+UQq67SK4qq6jihoVr2B2JMdI2MuaTSv70UAnfeXAgfXIoGHSIACAikkvmLN6QbvzvL/oZISKAx1OxkczbOZHvPLibW/6yjh+dOzTQQKBVQG3A932GDx/OsGHDGDlyJC+99FKTr9+xYwd33HFH+vvy8nImTpzY3sUUkQNk647dAFzqP0knIoQMfOKc6S/noaJbmB/6Md+r/R1/eugh+l//V0689Rle2ho54OUsuABQVlHFGb95lQHJP/SyiqpW37Nz586sWrWK1atXM3v2bG644YYmX58dAESkY+ndozMjbSPn+YnOoHNglvjl4zjOX890/2keKrqF34V+Ra+dq5n3ei0jfvxkm7RJzVVQASA1Lvf+rj04oGrHbm5YuKZN/8B37drFIYccAkB1dTXjx49n5MiRHH300Tz22GMAXH/99WzatInhw4dz7bXXpl974YUXMnjwYKZPn45zjqeffpqvfvWr6XsvWbKEyZMnt1lZRaR9XHvmIM4IrcJPNvqQCAKZgSAVDFJPBb8L/YoBu9e1eZvUlA41B3DLX9bxxtZdjf684t0d1Mbida7tjsS47uHXeWDZuw2+Z0jvg/nRuUOb/Nzdu3czfPhwampqeP/993nmmWcAKC4u5tFHH+Xggw/mww8/ZMyYMZx33nnceuutrF27llWrVgGJIaCKigrWrVtH7969OfHEE3nxxRc59dRT+cY3vsH27ds5/PDDueuuu7jiiita8CciIkE4f0Qflpd/Bjsh5owoPhWxL1Pqb8DPOIY3FRx8lwgEp/kreCo2inkPTWTO4mHtPllcUE8A2Y3/vq43V2oIaP369fz973/nsssuwzmHc45Zs2ZxzDHHcNppp1FVVcW2bdsavMfo0aPp27cvnucxfPhwNm/ejJlx6aWXcu+997Jjxw5efvllJkyY0KqyisgB8N4yRu58Cgd4vs8bw2/iqvBPuKj2R9wXG8+rscHEsHpPBZlPBHOrb6D6kW8y45bftNsTQYd6AthXT/3EW5+hKjk5k6lPj848+B/Ht0kZjj/+eD788EO2b9/OE088wfbt21mxYgXhcJj+/fs3umO5U6dO6a993ycajQJwxRVXcO6551JcXMyUKVMIhTrUX5lIh1NWUcX2x//AlS6OGcTjjhGHxVn1ozMoqxjKnMXDqNqxm5Gxjcz0H+c0f0X6qSDzieA4fz3HsZ6L48/w1CNlTH5oItu6t+1TQUG1JteeOSi9Njelc9jn2jMHtdlnrF+/nlgsRs+ePdm5cyef+9znCIfDPPvss7zzTiIja7du3fjkk0+adb/evXvTu3dvfvrTn7JkyZI2K6eItL3UPOP17mMslBj+qcVnWe1gTiExNJRqvMsqhnP9oqEcWbOu0UAAmcNDK7lp1xXcsLAWkvdqrYIKAKk/sF/8rZJ/7trTZrvzUnMAAM457rnnHnzfZ/r06Zx77rmUlpYyfPhwBg8eDEDPnj058cQTOeqoo5gwYQLnnHNOk/efPn0627dvZ8iQIa0qp4i0rzmLN1ASrWR60dM4B3E8bolcytLXOvPi+LqvTQWDsoqh/GTxMObtXF0vEEDmU0Gcn4TvZmNtP+YsLlIA2B/nj+jD+C8f3Ka5gGKxWIPXDzvsMF5++eUGf3b//ffX+T7zIIjbb7+9zs9eeOEF/v3f/711hRSRdrd1x24m+ZX4JIZ/zDkOter0voCGpALBz++P8d3KEkpqK5nsL+XLVNWZNE7cL84Yr5Lf7hjYJuUtuACQb0aNGkWXLl2YO3du0EURkX3o3aMzFTu/hBnEHEQI8Uq8hN49Ou/zvSf0DjOkZAj/tcjjpt2JBn7vPMFKcI4I4WbfrzkUAHLcihUrgi6CiDTTtWcO4v5H1gJQFjuJ+2KnURkqYXYz5xn3DgslEsmt3DGQq6PXMCK2kTFeJa/ES1p0v31RABARaSPnj+hD9LWdsBVui11IvPsXmb0f84x1J4urmLO4M7/dMZDePTrv1/0aowAgItKGSj59jT0uxPNXDca+2Prl5ZnBoK0V1EYwEZF29d4ySnY+T5FFsT9/Fd5bFnSJmqQAICLSRtw/nsdwGECsFjYvDbpITVIAEBFpIx8fdCQGiRDgF0H/k4MuUpMUAALU2P4BEclPW3YlduluH3gxzFgE/UYHXKKmFWQA8LaugKVz22R87rXXXuOYY46hpqaGTz/9lKFDh/L666/z9a9/naFDhzJx4kTOPvtsHn74YQD69+/Pj3/8Y0466SQWLFjQ6s8Xkdzhv/kkAF2GX5DzjT90tFVAf7se/rmm6dfs2cVB29aCi4N50Oso6HRw468/4miYcGujPz722GM577zzuOmmm9i9ezeXXHIJGzduZPPmzaxZs4YPPviAkpISvva1r6XfU1xczAsvvNDS2olILntvGSXvP4oDuiy8FLrpCSD31OxMNP6Q+L1mZ6tv+cMf/pAlS5awfPlyrrvuOl544QWmTJmC53kcccQRjBs3rs7rp06d2urPFJHcUVZRxe/uuQfPxTEgHs39CWDoaE8ATfTU095bBvecC7FIYpLmgj+0Okp//PHHVFdXE4lEqKmpwWUkcmpIly5dWvV5IpI7UhlAS6N9uKpT/QyguazwngD6jeazKQ/CqTe22STNzJkz+clPfsL06dP5wQ9+wEknncQjjzxCPB5n27ZtlJeXt77cIpKT5izewO5IDC+ZtXNB7BSm185i1mttk6+nPXWsJ4BmivceBYPGtsm9/vSnPxEKhZg2bRqxWIwTTjiByZMn07dvX4466igGDhzIcccdR/fu3dvk80Qkt6QyfU7wXiXuYFHseFa6gVgTGUBzRUEGgLZ02WWXcdlllwGJk7xeffVVIHHEY9euXfnoo48YPXo0Rx99NACbN28Oqqgi0g569+hMr52rmeI/hwF3Fs1leu0stnUfFnTR9kkBoJ1MnDiRHTt2UFtby80338wRRxwRdJFEpB1ce+YgNj96Px4OMwi7KCeF13PkmRcFXbR9UgBoJxr3FykM54/ow6uvhrB/JiaAoxZi1FfO45R2SuDWljpEAHDOYZmHaBaQfa04EpF29t4ySv/5EM6B5/sUnz2HU0qbPuY1V+T9KqDi4mI++uijgmwInXN89NFHFBcXB10UkcK1eSkesfQRkOz+KOgSNVtgTwBm1g/4E3AEEAfmOed+3dL79O3bly1btrB9+/Zmv6empibvG81UHYqLi+nbt2/QxREpSGUVVTxfHuFXQBwjbmFCOZ4ALlOQQ0BR4HvOuZVm1g1YYWZLnHNvtOQm4XCYAQMGtOiDy8vLGTFiRIvek2s6Qh1E8llqA9jM+AYIw9+jpfzZzmPqh304v1/QpWuewIaAnHPvO+dWJr/+BKgEcn/WRESExAawkmgl3wyV4RyM81ezJxpjzuINQRet2SwXxs7NrD/wPHCUc25X1s9mAjMBevXqNWr+/Pmt/rzq6mq6du3a6vsESXXIDapDbgiiDpf//VO+7pdxbeghzCDqPH4VncIdsUncfVbL0720Zx3GjRu3wjlXmn098FVAZtYVeAT4TnbjD+CcmwfMAygtLXVjx45t9WeWl5fTFvcJkuqQG1SH3BBEHfq88gxFn0Qwg7iDCCFeiZfQp0fn/SpLEHUIdBWQmYVJNP73OecWBlkWEZGW+Pmxu/lGaBEAMTxuiVxKZaiEa88cFHDJmi/IVUAG3AlUOud+FVQ5RET2xylF63GWONXPgP4H7WH2OUdzfh5sAEsJcgjoROBSYI2ZrUpem+WceyK4IomINE9Nr5F0cuDMCIU7cdWMGdAvfxp/CDAAOOdeIBE4RUTyznuVy/gXg496f4WeE27K+dO/GpL3O4FFRA60557+KwNW3opz0GXryzy3sfkbUXOJAoCISAuUVVSx4rlF+MQxg5CLsuK5RZRVVAVdtBZTABARaYE5izfgYnWXf74QGZxXG8BSAt8HICKST3rtXM23isqAvcs/8+UEsGx6AhARaYEzurxFiL3LPw+1aiBxMli+UQAQEWmBkX2K6+3+7Rz282oDWIqGgEREmuu9ZYx65y6cg5j53FJ7Kdu6D2P2mYPyagNYigKAiEgzlFVU8cHjd/LvLnH4iw/cOqEPnHxq0EXbbxoCEhHZh1Tu/+57tmKWOPt3j/N5rnZw0EVrFQUAEZF9SOX+n+IvxTmIJ1f/zHot/yZ+MykAiIjsw9Ydu7nS/yueucTZvzgOtWq25uHSz0wKACIi+3D6we9whr8c50hMAOPxSrwkL5d+ZlIAEBFpQllFFefU/h2fRO8/jrEgdkre5f5viAKAiEgjyiqqeHDhI5zjlmKW6P1H8VkSHsfsyfmV+78hCgAiIo2Ys3gD57hn8YkDe3v/b3UamveNPygAiIg0qtfO1VzkP1en978wdnLeT/6mKACIiDTijC5vEU7m/Un1/le6gXk/+ZuiACAi0ohR/bql8/7UEmZh7OS8zfvTEKWCEBHJUlZRxRN/e4w79vweB8TMy/u8Pw1RABARyZBK+/BD9yS+nzj1yxxcOqwbt07N37w/DdEQkIhIhlTahwszJn9jePx6U6+gi9bmFABERDJs3bGbyf5SQllLP5fs+mLAJWt7CgAiIhlOP/gdpvjl9ZZ+dpSVP5kUAEREMnz7S9vqLf3sCGkfGqJJYBGRpLKKKjZt/IQh7F36+fxBpzF7Qv6nfWiIAoCICHvz/tzj3Z089MVjdnwGZ0+Y1CEbf9AQkIgIkFj9M8k9lR7+Aega38WcxRsCLFX70hOAiBS0sooq5izewMm7HueicN2ln6/ESzpM3p+G6AlARApWatPXybse5+fhO/EscT0OHS7vT0MUAESkYM1ZvIFJsSfrNP6JM3/9Dpf3pyEaAhKRglRWUUWvnav5SdFdJNv+ZONv3By5vMPl/WlIoE8AZnaWmW0ws7fM7PogyyIihSM19JPa8Zsa949j3Bj5GksPnsiL15/aoRt/CPAJwMx84H+A04EtwGtmtsg590ZQZRKRji014Vu1YzcXe0/zr+Gn60z63hS5gsf8M5jdgYd9MgU5BDQaeMs59zaAmc0HJgEKACLS5lK9/pJoJTeHHucMf3l66CcOzI+NY358PLdN6ZibvhoSZADoA7yX8f0W4LiAyiIiHVhZRRXfe2g1U+wpflp0F35y2AfqTvr26dG5YBp/AHPOBfPBZlOAM51zVya/vxQY7Zz7ZtbrZgIzAXr16jVq/vz5rf7s6upqunbt2ur7BEl1yA2qQ25oqg4vbY1w99paJvM0Pwv/EQ+X1fgnxv0XMp7LjyrihN7hA1jyvdrz72HcuHErnHOl2deDfALYAvTL+L4vsDX7Rc65ecA8gNLSUjd27NhWf3B5eTltcZ8gqQ65QXXIDY3VoayiijsXr2aK1W38U/3e1Lj/0oMn8suAV/wE8fcQZAB4DfgXMxsAVAEXA9MCLI+IdCCpMf8p9hQ/D9+JQcZqH1gSK+Uem8TUKRdwawEN+2QKLAA456Jm9p/AYsAH/uicWxdUeUSkY8hc6TPdW8JPwndlbfJKDPkscKcx96JhBTXmny3QjWDOuSeAJ4Isg4h0HKle/+5IjOv8+7k69Hh6vD9lSWwUj/lnMHdy4az2aYx2AotIhzFn8QZ2R2LMDs3jYr88fT015h/B5w/xc5ldQEs9m6IAICIdxtYdu7nGfzDd+Dc05n/JlAvU+CcpAIhIh1BWUcVor5L/DD0G0GB6h2s7eG6fllIAEJG899LWCH9+eg33+g/Um/C9JX4lY6Z8t2BX+jRF6aBFJO89sjHCjNhCRnpv4dze3D43R/+NkZO/q15/IxQARCTvfXHPBq4LPwgkhn5SuX0eiHX8jJ6toSEgEclbqTX/PwstbPBAl458mldb0BOAiOSl1Jr/L+5cxle81+sO/UQupzJU0qFP82oLegIQkbyUWvN/U/jedO8/5oz5sXEsPXhihz/Nqy0oAIhI3imrqKJqx25uDP2ZEu+99EavaHLo58XrTw22gHlCAUBE8kJmjh8DZnqLuNL/G5CY+I05WBA7hW3dhwVb0DyiACAiOa2soor/WrSOHbsjjLSNXB1ayjG8xdH+O0Dmhi+fv9pYjfu3QJMBwMx27eP9BrzvnBvYdkUSkUKX3dsfYRuZGfoLp/sr8Njb6Gfu9r05cjlTleahRfb1BLDJOTeiqReYWUUblkdEClxmRs+RtpGr/EWM91emG/6UhlI9aLdvy+wrAFzQjHs05zUiIvuUOrt3GBv4eqiMU/1V6YNcUjJPsU2d6LWQ8fxSQz8t1mQAcM69nfm9mR2c+R7n3MfZrxERaYnM4Z6RtpE/+o9wsr+m0YY/hrE8Noi36MPC2Mls6z6My78Q09DPfmjWJLCZ/QfwY2A3kIq/DjiyncolIh1c9uTunNB8xvjrSbX5mQe3Q6Lhfyo2inmxifyj81B+dO5QFiYb/fLy8gNe/o6guauAvg8Mdc592J6FEZGOL9XwH1mzju/7z1MSfpcR3lt4GRO70HDDv637MKV0bkPNDQCbgM/asyAi0nFlD/P80l/E+KKVdVb0kPV1quH/Q/xcLpkyJd3bl7bT3ABwA/CSmb0K7ElddM59q11KJSIdQmZvf5JXSU9/B5eHFtdb0QMN9/grQyU6vrEdNTcA/B54BlhDItOqiEiDsnv7t/qPc1rRSjzijU7sRvB5JjaCD+nOwtjJrHQD6dOjs/L5tLPmBoCoc+6adi2JiOS1hnbsTvWfJZTsM2YO7zTU21/pBnLIQeE6k7vSvpobAJ41s5nAX6g7BPRxu5RKRPJC4zt2V+LhGuztRx2Ax1OxkXUa/tvOHare/gHW3AAwLfn7DRnXtAxUpIBl79i90H+Oi/xy/IyGP3vT1rzo2VTThVfiJelhnts0zBOYZgUA59yA9i6IiOSPsooq7l2wgCvsjeTE7pMN9vhjyZ5+5ti+hnlyx76SwY10zq1s7WtEJP9lDvdc7D3NA+G78LMmdhtK0/BgfDwO1NvPQft6ArjLzMYC1sRr7gSaTBgnIvkrcynn1f5SjgxtZYxfmW74Myd2s3v8/+g8lP/W2H7O2lcA6A6soOkAsL3tiiMiuSJ7x+7FReWJHn9WKmZITOw6fG6OXM78+Hj69OisHbt5YF/J4PofoHKISA7IXsN/nf8cU4qeI5Rs+DPFkoew/G90AtV0YZkbwiVTpiglcx7RiWAiUmcNP8AV3hPcFL6vzsRu5vh+BJ+HYmPTE7udwz6zJ2vHbr5RABApUHV7+w8xzavE+Y7zvJcY7G2p0/DHMaL4PBsbnh7fr3AD05O7Gu7JTwoAIgWmoR27F/nlhIjVm+xLTezOj41L9/YBDjkorMndDqC55wE8Dcx1zj2RcW2ec27m/nyomc0BzgVqSWQavcI5t2N/7iUi+9bQjt3JoeeZ6j+XaPgzevupyd3UGH9qYhdIr+FXw98xNPcJYADwAzM71jl3S/JaaSs+dwlwg3Muama/ILHD+AetuJ+INCJzxy7AdO9JbgnfU29835FY7hd1RgyfBbFT6iRm0zBPx9PcALADGA/8PzP7C3BJaz7UOfdkxrevABe25n4i0rDUGbsx5zjZVnNNaAHDvbezxvchSogFsVNYG+/PoVbNK/GSeqduScdjLnNqv7EXmVU450Ykv74c+B5wiHOub6sLkAgoDzrn7m3k5zOBmQC9evUaNX/+/NZ+JNXV1XTt2rXV9wmS6pAbcrEOL22N8MjGCB/VJP5vT/Re4pv+owz0quqt4W9ofL9nsXHBwDAn9A4HWY0WycW/h5ZqzzqMGzduhXOu3qhNc58Afpf6wjl3t5mtAb7R1BvM7CngiAZ+dKNz7rHka24EosB9jd3HOTcPmAdQWlrqxo4d28wiN668vJy2uE+QVIfckEt1yNy4NdV7g56hXZzpLaOv9zHZ/bzs8f18z8aZS38P+yuIOjQ3Gdzvs75fAXxtH+85ramfm9kMYCIw3jXnMURE6klN7vbauZoxXiXnumJuLrqXcHJi17nsid3Ecs7U+P627sOY+YUYs6adHnRVJACBLAM1s7NITPqe4pzTWcMiLZSdn+eioucIEa2XmK2hhr8s/pU6Z+yWl5cHVg8JVlD7AG4HOgFLLPGv9RXn3FUBlUUkb9TNz7OUi4r2rt83g7iDmPNIZOZx9Vb06IxdyRRIAHDOfTmIzxXJR9n5ea71n0/2+Ouu34872EMRt0Qu5VCr5mPXNb2iR2fsSkO0E1gkR2Xv2P126Bkm+y+kM3LC3vX7keQyzszVPIAOX5EmKQCI5JDs3v40r5KQH+VboUfrNfyp/DyZDX8qlYM2bklzKACI5ID6+XnqpmnInNBtqOEHpWmQllMAEAlIQ739aq8TN4XvSy/jhNSafQ9LbtxSwy9tRQFA5ADLzr0/0jZyf9HPKCJSbxlnHKjNmNh9JV6iNMzSZhQARA6A7GyciZ2PjmneU3w39AidiGQt42y8t680zNJWFABE2lF2bx+g1CqZ4T/JcO/NOmkaYs6oJazevhwwCgAi7SQ7DfNYq+A7oUcYlszGGU+u3fcskYL5xfhR/Dp6gXr7csAoAIi0g1Qa5mFs4NzQyxzF25T6b6ZX8kBiNY/DI+4cEULpxl+9fTlQFABE2khmYrbjvTf4kf8x0/2n8c3VOXAFEsM9kYzhnnVFR3PZRZO1YUsOKAUAkVbKzM/zLf9ZLihamt60lR7fz+jtZ07ubus+jGvPHMR1avglAAoAIvuh7qqeVZzlvcpvin5TZ7duKue+JYd4Mnv7Xz1fvX0JngKASAvUXdXjuNRfwmX+k3zJtqbX8Kd262av6FnNIOZOGabevuQMBQCRfchewz/a3mBGeDEj7C0+7/1f8mjFRNpl38UbXL/fOewzd7LSMEtuUQAQaUR2fp5p/koG2buM91fVW8aJMxbETmGrOyydflmJ2STXKQCIZGgoP09nv4arQ3/BJ17ntdnLOJWfR/KNAoAI9Xfsltp67i/6OWGidbJxRl3irK3sid3MQ1fU25d8oQAgBS274T/JXudrob8x2ltPkUWBxlfzpBp9Hboi+UoBQApO9qTuCNvIhNAyhvIPjvcrMUs0+hHnY9Rv9DPz88z8QoxZ004Pukoi+0UBQApG5oatSV4lr1gJx3ibuDl0b73dug6PB2Nj60zqQv38POXl5YHVR6S1FACkQ2uot39f0c/Tufe9ZMMPdXfrpiZ1lY1TOjIFAOmQGjp05TRvBSd7r1NMbXpi98VoCSP9TYRdtN5Qzz86D1U2TunQFACkw2iotz/Gr+QgPuOq0F/TyzjjGDioJczc2FSIwRivst6krhp+6egUAKRDyM69X2qV3Fc0u4FlnB7zY+PqbdhaGUss4bxNwzxSQBQAJO9l5t4f66/iEKtmkv9io8s4tWFLJEEBQPJSZu79MV4lX/eifDv0aDob55uxz9PJi+AT14YtkUYoAEheqZufZwMPFP2MonrDPMaj8a/wSrSkwbF9bdgSSVAAkJyXnZ/nMm8NXUI1XOg9T6c6wzwe5iBCKL1hS2P7Io1TAJCclb1xK+7B98ILCCWHed6LH0Y3txuvgWGe1Qziv6cOU6Mv0gQFAMk52WmYHyj6aWI1T/LnqWGeB2Kn8kp8SJ1hHlDufZHmUgCQnJA9zHOJtw7fj3Gx/2x6mCfuIFZnmGcIK5PDPMq9L9JyCgASqOwdu2NtFf9bNJcQMczgg/jB1Dofr4GkbFrKKdI6gQYAM/s+MAc43Dn3YZBlkQMn+0D1EbaRiaGX6ccHjPNXE7LEjt2YM+6OnVVvmCeRrE29fZHWCiwAmFk/4HTg3aDKIAdW/Wycgyn1NvKD0Px0Ns5XY4MZ7m8iRKzeMA+oty/SloJ8Avhv4DrgsQDLIO2s4WycP6MTiSEfL7l+HxLj+8+7Yfyy9uJ6E7tq+EXanrnU/74D+aFm5wHjnXPfNrPNQGljQ0BmNhOYCdCrV69R8+fPb/XnV1dX07Vr11bfJ0j5UIeXtka4e20ttfFENs4z/Nf4ir1OifdeeuPW0thRHOtvJEwiG+f02lnpRh+gZ7FxwcAwJ/QOB1iTxuXD38O+qA65oT3rMG7cuBXOudLs6+32BGBmTwFHNPCjG4FZwBnNuY9zbh4wD6C0tNSNHTu21WUrLy+nLe4TpFyvQ1lFFXcuXk3Mwdney/wmfDseic5GHMM5iBDmttiFeZ2NM9f/HppDdcgNQdSh3QKAc+60hq6b2dHAAGC1mQH0BVaa2Wjn3D/bqzzS/rKHe3yifN1/nGtCD+PhlI1TJMcc8DkA59wa4HOp7/c1BCT5ITsd81TvGa4JLeBz3k5ejJUwynuLkIspG6dIDtE+AGm1VDrmmHOEiDI79Acu9J8HoNaFmBudCtQd5unTozPn6EB1kUAFHgCcc/2DLoPsv1TPP+YcfW07vw7fzijvTeIuscLHc3HGeJXcEZvE6vgg5l40LJ2NUweqiwQr8AAg+SvzIJYZocWc5q8kjsevIhdwdegv6XN2X4mX0DnsM1v5eURyigKA7JdUz/9Y1nJv0WxCFifujK9HvsXf48fxQvzo9JDPtu7DmK3JXZGcowAgLZK50mewvcvt4d+kUzfEMY60xEKu1Qzisgsv4jo1+iI5SwFAmi3V6y+JVvJf4UWc4q2mms7scaH00Ysa7hHJHwoA0iyp8f4TqeCuojn45og545raq9lFl/Rwz2oGKRe/SJ5QAJBGZW/sGm2V3F70m/SOXocxxHuXO2KTWBkbqINYRPKMAoA0KHNjV6mt57uhRxjjrWObO4ROFsV3sfSQDyg1s0g+UgCQejI3dp1hr/HbotvwzRF1Ht+PXEUNndJDPpWhEm5Tr18kLykASB2ZG7vO9JZxW/h/0kM+AMO9t9NDPr6ZhnxE8pgCgAB7x/t77VzNld5ahoY2c1ZoOW/Ge9PPtqfz+KSGfLTSRyT/KQBIneWd9ycPazGDsugJXBu9iqPt7Xp5fDTeL5L/FAAKXGq8/yD3KTeE7083/jFnbHT9iBBipRtYL4+PiOQ/BYAClDncc7z3Bt/x93Cx/yw92UUMDxwa7hEpAAoABSR1KPuO3RFG2gYeKPoZRUQxgw3xPlwRuY4iohruESkQCgAdXGZvf4xXydHx/gz0t/AN/zE6WRRIDPc8FjuRte5IgPSmLi3vFOnYFAA6qLq9/Y08UPRTwkQxwAzeiPejq6vBS+fwGZJ+r2+mIR+RAqAA0MGUVVSx8LGFHBVZw+D4l+nnb+caf0G6tx93cHf0TG6JzmCkbawz3AMa7xcpJAoAHUBmzp6RtpH7i35GUSiS7u2/Gz+MWufj4YgQ4i+x4wFY6QayMpY4lB2UzkGk0CgA5LHM3v7n44MY4lUzK3QfxRYBEr39+dFx3BC9kpH2Zr3ePuhQdpFCpgCQZ+r29h9K9/YhcQbvh64bEedjyd7+gtgpgKm3LyL1KADkicxJXXCMtvXcEr47vXEr7uDR6Il8P3oVw2yTevsisk8KADksOx//CbaGGeEnGWTv8kVvO9WuU3LjliNCmD/HTieGn+7tp6jhF5GGKADkoLq9fRhk7/Id/2HO8pcne/vG7dFJ3B49nyH2Tr3evgEODfOISNMUAHJEdm+/1Cq5LLSEL3lbGeK9S9QlkjIbicPXP3PF1NBJvX0R2W8KAAHL7u335kOuCS3gAn9pemz/rugZPBsbzu+LbiPsonXy9Ki3LyL7SwEgIKmG/8iadUzz3qDGL+J4r5JTvZUYLqO377HdHcLzbjjTa2fVGe5Rb19EWkMB4ADKHuY52VZxZ9FcQsQwgx3uIH4bO4/XYwP4ddEd9Xr7Fcnhnj49OjPzCzFmTTs92AqJSF5TADgAMnv7k7w3+NA7mNH+Rs73XiBkcSCRkO3O6Nn8JjYZgOm1PZrs7ZeXlwdVHRHpIBQA2kl2b/84W8efim4lnOztf+aKeCo2krH+akIkjlt8MX5U+v2Zvf3bNLYvIu1AAaCNZU/qDrD3me4/xSX+UxRZDEj09udFJ3Jb7EJGxuonZNPYvogcCAoAbSC7tz/K1jMt9DQD7H1G+G8TcT6vxgdzrLcBP5l++fn4MQB1lnGq4ReRA0kBoBWye/uH8398N/QwF/vP4iWXcN4XPZXboheynR710i9rCaeIBCmwAGBm3wT+E4gCf3XOXRdUWVoiu7c/wjYyJfQc/ewDjvPWE7YYziVeG8ejyh3OdnoA6u2LSG4JJACY2ThgEnCMc26PmX0uiHK0RHZvvyuf8a3QQq70/4ZnDudgUex4FsdLmRv+vTZsiUjOC+oJ4GrgVufcHgDn3AcBlaNJ2QetTPMqed87lFJvI+f7L9DF9qR7+zE8Nrgv8ET8eP5Z21MbtkQk55lLtWAH8kPNVgGPAWcBNcD3nXOvNfLamcBMgF69eo2aP39+qz+/urqarl27Nvrzl7ZGuO+NWj5NnKLIaKvk3qLZiTN1DfY4n8diJ/FafCA/Dt9DmERvf3rtrDrpl3sWGxcMDHNC73Cry9zSOuQD1SE3qA65oT3rMG7cuBXOudLs6+32BGBmTwFHNPCjG5OfewgwBjgWeMjMjnQNRCPn3DxgHkBpaakbO3bsfpdpb4/e6NG5FjP4v88i+GbEnKNH5zC10RiDo+uZ4VXyln2eYf7bXO4vpih5pm7MGb+Lnsd/x6YAsKm2TyC9/fLyclrzZ5ELVIfcoDrkhiDq0G4BwDl3WmM/M7OrgYXJBn+ZmcWBw4Dt7VKY95ax7qW/8uDanlRFvgTAjt0RRtpG/tVPNN7r6E+n3R8zzlvHL4r+lxDRZC4eWB4fxHBvU3oJ53PxYelba8OWiOSroOYAyoBTgXIzGwgUAR+2yye9twzumsCQeJT7PGNLUU9q6EQxe+hrH2EkHjrM6r817mBedCK/iE2rt4QTtJJHRPJbUAHgj8AfzWwtUAvMaGj4p01sXgrxWHIVjuMzitnkevMl2wq4dMrl52NHszg+mh58wrdDC9O9/SXxxLDZSjeQithAreQRkQ4jkADgnKsFLjkgH9b/ZAgVE43uIUKIWZErWekGMtI2cl/Rz9PLNX8dvSDds38lPkS9fRHp8Dr+TuB+o2HGIja89Fd+urYnK11iDmClG1gvv36KNmyJSCHo+AEAoN9ohk4dzdSBVbybXNffo3OYf9hQVn42EN8SA0Q9OocTefk/i9Bbwzwi0sEVRgBIOn9EH84f0adDLBkTEWktL+gCiIhIMBQAREQKlAKAiEiBUgAQESlQCgAiIgUqkGyg+8vMtgPvtMGtDqO9Uk8cOKpDblAdcoPq0LQvOucOz76YVwGgrZjZ8oZSo+YT1SE3qA65QXXYPxoCEhEpUAoAIiIFqlADwLygC9AGVIfcoDrkBtVhPxTkHICIiBTuE4CISMFTABARKVAFFwDM7Cwz22Bmb5nZ9UGXpzFm9kcz+yB5alrq2qFmtsTM3kz+fkjGz25I1mmDmZ0ZTKn3MrN+ZvasmVWa2Toz+3byej7VodjMlpnZ6mQdbklez5s6pJiZb2YVZvZ48vu8qoOZbTazNWa2ysyWJ6/lWx16mNnDZrY++f/i+MDr4JwrmF+AD2wCjiRxDvFqYEjQ5WqkrF8BRgJrM679Erg++fX1wC+SXw9J1qUTMCBZRz/g8n8eGJn8uhuwMVnOfKqDAV2TX4eBV4Ex+VSHjLpcA9wPPJ5v/5aS5doMHJZ1Ld/qcA9wZfLrIqBH0HUotCeA0cBbzrm3XeJYyvnApIDL1CDn3PPAx1mXJ5H4R0Ty9/Mzrs93zu1xzv0DeItEXQPjnHvfObcy+fUnQCXQh/yqg3POVSe/DSd/OfKoDgBm1hc4B/hDxuW8qkMj8qYOZnYwiU7dnZA4Ftc5t4OA61BoAaAP8F7G91uS1/JFL+fc+5BoYIHPJa/ndL3MrD8wgkQPOq/qkBw6WQV8ACxxzuVdHYDbgOuAeMa1fKuDA540sxVmNjN5LZ/qcCSwHbgrORT3BzPrQsB1KLQAYA1c6wjrYHO2XmbWFXgE+I5zbldTL23gWuB1cM7FnHPDgb7AaDM7qomX51wdzGwi8IFzbkVz39LAtcD/HoATnXMjgQnAN8zsK028NhfrECIxpPtb59wI4FMSQz6NOSB1KLQAsAXol/F9X2BrQGXZH9vM7PMAyd8/SF7PyXqZWZhE43+fc25h8nJe1SEl+bheDpxFftXhROA8M9tMYsjzVDO7l/yqA865rcnfPwAeJTEckk912AJsST5BAjxMIiAEWodCCwCvAf9iZgPMrAi4GFgUcJlaYhEwI/n1DOCxjOsXm1knMxsA/AuwLIDypZmZkRjvrHTO/SrjR/lUh8PNrEfy687AacB68qgOzrkbnHN9nXP9Sfx7f8Y5dwl5VAcz62Jm3VJfA2cAa8mjOjjn/gm8Z2aDkpfGA28QdB2Cnhk/0L+As0msSNkE3Bh0eZoo5wPA+0CERG/g34CewNPAm8nfD814/Y3JOm0AJuRA+U8i8cj6OrAq+evsPKvDMUBFsg5rgR8mr+dNHbLqM5a9q4Dypg4kxs9XJ3+tS/2/zac6JMs0HFie/PdUBhwSdB2UCkJEpEAV2hCQiIgkKQCIiBQoBQARkQKlACAiUqAUAERECpQCgIhIgVIAENkPZtbfzHYn8wS15H1Tkyl+H2+nook0mwKAyP7b5BJ5gprNOfcgcGX7FEekZRQARLKY2bFm9nryQJguycNgmkoCl3oiWJ/M8rjWzO4zs9PM7MXkYR+5mlJZClgo6AKI5Brn3Gtmtgj4KdAZuNc5t3YfbwP4MjAFmEki79Q0EikxzgNmsTfXu0hOUAAQadiPSTTiNcC3mvmefzjn1gCY2TrgaeecM7M1QP92KaVIK2gISKRhhwJdSRxnWdzM9+zJ+Dqe8X0cdbYkBykAiDRsHnAzcB/wi4DLItIu1CsRyWJmlwFR59z9ZuYDL5nZqc65Z4Ium0hbUjpokf2QPOf4cedck6uDGnnvWOD7zrmJbVwskRbREJDI/okB3fdnIxhwB/B/7VEokZbQE4CISIHSE4CISIFSABARKVAKACIiBUoBQESkQP1/ZGD3u8Fp9LYAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "xgr,zgr = xgrid(x, bathy,dxmin=2)\n", "\n", "\n", "plt.figure()\n", "plt.plot(x,bathy,'-o')\n", "plt.plot(xgr,zgr,'.-')\n", "plt.legend(['Bathy','xgr'])\n", "plt.xlabel('x [m]')\n", "plt.ylabel('z [m]')\n", "plt.grid('on')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Interpolate\n", "Interpolate data to xb grid." ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYAAAAEWCAYAAABv+EDhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAq7UlEQVR4nO3de5hU1ZX38e+qU9WIkIiXDAlIgiYDNmoQ7CAmURvvRiOOijpegpnx5dHMTC5OdLxlEpMYkvCamDeZJMPEqIkoYiQtMZlBMbYaDUGgRcUGxHjBJvEyBg2KdNWp9f5Rp4rq6uqmoas5VV2/z/PwdPep6lN7C+61zzp7r2PujoiI1J9E3A0QEZF4KACIiNQpBQARkTqlACAiUqcUAERE6pQCgIhInVIAkLpmZq1mdlEf39tsZi9V8LMvNLPfVep8IjtKAUBkFzCzsWbmZpaMuy0ieQoAIiJ1SgFABjUz+6CZvW5mk6OfR5nZa2bWXPS2D5rZMjN7w8zuNrO9tnPOq6JzPG9m5xUdP9nM2szsTTPbYGZfKfq1h6Kvm8xss5kdXvR7/9fM/mJmz5nZSdGxGWa2ouRz/9XMWnbiP4NIWQoAMqi5+7PAvwHzzGx34CbgZndvLXrbp4B/AEYBGeD/9XLK9wL7AKOBmcBcMxsfvfZWdK4RwMnAJWZ2WvTakdHXEe4+3N1/H/18GLA2Oue3gRvNzIBFwH5m1lj02ecDP+9z50W2w1QLSOqBmS0C9gMc+Ii7b42OtwJL3f2K6OcJwOPAUHcPS87RDCwB9nD3t6JjC4An3f1rZT7zBsDd/QtmNhZ4Dki5eyZ6/ULgGnf/UPTz7uSCyPvc/c9m9iPgdXe/2swOBH4HvDffdpH+0hWA1Iv/Ag4Cvl9mAN1Q9P0LQIrcjLycv+QH/6L3jwIws8PM7AEze9XM3gAu7uU8eX/Of+Pub0ffDo++3gKcG10RXAAs0OAvlaQAIIOemQ0HbgBuBL5SJsc/puj79wNp4LUeTrenmQ0ref/G6PvbyKVuxrj7HsCPAYte2+FLbXdfCnQCRwDnovSPVJgCgNSD7wEr3P0i4NfkBuZi55vZhCgF81XgF6XpnxLXmlmDmR0BnALcGR1/F7mUzTtmNoXcoJ33KpAF9t/Btv8M+AGQcXftGZCKUgCQQc3MpgMnkkvHAFwKTC5evUNuZn0zuXTMbsBneznln4G/kJv1zwMudvc10WufAb5qZn8F/h1YkP+lKL1zHfCImW0ys6l97MLPyaWuNPuXitNNYJEqZmZDgVeAye7+TNztkcFFVwAi1e0S4DEN/jIQtC1dpEqZ2fPkbiKfFm9LZLBSCkhEpE4pBSQiUqdqKgW0zz77+NixY/t9nrfeeothw4Zt/41VTH2oDupDdVAferdixYrX3P09pcdrKgCMHTuW5cuX9/s8ra2tNDc3979BMVIfqoP6UB3Uh96Z2QvljisFJCJSpxQARETqlAKAiEidUgAQEalTCgAiInWqplYBiYjUipa2DuYsXsvGTVsYNWIol50wntMmjY67WV0oAIiIVFhLWwdXLnySLemQybaOqX9t52cLGvn8HeMYXUXBQAFARKTC5ixeWxj8b2u4jiRpIMGScDJz3ziFz9+xhWt/tZovf/LAWAOBAoCISIVt3LQFgAuCexlCGjNwz3JCsJxjgxUsD8ezvnM0P1twROGq4OT3hzTv4nbqJrCISIWNGjGUybaOU4NHAXAHs9yfAOewYA3nBfezoOFafpz8DiPfWMXcJzqZ9NV7aWnr2GXt1BWAiEiFXXbCeP688BaC6InQ7rk/kAsCeYF74apgSXgoc7ecwpULswC7JDWkKwARkQo7bdJomvZ6G4DQjU6S/CE8gBArBIPSq4ITguUsaLiW7/ocfrZgAR/75m8H/GpAVwAiIpW2YRmT31iCA4kg4OmDr+HiJw5k/3dWc3rwMB+ig6ZgLUF0WZC/Kii+Ili+eTzr7xrNzEXT+LtTTx+QKwIFABGRCmpp6+DVe37CRZ7FDLJZZ9I+WR7/8vG0tB3InMUT6di0hcnhOmYF93BssKJsIDgsWMNhrOGc7G9ZclcLpy84hZf3mFjRJaQKACIiFZJf/3+Fv44l8+mfgGWdB3AUudRQfvBuaTuEKxblrgp6CgRQfFWwkmve/DRXLuwEKnOPQPcAREQqZM7itTRm2jkvuB93yJLg2vQFXPXY0G7vPW3SaB7/8vF86qyz+Nrwqzmr88ssDpt6uU+Q5Wupm2nMtDNn8dqKtFdXACIiFbJx0xamB+0E5NI/5s5etrmwL6Cc/FXBN24L+UJ7I42d7T3eJzDPMjXRzo82jatIexUAREQqZNSIobS98UHMIHRIk2RptpFRI7pfAZT66KgUExon8JVFCa7Zkhvgt90nWAnupEn1+Xx9oQAgIlIhl50wntvuegqAlvDjzAuPpT3ZyOwTxvfp9/NXA/lCcis3jeOSzKVMCtcxNdHO0mzjDp1vexQAREQq5LRJo8k89gZshBvCM8nu8QFm78Sqna43izuYs3goP9o0jlEjhu7U+XqiACAiUkGNbz3GVk/y0MUHYB84vN/nKw4GlaZVQCIilbJhGY1vPESDZbCf/x1sWBZ3i3qlACAiUiH+3EMYjgGEnfD8w3E3qVcKACIiFfL67vtjkAsBQQOMPSLuJvVKAUBEpEJeejO3S/fVcefAzEUwZkrMLeqdAoCISIUEz9wLwLBDzqj6wR8UAEREKmPDMhr/9EscGLbwgqq/AQwKACIi/dbS1sGPb7mFhGcxIJup/hvAoAAgItIv+Qqgj7w1OioBYWz1gAc7D4i7aduljWAiIv0wZ/FatqRDEtF0+s7wKBaEzbz82FAeOSbetm2PAoCISD/kK32elPgDWYdF4eGs9HFYLxVAq4VSQCIi/TBqxFAm2zpmBA9iwI0N1zPZ1lWsYudAUgAQEemHy04Yz5Gp1SRwzCBFho+n1nBZhSp2DiQFABGRfjht0mgOH5Us3ADOWJJDjzp1wAq4VZLuAYiI9MeGZTT9eQHukAgCdvvEHI5qOjnuVvWJAoCISH88/zAJwtyD3N1hy//G3aI+iy0FZGZjzOwBM2s3s9Vm9rm42iIisjNa2jq4tDWd2/yFkbFU1ReAKxbnPYAM8K/u3ghMBf7JzCbE2B4RkT7LbwAb885aAP4n08QF6atpea36c/95sQUAd/+Tu6+Mvv8r0A7Uzn85EalrcxavpTHTzr8kW3CHacEqtmZC5ixeG3fT+szcPe42YGZjgYeAg9z9zZLXZgGzAEaOHHno/Pnz+/15mzdvZvjw4f0+T5zUh+qgPlSHOPpw4f+8xWeCFi5LLsAMMp7gO5kZ/DCczs0nDtvh8w1kH6ZNm7bC3ZtKj8d+E9jMhgN3AZ8vHfwB3H0uMBegqanJm5ub+/2Zra2tVOI8cVIfqoP6UB3i6MPopb+l4a9pzCDrkCbJ0mwjo0cM3am2xNGHWPcBmFmK3OA/z90XxtkWEZEd8Y2PbOGfkosACElwbfoC2pONNbEBLC+2KwAzM+BGoN3dvxNXO0REdsZRDWtwCwEwYOzuW5l98sE1sQEsL84U0MeAC4Anzezx6NhV7v6b+JokItI374yczBAHNyOZGsLFM2fCmNoZ/CHGAODuvyMXOEVEas6G9mX8rcH/jjqSvU+6piYeAVlKtYBERHbQg/f/mv1WfhN3GLbx9zy47tW4m7RTFABERHZAS1sHKx5cREAWM0h6hhUPLqKlrSPupu0wBQARkR0wZ/FaPOy6/PN36QNqagNYXuz7AEREasnIN1bx2YYWYNvyz1p5AlgpXQGIiOyA44etJ8m25Z972WaAmngCWCkFABGRHTB59G7ddv8OTQU1tQEsTykgEZG+2rCMQ1+4CXcILeDazgt4eY+JzD5hfE1tAMtTABAR6YOWtg5euedG/o/nHv4SAN88aTQccXTcTdtpSgGJiGxHvvb/Hls3Fp79u9UDHuw8IO6m9YsCgIjIduRr/88IHsYdstHqn6seq70bv8UUAEREtmPjpi1cFPyahDlmYDh72WY21uDSz2IKACIi23Hcu1/g+GA57rnnvockWJptrMmln8UUAEREetHS1sHJnf9DQG72n8W4Mzyq5mr/l6MAICLSg5a2Du5YeBcn+8OY5Wb/GQLuS01j9um1Vfu/HAUAEZEezFm8lpP9AQKywLbZ//ohB9b84A8KACIiPRr5xirOCh7sMvtfGB5R8zd/8xQARER6cPyw9aSiuj/52f9KH1fzN3/zFABERHpw6Jh3Fer+dJJiYXhEzdb9KUelIERESrS0dfCb/76bH279TxwILVHzdX/KUQAQESmSL/vw734vQZB76pc5XDDxXXzz7Nqt+1OOUkAiIkXyZR/OLLr5G5Lge8+OjLtpFacAICJSZOOmLZwePEyyZOnnfW9+IOaWVZ4CgIhIkePe/QIzgtZuSz8Hy8qfYgoAIiJFPvfBl7st/RwMZR/K0U1gEZFIS1sHz677KxPYtvTzod2PZfZJtV/2oRwFABERttX9uSVxc/TQlwSzszP5xEnTB+XgD0oBiYgAudU/031JIf0DMDz7JnMWr42xVQNLVwAiUtda2jqYs3gtR7x5D2elui79XJptHDR1f8rRFYCI1K38pq8j3ryHb6RuJGG541kYdHV/ylEAEJG6NWfxWqaH93YZ/HPP/A0GXd2fcpQCEpG61NLWwcg3VvG1hpuIxv5o8De+lL5w0NX9KSfWKwAzO9HM1prZejO7Is62iEj9yKd+8jt+83n/LMbV6X/g4XefwiNXHD2oB3+I8QrAzALgP4DjgJeAx8xskbs/HVebRGRwy9/w7di0hXMS9/P3qfu73PS9Jv1p7g6OZ/YgTvsUizMFNAVY7+5/BDCz+cB0QAFARCouP+tvzLTzpeQ9HB8sL6R+ssD8cBrzs8dww4zBuemrnDgDwGhgQ9HPLwGHxdQWERnEWto6+NcFq5hhS/h6w00EUdoHut70HT1iaN0M/gDm7vF8sNkM4AR3vyj6+QJgirv/S8n7ZgGzAEaOHHno/Pnz+/3ZmzdvZvjw4f0+T5zUh+qgPlSH3vrw6MY0Nz/Vyencz3Wpn5LASwb/XN5/Icdw4UENfHRUahe2fJuB/HuYNm3aCndvKj0e5xXAS8CYop/3BTaWvsnd5wJzAZqamry5ubnfH9za2kolzhMn9aE6qA/Voac+tLR1cOPiVcywroN/ft6bz/s//O5T+HbMK37i+HuIMwA8Bvytme0HdADnAOfG2B4RGUTyOf8ZtoRvpG7EoGi1D9wXNnGLTefsGWfwzTpK+xSLLQC4e8bM/hlYDATAT919dVztEZHBoXilz3mJ+/ha6qaSTV65lM+dfizXnzWxrnL+pWLdCObuvwF+E2cbRGTwyM/6t6RDLg9u45LkPYV8f9594aHcHRzP9afXz2qfnmgnsIgMGnMWr2VLOmR2ci7nBK2F4/mcf5qAn2Q/yew6WurZGwUAERk0Nm7awqXBHYXBv1zO//wZZ2jwjygAiMig0NLWwZREO/+cvBugbHmHywZ5bZ8dpQAgIjXv0Y1pfn7/k9wa3N7thu+12YuYOuMLdbvSpzcqBy0iNe+udWlmhguZnFiP+7baPl/K/COTT/+CZv09UAAQkZr3ga1ruTx1B5BL/eRr+9weDv6Knv2hFJCI1Kz8mv/rkgvLPtBlMD/NqxJ0BSAiNSm/5v8DbyzjyMQTXVM/6QtpTzYO6qd5VYKuAESkJuXX/F+TurUw+w/dmB9O4+F3nzLon+ZVCQoAIlJzWto66Ni0hauTP6cxsaGw0SsTpX4eueLoeBtYIxQARKQmFNf4MWBWYhEXBf8N5G78hg53hkfx8h4T421oDek1AJjZoj6c43V3v7AyzRER6aqlrYOvLFrNpi1pJts6Lkk+zIdZz8HBC0Dxhq+AX1uz8v47YHtXAI3ARb28buSe6ysiUjGls/1Jto5ZyV9xXLCCBNsG/eLdvl9KX8jZKvOwQ7YXAK529wd7e4OZXVvB9ohInSuu6DnZ1nFxsIhjgpWFgT+vXKkH7fbdMb0GAHdfsL0T9OU9IiJ9kX9270TW8plkC0cHjxce5JJX/BTb/BO9FnIM31bqZ4f16SawmTUBVwMfiH7HAHf3Dw9g20SkDhSneybbOn4a3MURwZM9DvwhxvJwPOsZzcLwCF7eYyIXvj9U6mcn9HUV0DzgMuBJcrusRUT6pfTm7pzkfKYGa8iP+cUPbofcwL8kPJS54Sk8N/RAvvzJA1kYDfqtra27vP2DQV8DwKvu3pcVQSIivcoP/Pu/s5ovBg/RmHqRSYn1JIpu7EL5gf/lPSaqpHMF9TUAfNnMfgLcD2zNH3T3hQPSKhEZVErTPN8OFnFMw8ouK3oo+T4/8P8k+0nOnzGjMNuXyulrAPg0cACQYlsKyAEFABHpUfFsf3qinb2DTVyYXNxtRQ+Un/G3Jxv1+MYB1NcAMNHdDx7QlojIoFA62/9mcA/HNqwkQbbHG7tpAn4bTuI19mBheAQrfRyjRwxVPZ8B1tcAsNTMJrj70wPaGhGpWeV27J4dPEAyShoUp3fKzfZX+jj23D3V5eauDKy+BoCPAzPN7Dly9wC0DFREetmxu5IEXna2n3GABEvCyV0G/hs+eaBm+7tYXwPAiQPaChGpOaU7ds8MHuSsoJWgaOAv3bQ1N/MJNjOMpdnGQprnBqV5YtOnAODuLwx0Q0SkdrS0dXDrnXfyaXs6urF7b9kZfxjN9Itz+0rzVI/tVQNd6e6T+/seEal9xemecxL3c3vqJoKSG7vlyjTckT0GB832q9B2q4Ga2RO9vG7AHhVsj4hUmeKlnJcED7N/ciNTg/bCwF98Y7d0xv/c0AP5rnL7VWt7AeCAPpwjrERDRKS6lO7YPaehNTfjLynFDLkbu07Al9IXMj97DKNHDNWO3RqwvWqgyv2L1JHSNfyXBw8yo+FBktHAXyyMHsLyX5mT2MwwlvkEzp8xQyWZa4geCSkiXdbwA3w68RuuSc3rcmO3OL+fJmBB2Fy4sTs0FTD7dO3YrTUKACJ1qutsfwHnJtrxwDk18SgHJF7qMvBnMTIEPBAeUsjvt/m4ws1dpXtqU1+fB9BtF7CZNbt764C0SkQGTLkdu2cFrSQJKcnyFG7szg+nFWb7AHvuntLN3UGgr1cAC8zs58C3gd2ir03A4TvzoWY2B/gk0Ak8C3za3TftzLlEZPvK7dg9PfkQZwcP5gb+otl+/uZuPsefv7ELFNbwa+AfHPoaAA4DvgU8CryL3ANiPtaPz70PuNLdM2b2LeBK4N/6cT4R6UHxjl2A8xL3cm3qlm75fSe3rjvjRkjAneFRXQqzKc0z+PQ1AKSBLcBQclcAz7n7Tj8ZzN3vLfpxKXDmzp5LRHqWf8Zu6M4RtopLk3dySOKPJfl9yJDkzvAonsqOZS/bzNJsY7enbsngY158a7+nN5mtAu4GvgbsDfwnkHb3fg/cZvYr4A53v7WH12cBswBGjhx56Pz58/v7kWzevJnhw4f3+zxxUh+qQzX24dGNae5al+Z/38n9v31K4lH+Jfgl4xId3dbwl8vv772bcca4FB8dlYqzGzukGv8edtRA9mHatGkr3L2p9HhfA0CTuy8vOXaBu/+8l99ZAry3zEtXu/vd0XuuJncv4XTvQ0Oampp8+fLl23vbdrW2ttLc3Nzv88RJfagO1dSH4o1bUxNPs7e9yQmJZeybeL3bE7eydM3v13puv5r+HnbWQPbBzMoGgL4Wg+s26vY2+EevH7udBs0ETgGO6cvgLyLd5W/ujnxjFVMT7XzSd+NLDbeSim7supfe2M0t58zn91/eYyKz3h9y1bnHxd0ViUEs+wDM7ERyN32Pcve342iDSC0rrc9zVsODJMl0K8xWbuBvyR7Z5Rm7ra2tsfVD4hXXRrAfAEOA+yz3r3Wpu18cU1tEakbX+jwPc1bDtvX7ZpB1CD1BrjKPd1vRo2fsSrFYAoC7fyiOzxWpRaX1eS4LHopm/F3X72cdttLAtekL2Ms287oPL6zo0TN2pRyVghCpUqU7dj+X/C2nB78rVOSEbev309EyzuLVPIAeviK9UgAQqSKls/1zE+0kgwyfTf6y28Cfr89TPPDnSzlo45b0hQKASBXoXp+na5mG4hu65QZ+UJkG2XEKACIxKTfb35wYwjWpeYVlnJCvyZPAoo1bGvilUhQARHax0tr7k20dtzVcRwPpbss4s0Bn0Y3dpdlGlWGWilEAENkFSqtx5nY+OucmlvCF5F0MIV2yjLPn2b7KMEulKACIDKDS2T5Ak7UzM7iXQxLPdCnTELrRSUqzfdllFABEBkhpGeZma+PzybuYGFXjzEZr9xOWK8H8SPYgvpc5Q7N92WUUAEQGQL4M80TW8snk7zmIP9IUPFNYyQO51TxOgqw7aZKFwV+zfdlVFABEKqS4MNvhiaf5cvA65wX3E5h3eeAK5NI96aJ0z+qGg/nUWadrw5bsUgoAIv1UXJ/ns8EDnNHwcGHTViG/XzTbL765+/IeE7nshPFcroFfYqAAILITuq7qeZwTE3/g+w3f77JbN/9MXYtSPMWz/b87TbN9iZ8CgMgO6Lqqx7kguI9PBffyQdtYWMOf361buqJnFeO5fsZEzfalaigAiGxH6Rr+KfY0M1OLmWTreV/iL9GjFXNllwPPll2/PzQVcP3pKsMs1UUBQKQHpfV5zg1WMt5e5Jjg8W7LOHHjzvAoNvo+hfLLKswm1U4BQKRIufo8Q4N3uCT5KwKyXd5buoxT9Xmk1igAiNB9x26TreG2hm+QItOlGmfGc8/aKr2xW/zQFc32pVYoAEhdKx34P25P8A/J/2ZKYg0NlgF6Xs2TH/T10BWpVQoAUndKb+pOsnWclFzGgTzH4UE7ZrlBP+0BRvdBv7g+z6z3h1x17nFxd0lkpygASN0o3rA1PdHOUmvkw4ln+VLy1m67dZ0Ed4TNXW7qQvf6PK2trbH1R6S/FABkUCs325/X8I1C7f1ENPBD1926+Zu6qsYpg5kCgAxK5R66cmxiBUcknmA3Ogs3dh/JNDI5eJaUZ7qlep4beqCqccqgpgAgg0a52f7UoJ3deZuLk78uLOPMYuDQSYrrw7MhhKmJ9m43dTXwy2CnACCDQmnt/SZrZ17D7DLLOBPMD6d127C1Mswt4bxBaR6pIwoAUvOKa+83B4+zp21mevBIj8s4tWFLJEcBQGpSce39qYl2PpPI8LnkLwvVOJ8J38eQRJqArDZsifRAAUBqStf6PGu5veE6GrqleYxfZo9kaaaxbG5fG7ZEchQApOqV1uf5VOJJhiXf4czEQwzpkuZJYA5pkoUNW8rti/RMAUCqVunGrWwC/jV1J8kozbMhuw/v8i0kyqR5VjGe7549UYO+SC8UAKTqlJZhvr3h67nVPNHr+TTP7eHRLM1O6JLmAdXeF+krBQCpCqVpnvMTqwmCkHOCBwppnqxD2CXNM4GVUZpHtfdFdpwCgMSqdMdusz3OfzVcT5IQM3gl+246PSBRpiiblnKK9E+sAcDMvgjMAd7j7q/F2RbZdUofqD7J1nFK8veM4RWmBatIWm7HbujGzeGJ3dI8uWJtmu2L9FdsAcDMxgDHAS/G1QbZtbpX4zyApsQ6/i05v1CN8w/hARwSPEuSsFuaBzTbF6mkOK8AvgtcDtwdYxtkgJWvxnkdQ8ilfBLR+n3I5fcf8ol8u/Ocbjd2NfCLVJ55/v++XfmhZqcCx7j758zseaCppxSQmc0CZgGMHDny0Pnz5/f78zdv3szw4cP7fZ441UIfHt2Y5uanOunM5qpxHh88xpH2BI2JDYWNWw+HB/GRYB0pctU4z+u8qjDoA+y9m3HGuBQfHZWKsSc9q4W/h+1RH6rDQPZh2rRpK9y9qfT4gF0BmNkS4L1lXroauAo4vi/ncfe5wFyApqYmb25u7nfbWltbqcR54lTtfWhp6+DGxasIHT6R+D3fT/2ABLnJRhbDHdKkuCE8s6arcVb730NfqA/VIY4+DFgAcPdjyx03s4OB/YBVZgawL7DSzKa4+58Hqj0y8ErTPQEZPhPcw6XJX5DAVY1TpMrs8nsA7v4k8Df5n7eXApLaUFqO+ezEb7k0eSd/k3iDR8JGDk2sJ+mhqnGKVBHtA5B+y5djDt1JkmF28iecGTwEQKcnuT5zNtA1zTN6xFBO1gPVRWIVewBw97Fxt0F2Xn7mH7qzr73K91I/4NDEM2Q9t8In4VmmJtr5YTidVdnxXH/WxEI1Tj1QXSResQcAqV3FD2KZmVzMscFKsiT4TvoMLkn+qvCc3aXZRoamAmarPo9IVVEAkJ2Sn/l/hKe4tWE2ScuSdeMz6c/yP9nD+F324ELK5+U9JjJbN3dFqo4CgOyQ4pU+B9iL/CD1/ULphizG/pZbyLWK8XzqzLO4XIO+SNVSAJA+y8/6GzPtfCW1iKMSq9jMULZ6svDoRaV7RGqHAoD0ST7f/zHauKlhDoE5oRuXdl7CmwwrpHtWMV61+EVqhAKA9Kh0Y9cUa+cHDd8v7Oh1jAmJF/lhOJ2V4Tg9iEWkxigASFnFG7uabA1fSN7F1MRqXvY9GWIZAg8LKR9QaWaRWqQAIN0Ub+w63h7jRw03EJiT8QRfTF/MOwwppHzak43coFm/SE1SAJAuijd2nZBYxg2p/yikfAAOSfyxkPIJzJTyEalhCgACbMv3j3xjFRclnuLA5POcmFzOM9lRjLFXC3V88ikfrfQRqX0KANJleedt0cNazKAl81Euy1zMwfbHbnV8lO8XqX0KAHUun+/f3d/iytRthcE/dGOdjyFNkpU+rlsdHxGpfQoAdag43XN44mk+H2zlnOAB9uZNQhLgKN0jUgcUAOpI/qHsm7akmWxrub3hOhrIYAZrs6P5dPpyGsgo3SNSJxQABrni2f7URDsHZ8cyLniJfwruZohlgFy65+7wYzzl+wMUNnVpeafI4KYAMEh1ne2v4/aGr5MigwFm8HR2DMP9HRKFGj4TCr8bmCnlI1IHFAAGmZa2DhbevZCD0k9yQPZDjAle5dLgzsJsP+twc+YErs3MZLKt65LuAeX7ReqJAsAgUFyzZ7Kt47aG62hIpguz/Rez+9DpAQmcNEl+FR4OwEofx8ow91B2UDkHkXqjAFDDimf778uOZ0JiM1cl57GbpYHcbH9+ZhpXZi5isj3TbbYPeii7SD1TAKgxXWf7Cwqzfcg9g/c1fxdpD7Botn9neBRgmu2LSDcKADWi+KYuOFNsDdembi5s3Mo6/DLzMb6YuZiJ9qxm+yKyXQoAVay0Hv9H7Ulmpu5lvL3IBxKvstmHRBu3nDQpfh4eR0hQmO3naeAXkXIUAKpQ19k+jLcX+XzwC04MlkezfeMHmen8IHMaE+yFbrN9AxyleUSkdwoAVaJ0tt9k7XwqeR8fTGxkQuJFMp4rymzkHr7+tu/GOwzRbF9EdpoCQMxKZ/ujeI1Lk3dyRvBwIbd/U+Z4HggP4T8bbiDlmS51ejTbF5GdpQAQk/zAv/87qzk38TTvBA0cnmjn6MRKDC+a7Sd41ffkIT+E8zqv6pLu0WxfRPpDAWAXKk3zHGGPc2PD9SQJMYNNvjs/Ck/liXA/vtfww26z/bYo3TN6xFBmvT/kqnOPi7dDIlLTFAB2geLZ/vTE07yWeDdTgnWclvgdScsCuYJsN2Y+wffD0wE4r3NEr7P91tbWuLojIoOEAsAAKZ3tH2ar+VnDN0lFs/23vYEl4WSag1UkyT1u8ZHsQYXfL57t36DcvogMAAWACiu9qbuf/YnzgiWcHyyhwUIgN9ufmzmFG8IzmRx2L8im3L6I7AoKABVQOts/1NZwbvJ+9rM/MSn4I2kP+EP2AD6SWEsQlV9+KPthgC7LODXwi8iupADQD6Wz/ffwF76Q/AXnBA+QiJZwzssczQ2ZM3mVEd3KL2sJp4jEKbYAYGb/AvwzkAF+7e6Xx9WWHVE6259k65iRfJAx9gqHJdaQshD33HuzJOjw9/AqIwDN9kWkusQSAMxsGjAd+LC7bzWzv4mjHTuidLY/nLf5bHIhFwX/TcIcd1gUHs7ibBPXp/5TG7ZEpOrFdQVwCfBNd98K4O6vxNSOXpU+aOXcRDt/SuxFU2IdpwW/Y5htLcz2QxKs9ffzm+zh/Llzb23YEpGqZ54fwXblh5o9DtwNnAi8A3zR3R/r4b2zgFkAI0eOPHT+/Pn9/vzNmzczfPjwHl9/dGOaeU938lbuKYpMsXZubZide6auwVYPuDv8OI9lx/HV1C2kyM32z+u8qkv55b13M84Yl+Kjo1L9bvOO9qEWqA/VQX2oDgPZh2nTpq1w96bS4wN2BWBmS4D3lnnp6uhz9wSmAh8BFpjZ/l4mGrn7XGAuQFNTkzc3N+90m7bN6I0RQzsxg7+8nSYwI3RnxNAUnZmQAzJrmJloZ729j4nBH7kwWExD9Ezd0I0fZ07lu+EMAJ7tHB3LbL+1tZX+/LeoBupDdVAfqkMcfRiwAODux/b0mpldAiyMBvxlZpYF9gFeHZDGbFjG6kd/zR1P7U1H+oMAbNqSZrKt4++D3OC9mrEM2fI60xKr+VbDf5EkE9XigeXZ8RySeLawhPPB7MTCqbVhS0RqVVz3AFqAo4FWMxsHNACvDcgnbVgGN53EhGyGeQnjpYa9eYch7MZW9rX/xchddJh1/9Wsw9zMKXwrPLfbEk7QSh4RqW1xBYCfAj81s6eATmBmufRPRTz/MGTDaBWO8za78ayP4oO2EfBCyeWHwoNZnJ3CCP7K55ILC7P9+7K5tNlKH0dbOE4reURk0IglALh7J3D+LvmwsUdAcjcyma2kSXJV+iJW+jgm2zrmNXyjsFzze5kzCjP7pdkJmu2LyKA3+HcCj5kCMxex9tFf8/Wn9mal5+4BrPRx3err52nDlojUg8EfAADGTOHAs6dw9rgOXozW9Y8YmuI5O5CVb48jsFyCaMTQVK4u/9tpRinNIyKDXH0EgMhpk0Zz2qTRg2LJmIhIfyXiboCIiMRDAUBEpE4pAIiI1CkFABGROqUAICJSp2KpBrqzzOxV4IUKnGofBqr0xK6jPlQH9aE6qA+9+4C7v6f0YE0FgEoxs+XlSqPWEvWhOqgP1UF92DlKAYmI1CkFABGROlWvAWBu3A2oAPWhOqgP1UF92Al1eQ9ARETq9wpARKTuKQCIiNSpugsAZnaima01s/VmdkXc7emJmf3UzF6JnpqWP7aXmd1nZs9EX/cseu3KqE9rzeyEeFq9jZmNMbMHzKzdzFab2eei47XUh93MbJmZrYr6cG10vGb6kGdmgZm1mdk90c811Qcze97MnjSzx81seXSs1vowwsx+YWZrov8vDo+9D+5eN3+AAHgW2J/cc4hXARPiblcPbT0SmAw8VXTs28AV0fdXAN+Kvp8Q9WUIsF/UxyDm9r8PmBx9/y5gXdTOWuqDAcOj71PAH4CptdSHor5cCtwG3FNr/5aidj0P7FNyrNb6cAtwUfR9AzAi7j7U2xXAFGC9u//Rc4+lnA9Mj7lNZbn7Q8DrJYenk/tHRPT1tKLj8919q7s/B6wn19fYuPuf3H1l9P1fgXZgNLXVB3f3zdGPqeiPU0N9ADCzfYGTgZ8UHa6pPvSgZvpgZu8mN6m7EXKPxXX3TcTch3oLAKOBDUU/vxQdqxUj3f1PkBtggb+Jjld1v8xsLDCJ3Ay6pvoQpU4eB14B7nP3musDcANwOZAtOlZrfXDgXjNbYWazomO11If9gVeBm6JU3E/MbBgx96HeAoCVOTYY1sFWbb/MbDhwF/B5d3+zt7eWORZ7H9w9dPdDgH2BKWZ2UC9vr7o+mNkpwCvuvqKvv1LmWOx/D8DH3H0ycBLwT2Z2ZC/vrcY+JMmldH/k7pOAt8ilfHqyS/pQbwHgJWBM0c/7AhtjasvOeNnM3gcQfX0lOl6V/TKzFLnBf567L4wO11Qf8qLL9VbgRGqrDx8DTjWz58mlPI82s1uprT7g7hujr68AvySXDqmlPrwEvBRdQQL8glxAiLUP9RYAHgP+1sz2M7MG4BxgUcxt2hGLgJnR9zOBu4uOn2NmQ8xsP+BvgWUxtK/AzIxcvrPd3b9T9FIt9eE9ZjYi+n4ocCywhhrqg7tf6e77uvtYcv/ef+vu51NDfTCzYWb2rvz3wPHAU9RQH9z9z8AGMxsfHToGeJq4+xD3nfFd/Qf4BLkVKc8CV8fdnl7aeTvwJyBNbjbwj8DewP3AM9HXvYref3XUp7XASVXQ/o+Tu2R9Ang8+vOJGuvDh4G2qA9PAf8eHa+ZPpT0p5ltq4Bqpg/k8ueroj+r8//f1lIfojYdAiyP/j21AHvG3QeVghARqVP1lgISEZGIAoCISJ1SABARqVMKACIidUoBQESkTikAiIjUKQUAkZ1gZmPNbEtUJ2hHfu/sqMTvPQPUNJE+UwAQ2XnPeq5OUJ+5+x3ARQPTHJEdowAgUsLMPmJmT0QPhBkWPQymtyJw+SuCNVGVx6fMbJ6ZHWtmj0QP+6jWkspSx5JxN0Ck2rj7Y2a2CPg6MBS41d2f2s6vAXwImAHMIld36lxyJTFOBa5iW613kaqgACBS3lfJDeLvAJ/t4+885+5PApjZauB+d3czexIYOyCtFOkHpYBEytsLGE7ucZa79fF3thZ9ny36OYsmW1KFFABEypsLfAmYB3wr5raIDAjNSkRKmNmngIy732ZmAfComR3t7r+Nu20ilaRy0CI7IXrO8T3u3uvqoB5+txn4orufUuFmiewQpYBEdk4I7LEzG8GAHwJ/GYhGiewIXQGIiNQpXQGIiNQpBQARkTqlACAiUqcUAERE6tT/B+jn1N/5RulZAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "zgr = np.interp(xgr, x, bathy)\n", "\n", "plt.figure()\n", "plt.plot(x,bathy,'-o')\n", "plt.plot(xgr,zgr,'.-')\n", "plt.xlabel('x [m]')\n", "plt.ylabel('x [m]')\n", "plt.title('xb bathy')\n", "plt.grid('on')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Seaward extend\n", "Extend the grid to the required offshore depth. In this case 20 m with a artifial slope of 1/20.\n", "\n", "If desired the function `offshore_depth` can also be used to determine the required offshore water depth." ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "Text(0, 0.5, 'z [m]')" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYcAAAEGCAYAAACO8lkDAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAAAhuklEQVR4nO3dd3hU55n38e+NhASiSKgAoqpRTDFNprpX4vh1SdaJC+CsHYPbOtVOHL95N7tZb9brbLxJnNjGJY4xuCS4EDuOE2zHMWCKRC8GVCiiqiAkBBKS5nn/mDEWlgABGp0pv8916WLmnDOaW0fD/PSc+8xzzDmHiIhIUx28LkBEREKPwkFERJpROIiISDMKBxERaUbhICIizcR6XUBbSE1NdRkZGV6XISISVvLz88ucc2ktrYuIcMjIyCAvL8/rMkREwoqZbT/ROh1WEhGRZhQOIiLSjMJBRESaUTiIiEgzCgcREWlG4SAiIs0oHEREpBmFg4hIiCg/VMfvl2xjSUEZtfWNntYSER+CExEJd7X1jdz41CcUldUAEBfbgbEDkpicncqUnBRG9UsiNqb9/p5XOIiIhICXl++gqKyGp6aNJT42hiWFZSwpLOfxhVv4xd+ge6dYLhiUxkWD07hwcBq9EzsFtR6Fg4iIx2rrG/nt3wuZmJXM1BHpAFwytCcAB2qOsqSwnI+27OejLaW8s24PAEN7d+OiwWlcMawXuRnJbV6TwkFExGML1uymtLqOX988ptm6Hl3i+PK56Xz53HScc2zeV81Hm0v5aEspzy8uZm9VrcJBRCQSLdy4jz6JnZiQefI3eTNjaO/uDO3dnVkXZVNT10B1bUNQatLZSiIiHqqtb2RRQRmXntMTMzutx3aJjw1a70HhICLioWXFFRw+2shl5/TyupTjKBxERDz0/qZ9dO4Yw6SsFK9LOY7CQUTEI4ePNvDmql1cdk5POnWM8bqc4ygcREQ8Mn/lLqpqG/jnKRlel9KMwkFExAM+n+OFxcWc2y+RsQN6eF1OMwoHEREP/GNrKYWlNdw+JfO0z1JqDwoHEREP/G7xNnp2i+fqkelel9IihYOISDsr2H+Ij7aUMm3iQOJiQ/NtODSrEhGJYC8sKSYutgO3TBjgdSknpHAQEWlHZYfqmJ+/i+tG9SG1a7zX5ZyQwkFEpB3978It1Df6uOvibK9LOSmFg4hIO9m6r5qXl+/k1gkDyE7r6nU5J6VwEBFpBw2NPr7/x7V0jY/l/ssGeV3OKWnKbhGRdvCbDwtZs7OSJ24ZQ0oI9xo+o5GDiEiQrdlZya8+2Mr1o/twzbl9vC6nVRQOIiJBdORoI995dTW9usXzb9eN8LqcVtNhJRGRIPqvdzdRVFbDvG9OILFzR6/LaTWNHEREgmTR1jJ+/8l2bp+SyeScVK/LOS0KBxGRIDh4pJ4H/riG7LQuPDh1iNflnDYdVhIRaWPOOf7vm+vZX13H63dPDrkL+bRGyI4czGyqmW02swIz+6HX9YiItEZFzVHue3kVf1qzm+9eMZhR/ZO8LumMhOTIwcxigN8AVwAlwAozW+Cc2+htZSIiJ7ZlXzUznltOeU0dD1w1hHtCfIqMkwnJcADGAwXOuSIAM3sFuA5QOIhIyHHO8c66Pfzo9XV06hjDG/dMYUTfRK/LOiuhGg59gZ1N7pcAEzyqRUSkRY0+x7vr9/CbDwvZtKeKEX278+St4+ifnOB1aWctVMOhpWvmueM2MJsJzAQYMCB050QXkchz5Ggjb6zaxTMfF1FcVkNWahf+58ZRXDe6D7ExIdvKPS2hGg4lQP8m9/sBu5tu4JybDcwGyM3NPS44RESCoaaugV9/UMC8Zdupqm1gZN9Enrx1LFcO701Mh9C7DvTZCNVwWAEMMrNMYBdwE3CLtyWJSDSrqDnKtGeXsWlvFVePTOcbkzPIHdgDs8gKhc+EZDg45xrM7D7gPSAGeN45t8HjskQkSpVW1zHt2WVsK6/h+dvO45KhPb0uKehCMhwAnHN/Bv7sdR0iEt32HDzCtGeXsavyCM9/4zymhNk0GGcqZMNBRMRr+dsrmDVnJbX1jfz+n8czISvF65LajcJBRKSJkgOH+XhrGR9vLWXhxv2kJ3Xi5TsnMKhXN69La1cKBxGJatW19SwtquDjraUs2lpGUVkNAL27d+Kr4/rxg6lDSEqI87jK9qdwEJGo0uhzrC2p5OOtZSzaWsbKHQdo8Dk6d4xhYlYy0yYO5IJBqeT07BqxZyK1hsJBRCLezorDLCrwHypaXFDOwSP1mMGIPonMvDCLCwalMXZgEvGx4Td7arAoHEQk4jQ9VPTx1jKKA4eK0hM7cdXwXpw/KI0p2SmkdI33uNLQpXAQkbDX9FDRx1tLWbWj8tihoknZKUyfOJALB6eSnRbdh4pOh8JBRMLSzorPzypaXFBGVW2DDhW1IYWDiISF6tp6Piks9zeSC44/VDR1RG8uGJTGlJxUkrtE35lFwaBwEJGQdLTBx6odB1hcWM6SgjJW7ayk0edIiIthYlYKMyYN5IJBaWSnddGhoiBQOIhISPD5HJv2VrG4oIzFBeUsL67gSH0jHQxG9k1k1oVZXDg4jbEDehAXGxnTYocyhYOIeMI5x/bywywuLGNJQTmfFJVTUXMUgOy0Lnwttx+Tc1KZmJVCYueOHlcbfRQOItJuSqvrWFJYdmx0sKvyCOD/NPIlQ3oyJSeFydmp9E7s5HGlonAQkaCprq1nWVHFsdHB5n3VAHTvFMvk7FTuuiiLyTmpZKWqbxBqFA4i0mbqGhpZub3y2OhgTclBGn2O+NgOjM9M5voxfZmSk8LwPokRd+W0SKNwEJEz1uhzbNxdxeJAGKzYVkFtvY8OBqP6J3H3RdlMzklh7IAedOqozxuEE4WDiLSac47isppjp5cuKfTPUwQwuFdXbjpvAFNyUpmQlUz3TmoihzOFg4ic1P6q2sDIoJzFBWXsOVgLQN+kzlw5rBdTclKZnJ1Cz+5qIkcShYOIHKeqtp6lheUsKSxnUUEZBfsPAZCU0JHJ2Sncm53K+TmpDExJUBM5gikcRKJcbX0jK7cfYHFhGYsKyllXUonPQeeOMZyXmcyN4/oxJSeVYend6aAmctRQOIhEmUafY/2ug8eayHnbDlDX4COmgzG6fxL3XZLDlJxURg/QpHXRTOEgEuGccxSW1rCk0H/ls6VF5VTVNgAwtHc3pk0cyJScFMZnptA1Xm8J4qdXgkgE2nuwNvAp5DIWF5axr6oOgH49OnP1yHQmB5rIqbrYjZyAwkEkAhw8XM8nReXHwqCo1D+ddXKXOCZnpzAlJ5Up2akMSEnwuFIJFwoHkTBUW99I3rYDLCooY0lhGet3HcTnICEuhgmZydwyfgCTs1MZ2rubmshyRhQOImGgodHH2l0HWRKYsC5/xwGONviI7WCMHdCD+y8bxJScVEb1S9J01tImFA4iIcg5R8H+QywKhMGyonKq6/xN5GHp3blt0kAm56QyPiOZLmoiSxDoVSUSInZVHmFxQdmxaSn2V/ubyANTErhmVB+m5KQwKSuFFDWRpR0oHEQ84pwjb/sBFqzefdw1kVO7xjE5O/XYtQ36J6uJLO1P4SDSzg7UHOWVFTt5LW8nxWU1JMTFMCkr5djnDYb06qZpKcRzCgeRdrJlXzW/W7yN11eWUNfgY3xGMvdcnM3VI9PVN5CQo1ekSBD5fI6PtpTy/OJiPt5aRnxsB74yti/fmJzJkN7dvC5P5IQUDiJBUFPXwOsrS/jd4m0UldXQq3s8D1w1hJvHDyC5S5zX5YmcksJBpA2VHDjMi59s55XlO6iqbWBUv0R+edNovjQiXZ8/kLAScuFgZj8B7gRKA4t+5Jz7s3cViZycc44V2w7wwpJi/rJ+L2bG1BG9uX1KJmMHJKm5LGEp5MIh4HHn3M+9LkLkRJxzbN1/iLdW7+Kt1bspOXCExM4dufPCLGZMyqBvUmevSxQ5K6EaDiIhxTnHtvLDLCsqZ2lROcuKK9hzsJYOBlNyUvnO5YP50sjeJMTpv5REhlB9Jd9nZjOAPOB7zrkDX9zAzGYCMwEGDBjQzuVJpHPOUVRWw7KiikAYlB+b9jq1axwTMlOYmJ3CVcN70bObrp0skcecc+3/pGYLgd4trHoYWAqUAQ74KZDunLv9ZN8vNzfX5eXltXmdEj38F8Q5xNJjYVBBaWD6irRu8UzMSmFCZjITs1LITuuiPoJEBDPLd87ltrTOk5GDc+7y1mxnZs8Abwe5HIlCn/UM/IeJKlhWXE7ZoaMA9Ooez+TsFP/oICuZzFSFgUSfkDusZGbpzrk9gbs3AOu9rEcig8/n2LK/+thhouXFFZTX+MMgPbETFwxKOzYyGJiSoDCQqBdy4QD8t5mNxn9YaRswy9NqJCz5fI5P91azrLj8WBgcOFwPQN+kzlw0JI2JmSlMzEqhf3JnhYHIF4RcODjnpntdg4SfRp9j054qlhV/PjI4eMQfBv16dObSob2YmJUcCAPNcipyKiEXDiKt0ehzbNxddexMouXFFVTV+i+GMyA5gauG92JCZgoTspLp10NhIHK6FA4SFhoafWw4FgYVrCiuOHZltIyUBK4emc6ErGQmZKbQRx9AEzlrCgcJOT6f40h9I1v2VR87TJS37QCHAmGQldqFa0b1YWIgDHon6nMGIm1N4SDtrra+kUfe2cTCTfuOW17f6KOmrpEj9Y3HLc9O68J1o/swISuFiZnJ9OyuMBAJNoWDtKu9B2uZ9VI+a3ZW8qURvenW6fOXYGxMBxI6xpAQH0tCXAz9enRmfGayPoEs4gGFg7SbFdsquPullRw+2sBT08YydUS61yWJyAkoHKRdLFizm+++upp+PToz784JDO6lq6CJhDKFgwTdmp2VfP8Paxg7sAfPzMglsXNHr0sSkVPQpakkqEqr65g1J5+e3eJ5eto4BYNImNDIQYLmaIOPe+bmU3nkKK/fPYUeunaySNhQOEjQ/PTtjazYdoBf3TyGYX26e12OiJyGk4aDmVWd4vEG7HHODW67kiQSvLpiB3OWbmfWRVlcO6qP1+WIyGk61cih0Dk35mQbmNmqNqxHIsDKHQf48ZsbuGBQKg9eNdTrckTkDJyqIf3VVnyP1mwjUWJ/VS13zcmnd2Infn3zGGI6aCpskXB00pGDc66o6X0z6970Mc65ii9uI9GrrqGRu17K51BdAy/eMZ6kBDWgRcJVqxrSZjYL+HfgCP6L8BD4NytIdUkY+smCjazcUclvbx3L0N5qQIuEs9aerfR9YLhzriyYxUj4mrtsOy8v38E9F2dz9UhNiyES7lr7IbhC4HAwC5Hwlbetgp8s2MDFQ9L43pVDvC5HRNpAa0cODwFLzGwZUPfZQufc/UGpSsLG3oO13D13JX2TOvPLm9SAFokUrQ2Hp4EPgHWAL3jlSDiprW9k1kv5HK5rYO43J2hqDJEI0tpwaHDOfTeolUhYcc7x/95az5qdlTw1bZxmWRWJMK3tOXxoZjPNLN3Mkj/7CmplEtJeWrqd1/JKuP/SHKaO6O11OSLSxlo7crgl8O9DTZbpVNYotayonH/700YuG9qTb1+umVNEIlGrwsE5lxnsQiQ87K48wr3zVjIgJYHHbxpNBzWgRSLSSQ8rmdnYU32D1mwjkaG23v8J6Np6H7On59K9kxrQIpHqVCOH35nZxfhnXz2R54CTTs4n4c85x8NvrGdtyUGemZFLTs+uXpckIkF0qnBIBPI5eTiUtl05EqpeWLKN+StL+Pblg7hiWC+vyxGRIDvVxHsZ7VSHhLBPCsv5j3c2ccWwXtx/6SCvyxGRdqBrSMtJlRw4zL3zVpKRksAvvjZKDWiRKKFwkBM6crSRWXPyqW/w8cyMXLqpAS0SNXQNaWmRc46HXl/Lxj1VPHdbLllpakCLRJNWjRzM7H0zu/oLy2YHpyQJBc8tKubN1bv53hWDuXSoGtAi0aa1h5UygR+Y2b82WZYbhHokBCwuKOM//7yJqcN7c+8lOV6XIyIeaG04VAKXAb3M7E9mlng2T2pmN5rZBjPzmVnuF9Y9ZGYFZrbZzK46m+eR07ez4jD3zVtJTs+u/PxrozBTA1okGrU2HMw51+CcuweYDywCep7F864HvgL847gnMRsG3AQMB6YCvzWzmLN4HjkNR442MnNOPo0+x+zpuXSNV0tKJFq1Nhye+uyGc+4F4BvAX8/0SZ1zm5xzm1tYdR3winOuzjlXDBQA48/0eaT1nHM8OH8tn+6t4lc3jyEjtYvXJYmIh1oVDs65p79wP985d3sQ6ukL7GxyvySwrJnAFOJ5ZpZXWqoPaZ+t2f8o4k9rdvPAVUO4eMjZDApFJBIE7biBmS0EWpro/2Hn3FsnelgLy1xLGzrnZgOzAXJzc1vcRlrnH1tKefQvn/LlkencfVG21+WISAgIWjg45y4/g4eVAP2b3O8H7G6biqQl28tr+JeXVzG4Vzceu/FcNaBFBAi9T0gvAG4ys3gzywQGAcs9rili1dQ1MGtOPgCzp+eSEKcGtIj4eRIOZnaDmZUAk4B3zOw9AOfcBuA1YCPwF+Be51yjFzVGOuccD/xxDVv2VfPELWMYkJLgdUkiEkI8+VPROfcG8MYJ1j0CPNK+FUWfJz8q5M/r9vKjq4dywaA0r8sRkRATaoeVpB18uHk/j723mWtH9eHOC3QZcBFpTuEQZYrLavjWy6s4p3d3Hv2qGtAi0jKFQxQ5VNfAzBfziOlgPD19HJ3j9OFzEWmZTk+JEj6f43uvraaorIYXbx9P/2Q1oEXkxDRyiBK/+bCA9zbs46EvDWVKTqrX5YhIiFM4RIH3N+3jFwu3cMOYvtxxfqbX5YhIGFA4RLjC0kN8+5XVDO/TnZ99ZaQa0CLSKgqHCFZdW8/MF/OIi+3A09Nz6dRRDWgRaR01pCOUz+f4zqtr2FZ+mJfumEDfpM5elyQiYUQjhwj1y/e3snDTPn785XOYlJ3idTkiEmYUDhHorxv28sv3t/LVsf24bXKG1+WISBhSOESYgv3VfOfV1Yzql8gjN4xQA1pEzojCIYIcPFLPnS/m0zkuhqemj1MDWkTOmBrSEcLfgF7NzorDzLtzIumJakCLyJnTyCFCPL5wCx98up9//T/DGJ+Z7HU5IhLmFA4R4C/r9/DrDwr4em5/pk0c6HU5IhIBFA5hbvPear772hpG90/i368frga0iLQJhUMYO3i4nplz8ugSH8tT08YRH6sGtIi0DYVDmGr0Oe5/ZRW7K4/w1LSx9E7s5HVJIhJBdLZSmPr5Xzfz0ZZS/vOGkYwbqAa0iLQtjRzC0Ntrd/Pk3wu5efwAbpkwwOtyRCQCKRzCzKY9VTzwh7WMG9iDn1w7zOtyRCRCKRzCSOXho8yck0f3zrE8eetYNaBFJGjUcwgTDY0+/uXlVew7WMcrsybSs7sa0CISPAqHMPHYe5v5eGsZj351JGMH9PC6HBGJcDqsFAbeWr2Lp/9RxPSJA/n6eWpAi0jwKRxC3IbdB/nB/LWcl9GDH1+jBrSItA+FQwirqDnKzBfzSeocx29vHUdcrH5dItI+1HMIUQ2NPu6bt5LSQ3X8YdYk0rrFe12SiEQR/Skaon727qcsKSznketHMKp/ktfliEiUUTiEoDdWlfDcomK+MTmDG3P7e12OiEQhhUOIWb/rID+cv44Jmck8/OVzvC5HRKKUwiGElB2qY+aLeaR0ieM3t46lY4x+PSLiDU/efczsRjPbYGY+M8ttsjzDzI6Y2erA11Ne1OeF+kYf985dSXnNUZ6enktqVzWgRcQ7Xp2ttB74CvB0C+sKnXOj27cc7z3yziaWFVfw+NdHMbJfotfliEiU8yQcnHObAF3SMuAPeTt5Yck27jg/kxvG9PO6HBGRkOw5ZJrZKjP7yMwuONFGZjbTzPLMLK+0tLQ962tTa3ZW8vCb65mcncJDXxrqdTkiIkAQRw5mthDo3cKqh51zb53gYXuAAc65cjMbB7xpZsOdc1Vf3NA5NxuYDZCbm+vaqu72VFpdx6w5+aR1jeeJW8YSqwa0iISIoIWDc+7yM3hMHVAXuJ1vZoXAYCCvjcvz3NEGH/fMzafyyFHm3z2Z5C5xXpckInJMSP2pamZpZhYTuJ0FDAKKvK0qOH769kZWbDvAo189l+F91IAWkdDi1amsN5hZCTAJeMfM3gusuhBYa2ZrgD8CdznnKryoMZheXbGDOUu3M/PCLK4b3dfrckREmvHqbKU3gDdaWD4fmN/+FbWflTsO8OM3N3DBoFQevGqI1+WIiLQopA4rRbr9VbXcNSefXonx/PrmMWpAi0jI0rtTO6lraOSul/Kprm1g9vRckhLUgBaR0KXrObSTnyzYyModlTxxyxjOSe/udTkiIielkUM7mLtsOy8v38HdF2dzzbl9vC5HROSUFA5Blretgp8s2MBFg9P4/pVqQItIeFA4BNHeg7XcPXclfZI686ubxhDTQXNJiUh4UM8hSD5rQNfUNfDSHRNITOjodUkiIq2mcAgC5xw/fnM9q3dW8tS0sQzp3c3rkkRETosOKwXBS0u381peCfddksPUEelelyMictoUDm1seXEF//anjVw6tCffuWKw1+WIiJwRhUMb2l15hHvm5tM/OYHHvz5aDWgRCVvqObSR2np/A7q23scrM8eR2FkNaBEJXwqHNuCc4+E31rO25CCzp48jp6ca0CIS3nRYqQ38fsk25q8s4VuXDeLK4S1d/E5EJLwoHM7SJ4Xl/PSdTVx+Ti++ddkgr8sREWkTCoezsKvyCPfOW0lGSgKPf30UHdSAFpEIoXA4Q7X1jcyak0d9g4/ZM3Lp1kkNaBGJHGpInwHnHA+9vo4Nu6t4dkYu2WldvS5JRKRNaeRwBp5bVMwbq3bxncsHc9k5vbwuR0SkzSkcTtPigjJ+9u6nXDW8F/ddkuN1OSIiQaFwOA07Kw5z37yVZKV24X++NloNaBGJWAqHVjpytJGZc/Jp8Dlmz8ila7zaNSISufQO1wrOOR6cv5ZP91bx/DfOIzO1i9cliYgElUYOrfDMx0X8ac1uvn/lEC4Z0tPrckREgk7hcAofby3lv979lKtH9uaei7O9LkdEpF0oHE5iR/lh7pu3ikE9u/HYP43CTA1oEYkOCocTOHy0gZlz8gCYPWMcXdSAFpEoonBogXOOB/6wli37qvnVzWMYmKIGtIhEF4VDC578qJB31u3hwalDuWhwmtfliIi0O4XDF/x9834ee28z15ybzqwLs7wuR0TEEwqHJraV1XD/y6sY0qsb//1P56oBLSJRS+EQcKjO34Du0MF4ZkYuCXFqQItI9FI4AD6f43uvraZg/yGeuHks/ZMTvC5JRMRTnoSDmT1mZp+a2Voze8PMkpqse8jMCsxss5ld1R71/ObDAt7bsI8fXX0O5w9KbY+nFBEJaV6NHP4GjHDOnQtsAR4CMLNhwE3AcGAq8FsziwlmIe9v2scvFm7h+tF9uOP8zGA+lYhI2PAkHJxzf3XONQTuLgX6BW5fB7zinKtzzhUDBcD4YNVRWHqIb7+ymmHp3fnZV9SAFhH5TCj0HG4H3g3c7gvsbLKuJLCsGTObaWZ5ZpZXWlp6Rk8cH9uB0QOSeHr6ODrHBXWAIiISVoJ2So6ZLQR6t7DqYefcW4FtHgYagLmfPayF7V1L3985NxuYDZCbm9viNqfSr0cCc+6YcCYPFRGJaEELB+fc5Sdbb2a3AdcAlznnPntzLwH6N9msH7A7OBWKiMiJeHW20lTgB8C1zrnDTVYtAG4ys3gzywQGAcu9qFFEJJp59UmvJ4B44G+BJvBS59xdzrkNZvYasBH/4aZ7nXONHtUoIhK1PAkH51zOSdY9AjzSjuWIiMgXhMLZSiIiEmIUDiIi0ozCQUREmlE4iIhIM/b5RwzCl5mVAtu9ruM0pQJlXhcRIrQvjqf98Tnti+O19f4Y6Jxr8XKXEREO4cjM8pxzuV7XEQq0L46n/fE57Yvjtef+0GElERFpRuEgIiLNKBy8M9vrAkKI9sXxtD8+p31xvHbbH+o5iIhIMxo5iIhIMwoHERFpRuEQZGb2mJl9amZrzewNM0tqsu4hMysws81mdlWT5ePMbF1g3a8sgq9famZTAz9/gZn90Ot6gs3M+pvZh2a2ycw2mNm3AsuTzexvZrY18G+PJo9p8XUSKcwsxsxWmdnbgfvRvC+SzOyPgfeMTWY2ybP94ZzTVxC/gCuB2MDtR4FHA7eHAWvwT12eCRQCMYF1y4FJ+K+M9y7wJa9/jiDtm5jAz50FxAX2xzCv6wryz5wOjA3c7gZsCbwW/hv4YWD5D1vzOomUL+C7wDzg7cD9aN4Xvwe+GbgdByR5tT80cggy59xfnXMNgbtL8V/dDuA64BXnXJ1zrhgoAMabWTrQ3Tn3ifO/Al4Erm/vutvJeKDAOVfknDsKvIJ/v0Qs59we59zKwO1qYBP+66Rfh/+NgcC/1wdut/g6adeig8jM+gFfBp5tsjha90V34ELgOQDn3FHnXCUe7Q+FQ/u6Hf9IAPxvCDubrCsJLOsbuP3F5ZHoRPsgKphZBjAGWAb0cs7tAX+AAD0Dm0X6Pvpf4EHA12RZtO6LLKAU+F3gMNuzZtYFj/aHwqENmNlCM1vfwtd1TbZ5GP/V7eZ+tqiFb+VOsjwSRdPPehwz6wrMB77tnKs62aYtLIuIfWRm1wD7nXP5rX1IC8siYl8ExAJjgSedc2OAGvyHkU4kqPvDq8uERhTn3OUnW29mtwHXAJcFDhWBP+X7N9msH7A7sLxfC8sj0Yn2QUQzs474g2Guc+71wOJ9ZpbunNsTOLS4P7A8kvfRFOBaM7sa6AR0N7OXiM59Af6fr8Q5tyxw/4/4w8GT/aGRQ5CZ2VTgB8C1zrnDTVYtAG4ys3gzywQGAcsDw8ZqM5sYOEtpBvBWuxfePlYAg8ws08zigJvw75eIFfidPgdscs79osmqBcBtgdu38fnvvMXXSXvVG0zOuYecc/2ccxn4f/cfOOemEYX7AsA5txfYaWZDAosuAzbi0f7QyCH4nsB/NsHfAmekLnXO3eWc22Bmr+H/5TcA9zrnGgOPuRt4AeiMv0fxbrPvGgGccw1mdh/wHv4zl553zm3wuKxgmwJMB9aZ2erAsh8B/wW8ZmZ3ADuAGwFO8TqJVNG8L/4FmBv4Y6kI+Gf8f8S3+/7Q9BkiItKMDiuJiEgzCgcREWlG4SAiIs0oHEREpBmFg4iINKNwEBGRZhQOIm3IzDLM7EiTzzC09nFfD0y9/HaQShM5LQoHkbZX6JwbfToPcM69CnwzOOWInD6Fg0grmdl5gYs2dTKzLoGL9Yw4xWMyAhdueTYwGeNcM7vczBYHLt4SMVNOS2TR9BkireScW2FmC4D/wD+1yUvOufWteGgO/ikPZuKfT+oW4HzgWvxTZ1wflIJFzoLCQeT0/Dv+N/ha4P5WPqbYObcOwMw2AO8755yZrQMyglKlyFnSYSWR05MMdMV/ic9OrXxMXZPbvib3fegPNAlRCgeR0zMb+DH+izY96nEtIkGjv1pEWsnMZgANzrl5ZhYDLDGzS51zH3hdm0hb05TdIm0ocF3ot51zJz2L6QSPvRj4vnPumjYuS+S06bCSSNtqBBLP5ENwwG+BA8EoSuR0aeQgIiLNaOQgIiLNKBxERKQZhYOIiDSjcBARkWb+P35o6dvHpYtrAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "#d_start, slope, Hm0_shoal = offshore_depth(Hm0=9, Tp=15, depth_offshore_profile=15, depth_boundary_conditions=20)\n", "\n", "xgr, ygr, zgr = seaward_extend(xgr,[0],zgr,slope=1/20,depth=-20)\n", "\n", "\n", "\n", "plt.figure()\n", "plt.plot(xgr.T,zgr[:,:].T)\n", "plt.xlabel('x [m]')\n", "plt.ylabel('z [m]')\n" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "### Create model setup\n", "Create the model setup class" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Test som 2\n" ] } ], "source": [ "xb_setup = XBeachModelSetup('Test som 2')\n", "\n", "print(xb_setup)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Add the grid, wave boundary conditions and parameter to the model.\n", "\n", "Note that the following grid settings can be specified:\n", "\n", "`set_grid(self,xgr,ygr,zgr, posdwn=1, xori=0, yori=0,alfa=0, thetamin=-90, thetamax = 90, dtheta=10, dtheta_s=10)`\n", "\n", "Two different boundary conditions can be specified: `jonstable` or `parametric`.\n" ] }, { "cell_type": "code", "execution_count": 8, "metadata": {}, "outputs": [], "source": [ "xb_setup.set_grid(xgr,None,zgr)\n", "\n", "xb_setup.set_waves('parametric',{'Hm0':2, 'Tp':5, 'mainang':270, 'gammajsp':3.3, 's' : 10000, 'fnyq':1})\n", "#xb_setup.set_waves('jonstable',{'Hm0':[1.5, 2, 1.5],'Tp':[4, 5, 4],'gammajsp':[3.3, 3.3, 3.3], 's' : [20,20,20], 'mainang':[270,280, 290],'duration':[3600, 3600, 3600],'dtbc':[1,1,1]})\n", "\n", "xb_setup.set_params({'Wavemodel':'surfbeat',\n", " 'morphology':0,\n", " 'bedfriccoef':0.01,\n", " 'tstop':3600,\n", " 'zs0':0,\n", " 'nglobalvar':['zb','zs','H'],\n", " 'npointvar':['zs','zb'],\n", " 'nmeanvar':['zb'],\n", " 'npoints':['1 0', '6 0', '10 0', '12 0']})" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Write the model setup" ] }, { "cell_type": "code", "execution_count": 10, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Not possible to plot wave boundary\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAfYAAAIZCAYAAABUER12AAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8/fFQqAAAACXBIWXMAAAsTAAALEwEAmpwYAABEc0lEQVR4nO3dd3yV9f3//8crg4QsZkjCRmUPGWHX2lZbR60b9x6Ao7X9+PlVre2309r1aatVQFygWBG31qpV1FoZShBE9hCQQBY7YSQkef3+yBGjDchIcl3n5Hm/3XIj5zrXOXn65pgn13Xe1/uYuyMiIiKxIS7oACIiIlJ/VOwiIiIxRMUuIiISQ1TsIiIiMUTFLiIiEkNU7CIiIjFExS4iIhJDVOwiIWdmZbW+qs1sT63blx7B871jZtc1RNYjZWZJZvawma03s1IzW2BmpwWdSyQaJQQdQEQOzt3TPvvezNYB17n7m8ElahAJwAbgROBT4HRghpn1d/d1QQYTiTY6YheJUmYWZ2a3m9kaM9tiZjPMrHXkvmQzmxbZvt3M5plZlpndBZwA3Bc54r+vjuet87GR+9qb2UtmttXMVpvZ9bUe9wszezry2FIz+9jMepjZHWZWbGYbzOw7df23uPsud/+Fu69z92p3/wewFhjSEGMnEstU7CLR6wfA2dQc5bYHtgH3R+67EmgBdALaAOOBPe5+J/Af4GZ3T3P3m+t43jofG7nvSSA/8vPOB35rZifVeuz3gMeBVsAC4HVqfs90AH4FPHAo/2GRf0j0AJYcyv4i8jkVu0j0Ggfc6e757l4O/AI438wSgH3UlPJx7l7l7vPdfechPm+djzWzTsDXgNvcfa+7LwQeAi6v9dj/uPvr7l4JPA1kAr9z933AdKCrmbU82A83s0TgCWCquy8/xMwiEqFiF4leXYDnI6fLtwPLgCogi5qj5teB6Wa2ycz+ECnMQ3Ggx7YHtrp7aa1911NzNP6Zolrf7wE2u3tVrdsAaRyAmcVFfn4FUNfZBBH5Cip2kei1ATjN3VvW+kp2943uvs/df+nufYBRwBnAFZHHHfQjHQ/y2E1AazNLr7V7Z2BjffzHmJkBD1PzD5PzIkf5InKYVOwi0WsScJeZdQEws0wzOyvy/TfNrL+ZxQM7qTm9/tmRcxFwzIGe9ECPdfcNwGzg7sgEuwHAtdScNq8PE4HewPfcfc9X7SwidVOxi0Sve4CXgH+ZWSkwFxgeuS8beIaaYl4G/BuYVutx55vZNjO7t47nPdhjLwa6UnP0/jzwc3d/42j/QyL/OBkHDAQKj+Y6fZGmztwPelZOREREooiO2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhqjYRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhiQEHaA+tG3b1rt27Rp0DBERkUYzf/78ze6e+eXtMVHsXbt2JS8vL+gYIiIijcbM1te1XafiRUREYoiKXUREJIao2EVERGKIil1ERCSGqNhFRERiiIpdREQkhsTE5W4iIiINZcfuffzxX8tZvHEn3dul0TM7ff9XZloSZhZ0xC9QsYuIiBxAVbUzbloe89ZtY3Dnlry9ooSn5+fvv791ajN6ZqXTp30Gfdtn0Ld9C47NTCUhPrgT4ip2ERGRA/jbW6uY+8lW/nj+AMbkdgJgS1k5KwpLWV5YysqiUpYVljJt7nrKK6sBSEqIo1fOZ0VfU/a9stNJToxvlMwqdhERkTos3LCde2eu4txBHfaXOkCbtCRGHZfEqOPa7t9WWVXNJ5t3sXjjDpZs2smSTTt4+aNN/P39TwGIjzOeu2EUx3dq2eC5VewiIiJ1mPTOGtKTE/nV2f2+ct+E+Dh6ZKXTIyudcwfXbHN38rftYcmmmrI/JjO1gRNHsjTKTxEREYkiG7bu5l9LCxl34rGkJR1ZVZoZnVqn0Kl1Cqf2y6nnhAemy91ERES+ZOrsdZgZV4zsEnSUw6ZiFxERqaV07z6emreB0/vnkNOiedBxDpuKXUREpJYH/7OW0vJKxp5wTNBRjoiKXUREJCJ/224e+s8nfLd/Dv07tgg6zhFRsYuIiADV1c6Pn1mEAbef1ivoOEdMxS4iIgJMe389s9ds4adn9KFT65Sg4xwxFbuIiDR56zbv4u5/LufEHplcNLTTVz8gxFTsIiLSpFVWVfO/T39EQrzxu/P6h+5DXQ6XFqgREZEmbdK/15C3fht/vXBgVF7e9mU6YhcRkSZr4Ybt/PXNVZx5fHvOHtQh6Dj1QsUuIiJN0qdbdnP9Y3lkZSTz60NYDz5a6FS8iIg0Oes27+LiB+dSUVnNk9cPp0XzxKAj1RsVu4iINClrSsq4ZH+pj+C4dulBR6pXKnYREWkyVhWVcslD7+PuTB87kp7ZsVXqoGIXEZEmYkVhKZc+NBczY/rY2DtS/4wmz4mISMxbUVjKJQ/OJS7GSx10xC4iIjFueeFOLnnwfRLjjeljR9KtbWrQkRqUjthFRCRmLd64g4snz6VZfFyTKHXQEbuIiMSgrbsqeHzOeia8s5rWqc148voRdG0CpQ4hLnYzOxW4B4gHHnL33wUcSUREQq68sopJ73zCxH+vZu++ar7TJ4u7zulPZnpS0NEaTSiL3czigfuBbwP5wDwze8ndlwabTEREwqh07z5eXVzIhLdXs27Lbs4YkMMtJ3Wne1bsTpI7kFAWOzAMWO3unwCY2XTgLEDFLiIiAOzYvY+Zy4t4bXEh76wsoaKyml7Z6Tx+7TBO6J4ZdLzAhLXYOwAbat3OB4bX3sHMxgJjATp37tx4yUREJDDrt+zinRUlvLmsiDlrtlBZ7WRnJHPJsM6cMSCHIV1aRf3Hrh6tsBZ7XX8r/oUb7pOByQC5ublex/4iIhIDKquqefKDT5kyex1rSnYB0K1tKtee0I3T+uUwoEML4uKadpnXFtZizwc61brdEdgUUBYREQnIyqJS/r+nP+Kj/B0M6tySn3+vD9/o2a5JXLZ2pMJa7POA7mbWDdgIXARcEmwkERFpLJVV1Tzw7ifc8+Yq0pITuO+SQXy3f06TP81+KEJZ7O5eaWY3A69Tc7nbI+6+JOBYIiLSCN5dWcJv/7mM5YWlnN4/m1+d1Y+2aU3ncrWjFcpiB3D3fwL/DDqHiIg0rL37qlhRWMrSgp28uriQd1eW0Kl1cyZdNphT++UEHS/qhLbYRUQk9hSX7mXppp0sK6gp8mUFO/mkpIzqyBToVimJ3Hl6b64Y1YWkhPhgw0YpFbuIiNS7yqpq1m7exdKCnTVfm2pKfHNZxf59OrRsTu+cDE7vn0OfnHT65LSgY6vmmuF+lFTsIiJyVEr37mN5Yen+8l5asJMVhaWUV1YD0Cw+ju5ZaXyjZzv65GTQp30GvbMzaJGSGHDy2KRiFxGRQ+LubNy+p+Y0+qadLC3YwbKCUj7dunv/Pq1SEunTPoPLR3ShT/uaEj82M43EeH2YaGNRsYuIyH8pr6xiVVHZ/vfBPzsa37m3EgAz6Nomlf4dWnDh0E70jpxKz8pI0iVpAVOxi4g0cVt3VXyhvJcW7GR1cRmVkRltzRPj6ZmdzhnHt6dPTga9czLolZ1OapIqJIz0tyIi0kRUVlWzbssulhWUsrxw5/5T6oU79+7fJysjiT45GXyrV7ua98JzMujaJpV4TWiLGip2EZEYtLmsnOWRAl9eWPPnyqIyKiIT2uLjjOMy0xh5bJv9p9F756TTRgvBRD0Vu4hIFCuvrGJ1cdkXSnxZQSmby8r375OZnkSv7HSuGtWVnlnp9MpJ57h2abpOPEap2EVEooC7U7Bj7/5T6MsLS1lRuJM1JbuoirwX3iwhjh5ZaXyjZya9stPpnZNBz+x0LcfaxKjYRURCZld5JSuKSllRWMrygp0si/z52Yx0+Gxxl3S+0yebXjnp9MrOoGubFBJ0WVmTp2IXEQlIdbXz6dbdtY7Ca06lr9/y+XXhqc3i6ZWTwfeOb0+vnAx6Z6fTIzudjGQt7iJ1U7GLiDSC7bsraiaxFUTeBy8sZWVhKXv2VQEQZ9C1bSr92rfgvMEd959K79BSS6zK4VGxi4jUo31V1XxSsmv/UfiKyFF4wY7PLylrlZJIr+wMLhrWid7ZGfTKSad7u3SaN9NkNjl6KnYRkSPg7pSUlu+/lGx5Qc1R+OriUvZV1UxmS4w3js1MY8QxbeiVnb7/VHpmulZnk4ajYhcR+Qp799Usr7osUuCfvRe+ddfnn1SWnZFMr5x0TuyRGSnxdI5pm0azBE1mk8alYhcRiXB38rft+cJ74csLd7J28679nxeenBhHz+wMvt07a/9s9F7Z6bRKbRZseJEIFbuINEmle/exsqj089noBTWXl5WWf35JWefWKfTKTue7A9rTO3IqvXPrFC2vKqGmYheRmFZV7azbsmv/KfTPijx/2579+6QnJ9A7O4NzBneoOQLPSadHVjpp+pATiUJ61YpIzNi6q+ILC7osLyxlZVEp5bXWRz+mbSqDOrfi4mGd909oa98iWZPZJGao2EUk6pRXVrGmuOaSshWRa8KXF+ykuPTz9dHbptWsj375iC70inzM6HHt0khO1CVlEttU7CISWu5O0c7yL85GLyhlTcnnnxXeLD6O7llpnNA9k96RyWw9I5eUiTRFKnYRCYXdFZWsLCr7fGW2yJ879uzbv0+Hls3plZ3OyX3a7Z+N3q1tqtZHF6lFxS4ijaq62tmwbff+SWwrCms+qWzdll145JKy1Gbx9MxO5/T+OV84Cm/RXOuji3wVFbuINKiy8krmrd3K7DWbyVu/jRWFpeyuqFkf3Qy6tkmlV3Y6Zw/sQK+cdHpnZ9CxldZHFzlSKnYRqVfuzvLCUt5aXszby4tZuGE7ldVOs4Q4BnZsyQW5nfYfhXfPSiOlmX4NidQn/R8lIkdt774q5nyyhZnLinhrWTGbIh940r9DC8adeAyjjm3LkC6tNCNdpBGo2EXkiBTt3Mtby4uZuayYWas3s2dfFc0T4zmhe1tuObk73+zZjnYZyUHHFGlyVOwickiqq53Fm3Ywc1kxM5cXsXjjTqBmpvqY3I58q1c7RhzTRkflIgELXbGb2S+A64GSyKafuPs/g0sk0nTtrqjkvVWbmbmsmLdWFFNSWk6cweDOrfjxqT05qVcWPbLStGqbSIiErtgj/uLufwo6hEhTU1ZeyYfrt/HB2q28v3YLH23YQUVVNelJCXy9ZyYn9WrHN3q2o7U+yUwktMJa7CLSgCqrqlm/dTerikpZWVTGyqKaNdXXlOyiqtqJjzP6dWjBVaO78o0emQzt1ppELQIjEhXCWuw3m9kVQB5wq7tv+/IOZjYWGAvQuXPnRo4nEh0+WwymdnmvLCpjTUkZFZEPRjGDTq1S6JGVxil9sxnWrTWDO7ciVZ9sJhKVzD9b6qkxf6jZm0B2HXfdCcwFNgMO/BrIcfdrDvZ8ubm5npeXV+85RaKFu7Nx+579xb2yqJRVRWWsKi5l777q/ft1aNmc7llp9MxKp3tWOj2y0jiuna4lF4lGZjbf3XO/vD2Q/5vd/eRD2c/MHgT+0cBxRKKGu1O4cy8ri8oip9FLWVFUxuqiUnZFVnMDyMpIokdWOpcO70KPrDS6Z6XTvV0a6claklUk1oXun+lmluPuBZGb5wCLg8wjEgR3p6SsnFVfOoW+sqiU0r2V+/drm9aMHlnpjMnt9PmReLt0WqSowEWaqtAVO/AHMxtIzan4dcC4QNOINLCtuypqlffnBb599+efatYyJZEeWemcNbB9rdPo6ZqdLiL/JXTF7u6XB51BpCHs2L2PlcWl+9//XlFYyqriUjaXVezfJz05gR5Z6ZzWL5sekfLunpVGZlqSrhUXkUMSumIXiXale/exqrjmPfAVhTUT2FYWlVK0s3z/PqnN4jkuK51v9WoXKe90emalk5WhAheRo6NiFzlCuysqWV382ZH35zPRN27fs3+f5MQ4urdLZ/RxbekRKe/uWWm0b6GPJRWRhqFiFwHyt+3mof+sZdvuii9sd4eqaqe8soryyur9X1t3lZO/bQ+fXS3aLCGOYzPTyO3aikuyOkdOo6fRsVUK8SpwEWlEKnZp0rbvruD+t1czdfZ6MGjf4r8/jSw+zkhKiCcpMY6khDhaNk+kU6vmjBnSaf+lZF1ap5CgldlEJARU7NIk7d1XxWNz1nHfW6spLa/k/MEd+dG3e9C+ZfOgo4mIHBUVuzQ5eeu2csv0hWzcvodv9szkttN60Ss7I+hYIiL1QsUuTcoLCzby42cW0aFVc/5+3XBGHdc26EgiIvVKxS5Ngrvz1zdXcc/MVYw4pjWTLhtCyxQt7iIisUfFLjFv774qfvzMIl76aBNjhnTkrnP60yxBE91EJDap2CWmbS4rZ+xjeXz46XZuO7UX4088RgvAiEhMU7FLzFpVVMrVU+axuayciZcO5rT+OUFHEhFpcCp2iUnvrizhpic+JLlZPDPGjWRAx5ZBRxIRaRQqdok50+au5+cvLaF7uzQeuWqork0XkSZFxS4xo6raueuVZTwyay3f6tWOey8eRFqSXuIi0rTot57EhF3llfzgyQXMXF7MNaO7ced3e2uNdhFpklTsEvU2bd/DtVPzWFlUyq/P7sflI7oEHUlEJDAqdolqH+fv4Nqp89hTUcUjVw3lxB6ZQUcSEQmUil2i1muLC/jhUwtpk5rEtBuH0yMrPehIIiKBO6Tlt8xsppmd/qVtkxsmksjBuTuT/r2G8dM+pHdOBi/cNFqlLiIScahH7N2A28xsqLv/MrItt4EyiRxQRWU1P33hY2bk5XPGgBz+NOZ4khPjg44lIhIah7pg9nbgJCDLzF42sxYNF0mkbjt27+PKRz5gRl4+P/jWcdx70SCVuojIlxzqEbu5eyVwo5ldBbwHtGqwVCJfsm7zLq6ZMo/8bXv4y4XHc86gjkFHEhEJpUMt9kmffePuU8zsY+Cmhokk8kXvf7KFcdPmY8AT1w9naNfWQUcSEQmtQyp2d3/gS7fnA9c0SCKRWp6dn8/tzy2ic+sUHrlqKF3apAYdSUQk1HS5m4RSdbXz5zdWct/bqxl1bBsmXjqEFimJQccSEQk9FbuEzt59Vdw64yNe+biAi4d14ldn9SMx/lDneYqING0qdgmV4tK9XP/YfBblb+fO03tz3QndMNOa7yIih0rFLqGxvHAn107JY+uuCiZdNoRT+mYHHUlEJOqo2CUU3llRzM1/X0BqUjxPjx9Jvw5aKkFE5EgE8salmY0xsyVmVm1muV+67w4zW21mK8zslCDySeOaOnsd10yZR5c2Kbx409dU6iIiRyGoI/bFwLnAFy6jM7M+wEVAX6A98KaZ9XD3qsaPKA2tsqqa37yyjCmz13Fy7yzuuWggqUk6iSQicjQC+S3q7suAuiZFnQVMd/dyYK2ZrQaGAXMaN6E0tNK9+/j+kwt4Z0UJ15/QjdtP6018nCbJiYgcrbAdHnUA5ta6nR/Z9l/MbCwwFqBz584Nn0zqzcbte7h2yjxWFZfx23P6c8lw/f2JiNSXBit2M3sTqGta853u/uKBHlbHNq9rR3efDEwGyM3NrXMfCZ+FG7Zz3dQ8yiurmHL1UE7onhl0JBGRmNJgxe7uJx/Bw/KBTrVudwQ21U8iCdoriwr4nxkLaZeRxPSxwzmunT5DXUSkvoVtOa+XgIvMLMnMugHdgQ8CziRHyd25/+3V3PT3D+nfoQUv3DhapS4i0kACeY/dzM4B/gZkAq+Y2UJ3P8Xdl5jZDGApUAncpBnx0a2ispo7nvuYZz/M56yB7fn9eQP0GeoiIg0oqFnxzwPPH+C+u4C7GjeRNIRtuyoYN20+H6zdyo9O7sEPTjpOy8OKiDSwsM2KlxjxSUkZ10yZx6Yde7nnooGcNbDOixtERKSeqdil3s1Zs4Xx0+aTEGc8ef1whnRpHXQkEZEmQ8Uu9WpG3gZ+8tzHdGubyiNXDaVT65SgI4mINCkqdqkX1dXOH15fwaR/r+GE7m25/9LBZCQnBh1LRKTJUbHLUdtTUcWPnlrIa0sKuXR4Z35xZl8S48N2JaWISNOgYpejUrxzL9c9lsfHG3fwszP6cM3orpr5LiISIBW7HLGlm3Zy7dR57Nizjwcvz+XkPllBRxIRafJU7HJEZi4r4gdPLiA9OZGnx4+kb3t9hrqISBio2OWwuDuPzlrHb15ZSt/2LXjoylyyMpKDjiUiIhEqdjlklVXV/PLlpTw+dz2n9M3iLxcOJKWZXkIiImGi38pySHbu3cfNf1/AuytLGHfiMdx2Si/i4jRJTkQkbFTs8pU2bN3NNVPmsXbzLn5/Xn8uHNo56EgiInIAKnY5qPnrtzHu8TwqKqt57JphjDqubdCRRETkIFTsckAvfbSJ/336I3JaJPPUuKEcm5kWdCQREfkKKnb5L+7OvTNX85c3VzKsa2smXT6E1qnNgo4lIiKHQMUuX1BeWcXtz37M8ws2cu7gDtx9bn+SEuKDjiUiIodIxS77bd1VwbjH85i3bhv/+50e3PTN47Q8rIhIlFGxCwCri8u4Zso8inbu5b5LBnHGgPZBRxIRkSOgYhdmrd7M+GnzSUqIY/rYEQzq3CroSCIicoRU7E3ckx98ys9eWMyxmWk8dGUunVqnBB1JRESOgoq9iaqqdn7/2nImv/sJJ/bI5L5LBpGenBh0LBEROUoq9iZod0Ult0xfyBtLi7hiZBf+3xl9SIiPCzqWiIjUAxV7E1O4Yy/XTp3HsoKd/OJ7fbhqdLegI4mISD1SsTchizfu4Nqp8yjbW8nDVw7lm73aBR1JRETqmYq9iXhjaRE/eHIBrVISeeaGUfTOyQg6koiINAAVe4xzdx76z1p+++oyBnRowYNX5tIuPTnoWCIi0kBU7DFsX1U1/+/FJTz5waec3j+b/xszkObNtDysiEgsU7HHqB179nHjE/OZtXoLN37jWP73Oz2Ji9PysCIisU7FHoPWb9nFNVPm8enW3fzx/AGMye0UdCQREWkkgVy8bGZjzGyJmVWbWW6t7V3NbI+ZLYx8TQoiXzTLW7eVcybMZsuuCh6/drhKXUSkiQnqiH0xcC7wQB33rXH3gY0bJza8sGAjP35mER1aNeeRq4bSrW1q0JFERKSRBVLs7r4M0EeC1hN35y9vruLemasY3q01D1w+hJYpzYKOJSIiAQjjOqLdzGyBmf3bzE440E5mNtbM8swsr6SkpDHzhcrefVXcMn0h985cxZghHXn82uEqdRGRJqzBjtjN7E0gu4677nT3Fw/wsAKgs7tvMbMhwAtm1tfdd355R3efDEwGyM3N9frKHU02l5Uz9rE8Pvx0Oz8+tSc3nHiszoKIiDRxDVbs7n7yETymHCiPfD/fzNYAPYC8eo4X9VYWlXLNlHlsLitn4qWDOa1/TtCRREQkBEJ1uZuZZQJb3b3KzI4BugOfBBwrdN5dWcJNT3xIcrN4nho7kuM7tQw6koiIhERQl7udY2b5wEjgFTN7PXLX14FFZvYR8Aww3t23BpExrKbNXc/VU+bRoVVzXrhptEpdRES+IKhZ8c8Dz9ex/Vng2cZPFH5V1c5dryzjkVlr+Vavdtx78SDSkkJ1wkVEREJAzRAFysorueXJBcxcXszVo7vy0+/2IV7Lw4qISB1U7CG3afserp2ax8qiUn59Vl8uH9k16EgiIhJiKvYQW5S/neum5rG7ooqHr8zlGz3bBR1JRERCTsUeUq8tLuCHTy2kTWoSz94wnJ7Z6UFHEhGRKKBiDxl354F3P+F3ry5nUOeWTL48l8z0pKBjiYhIlFCxh0hFZTU/feFjZuTlc8aAHP405niSE+ODjiUiIlFExR4SO3bvY/y0+cz5ZAvf/9Zx/OjkHsRp5ruIiBwmFXsIrNu8i2umzCN/2x7+fMHxnDu4Y9CRREQkSqnYA/b+J1sYN20+Bky7bjjDurUOOpKIiEQxFXuAnp2fz+3PLaJT6xQevWooXdqkBh1JRESinIo9ANXVzp/fWMl9b69m1LFtmHjpEFqkJAYdS0REYoCKvZHt3VfFrTM+4pWPC7hoaCd+fXY/EuMD+SweERGJQSr2RlRcupfrH5vPovzt/OT0Xlx/wjGYaea7iIjUHxV7I1leuJNrp+SxdVcFky4bwil9s4OOJCIiMUjF3gjeXlHM9/++gNSkeJ4eP5J+HVoEHUlERGKUir2BTZ29jl++vIRe2Rk8fFUuOS2aBx1JRERimIq9gVRWVfPrfyxl6pz1nNy7HfdcNIjUJA23iIg0LDVNAyjdu4/vP7mAd1aUcN3XunHH6b2J1/KwIiLSCFTs9Wzj9j1cO2Ueq4rLuOucflw6vEvQkUREpAlRsdejhRu2c93UPMorq5hy9VBO6J4ZdCQREWliVOz15JVFBfzPjIW0y0jiyeuH0z0rPehIIiLSBKnYj5K7M+GdNfzx9RUM6dKKyZcPoU1aUtCxRESkiVKxH4Xyyip+8txinv0wn7MGtuf35w0gOTE+6FgiItKEqdiP0LZdFYybNp8P1m7lhyd355aTumt5WBERCZyK/QisKSnj2inz2LR9L/dcNJCzBnYIOpKIiAigYj9sc9ZsYfy0+STEGU+OHc6QLq2DjiQiIrKfiv0wzJi3gZ88/zFd26by6FVD6dQ6JehIIiIiX6BiPwTV1c4fXl/BpH+v4YTubbnvksG0aJ4YdCwREZH/omL/CnsqqvjRUwt5bUkhlwzvzC/P7EtifFzQsUREROoUSEOZ2R/NbLmZLTKz582sZa377jCz1Wa2wsxOCSLfZ4p37uXCyXN4fWkhP/1ub+46u59KXUREQi2olnoD6OfuA4CVwB0AZtYHuAjoC5wKTDCzQC4MX7ppJ2fdP4vVxWVMvjyX6044RpeziYhI6AVS7O7+L3evjNycC3SMfH8WMN3dy919LbAaGNbY+WYuK2LMpNm4w9PjR/LtPlmNHUFEROSIhOG88jXAq5HvOwAbat2XH9nWaD7O38H1j+XRLTOVF28eTd/2LRrzx4uIiByVBps8Z2ZvAtl13HWnu78Y2edOoBJ44rOH1bG/H+D5xwJjATp37nzUeT/Tr0MGvzqrH+cO7kBKM80tFBGR6NJgzeXuJx/sfjO7EjgDOMndPyvvfKBTrd06ApsO8PyTgckAubm5dZb/kTAzLhuhz1AXEZHoFNSs+FOB24Az3X13rbteAi4ysyQz6wZ0Bz4IIqOIiEg0Cupc831AEvBGZKb5XHcf7+5LzGwGsJSaU/Q3uXtVQBlFRESiTiDF7u7HHeS+u4C7GjGOiIhIzAjDrHgRERGpJyp2ERGRGKJiFxERiSEqdhERkRhin19CHr3MrARYH3SORtIW2Bx0iCZCY904NM6NQ+PcOBpznLu4e+aXN8ZEsTclZpbn7rlB52gKNNaNQ+PcODTOjSMM46xT8SIiIjFExS4iIhJDVOzRZ3LQAZoQjXXj0Dg3Do1z4wh8nPUeu4iISAzREbuIiEgMUbGHmJn90cyWm9kiM3vezFrWuu8OM1ttZivM7JRa24eY2ceR++61yKfsyKEzs1Mj47razG4POk80M7NOZva2mS0zsyVmdktke2sze8PMVkX+bFXrMXW+tuWrmVm8mS0ws39Ebmuc65mZtTSzZyK/m5eZ2ciwjbOKPdzeAPq5+wBgJXAHgJn1AS4C+gKnAhPMLD7ymInAWGo+8rZ75H45RJFxvB84DegDXBwZbzkylcCt7t4bGAHcFBnP24GZ7t4dmBm5/VWvbflqtwDLat3WONe/e4DX3L0XcDw14x2qcVaxh5i7/8vdKyM35wIdI9+fBUx393J3XwusBoaZWQ6Q4e5zvGbyxGPA2Y2dO8oNA1a7+yfuXgFMp2a85Qi4e4G7fxj5vpSaX4IdqBnTqZHdpvL567TO13ajho5SZtYR+C7wUK3NGud6ZGYZwNeBhwHcvcLdtxOycVaxR49rgFcj33cANtS6Lz+yrUPk+y9vl0N3oLGVo2RmXYFBwPtAlrsXQE35A+0iu2n8j9xfgR8D1bW2aZzr1zFACfBo5C2Ph8wslZCNs4o9YGb2ppktruPrrFr73EnNKc0nPttUx1P5QbbLodMYNgAzSwOeBX7o7jsPtmsd2zT+X8HMzgCK3X3+oT6kjm0a56+WAAwGJrr7IGAXkdPuBxDIOCc09A+Qg3P3kw92v5ldCZwBnOSfX5uYD3SqtVtHYFNke8c6tsuhO9DYyhEys0RqSv0Jd38usrnIzHLcvSDyFlJxZLvG/8iMBs40s9OBZCDDzKahca5v+UC+u78fuf0MNcUeqnHWEXuImdmpwG3Ame6+u9ZdLwEXmVmSmXWjZpLcB5FTQKVmNiIyG/4K4MVGDx7d5gHdzaybmTWjZuLLSwFnilqR1+HDwDJ3/3Otu14Crox8fyWfv07rfG03Vt5o5e53uHtHd+9KzWv2LXe/DI1zvXL3QmCDmfWMbDoJWErIxllH7OF2H5AEvBG5am2uu4939yVmNoOaF1QlcJO7V0UecwMwBWhOzXvyr/7Xs8oBuXulmd0MvA7EA4+4+5KAY0Wz0cDlwMdmtjCy7SfA74AZZnYt8CkwBuArXtty+DTO9e/7wBORf/h/AlxNzUFyaMZZK8+JiIjEEJ2KFxERiSEqdhERkRiiYhcREYkhKnYREZEYomIXERGJISp2ERGRGKJiFxERiSEqdhERkRiiYhcREYkhKnYREZEYomIXERGJISp2ERGRGKJiFxERiSEqdhERkRiiYhcREYkhKnYREZEYomIXERGJISp2ERGRGKJiFxERiSEqdhERkRiiYhcREYkhKnYREZEYomIXERGJIQlB/nAzewQ4Ayh2936RbX8EvgdUAGuAq919+8Gep23btt61a9eGDSsiIhIi8+fP3+zumV/ebu4eRJ6aH272daAMeKxWsX8HeMvdK83s9wDuftvBnic3N9fz8vIaPK+IiEhYmNl8d8/98vZAT8W7+7vA1i9t+5e7V0ZuzgU6NnowERGRKBX299ivAV6t6w4zG2tmeWaWV1JS0sixREREwim0xW5mdwKVwBN13e/uk909191zMzP/6y0GERGRJinQyXMHYmZXUjOp7iQPchKAiIhIlAldsZvZqcBtwInuvjvoPCIiItEk0FPxZvYkMAfoaWb5ZnYtcB+QDrxhZgvNbFKQGUVERKJJoEfs7n5xHZsfbvQgtWwuK+euV5YFGUECZgbpSQmkJyeSnlzzZ4vmieR2bUVWRnLQ8UREDip0p+KDVl5Zzfz124KOIQGqqnbKyisp3buP6lozPMxgaJfWnN4/m9P656jkRSSUAl2gpr5ogRppCO7O7ooqSvdWsrmsnLeWF/PKogJWFJXWlHzX1ny3fw6n9cumnUpeRBrZgRaoUbGLHKbVxaW8sqiQVz7exMqisv0lf8aAHE7tl027dJW8iDQ8FbtIA1hVVMorHxfwyqICVhXXlPywSMmfopIXkQakYhdpYCuLSnllUQGvfFzA6kjJD+9Wc7r+1H45ZKYnBR1RRGKIil2kkbg7K4vKIkfym1hTsos4g1HHtuXMge05pW82LZonBh1TRKKcil0kAJ+V/MsfbeKljzbx6dbdNIuP45u9Mjnz+A6c1LsdyYnxQccUkSikYhcJmLvzUf4OXly4kX8sKqCktJzUZvGc0jebMwe2Z/RxbUmMD+3HN4hIyKjYRUKkqtp5/5MtvLhwE68uLmDn3kpapzbju/1zuP6EY+jcJiXoiCIScip2kZAqr6zi3ZWbeXHhRt5cVkS1w/UndOPGbxxHapLWkBKRuqnYRaJA0c69/P7V5Ty3YCNZGUncflovzh7YATMLOpqIhMyBil1v6ImESFZGMn++cCDP3TiK7IxkfvTUR5w3cTari0uDjiYiUULFLhJCgzu34vkbR/OnMcezfstuzr5/Nm8uLQo6lohEARW7SEjFxRnnD+nIy9//Gl3bpnD943nc//ZqYuHtMxFpOCp2kZBr37I5z4wfxZnHt+ePr6/g5r8vYHdFZdCxRCSkVOwiUSA5MZ6/XjiQn5zei1cXF3DuhNls2Lo76FgiEkIqdpEoYWaM/fqxPHLVUDZu38OZ973H7DWbg44lIiGjYheJMt/o2Y6Xbv4abdKSuPzhD5gya63edxeR/VTsIlGoW9tUnr9xFN/smckvXl7Kbc8uoryyKuhYIhICKnaRKJWenMjky3P5wbeOY0ZePhdPnkvxzr1BxxKRgKnYRaJYXJzxP9/pyYRLB7OsoJTv3fceCzdsDzqWiARIxS4SA07vn8NzN44iMT6OCx6Yw7Pz84OOJCIBUbGLxIjeORm8dPPXGNK5Fbc+/RG//sdSKquqg44lIo1MxS4SQ1qnNuOxa4dx1aiuPPzeWq56dB7bd1cEHUtEGpGKXSTGJMbH8Ysz+/KH8wfwwdqtnHnfLFYU6kNkRJoKFbtIjLogtxPTx41gz74qzpkwi9cWFwYdSUQagYpdJIYN7tyKl2/+Gt2z0hk/bT5/eWMl1dVazEYklqnYRWJcdotknho7gvMGd+SemasYP20+ZeX6EBmRWBVosZvZI2ZWbGaLa21rbWZvmNmqyJ+tgswoEguSE+P505gB/OyMPsxcXsy5E2axfsuuoGOJSAMI+oh9CnDql7bdDsx09+7AzMhtETlKZsa1X+vG1KuHUbSznDPvm8V/VpUEHUtE6lmgxe7u7wJbv7T5LGBq5PupwNmNmUkk1n2te1teunk02RnJXPnIBzz0n0+CjiQi9SjoI/a6ZLl7AUDkz3Z17WRmY80sz8zySkp01CFyOLq0SeW5G0fx7T5Z/OaVZUydvS7oSCJST8JY7IfE3Se7e66752ZmZgYdRyTqpCYlMOHSIZzcO4tfvryEmcuKgo4kIvUgjMVeZGY5AJE/iwPOIxKz4uOMey8eSJ/2GXz/yQUs3rgj6EgicpTCWOwvAVdGvr8SeDHALCIxL6VZAg9fOZSWzRO5duo8CnbsCTqSiByFoC93exKYA/Q0s3wzuxb4HfBtM1sFfDtyW0QaUFZGMo9cPZRd5VVc/eg8XecuEsWCnhV/sbvnuHuiu3d094fdfYu7n+Tu3SN/fnnWvIg0gF7ZGdx/6WBWFZdx27OLcNcKdSLRKIyn4kUkICf2yORHJ3fnlUUFvLBwY9BxROQIqNhF5Atu+MZx5HZpxf97YQn523YHHUdEDpOKXUS+ID7O+MuFA6l259YZH1GlD40RiSoqdhH5L51ap/DzM/vy/tqtWplOJMqo2EWkTmOGdOSUvln86V8rWLppZ9BxROQQqdhFpE5mxt3nDqBlSjN++NQC9u6rCjqSiBwCFbuIHFDr1Gb84fwBrCwq4w+vrQg6jogcAhW7iBzUN3u24/IRXXhk1lreW7U56Dgi8hVU7CLylX5yem+OyUzlf5/+iO27K4KOIyIHoWIXka/UvFk891w4iM1l5fz0hcValU4kxFTsInJI+ndswQ9P7s4/FhXw4sJNQccRkQNQsYvIIRt/4rEM6dKKn72wWKvSiYSUil1EDllCfBx/uUCr0omEmYpdRA5L5zYp/Px7WpVOJKxU7CJy2MbkalU6kbBSsYvIYdOqdCLhpWIXkSNSe1W6P76uVelEwkLFLiJH7LNV6R5+T6vSiYSFil1EjopWpRMJFxW7iByV5s3i+euFA7UqnUhIqNhF5KgN6NhSq9KJhISKXUTqxf5V6V5czMbte4KOI9JkqdhFpF7sX5Wu2vnh9AVUVFYHHUmkSVKxi0i96dwmhd+e259567bxm1eWBh1HpElKCDqAiMSWswZ2YPHGHTz4n7X0bZ/BhUM7Bx1JpEnREbuI1LvbTu3F145ry89eWMKHn24LOo5Ik6JiF5F6lxAfx98uHkRWiyTGPz6fop17g44k0mSo2EWkQbRKbcbky3MpK69kzKQ5rC4uCzqSSJMQ2mI3sx+Z2RIzW2xmT5pZctCZROTw9M7JYNp1w9ldUcm5E2Yxe42WnRVpaKEsdjPrAPwAyHX3fkA8cFGwqUTkSAzu3IrnbxxNu4xkrnj4A2bkbQg6kkhMC2WxRyQAzc0sAUgBtJyVSJTq1DqFZ28YxfBjWvPjZxbxh9eWU12tpWdFGkIoi93dNwJ/Aj4FCoAd7v6v2vuY2VgzyzOzvJKSkiBiishhaNE8kSlXD+OioZ2Y8M4avj9dn+Mu0hBCWexm1go4C+gGtAdSzeyy2vu4+2R3z3X33MzMzCBiishhSoyP4+5z+3PHab14ZVEBFz84l81l5UHHEokpoSx24GRgrbuXuPs+4DlgVMCZRKQemBnjTjyWiZcOZummnZx9/yxWFZUGHUskZoS12D8FRphZipkZcBKwLOBMIlKPTuufw1PjRrJ3XzXnTpzNe6s0Y16kPoSy2N39feAZ4EPgY2pyTg40lIjUu4GdWvLCTaNo36I5Vz36AdM/+DToSCJRL5TFDuDuP3f3Xu7ez90vd3e9EScSgzq2SuHpG0Yy8tg23P7cx/zuVc2YFzkaoS12EWk6MpITefSqoVw6vDOT/r2Gm/7+IXsqNGNe5Eio2EUkFBLi4/jN2f346Xd789qSQi56cC7FpVpjXuRwqdhFJDTMjOtOOIZJlw1hZWEp59w/mxWFmjEvcjhU7CISOqf0zWbGuJHsq6rm/ImzeXelFqESOVQqdhEJpf4dW/DCTaPp0Ko5V0+ZxxPvrw86kkhUULGLSGi1b9mcZ24YxQnd23Ln84u565WlmjEv8hVU7CISamlJCTx0RS5XjOzCg/9Zy/hp89ldURl0LJHQUrGLSOglxMfxq7P68fPv9eGNZUVc+MBcindqxrxIXVTsIhI1rh7djQcvz2VNSRln3z+LZQU7g44kEjoqdhGJKif3yWLGuJFUuTNm0hzeXlEcdCSRUFGxi0jU6dehBS/e9DU6t07h2inzeHzOuqAjiYSGil1EolJ2i2SeHj+Sb/Zsx89eXMKvXl5KlWbMi6jYRSR6pSYlMPmKXK4e3ZVHZq1l3OPz2VWuGfPStNVLsZvZr80sodbtDDN7tD6eW0TkYOLjjJ9/ry+/Oqsvby0v4oIH5lC4QzPmpemqryP2BOB9MxtgZt8B5gHz6+m5RUS+0hUju/LwlUNZt3kXZ98/iyWbdgQdSSQQ9VLs7n4HcBvwPjAF+K6731cfzy0icqi+2asdT48fhRmMmTSHt5YXBR1JpNHV16n4rwP3Ar8C/g3cZ2bt6+O5RUQOR5/2Gbxw02iOyUzluql5TJm1NuhIIo2qvk7F/wk4z93vdveLgcnAW/X03CIihyUrI5kZ40ZyUu8sfvHyUn7x0hLNmJcmI+GrdzkwM/ufyLdPAaeZ2Wm17n7waJ5bRORopDRLYNJlQ7j7n8t46L21fLp1N/dePIi0pKP6tScSekd7xJ4e+RoC3AB0iHyNB3of5XOLiByV+Djjp2f04ddn9+PfK0sYM2kOBTv2BB1LpEEdVbG7+y/d/ZdAW2Cwu9/q7rdSU/Qd6yOgiMjRunxEFx65aigbtu7mrPtmsXijZsxL7Kqv99g7AxW1blcAXevpuUVEjtqJPTJ55oaRJMbHMWbSHN5YqhnzEpvqq9gfBz4ws1+Y2c+puextaj09t4hIveiVncHzN42iR1YaYx/P46H/fIK7JtVJbKmv69jvAq4GtgHbgavd/e76eG4RkfrULj2Z6WNHckqfbH7zyjJ+9uJiKquqg44lUm/qbXqou38IfFhfzyci0lCaN4tnwqWD+f1ry3ng3U/YsHUP910yiPTkxKCjiRw1fQiMiDRJcXHGHaf35u5z+/Pe6s2MmTSHjds1Y16in4pdRJq0i4d1ZurVw9i4bQ9n3z+LRfnbg44kclRU7CLS5H2te1uevXEUzeLjuOCBOby2uDDoSCJHLLTFbmYtzewZM1tuZsvMbGTQmUQkdvXISueFm0bTKzuDG56Yz+R312jGvESl0BY7cA/wmrv3Ao4HlgWcR0RiXGZ6EtPHjuD0fjn89p/LufOFxezTjHmJMqFcNNnMMoCvA1cBuHsFX1wAR0SkQSQnxvO3iwfRpU0KE95Zw4atu7n/0sFkaMa8RImwHrEfA5QAj5rZAjN7yMxSa+9gZmPNLM/M8kpKSoJJKSIxKS7O+PGpvfjDeQOYs2YL50+czYatu4OOJXJIwlrsCcBgYKK7DwJ2AbfX3sHdJ7t7rrvnZmZmBpFRRGLcBUM78dg1wyjcsZdzJsxi4YbtQUcS+UphLfZ8IN/d34/cfoaaohcRaVSjjmvLczeOpnmzeC58YA7//Lgg6EgiBxXKYnf3QmCDmfWMbDoJWBpgJBFpwo5rl8YLN46mX4cW3PjEh0x8RzPmJbxCWewR3weeMLNFwEDgt8HGEZGmrE1aEk9cN5zvHd+e37+2nNuf/Vgz5iWUQjkrHsDdFwK5QecQEflMcmI891w4kK5tUvjbW6vJ376bCZcOoUVzzZiX8AjzEbuISOjExRm3fqcnfxpzPB+s3cp5mjEvIaNiFxE5AucP6chj1wynpLScs++fxfz124KOJAKo2EVEjtjIY9vw3I2jSEtO4OIH5/LyR5uCjiSiYhcRORrHZqbx/I2jOb5jC77/5ALuf3u1ZsxLoFTsIiJHqXVqM6ZdN5yzB7bnj6+v4P97ZhEVlZoxL8EI7ax4EZFokpQQz18uHEjXtqn89c1VbNy2h0mXDaFFimbMS+PSEbuISD0xM354cg/+cuHxzF+/jXMmzmL9ll1Bx5ImRsUuIlLPzhnUkWnXDWfbrgrOvn8Weeu2Bh1JmhAVu4hIAxjWrTXP3zialinNuOTB93lx4cagI0kToWIXEWkgXdum8twNoxjYuSW3TF/IvTNXaca8NDgVu4hIA2qV2ozHrx3GuYM78Oc3VnLrjI8or6wKOpbEMM2KFxFpYEkJ8fzfmOPp1iaV/3tjJfnb9/DAZUNoldos6GgSg3TELiLSCMyM75/UnXsuGsjCDds5d+Js1m7WjHmpfyp2EZFGdNbADvz9uuHs2LOPcybM4v1PtgQdSWKMil1EpJHldm3N8zeOonVqMy57+H2eX5AfdCSJISp2EZEAdGmTyvM3jCa3S2t+9NRH/OWNlZoxL/VCxS4iEpAWKYlMvWYYY4Z05J6Zq/jhUwvZu08z5uXoaFa8iEiAmiXE8YfzB9C1bSp/fH0FG7ftYfIVubTWjHk5QjpiFxEJmJlx0zeP475LBrFo4w7OmTCLNSVlQceSKKViFxEJiTMGtGf62BGU7a3k3AmzmbNGM+bl8KnYRURCZHDnVrxw02gy05O44pH3eWa+ZszL4VGxi4iETKfWKTx7wyiGdWvN/z79EX96fQXV1ZoxL4dGxS4iEkItmicy5ephXDS0E/e9vZofTF+gGfNySDQrXkQkpBLj47j73P50a5vK3a8uZ9P2PTx4RS5t0pKCjiYhpiN2EZEQMzPGnXgsEy8dzJJNOzl7wixWF5cGHUtCTMUuIhIFTuufw1PjRrKnoppzJsxm9urNQUeSkFKxi4hEiYGdWvLCTaPIaZHMFY98wIx5G4KOJCGkYhcRiSIdW6XwzA2jGHlsG3787CJ+/9pyzZiXL1Cxi4hEmYzkRB69aiiXDO/MxHfWcPOTH2rGvOwX2mI3s3gzW2Bm/wg6i4hI2CTEx3HX2f248/TevLq4kCse/oAdu/cFHUtCILTFDtwCLAs6hIhIWJkZ13/9GP528SAWbtjO+ZNms3H7nqBjScBCWexm1hH4LvBQ0FlERMLujAHtmXrNMAp37uXcCbNYVrAz6EgSoFAWO/BX4MdA9YF2MLOxZpZnZnklJSWNFkxEJIxGHtuGZ8aPwjAumDSH2Wt0OVxTFbpiN7MzgGJ3n3+w/dx9srvnuntuZmZmI6UTEQmvntnpPHfjKHJaJnPVI/N4+aNNQUeSAISu2IHRwJlmtg6YDnzLzKYFG0lEJDq0b9mcp8eNYmDnlnz/yQU89J9Pgo4kjSx0xe7ud7h7R3fvClwEvOXulwUcS0QkarRISeSxa4Zxev9sfvPKMn7zj6W61r0J0YfAiIjEoOTEeP528WDapS/loffWUrhzL/93wfEkJcQHHU0aWKiL3d3fAd4JOIaISFSKjzN+/r0+5LRI5u5Xl7O5rJzJV+SSkZwYdDRpQKE7FS8iIvXns0+H++uFA5m/fhsXTJpD4Y69QceSBqRiFxFpAs4e1IFHrxpG/rY9nDthFquK9NGvsUrFLiLSRHyte1ueGjeCfdXOeRNn88HarUFHkgagYhcRaUL6tm/BczeMom16Epc9/D6vflwQdCSpZyp2EZEmplPrFJ4dP4p+7TO48e8fMnX2uqAjST1SsYuINEGtUpvxxHUjOLl3Fj9/aQm/f2057rrWPRao2EVEmqjmzeKZdNkQLo18rvutMz6iovKAH9EhUSLU17GLiEjDio8zfnN2P3JaJPOnf62kpKyciZcNIS1J9RCtdMQuItLEmRk3f6s7fzx/ALPXbOHCB+ZQXKpr3aOVil1ERAAYk9uJh67MZe3mXZw7YTZrSsqCjiRHQMUuIiL7fbNnO568fgR7Kqo4f+JsPvx0W9CR5DCp2EVE5AuO79SS524cRUbzRC55cC5vLC0KOpIcBhW7iIj8ly5tUnn2hlH0zEpn3ON5/P39T4OOJIdIxS4iInVqm5bEk2NHcGKPTH7y/Mf8+V8rdK17FFCxi4jIAaU0S+DBK3K5ILcj9761mtueXcS+Kl3rHma6UFFERA4qIT6O3583gOwWzbl35ipKSsu5/9LBpDRThYSRjthFROQrmRn/8+0e/Pac/vx7ZQkXT57L5rLyoGNJHVTsIiJyyC4Z3pkHLs9lRVEp50+czfotu4KOJF+iYhcRkcPy7T5ZPHHdCHbs2ce5E2azKH970JGkFhW7iIgctiFdWvHMDaNo3iyeiybP5e0VxUFHkggVu4iIHJFjM9N47sZRdGubynVT85iRtyHoSIKKXUREjkK79GSeGjeSUce24cfPLOJvM1fpWveAqdhFROSopCUl8PCVQzlnUAf+742V/PSFxVRVq9yDoosQRUTkqDVLiOPPFxxPdotkJr6zhuLScu69aBDNm8UHHa3J0RG7iIjUCzPjtlN78csz+/LmsiIufWgu23ZVBB2ryVGxi4hIvbpyVFcmXjqYxZt2ct6k2WzYujvoSE2Kil1EROrdqf1yeOK64WwuLefcibNZvHFH0JGaDBW7iIg0iKFdW/PsDaNIjDMumjyX91ZtDjpSkxDKYjezTmb2tpktM7MlZnZL0JlEROTwdc9K57kbR9OxVXOuevQDXliwMehIMS+UxQ5UAre6e29gBHCTmfUJOJOIiByB7BbJzBg/ktyurfjhUwu5/+3V+ujXBmTRsJCAmb0I3Ofub9R1f25urufl5TVyKhERORzllVXcOuMj/rGogOaJ8STGW9CRGtXT40fRMzu93p7PzOa7e+6Xt4f+OnYz6woMAt7/0vaxwFiAzp07N34wERE5LEkJ8fzt4kGcM6gD763eTBQcV9arFs0TG+XnhPqI3czSgH8Dd7n7cwfaT0fsIiLS1BzoiD2s77FjZonAs8ATByt1ERER+Vwoi93MDHgYWObufw46j4iISLQIZbEDo4HLgW+Z2cLI1+lBhxIREQm7UE6ec/f3gKY1XVJERKQehPWIXURERI6Ail1ERCSGqNhFRERiSKivYz9UZlYCrA86RwNpC+iTExqHxrrxaKwbj8a6cTXmeHdx98wvb4yJYo9lZpZX1wIEUv801o1HY914NNaNKwzjrVPxIiIiMUTFLiIiEkNU7OE3OegATYjGuvForBuPxrpxBT7eeo9dREQkhuiIXUREJIao2EVERGKIij0kzOyPZrbczBaZ2fNm1rLWfXeY2WozW2Fmp9TaPsTMPo7cd2/kU/HkMJnZqZGxXW1mtwedJ9qZWScze9vMlpnZEjO7JbK9tZm9YWarIn+2qvWYOl/jcmjMLN7MFpjZPyK3NdYNxMxamtkzkd/Xy8xsZNjGW8UeHm8A/dx9ALASuAPAzPoAFwF9gVOBCWYWH3nMRGAs0D3ydWpjh452kbG8HzgN6ANcHBlzOXKVwK3u3hsYAdwUGdPbgZnu3h2YGbn9Va9xOTS3AMtq3dZYN5x7gNfcvRdwPDXjHqrxVrGHhLv/y90rIzfnAh0j358FTHf3cndfC6wGhplZDpDh7nO8ZgbkY8DZjZ07BgwDVrv7J+5eAUynZszlCLl7gbt/GPm+lJpffB2oGdepkd2m8vnrtc7XeKOGjmJm1hH4LvBQrc0a6wZgZhnA14GHAdy9wt23E7LxVrGH0zXAq5HvOwAbat2XH9nWIfL9l7fL4TnQ+Eo9MLOuwCDgfSDL3QugpvyBdpHd9HdwdP4K/BiorrVNY90wjgFKgEcjb308ZGaphGy8VeyNyMzeNLPFdXydVWufO6k5lfnEZ5vqeCo/yHY5PBrHBmJmacCzwA/dfefBdq1jm/4ODoGZnQEUu/v8Q31IHds01ocuARgMTHT3QcAuIqfdDyCQ8U5o6B8gn3P3kw92v5ldCZwBnOSfLzCQD3SqtVtHYFNke8c6tsvhOdD4ylEws0RqSv0Jd38usrnIzHLcvSDyVlJxZLv+Do7caOBMMzsdSAYyzGwaGuuGkg/ku/v7kdvPUFPsoRpvHbGHhJmdCtwGnOnuu2vd9RJwkZklmVk3aibJfRA53VNqZiMis+GvAF5s9ODRbx7Q3cy6mVkzaia6vBRwpqgWeT0+DCxz9z/Xuusl4MrI91fy+eu1ztd4Y+WNZu5+h7t3dPeu1Lx233L3y9BYNwh3LwQ2mFnPyKaTgKWEbLx1xB4e9wFJwBuRq9bmuvt4d19iZjOoefFUAje5e1XkMTcAU4Dm1Lwn/+p/PasclLtXmtnNwOtAPPCIuy8JOFa0Gw1cDnxsZgsj234C/A6YYWbXAp8CYwC+4jUuR0Zj3XC+DzwRORD4BLiamoPk0Iy3lpQVERGJIToVLyIiEkNU7CIiIjFExS4iIhJDVOwiIiIxRMUuIiISQ1TsIiIiMUTFLiIiEkNU7CJy2MxsqJktMrNkM0uNfO56v6BziYgWqBGRI2Rmv6FmffLm1KyffXfAkUQEFbuIHKHIkprzgL3AKC1NKhIOOhUvIkeqNZAGpFNz5C4iIaAjdhE5Imb2EjAd6AbkuPvNAUcSEfTpbiJyBMzsCqDS3f9uZvHAbDP7lru/FXQ2kaZOR+wiIiIxRO+xi4iIxBAVu4iISAxRsYuIiMQQFbuIiEgMUbGLiIjEEBW7iIhIDFGxi4iIxJD/H8LN1DDwQEnxAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sim_path = os.path.join('xb-1D')\n", "if not os.path.exists(sim_path):\n", " os.mkdir(sim_path)\n", "xb_setup.write_model(sim_path)" ] }, { "cell_type": "code", "execution_count": null, "id": "72bcee2a", "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.8.12" } }, "nbformat": 4, "nbformat_minor": 5 }