tsl_calc_checkup.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "a23560a8-1945-48f8-bb7e-18add8188fbc",
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/niki/anaconda3/envs/cosipy_test/lib/python3.10/site-packages/dask_jobqueue/core.py:20: FutureWarning: tmpfile is deprecated and will be removed in a future release. Please use dask.utils.tmpfile instead.\n",
      "  from distributed.utils import tmpfile\n"
     ]
    }
   ],
   "source": [
    "from config import *\n",
    "from constants import *\n",
    "import xarray as xr\n",
    "from COSIPY import prereq_res, construct_resampled_ds\n",
    "from cosipy.modules.evaluation import * \n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import sys\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "0c18fb5a-1647-40a7-8e70-a1bc3aee75eb",
   "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'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:       (lat: 14, lon: 30, time: 96444)\n",
       "Coordinates:\n",
       "  * time          (time) datetime64[ns] 1999-01-01 1999-01-01 ... 2010-01-01\n",
       "  * lat           (lat) float64 46.79 46.79 46.79 46.79 ... 46.82 46.82 46.82\n",
       "  * lon           (lon) float64 10.72 10.73 10.73 10.73 ... 10.8 10.8 10.8 10.8\n",
       "Data variables: (12/37)\n",
       "    HGT           (lat, lon) float32 3.227e+03 3.244e+03 ... 2.491e+03 2.439e+03\n",
       "    MASK          (lat, lon) float32 nan nan nan nan nan ... nan nan nan nan nan\n",
       "    SLOPE         (lat, lon) float32 18.99 18.21 16.33 ... 26.08 21.17 11.44\n",
       "    ASPECT        (lat, lon) float32 197.8 204.6 186.6 ... 156.9 158.5 153.8\n",
       "    T2            (time, lat, lon) float64 ...\n",
       "    RH2           (time, lat, lon) float64 ...\n",
       "    ...            ...\n",
       "    CONDENSATION  (time, lat, lon) float64 ...\n",
       "    DEPOSITION    (time, lat, lon) float64 ...\n",
       "    REFREEZE      (time, lat, lon) float64 ...\n",
       "    subM          (time, lat, lon) float64 ...\n",
       "    Z0            (time, lat, lon) float64 ...\n",
       "    surfM         (time, lat, lon) float64 ...\n",
       "Attributes: (12/63)\n",
       "    Start_from_restart_file:                    False\n",
       "    Stake_evaluation:                           False\n",
       "    WRF_simulation:                             False\n",
       "    Compression_level:                          2\n",
       "    Slurm_use:                                  True\n",
       "    Full_fiels:                                 False\n",
       "    ...                                         ...\n",
       "    water_density:                              1000.0\n",
       "    ice_density:                                917.0\n",
       "    air_density:                                1.1\n",
       "    sigma:                                      5.67e-08\n",
       "    zero_temperature:                           273.16\n",
       "    Surface_emission_coeff:                     0.99</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-8219def2-3257-4313-b451-b7b0d23902ba' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-8219def2-3257-4313-b451-b7b0d23902ba' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>lat</span>: 14</li><li><span class='xr-has-index'>lon</span>: 30</li><li><span class='xr-has-index'>time</span>: 96444</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-d9c1c6ea-31c0-47d4-9333-1dd0ea6bc8dd' class='xr-section-summary-in' type='checkbox'  checked><label for='section-d9c1c6ea-31c0-47d4-9333-1dd0ea6bc8dd' class='xr-section-summary' >Coordinates: <span>(3)</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'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>1999-01-01 ... 2010-01-01</div><input id='attrs-85550568-5c27-4d5a-9dff-94dde4c20727' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-85550568-5c27-4d5a-9dff-94dde4c20727' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f8a0407c-ff34-44b8-8d98-886c84dd3e5f' class='xr-var-data-in' type='checkbox'><label for='data-f8a0407c-ff34-44b8-8d98-886c84dd3e5f' 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([&#x27;1999-01-01T00:00:00.000000000&#x27;, &#x27;1999-01-01T00:00:00.000000000&#x27;,\n",
       "       &#x27;1999-01-01T01:00:00.000000000&#x27;, ..., &#x27;2009-12-31T22:00:00.000000000&#x27;,\n",
       "       &#x27;2009-12-31T23:00:00.000000000&#x27;, &#x27;2010-01-01T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>46.79 46.79 46.79 ... 46.82 46.82</div><input id='attrs-3476b76d-0304-45fd-a81f-9627fa818882' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3476b76d-0304-45fd-a81f-9627fa818882' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-79e7d052-8601-4868-b832-efd450ba4096' class='xr-var-data-in' type='checkbox'><label for='data-79e7d052-8601-4868-b832-efd450ba4096' 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'><dt><span>standard_name :</span></dt><dd>lat</dd><dt><span>long_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd></dl></div><div class='xr-var-data'><pre>array([46.785139, 46.787917, 46.790694, 46.793472, 46.79625 , 46.799028,\n",
       "       46.801806, 46.804583, 46.807361, 46.810139, 46.812917, 46.815694,\n",
       "       46.818472, 46.82125 ])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>10.72 10.73 10.73 ... 10.8 10.8</div><input id='attrs-7e8230fb-3611-46d6-95a7-f602f951cc93' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7e8230fb-3611-46d6-95a7-f602f951cc93' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9cb03a33-27f4-4c90-9aac-29cbc98fc337' class='xr-var-data-in' type='checkbox'><label for='data-9cb03a33-27f4-4c90-9aac-29cbc98fc337' 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'><dt><span>standard_name :</span></dt><dd>lon</dd><dt><span>long_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>array([10.724306, 10.727083, 10.729861, 10.732639, 10.735417, 10.738194,\n",
       "       10.740972, 10.74375 , 10.746528, 10.749306, 10.752083, 10.754861,\n",
       "       10.757639, 10.760417, 10.763194, 10.765972, 10.76875 , 10.771528,\n",
       "       10.774306, 10.777083, 10.779861, 10.782639, 10.785417, 10.788195,\n",
       "       10.790972, 10.79375 , 10.796528, 10.799306, 10.802083, 10.804861])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-77703ae5-6ed1-43af-8dbf-4ab9d434cbf4' class='xr-section-summary-in' type='checkbox'  ><label for='section-77703ae5-6ed1-43af-8dbf-4ab9d434cbf4' class='xr-section-summary' >Data variables: <span>(37)</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>HGT</span></div><div class='xr-var-dims'>(lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d0fdb9c2-a133-423a-8188-b14262517b9f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d0fdb9c2-a133-423a-8188-b14262517b9f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3e01fdae-830f-4de9-a93b-f7a1472855a4' class='xr-var-data-in' type='checkbox'><label for='data-3e01fdae-830f-4de9-a93b-f7a1472855a4' 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'><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Elevation</dd></dl></div><div class='xr-var-data'><pre>array([[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833, 2761.3525],\n",
       "       [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427, 2740.5234],\n",
       "       [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124, 2793.9417],\n",
       "       ...,\n",
       "       [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153, 2687.6873],\n",
       "       [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822, 2485.7278],\n",
       "       [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 , 2438.5393]],\n",
       "      dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>MASK</span></div><div class='xr-var-dims'>(lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-aa98c1d1-0b29-4443-898c-c4bb6f2c2c18' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-aa98c1d1-0b29-4443-898c-c4bb6f2c2c18' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3546ce31-81ad-4116-868b-468c534e49f1' class='xr-var-data-in' type='checkbox'><label for='data-3546ce31-81ad-4116-868b-468c534e49f1' 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'><dt><span>units :</span></dt><dd>boolean</dd><dt><span>long_name :</span></dt><dd>Glacier mask</dd></dl></div><div class='xr-var-data'><pre>array([[nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan,  1., ..., nan, nan, nan],\n",
       "       ...,\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ...,  1.,  1., nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>SLOPE</span></div><div class='xr-var-dims'>(lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-c0d46b43-5cf6-4d2d-a8b6-18706eba104e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c0d46b43-5cf6-4d2d-a8b6-18706eba104e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0642fb3e-fc68-4a78-9b25-4b1cc64de3ce' class='xr-var-data-in' type='checkbox'><label for='data-0642fb3e-fc68-4a78-9b25-4b1cc64de3ce' 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'><dt><span>units :</span></dt><dd>degrees</dd><dt><span>long_name :</span></dt><dd>Terrain slope</dd></dl></div><div class='xr-var-data'><pre>array([[18.992346, 18.207712, 16.334656, ...,  8.222678,  2.090747,  5.847179],\n",
       "       [18.161448, 16.73097 , 10.739418, ..., 18.161112, 14.222764,  6.326509],\n",
       "       [17.340193, 15.885959, 13.283587, ..., 16.498058, 18.527357, 16.452147],\n",
       "       ...,\n",
       "       [19.914795, 20.332918, 18.758577, ..., 23.347664, 28.853716, 29.638905],\n",
       "       [21.125399, 18.110556, 14.433165, ...,  3.889344, 11.088965, 21.613993],\n",
       "       [ 7.031674,  3.742292,  5.892725, ..., 26.075985, 21.169958, 11.443837]],\n",
       "      dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ASPECT</span></div><div class='xr-var-dims'>(lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e5284536-80a7-4320-a71a-badb5fc10b70' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e5284536-80a7-4320-a71a-badb5fc10b70' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-86c397f6-00a1-4924-bb49-1185f9df94e9' class='xr-var-data-in' type='checkbox'><label for='data-86c397f6-00a1-4924-bb49-1185f9df94e9' 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'><dt><span>units :</span></dt><dd>degrees</dd><dt><span>long_name :</span></dt><dd>Aspect of slope</dd></dl></div><div class='xr-var-data'><pre>array([[197.84161 , 204.59032 , 186.60773 , ..., 136.02306 ,  79.32202 ,\n",
       "        350.1807  ],\n",
       "       [195.93721 , 198.19664 , 171.0687  , ..., 148.60442 , 141.58025 ,\n",
       "        123.47226 ],\n",
       "       [208.8029  , 182.87631 , 129.41533 , ..., 143.67503 , 138.62671 ,\n",
       "        133.12325 ],\n",
       "       ...,\n",
       "       [  5.427312,  11.55339 , 358.742   , ..., 342.4476  , 340.5044  ,\n",
       "        342.51035 ],\n",
       "       [  4.887404,   3.41128 , 348.73535 , ..., 119.23053 , 348.03845 ,\n",
       "        345.51755 ],\n",
       "       [344.67258 , 316.44064 , 272.25903 , ..., 156.92636 , 158.54759 ,\n",
       "        153.81514 ]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>T2</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-1e5ee33b-a0fd-4473-825c-aa2e07e2b5d7' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-1e5ee33b-a0fd-4473-825c-aa2e07e2b5d7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f04df761-955d-4d43-952e-bd08cfa53fe6' class='xr-var-data-in' type='checkbox'><label for='data-f04df761-955d-4d43-952e-bd08cfa53fe6' 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'><dt><span>units :</span></dt><dd>K</dd><dt><span>long_name :</span></dt><dd>Air temperature at 2 m</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>RH2</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-731c0ac0-bf6f-4391-82fb-6b74ecf9dd0a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-731c0ac0-bf6f-4391-82fb-6b74ecf9dd0a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4b5d50d6-fc7b-4d2b-ab62-2a5b09027d62' class='xr-var-data-in' type='checkbox'><label for='data-4b5d50d6-fc7b-4d2b-ab62-2a5b09027d62' 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'><dt><span>units :</span></dt><dd>%</dd><dt><span>long_name :</span></dt><dd>Relative humidity at 2 m</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>U2</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-ee41db6b-7b58-4cc2-b4c4-bc180ef5c4c5' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-ee41db6b-7b58-4cc2-b4c4-bc180ef5c4c5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9b800dbe-311e-444e-82fe-b18a91c627cf' class='xr-var-data-in' type='checkbox'><label for='data-9b800dbe-311e-444e-82fe-b18a91c627cf' 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'><dt><span>units :</span></dt><dd>m s⁻¹</dd><dt><span>long_name :</span></dt><dd>Wind velocity at 2 m</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>PRES</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-c80d9344-161f-48a5-8ea9-bee91e9854bd' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c80d9344-161f-48a5-8ea9-bee91e9854bd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a6edeb94-7f70-4cb1-be54-b60c29ae316e' class='xr-var-data-in' type='checkbox'><label for='data-a6edeb94-7f70-4cb1-be54-b60c29ae316e' 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'><dt><span>units :</span></dt><dd>hPa</dd><dt><span>long_name :</span></dt><dd>Atmospheric pressure</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>G</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-03d23490-7a85-46b1-913b-07e15cc946f2' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-03d23490-7a85-46b1-913b-07e15cc946f2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-294c2234-fc17-4193-a9bb-2718d63343af' class='xr-var-data-in' type='checkbox'><label for='data-294c2234-fc17-4193-a9bb-2718d63343af' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Incoming shortwave radiation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>RRR</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-9d5ea74a-7f67-45e5-bda6-15d5baf292f3' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9d5ea74a-7f67-45e5-bda6-15d5baf292f3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f33ca9c3-da5a-4bb5-9220-0a956aef5f4c' class='xr-var-data-in' type='checkbox'><label for='data-f33ca9c3-da5a-4bb5-9220-0a956aef5f4c' 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'><dt><span>units :</span></dt><dd>mm</dd><dt><span>long_name :</span></dt><dd>Total precipiation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>N</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-33ae21bd-56b3-4a14-b2af-2e3c558aa1d2' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-33ae21bd-56b3-4a14-b2af-2e3c558aa1d2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4958d1d4-4fe2-4238-9441-26154f905b09' class='xr-var-data-in' type='checkbox'><label for='data-4958d1d4-4fe2-4238-9441-26154f905b09' 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'><dt><span>units :</span></dt><dd>-</dd><dt><span>long_name :</span></dt><dd>Cloud fraction</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>LWin</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-6b44599e-2e42-464a-beb4-08a6612671aa' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6b44599e-2e42-464a-beb4-08a6612671aa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5bfe522b-3c3b-46ee-9fa3-afaebf464719' class='xr-var-data-in' type='checkbox'><label for='data-5bfe522b-3c3b-46ee-9fa3-afaebf464719' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Incoming longwave radiation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>RAIN</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-070f4dd8-af61-4434-b5a0-be4a8f994dce' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-070f4dd8-af61-4434-b5a0-be4a8f994dce' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-44b1dde7-b22e-4b0d-a3f0-81a1c214eefd' class='xr-var-data-in' type='checkbox'><label for='data-44b1dde7-b22e-4b0d-a3f0-81a1c214eefd' 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'><dt><span>units :</span></dt><dd>mm</dd><dt><span>long_name :</span></dt><dd>Liquid precipitation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>SNOWFALL</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f37cae8f-f2af-4c17-a717-acc2c591332f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f37cae8f-f2af-4c17-a717-acc2c591332f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f386f9ad-4ea0-43e8-82a0-753b7b8ee79e' class='xr-var-data-in' type='checkbox'><label for='data-f386f9ad-4ea0-43e8-82a0-753b7b8ee79e' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Snowfall</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>LWout</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-1d47fc23-48ff-4890-a33a-0663cb8f8e87' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-1d47fc23-48ff-4890-a33a-0663cb8f8e87' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3e763f51-b22f-4a3e-a72b-244ea737b3b0' class='xr-var-data-in' type='checkbox'><label for='data-3e763f51-b22f-4a3e-a72b-244ea737b3b0' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Outgoing longwave radiation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>H</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-5611c80d-d1a0-4032-a1b4-1b5c71f81c05' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-5611c80d-d1a0-4032-a1b4-1b5c71f81c05' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9537cb7b-737b-4933-82d0-d1d9755dc90c' class='xr-var-data-in' type='checkbox'><label for='data-9537cb7b-737b-4933-82d0-d1d9755dc90c' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Sensible heat flux</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>LE</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e5d1beb6-e96e-472d-a10e-8ce7f34343c5' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e5d1beb6-e96e-472d-a10e-8ce7f34343c5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-647576c0-4436-499a-9ea6-c539155346a5' class='xr-var-data-in' type='checkbox'><label for='data-647576c0-4436-499a-9ea6-c539155346a5' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Latent heat flux</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>B</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-0f58d57e-4ee3-4d10-b1d8-2584d1b36da1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-0f58d57e-4ee3-4d10-b1d8-2584d1b36da1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c026237c-0484-4eda-ac30-571ac45ead97' class='xr-var-data-in' type='checkbox'><label for='data-c026237c-0484-4eda-ac30-571ac45ead97' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Ground heat flux</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>QRR</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-fd15d6b6-0a70-4735-8788-9687a4b45ae5' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-fd15d6b6-0a70-4735-8788-9687a4b45ae5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e424c19a-405c-4881-b189-af5a1eaf7bf2' class='xr-var-data-in' type='checkbox'><label for='data-e424c19a-405c-4881-b189-af5a1eaf7bf2' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Rain heat flux</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>surfMB</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-18a2d6d4-4ab4-4303-8486-b77e8acc7606' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-18a2d6d4-4ab4-4303-8486-b77e8acc7606' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c5b83319-5bb5-4b76-98fd-76a340d45440' class='xr-var-data-in' type='checkbox'><label for='data-c5b83319-5bb5-4b76-98fd-76a340d45440' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Surface mass balance</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>MB</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-47e53d07-feb9-4871-bcee-a9d7de4a9dd4' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-47e53d07-feb9-4871-bcee-a9d7de4a9dd4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3ebd958b-db06-4e90-abfe-7fef74e94c2e' class='xr-var-data-in' type='checkbox'><label for='data-3ebd958b-db06-4e90-abfe-7fef74e94c2e' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Mass balance</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>Q</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-99a40dd4-e081-4b29-87e2-d55d99079e6e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-99a40dd4-e081-4b29-87e2-d55d99079e6e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-07b1c1a6-1bd0-4f24-a2f5-6afe62b918c1' class='xr-var-data-in' type='checkbox'><label for='data-07b1c1a6-1bd0-4f24-a2f5-6afe62b918c1' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Runoff</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>SNOWHEIGHT</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-6f14fb45-f3aa-485f-b7c5-35d12c156f15' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6f14fb45-f3aa-485f-b7c5-35d12c156f15' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-988d1174-346b-480c-8d3a-1cc5ce28a27f' class='xr-var-data-in' type='checkbox'><label for='data-988d1174-346b-480c-8d3a-1cc5ce28a27f' 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'><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Snowheight</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>TOTALHEIGHT</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4b1c5cc6-aa5a-4b28-923f-d46568c1b942' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4b1c5cc6-aa5a-4b28-923f-d46568c1b942' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4c342a96-9662-4a5a-ae89-f37af65a37a9' class='xr-var-data-in' type='checkbox'><label for='data-4c342a96-9662-4a5a-ae89-f37af65a37a9' 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'><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Total domain height</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>TS</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-03792116-c4ab-4ae3-99e6-5152e131b080' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-03792116-c4ab-4ae3-99e6-5152e131b080' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cb1680d5-9b6a-4516-b32c-cd07aa68f519' class='xr-var-data-in' type='checkbox'><label for='data-cb1680d5-9b6a-4516-b32c-cd07aa68f519' 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'><dt><span>units :</span></dt><dd>K</dd><dt><span>long_name :</span></dt><dd>Surface temperature</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ALBEDO</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-dd39c0d3-c99e-4c82-ad87-2348b7255648' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-dd39c0d3-c99e-4c82-ad87-2348b7255648' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e5b2554f-8967-4638-90cc-362bb671495c' class='xr-var-data-in' type='checkbox'><label for='data-e5b2554f-8967-4638-90cc-362bb671495c' 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'><dt><span>units :</span></dt><dd>-</dd><dt><span>long_name :</span></dt><dd>Albedo</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>LAYERS</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-44f2be7c-b6a1-46a7-9ce8-eb4daa3a604f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-44f2be7c-b6a1-46a7-9ce8-eb4daa3a604f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-24e9edab-efdf-4c95-acc1-fbf27f58524f' class='xr-var-data-in' type='checkbox'><label for='data-24e9edab-efdf-4c95-acc1-fbf27f58524f' 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'><dt><span>units :</span></dt><dd>-</dd><dt><span>long_name :</span></dt><dd>Number of layers</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ME</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-cee47421-aef3-4dce-b64a-cc292d43084c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-cee47421-aef3-4dce-b64a-cc292d43084c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-811ca873-b9f3-4c2b-a23c-c889c0087ba3' class='xr-var-data-in' type='checkbox'><label for='data-811ca873-b9f3-4c2b-a23c-c889c0087ba3' 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'><dt><span>units :</span></dt><dd>W m⁻²</dd><dt><span>long_name :</span></dt><dd>Available melt energy</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>intMB</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-8bb8f8d1-c607-468c-b897-7e38cebd8cfd' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-8bb8f8d1-c607-468c-b897-7e38cebd8cfd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1bed3532-3158-4d4e-b846-59ee42f09f92' class='xr-var-data-in' type='checkbox'><label for='data-1bed3532-3158-4d4e-b846-59ee42f09f92' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Internal mass balance</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>EVAPORATION</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e2440b3f-793f-4ef8-afab-e0ba96adcd49' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e2440b3f-793f-4ef8-afab-e0ba96adcd49' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-97eb5504-b714-4cce-969f-12fc1fa5b72d' class='xr-var-data-in' type='checkbox'><label for='data-97eb5504-b714-4cce-969f-12fc1fa5b72d' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Evaporation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>SUBLIMATION</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-580706de-cfad-4126-905e-2c81d14baccb' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-580706de-cfad-4126-905e-2c81d14baccb' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d3c11d03-b3d7-41ed-a34d-1b1106eb8da6' class='xr-var-data-in' type='checkbox'><label for='data-d3c11d03-b3d7-41ed-a34d-1b1106eb8da6' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Sublimation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>CONDENSATION</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-690c75fd-29a2-4e92-b052-9d8ee6133068' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-690c75fd-29a2-4e92-b052-9d8ee6133068' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4275497a-24a5-4671-9ba5-a077d73123ed' class='xr-var-data-in' type='checkbox'><label for='data-4275497a-24a5-4671-9ba5-a077d73123ed' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Condensation</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>DEPOSITION</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-33877f66-d395-4c2d-8014-e958f9c4be58' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-33877f66-d395-4c2d-8014-e958f9c4be58' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3688a679-c2a5-43f7-8159-fa718e750593' class='xr-var-data-in' type='checkbox'><label for='data-3688a679-c2a5-43f7-8159-fa718e750593' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Deposition</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>REFREEZE</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-b4ff6193-f3b6-44ec-9f66-ddb624b5d95a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-b4ff6193-f3b6-44ec-9f66-ddb624b5d95a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-63b7b4f5-7122-45ae-8ef0-8e81555570b9' class='xr-var-data-in' type='checkbox'><label for='data-63b7b4f5-7122-45ae-8ef0-8e81555570b9' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Refreezing</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>subM</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e5d594cc-a6be-44f0-b4ca-33ab3f73b09a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e5d594cc-a6be-44f0-b4ca-33ab3f73b09a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c29f4d88-aca8-420d-aec5-37135d984469' class='xr-var-data-in' type='checkbox'><label for='data-c29f4d88-aca8-420d-aec5-37135d984469' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Subsurface melt</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>Z0</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-673a1039-83c5-4248-9a81-493bc9c001d3' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-673a1039-83c5-4248-9a81-493bc9c001d3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-87aeb94c-d9c2-4a6c-9822-a4b5ec48ed6b' class='xr-var-data-in' type='checkbox'><label for='data-87aeb94c-d9c2-4a6c-9822-a4b5ec48ed6b' 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'><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Roughness length</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>surfM</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-c0a171f9-dc1b-48d9-a8db-75a9bbc2dc36' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c0a171f9-dc1b-48d9-a8db-75a9bbc2dc36' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-43644b5d-69c6-4f01-9120-a9fab5a5fac2' class='xr-var-data-in' type='checkbox'><label for='data-43644b5d-69c6-4f01-9120-a9fab5a5fac2' 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'><dt><span>units :</span></dt><dd>m w.e.</dd><dt><span>long_name :</span></dt><dd>Surface melt</dd></dl></div><div class='xr-var-data'><pre>[40506480 values with dtype=float64]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-9e679907-d8bf-48cc-aea4-7962504243ec' class='xr-section-summary-in' type='checkbox'  ><label for='section-9e679907-d8bf-48cc-aea4-7962504243ec' class='xr-section-summary' >Attributes: <span>(63)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>Start_from_restart_file :</span></dt><dd>False</dd><dt><span>Stake_evaluation :</span></dt><dd>False</dd><dt><span>WRF_simulation :</span></dt><dd>False</dd><dt><span>Compression_level :</span></dt><dd>2</dd><dt><span>Slurm_use :</span></dt><dd>True</dd><dt><span>Full_fiels :</span></dt><dd>False</dd><dt><span>Force_use_TP :</span></dt><dd>True</dd><dt><span>Force_use_N :</span></dt><dd>False</dd><dt><span>Tile_of_glacier_of_interest :</span></dt><dd>False</dd><dt><span>Time_step_input_file_seconds :</span></dt><dd>3600</dd><dt><span>Max_layers :</span></dt><dd>50</dd><dt><span>Z_measurment_height :</span></dt><dd>2.0</dd><dt><span>Stability_correction :</span></dt><dd>MO</dd><dt><span>Albedo_method :</span></dt><dd>Oerlemans98</dd><dt><span>Densification_method :</span></dt><dd>Boone</dd><dt><span>Penetrating_method :</span></dt><dd>Bintanja95</dd><dt><span>Roughness_method :</span></dt><dd>Moelg12</dd><dt><span>Saturation_water_vapour_method :</span></dt><dd>Sonntag90</dd><dt><span>Initial_snowheight :</span></dt><dd>0.3</dd><dt><span>Initial_snow_layer_heights :</span></dt><dd>0.1</dd><dt><span>Initial_glacier_height :</span></dt><dd>191.366</dd><dt><span>Initial_glacier_layer_heights :</span></dt><dd>0.5</dd><dt><span>Initial_top_density_snowpack :</span></dt><dd>300.0</dd><dt><span>Initial_bottom_density_snowpack :</span></dt><dd>600.0</dd><dt><span>Temperature_bottom :</span></dt><dd>270.16</dd><dt><span>Const_init_temp :</span></dt><dd>0.1</dd><dt><span>Center_snow_transfer_function :</span></dt><dd>-0.3193</dd><dt><span>Spread_snow_transfer_function :</span></dt><dd>2.705</dd><dt><span>Multiplication_factor_for_RRR_or_SNOWFALL :</span></dt><dd>1.092</dd><dt><span>Minimum_snow_layer_height :</span></dt><dd>0.001</dd><dt><span>Minimum_snowfall :</span></dt><dd>0.001</dd><dt><span>Remesh_method :</span></dt><dd>log_profile</dd><dt><span>First_layer_height_log_profile :</span></dt><dd>0.01</dd><dt><span>Layer_stretching_log_profile :</span></dt><dd>1.2</dd><dt><span>Merge_max :</span></dt><dd>1</dd><dt><span>Density_threshold_merging :</span></dt><dd>5</dd><dt><span>Temperature_threshold_merging :</span></dt><dd>0.01</dd><dt><span>Density_fresh_snow :</span></dt><dd>300.0</dd><dt><span>Albedo_fresh_snow :</span></dt><dd>0.925</dd><dt><span>Albedo_firn :</span></dt><dd>0.4636</dd><dt><span>Albedo_ice :</span></dt><dd>0.3708</dd><dt><span>Albedo_mod_snow_aging :</span></dt><dd>5.08</dd><dt><span>Albedo_mod_snow_depth :</span></dt><dd>5.195</dd><dt><span>Roughness_fresh_snow :</span></dt><dd>3.418</dd><dt><span>Roughness_ice :</span></dt><dd>5.215</dd><dt><span>Roughness_firn :</span></dt><dd>7.0</dd><dt><span>Aging_factor_roughness :</span></dt><dd>0.0026</dd><dt><span>Snow_ice_threshold :</span></dt><dd>900.0</dd><dt><span>lat_heat_melting :</span></dt><dd>334000.0</dd><dt><span>lat_heat_vaporize :</span></dt><dd>2500000.0</dd><dt><span>lat_heat_sublimation :</span></dt><dd>2834000.0</dd><dt><span>spec_heat_air :</span></dt><dd>1004.67</dd><dt><span>spec_heat_ice :</span></dt><dd>2050.0</dd><dt><span>spec_heat_water :</span></dt><dd>4217.0</dd><dt><span>k_i :</span></dt><dd>2.22</dd><dt><span>k_w :</span></dt><dd>0.55</dd><dt><span>k_a :</span></dt><dd>0.024</dd><dt><span>water_density :</span></dt><dd>1000.0</dd><dt><span>ice_density :</span></dt><dd>917.0</dd><dt><span>air_density :</span></dt><dd>1.1</dd><dt><span>sigma :</span></dt><dd>5.67e-08</dd><dt><span>zero_temperature :</span></dt><dd>273.16</dd><dt><span>Surface_emission_coeff :</span></dt><dd>0.99</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:       (lat: 14, lon: 30, time: 96444)\n",
       "Coordinates:\n",
       "  * time          (time) datetime64[ns] 1999-01-01 1999-01-01 ... 2010-01-01\n",
       "  * lat           (lat) float64 46.79 46.79 46.79 46.79 ... 46.82 46.82 46.82\n",
       "  * lon           (lon) float64 10.72 10.73 10.73 10.73 ... 10.8 10.8 10.8 10.8\n",
       "Data variables: (12/37)\n",
       "    HGT           (lat, lon) float32 ...\n",
       "    MASK          (lat, lon) float32 ...\n",
       "    SLOPE         (lat, lon) float32 ...\n",
       "    ASPECT        (lat, lon) float32 ...\n",
       "    T2            (time, lat, lon) float64 ...\n",
       "    RH2           (time, lat, lon) float64 ...\n",
       "    ...            ...\n",
       "    CONDENSATION  (time, lat, lon) float64 ...\n",
       "    DEPOSITION    (time, lat, lon) float64 ...\n",
       "    REFREEZE      (time, lat, lon) float64 ...\n",
       "    subM          (time, lat, lon) float64 ...\n",
       "    Z0            (time, lat, lon) float64 ...\n",
       "    surfM         (time, lat, lon) float64 ...\n",
       "Attributes: (12/63)\n",
       "    Start_from_restart_file:                    False\n",
       "    Stake_evaluation:                           False\n",
       "    WRF_simulation:                             False\n",
       "    Compression_level:                          2\n",
       "    Slurm_use:                                  True\n",
       "    Full_fiels:                                 False\n",
       "    ...                                         ...\n",
       "    water_density:                              1000.0\n",
       "    ice_density:                                917.0\n",
       "    air_density:                                1.1\n",
       "    sigma:                                      5.67e-08\n",
       "    zero_temperature:                           273.16\n",
       "    Surface_emission_coeff:                     0.99"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Load sample file\n",
    "result = xr.open_dataset(\"/home/niki/Dokumente/cosipy_test_space/data/output/\"+\n",
    "                         \"HEF_COSMO_300m_1999_2010_DREAM_19990101-20100101_num9.nc\")\n",
    "\n",
    "result"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 25,
   "id": "f59ff635-45a3-48b9-9904-bc87606ff65b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "DatetimeIndex(['2000-01-01', '2000-01-02', '2000-01-03', '2000-01-04',\n",
      "               '2000-01-05', '2000-01-06', '2000-01-07', '2000-01-08',\n",
      "               '2000-01-09', '2000-01-10',\n",
      "               ...\n",
      "               '2009-12-22', '2009-12-23', '2009-12-24', '2009-12-25',\n",
      "               '2009-12-26', '2009-12-27', '2009-12-28', '2009-12-29',\n",
      "               '2009-12-30', '2009-12-31'],\n",
      "              dtype='datetime64[ns]', length=3653, freq=None) [ 946684800000000000  946771200000000000  946857600000000000 ...\n",
      " 1262044800000000000 1262131200000000000 1262217600000000000] [ 946684800000000000  946684800000000000  946688400000000000 ...\n",
      " 1262293200000000000 1262296800000000000 1262300400000000000] [[[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",
      "\n",
      " [[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",
      "\n",
      " [[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",
      "\n",
      " ...\n",
      "\n",
      " [[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",
      "\n",
      " [[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",
      "\n",
      " [[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"
     ]
    }
   ],
   "source": [
    "## Start resampling by hand routine\n",
    "\n",
    "tsla_observations = pd.read_csv(tsl_data_file)\n",
    "tsla_observations\n",
    "\n",
    "dates, clean_day_vals, secs, holder = prereq_res(result.sel(time=slice(\"2000-01-01\",\"2009-12-31\")))\n",
    "print(dates, clean_day_vals, secs, holder)\n",
    "\n",
    "resampled_array = resample_by_hand(holder, result.sel(time=slice(\"2000-01-01\",\"2009-12-31\")).SNOWHEIGHT.values, secs, clean_day_vals)\n",
    "resampled_out = construct_resampled_ds(result.sel(time=slice(\"2000-01-01\",\"2009-12-31\")), resampled_array, dates.values)\n",
    "resampled_out['HGT'] = (('lat','lon'), result['HGT'].data)\n",
    "resampled_out['MASK'] = (('lat','lon'), result['MASK'].data)\n",
    "\n",
    "tsl_out = create_tsl_df(resampled_out, min_snowheight, tsl_method, tsl_normalize)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 73,
   "id": "a050695f-b60d-464c-8329-87779bac2f99",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>time</th>\n",
       "      <th>Med_TSL</th>\n",
       "      <th>Mean_TSL</th>\n",
       "      <th>Std_TSL</th>\n",
       "      <th>Max_TSL</th>\n",
       "      <th>Min_TSL</th>\n",
       "      <th>Flag</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>2000-01-27</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-2.013260</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>58</th>\n",
       "      <td>2000-02-28</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-2.013260</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>170</th>\n",
       "      <td>2000-06-19</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-2.013260</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>186</th>\n",
       "      <td>2000-07-05</td>\n",
       "      <td>0.058158</td>\n",
       "      <td>0.058158</td>\n",
       "      <td>-1.998620</td>\n",
       "      <td>0.072798</td>\n",
       "      <td>0.043518</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>266</th>\n",
       "      <td>2000-09-23</td>\n",
       "      <td>0.555630</td>\n",
       "      <td>0.555630</td>\n",
       "      <td>-1.956499</td>\n",
       "      <td>0.612391</td>\n",
       "      <td>0.498869</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3482</th>\n",
       "      <td>2009-07-14</td>\n",
       "      <td>0.044016</td>\n",
       "      <td>0.044016</td>\n",
       "      <td>-2.012763</td>\n",
       "      <td>0.044514</td>\n",
       "      <td>0.043518</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3498</th>\n",
       "      <td>2009-07-30</td>\n",
       "      <td>0.297253</td>\n",
       "      <td>0.297253</td>\n",
       "      <td>-2.003701</td>\n",
       "      <td>0.306813</td>\n",
       "      <td>0.287693</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3514</th>\n",
       "      <td>2009-08-15</td>\n",
       "      <td>0.482437</td>\n",
       "      <td>0.482437</td>\n",
       "      <td>-1.936264</td>\n",
       "      <td>0.559432</td>\n",
       "      <td>0.405441</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3530</th>\n",
       "      <td>2009-08-31</td>\n",
       "      <td>0.613547</td>\n",
       "      <td>0.613547</td>\n",
       "      <td>-1.943485</td>\n",
       "      <td>0.683323</td>\n",
       "      <td>0.543771</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3562</th>\n",
       "      <td>2009-10-02</td>\n",
       "      <td>0.281018</td>\n",
       "      <td>0.281018</td>\n",
       "      <td>-1.732242</td>\n",
       "      <td>0.562037</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>62 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           time   Med_TSL  Mean_TSL   Std_TSL   Max_TSL   Min_TSL  Flag\n",
       "26   2000-01-27  0.000000  0.000000 -2.013260  0.000000  0.000000   1.0\n",
       "58   2000-02-28  0.000000  0.000000 -2.013260  0.000000  0.000000   1.0\n",
       "170  2000-06-19  0.000000  0.000000 -2.013260  0.000000  0.000000   1.0\n",
       "186  2000-07-05  0.058158  0.058158 -1.998620  0.072798  0.043518   0.0\n",
       "266  2000-09-23  0.555630  0.555630 -1.956499  0.612391  0.498869   0.0\n",
       "...         ...       ...       ...       ...       ...       ...   ...\n",
       "3482 2009-07-14  0.044016  0.044016 -2.012763  0.044514  0.043518   0.0\n",
       "3498 2009-07-30  0.297253  0.297253 -2.003701  0.306813  0.287693   0.0\n",
       "3514 2009-08-15  0.482437  0.482437 -1.936264  0.559432  0.405441   0.0\n",
       "3530 2009-08-31  0.613547  0.613547 -1.943485  0.683323  0.543771   0.0\n",
       "3562 2009-10-02  0.281018  0.281018 -1.732242  0.562037  0.000000   0.0\n",
       "\n",
       "[62 rows x 7 columns]"
      ]
     },
     "execution_count": 73,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tsl_out.loc[tsl_out['time'].isin(tsla_observations['LS_DATE'])]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 74,
   "id": "81214b5c-94f5-4b47-b930-b4164dd188ec",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>LS_DATE</th>\n",
       "      <th>CC_area</th>\n",
       "      <th>CC_total_port</th>\n",
       "      <th>DC_area</th>\n",
       "      <th>DIC_area</th>\n",
       "      <th>DIC_max</th>\n",
       "      <th>IC_area</th>\n",
       "      <th>LS_ID</th>\n",
       "      <th>LS_SAT</th>\n",
       "      <th>RGI_ID</th>\n",
       "      <th>...</th>\n",
       "      <th>class_coverage</th>\n",
       "      <th>glacier_DEM_max</th>\n",
       "      <th>glacier_DEM_min</th>\n",
       "      <th>glacier_area</th>\n",
       "      <th>tool_version</th>\n",
       "      <th>scene</th>\n",
       "      <th>TSL_stddev_norm</th>\n",
       "      <th>TSL_max_norm</th>\n",
       "      <th>TSL_min_norm</th>\n",
       "      <th>TSL_normalized</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>2000-01-27</td>\n",
       "      <td>0.128650</td>\n",
       "      <td>2.635219</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.771875</td>\n",
       "      <td>3681.0</td>\n",
       "      <td>0.771875</td>\n",
       "      <td>LT05_L1TP_193027_20000127_20171212_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>60.910861</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.860917</td>\n",
       "      <td>0.054432</td>\n",
       "      <td>0.001555</td>\n",
       "      <td>0.031882</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>2000-02-28</td>\n",
       "      <td>0.027090</td>\n",
       "      <td>0.371488</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.128655</td>\n",
       "      <td>2926.0</td>\n",
       "      <td>0.128655</td>\n",
       "      <td>LT05_L1TP_193027_20000228_20171211_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>90.986062</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.854551</td>\n",
       "      <td>0.080093</td>\n",
       "      <td>0.005443</td>\n",
       "      <td>0.044323</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>2000-06-19</td>\n",
       "      <td>0.291125</td>\n",
       "      <td>3.106597</td>\n",
       "      <td>0.331701</td>\n",
       "      <td>2.193678</td>\n",
       "      <td>3682.0</td>\n",
       "      <td>1.861979</td>\n",
       "      <td>LT05_L1TP_193027_20000619_20171211_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.865007</td>\n",
       "      <td>0.210731</td>\n",
       "      <td>0.171073</td>\n",
       "      <td>0.195956</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>2000-07-05</td>\n",
       "      <td>0.390600</td>\n",
       "      <td>4.593407</td>\n",
       "      <td>0.269377</td>\n",
       "      <td>2.472011</td>\n",
       "      <td>3682.0</td>\n",
       "      <td>2.202634</td>\n",
       "      <td>LT05_L1TP_193027_20000705_20171210_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.862746</td>\n",
       "      <td>0.293157</td>\n",
       "      <td>0.235614</td>\n",
       "      <td>0.280715</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>2000-09-23</td>\n",
       "      <td>0.051182</td>\n",
       "      <td>0.696388</td>\n",
       "      <td>0.004490</td>\n",
       "      <td>0.550438</td>\n",
       "      <td>3651.0</td>\n",
       "      <td>0.545948</td>\n",
       "      <td>LT05_L1TP_193027_20000923_20171211_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>91.700154</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.851573</td>\n",
       "      <td>0.087092</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.055210</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>63</th>\n",
       "      <td>2009-07-14</td>\n",
       "      <td>1.659006</td>\n",
       "      <td>17.703254</td>\n",
       "      <td>0.595720</td>\n",
       "      <td>2.985878</td>\n",
       "      <td>3682.0</td>\n",
       "      <td>2.390158</td>\n",
       "      <td>LT05_L1TP_193027_20090714_20161023_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.870814</td>\n",
       "      <td>0.255832</td>\n",
       "      <td>0.237170</td>\n",
       "      <td>0.248834</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>64</th>\n",
       "      <td>2009-07-30</td>\n",
       "      <td>3.961122</td>\n",
       "      <td>42.794823</td>\n",
       "      <td>0.751441</td>\n",
       "      <td>3.690090</td>\n",
       "      <td>3681.0</td>\n",
       "      <td>2.938649</td>\n",
       "      <td>LT05_L1TP_193027_20090730_20161027_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.869126</td>\n",
       "      <td>0.314930</td>\n",
       "      <td>0.291602</td>\n",
       "      <td>0.309487</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>65</th>\n",
       "      <td>2009-08-15</td>\n",
       "      <td>1.184912</td>\n",
       "      <td>13.088718</td>\n",
       "      <td>1.218616</td>\n",
       "      <td>4.678650</td>\n",
       "      <td>3649.0</td>\n",
       "      <td>3.460037</td>\n",
       "      <td>LT05_L1TP_193027_20090815_20161022_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.865531</td>\n",
       "      <td>0.380249</td>\n",
       "      <td>0.342924</td>\n",
       "      <td>0.361586</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>66</th>\n",
       "      <td>2009-08-31</td>\n",
       "      <td>1.272939</td>\n",
       "      <td>14.641502</td>\n",
       "      <td>1.557160</td>\n",
       "      <td>5.734950</td>\n",
       "      <td>3696.0</td>\n",
       "      <td>4.177785</td>\n",
       "      <td>LT05_L1TP_193027_20090831_20161021_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>100.000000</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.830984</td>\n",
       "      <td>0.486003</td>\n",
       "      <td>0.372473</td>\n",
       "      <td>0.481337</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>67</th>\n",
       "      <td>2009-10-02</td>\n",
       "      <td>1.151016</td>\n",
       "      <td>14.833442</td>\n",
       "      <td>1.340535</td>\n",
       "      <td>5.064627</td>\n",
       "      <td>3651.0</td>\n",
       "      <td>3.724104</td>\n",
       "      <td>LT05_L1TP_193027_20091002_20180201_01_T1</td>\n",
       "      <td>LANDSAT_5</td>\n",
       "      <td>RGI60-11.00897</td>\n",
       "      <td>...</td>\n",
       "      <td>96.815050</td>\n",
       "      <td>3698</td>\n",
       "      <td>2412</td>\n",
       "      <td>8.014875</td>\n",
       "      <td>0.8.1</td>\n",
       "      <td>193027</td>\n",
       "      <td>-1.858797</td>\n",
       "      <td>0.532659</td>\n",
       "      <td>0.487558</td>\n",
       "      <td>0.520995</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>62 rows × 27 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "       LS_DATE   CC_area  CC_total_port   DC_area  DIC_area  DIC_max  \\\n",
       "6   2000-01-27  0.128650       2.635219  0.000000  0.771875   3681.0   \n",
       "7   2000-02-28  0.027090       0.371488  0.000000  0.128655   2926.0   \n",
       "8   2000-06-19  0.291125       3.106597  0.331701  2.193678   3682.0   \n",
       "9   2000-07-05  0.390600       4.593407  0.269377  2.472011   3682.0   \n",
       "10  2000-09-23  0.051182       0.696388  0.004490  0.550438   3651.0   \n",
       "..         ...       ...            ...       ...       ...      ...   \n",
       "63  2009-07-14  1.659006      17.703254  0.595720  2.985878   3682.0   \n",
       "64  2009-07-30  3.961122      42.794823  0.751441  3.690090   3681.0   \n",
       "65  2009-08-15  1.184912      13.088718  1.218616  4.678650   3649.0   \n",
       "66  2009-08-31  1.272939      14.641502  1.557160  5.734950   3696.0   \n",
       "67  2009-10-02  1.151016      14.833442  1.340535  5.064627   3651.0   \n",
       "\n",
       "     IC_area                                     LS_ID     LS_SAT  \\\n",
       "6   0.771875  LT05_L1TP_193027_20000127_20171212_01_T1  LANDSAT_5   \n",
       "7   0.128655  LT05_L1TP_193027_20000228_20171211_01_T1  LANDSAT_5   \n",
       "8   1.861979  LT05_L1TP_193027_20000619_20171211_01_T1  LANDSAT_5   \n",
       "9   2.202634  LT05_L1TP_193027_20000705_20171210_01_T1  LANDSAT_5   \n",
       "10  0.545948  LT05_L1TP_193027_20000923_20171211_01_T1  LANDSAT_5   \n",
       "..       ...                                       ...        ...   \n",
       "63  2.390158  LT05_L1TP_193027_20090714_20161023_01_T1  LANDSAT_5   \n",
       "64  2.938649  LT05_L1TP_193027_20090730_20161027_01_T1  LANDSAT_5   \n",
       "65  3.460037  LT05_L1TP_193027_20090815_20161022_01_T1  LANDSAT_5   \n",
       "66  4.177785  LT05_L1TP_193027_20090831_20161021_01_T1  LANDSAT_5   \n",
       "67  3.724104  LT05_L1TP_193027_20091002_20180201_01_T1  LANDSAT_5   \n",
       "\n",
       "            RGI_ID  ...  class_coverage  glacier_DEM_max  glacier_DEM_min  \\\n",
       "6   RGI60-11.00897  ...       60.910861             3698             2412   \n",
       "7   RGI60-11.00897  ...       90.986062             3698             2412   \n",
       "8   RGI60-11.00897  ...      100.000000             3698             2412   \n",
       "9   RGI60-11.00897  ...      100.000000             3698             2412   \n",
       "10  RGI60-11.00897  ...       91.700154             3698             2412   \n",
       "..             ...  ...             ...              ...              ...   \n",
       "63  RGI60-11.00897  ...      100.000000             3698             2412   \n",
       "64  RGI60-11.00897  ...      100.000000             3698             2412   \n",
       "65  RGI60-11.00897  ...      100.000000             3698             2412   \n",
       "66  RGI60-11.00897  ...      100.000000             3698             2412   \n",
       "67  RGI60-11.00897  ...       96.815050             3698             2412   \n",
       "\n",
       "    glacier_area  tool_version   scene  TSL_stddev_norm  TSL_max_norm  \\\n",
       "6       8.014875         0.8.1  193027        -1.860917      0.054432   \n",
       "7       8.014875         0.8.1  193027        -1.854551      0.080093   \n",
       "8       8.014875         0.8.1  193027        -1.865007      0.210731   \n",
       "9       8.014875         0.8.1  193027        -1.862746      0.293157   \n",
       "10      8.014875         0.8.1  193027        -1.851573      0.087092   \n",
       "..           ...           ...     ...              ...           ...   \n",
       "63      8.014875         0.8.1  193027        -1.870814      0.255832   \n",
       "64      8.014875         0.8.1  193027        -1.869126      0.314930   \n",
       "65      8.014875         0.8.1  193027        -1.865531      0.380249   \n",
       "66      8.014875         0.8.1  193027        -1.830984      0.486003   \n",
       "67      8.014875         0.8.1  193027        -1.858797      0.532659   \n",
       "\n",
       "    TSL_min_norm  TSL_normalized  \n",
       "6       0.001555        0.031882  \n",
       "7       0.005443        0.044323  \n",
       "8       0.171073        0.195956  \n",
       "9       0.235614        0.280715  \n",
       "10      0.000000        0.055210  \n",
       "..           ...             ...  \n",
       "63      0.237170        0.248834  \n",
       "64      0.291602        0.309487  \n",
       "65      0.342924        0.361586  \n",
       "66      0.372473        0.481337  \n",
       "67      0.487558        0.520995  \n",
       "\n",
       "[62 rows x 27 columns]"
      ]
     },
     "execution_count": 74,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tsla_observations.iloc[6:None]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 40,
   "id": "a930db0e-d449-4169-9669-5814a9c51874",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Required time for resample only:  0:00:11.498052\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/home/niki/anaconda3/envs/cosipy_test/lib/python3.10/site-packages/xarray/core/common.py:1127: UserWarning: Passing ``keep_attrs`` to ``resample`` has no effect and will raise an error in xarray 0.20. Pass ``keep_attrs`` directly to the applied function, e.g. ``resample(...).mean(keep_attrs=True)``.\n",
      "  warnings.warn(\n"
     ]
    },
    {
     "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'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:     (lat: 14, lon: 30, time: 3653)\n",
       "Coordinates:\n",
       "  * time        (time) datetime64[ns] 2000-01-01 2000-01-02 ... 2009-12-31\n",
       "  * lat         (lat) float64 46.79 46.79 46.79 46.79 ... 46.82 46.82 46.82\n",
       "  * lon         (lon) float64 10.72 10.73 10.73 10.73 ... 10.8 10.8 10.8 10.8\n",
       "Data variables:\n",
       "    SNOWHEIGHT  (time, lat, lon) float64 nan nan nan nan nan ... nan nan nan nan\n",
       "    HGT         (time, lat, lon) float32 3.227e+03 3.244e+03 ... 2.439e+03\n",
       "    MASK        (time, lat, lon) float32 nan nan nan nan nan ... nan nan nan nan</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-def85bd0-0240-46c4-a3d8-4052940c4daf' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-def85bd0-0240-46c4-a3d8-4052940c4daf' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>lat</span>: 14</li><li><span class='xr-has-index'>lon</span>: 30</li><li><span class='xr-has-index'>time</span>: 3653</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-c69a7224-ba04-43e7-9c8c-c1c2a093fa6e' class='xr-section-summary-in' type='checkbox'  checked><label for='section-c69a7224-ba04-43e7-9c8c-c1c2a093fa6e' class='xr-section-summary' >Coordinates: <span>(3)</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'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2000-01-01 ... 2009-12-31</div><input id='attrs-b7d392ce-0b6d-46d9-b5b2-c2f21201d59e' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b7d392ce-0b6d-46d9-b5b2-c2f21201d59e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b6774eb1-aaa2-452d-9bf7-5180061b7bfc' class='xr-var-data-in' type='checkbox'><label for='data-b6774eb1-aaa2-452d-9bf7-5180061b7bfc' 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([&#x27;2000-01-01T00:00:00.000000000&#x27;, &#x27;2000-01-02T00:00:00.000000000&#x27;,\n",
       "       &#x27;2000-01-03T00:00:00.000000000&#x27;, ..., &#x27;2009-12-29T00:00:00.000000000&#x27;,\n",
       "       &#x27;2009-12-30T00:00:00.000000000&#x27;, &#x27;2009-12-31T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>46.79 46.79 46.79 ... 46.82 46.82</div><input id='attrs-e7f6c77f-fd05-413e-ae3d-194181ac2917' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e7f6c77f-fd05-413e-ae3d-194181ac2917' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e33a8302-7967-4a12-8a15-e08fa7f56add' class='xr-var-data-in' type='checkbox'><label for='data-e33a8302-7967-4a12-8a15-e08fa7f56add' 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'><dt><span>standard_name :</span></dt><dd>lat</dd><dt><span>long_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd></dl></div><div class='xr-var-data'><pre>array([46.785139, 46.787917, 46.790694, 46.793472, 46.79625 , 46.799028,\n",
       "       46.801806, 46.804583, 46.807361, 46.810139, 46.812917, 46.815694,\n",
       "       46.818472, 46.82125 ])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>10.72 10.73 10.73 ... 10.8 10.8</div><input id='attrs-9a71197f-d907-4f19-90e5-03681d4f8647' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9a71197f-d907-4f19-90e5-03681d4f8647' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4fbc2c3b-24dc-452d-a8b9-f13b676476d1' class='xr-var-data-in' type='checkbox'><label for='data-4fbc2c3b-24dc-452d-a8b9-f13b676476d1' 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'><dt><span>standard_name :</span></dt><dd>lon</dd><dt><span>long_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>array([10.724306, 10.727083, 10.729861, 10.732639, 10.735417, 10.738194,\n",
       "       10.740972, 10.74375 , 10.746528, 10.749306, 10.752083, 10.754861,\n",
       "       10.757639, 10.760417, 10.763194, 10.765972, 10.76875 , 10.771528,\n",
       "       10.774306, 10.777083, 10.779861, 10.782639, 10.785417, 10.788195,\n",
       "       10.790972, 10.79375 , 10.796528, 10.799306, 10.802083, 10.804861])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-4a446cdd-a40d-4936-ad5e-8395a1edf2e3' class='xr-section-summary-in' type='checkbox'  checked><label for='section-4a446cdd-a40d-4936-ad5e-8395a1edf2e3' class='xr-section-summary' >Data variables: <span>(3)</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>SNOWHEIGHT</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan nan nan nan ... nan nan nan nan</div><input id='attrs-3f6dd500-44d2-4638-9b3d-8ccaff855710' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-3f6dd500-44d2-4638-9b3d-8ccaff855710' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-330293fc-e951-4da7-83fb-d4111b8ca453' class='xr-var-data-in' type='checkbox'><label for='data-330293fc-e951-4da7-83fb-d4111b8ca453' 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([[[        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,  4.44892803, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        ...,\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,  1.89278899,\n",
       "          2.01666253,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan]],\n",
       "\n",
       "       [[        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,  4.40093329, ...,         nan,\n",
       "                 nan,         nan],\n",
       "...\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,  1.56140241,\n",
       "          1.6120707 ,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan]],\n",
       "\n",
       "       [[        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan, 23.32256898, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        ...,\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,  1.6391997 ,\n",
       "          1.68981268,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>HGT</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>3.227e+03 3.244e+03 ... 2.439e+03</div><input id='attrs-36e24832-928e-434b-9ffa-422a966de223' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-36e24832-928e-434b-9ffa-422a966de223' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e25039d3-69d3-4e2d-b84c-a529e3141c01' class='xr-var-data-in' type='checkbox'><label for='data-e25039d3-69d3-4e2d-b84c-a529e3141c01' 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([[[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833,\n",
       "         2761.3525],\n",
       "        [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427,\n",
       "         2740.5234],\n",
       "        [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124,\n",
       "         2793.9417],\n",
       "        ...,\n",
       "        [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153,\n",
       "         2687.6873],\n",
       "        [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822,\n",
       "         2485.7278],\n",
       "        [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 ,\n",
       "         2438.5393]],\n",
       "\n",
       "       [[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833,\n",
       "         2761.3525],\n",
       "        [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427,\n",
       "         2740.5234],\n",
       "        [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124,\n",
       "         2793.9417],\n",
       "...\n",
       "        [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153,\n",
       "         2687.6873],\n",
       "        [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822,\n",
       "         2485.7278],\n",
       "        [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 ,\n",
       "         2438.5393]],\n",
       "\n",
       "       [[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833,\n",
       "         2761.3525],\n",
       "        [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427,\n",
       "         2740.5234],\n",
       "        [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124,\n",
       "         2793.9417],\n",
       "        ...,\n",
       "        [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153,\n",
       "         2687.6873],\n",
       "        [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822,\n",
       "         2485.7278],\n",
       "        [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 ,\n",
       "         2438.5393]]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>MASK</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>nan nan nan nan ... nan nan nan nan</div><input id='attrs-761f5120-9ebc-4968-a33a-d141291809fa' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-761f5120-9ebc-4968-a33a-d141291809fa' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fc478883-8f77-4b4c-a456-c22693976ff4' class='xr-var-data-in' type='checkbox'><label for='data-fc478883-8f77-4b4c-a456-c22693976ff4' 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([[[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]]], dtype=float32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-d09c993d-0205-4037-9029-66a84c17f31a' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-d09c993d-0205-4037-9029-66a84c17f31a' 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.Dataset>\n",
       "Dimensions:     (lat: 14, lon: 30, time: 3653)\n",
       "Coordinates:\n",
       "  * time        (time) datetime64[ns] 2000-01-01 2000-01-02 ... 2009-12-31\n",
       "  * lat         (lat) float64 46.79 46.79 46.79 46.79 ... 46.82 46.82 46.82\n",
       "  * lon         (lon) float64 10.72 10.73 10.73 10.73 ... 10.8 10.8 10.8 10.8\n",
       "Data variables:\n",
       "    SNOWHEIGHT  (time, lat, lon) float64 nan nan nan nan nan ... nan nan nan nan\n",
       "    HGT         (time, lat, lon) float32 3.227e+03 3.244e+03 ... 2.439e+03\n",
       "    MASK        (time, lat, lon) float32 nan nan nan nan nan ... nan nan nan nan"
      ]
     },
     "execution_count": 40,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## other version of resampling\n",
    "a_resampled_out = resample_output(result.sel(time=slice(\"2000-01-01\", \"2019-12-31\"))).isel(time=slice(None, -1)) # for some reason it features 2010-01-01 \n",
    "#a_tsl_out = create_tsl_df(a_resampled_out, min_snowheight, tsl_method, tsl_normalize)\n",
    "\n",
    "## xarray version of resampling\n",
    "def resample_xarray(cos_output):\n",
    "    ds = cos_output\n",
    "    ds = ds[['SNOWHEIGHT','HGT','MASK']]\n",
    "    ds_daily = ds.resample(time=\"1d\", keep_attrs=True).mean(dim='time')\n",
    "    return ds_daily\n",
    "\n",
    "test = resample_xarray(result.sel(time=slice(\"2000-01-01\", \"2019-12-31\"))).isel(time=slice(None, -1))\n",
    "test\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "id": "03afde79-ed1c-460d-8589-b4e02ddc4315",
   "metadata": {},
   "outputs": [],
   "source": [
    "## Start testing if resampling worked\n",
    "np.testing.assert_allclose(a_resampled_out.SNOWHEIGHT.values, test.SNOWHEIGHT.values, err_msg=\"These two are not equal\")\n",
    "#these two are equivalent "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "id": "09a9612b-d18d-4491-befc-7bef500c9cf8",
   "metadata": {},
   "outputs": [],
   "source": [
    "np.testing.assert_allclose(a_resampled_out.SNOWHEIGHT.values, resampled_out.SNOWHEIGHT.values, err_msg=\"These two are not equal\")\n",
    "#there seems to be a difference between the two versions - date time "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 44,
   "id": "54e4bb66-44f0-406b-8393-92516a041884",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(array([], dtype=int64), array([], dtype=int64), array([], dtype=int64))"
      ]
     },
     "execution_count": 44,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "checkup = ~np.isclose(a_resampled_out.SNOWHEIGHT.values, resampled_out.SNOWHEIGHT.values, rtol=1e-4, equal_nan=True)\n",
    "checkup\n",
    "\n",
    "\n",
    "#get indices where it is false\n",
    "indices = np.where(checkup)\n",
    "indices"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 49,
   "id": "f4af217d-4a25-407c-a143-10879186074b",
   "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'>&lt;xarray.Dataset&gt;\n",
       "Dimensions:     (lat: 14, lon: 30, time: 3653)\n",
       "Coordinates:\n",
       "  * time        (time) datetime64[ns] 2000-01-01 2000-01-02 ... 2009-12-31\n",
       "  * lat         (lat) float64 46.79 46.79 46.79 46.79 ... 46.82 46.82 46.82\n",
       "  * lon         (lon) float64 10.72 10.73 10.73 10.73 ... 10.8 10.8 10.8 10.8\n",
       "Data variables:\n",
       "    SNOWHEIGHT  (time, lat, lon) float64 nan nan nan nan nan ... nan nan nan nan\n",
       "    HGT         (time, lat, lon) float32 3.227e+03 3.244e+03 ... 2.439e+03\n",
       "    MASK        (time, lat, lon) float32 nan nan nan nan nan ... nan nan nan nan</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-9e5394bd-6439-4d2b-bde6-11676bd9c8a5' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-9e5394bd-6439-4d2b-bde6-11676bd9c8a5' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>lat</span>: 14</li><li><span class='xr-has-index'>lon</span>: 30</li><li><span class='xr-has-index'>time</span>: 3653</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-aa4275b1-6379-41d8-9bd9-61fa19af7021' class='xr-section-summary-in' type='checkbox'  checked><label for='section-aa4275b1-6379-41d8-9bd9-61fa19af7021' class='xr-section-summary' >Coordinates: <span>(3)</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'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2000-01-01 ... 2009-12-31</div><input id='attrs-6cc13b0a-92e7-41e8-97c4-dc0b3289f6b8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-6cc13b0a-92e7-41e8-97c4-dc0b3289f6b8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1f0f2381-1252-4624-8f41-c4b90c58e12d' class='xr-var-data-in' type='checkbox'><label for='data-1f0f2381-1252-4624-8f41-c4b90c58e12d' 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([&#x27;2000-01-01T00:00:00.000000000&#x27;, &#x27;2000-01-02T00:00:00.000000000&#x27;,\n",
       "       &#x27;2000-01-03T00:00:00.000000000&#x27;, ..., &#x27;2009-12-29T00:00:00.000000000&#x27;,\n",
       "       &#x27;2009-12-30T00:00:00.000000000&#x27;, &#x27;2009-12-31T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>46.79 46.79 46.79 ... 46.82 46.82</div><input id='attrs-1313966b-bd17-4efd-b0c2-664379c53e39' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1313966b-bd17-4efd-b0c2-664379c53e39' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d51627be-eed6-473f-aa64-ab97bc4aa552' class='xr-var-data-in' type='checkbox'><label for='data-d51627be-eed6-473f-aa64-ab97bc4aa552' 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([46.785139, 46.787917, 46.790694, 46.793472, 46.79625 , 46.799028,\n",
       "       46.801806, 46.804583, 46.807361, 46.810139, 46.812917, 46.815694,\n",
       "       46.818472, 46.82125 ])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>10.72 10.73 10.73 ... 10.8 10.8</div><input id='attrs-e318bb63-85a9-4571-b891-eeed80333a46' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-e318bb63-85a9-4571-b891-eeed80333a46' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9e8ab8cf-bcf5-4568-8d69-8168d83b61e9' class='xr-var-data-in' type='checkbox'><label for='data-9e8ab8cf-bcf5-4568-8d69-8168d83b61e9' 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([10.724306, 10.727083, 10.729861, 10.732639, 10.735417, 10.738194,\n",
       "       10.740972, 10.74375 , 10.746528, 10.749306, 10.752083, 10.754861,\n",
       "       10.757639, 10.760417, 10.763194, 10.765972, 10.76875 , 10.771528,\n",
       "       10.774306, 10.777083, 10.779861, 10.782639, 10.785417, 10.788195,\n",
       "       10.790972, 10.79375 , 10.796528, 10.799306, 10.802083, 10.804861])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-1c20a362-5e8e-44b8-8f8a-af342460a416' class='xr-section-summary-in' type='checkbox'  checked><label for='section-1c20a362-5e8e-44b8-8f8a-af342460a416' class='xr-section-summary' >Data variables: <span>(3)</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>SNOWHEIGHT</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>nan nan nan nan ... nan nan nan nan</div><input id='attrs-bad7479e-6b04-4f09-b17c-a88691bf1bcf' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bad7479e-6b04-4f09-b17c-a88691bf1bcf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6115c587-a6a8-43b3-aae4-9d7f7a1530e3' class='xr-var-data-in' type='checkbox'><label for='data-6115c587-a6a8-43b3-aae4-9d7f7a1530e3' 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([[[        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,  4.44892803, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        ...,\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,  1.89278899,\n",
       "          2.01666253,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan]],\n",
       "\n",
       "       [[        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,  4.40093329, ...,         nan,\n",
       "                 nan,         nan],\n",
       "...\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,  1.56140241,\n",
       "          1.6120707 ,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan]],\n",
       "\n",
       "       [[        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan, 23.32256898, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        ...,\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan],\n",
       "        [        nan,         nan,         nan, ...,  1.6391997 ,\n",
       "          1.68981268,         nan],\n",
       "        [        nan,         nan,         nan, ...,         nan,\n",
       "                 nan,         nan]]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>HGT</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>3.227e+03 3.244e+03 ... 2.439e+03</div><input id='attrs-57ca4f2b-793e-48b1-8f08-9cfe8b838f41' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-57ca4f2b-793e-48b1-8f08-9cfe8b838f41' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a722840d-9772-4538-ac0f-1fc6abb6db94' class='xr-var-data-in' type='checkbox'><label for='data-a722840d-9772-4538-ac0f-1fc6abb6db94' 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([[[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833,\n",
       "         2761.3525],\n",
       "        [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427,\n",
       "         2740.5234],\n",
       "        [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124,\n",
       "         2793.9417],\n",
       "        ...,\n",
       "        [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153,\n",
       "         2687.6873],\n",
       "        [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822,\n",
       "         2485.7278],\n",
       "        [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 ,\n",
       "         2438.5393]],\n",
       "\n",
       "       [[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833,\n",
       "         2761.3525],\n",
       "        [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427,\n",
       "         2740.5234],\n",
       "        [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124,\n",
       "         2793.9417],\n",
       "...\n",
       "        [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153,\n",
       "         2687.6873],\n",
       "        [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822,\n",
       "         2485.7278],\n",
       "        [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 ,\n",
       "         2438.5393]],\n",
       "\n",
       "       [[3227.291 , 3244.419 , 3298.0408, ..., 2770.815 , 2751.1833,\n",
       "         2761.3525],\n",
       "        [3319.9644, 3344.6167, 3412.313 , ..., 2847.6636, 2781.2427,\n",
       "         2740.5234],\n",
       "        [3419.6213, 3451.8706, 3395.9763, ..., 2959.1333, 2885.9124,\n",
       "         2793.9417],\n",
       "        ...,\n",
       "        [2840.8828, 2786.694 , 2772.5042, ..., 2536.203 , 2618.2153,\n",
       "         2687.6873],\n",
       "        [2673.9604, 2654.9146, 2665.6245, ..., 2447.7021, 2442.7822,\n",
       "         2485.7278],\n",
       "        [2576.956 , 2594.846 , 2613.9517, ..., 2565.7695, 2491.275 ,\n",
       "         2438.5393]]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>MASK</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>nan nan nan nan ... nan nan nan nan</div><input id='attrs-c3125944-2a30-4f53-862a-a14076ecca03' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-c3125944-2a30-4f53-862a-a14076ecca03' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fc623c96-66e7-4920-a5dc-917b5179f12b' class='xr-var-data-in' type='checkbox'><label for='data-fc623c96-66e7-4920-a5dc-917b5179f12b' 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([[[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "...\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]],\n",
       "\n",
       "       [[nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan,  1., ..., nan, nan, nan],\n",
       "        ...,\n",
       "        [nan, nan, nan, ..., nan, nan, nan],\n",
       "        [nan, nan, nan, ...,  1.,  1., nan],\n",
       "        [nan, nan, nan, ..., nan, nan, nan]]], dtype=float32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-c8241139-4958-4404-af1f-659add38f3d6' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-c8241139-4958-4404-af1f-659add38f3d6' 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.Dataset>\n",
       "Dimensions:     (lat: 14, lon: 30, time: 3653)\n",
       "Coordinates:\n",
       "  * time        (time) datetime64[ns] 2000-01-01 2000-01-02 ... 2009-12-31\n",
       "  * lat         (lat) float64 46.79 46.79 46.79 46.79 ... 46.82 46.82 46.82\n",
       "  * lon         (lon) float64 10.72 10.73 10.73 10.73 ... 10.8 10.8 10.8 10.8\n",
       "Data variables:\n",
       "    SNOWHEIGHT  (time, lat, lon) float64 nan nan nan nan nan ... nan nan nan nan\n",
       "    HGT         (time, lat, lon) float32 3.227e+03 3.244e+03 ... 2.439e+03\n",
       "    MASK        (time, lat, lon) float32 nan nan nan nan nan ... nan nan nan nan"
      ]
     },
     "execution_count": 49,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "## Resampling led to same values, now check the next routine which is creating the TSL data\n",
    "a_resampled_out"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 75,
   "id": "2762991b-fc4a-4f2e-9a6e-a4a0a1908821",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Calculating TSLA using conservative. Normalization is set to True.\n",
      "Max elev. 3656.1287 .\n",
      " Min elev. 2442.7822\n",
      "Numba issue exists.\n"
     ]
    }
   ],
   "source": [
    "try:\n",
    "    a_tsl_out = create_tsl_df(a_resampled_out, min_snowheight, tsl_method, tsl_normalize)\n",
    "except:\n",
    "    print(\"Numba issue exists.\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 57,
   "id": "c9cca32f-b887-4c10-8fb1-2a4797257d52",
   "metadata": {},
   "outputs": [],
   "source": [
    "#numba causes issues here, for whatever a reason. instead just use a copy of the function without numba\n",
    "def tsl_method_conservative_nojit(snowheights, hgts, mask, min_snowheight):\n",
    "    amed = np.zeros(snowheights.shape[0])\n",
    "    amean = np.zeros(snowheights.shape[0])\n",
    "    astd = np.zeros(snowheights.shape[0])\n",
    "    amax = np.zeros(snowheights.shape[0])\n",
    "    amin = np.zeros(snowheights.shape[0])\n",
    "    flag = np.zeros(snowheights.shape[0])\n",
    "\n",
    "    for n in np.arange(0, snowheights.shape[0]):\n",
    "        #min altitude of snow\n",
    "        #print(\"Processing \", n, \"out of \", snowheights.shape[0])\n",
    "        filtered_elev_snow = np.nanmin(np.where(snowheights[n,:,:]>=min_snowheight, hgts, np.nan).ravel())\n",
    "        #print(filtered_elev_snow)\n",
    "        #this line is basically redudant, ensures that snowline altitude cannot fall below glacier altitude\n",
    "        filtered_elev_snow = np.nanmax(np.append(filtered_elev_snow, np.nanmin(np.where(mask==1, hgts, np.nan)))) #numba does not support np.maximum\n",
    "        #print(filtered_elev_snow)\n",
    "        #now for snow-free surfaces\n",
    "        filtered_elev_nosnow = np.nanmax(np.where(snowheights[n,:,:]<min_snowheight, hgts, np.nan))\n",
    "        #print(filtered_elev_nosnow)\n",
    "        if np.isnan(filtered_elev_nosnow):\n",
    "            #print(\"Glacier seems to be fully snow-covered. Assigning minimum elevation.\")\n",
    "            filtered_elev_nosnow = np.nanmin(np.where(mask==1, hgts, np.nan))\n",
    "            #print(filtered_elev_nosnow,\"m a.s.l.\")\n",
    "            flag[n] = 1\n",
    "        else:\n",
    "            #technically also redudant, ascertain that max. elevation of glacier is not exceeded\n",
    "            filtered_elev_nosnow = np.nanmin(np.append(filtered_elev_nosnow, np.nanmax(np.where(mask==1, hgts, np.nan)))) #numba does not support minimum/maximum of numpy\n",
    "        amed[n] = np.nanmedian(np.append(filtered_elev_snow, filtered_elev_nosnow))\n",
    "        amean[n] = np.nanmean(np.append(filtered_elev_snow, filtered_elev_nosnow))\n",
    "        astd[n] = np.nanstd(np.append(filtered_elev_snow, filtered_elev_nosnow)) #might be a bit large.. std of all values essentially function of resolution  \n",
    "        amax[n] = np.nanmax(np.append(filtered_elev_snow, filtered_elev_nosnow))\n",
    "        amin[n] = np.nanmin(np.append(filtered_elev_snow, filtered_elev_nosnow))\n",
    "\n",
    "    return (amed, amean, astd, amax, amin, flag)\n",
    "\n",
    "def calculate_tsl_byhand_nojit(snowheights,hgts,mask, min_snowheight,tsl_method, tsl_normalize):\n",
    "    max_elev = np.nanmax(np.where(mask==1, hgts, np.nan))\n",
    "    min_elev = np.nanmin(np.where(mask==1, hgts, np.nan))\n",
    "    print(\"Calculating TSLA using {}. Normalization is set to {}.\".format(tsl_method, tsl_normalize))\n",
    "    print(\"Max elev.\", max_elev,\".\\n Min elev.\", min_elev)\n",
    "    if tsl_method == 'mantra':\n",
    "        amed, amean, astd, amax, amin, flag = tsl_method_mantra(snowheights, hgts, mask, min_snowheight)\n",
    "    elif tsl_method == 'conservative':\n",
    "        amed, amean, astd, amax, amin, flag = tsl_method_conservative_nojit(snowheights, hgts, mask, min_snowheight)\n",
    "    else:\n",
    "        amed, amean, astd, amax, amin, flag = tsl_method_gridsearch(snowheights, hgts, mask, min_snowheight)\n",
    "    if tsl_normalize:\n",
    "        #Start normalizing:\n",
    "        amed_norm = (amed - min_elev) / (max_elev - min_elev)\n",
    "        amean_norm = (amean - min_elev) / (max_elev - min_elev)\n",
    "        astd_norm = (astd - min_elev) / (max_elev - min_elev)\n",
    "        amax_norm = (amax - min_elev) / (max_elev - min_elev)\n",
    "        amin_norm = (amin - min_elev) / (max_elev - min_elev)\n",
    "        return (amed_norm, amean_norm, astd_norm, amax_norm, amin_norm, flag)\n",
    "    else:\n",
    "        return (amed, amean, astd, amax, amin, flag)\n",
    "\n",
    "def create_tsl_df_nojit(cos_output,min_snowheight, tsl_method, tsl_normalize):\n",
    "    times = datetime.now()\n",
    "    amed,amean,astd,amax,amin,flag = calculate_tsl_byhand_nojit(cos_output.SNOWHEIGHT.values, cos_output.HGT.values, cos_output.MASK.values, min_snowheight, tsl_method, tsl_normalize)\n",
    "    print(amed)\n",
    "    tsl_df = pd.DataFrame({'time':pd.to_datetime(cos_output.time.values),\n",
    "                            'Med_TSL':amed,\n",
    "                            'Mean_TSL':amean,\n",
    "                            'Std_TSL':astd,\n",
    "                            'Max_TSL':amax,\n",
    "                            'Min_TSL':amin,\n",
    "                            'Flag': flag})\n",
    "    print(\"Time required for calculating TSL only :\", datetime.now()-times)\n",
    "    return tsl_df"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 58,
   "id": "6a0cf5fe-3fe9-4682-bec4-2fc89470943e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Calculating TSLA using conservative. Normalization is set to True.\n",
      "Max elev. 3656.1287 .\n",
      " Min elev. 2442.7822\n"
     ]
    },
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "/tmp/ipykernel_19149/3061086146.py:19: RuntimeWarning: All-NaN slice encountered\n",
      "  filtered_elev_nosnow = np.nanmax(np.where(snowheights[n,:,:]<min_snowheight, hgts, np.nan))\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[0. 0. 0. ... 0. 0. 0.]\n",
      "Time required for calculating TSL only : 0:00:43.580121\n"
     ]
    }
   ],
   "source": [
    "a_tsl_out = create_tsl_df_nojit(a_resampled_out, min_snowheight, tsl_method, tsl_normalize)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 59,
   "id": "0cc31017-519f-4be8-a5b7-1942b3eab0b4",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>time</th>\n",
       "      <th>Med_TSL</th>\n",
       "      <th>Mean_TSL</th>\n",
       "      <th>Std_TSL</th>\n",
       "      <th>Max_TSL</th>\n",
       "      <th>Min_TSL</th>\n",
       "      <th>Flag</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2000-01-01</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2000-01-02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2000-01-03</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2000-01-04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2000-01-05</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3648</th>\n",
       "      <td>2009-12-27</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3649</th>\n",
       "      <td>2009-12-28</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3650</th>\n",
       "      <td>2009-12-29</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3651</th>\n",
       "      <td>2009-12-30</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3652</th>\n",
       "      <td>2009-12-31</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3653 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           time  Med_TSL  Mean_TSL  Std_TSL  Max_TSL  Min_TSL  Flag\n",
       "0    2000-01-01      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "1    2000-01-02      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "2    2000-01-03      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3    2000-01-04      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "4    2000-01-05      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "...         ...      ...       ...      ...      ...      ...   ...\n",
       "3648 2009-12-27      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3649 2009-12-28      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3650 2009-12-29      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3651 2009-12-30      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3652 2009-12-31      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "\n",
       "[3653 rows x 7 columns]"
      ]
     },
     "execution_count": 59,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "a_tsl_out"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 60,
   "id": "44a8f599-8f4f-4ad2-a594-96f52d94c58c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>time</th>\n",
       "      <th>Med_TSL</th>\n",
       "      <th>Mean_TSL</th>\n",
       "      <th>Std_TSL</th>\n",
       "      <th>Max_TSL</th>\n",
       "      <th>Min_TSL</th>\n",
       "      <th>Flag</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>2000-01-01</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>2000-01-02</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>2000-01-03</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>2000-01-04</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>2000-01-05</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>...</th>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "      <td>...</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3648</th>\n",
       "      <td>2009-12-27</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3649</th>\n",
       "      <td>2009-12-28</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3650</th>\n",
       "      <td>2009-12-29</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3651</th>\n",
       "      <td>2009-12-30</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3652</th>\n",
       "      <td>2009-12-31</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>-2.01326</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>1.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>3653 rows × 7 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "           time  Med_TSL  Mean_TSL  Std_TSL  Max_TSL  Min_TSL  Flag\n",
       "0    2000-01-01      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "1    2000-01-02      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "2    2000-01-03      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3    2000-01-04      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "4    2000-01-05      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "...         ...      ...       ...      ...      ...      ...   ...\n",
       "3648 2009-12-27      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3649 2009-12-28      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3650 2009-12-29      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3651 2009-12-30      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "3652 2009-12-31      0.0       0.0 -2.01326      0.0      0.0   1.0\n",
       "\n",
       "[3653 rows x 7 columns]"
      ]
     },
     "execution_count": 60,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "tsl_out"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 61,
   "id": "290f7f1a-cb91-404e-a8c9-9b28aa6bc0ae",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "True"
      ]
     },
     "execution_count": 61,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "np.allclose(a_tsl_out['Med_TSL'], tsl_out['Med_TSL'])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "id": "4d1e2386-0221-4051-a213-dd025c63b4ba",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Column  Med_TSL are equal?\n",
      "True\n",
      "Column  Mean_TSL are equal?\n",
      "True\n",
      "Column  Std_TSL are equal?\n",
      "True\n",
      "Column  Max_TSL are equal?\n",
      "True\n",
      "Column  Min_TSL are equal?\n",
      "True\n",
      "Column  Flag are equal?\n",
      "True\n"
     ]
    }
   ],
   "source": [
    "for column in tsl_out.columns:\n",
    "    if column != \"time\":\n",
    "        print(\"Column \", column, \"are equal?\")\n",
    "        print(np.allclose(tsl_out[column],a_tsl_out[column]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "bb4eb471-6210-453d-8a73-1366ff9167ff",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}