Created
June 18, 2022 02:48
-
-
Save julia-neme/8e1e2528114753d1027f097bc1e7f179 to your computer and use it in GitHub Desktop.
Checking mask, kmt and topog
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| { | |
| "cells": [ | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "e328f7b3-60c9-4f44-96f9-4e23b2a042ba", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "import xarray as xr\n", | |
| "import matplotlib.pyplot as plt\n", | |
| "import numpy as np" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "88405ba6-e432-44e7-a23c-8daed3a804a5", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "mask_orig = xr.open_dataset('/g/data/ik11/inputs/access-om2/input_20200530/mom_025deg/ocean_mask.nc')['mask']\n", | |
| "kmt_orig = xr.open_dataset('/g/data/ik11/inputs/access-om2/input_20200530/cice_025deg/kmt.nc')['kmt']\n", | |
| "topog = xr.open_dataset('/scratch/e14/jn8053/access-om2/input/025deg_jra55_ryf_nogap/mom_025deg/topog.nc')['depth']" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 9, | |
| "id": "5767ce15-506e-4a66-8c7f-8f36abbbf635", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "<matplotlib.collections.QuadMesh at 0x146098dda790>" | |
| ] | |
| }, | |
| "execution_count": 9, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAwvklEQVR4nO2deRQlVXXufx+ziAgIkmZIQFdLHhgRbRmMy6hEAaO2MdGHcUBECQYV9SUBnms55ZFgjESNCvbDAacgD4l2ErBR4vCMyqBhRqEFhRYCAoI8iUB3f++PqktX375D3Xurbp2qu39r1bq35l2nqs5X55x99pFtgiAIgmBWNmvagCAIgqAbhKAEQRAElRCCEgRBEFRCCEoQBEFQCSEoQRAEQSWEoARBEASVEIISBEGQEJI+IekOSVcPWS9JH5K0WtKVkp4ybxuHEYISBEGQFp8CDh+x/ghgaT4dC5w+B5tKEYISBEGQELa/Bdw9YpPlwKed8T1gB0lL5mPdaLZo2oC62EpbexseOdW+T3jS/RVbExS5/sptN0nj66/cdqP54vr+dSkwyr5B6wZdbzxn1THrM3Ifv7jT9i7T7n/Ysx/pu+5eV2rb71/5wDXArwuLVtheMcHpdgduKcyvyZfdNsExaqGzgrINj+QgHTrdzldlP6tuvaI6gwIO221/AO65bXC69taPSvfeNk2zatUGG/ttKq4rbhPPUz0cttv+HKTZjvE1n/vTWfa/8+51XLxqj1Lbbrnkx7+2vWyG0w262iRiaHVWUKogMoFqKZOWqYvJqluv2FRA+mzuf24GzQ86bjAdg+7J/DHrvH5eJ1sD7FmY3wO4dV4nH0W0oQTJ0MtUh2UOTWe6vYxrkIAM+l/cb9i6/vXNZ4zBNBhYj0tNFbASeHXu7XUwcK/txqu7IEooQUKUqfJqmlFi0ls/SBR6y4aVbgaJUsrpEGzKeqopoUj6R+BZwM6S1gDvBLYEsH0GcD7wfGA1cD9wdCUnroAQlJLMUv0VVRwbGJVZ9jLjYVVEKaRZ8V6uuvWKsQJTxblSuO5ZqbP6OIVSnTEPVVTlZfvlY9YbOL6Sk1VMVHmNoVgNMc2DO+iLtAsZxLSMqv4Zl75NZxzDzt97RkbZN+yel7mmpq87dVJ4nwysw6WmLhOCMiHTvtyLLiRFhqVFb1n/ujalW+/ayohMj66LyjSlkzZe7xzbUJIlBGUKJn3Y25QhNsWgqq7UmMW2KjLINmayMH26lb3eFNLFwDq71NRlQlCmZBKPnBQe+NQY5Rk1rM0pBbHp2THOlnFuwtOeO9hAal5x60tOXSYEpUa61KhaJdO0SaWWeRSZRTzKCFMqYpoaKaWJS7afdL0NJby8ZmRUNU0aHa7Soz9d+vt3jPIAa5Ky3nq9aynaPKqjY5mMcdTHSXgRNo8ND3VbK0oRglIBozqsxYtdjtR7yI9jEsEb90yMOs4gIS6zXVdJ4UMjQ6wbGBFlsahNUCR9AngBcIftJ+bLdgK+AOwF/AR4me1f5OtOBo4B1gFvtr0qX/5UsnDOjyDr0HNC7oedNF1+iauizWlU1vZJ3IUnzRwnacNrc1q3AQPrk8+V6qfONpRPsWlM/5OAi2wvBS7K55G0L3AksF++z0clbZ7vczpZzP9e/P9R4wQELWFURjtr35+6KboF9y8v/k5KitcalGddXkoZN3WZ2gRlSEz/5cBZ+f+zgBcXlp9t+wHbN5GFFDgwj/G/ve3v5qWSTxf2CTpOGzLYfmEp2x4SJYZukXVsDEGZdxvKrr0gZrZvk/TYfPnuwPcK2/Xi+z+U/+9fPhBJx5KVZtiGZsfQiExjMtogHj3a2KdkXITkWY5Xx3M+LAxPqhh4yOE0m0oKDIvvP1Hcf9srbC+zvWxLtq7MuKAZUsxIZgmhMsl2dVKlmAQZRqxjs1JTl5l3CeV2SUvy0skS4I58+bD4/mvy//3Lg46RjrfOaGa1cVwo+7qpord/ikKfAuvd7eqsMsxbLlcCR+X/jwK+XFh+pKStJe1N1vh+SV49dp+kgyUJeHVhnyBoLU24lFcVOmaSOGVV0IYPjWhDyajTbXhQTP9TgXMkHQPcDLwUwPY1ks4BrgXWAsfb7g3Q/AY2uA1fkE9JEl9u1dP0F31ZUr/3dbZz1HXtZUqt6TwXYl20odTq5fVy20tsb2l7D9sft32X7UNtL81/7y5sf4rtx9vex/YFheWX2X5ivu6NKfVBKYbFaEvGlzJ1ji1SN7NEoW4Do1yi23Sf6iIbsXGzUlOXiZ7yMxBjnQSpU+XzOKjEEM97hi0e9ObjN+w43ZbLORNfatWyCJlV3c9M1cdvqt2nDe/WelRq6jJRQqmY8ISZjTZlIItKU892yn1Sskb5+D6PFKiJaTPElMO0z5NUM4620bV0TPfdyBrly0xdpttX1yBVdB4LcWkXca/qYdbOpPMgGuUzosqrJopVX5NWgxWrfabZv2uklHEEzdCGZ39ddGwMQambQUPdjno5Bg281C8wi0SISdAGjHjIkZ12u/zVMRZNTNpIHeOZBOnTa5Rf9Fhe3b66lrLobScpj4XSJuIDZH4Ysc7lpi4TZbQ5U+Ylj4ygvWIS925x6XqDexkiBRpiUIZZ9Oxqa4Y6K22/7knsr1N8ui5sqT0nNuE2TAhK44x6MRZJXBblOvvpesZfBW14LrJG+c1LTV0mqrzmzCCvrza8MMOY1fOszdc+iBQEIgUbilTp9p6yp2PXG9zLEILSElJ8kaYVxNSuownqFNJUnpX+a5zWrjZEoTaKAbaIKq9WMWiQo6qON+992zRe+CQ0Ka6pRLuuq/oyRSEpEm7DUUJpHVV99cH0mVmVL3bqmURddLFBvsyzOOh+jxpjpS3BQg2s73iDexkiBRaAEID50lQapeDYMM35m7a5GsoN/xtDAAfJUybmV9kvw2C+VJmZplwymfQYxf9tEBxD5z24yhAllI4wqH2lyvG4U/j6bRNl0qpL6VnFtQzygGwLtljvzUpNXabbVxdUJipRmqmWWTPMYuDQphvj6zj3oOGGU6fKjo2SDpf0I0mrJZ00YP2jJf2zpCskXSPp6MovaApCUILSpY8QlcmoK0Psj0LdNHW7QLeBbDyUaoYAlrQ58BHgCGBf4OWS9u3b7HjgWtv7A88C3i9pq0ovagpCUIKHCVGZnUnHvFlUmi5VVU+lIzYeCKy2faPtB4GzgeV92xh4lCQB2wF3A2urvKJpCEEJNmGcsHQrI6ieMunTli/vccwy1DVsLCxFF+G2pU/mNqxSE7CzpMsK07F9h9sduKUwvyZfVuTDwH8DbgWuAk6wvb6eqytPCEqwEaN6v7ftJZ83k/SZmHffoRQpplNb+psMY8JYXnfaXlaYVvQdblC9mPvmDwMuB3YDngx8WNL2FV/WxITb8IJTpvGzzd4386Tn+tpL01GZf1096kd1CKxajMqWZBflmakwfP0aYM/C/B5kJZEiRwOn2jawWtJNwG8Dl1RlxDRECSUYWJ/d3yegS1/GKTBJW8skaT+q2mieGfuiPS9Z+PrKBti6FFgqae+8of1IYGXfNjcDhwJI2hXYB7ixwkuaihCUBWSYq2nZL+quZxaD0qVsxl5lZN2qacKmReqPM0EbykhsrwXeCKwCrgPOsX2NpOMkHZdv9lfA0yVdBVwEnGj7zpourTRR5bWAjKqOmXR51xh3nVVGvq2rr8Wg4867umuS47Sxz0k/WbTh6r7PbZ8PnN+37IzC/1uB51V2woqIEsqC0vYXuGqGldjGZXZlw9zMm2FhTGaljAfWpG0nKaTXrGShVzYrNXWZbl9dMJJZXT67xKgv+rraO8ocb9b9512yXLSG+A1E6BUIQekcdWYgbewfUIZRGe+k8c6mOXdX6NK1TENVPeXbTCNtKJLeCryOrKR4FZkL3LbAF4C9gJ8AL7P9i3z7k4FjgHXAm22vmr/V6VOmLnrRX/p+yraZDKvrn9S5oa108UOiSnpeXovO3EsoknYH3gwss/1EYHMyt7iTgItsLyXzWjgp337ffP1+wOHAR/NYNwvPtC6lweRUUYKZ5Nhxr9pHalVekp4yaqrjnE15eW0BPELSQ2Qlk1uBk8mCnAGcBXwDOJEshs3Zth8AbpK0mizWzXfnbHNSTPpVPM8Mqn9slq5/3VYxHsjitj10g0THlH//iHUGnlP1CecuKLZ/JunvyDrm/Bdwoe0LJe1q+7Z8m9skPTbfZXfge4VDDIprs1AME5M6+0CUde8c5bIamWW3KQrrsGgBXX0GDKxNrMHd9rPnfc65C4qkHclKHXsD9wD/R9IrR+0yYFl/XJvesY8FjgXYhm1nMzQhJgmPUldJZFAAvzLbD1vW1YwFZiuxdKmqa1HEpEdqHlySXjJqve3zqj5nE1Vevw/cZPvnAJLOA54O3C5pSV46WQLckW9fJq4NAHmQtRUA22ungaKTOtN0LGyid/YokStjR9dLLVVUg6XKsHufap+cuVCyF/yceeGIdQYqF5QmJPVm4GBJ2+ax/A8lCy+wEjgq3+Yo4Mv5/5XAkZK2lrQ3sJSGA6DVxTQZUJMvbL9TwDTtJU3HCauyx3eVx5s3ZUuc/b/968ser2tUOcBWZTbZR4+YXlvHOecuKLYvBs4FfkDmMrwZWaniVOC5km4AnpvPY/sa4BzgWuArwPG2183b7rppm5j0M0vje0rXMY5JIgi36bpgfD+j8CgcTVWxvKpG0q6SPi7pgnx+X0nH1HGuRir9bL/T9m/bfqLtV9l+wPZdtg+1vTT/vbuw/Sm2H297H9sXNGFznbT55euv2mnrtbQt86+LcfevbOfWRUvPCQfYmjefIgs0uVs+fz3wljpOlFYrUlCaFAbAqkpMUqn26mokgLooIz6LghFr129WamqAnW2fA6yHh6MZ11LLE4LSMCnFfipLL+Nta2lkHFVmhG3JVCdxkuiv3uzqczApqbWhFPiVpMeQe8dKOhi4t44TRfj6Bhnmq586ddqbSkfIYRnmNLalfI+nHTlyWOP8NMfsBCZFL68ebyNzbnq8pH8HdgH+uI4TRQmlQaoaC6LpKqOqSP0aZrlPKWay/aJZtpQyyKOtV2pN8TrnQcptKLZ/APweWfeMPwX2s31lHecKQWmY/pdzUV9I6P61p3p9ZYW8KBjR7rQpqQqKpOOB7WxfY/tqYDtJf1bHuUJQWk7qX/WTMKoaJaif/hLHsHXBphixbv1mpaYGeL3tex62NYvi/vo6ThSCkhCTvLRdz3RTur4uZqb90RW6eI3zJuFG+c3yTuQA5NHat6rjRNEonxCTtqeklOl2nSrauXrHSeG+DerpnoIreltx2o3yq4BzJJ1B1txzHFkn8coJQQmCiun3CktNTIax6OFTZsXpCsqJZEFz30AWbPdC4Mw6ThRVXgkSL3JGypnvMIaJSduIZ3BSyjXIN+Tltd72Gbb/GDjF9sfqCl8VgpIYXXEBroq2pEXvvhW9ngb1X0nZKyo8DWfDVqmpYWopmfSIKq9EaEvGGQymKz3HQ0ymw4Z16xsXi42QtEUeZmWjxXWeM0ooCbAoA09NS5sz6B69EkyxJJMaUTqejQS9vAYN8/HuOk8YJZTEiBe6/fRHQBg3euWi0UVXZZNko/wmBtn+0sgdpLeVOO6vbH9s0IoQlIZZ5Iylq3QhnH8Zph1QrZskOWLjLqMEwvZpAxb/BXA6o6vGjgNCUFJjWGC97r503abfPbiL93FQpOGyotL1UR6d3qDjmwPbMVm7yWdsv2fUBpIeOWxdtKEkRHjZtJe2lkrKeJ0NizA8qkpvGtqUboNI0MvrNtvvsf3uQdPga/BfjjvoqG2ihDJHxg0f29Zw9otOF+7XMFEpG1+tqnaRVB0WxpF5eSX3fT61eknaFfhrYHfbh0vaFzjE9sdH7ZdcCiw6Xcicgu5TFJBpXKa7GO7eLjfNkUNn2PdTZCFbluTzpYYNDkGpmGEvVAhFN+nyfR03eFaZMVQGVZMNimTcBYFJrcrL9t0z7D7VsMFR5VUhs4pJVHe1iy7fq0mubdxz39ZqrEkwSfSCr5Kphg2OEkpF9DfKdjmzCbohJtNk8sOuuyuRAmbBJad5I+m9ZZb10T9s8KeBN407VwhKRYx6Oad1qww2kFLapGTLLFRRop50n85i8HqVmsog6XBJP5K0WtJJQ7Z5lqTLJV0j6ZsjDvfcAcuOGHk5Uw4bHIJSIbM0UAZBE0wjBGXchcf1OekiVbWh5ANgfYQs098XeHnuZVXcZgfgo8CLbO8HvHTAcd4g6SpgH0lXFqabgDJjyh8I7A88Jbfh1eN2iDaUipnkxQkhmYxFqItvgnGN73Wep65zNUGFHlwHAqtt3wgg6WxgOXBtYZs/Ac6zfXN2bt8x4DifBy4A/gYolnLuG9dgL+kzwOOBy9nQGG+yqq+hhKBUyCwlkkX4ggvSpIkMvWsOKBPG8tpZ0mWF+RW2VxTmdwduKcyvAQ7qO8YTgC0lfQN4FPBB2xtl9rbvJWtIfzmApMcC2wDbSdquJ0ZDWAbsa08mkyEoFTLoJRk24FKXXqZ50mRgwbhnGVWkQ+fS0kB5QbnT9rIR6wcdqD9j3wJ4Kllfk0cA35X0PdvXb3Iw6YXAacBuwB3AbwHXAfuNsOFq4DeA20ZsswnRhtIAnXuZGiDSsBqmEeYq+4t0pboLKu3YuAbYszC/B3DrgG2+YvtXtu8EvkXW3jGI/wUcDFxve28yEfr3MTbsDFwraZWklb1pnOFRQgmCEnRRwGbNzIc1zvcft4tptynlPbhKcCmwVNLewM+AI8naTIp8GfiwpC2ArciqxP5+yPEesn2XpM0kbWb76yXcht81jeEhKHNkMV6s+VGMfVbnl27ct9GM8/oqru9SiWQTKmqUt71W0hvJQp9sDnzC9jWSjsvXn2H7OklfIfPWWg+cafvqIYe8R9J2ZKWYz0m6A+gfybHfhm/m8byeli+6ZEjD/0aEoFRM1xobu05/yPl+etU7cU/HM04sZhGT5AflcrUDbNk+Hzi/b9kZffPvA95X4nDLgf8C3gq8Ang0MC5E/cvyY3+DrE3nHyT9he1zR+0XghK0mmlLKYvYT6JIshlzm0lvPBQAbP8q/7seOKvkbm8HntYrlUjaBfgaMFJQGmmUl7SDpHMl/VDSdZIOkbSTpK9KuiH/3bGw/cl5j9EfSTqsCZvLsmgZU0qUTfsq4lS1lZQCME4Toijt+6GSUyvYrK+K6y5K6EVTXl4fJPNQ+G0yz4TryDreXGR7KXBRPk/eQ/RIMhe3w4GP5j1JkyPth727jBu/PWg/g6IcJxczb33JqR18Jffweo2k1wD/Sl8V3CDmLiiStgeeCXwcwPaDtu8hq+frFcfOAl6c/18OnG37Ads3AavJepIGwcPUmbEklWnNQColkx6TfgQMCnvf239QcNa53rdeP5QyUwuw/RfACuBJZB/9K2yfOG6/JkoojwN+DnxS0n9IOjMfo3hX27cB5L+Pzbcf1Gt090EHlnSspMskXfYQD9R3BQNIvtFwAYhSynC68lyOGia76RJLggNsASDpd/NmhOsl3SjpJkk3jtvP9hdtv832W23/U5lzNSEoW5AFGzvd9gHAr9g4zkw/ZXqNZgvtFbaX2V62JVvPbmlJ2jqeeNeo6z50JTNOkWk+Aor3Y1yb0FxFJtX49Vlt0GnAM8jcgJexwR14IJJekrdn3yvpl5Luk/TLcSdqwstrDbDG9sX5/LlkgnK7pCW2b5O0hCxEQG/7cb1GSzGuA9a0xwnSpKp+D+Nci4MASLk6617bF0y4z98CL7R93SQ7zb2EYvs/gVsk7ZMvOpQsiuZK4Kh82VFkPUHJlx8paeu85+hS4JI5mpxuI2AwFcNEpszX7qj9g+mYR3rO492Vy00N8HVJ78u9aZ/Sm8bsc/ukYgLN9UN5E1mPza2AG4GjycTtHEnHADeTx/fPe4ieQyY6a4HjbY8d27jIqIGEJnmYo45+OIv09d4f8LN/WTAd0wp2EmlvQXWhV6qmF6m4GJDSwHNG7HOZpC8AX4INDdK2zxt1okYExfblbHxxPQ4dsv0pwCl12lSGMg/uImWsRVK95mEfDcMyrXHXUby/g8RkUe9/UySV1ul2bHz2FLttD9wPPK94KCA9QZkXk3ZgG1UVktSDG0xElVVV4xr+4zmZjv6hHcrcq2nSutZ4YokKiqRHA+8k664B8E3gPfl4KQOxffQ05+ps+PonPOn+UtuVrVst8xBG3Xra9LeDRebfTpJty0zXy+sTwH3Ay/Lpl8AnB20o6dhxBxu1TadLKPOm89FUO8Swr+GqQ6+nULptyzM5Lq2aTseR9Do2psnjbf9RYf7dki4fsu1Jku4ccSwBJ5B1etyEUoKSD1f5SeDztn9RZp9FpRisMGgHgxrZ6zh2U7TpQ2fcqKdVUFd6NOTBVYb/kvQM29+GrKMjWfThQXwTeOGY43112IqyJZQjyTyxLi2Iy4WTjjecOlWNr5FCJhJMR5sy367S2g+ydHPDNwBn5W0pAu4GXjNow2nbTnqUakOxvdr224EnAJ8nq5O7WdK7Je00iwGpUUX9emRI7WaS2FLFKWVamUEXqDqN64iqkGo/FNuX296fLC7X79g+wHYtD0TpNhRJTwJeCxwBfBH4HFlX/n8DnlyHcfOm/8tokuqr/u1a+5UVAONLKmVdkVN6BtpU+mqLnRuRWBuKpFfa/qykt/UtB8D2aVWfs2wbyveBe4AzgRNt9zq6XJzXxyXH9VduW8lxphGTCNXRDfoz4Gky5HgGqqOOdKzsmM15cI3ikfnvowasq8XasiWUVwEHAHsDJxYU7j22X1KHYU0wbYN6ZBjdpYp7m5KotKmUUjcb90+6YfYDJiYotj+W//2a7X8vrhtXEJB0Allb+X1kBYkDgJNsXzhqv7L9UE4ja/lfSxYduDd1jmm/QoPFYFphSOkZmbTDbwq0YbwbrS83NcA/lFxW5LW2f0nWU34XMqesU8edqGwJZQ/bh5fcNgnKdmyE2UJoFMdnqNvdMWg3VfdxmYVB1Xkw/DlOSRCTJbESiqRDgKcDu/S1o2wPjBv1ttcg9Hzgk7avUK9qagRlBeU7kn7H9lUlt28V074wg17IIChL01Vhk4SOqbOvzjja8G41GEl4FFsB25Hl88V2lF8Cfzxm3+9LupCsmeNkSY+ixADGZQXlGcBrJN1EFnlSgG0/qeT+yVMUlWk8fILFIO79fGmDmDxMYl5etr8JfFPSp2z/tOx+eUnkHWRVXTfavl/SY8iqvUZStg3lCLJxSJ5H1pbyAsb3pmwd47yzpvWFj0yoG9RxH+PZGE6rxARSjuV1pqQdejOSdpS0atjGeYf1L9n+ge178mV32b5y3InKdmz86aCpzL5tpa5OVEFQJJ6NwbQxXVLt2Ajs3BMGgDx81mPH7PM9SSOHCR5EZ6MNT8uouuJx9cvjek33r4+v0/ZQ171qS0/7edJGMcFJe3mtl/SbvRlJv8X4stKzyUTlx5KulHSVpLEllIg2XANlhKc3HyP+tYN5eDrFs9Dya0+vUb7H24FvS/pmPv9MYFyY+iOmOVEIygCKmcewXvBlKeuGGR5jaTPPEkRK7sWDCDfiISQqKLa/ko8hfzCZQ9VbbY8KUY/tn0p6BrDU9icl7ULmMTaSqPIawrQv8TD//kHbDVsXVSBBkRSfhaTDoDREwm0oAOuAO4B7gX0lPXPUxpLeCZwInJwv2hL47LiThKCMYNKR4aaJVDxq22hrCXrER0YwLZJeB3wLWAW8O/9915jd/hB4EXlEFNu3Mjgm2EaEoJSkyVDlbf9y6wJxDzYl0qSPdN2GTwCeBvzU9rPJ4nL9fMw+D+buwwaQ9Mgx2wPRhlKaMiWJQduW6Q09qLF+0HGDIEgUN+bBVYZf2/61JCRtbfuHkvYZs885kj4G7CDp9WRDl5w57kRRQqmAYtXYMHGY9Di9+TIlohCe+dD0F3nT56+b1j/H6ZZQ1uQdG78EfFXSl4FbR+1g+++Ac8nGvtoHeIftD407UZRQKqSOUeDKHKfrGU0Q9zh1RJKxvACw/Yf533dJ+jrwaOAro/aR9F7bJ1IYP76wbChRQmmAcZ0fg3SJjD0YSoIlFEmbSbr6YRPtb9peafvBMbs+d8CysX1TQlDmSHjqBNMQItaCD62SLsNlSzGSDpf0I0mrJZ00YrunSVonaWD0YNvrgSuKPeXHnPcNkq4C9sl7yPemm4DoKd8k0/Z8bjqseTCa6NiXEWnQR0WN8pI2Bz5CVkpYA1wqaaXtawds914yN+BRLAGukXQJhYERbb9owLafBy4A/gYoCtl9tu8eZ3sIyhwYVb3V5DgTwfSEqMyHNqVxhW0oBwKrbd8IIOlsYDlwbd92byJrNB8YxDH36HqArO9JKWzfS9b58eW5YO1KphPbSdrO9s2j9g9BqYFJXoJhpZgopaRPuHdXx7DhIlpFeUHZWdJlhfkVtlcU5ncHbinMrwEOKh5A0u5knQ+fwxBBAb4LPAV4ne1XlbYuO/4byTo/3s6GspeBkWNghaAkQutenuBh6iytLMJHRSeucbIG9zttLxuxftBIXf1H/wBwou11I0bm3UrSUcDTJb1kkwPa542w4S3APrbvGrHNJoSgVEyVGUuUUoKUqUJIyz7fbXgPKqzyWgPsWZjfg037jSwDzs7FZGfg+ZLW2v5SYZvjgFcAO7DpgIgGRgnKLWRVXxPRmKDk9XOXAT+z/QJJOwFfAPYCfgK8LB8IBkknA8eQBTh7s+1xjVBzp87xMtrwMgXBwlOdoFwKLJW0N/Az4EjgTzY6lb1377+kTwH/0icm2P42Wdj6y2x/fEIbbgS+IelfyYZ97x3ztFE7Nek2fAJwXWH+JOAi20uBi/J5JO1LlqD7AYcDH83FKAlSdAVOzZ4gmJQ2PsNVDbBley3wRjLvreuAc2xfI+k4ScdNatcUYgJwM1mnxq3IgkL2ppE0UkKRtAfwB8ApwNvyxcuBZ+X/zwK+QRY+eTlwdu6tcJOk1WReEN+do8kDmddDP2kpJUo07aZN9y+83XIq7rRo+3zg/L5lZwzZ9jXVnfnhY74bsqCQtn81bvseTZVQPgD8JRt7bu9q+zaA/Lc35vEgj4fdBx1U0rGSLpN02c/vWle50U0SL233mXS4hC5TRxijOtEEUxuQdIika8lrkSTtL+mj4/abu6BIegFwh+3vl91lwLKB3wK2V9heZnvZLo9JplasMkJUuksbMs2q6VycugRDr/SQtLukp0t6Zm8as8sHgMOAuwBsX0E2dPBImqjy+l3gRZKeD2wDbC/ps8DtkpbYvk3SErLRxaCcx0MQNEYI/eR0TkxINzikpPcC/52sY2Sv6sZkg24NxfYtfS7JY6t95l5CsX2y7T1s70XW2P5vtl8JrASOyjc7Cvhy/n8lcKSkrXOvh6XAJXM2OwiCIbQt46+NdEsoLybrU/J82y/Mp0FhV4rcIunpgCVtJenP2diJaiApBYc8FXiupBvIYticCmD7GuAcMnX9CnC87SQaSOJFCoIAeHiArSq8vGrgRrIx4SfhOOB4svbqNcCT8/mRNNqx0fY3yLy5yHtkHjpku1PIPMKCoHN04cNkEm+vLlzvQBKt8gLuBy6XdBEb9yl587AdbN9J1ilyIqKnfBA0SGcz1wUk1TYUsmaDlZPsIOks4ATb9+TzOwLvt/3aUfuFoMyJCCQY9NM1MSlTSunaNW9EooJi+6wpdntST0zyY/xC0gHjdgpBqZFBL88swtLplzHoPF1/flMtoUhaSja+yb5knrUA2H7ciN02k7RjIfzVTpTQi5Qa5TvFuJent77rL1kwmEW67wvRYdNk3bTLTPPnk8DpwFrg2cCngc+M2ef9wHck/ZWk9wDfAf523ImihFIR07ww0ZDZfiYpaS7CPVyEaxyESLeEAjzC9kWSZPunwLsk/V/gncN2sP3pfMyW55Bd3kv6R4wcRAjKjMz6As2zbSUiFwdBjaQrKL+WtBlwQz5w1s/YENpqKLmAjBWRIlHltUBEIL9mibTvNrJLTQ3wFmBb4M3AU4FXsqETeaVECSURhmX2dZUoorQyO5MKRKR3h2kwTtc4bF8KkNV4+eg6zxWCkhDzynBCTKqhbIkv0noxSLUNRdIhwMeB7YDflLQ/8Ke2/6zqc0WVVxDUSIjJ4pBw6JUPMEXk4GmIEsqCERlcdUSbSLARiZZQYLrIwdMQgpIovcyqKgGo+njBaCKdFwynW+VFX+Rgssb5sZGDpyGqvBaIyOSqZVR6RullAUk3fP1UkYOnIUooiVJ1ySSYHyHci0fKHRunjRw8DVFC6TBFMYlMLgjqRetdapq7XdJZknYozO8o6RN1nCtKKB0lxKQZIq0XlIT7oTBl5OBpCEEJgimIqsSgn4ZcgsswVeTgaQhB6SjxpVwvEcYm2IR0Syi9yMHn5vMvpaYRcKMNJQimYJiYhMgsLnK5ad7Y/jTwR8Dt+fQS2+PC109FlFCCoEKiZLigGGgm8GNZtiRzRnP+vxaihBIEQaeZV6kx1dArkk4APgfsTBa2/rOS3lTHuUJQgmAKRg3vHKTFPO5Lrx9KilVewDHAQbbfafsdwMHA6+s4UVR5BcEU9H/1hpgsOHbKVV5i49hd6/JllRMllCCYgn4BGVWtEg31i0HCJZRPAhdLepekdwHfIwtnXzlRQgmCChhVQonSy4KQaAHF9mmSvgE8g6xkcrTt/6jjXCEoQVABMWhZkGosLwDbPwB+UPd5osorCKZkkmqvoOMYWOdyU4cJQQmCGQhRCXok3IYyN6LKKwhmJKq6AqBSLy9JhwMfBDYHzrR9at/6VwAn5rP/D3hDPrRvo4SgBEEQVEBVpQ9JmwMfAZ5LNiDWpZJW2r62sNlNwO/lkYOPAFYAB1VjwfRElVcQBMGslB2tsZzoHAistn2j7QeBs4HlG53O/k4vejCZG/Aes1/E7MxdUCTtKenrkq6TdE0eFgBJO0n6qqQb8t8dC/ucLGm1pB9JOmzeNgdBEIxCgNa51ATsLOmywnRs3+F2B24pzK/Jlw3jGOCCSi9oSpqo8loL/A/bP5D0KOD7kr4KvAa4yPapkk4CTgJOlLQvcCSwH7Ab8DVJT7C9bsjxgyAI5o7Kt6HcaXvZqEMNWDbw4JKeTSYozyh78jqZewnF9m25TzS27wOuI1Pf5cBZ+WZnAS/O/y8Hzrb9gO2bgNVkRcIgCII0qLbKaw2wZ2F+D+DW/o0kPQk4E1hu+67pja+ORttQJO0FHABcDOxq+zbIRIcsKiZMUPyTdGyvGPnzu6IAEwTBvPCGeF7jpvFcCiyVtLekrchqaFYWN5D0m8B5wKtsX1/55UxJY15ekrYDvgi8xfYvpaGxykoX/2yvIPN2YNn+23Tc4zsIgpSoysvL9lpJbwRWkbkNf8L2NZKOy9efAbwDeAzw0TzvXDumGm0uNCIokrYkE5PP2T4vX3y7pCW2b5O0BLgjX16q+BcEQdAoFfZDsX0+cH7fsjMK/18HvK6yE1ZEE15eIot0eZ3t0wqrVgJH5f+PAr5cWH6kpK0l7Q0sBS6Zl71BEARj8UReXp2liRLK7wKvAq6SdHm+7H8CpwLnSDoGuBl4KUBe1DsHuJbMQ+z48PAKgiA5uq0VpZi7oNj+NsMHdzl0yD6nAKfUZlQQBMGMTOA23Fki9EoQBEEVhKBE6JUgmAcRhbjjGFhfcuowUUIJgpoJMek+wlHlRQhKENRCbwTHGMlxgVjf8eJHCaLKi/iCDOohxGSBiCovIATlYYqiEgITVEU8S4uD7FJTl4kqrwLDRKX4lTls+aTHn2b/oB2EiNRH0u9Qx8WiDAshKINEYJKXvrftsPHDx40rPuyhT/rlCCZmUBVX3NPqSFuoSwd+7DRyRxNh2f7beMerXti0GTMRmVH7iI+EehgnJrOm8+ZLbvj+LMEVH/2IJT7kca8tte2qa/96pnOlzEKUUNpMNOwGi8o0tQj9zPPd6Xr7SBk6KyjXX7ktBw2NiN8Oei/JsKq1IOgiVVZtlXl3sm1umP1kISjh5dUmDttt/4enIB3ifkzHoHSrKy2HHbey8xlY73JTh+lsCaXrzOJtFlRHuJtPRn+JYZ5VuoPOXR3RKA8hKJ0g2lnmT4hHxqB0GPQs9m/XpBDXdr4QlBCUrhCiMl96YVV6jKpSaeN9mSXTXUixNbCu493gSxCC0iHCZTWYhIXM+GvD4BCUaJTvMJFh1Efb07bt9ieJXW7qMFFC6ThtrXKpmqpLb/1VXsGC0/PyWnBCUDpOiMlgqhDatolKm2xtJR0vfZQhBCVYWKbpXd3WTLmtdreKEJQQlK4SJZONmaQ0MWsw0ZRoq92tw4Z165q2onFCUDpIV8SkTGZY9lqnzVgjQw5KEyWUEJSu0WYxmSbzjgx/OJE2cyYEJQSlS7RNTCLDq4dI1ybofpyuMoSgBHNjUTO6JuJVBXPG4OjYGIIS1E9kcvMZgiDSuWEi9EoISlAvkcltTNWllUjfRLBhfQhKhF7pEE1mLvMc26LtVDWmTaRvYkTolSihdI3URCUYTqRXt3CUUEJQgiAIZqf7pY8yhKAEQRDMSgSHBEJQgiAIZsaAI/RKexrlJR0u6UeSVks6qWl7giAIHsb5AFtlphKMy++U8aF8/ZWSnlL5NU1BKwRF0ubAR4AjgH2Bl0vat1mrgiAINuD1LjWNo2R+dwSwNJ+OBU6v9mqmoxWCAhwIrLZ9o+0HgbOB5Q3bFARBsIHqSihl8rvlwKed8T1gB0lLqr2gyWlLG8ruwC2F+TXAQf0bSTqWTK0BHviaz716DrZNws7AnU0b0UfYVJ4U7QqbyjHOpt+a5eD38YtVX/O5O5fcfBtJlxXmV9heUZgvk98N2mZ34LaSNtRCWwRFA5ZtUnbMb8oKAEmX2V5Wt2GTEDaVI0WbIE27wqZy1G2T7cMrPFyZ/K5Unjhv2lLltQbYszC/B3BrQ7YEQRDUSZn8Lsk8sS2CcimwVNLekrYCjgRWNmxTEARBHZTJ71YCr869vQ4G7rXdaHUXtKTKy/ZaSW8EVgGbA5+wfc2Y3VaMWd8EYVM5UrQJ0rQrbCpHijYNZFh+J+m4fP0ZwPnA84HVwP3A0U3ZW0SOcAFBEARBBbSlyisIgiBInBCUIAiCoBI6JyhNhWiRtKekr0u6TtI1kk7Il+8k6auSbsh/dyzsc3Ju548kHVajbZtL+g9J/5KQTTtIOlfSD/M0O6RpuyS9Nb93V0v6R0nbzNsmSZ+QdIekqwvLJrZB0lMlXZWv+5CkQW6ms9j0vvzeXSnpnyTtME+bhtlVWPfnkixp58Kyudi10NjuzETWgPVj4HHAVsAVwL5zOvcS4Cn5/0cB15OFTfhb4KR8+UnAe/P/++b2bQ3sndu9eU22vQ34PPAv+XwKNp0FvC7/vxWwQ5N2kXUKuwl4RD5/DvCaedsEPBN4CnB1YdnENgCXAIeQ9Ve4ADiiYpueB2yR/3/vvG0aZle+fE+yBu2fAjvP265FnrpWQmksRIvt22z/IP9/H3AdWSa1nCzzJP99cf5/OXC27Qds30TmrXFg1XZJ2gP4A+DMwuKmbdqeLDP4OIDtB23f07RdZF6Pj5C0BbAtmV//XG2y/S3g7r7FE9mgLATH9ra/6yzH/HRhn0pssn2h7bX57PfI+kHMzaZhduX8PfCXbNzRb252LTJdE5Rh4QjmiqS9gAOAi4FdnfuH57+PzTebl60fIHu5ikGEmrbpccDPgU/mVXFnSnpkk3bZ/hnwd8DNZOEr7rV9YZM2FZjUht3z//OwDeC1ZF/2jdsk6UXAz2xf0bcqlbTqNF0TlMbDEUjaDvgi8Bbbvxy16YBlldoq6QXAHba/X3aXAcvqSL8tyKoqTrd9APArsqqcxuzK2yWWk1WH7AY8UtIrm7SpBMNsmJttkt4OrAU+17RNkrYF3g68Y9DqpuxaJLomKI2GI5C0JZmYfM72efni2/NiNfnvHXO09XeBF0n6CVn133MkfbZhm3rnWWP74nz+XDKBadKu3wdusv1z2w8B5wFPb9imHpPasIYNVVC12SbpKOAFwCvy6qKmbXo82QfBFfkzvwfwA0m/0bBdC0PXBKWxEC25Z8jHgetsn1ZYtRI4Kv9/FPDlwvIjJW0taW+ycQ0uqdIm2yfb3sP2XmRp8W+2X9mkTbld/wncImmffNGhwLUN23UzcLCkbfN7eShZO1ijaVU4V2kb8mqx+yQdnF/Lqwv7VIKkw4ETgRfZvr/P1kZssn2V7cfa3it/5teQOcr8Z5N2LRRNewVUPZGFI7iezIvj7XM87zPIispXApfn0/OBxwAXATfkvzsV9nl7buePqNmzBHgWG7y8GrcJeDJwWZ5eXwJ2bNou4N3AD4Grgc+QeQTN1SbgH8nacB4iyxCPmcYGYFl+HT8GPkweFaNCm1aTtUn0nvUz5mnTMLv61v+E3MtrnnYt8hShV4IgCIJK6FqVVxAEQdAQIShBEARBJYSgBEEQBJUQghIEQRBUQghKEARBUAkhKEEQBEElhKAEQRAElRCCEnQSSXspG2flfysb4+RCSY+SdKmkZ+Xb/I2kU5q1NAi6QwhK0GWWAh+xvR9wD1nwx9cAp0t6LnA4We/4IAgqIAQl6DI32b48//99YC/b15CFVfln4LXOxs0JgqACQlCCLvNA4f86srD5AL9DVmLZdd4GBUGXCUEJFgpJLyELtvhM4EPFsdCDIJiNEJRgkdgZOJUsKu31ZJFlP9isSUHQHSLacBAEQVAJUUIJgiAIKiEEJQiCIKiEEJQgCIKgEkJQgiAIgkoIQQmCIAgqIQQlCIIgqIQQlCAIgqAS/j99s5k/MulhfAAAAABJRU5ErkJggg==\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 2 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "mask_orig['mask'].plot()" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 10, | |
| "id": "4af0fa2c-5bf5-418b-8909-f8778d66e8d6", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/plain": [ | |
| "<matplotlib.collections.QuadMesh at 0x146098d0f3d0>" | |
| ] | |
| }, | |
| "execution_count": 10, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| }, | |
| { | |
| "data": { | |
| "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZQAAAEKCAYAAAA1qaOTAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8qNh9FAAAACXBIWXMAAAsTAAALEwEAmpwYAAAwvklEQVR4nO2deRQlVXXufx+ziAgIkmZIQFdLHhgRbRmMy6hEAaO2MdGHcUBECQYV9SUBnms55ZFgjESNCvbDAacgD4l2ErBR4vCMyqBhRqEFhRYCAoI8iUB3f++PqktX375D3Xurbp2qu39r1bq35l2nqs5X55x99pFtgiAIgmBWNmvagCAIgqAbhKAEQRAElRCCEgRBEFRCCEoQBEFQCSEoQRAEQSWEoARBEASVEIISBEGQEJI+IekOSVcPWS9JH5K0WtKVkp4ybxuHEYISBEGQFp8CDh+x/ghgaT4dC5w+B5tKEYISBEGQELa/Bdw9YpPlwKed8T1gB0lL5mPdaLZo2oC62EpbexseOdW+T3jS/RVbExS5/sptN0nj66/cdqP54vr+dSkwyr5B6wZdbzxn1THrM3Ifv7jT9i7T7n/Ysx/pu+5eV2rb71/5wDXArwuLVtheMcHpdgduKcyvyZfdNsExaqGzgrINj+QgHTrdzldlP6tuvaI6gwIO221/AO65bXC69taPSvfeNk2zatUGG/ttKq4rbhPPUz0cttv+HKTZjvE1n/vTWfa/8+51XLxqj1Lbbrnkx7+2vWyG0w262iRiaHVWUKogMoFqKZOWqYvJqluv2FRA+mzuf24GzQ86bjAdg+7J/DHrvH5eJ1sD7FmY3wO4dV4nH0W0oQTJ0MtUh2UOTWe6vYxrkIAM+l/cb9i6/vXNZ4zBNBhYj0tNFbASeHXu7XUwcK/txqu7IEooQUKUqfJqmlFi0ls/SBR6y4aVbgaJUsrpEGzKeqopoUj6R+BZwM6S1gDvBLYEsH0GcD7wfGA1cD9wdCUnroAQlJLMUv0VVRwbGJVZ9jLjYVVEKaRZ8V6uuvWKsQJTxblSuO5ZqbP6OIVSnTEPVVTlZfvlY9YbOL6Sk1VMVHmNoVgNMc2DO+iLtAsZxLSMqv4Zl75NZxzDzt97RkbZN+yel7mmpq87dVJ4nwysw6WmLhOCMiHTvtyLLiRFhqVFb1n/ujalW+/ayohMj66LyjSlkzZe7xzbUJIlBGUKJn3Y25QhNsWgqq7UmMW2KjLINmayMH26lb3eFNLFwDq71NRlQlCmZBKPnBQe+NQY5Rk1rM0pBbHp2THOlnFuwtOeO9hAal5x60tOXSYEpUa61KhaJdO0SaWWeRSZRTzKCFMqYpoaKaWJS7afdL0NJby8ZmRUNU0aHa7Soz9d+vt3jPIAa5Ky3nq9aynaPKqjY5mMcdTHSXgRNo8ND3VbK0oRglIBozqsxYtdjtR7yI9jEsEb90yMOs4gIS6zXVdJ4UMjQ6wbGBFlsahNUCR9AngBcIftJ+bLdgK+AOwF/AR4me1f5OtOBo4B1gFvtr0qX/5UsnDOjyDr0HNC7oedNF1+iauizWlU1vZJ3IUnzRwnacNrc1q3AQPrk8+V6qfONpRPsWlM/5OAi2wvBS7K55G0L3AksF++z0clbZ7vczpZzP9e/P9R4wQELWFURjtr35+6KboF9y8v/k5KitcalGddXkoZN3WZ2gRlSEz/5cBZ+f+zgBcXlp9t+wHbN5GFFDgwj/G/ve3v5qWSTxf2CTpOGzLYfmEp2x4SJYZukXVsDEGZdxvKrr0gZrZvk/TYfPnuwPcK2/Xi+z+U/+9fPhBJx5KVZtiGZsfQiExjMtogHj3a2KdkXITkWY5Xx3M+LAxPqhh4yOE0m0oKDIvvP1Hcf9srbC+zvWxLtq7MuKAZUsxIZgmhMsl2dVKlmAQZRqxjs1JTl5l3CeV2SUvy0skS4I58+bD4/mvy//3Lg46RjrfOaGa1cVwo+7qpord/ikKfAuvd7eqsMsxbLlcCR+X/jwK+XFh+pKStJe1N1vh+SV49dp+kgyUJeHVhnyBoLU24lFcVOmaSOGVV0IYPjWhDyajTbXhQTP9TgXMkHQPcDLwUwPY1ks4BrgXWAsfb7g3Q/AY2uA1fkE9JEl9u1dP0F31ZUr/3dbZz1HXtZUqt6TwXYl20odTq5fVy20tsb2l7D9sft32X7UNtL81/7y5sf4rtx9vex/YFheWX2X5ivu6NKfVBKYbFaEvGlzJ1ji1SN7NEoW4Do1yi23Sf6iIbsXGzUlOXiZ7yMxBjnQSpU+XzOKjEEM97hi0e9ObjN+w43ZbLORNfatWyCJlV3c9M1cdvqt2nDe/WelRq6jJRQqmY8ISZjTZlIItKU892yn1Sskb5+D6PFKiJaTPElMO0z5NUM4620bV0TPfdyBrly0xdpttX1yBVdB4LcWkXca/qYdbOpPMgGuUzosqrJopVX5NWgxWrfabZv2uklHEEzdCGZ39ddGwMQambQUPdjno5Bg281C8wi0SISdAGjHjIkZ12u/zVMRZNTNpIHeOZBOnTa5Rf9Fhe3b66lrLobScpj4XSJuIDZH4Ysc7lpi4TZbQ5U+Ylj4ygvWIS925x6XqDexkiBRpiUIZZ9Oxqa4Y6K22/7knsr1N8ui5sqT0nNuE2TAhK44x6MRZJXBblOvvpesZfBW14LrJG+c1LTV0mqrzmzCCvrza8MMOY1fOszdc+iBQEIgUbilTp9p6yp2PXG9zLEILSElJ8kaYVxNSuownqFNJUnpX+a5zWrjZEoTaKAbaIKq9WMWiQo6qON+992zRe+CQ0Ka6pRLuuq/oyRSEpEm7DUUJpHVV99cH0mVmVL3bqmURddLFBvsyzOOh+jxpjpS3BQg2s73iDexkiBRaAEID50lQapeDYMM35m7a5GsoN/xtDAAfJUybmV9kvw2C+VJmZplwymfQYxf9tEBxD5z24yhAllI4wqH2lyvG4U/j6bRNl0qpL6VnFtQzygGwLtljvzUpNXabbVxdUJipRmqmWWTPMYuDQphvj6zj3oOGGU6fKjo2SDpf0I0mrJZ00YP2jJf2zpCskXSPp6MovaApCUILSpY8QlcmoK0Psj0LdNHW7QLeBbDyUaoYAlrQ58BHgCGBf4OWS9u3b7HjgWtv7A88C3i9pq0ovagpCUIKHCVGZnUnHvFlUmi5VVU+lIzYeCKy2faPtB4GzgeV92xh4lCQB2wF3A2urvKJpCEEJNmGcsHQrI6ieMunTli/vccwy1DVsLCxFF+G2pU/mNqxSE7CzpMsK07F9h9sduKUwvyZfVuTDwH8DbgWuAk6wvb6eqytPCEqwEaN6v7ftJZ83k/SZmHffoRQpplNb+psMY8JYXnfaXlaYVvQdblC9mPvmDwMuB3YDngx8WNL2FV/WxITb8IJTpvGzzd4386Tn+tpL01GZf1096kd1CKxajMqWZBflmakwfP0aYM/C/B5kJZEiRwOn2jawWtJNwG8Dl1RlxDRECSUYWJ/d3yegS1/GKTBJW8skaT+q2mieGfuiPS9Z+PrKBti6FFgqae+8of1IYGXfNjcDhwJI2hXYB7ixwkuaihCUBWSYq2nZL+quZxaD0qVsxl5lZN2qacKmReqPM0EbykhsrwXeCKwCrgPOsX2NpOMkHZdv9lfA0yVdBVwEnGj7zpourTRR5bWAjKqOmXR51xh3nVVGvq2rr8Wg4867umuS47Sxz0k/WbTh6r7PbZ8PnN+37IzC/1uB51V2woqIEsqC0vYXuGqGldjGZXZlw9zMm2FhTGaljAfWpG0nKaTXrGShVzYrNXWZbl9dMJJZXT67xKgv+rraO8ocb9b9512yXLSG+A1E6BUIQekcdWYgbewfUIZRGe+k8c6mOXdX6NK1TENVPeXbTCNtKJLeCryOrKR4FZkL3LbAF4C9gJ8AL7P9i3z7k4FjgHXAm22vmr/V6VOmLnrRX/p+yraZDKvrn9S5oa108UOiSnpeXovO3EsoknYH3gwss/1EYHMyt7iTgItsLyXzWjgp337ffP1+wOHAR/NYNwvPtC6lweRUUYKZ5Nhxr9pHalVekp4yaqrjnE15eW0BPELSQ2Qlk1uBk8mCnAGcBXwDOJEshs3Zth8AbpK0mizWzXfnbHNSTPpVPM8Mqn9slq5/3VYxHsjitj10g0THlH//iHUGnlP1CecuKLZ/JunvyDrm/Bdwoe0LJe1q+7Z8m9skPTbfZXfge4VDDIprs1AME5M6+0CUde8c5bIamWW3KQrrsGgBXX0GDKxNrMHd9rPnfc65C4qkHclKHXsD9wD/R9IrR+0yYFl/XJvesY8FjgXYhm1nMzQhJgmPUldJZFAAvzLbD1vW1YwFZiuxdKmqa1HEpEdqHlySXjJqve3zqj5nE1Vevw/cZPvnAJLOA54O3C5pSV46WQLckW9fJq4NAHmQtRUA22ungaKTOtN0LGyid/YokStjR9dLLVVUg6XKsHufap+cuVCyF/yceeGIdQYqF5QmJPVm4GBJ2+ax/A8lCy+wEjgq3+Yo4Mv5/5XAkZK2lrQ3sJSGA6DVxTQZUJMvbL9TwDTtJU3HCauyx3eVx5s3ZUuc/b/968ser2tUOcBWZTbZR4+YXlvHOecuKLYvBs4FfkDmMrwZWaniVOC5km4AnpvPY/sa4BzgWuArwPG2183b7rppm5j0M0vje0rXMY5JIgi36bpgfD+j8CgcTVWxvKpG0q6SPi7pgnx+X0nH1HGuRir9bL/T9m/bfqLtV9l+wPZdtg+1vTT/vbuw/Sm2H297H9sXNGFznbT55euv2mnrtbQt86+LcfevbOfWRUvPCQfYmjefIgs0uVs+fz3wljpOlFYrUlCaFAbAqkpMUqn26mokgLooIz6LghFr129WamqAnW2fA6yHh6MZ11LLE4LSMCnFfipLL+Nta2lkHFVmhG3JVCdxkuiv3uzqczApqbWhFPiVpMeQe8dKOhi4t44TRfj6Bhnmq586ddqbSkfIYRnmNLalfI+nHTlyWOP8NMfsBCZFL68ebyNzbnq8pH8HdgH+uI4TRQmlQaoaC6LpKqOqSP0aZrlPKWay/aJZtpQyyKOtV2pN8TrnQcptKLZ/APweWfeMPwX2s31lHecKQWmY/pdzUV9I6P61p3p9ZYW8KBjR7rQpqQqKpOOB7WxfY/tqYDtJf1bHuUJQWk7qX/WTMKoaJaif/hLHsHXBphixbv1mpaYGeL3tex62NYvi/vo6ThSCkhCTvLRdz3RTur4uZqb90RW6eI3zJuFG+c3yTuQA5NHat6rjRNEonxCTtqeklOl2nSrauXrHSeG+DerpnoIreltx2o3yq4BzJJ1B1txzHFkn8coJQQmCiun3CktNTIax6OFTZsXpCsqJZEFz30AWbPdC4Mw6ThRVXgkSL3JGypnvMIaJSduIZ3BSyjXIN+Tltd72Gbb/GDjF9sfqCl8VgpIYXXEBroq2pEXvvhW9ngb1X0nZKyo8DWfDVqmpYWopmfSIKq9EaEvGGQymKz3HQ0ymw4Z16xsXi42QtEUeZmWjxXWeM0ooCbAoA09NS5sz6B69EkyxJJMaUTqejQS9vAYN8/HuOk8YJZTEiBe6/fRHQBg3euWi0UVXZZNko/wmBtn+0sgdpLeVOO6vbH9s0IoQlIZZ5Iylq3QhnH8Zph1QrZskOWLjLqMEwvZpAxb/BXA6o6vGjgNCUFJjWGC97r503abfPbiL93FQpOGyotL1UR6d3qDjmwPbMVm7yWdsv2fUBpIeOWxdtKEkRHjZtJe2lkrKeJ0NizA8qkpvGtqUboNI0MvrNtvvsf3uQdPga/BfjjvoqG2ihDJHxg0f29Zw9otOF+7XMFEpG1+tqnaRVB0WxpF5eSX3fT61eknaFfhrYHfbh0vaFzjE9sdH7ZdcCiw6Xcicgu5TFJBpXKa7GO7eLjfNkUNn2PdTZCFbluTzpYYNDkGpmGEvVAhFN+nyfR03eFaZMVQGVZMNimTcBYFJrcrL9t0z7D7VsMFR5VUhs4pJVHe1iy7fq0mubdxz39ZqrEkwSfSCr5Kphg2OEkpF9DfKdjmzCbohJtNk8sOuuyuRAmbBJad5I+m9ZZb10T9s8KeBN407VwhKRYx6Oad1qww2kFLapGTLLFRRop50n85i8HqVmsog6XBJP5K0WtJJQ7Z5lqTLJV0j6ZsjDvfcAcuOGHk5Uw4bHIJSIbM0UAZBE0wjBGXchcf1OekiVbWh5ANgfYQs098XeHnuZVXcZgfgo8CLbO8HvHTAcd4g6SpgH0lXFqabgDJjyh8I7A88Jbfh1eN2iDaUipnkxQkhmYxFqItvgnGN73Wep65zNUGFHlwHAqtt3wgg6WxgOXBtYZs/Ac6zfXN2bt8x4DifBy4A/gYolnLuG9dgL+kzwOOBy9nQGG+yqq+hhKBUyCwlkkX4ggvSpIkMvWsOKBPG8tpZ0mWF+RW2VxTmdwduKcyvAQ7qO8YTgC0lfQN4FPBB2xtl9rbvJWtIfzmApMcC2wDbSdquJ0ZDWAbsa08mkyEoFTLoJRk24FKXXqZ50mRgwbhnGVWkQ+fS0kB5QbnT9rIR6wcdqD9j3wJ4Kllfk0cA35X0PdvXb3Iw6YXAacBuwB3AbwHXAfuNsOFq4DeA20ZsswnRhtIAnXuZGiDSsBqmEeYq+4t0pboLKu3YuAbYszC/B3DrgG2+YvtXtu8EvkXW3jGI/wUcDFxve28yEfr3MTbsDFwraZWklb1pnOFRQgmCEnRRwGbNzIc1zvcft4tptynlPbhKcCmwVNLewM+AI8naTIp8GfiwpC2ArciqxP5+yPEesn2XpM0kbWb76yXcht81jeEhKHNkMV6s+VGMfVbnl27ct9GM8/oqru9SiWQTKmqUt71W0hvJQp9sDnzC9jWSjsvXn2H7OklfIfPWWg+cafvqIYe8R9J2ZKWYz0m6A+gfybHfhm/m8byeli+6ZEjD/0aEoFRM1xobu05/yPl+etU7cU/HM04sZhGT5AflcrUDbNk+Hzi/b9kZffPvA95X4nDLgf8C3gq8Ang0MC5E/cvyY3+DrE3nHyT9he1zR+0XghK0mmlLKYvYT6JIshlzm0lvPBQAbP8q/7seOKvkbm8HntYrlUjaBfgaMFJQGmmUl7SDpHMl/VDSdZIOkbSTpK9KuiH/3bGw/cl5j9EfSTqsCZvLsmgZU0qUTfsq4lS1lZQCME4Toijt+6GSUyvYrK+K6y5K6EVTXl4fJPNQ+G0yz4TryDreXGR7KXBRPk/eQ/RIMhe3w4GP5j1JkyPth727jBu/PWg/g6IcJxczb33JqR18Jffweo2k1wD/Sl8V3CDmLiiStgeeCXwcwPaDtu8hq+frFcfOAl6c/18OnG37Ads3AavJepIGwcPUmbEklWnNQColkx6TfgQMCnvf239QcNa53rdeP5QyUwuw/RfACuBJZB/9K2yfOG6/JkoojwN+DnxS0n9IOjMfo3hX27cB5L+Pzbcf1Gt090EHlnSspMskXfYQD9R3BQNIvtFwAYhSynC68lyOGia76RJLggNsASDpd/NmhOsl3SjpJkk3jtvP9hdtv832W23/U5lzNSEoW5AFGzvd9gHAr9g4zkw/ZXqNZgvtFbaX2V62JVvPbmlJ2jqeeNeo6z50JTNOkWk+Aor3Y1yb0FxFJtX49Vlt0GnAM8jcgJexwR14IJJekrdn3yvpl5Luk/TLcSdqwstrDbDG9sX5/LlkgnK7pCW2b5O0hCxEQG/7cb1GSzGuA9a0xwnSpKp+D+Nci4MASLk6617bF0y4z98CL7R93SQ7zb2EYvs/gVsk7ZMvOpQsiuZK4Kh82VFkPUHJlx8paeu85+hS4JI5mpxuI2AwFcNEpszX7qj9g+mYR3rO492Vy00N8HVJ78u9aZ/Sm8bsc/ukYgLN9UN5E1mPza2AG4GjycTtHEnHADeTx/fPe4ieQyY6a4HjbY8d27jIqIGEJnmYo45+OIv09d4f8LN/WTAd0wp2EmlvQXWhV6qmF6m4GJDSwHNG7HOZpC8AX4INDdK2zxt1okYExfblbHxxPQ4dsv0pwCl12lSGMg/uImWsRVK95mEfDcMyrXHXUby/g8RkUe9/UySV1ul2bHz2FLttD9wPPK94KCA9QZkXk3ZgG1UVktSDG0xElVVV4xr+4zmZjv6hHcrcq2nSutZ4YokKiqRHA+8k664B8E3gPfl4KQOxffQ05+ps+PonPOn+UtuVrVst8xBG3Xra9LeDRebfTpJty0zXy+sTwH3Ay/Lpl8AnB20o6dhxBxu1TadLKPOm89FUO8Swr+GqQ6+nULptyzM5Lq2aTseR9Do2psnjbf9RYf7dki4fsu1Jku4ccSwBJ5B1etyEUoKSD1f5SeDztn9RZp9FpRisMGgHgxrZ6zh2U7TpQ2fcqKdVUFd6NOTBVYb/kvQM29+GrKMjWfThQXwTeOGY43112IqyJZQjyTyxLi2Iy4WTjjecOlWNr5FCJhJMR5sy367S2g+ydHPDNwBn5W0pAu4GXjNow2nbTnqUakOxvdr224EnAJ8nq5O7WdK7Je00iwGpUUX9emRI7WaS2FLFKWVamUEXqDqN64iqkGo/FNuX296fLC7X79g+wHYtD0TpNhRJTwJeCxwBfBH4HFlX/n8DnlyHcfOm/8tokuqr/u1a+5UVAONLKmVdkVN6BtpU+mqLnRuRWBuKpFfa/qykt/UtB8D2aVWfs2wbyveBe4AzgRNt9zq6XJzXxyXH9VduW8lxphGTCNXRDfoz4Gky5HgGqqOOdKzsmM15cI3ikfnvowasq8XasiWUVwEHAHsDJxYU7j22X1KHYU0wbYN6ZBjdpYp7m5KotKmUUjcb90+6YfYDJiYotj+W//2a7X8vrhtXEJB0Allb+X1kBYkDgJNsXzhqv7L9UE4ja/lfSxYduDd1jmm/QoPFYFphSOkZmbTDbwq0YbwbrS83NcA/lFxW5LW2f0nWU34XMqesU8edqGwJZQ/bh5fcNgnKdmyE2UJoFMdnqNvdMWg3VfdxmYVB1Xkw/DlOSRCTJbESiqRDgKcDu/S1o2wPjBv1ttcg9Hzgk7avUK9qagRlBeU7kn7H9lUlt28V074wg17IIChL01Vhk4SOqbOvzjja8G41GEl4FFsB25Hl88V2lF8Cfzxm3+9LupCsmeNkSY+ixADGZQXlGcBrJN1EFnlSgG0/qeT+yVMUlWk8fILFIO79fGmDmDxMYl5etr8JfFPSp2z/tOx+eUnkHWRVXTfavl/SY8iqvUZStg3lCLJxSJ5H1pbyAsb3pmwd47yzpvWFj0yoG9RxH+PZGE6rxARSjuV1pqQdejOSdpS0atjGeYf1L9n+ge178mV32b5y3InKdmz86aCpzL5tpa5OVEFQJJ6NwbQxXVLt2Ajs3BMGgDx81mPH7PM9SSOHCR5EZ6MNT8uouuJx9cvjek33r4+v0/ZQ171qS0/7edJGMcFJe3mtl/SbvRlJv8X4stKzyUTlx5KulHSVpLEllIg2XANlhKc3HyP+tYN5eDrFs9Dya0+vUb7H24FvS/pmPv9MYFyY+iOmOVEIygCKmcewXvBlKeuGGR5jaTPPEkRK7sWDCDfiISQqKLa/ko8hfzCZQ9VbbY8KUY/tn0p6BrDU9icl7ULmMTaSqPIawrQv8TD//kHbDVsXVSBBkRSfhaTDoDREwm0oAOuAO4B7gX0lPXPUxpLeCZwInJwv2hL47LiThKCMYNKR4aaJVDxq22hrCXrER0YwLZJeB3wLWAW8O/9915jd/hB4EXlEFNu3Mjgm2EaEoJSkyVDlbf9y6wJxDzYl0qSPdN2GTwCeBvzU9rPJ4nL9fMw+D+buwwaQ9Mgx2wPRhlKaMiWJQduW6Q09qLF+0HGDIEgUN+bBVYZf2/61JCRtbfuHkvYZs885kj4G7CDp9WRDl5w57kRRQqmAYtXYMHGY9Di9+TIlohCe+dD0F3nT56+b1j/H6ZZQ1uQdG78EfFXSl4FbR+1g+++Ac8nGvtoHeIftD407UZRQKqSOUeDKHKfrGU0Q9zh1RJKxvACw/Yf533dJ+jrwaOAro/aR9F7bJ1IYP76wbChRQmmAcZ0fg3SJjD0YSoIlFEmbSbr6YRPtb9peafvBMbs+d8CysX1TQlDmSHjqBNMQItaCD62SLsNlSzGSDpf0I0mrJZ00YrunSVonaWD0YNvrgSuKPeXHnPcNkq4C9sl7yPemm4DoKd8k0/Z8bjqseTCa6NiXEWnQR0WN8pI2Bz5CVkpYA1wqaaXtawds914yN+BRLAGukXQJhYERbb9owLafBy4A/gYoCtl9tu8eZ3sIyhwYVb3V5DgTwfSEqMyHNqVxhW0oBwKrbd8IIOlsYDlwbd92byJrNB8YxDH36HqArO9JKWzfS9b58eW5YO1KphPbSdrO9s2j9g9BqYFJXoJhpZgopaRPuHdXx7DhIlpFeUHZWdJlhfkVtlcU5ncHbinMrwEOKh5A0u5knQ+fwxBBAb4LPAV4ne1XlbYuO/4byTo/3s6GspeBkWNghaAkQutenuBh6iytLMJHRSeucbIG9zttLxuxftBIXf1H/wBwou11I0bm3UrSUcDTJb1kkwPa542w4S3APrbvGrHNJoSgVEyVGUuUUoKUqUJIyz7fbXgPKqzyWgPsWZjfg037jSwDzs7FZGfg+ZLW2v5SYZvjgFcAO7DpgIgGRgnKLWRVXxPRmKDk9XOXAT+z/QJJOwFfAPYCfgK8LB8IBkknA8eQBTh7s+1xjVBzp87xMtrwMgXBwlOdoFwKLJW0N/Az4EjgTzY6lb1377+kTwH/0icm2P42Wdj6y2x/fEIbbgS+IelfyYZ97x3ztFE7Nek2fAJwXWH+JOAi20uBi/J5JO1LlqD7AYcDH83FKAlSdAVOzZ4gmJQ2PsNVDbBley3wRjLvreuAc2xfI+k4ScdNatcUYgJwM1mnxq3IgkL2ppE0UkKRtAfwB8ApwNvyxcuBZ+X/zwK+QRY+eTlwdu6tcJOk1WReEN+do8kDmddDP2kpJUo07aZN9y+83XIq7rRo+3zg/L5lZwzZ9jXVnfnhY74bsqCQtn81bvseTZVQPgD8JRt7bu9q+zaA/Lc35vEgj4fdBx1U0rGSLpN02c/vWle50U0SL233mXS4hC5TRxijOtEEUxuQdIika8lrkSTtL+mj4/abu6BIegFwh+3vl91lwLKB3wK2V9heZnvZLo9JplasMkJUuksbMs2q6VycugRDr/SQtLukp0t6Zm8as8sHgMOAuwBsX0E2dPBImqjy+l3gRZKeD2wDbC/ps8DtkpbYvk3SErLRxaCcx0MQNEYI/eR0TkxINzikpPcC/52sY2Sv6sZkg24NxfYtfS7JY6t95l5CsX2y7T1s70XW2P5vtl8JrASOyjc7Cvhy/n8lcKSkrXOvh6XAJXM2OwiCIbQt46+NdEsoLybrU/J82y/Mp0FhV4rcIunpgCVtJenP2diJaiApBYc8FXiupBvIYticCmD7GuAcMnX9CnC87SQaSOJFCoIAeHiArSq8vGrgRrIx4SfhOOB4svbqNcCT8/mRNNqx0fY3yLy5yHtkHjpku1PIPMKCoHN04cNkEm+vLlzvQBKt8gLuBy6XdBEb9yl587AdbN9J1ilyIqKnfBA0SGcz1wUk1TYUsmaDlZPsIOks4ATb9+TzOwLvt/3aUfuFoMyJCCQY9NM1MSlTSunaNW9EooJi+6wpdntST0zyY/xC0gHjdgpBqZFBL88swtLplzHoPF1/flMtoUhaSja+yb5knrUA2H7ciN02k7RjIfzVTpTQi5Qa5TvFuJent77rL1kwmEW67wvRYdNk3bTLTPPnk8DpwFrg2cCngc+M2ef9wHck/ZWk9wDfAf523ImihFIR07ww0ZDZfiYpaS7CPVyEaxyESLeEAjzC9kWSZPunwLsk/V/gncN2sP3pfMyW55Bd3kv6R4wcRAjKjMz6As2zbSUiFwdBjaQrKL+WtBlwQz5w1s/YENpqKLmAjBWRIlHltUBEIL9mibTvNrJLTQ3wFmBb4M3AU4FXsqETeaVECSURhmX2dZUoorQyO5MKRKR3h2kwTtc4bF8KkNV4+eg6zxWCkhDzynBCTKqhbIkv0noxSLUNRdIhwMeB7YDflLQ/8Ke2/6zqc0WVVxDUSIjJ4pBw6JUPMEXk4GmIEsqCERlcdUSbSLARiZZQYLrIwdMQgpIovcyqKgGo+njBaCKdFwynW+VFX+Rgssb5sZGDpyGqvBaIyOSqZVR6RullAUk3fP1UkYOnIUooiVJ1ySSYHyHci0fKHRunjRw8DVFC6TBFMYlMLgjqRetdapq7XdJZknYozO8o6RN1nCtKKB0lxKQZIq0XlIT7oTBl5OBpCEEJgimIqsSgn4ZcgsswVeTgaQhB6SjxpVwvEcYm2IR0Syi9yMHn5vMvpaYRcKMNJQimYJiYhMgsLnK5ad7Y/jTwR8Dt+fQS2+PC109FlFCCoEKiZLigGGgm8GNZtiRzRnP+vxaihBIEQaeZV6kx1dArkk4APgfsTBa2/rOS3lTHuUJQgmAKRg3vHKTFPO5Lrx9KilVewDHAQbbfafsdwMHA6+s4UVR5BcEU9H/1hpgsOHbKVV5i49hd6/JllRMllCCYgn4BGVWtEg31i0HCJZRPAhdLepekdwHfIwtnXzlRQgmCChhVQonSy4KQaAHF9mmSvgE8g6xkcrTt/6jjXCEoQVABMWhZkGosLwDbPwB+UPd5osorCKZkkmqvoOMYWOdyU4cJQQmCGQhRCXok3IYyN6LKKwhmJKq6AqBSLy9JhwMfBDYHzrR9at/6VwAn5rP/D3hDPrRvo4SgBEEQVEBVpQ9JmwMfAZ5LNiDWpZJW2r62sNlNwO/lkYOPAFYAB1VjwfRElVcQBMGslB2tsZzoHAistn2j7QeBs4HlG53O/k4vejCZG/Aes1/E7MxdUCTtKenrkq6TdE0eFgBJO0n6qqQb8t8dC/ucLGm1pB9JOmzeNgdBEIxCgNa51ATsLOmywnRs3+F2B24pzK/Jlw3jGOCCSi9oSpqo8loL/A/bP5D0KOD7kr4KvAa4yPapkk4CTgJOlLQvcCSwH7Ab8DVJT7C9bsjxgyAI5o7Kt6HcaXvZqEMNWDbw4JKeTSYozyh78jqZewnF9m25TzS27wOuI1Pf5cBZ+WZnAS/O/y8Hzrb9gO2bgNVkRcIgCII0qLbKaw2wZ2F+D+DW/o0kPQk4E1hu+67pja+ORttQJO0FHABcDOxq+zbIRIcsKiZMUPyTdGyvGPnzu6IAEwTBvPCGeF7jpvFcCiyVtLekrchqaFYWN5D0m8B5wKtsX1/55UxJY15ekrYDvgi8xfYvpaGxykoX/2yvIPN2YNn+23Tc4zsIgpSoysvL9lpJbwRWkbkNf8L2NZKOy9efAbwDeAzw0TzvXDumGm0uNCIokrYkE5PP2T4vX3y7pCW2b5O0BLgjX16q+BcEQdAoFfZDsX0+cH7fsjMK/18HvK6yE1ZEE15eIot0eZ3t0wqrVgJH5f+PAr5cWH6kpK0l7Q0sBS6Zl71BEARj8UReXp2liRLK7wKvAq6SdHm+7H8CpwLnSDoGuBl4KUBe1DsHuJbMQ+z48PAKgiA5uq0VpZi7oNj+NsMHdzl0yD6nAKfUZlQQBMGMTOA23Fki9EoQBEEVhKBE6JUgmAcRhbjjGFhfcuowUUIJgpoJMek+wlHlRQhKENRCbwTHGMlxgVjf8eJHCaLKi/iCDOohxGSBiCovIATlYYqiEgITVEU8S4uD7FJTl4kqrwLDRKX4lTls+aTHn2b/oB2EiNRH0u9Qx8WiDAshKINEYJKXvrftsPHDx40rPuyhT/rlCCZmUBVX3NPqSFuoSwd+7DRyRxNh2f7beMerXti0GTMRmVH7iI+EehgnJrOm8+ZLbvj+LMEVH/2IJT7kca8tte2qa/96pnOlzEKUUNpMNOwGi8o0tQj9zPPd6Xr7SBk6KyjXX7ktBw2NiN8Oei/JsKq1IOgiVVZtlXl3sm1umP1kISjh5dUmDttt/4enIB3ifkzHoHSrKy2HHbey8xlY73JTh+lsCaXrzOJtFlRHuJtPRn+JYZ5VuoPOXR3RKA8hKJ0g2lnmT4hHxqB0GPQs9m/XpBDXdr4QlBCUrhCiMl96YVV6jKpSaeN9mSXTXUixNbCu493gSxCC0iHCZTWYhIXM+GvD4BCUaJTvMJFh1Efb07bt9ieJXW7qMFFC6ThtrXKpmqpLb/1VXsGC0/PyWnBCUDpOiMlgqhDatolKm2xtJR0vfZQhBCVYWKbpXd3WTLmtdreKEJQQlK4SJZONmaQ0MWsw0ZRoq92tw4Z165q2onFCUDpIV8SkTGZY9lqnzVgjQw5KEyWUEJSu0WYxmSbzjgx/OJE2cyYEJQSlS7RNTCLDq4dI1ybofpyuMoSgBHNjUTO6JuJVBXPG4OjYGIIS1E9kcvMZgiDSuWEi9EoISlAvkcltTNWllUjfRLBhfQhKhF7pEE1mLvMc26LtVDWmTaRvYkTolSihdI3URCUYTqRXt3CUUEJQgiAIZqf7pY8yhKAEQRDMSgSHBEJQgiAIZsaAI/RKexrlJR0u6UeSVks6qWl7giAIHsb5AFtlphKMy++U8aF8/ZWSnlL5NU1BKwRF0ubAR4AjgH2Bl0vat1mrgiAINuD1LjWNo2R+dwSwNJ+OBU6v9mqmoxWCAhwIrLZ9o+0HgbOB5Q3bFARBsIHqSihl8rvlwKed8T1gB0lLqr2gyWlLG8ruwC2F+TXAQf0bSTqWTK0BHviaz716DrZNws7AnU0b0UfYVJ4U7QqbyjHOpt+a5eD38YtVX/O5O5fcfBtJlxXmV9heUZgvk98N2mZ34LaSNtRCWwRFA5ZtUnbMb8oKAEmX2V5Wt2GTEDaVI0WbIE27wqZy1G2T7cMrPFyZ/K5Unjhv2lLltQbYszC/B3BrQ7YEQRDUSZn8Lsk8sS2CcimwVNLekrYCjgRWNmxTEARBHZTJ71YCr869vQ4G7rXdaHUXtKTKy/ZaSW8EVgGbA5+wfc2Y3VaMWd8EYVM5UrQJ0rQrbCpHijYNZFh+J+m4fP0ZwPnA84HVwP3A0U3ZW0SOcAFBEARBBbSlyisIgiBInBCUIAiCoBI6JyhNhWiRtKekr0u6TtI1kk7Il+8k6auSbsh/dyzsc3Ju548kHVajbZtL+g9J/5KQTTtIOlfSD/M0O6RpuyS9Nb93V0v6R0nbzNsmSZ+QdIekqwvLJrZB0lMlXZWv+5CkQW6ms9j0vvzeXSnpnyTtME+bhtlVWPfnkixp58Kyudi10NjuzETWgPVj4HHAVsAVwL5zOvcS4Cn5/0cB15OFTfhb4KR8+UnAe/P/++b2bQ3sndu9eU22vQ34PPAv+XwKNp0FvC7/vxWwQ5N2kXUKuwl4RD5/DvCaedsEPBN4CnB1YdnENgCXAIeQ9Ve4ADiiYpueB2yR/3/vvG0aZle+fE+yBu2fAjvP265FnrpWQmksRIvt22z/IP9/H3AdWSa1nCzzJP99cf5/OXC27Qds30TmrXFg1XZJ2gP4A+DMwuKmbdqeLDP4OIDtB23f07RdZF6Pj5C0BbAtmV//XG2y/S3g7r7FE9mgLATH9ra/6yzH/HRhn0pssn2h7bX57PfI+kHMzaZhduX8PfCXbNzRb252LTJdE5Rh4QjmiqS9gAOAi4FdnfuH57+PzTebl60fIHu5ikGEmrbpccDPgU/mVXFnSnpkk3bZ/hnwd8DNZOEr7rV9YZM2FZjUht3z//OwDeC1ZF/2jdsk6UXAz2xf0bcqlbTqNF0TlMbDEUjaDvgi8Bbbvxy16YBlldoq6QXAHba/X3aXAcvqSL8tyKoqTrd9APArsqqcxuzK2yWWk1WH7AY8UtIrm7SpBMNsmJttkt4OrAU+17RNkrYF3g68Y9DqpuxaJLomKI2GI5C0JZmYfM72efni2/NiNfnvHXO09XeBF0n6CVn133MkfbZhm3rnWWP74nz+XDKBadKu3wdusv1z2w8B5wFPb9imHpPasIYNVVC12SbpKOAFwCvy6qKmbXo82QfBFfkzvwfwA0m/0bBdC0PXBKWxEC25Z8jHgetsn1ZYtRI4Kv9/FPDlwvIjJW0taW+ycQ0uqdIm2yfb3sP2XmRp8W+2X9mkTbld/wncImmffNGhwLUN23UzcLCkbfN7eShZO1ijaVU4V2kb8mqx+yQdnF/Lqwv7VIKkw4ETgRfZvr/P1kZssn2V7cfa3it/5teQOcr8Z5N2LRRNewVUPZGFI7iezIvj7XM87zPIispXApfn0/OBxwAXATfkvzsV9nl7buePqNmzBHgWG7y8GrcJeDJwWZ5eXwJ2bNou4N3AD4Grgc+QeQTN1SbgH8nacB4iyxCPmcYGYFl+HT8GPkweFaNCm1aTtUn0nvUz5mnTMLv61v+E3MtrnnYt8hShV4IgCIJK6FqVVxAEQdAQIShBEARBJYSgBEEQBJUQghIEQRBUQghKEARBUAkhKEEQBEElhKAEQRAElRCCEnQSSXspG2flfysb4+RCSY+SdKmkZ+Xb/I2kU5q1NAi6QwhK0GWWAh+xvR9wD1nwx9cAp0t6LnA4We/4IAgqIAQl6DI32b48//99YC/b15CFVfln4LXOxs0JgqACQlCCLvNA4f86srD5AL9DVmLZdd4GBUGXCUEJFgpJLyELtvhM4EPFsdCDIJiNEJRgkdgZOJUsKu31ZJFlP9isSUHQHSLacBAEQVAJUUIJgiAIKiEEJQiCIKiEEJQgCIKgEkJQgiAIgkoIQQmCIAgqIQQlCIIgqIQQlCAIgqAS/j99s5k/MulhfAAAAABJRU5ErkJggg==\n", | |
| "text/plain": [ | |
| "<Figure size 432x288 with 2 Axes>" | |
| ] | |
| }, | |
| "metadata": { | |
| "needs_background": "light" | |
| }, | |
| "output_type": "display_data" | |
| } | |
| ], | |
| "source": [ | |
| "kmt_orig['kmt'].plot()" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "d633d24c-882a-4cfe-90c6-4c8777024d13", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/html": [ | |
| "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", | |
| "<defs>\n", | |
| "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", | |
| "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "</symbol>\n", | |
| "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", | |
| "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "</symbol>\n", | |
| "</defs>\n", | |
| "</svg>\n", | |
| "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", | |
| " *\n", | |
| " */\n", | |
| "\n", | |
| ":root {\n", | |
| " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", | |
| " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", | |
| " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", | |
| " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", | |
| " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", | |
| " --xr-background-color: var(--jp-layout-color0, white);\n", | |
| " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", | |
| " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", | |
| "}\n", | |
| "\n", | |
| "html[theme=dark],\n", | |
| "body.vscode-dark {\n", | |
| " --xr-font-color0: rgba(255, 255, 255, 1);\n", | |
| " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", | |
| " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", | |
| " --xr-border-color: #1F1F1F;\n", | |
| " --xr-disabled-color: #515151;\n", | |
| " --xr-background-color: #111111;\n", | |
| " --xr-background-color-row-even: #111111;\n", | |
| " --xr-background-color-row-odd: #313131;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-wrap {\n", | |
| " display: block;\n", | |
| " min-width: 300px;\n", | |
| " max-width: 700px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-text-repr-fallback {\n", | |
| " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header {\n", | |
| " padding-top: 6px;\n", | |
| " padding-bottom: 6px;\n", | |
| " margin-bottom: 4px;\n", | |
| " border-bottom: solid 1px var(--xr-border-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header > div,\n", | |
| ".xr-header > ul {\n", | |
| " display: inline;\n", | |
| " margin-top: 0;\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type,\n", | |
| ".xr-array-name {\n", | |
| " margin-left: 2px;\n", | |
| " margin-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-sections {\n", | |
| " padding-left: 0 !important;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input + label {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label {\n", | |
| " cursor: pointer;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label:hover {\n", | |
| " color: var(--xr-font-color0);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary {\n", | |
| " grid-column: 1;\n", | |
| " color: var(--xr-font-color2);\n", | |
| " font-weight: 500;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary > span {\n", | |
| " display: inline-block;\n", | |
| " padding-left: 0.5em;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in + label:before {\n", | |
| " display: inline-block;\n", | |
| " content: '►';\n", | |
| " font-size: 11px;\n", | |
| " width: 15px;\n", | |
| " text-align: center;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label:before {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label:before {\n", | |
| " content: '▼';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label > span {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary,\n", | |
| ".xr-section-inline-details {\n", | |
| " padding-top: 4px;\n", | |
| " padding-bottom: 4px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-inline-details {\n", | |
| " grid-column: 2 / -1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-details {\n", | |
| " display: none;\n", | |
| " grid-column: 1 / -1;\n", | |
| " margin-bottom: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked ~ .xr-section-details {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap {\n", | |
| " grid-column: 1 / -1;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 20px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap > label {\n", | |
| " grid-column: 1;\n", | |
| " vertical-align: top;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-preview {\n", | |
| " color: var(--xr-font-color3);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-preview,\n", | |
| ".xr-array-data {\n", | |
| " padding: 0 5px !important;\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-data,\n", | |
| ".xr-array-in:checked ~ .xr-array-preview {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-in:checked ~ .xr-array-data,\n", | |
| ".xr-array-preview {\n", | |
| " display: inline-block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list {\n", | |
| " display: inline-block !important;\n", | |
| " list-style: none;\n", | |
| " padding: 0 !important;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li {\n", | |
| " display: inline-block;\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:before {\n", | |
| " content: '(';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:after {\n", | |
| " content: ')';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li:not(:last-child):after {\n", | |
| " content: ',';\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-has-index {\n", | |
| " font-weight: bold;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list,\n", | |
| ".xr-var-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > div,\n", | |
| ".xr-var-item label,\n", | |
| ".xr-var-item > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-even);\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > .xr-var-name:hover span {\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list > li:nth-child(odd) > div,\n", | |
| ".xr-var-list > li:nth-child(odd) > label,\n", | |
| ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-odd);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name {\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dims {\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dtype {\n", | |
| " grid-column: 3;\n", | |
| " text-align: right;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-preview {\n", | |
| " grid-column: 4;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name,\n", | |
| ".xr-var-dims,\n", | |
| ".xr-var-dtype,\n", | |
| ".xr-preview,\n", | |
| ".xr-attrs dt {\n", | |
| " white-space: nowrap;\n", | |
| " overflow: hidden;\n", | |
| " text-overflow: ellipsis;\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name:hover,\n", | |
| ".xr-var-dims:hover,\n", | |
| ".xr-var-dtype:hover,\n", | |
| ".xr-attrs dt:hover {\n", | |
| " overflow: visible;\n", | |
| " width: auto;\n", | |
| " z-index: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " display: none;\n", | |
| " background-color: var(--xr-background-color) !important;\n", | |
| " padding-bottom: 5px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", | |
| ".xr-var-data-in:checked ~ .xr-var-data {\n", | |
| " display: block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-data > table {\n", | |
| " float: right;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name span,\n", | |
| ".xr-var-data,\n", | |
| ".xr-attrs {\n", | |
| " padding-left: 25px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs,\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " grid-column: 1 / -1;\n", | |
| "}\n", | |
| "\n", | |
| "dl.xr-attrs {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 125px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt,\n", | |
| ".xr-attrs dd {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " float: left;\n", | |
| " padding-right: 10px;\n", | |
| " width: auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt {\n", | |
| " font-weight: normal;\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt:hover span {\n", | |
| " display: inline-block;\n", | |
| " background: var(--xr-background-color);\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dd {\n", | |
| " grid-column: 2;\n", | |
| " white-space: pre-wrap;\n", | |
| " word-break: break-all;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-icon-database,\n", | |
| ".xr-icon-file-text2 {\n", | |
| " display: inline-block;\n", | |
| " vertical-align: middle;\n", | |
| " width: 1em;\n", | |
| " height: 1.5em !important;\n", | |
| " stroke-width: 0;\n", | |
| " stroke: currentColor;\n", | |
| " fill: currentColor;\n", | |
| "}\n", | |
| "</style><pre class='xr-text-repr-fallback'><xarray.DataArray (ny: 1080, nx: 1440)>\n", | |
| "array([[0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " ...,\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.]])\n", | |
| "Dimensions without coordinates: ny, nx</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span>ny</span>: 1080</li><li><span>nx</span>: 1440</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-0ac77133-9579-49fe-b797-94cc1319c9de' class='xr-array-in' type='checkbox' checked><label for='section-0ac77133-9579-49fe-b797-94cc1319c9de' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 ... 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0</span></div><div class='xr-array-data'><pre>array([[0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " ...,\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.]])</pre></div></div></li><li class='xr-section-item'><input id='section-6c99df38-4b6b-4f19-8ba2-f11c890c1b1a' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-6c99df38-4b6b-4f19-8ba2-f11c890c1b1a' class='xr-section-summary' title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-2a8a46eb-eb71-4f7d-a06e-778193d3d407' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-2a8a46eb-eb71-4f7d-a06e-778193d3d407' class='xr-section-summary' title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>" | |
| ], | |
| "text/plain": [ | |
| "<xarray.DataArray (ny: 1080, nx: 1440)>\n", | |
| "array([[0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " ...,\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.],\n", | |
| " [0., 0., 0., ..., 0., 0., 0.]])\n", | |
| "Dimensions without coordinates: ny, nx" | |
| ] | |
| }, | |
| "execution_count": 3, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| } | |
| ], | |
| "source": [ | |
| "mask_orig-kmt_orig" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 4, | |
| "id": "01aeb457-7787-47e3-ac61-e61130b70095", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [ | |
| "diff = (mask_orig-kmt_orig).stack(xy=('nx','ny'))" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 5, | |
| "id": "be50a63f-d769-440a-b21a-2eef827a2192", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/html": [ | |
| "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", | |
| "<defs>\n", | |
| "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", | |
| "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "</symbol>\n", | |
| "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", | |
| "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "</symbol>\n", | |
| "</defs>\n", | |
| "</svg>\n", | |
| "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", | |
| " *\n", | |
| " */\n", | |
| "\n", | |
| ":root {\n", | |
| " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", | |
| " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", | |
| " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", | |
| " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", | |
| " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", | |
| " --xr-background-color: var(--jp-layout-color0, white);\n", | |
| " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", | |
| " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", | |
| "}\n", | |
| "\n", | |
| "html[theme=dark],\n", | |
| "body.vscode-dark {\n", | |
| " --xr-font-color0: rgba(255, 255, 255, 1);\n", | |
| " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", | |
| " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", | |
| " --xr-border-color: #1F1F1F;\n", | |
| " --xr-disabled-color: #515151;\n", | |
| " --xr-background-color: #111111;\n", | |
| " --xr-background-color-row-even: #111111;\n", | |
| " --xr-background-color-row-odd: #313131;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-wrap {\n", | |
| " display: block;\n", | |
| " min-width: 300px;\n", | |
| " max-width: 700px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-text-repr-fallback {\n", | |
| " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header {\n", | |
| " padding-top: 6px;\n", | |
| " padding-bottom: 6px;\n", | |
| " margin-bottom: 4px;\n", | |
| " border-bottom: solid 1px var(--xr-border-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header > div,\n", | |
| ".xr-header > ul {\n", | |
| " display: inline;\n", | |
| " margin-top: 0;\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type,\n", | |
| ".xr-array-name {\n", | |
| " margin-left: 2px;\n", | |
| " margin-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-sections {\n", | |
| " padding-left: 0 !important;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input + label {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label {\n", | |
| " cursor: pointer;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label:hover {\n", | |
| " color: var(--xr-font-color0);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary {\n", | |
| " grid-column: 1;\n", | |
| " color: var(--xr-font-color2);\n", | |
| " font-weight: 500;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary > span {\n", | |
| " display: inline-block;\n", | |
| " padding-left: 0.5em;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in + label:before {\n", | |
| " display: inline-block;\n", | |
| " content: '►';\n", | |
| " font-size: 11px;\n", | |
| " width: 15px;\n", | |
| " text-align: center;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label:before {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label:before {\n", | |
| " content: '▼';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label > span {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary,\n", | |
| ".xr-section-inline-details {\n", | |
| " padding-top: 4px;\n", | |
| " padding-bottom: 4px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-inline-details {\n", | |
| " grid-column: 2 / -1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-details {\n", | |
| " display: none;\n", | |
| " grid-column: 1 / -1;\n", | |
| " margin-bottom: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked ~ .xr-section-details {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap {\n", | |
| " grid-column: 1 / -1;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 20px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap > label {\n", | |
| " grid-column: 1;\n", | |
| " vertical-align: top;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-preview {\n", | |
| " color: var(--xr-font-color3);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-preview,\n", | |
| ".xr-array-data {\n", | |
| " padding: 0 5px !important;\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-data,\n", | |
| ".xr-array-in:checked ~ .xr-array-preview {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-in:checked ~ .xr-array-data,\n", | |
| ".xr-array-preview {\n", | |
| " display: inline-block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list {\n", | |
| " display: inline-block !important;\n", | |
| " list-style: none;\n", | |
| " padding: 0 !important;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li {\n", | |
| " display: inline-block;\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:before {\n", | |
| " content: '(';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:after {\n", | |
| " content: ')';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li:not(:last-child):after {\n", | |
| " content: ',';\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-has-index {\n", | |
| " font-weight: bold;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list,\n", | |
| ".xr-var-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > div,\n", | |
| ".xr-var-item label,\n", | |
| ".xr-var-item > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-even);\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > .xr-var-name:hover span {\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list > li:nth-child(odd) > div,\n", | |
| ".xr-var-list > li:nth-child(odd) > label,\n", | |
| ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-odd);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name {\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dims {\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dtype {\n", | |
| " grid-column: 3;\n", | |
| " text-align: right;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-preview {\n", | |
| " grid-column: 4;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name,\n", | |
| ".xr-var-dims,\n", | |
| ".xr-var-dtype,\n", | |
| ".xr-preview,\n", | |
| ".xr-attrs dt {\n", | |
| " white-space: nowrap;\n", | |
| " overflow: hidden;\n", | |
| " text-overflow: ellipsis;\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name:hover,\n", | |
| ".xr-var-dims:hover,\n", | |
| ".xr-var-dtype:hover,\n", | |
| ".xr-attrs dt:hover {\n", | |
| " overflow: visible;\n", | |
| " width: auto;\n", | |
| " z-index: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " display: none;\n", | |
| " background-color: var(--xr-background-color) !important;\n", | |
| " padding-bottom: 5px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", | |
| ".xr-var-data-in:checked ~ .xr-var-data {\n", | |
| " display: block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-data > table {\n", | |
| " float: right;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name span,\n", | |
| ".xr-var-data,\n", | |
| ".xr-attrs {\n", | |
| " padding-left: 25px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs,\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " grid-column: 1 / -1;\n", | |
| "}\n", | |
| "\n", | |
| "dl.xr-attrs {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 125px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt,\n", | |
| ".xr-attrs dd {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " float: left;\n", | |
| " padding-right: 10px;\n", | |
| " width: auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt {\n", | |
| " font-weight: normal;\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt:hover span {\n", | |
| " display: inline-block;\n", | |
| " background: var(--xr-background-color);\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dd {\n", | |
| " grid-column: 2;\n", | |
| " white-space: pre-wrap;\n", | |
| " word-break: break-all;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-icon-database,\n", | |
| ".xr-icon-file-text2 {\n", | |
| " display: inline-block;\n", | |
| " vertical-align: middle;\n", | |
| " width: 1em;\n", | |
| " height: 1.5em !important;\n", | |
| " stroke-width: 0;\n", | |
| " stroke: currentColor;\n", | |
| " fill: currentColor;\n", | |
| "}\n", | |
| "</style><pre class='xr-text-repr-fallback'><xarray.DataArray (xy: 0)>\n", | |
| "array([], dtype=float64)\n", | |
| "Coordinates:\n", | |
| " * xy (xy) MultiIndex\n", | |
| " - nx (xy) int64 \n", | |
| " - ny (xy) int64 </pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'></div><ul class='xr-dim-list'><li><span class='xr-has-index'>xy</span>: 0</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-21dbc41c-4c99-42b7-b86d-13f36dfa49b4' class='xr-array-in' type='checkbox' checked><label for='section-21dbc41c-4c99-42b7-b86d-13f36dfa49b4' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span></span></div><div class='xr-array-data'><pre>array([], dtype=float64)</pre></div></div></li><li class='xr-section-item'><input id='section-58349c49-5598-4e08-a30d-753bd2f6702c' class='xr-section-summary-in' type='checkbox' checked><label for='section-58349c49-5598-4e08-a30d-753bd2f6702c' class='xr-section-summary' >Coordinates: <span>(1)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>xy</span></div><div class='xr-var-dims'>(xy)</div><div class='xr-var-dtype'>MultiIndex</div><div class='xr-var-preview xr-preview'>(nx, ny)</div><input id='attrs-a4b7ba95-4721-4694-8e4d-a63403562990' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a4b7ba95-4721-4694-8e4d-a63403562990' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5d36ccb8-898a-4433-9893-cc31371bcf5b' class='xr-var-data-in' type='checkbox'><label for='data-5d36ccb8-898a-4433-9893-cc31371bcf5b' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([], dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>nx</span></div><div class='xr-var-dims'>(xy)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'></div><input id='attrs-ffd2ad81-513b-4cd3-bda8-c6693104a001' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ffd2ad81-513b-4cd3-bda8-c6693104a001' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e379b66c-4073-4659-94c0-faac349e00c1' class='xr-var-data-in' type='checkbox'><label for='data-e379b66c-4073-4659-94c0-faac349e00c1' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([], dtype=int64)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ny</span></div><div class='xr-var-dims'>(xy)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'></div><input id='attrs-2328bc3a-3e9d-4ec0-91bc-586a99bb7305' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-2328bc3a-3e9d-4ec0-91bc-586a99bb7305' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6cad1f42-065e-498c-aed1-1288051e09d2' class='xr-var-data-in' type='checkbox'><label for='data-6cad1f42-065e-498c-aed1-1288051e09d2' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([], dtype=int64)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-ec70c78e-8636-48e1-960d-33dfe9b475aa' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-ec70c78e-8636-48e1-960d-33dfe9b475aa' class='xr-section-summary' title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>" | |
| ], | |
| "text/plain": [ | |
| "<xarray.DataArray (xy: 0)>\n", | |
| "array([], dtype=float64)\n", | |
| "Coordinates:\n", | |
| " * xy (xy) MultiIndex\n", | |
| " - nx (xy) int64 \n", | |
| " - ny (xy) int64 " | |
| ] | |
| }, | |
| "execution_count": 5, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| } | |
| ], | |
| "source": [ | |
| "diff.where(diff != 0, drop=True)" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 10, | |
| "id": "e695913d-e933-4d90-ba7d-1c880f296a98", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/html": [ | |
| "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", | |
| "<defs>\n", | |
| "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", | |
| "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "</symbol>\n", | |
| "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", | |
| "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "</symbol>\n", | |
| "</defs>\n", | |
| "</svg>\n", | |
| "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", | |
| " *\n", | |
| " */\n", | |
| "\n", | |
| ":root {\n", | |
| " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", | |
| " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", | |
| " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", | |
| " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", | |
| " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", | |
| " --xr-background-color: var(--jp-layout-color0, white);\n", | |
| " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", | |
| " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", | |
| "}\n", | |
| "\n", | |
| "html[theme=dark],\n", | |
| "body.vscode-dark {\n", | |
| " --xr-font-color0: rgba(255, 255, 255, 1);\n", | |
| " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", | |
| " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", | |
| " --xr-border-color: #1F1F1F;\n", | |
| " --xr-disabled-color: #515151;\n", | |
| " --xr-background-color: #111111;\n", | |
| " --xr-background-color-row-even: #111111;\n", | |
| " --xr-background-color-row-odd: #313131;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-wrap {\n", | |
| " display: block;\n", | |
| " min-width: 300px;\n", | |
| " max-width: 700px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-text-repr-fallback {\n", | |
| " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header {\n", | |
| " padding-top: 6px;\n", | |
| " padding-bottom: 6px;\n", | |
| " margin-bottom: 4px;\n", | |
| " border-bottom: solid 1px var(--xr-border-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header > div,\n", | |
| ".xr-header > ul {\n", | |
| " display: inline;\n", | |
| " margin-top: 0;\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type,\n", | |
| ".xr-array-name {\n", | |
| " margin-left: 2px;\n", | |
| " margin-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-sections {\n", | |
| " padding-left: 0 !important;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input + label {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label {\n", | |
| " cursor: pointer;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label:hover {\n", | |
| " color: var(--xr-font-color0);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary {\n", | |
| " grid-column: 1;\n", | |
| " color: var(--xr-font-color2);\n", | |
| " font-weight: 500;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary > span {\n", | |
| " display: inline-block;\n", | |
| " padding-left: 0.5em;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in + label:before {\n", | |
| " display: inline-block;\n", | |
| " content: '►';\n", | |
| " font-size: 11px;\n", | |
| " width: 15px;\n", | |
| " text-align: center;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label:before {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label:before {\n", | |
| " content: '▼';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label > span {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary,\n", | |
| ".xr-section-inline-details {\n", | |
| " padding-top: 4px;\n", | |
| " padding-bottom: 4px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-inline-details {\n", | |
| " grid-column: 2 / -1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-details {\n", | |
| " display: none;\n", | |
| " grid-column: 1 / -1;\n", | |
| " margin-bottom: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked ~ .xr-section-details {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap {\n", | |
| " grid-column: 1 / -1;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 20px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap > label {\n", | |
| " grid-column: 1;\n", | |
| " vertical-align: top;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-preview {\n", | |
| " color: var(--xr-font-color3);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-preview,\n", | |
| ".xr-array-data {\n", | |
| " padding: 0 5px !important;\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-data,\n", | |
| ".xr-array-in:checked ~ .xr-array-preview {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-in:checked ~ .xr-array-data,\n", | |
| ".xr-array-preview {\n", | |
| " display: inline-block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list {\n", | |
| " display: inline-block !important;\n", | |
| " list-style: none;\n", | |
| " padding: 0 !important;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li {\n", | |
| " display: inline-block;\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:before {\n", | |
| " content: '(';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:after {\n", | |
| " content: ')';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li:not(:last-child):after {\n", | |
| " content: ',';\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-has-index {\n", | |
| " font-weight: bold;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list,\n", | |
| ".xr-var-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > div,\n", | |
| ".xr-var-item label,\n", | |
| ".xr-var-item > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-even);\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > .xr-var-name:hover span {\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list > li:nth-child(odd) > div,\n", | |
| ".xr-var-list > li:nth-child(odd) > label,\n", | |
| ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-odd);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name {\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dims {\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dtype {\n", | |
| " grid-column: 3;\n", | |
| " text-align: right;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-preview {\n", | |
| " grid-column: 4;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name,\n", | |
| ".xr-var-dims,\n", | |
| ".xr-var-dtype,\n", | |
| ".xr-preview,\n", | |
| ".xr-attrs dt {\n", | |
| " white-space: nowrap;\n", | |
| " overflow: hidden;\n", | |
| " text-overflow: ellipsis;\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name:hover,\n", | |
| ".xr-var-dims:hover,\n", | |
| ".xr-var-dtype:hover,\n", | |
| ".xr-attrs dt:hover {\n", | |
| " overflow: visible;\n", | |
| " width: auto;\n", | |
| " z-index: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " display: none;\n", | |
| " background-color: var(--xr-background-color) !important;\n", | |
| " padding-bottom: 5px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", | |
| ".xr-var-data-in:checked ~ .xr-var-data {\n", | |
| " display: block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-data > table {\n", | |
| " float: right;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name span,\n", | |
| ".xr-var-data,\n", | |
| ".xr-attrs {\n", | |
| " padding-left: 25px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs,\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " grid-column: 1 / -1;\n", | |
| "}\n", | |
| "\n", | |
| "dl.xr-attrs {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 125px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt,\n", | |
| ".xr-attrs dd {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " float: left;\n", | |
| " padding-right: 10px;\n", | |
| " width: auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt {\n", | |
| " font-weight: normal;\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt:hover span {\n", | |
| " display: inline-block;\n", | |
| " background: var(--xr-background-color);\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dd {\n", | |
| " grid-column: 2;\n", | |
| " white-space: pre-wrap;\n", | |
| " word-break: break-all;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-icon-database,\n", | |
| ".xr-icon-file-text2 {\n", | |
| " display: inline-block;\n", | |
| " vertical-align: middle;\n", | |
| " width: 1em;\n", | |
| " height: 1.5em !important;\n", | |
| " stroke-width: 0;\n", | |
| " stroke: currentColor;\n", | |
| " fill: currentColor;\n", | |
| "}\n", | |
| "</style><pre class='xr-text-repr-fallback'><xarray.DataArray 'mask' ()>\n", | |
| "array(0.)</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'mask'</div></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-f5b379a1-bfd6-401c-aa7c-861243d76561' class='xr-array-in' type='checkbox' checked><label for='section-f5b379a1-bfd6-401c-aa7c-861243d76561' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>0.0</span></div><div class='xr-array-data'><pre>array(0.)</pre></div></div></li><li class='xr-section-item'><input id='section-be8d9b13-4e82-48d0-a266-c4bae438eade' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-be8d9b13-4e82-48d0-a266-c4bae438eade' class='xr-section-summary' title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-8a823e4b-d6c6-41f2-985f-f86c5a385aa4' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-8a823e4b-d6c6-41f2-985f-f86c5a385aa4' class='xr-section-summary' title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>" | |
| ], | |
| "text/plain": [ | |
| "<xarray.DataArray 'mask' ()>\n", | |
| "array(0.)" | |
| ] | |
| }, | |
| "execution_count": 10, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| } | |
| ], | |
| "source": [ | |
| "mask_orig.where(np.isnan(topog), drop=True).sum(dim='nx').sum('ny')" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": 11, | |
| "id": "d4e947a7-c587-4e1a-a4c0-6747c846f65a", | |
| "metadata": {}, | |
| "outputs": [ | |
| { | |
| "data": { | |
| "text/html": [ | |
| "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n", | |
| "<defs>\n", | |
| "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n", | |
| "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n", | |
| "</symbol>\n", | |
| "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n", | |
| "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n", | |
| "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n", | |
| "</symbol>\n", | |
| "</defs>\n", | |
| "</svg>\n", | |
| "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n", | |
| " *\n", | |
| " */\n", | |
| "\n", | |
| ":root {\n", | |
| " --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n", | |
| " --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n", | |
| " --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n", | |
| " --xr-border-color: var(--jp-border-color2, #e0e0e0);\n", | |
| " --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n", | |
| " --xr-background-color: var(--jp-layout-color0, white);\n", | |
| " --xr-background-color-row-even: var(--jp-layout-color1, white);\n", | |
| " --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n", | |
| "}\n", | |
| "\n", | |
| "html[theme=dark],\n", | |
| "body.vscode-dark {\n", | |
| " --xr-font-color0: rgba(255, 255, 255, 1);\n", | |
| " --xr-font-color2: rgba(255, 255, 255, 0.54);\n", | |
| " --xr-font-color3: rgba(255, 255, 255, 0.38);\n", | |
| " --xr-border-color: #1F1F1F;\n", | |
| " --xr-disabled-color: #515151;\n", | |
| " --xr-background-color: #111111;\n", | |
| " --xr-background-color-row-even: #111111;\n", | |
| " --xr-background-color-row-odd: #313131;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-wrap {\n", | |
| " display: block;\n", | |
| " min-width: 300px;\n", | |
| " max-width: 700px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-text-repr-fallback {\n", | |
| " /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header {\n", | |
| " padding-top: 6px;\n", | |
| " padding-bottom: 6px;\n", | |
| " margin-bottom: 4px;\n", | |
| " border-bottom: solid 1px var(--xr-border-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-header > div,\n", | |
| ".xr-header > ul {\n", | |
| " display: inline;\n", | |
| " margin-top: 0;\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type,\n", | |
| ".xr-array-name {\n", | |
| " margin-left: 2px;\n", | |
| " margin-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-obj-type {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-sections {\n", | |
| " padding-left: 0 !important;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 150px auto auto 1fr 20px 20px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input + label {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label {\n", | |
| " cursor: pointer;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-item input:enabled + label:hover {\n", | |
| " color: var(--xr-font-color0);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary {\n", | |
| " grid-column: 1;\n", | |
| " color: var(--xr-font-color2);\n", | |
| " font-weight: 500;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary > span {\n", | |
| " display: inline-block;\n", | |
| " padding-left: 0.5em;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label {\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in + label:before {\n", | |
| " display: inline-block;\n", | |
| " content: '►';\n", | |
| " font-size: 11px;\n", | |
| " width: 15px;\n", | |
| " text-align: center;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:disabled + label:before {\n", | |
| " color: var(--xr-disabled-color);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label:before {\n", | |
| " content: '▼';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked + label > span {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary,\n", | |
| ".xr-section-inline-details {\n", | |
| " padding-top: 4px;\n", | |
| " padding-bottom: 4px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-inline-details {\n", | |
| " grid-column: 2 / -1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-details {\n", | |
| " display: none;\n", | |
| " grid-column: 1 / -1;\n", | |
| " margin-bottom: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-section-summary-in:checked ~ .xr-section-details {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap {\n", | |
| " grid-column: 1 / -1;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 20px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-wrap > label {\n", | |
| " grid-column: 1;\n", | |
| " vertical-align: top;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-preview {\n", | |
| " color: var(--xr-font-color3);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-preview,\n", | |
| ".xr-array-data {\n", | |
| " padding: 0 5px !important;\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-data,\n", | |
| ".xr-array-in:checked ~ .xr-array-preview {\n", | |
| " display: none;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-array-in:checked ~ .xr-array-data,\n", | |
| ".xr-array-preview {\n", | |
| " display: inline-block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list {\n", | |
| " display: inline-block !important;\n", | |
| " list-style: none;\n", | |
| " padding: 0 !important;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li {\n", | |
| " display: inline-block;\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:before {\n", | |
| " content: '(';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list:after {\n", | |
| " content: ')';\n", | |
| "}\n", | |
| "\n", | |
| ".xr-dim-list li:not(:last-child):after {\n", | |
| " content: ',';\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-has-index {\n", | |
| " font-weight: bold;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list,\n", | |
| ".xr-var-item {\n", | |
| " display: contents;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > div,\n", | |
| ".xr-var-item label,\n", | |
| ".xr-var-item > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-even);\n", | |
| " margin-bottom: 0;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-item > .xr-var-name:hover span {\n", | |
| " padding-right: 5px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-list > li:nth-child(odd) > div,\n", | |
| ".xr-var-list > li:nth-child(odd) > label,\n", | |
| ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n", | |
| " background-color: var(--xr-background-color-row-odd);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name {\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dims {\n", | |
| " grid-column: 2;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-dtype {\n", | |
| " grid-column: 3;\n", | |
| " text-align: right;\n", | |
| " color: var(--xr-font-color2);\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-preview {\n", | |
| " grid-column: 4;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name,\n", | |
| ".xr-var-dims,\n", | |
| ".xr-var-dtype,\n", | |
| ".xr-preview,\n", | |
| ".xr-attrs dt {\n", | |
| " white-space: nowrap;\n", | |
| " overflow: hidden;\n", | |
| " text-overflow: ellipsis;\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name:hover,\n", | |
| ".xr-var-dims:hover,\n", | |
| ".xr-var-dtype:hover,\n", | |
| ".xr-attrs dt:hover {\n", | |
| " overflow: visible;\n", | |
| " width: auto;\n", | |
| " z-index: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " display: none;\n", | |
| " background-color: var(--xr-background-color) !important;\n", | |
| " padding-bottom: 5px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n", | |
| ".xr-var-data-in:checked ~ .xr-var-data {\n", | |
| " display: block;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-data > table {\n", | |
| " float: right;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-var-name span,\n", | |
| ".xr-var-data,\n", | |
| ".xr-attrs {\n", | |
| " padding-left: 25px !important;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs,\n", | |
| ".xr-var-attrs,\n", | |
| ".xr-var-data {\n", | |
| " grid-column: 1 / -1;\n", | |
| "}\n", | |
| "\n", | |
| "dl.xr-attrs {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " display: grid;\n", | |
| " grid-template-columns: 125px auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt,\n", | |
| ".xr-attrs dd {\n", | |
| " padding: 0;\n", | |
| " margin: 0;\n", | |
| " float: left;\n", | |
| " padding-right: 10px;\n", | |
| " width: auto;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt {\n", | |
| " font-weight: normal;\n", | |
| " grid-column: 1;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dt:hover span {\n", | |
| " display: inline-block;\n", | |
| " background: var(--xr-background-color);\n", | |
| " padding-right: 10px;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-attrs dd {\n", | |
| " grid-column: 2;\n", | |
| " white-space: pre-wrap;\n", | |
| " word-break: break-all;\n", | |
| "}\n", | |
| "\n", | |
| ".xr-icon-database,\n", | |
| ".xr-icon-file-text2 {\n", | |
| " display: inline-block;\n", | |
| " vertical-align: middle;\n", | |
| " width: 1em;\n", | |
| " height: 1.5em !important;\n", | |
| " stroke-width: 0;\n", | |
| " stroke: currentColor;\n", | |
| " fill: currentColor;\n", | |
| "}\n", | |
| "</style><pre class='xr-text-repr-fallback'><xarray.DataArray 'kmt' ()>\n", | |
| "array(0.)</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'kmt'</div></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-267044de-6a54-4f95-a260-8f4c3b22bb58' class='xr-array-in' type='checkbox' checked><label for='section-267044de-6a54-4f95-a260-8f4c3b22bb58' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>0.0</span></div><div class='xr-array-data'><pre>array(0.)</pre></div></div></li><li class='xr-section-item'><input id='section-d221f49b-c9ed-438f-b0ef-5fcacf42b07d' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-d221f49b-c9ed-438f-b0ef-5fcacf42b07d' class='xr-section-summary' title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-c52dc229-24d3-4d84-8db4-3750d35b6452' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-c52dc229-24d3-4d84-8db4-3750d35b6452' class='xr-section-summary' title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>" | |
| ], | |
| "text/plain": [ | |
| "<xarray.DataArray 'kmt' ()>\n", | |
| "array(0.)" | |
| ] | |
| }, | |
| "execution_count": 11, | |
| "metadata": {}, | |
| "output_type": "execute_result" | |
| } | |
| ], | |
| "source": [ | |
| "kmt_orig.where(np.isnan(topog), drop=True).sum(dim='nx').sum('ny')" | |
| ] | |
| }, | |
| { | |
| "cell_type": "code", | |
| "execution_count": null, | |
| "id": "59198164-101f-4fa8-b443-4accfc117990", | |
| "metadata": {}, | |
| "outputs": [], | |
| "source": [] | |
| } | |
| ], | |
| "metadata": { | |
| "kernelspec": { | |
| "display_name": "Python [conda env:analysis3-21.07]", | |
| "language": "python", | |
| "name": "conda-env-analysis3-21.07-py" | |
| }, | |
| "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.7" | |
| } | |
| }, | |
| "nbformat": 4, | |
| "nbformat_minor": 5 | |
| } |
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment