5Generate product_glacier_level_mb.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "09d0b8e1-9657-48a4-bf7c-860a4e28af7a",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'D:\\\\Py_works\\\\00Dh_process'"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pandas as pd\n",
    "import os\n",
    "import numpy as np\n",
    "import proplot as plot\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib as mpl\n",
    "import geopandas as gpd\n",
    "import xarray as xr\n",
    "import rioxarray as rioxr\n",
    "os.getcwd()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "6446aa27-6be8-4455-9555-a4a4052365dd",
   "metadata": {},
   "outputs": [],
   "source": [
    "pth_nc='D:\\\\Key_works\\\\20231212dh_calculate\\\\xnhl_test\\\\gla_level_dh_final'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "8c8a5b13-5a6a-479e-8537-789f7776a06f",
   "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>Unnamed: 0</th>\n",
       "      <th>group_id</th>\n",
       "      <th>mean</th>\n",
       "      <th>median</th>\n",
       "      <th>mean_fill</th>\n",
       "      <th>median_fill</th>\n",
       "      <th>Glc_Area</th>\n",
       "      <th>Mean_Elev</th>\n",
       "      <th>Mean_Slp</th>\n",
       "      <th>Mean_Asp</th>\n",
       "      <th>MA_Elev</th>\n",
       "      <th>basin</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>G101096E33298N</td>\n",
       "      <td>g0018</td>\n",
       "      <td>-2.247736</td>\n",
       "      <td>-4.501035</td>\n",
       "      <td>-2.532606</td>\n",
       "      <td>-4.124494</td>\n",
       "      <td>293031.2539</td>\n",
       "      <td>4942.3</td>\n",
       "      <td>39.4</td>\n",
       "      <td>315.8</td>\n",
       "      <td>4939.9</td>\n",
       "      <td>Yellow</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "       Unnamed: 0 group_id      mean    median  mean_fill  median_fill  \\\n",
       "0  G101096E33298N    g0018 -2.247736 -4.501035  -2.532606    -4.124494   \n",
       "\n",
       "      Glc_Area  Mean_Elev  Mean_Slp  Mean_Asp  MA_Elev   basin  \n",
       "0  293031.2539     4942.3      39.4     315.8   4939.9  Yellow  "
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pth='D:/Key_works/20231212dh_calculate/xnhl_test'\n",
    "df_log=pd.read_csv(pth+'/dh-effe-stat.csv')\n",
    "df_log.head(1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "744c1def-194e-4525-b656-9fd8a806b654",
   "metadata": {},
   "outputs": [],
   "source": [
    "stat_df=list(df_log.iloc[:,0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "8afaab42-e775-4410-aa00-d761a0edfeb3",
   "metadata": {},
   "outputs": [],
   "source": [
    "my_prj='''PROJCS[\"Albers_Conic_Equal_Area\",GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]],PROJECTION[\"Albers_Conic_Equal_Area\"],PARAMETER[\"latitude_of_center\",0],PARAMETER[\"longitude_of_center\",105],PARAMETER[\"standard_parallel_1\",25],PARAMETER[\"standard_parallel_2\",47],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"Easting\",EAST],AXIS[\"Northing\",NORTH]]'''"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "id": "85d9aafc-3780-4240-b5f4-b1be07099292",
   "metadata": {},
   "outputs": [],
   "source": [
    "newp0=pth_nc+'\\\\nc_yr'\n",
    "newp1=pth_nc+'\\\\nc_yr_mb'\n",
    "if os.path.exists(newp1):pass\n",
    "else:os.makedirs(newp1)\n",
    "\n",
    "aaa=[]\n",
    "for glai in stat_df: #G075167E36993N stat_df\n",
    "    innm=newp0+'/'+glai+'.nc'\n",
    "    outnm=newp1+'/'+glai+'.nc'\n",
    "    ds_gla=xr.open_dataset(innm)\n",
    "\n",
    "    ds_gla['mb']=(ds_gla['topo_dh']/(2000-ds_gla['dxt_id'])*0.85)\n",
    "    t=(2000-ds_gla['dxt_id']).mean()\n",
    "    pi=0.85;pw=1;dp=0.06\n",
    "    mb_err=(((ds_gla['topo_dh'].mean()/t)*(dp/pw))**2+((ds_gla.attrs['dh error']/t)*(pi/pw))**2)**0.5\n",
    "    ds_gla['mb'].attrs={'long_name' :'Mass balance (by using [topo_dh])',\n",
    "                        'units' :'m w.e. yr-1'}\n",
    "    ds_gla['dxt_id'].attrs={'long_name' :'Photography time of aerial photos (Topo maps time)',}\n",
    "    ds_gla.attrs['mb error']=round(float(mb_err.values),4)\n",
    "    aaa.append([ds_gla['mb'].mean(),mb_err.values])\n",
    "    ds_gla0=ds_gla.rio.write_crs(my_prj)\n",
    "    encd={'topo_dh':{\"zlib\": True, \"complevel\": 8},'topo_dh_fill':{\"zlib\": True, \"complevel\": 8},\n",
    "         'elevation':{\"zlib\": True, \"complevel\": 8},'dxt_id':{\"zlib\": True, \"complevel\": 8},\n",
    "         'mb':{\"zlib\": True, \"complevel\": 8},}\n",
    "    ds_gla0.to_netcdf(outnm)\n",
    "    #break"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "c50a6e17-7f72-43c0-92b3-20f75d3130c5",
   "metadata": {},
   "source": [
    "### 重投影"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "2e57c018-61e0-45f9-8b24-854607b998e3",
   "metadata": {},
   "outputs": [],
   "source": [
    "newp1=pth_nc+'\\\\nc_yr_mb'\n",
    "newp2=pth_nc+'\\\\Glacier-level mass change'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2fb41e94-acc1-4b90-9b8b-00249982c067",
   "metadata": {},
   "outputs": [],
   "source": [
    "flst=[i for i in os.listdir(newp1) if i[-3:]=='.nc']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "f54c3f62-a78f-4119-bb77-73680dce5743",
   "metadata": {},
   "outputs": [],
   "source": [
    "for glai in flst: #G075167E36993N stat_df\n",
    "    innm=newp1+'/'+glai\n",
    "    outnm=newp2+'/'+glai\n",
    "    ds_gla=xr.open_dataset(innm)\n",
    "    ds_gla0=ds_gla.rio.write_crs(my_prj)\n",
    "    ds_gla1=ds_gla0.rio.reproject('epsg:4326')\n",
    "    encd={'topo_dh':{\"zlib\": True, \"complevel\": 8},'topo_dh_fill':{\"zlib\": True, \"complevel\": 8},\n",
    "         'elevation':{\"zlib\": True, \"complevel\": 8},'dxt_id':{\"zlib\": True, \"complevel\": 8},\n",
    "         'mb':{\"zlib\": True, \"complevel\": 8},}\n",
    "    ds_gla1.to_netcdf(outnm)\n",
    "    #aaa"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "1c953955-495a-420d-b5e3-2b66b4607de1",
   "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[data-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 !important;\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-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\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",
       ".xr-index-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",
       ".xr-index-data-in:checked ~ .xr-index-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-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-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",
       ".xr-no-icon {\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; Size: 74kB\n",
       "Dimensions:       (x: 56, y: 41)\n",
       "Coordinates:\n",
       "  * x             (x) float64 448B 78.74 78.74 78.74 78.74 ... 78.76 78.76 78.76\n",
       "  * y             (y) float64 328B 31.19 31.19 31.19 31.19 ... 31.18 31.18 31.18\n",
       "    spatial_ref   int32 4B 0\n",
       "Data variables:\n",
       "    topo_dh       (y, x) float64 18kB nan nan nan nan nan ... nan nan nan nan\n",
       "    elevation     (y, x) float32 9kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    topo_dh_fill  (y, x) float64 18kB nan nan nan nan nan ... nan nan nan nan\n",
       "    dxt_id        (y, x) float32 9kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    mb            (y, x) float64 18kB nan nan nan nan nan ... nan nan nan nan\n",
       "Attributes:\n",
       "    transform:           [ 3.00000000e+01  0.00000000e+00 -2.43928035e+06  0....\n",
       "    res:                 [30. 30.]\n",
       "    nodatavals:          nan\n",
       "    scales:              1.0\n",
       "    offsets:             0.0\n",
       "    topo process group:  g0309\n",
       "    effective samples:   17.4075\n",
       "    dh error:            2.5313\n",
       "    mb error:            0.2098</pre><div class='xr-wrap' style='display:none'><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-1326a971-3aa1-4493-b601-c2f91cd2f97c' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-1326a971-3aa1-4493-b601-c2f91cd2f97c' 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'>x</span>: 56</li><li><span class='xr-has-index'>y</span>: 41</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-bba61f04-c47b-4212-8329-6eb7f5782a5c' class='xr-section-summary-in' type='checkbox'  checked><label for='section-bba61f04-c47b-4212-8329-6eb7f5782a5c' 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'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>78.74 78.74 78.74 ... 78.76 78.76</div><input id='attrs-f0aebe42-6571-487c-97e9-137f31b3c16c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f0aebe42-6571-487c-97e9-137f31b3c16c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8d321c2c-01a2-4001-b9c0-83d17e63c97a' class='xr-var-data-in' type='checkbox'><label for='data-8d321c2c-01a2-4001-b9c0-83d17e63c97a' 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>axis :</span></dt><dd>X</dd><dt><span>long_name :</span></dt><dd>longitude</dd><dt><span>standard_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>array([78.743441, 78.743752, 78.744064, 78.744376, 78.744687, 78.744999,\n",
       "       78.745311, 78.745622, 78.745934, 78.746246, 78.746558, 78.746869,\n",
       "       78.747181, 78.747493, 78.747804, 78.748116, 78.748428, 78.748739,\n",
       "       78.749051, 78.749363, 78.749675, 78.749986, 78.750298, 78.75061 ,\n",
       "       78.750921, 78.751233, 78.751545, 78.751856, 78.752168, 78.75248 ,\n",
       "       78.752792, 78.753103, 78.753415, 78.753727, 78.754038, 78.75435 ,\n",
       "       78.754662, 78.754973, 78.755285, 78.755597, 78.755908, 78.75622 ,\n",
       "       78.756532, 78.756844, 78.757155, 78.757467, 78.757779, 78.75809 ,\n",
       "       78.758402, 78.758714, 78.759025, 78.759337, 78.759649, 78.759961,\n",
       "       78.760272, 78.760584])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>31.19 31.19 31.19 ... 31.18 31.18</div><input id='attrs-a25a6d2c-c570-4a32-a393-02b6e72835f4' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-a25a6d2c-c570-4a32-a393-02b6e72835f4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8fda4ed5-2623-4773-90ab-f32a1dc460a4' class='xr-var-data-in' type='checkbox'><label for='data-8fda4ed5-2623-4773-90ab-f32a1dc460a4' 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>axis :</span></dt><dd>Y</dd><dt><span>long_name :</span></dt><dd>latitude</dd><dt><span>standard_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd></dl></div><div class='xr-var-data'><pre>array([31.187515, 31.187203, 31.186891, 31.18658 , 31.186268, 31.185956,\n",
       "       31.185644, 31.185333, 31.185021, 31.184709, 31.184398, 31.184086,\n",
       "       31.183774, 31.183463, 31.183151, 31.182839, 31.182527, 31.182216,\n",
       "       31.181904, 31.181592, 31.181281, 31.180969, 31.180657, 31.180346,\n",
       "       31.180034, 31.179722, 31.17941 , 31.179099, 31.178787, 31.178475,\n",
       "       31.178164, 31.177852, 31.17754 , 31.177229, 31.176917, 31.176605,\n",
       "       31.176293, 31.175982, 31.17567 , 31.175358, 31.175047])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-9c79cd9d-bdf1-41cd-8505-2c51004e938b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9c79cd9d-bdf1-41cd-8505-2c51004e938b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ab4b7c94-8dba-4376-b908-093410b0b8fc' class='xr-var-data-in' type='checkbox'><label for='data-ab4b7c94-8dba-4376-b908-093410b0b8fc' 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>crs_wkt :</span></dt><dd>GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314245179</dd><dt><span>inverse_flattening :</span></dt><dd>298.257223563</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>WGS 84</dd><dt><span>longitude_of_prime_meridian :</span></dt><dd>0.0</dd><dt><span>prime_meridian_name :</span></dt><dd>Greenwich</dd><dt><span>geographic_crs_name :</span></dt><dd>WGS 84</dd><dt><span>horizontal_datum_name :</span></dt><dd>World Geodetic System 1984</dd><dt><span>grid_mapping_name :</span></dt><dd>latitude_longitude</dd><dt><span>spatial_ref :</span></dt><dd>GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0,AUTHORITY[&quot;EPSG&quot;,&quot;8901&quot;]],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AXIS[&quot;Latitude&quot;,NORTH],AXIS[&quot;Longitude&quot;,EAST],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]]</dd><dt><span>GeoTransform :</span></dt><dd>78.74328468830531 0.00031169886301098257 0.0 31.18767046219548 0.0 -0.00031169886301098257</dd></dl></div><div class='xr-var-data'><pre>array(0)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-b2148172-631d-4705-9a77-813c778a12e2' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b2148172-631d-4705-9a77-813c778a12e2' class='xr-section-summary' >Data variables: <span>(5)</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>topo_dh</span></div><div class='xr-var-dims'>(y, x)</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-f74dd4e1-45ce-4a58-a46c-0a36f66c6741' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f74dd4e1-45ce-4a58-a46c-0a36f66c6741' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cf33aeab-abd0-4aae-90e3-408921079880' class='xr-var-data-in' type='checkbox'><label for='data-cf33aeab-abd0-4aae-90e3-408921079880' 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>long_name :</span></dt><dd>Elevation difference from Topo DEM and SRTM DEM (viod version, v1)</dd><dt><span>units :</span></dt><dd>m</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, nan, ..., nan, nan, nan],\n",
       "       ...,\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>elevation</span></div><div class='xr-var-dims'>(y, x)</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-560b1097-7be2-4bc7-8c95-5cd5bb9931e9' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-560b1097-7be2-4bc7-8c95-5cd5bb9931e9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a49ea7ff-12b2-46c7-a15e-dba98f7a9ccc' class='xr-var-data-in' type='checkbox'><label for='data-a49ea7ff-12b2-46c7-a15e-dba98f7a9ccc' 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>long_name :</span></dt><dd>elevation from viod filled SRTM DEM</dd><dt><span>units :</span></dt><dd>m a.s.l</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, nan, ..., nan, nan, nan],\n",
       "       ...,\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>topo_dh_fill</span></div><div class='xr-var-dims'>(y, x)</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-6e88c5bd-6cc3-40f4-8d70-6d3d53df8ff4' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6e88c5bd-6cc3-40f4-8d70-6d3d53df8ff4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f07db405-f97e-48a4-8e26-df28a047f393' class='xr-var-data-in' type='checkbox'><label for='data-f07db405-f97e-48a4-8e26-df28a047f393' 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>long_name :</span></dt><dd>Filled elevation difference from Topo DEM and SRTM DEM</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>Fill method :</span></dt><dd>Normalized regional hypsometric interpolation</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, nan, ..., nan, nan, nan],\n",
       "       ...,\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>dxt_id</span></div><div class='xr-var-dims'>(y, x)</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-4954b516-cc1d-4840-a73f-9e8ea72a1c59' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4954b516-cc1d-4840-a73f-9e8ea72a1c59' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-32e7c37d-3a8b-4384-959f-27b94168cfce' class='xr-var-data-in' type='checkbox'><label for='data-32e7c37d-3a8b-4384-959f-27b94168cfce' 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>long_name :</span></dt><dd>Photography time of aerial photos (Topo maps time)</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, nan, ..., nan, nan, nan],\n",
       "       ...,\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>mb</span></div><div class='xr-var-dims'>(y, x)</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-e89afa4a-e861-4807-b369-5b3c28e04d46' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e89afa4a-e861-4807-b369-5b3c28e04d46' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1aa95ce1-794e-46b5-bde7-a15ca6e9889e' class='xr-var-data-in' type='checkbox'><label for='data-1aa95ce1-794e-46b5-bde7-a15ca6e9889e' 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>long_name :</span></dt><dd>Mass balance (by using [topo_dh])</dd><dt><span>units :</span></dt><dd>m w.e. yr-1</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, nan, ..., nan, nan, nan],\n",
       "       ...,\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan],\n",
       "       [nan, nan, nan, ..., nan, nan, nan]])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-cbf5bb9c-ea2f-4f7d-ab3d-620986ba7d99' class='xr-section-summary-in' type='checkbox'  ><label for='section-cbf5bb9c-ea2f-4f7d-ab3d-620986ba7d99' class='xr-section-summary' >Indexes: <span>(2)</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-index-name'><div>x</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-5c721bef-5469-409a-9c90-2a2111153f15' class='xr-index-data-in' type='checkbox'/><label for='index-5c721bef-5469-409a-9c90-2a2111153f15' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([78.74344053773682, 78.74375223659983, 78.74406393546283,\n",
       "       78.74437563432585, 78.74468733318886, 78.74499903205187,\n",
       "       78.74531073091488,  78.7456224297779, 78.74593412864091,\n",
       "       78.74624582750391, 78.74655752636693, 78.74686922522994,\n",
       "       78.74718092409294, 78.74749262295596, 78.74780432181898,\n",
       "       78.74811602068198, 78.74842771954499, 78.74873941840801,\n",
       "       78.74905111727101, 78.74936281613402, 78.74967451499704,\n",
       "       78.74998621386004, 78.75029791272306, 78.75060961158607,\n",
       "       78.75092131044909, 78.75123300931209,  78.7515447081751,\n",
       "       78.75185640703812, 78.75216810590112, 78.75247980476414,\n",
       "       78.75279150362715, 78.75310320249015, 78.75341490135317,\n",
       "       78.75372660021618, 78.75403829907918,  78.7543499979422,\n",
       "       78.75466169680521, 78.75497339566822, 78.75528509453123,\n",
       "       78.75559679339425, 78.75590849225726, 78.75622019112026,\n",
       "       78.75653188998328,  78.7568435888463,  78.7571552877093,\n",
       "       78.75746698657231, 78.75777868543533, 78.75809038429833,\n",
       "       78.75840208316134, 78.75871378202436, 78.75902548088736,\n",
       "       78.75933717975037, 78.75964887861339, 78.75996057747639,\n",
       "        78.7602722763394, 78.76058397520242],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;x&#x27;))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>y</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-0a4ff80e-e6b3-4c91-80f4-d30316690045' class='xr-index-data-in' type='checkbox'/><label for='index-0a4ff80e-e6b3-4c91-80f4-d30316690045' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([31.187514612763977, 31.187202913900965, 31.186891215037953,\n",
       "       31.186579516174945, 31.186267817311933,  31.18595611844892,\n",
       "        31.18564441958591,   31.1853327207229,  31.18502102185989,\n",
       "       31.184709322996877, 31.184397624133865, 31.184085925270857,\n",
       "       31.183774226407845, 31.183462527544833,  31.18315082868182,\n",
       "       31.182839129818813,   31.1825274309558,  31.18221573209279,\n",
       "       31.181904033229777,  31.18159233436677, 31.181280635503757,\n",
       "       31.180968936640745, 31.180657237777734, 31.180345538914725,\n",
       "       31.180033840051713,   31.1797221411887,  31.17941044232569,\n",
       "        31.17909874346268,  31.17878704459967, 31.178475345736658,\n",
       "       31.178163646873646, 31.177851948010638, 31.177540249147626,\n",
       "       31.177228550284614, 31.176916851421602, 31.176605152558594,\n",
       "       31.176293453695582,  31.17598175483257, 31.175670055969558,\n",
       "        31.17535835710655, 31.175046658243538],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;y&#x27;))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-5f2d7e7d-ed95-41e7-a10c-53f25ae7d2b1' class='xr-section-summary-in' type='checkbox'  checked><label for='section-5f2d7e7d-ed95-41e7-a10c-53f25ae7d2b1' class='xr-section-summary' >Attributes: <span>(9)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>transform :</span></dt><dd>[ 3.00000000e+01  0.00000000e+00 -2.43928035e+06  0.00000000e+00\n",
       " -3.00000000e+01  3.62912094e+06]</dd><dt><span>res :</span></dt><dd>[30. 30.]</dd><dt><span>nodatavals :</span></dt><dd>nan</dd><dt><span>scales :</span></dt><dd>1.0</dd><dt><span>offsets :</span></dt><dd>0.0</dd><dt><span>topo process group :</span></dt><dd>g0309</dd><dt><span>effective samples :</span></dt><dd>17.4075</dd><dt><span>dh error :</span></dt><dd>2.5313</dd><dt><span>mb error :</span></dt><dd>0.2098</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset> Size: 74kB\n",
       "Dimensions:       (x: 56, y: 41)\n",
       "Coordinates:\n",
       "  * x             (x) float64 448B 78.74 78.74 78.74 78.74 ... 78.76 78.76 78.76\n",
       "  * y             (y) float64 328B 31.19 31.19 31.19 31.19 ... 31.18 31.18 31.18\n",
       "    spatial_ref   int32 4B 0\n",
       "Data variables:\n",
       "    topo_dh       (y, x) float64 18kB nan nan nan nan nan ... nan nan nan nan\n",
       "    elevation     (y, x) float32 9kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    topo_dh_fill  (y, x) float64 18kB nan nan nan nan nan ... nan nan nan nan\n",
       "    dxt_id        (y, x) float32 9kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    mb            (y, x) float64 18kB nan nan nan nan nan ... nan nan nan nan\n",
       "Attributes:\n",
       "    transform:           [ 3.00000000e+01  0.00000000e+00 -2.43928035e+06  0....\n",
       "    res:                 [30. 30.]\n",
       "    nodatavals:          nan\n",
       "    scales:              1.0\n",
       "    offsets:             0.0\n",
       "    topo process group:  g0309\n",
       "    effective samples:   17.4075\n",
       "    dh error:            2.5313\n",
       "    mb error:            0.2098"
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds_gla1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "91f38f25-3a35-4a6f-bca5-6e70614bfe84",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3e168c86-3033-4338-8a89-15d6e2463adc",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f391b8d6-d638-44be-b1ef-b4efdd5808c9",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "f5fc5c41-0e64-4699-8272-b79ff0ab59ce",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "32f62251-e14e-4924-9aa3-2f029a535d73",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "56ac7976-8a78-47f6-8098-750b35d52803",
   "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[data-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 !important;\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-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\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",
       ".xr-index-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",
       ".xr-index-data-in:checked ~ .xr-index-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-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-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",
       ".xr-no-icon {\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; Size: 15kB\n",
       "Dimensions:       (x: 22, y: 21)\n",
       "Coordinates:\n",
       "  * x             (x) float64 176B -3.58e+05 -3.58e+05 ... -3.574e+05 -3.574e+05\n",
       "  * y             (y) float64 168B 3.55e+06 3.55e+06 ... 3.549e+06 3.549e+06\n",
       "    spatial_ref   int32 4B 0\n",
       "Data variables:\n",
       "    topo_dh       (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    elevation     (y, x) float32 2kB ...\n",
       "    topo_dh_fill  (y, x) float64 4kB ...\n",
       "    dxt_id        (y, x) float32 2kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    mb            (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "Attributes:\n",
       "    transform:           [ 3.00000000e+01  0.00000000e+00 -3.58000351e+05  0....\n",
       "    res:                 [30. 30.]\n",
       "    nodatavals:          nan\n",
       "    scales:              1.0\n",
       "    offsets:             0.0\n",
       "    topo process group:  g0018\n",
       "    effective samples:   10.5176\n",
       "    dh error:            3.6774\n",
       "    mb error:            0.0948</pre><div class='xr-wrap' style='display:none'><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-38205f22-93e7-4cda-bf54-4c6573daf972' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-38205f22-93e7-4cda-bf54-4c6573daf972' 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'>x</span>: 22</li><li><span class='xr-has-index'>y</span>: 21</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-2e693a62-9ad7-4ee3-90d6-51cafc1c8676' class='xr-section-summary-in' type='checkbox'  checked><label for='section-2e693a62-9ad7-4ee3-90d6-51cafc1c8676' 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'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-3.58e+05 -3.58e+05 ... -3.574e+05</div><input id='attrs-a2f435cc-ac1f-4725-95fa-3557ef371e30' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-a2f435cc-ac1f-4725-95fa-3557ef371e30' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-945a81d6-1cbc-444a-900d-4579644e2422' class='xr-var-data-in' type='checkbox'><label for='data-945a81d6-1cbc-444a-900d-4579644e2422' 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([-357985.350904, -357955.350904, -357925.350904, -357895.350904,\n",
       "       -357865.350904, -357835.350904, -357805.350904, -357775.350904,\n",
       "       -357745.350904, -357715.350904, -357685.350904, -357655.350904,\n",
       "       -357625.350904, -357595.350904, -357565.350904, -357535.350904,\n",
       "       -357505.350904, -357475.350904, -357445.350904, -357415.350904,\n",
       "       -357385.350904, -357355.350904])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>3.55e+06 3.55e+06 ... 3.549e+06</div><input id='attrs-824ee102-3ae1-40db-baa7-419b65d9b3a7' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-824ee102-3ae1-40db-baa7-419b65d9b3a7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3cba9125-a302-4468-a1f2-62235ef912f9' class='xr-var-data-in' type='checkbox'><label for='data-3cba9125-a302-4468-a1f2-62235ef912f9' 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([3550055.936844, 3550025.936844, 3549995.936844, 3549965.936844,\n",
       "       3549935.936844, 3549905.936844, 3549875.936844, 3549845.936844,\n",
       "       3549815.936844, 3549785.936844, 3549755.936844, 3549725.936844,\n",
       "       3549695.936844, 3549665.936844, 3549635.936844, 3549605.936844,\n",
       "       3549575.936844, 3549545.936844, 3549515.936844, 3549485.936844,\n",
       "       3549455.936844])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>spatial_ref</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>0</div><input id='attrs-aece518d-6a8f-43ed-8bdf-729e0717e422' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-aece518d-6a8f-43ed-8bdf-729e0717e422' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cf09e81b-dd29-47ce-813c-e0a09f139d8e' class='xr-var-data-in' type='checkbox'><label for='data-cf09e81b-dd29-47ce-813c-e0a09f139d8e' 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>crs_wkt :</span></dt><dd>PROJCS[&quot;Albers_Conic_Equal_Area&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Albers_Conic_Equal_Area&quot;],PARAMETER[&quot;latitude_of_center&quot;,0],PARAMETER[&quot;longitude_of_center&quot;,105],PARAMETER[&quot;standard_parallel_1&quot;,25],PARAMETER[&quot;standard_parallel_2&quot;,47],PARAMETER[&quot;false_easting&quot;,0],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH]]</dd><dt><span>semi_major_axis :</span></dt><dd>6378137.0</dd><dt><span>semi_minor_axis :</span></dt><dd>6356752.314245179</dd><dt><span>inverse_flattening :</span></dt><dd>298.257223563</dd><dt><span>reference_ellipsoid_name :</span></dt><dd>WGS 84</dd><dt><span>longitude_of_prime_meridian :</span></dt><dd>0.0</dd><dt><span>prime_meridian_name :</span></dt><dd>Greenwich</dd><dt><span>geographic_crs_name :</span></dt><dd>WGS 84</dd><dt><span>horizontal_datum_name :</span></dt><dd>World Geodetic System 1984</dd><dt><span>projected_crs_name :</span></dt><dd>Albers_Conic_Equal_Area</dd><dt><span>grid_mapping_name :</span></dt><dd>albers_conical_equal_area</dd><dt><span>standard_parallel :</span></dt><dd>(25.0, 47.0)</dd><dt><span>latitude_of_projection_origin :</span></dt><dd>0.0</dd><dt><span>longitude_of_central_meridian :</span></dt><dd>105.0</dd><dt><span>false_easting :</span></dt><dd>0.0</dd><dt><span>false_northing :</span></dt><dd>0.0</dd><dt><span>spatial_ref :</span></dt><dd>PROJCS[&quot;Albers_Conic_Equal_Area&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433,AUTHORITY[&quot;EPSG&quot;,&quot;9122&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Albers_Conic_Equal_Area&quot;],PARAMETER[&quot;latitude_of_center&quot;,0],PARAMETER[&quot;longitude_of_center&quot;,105],PARAMETER[&quot;standard_parallel_1&quot;,25],PARAMETER[&quot;standard_parallel_2&quot;,47],PARAMETER[&quot;false_easting&quot;,0],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AXIS[&quot;Easting&quot;,EAST],AXIS[&quot;Northing&quot;,NORTH]]</dd><dt><span>GeoTransform :</span></dt><dd>-358000.35090433015 30.0 0.0 3550070.936844331 0.0 -30.0</dd></dl></div><div class='xr-var-data'><pre>array(0)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-52cf78c2-3dc3-4106-ac87-c6ae2888b0bd' class='xr-section-summary-in' type='checkbox'  checked><label for='section-52cf78c2-3dc3-4106-ac87-c6ae2888b0bd' class='xr-section-summary' >Data variables: <span>(5)</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>topo_dh</span></div><div class='xr-var-dims'>(y, x)</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-6d378fa8-31d4-47d2-947b-ef9ebaedce9a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6d378fa8-31d4-47d2-947b-ef9ebaedce9a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-76e16d37-2663-4fea-ac22-e8bee81983dc' class='xr-var-data-in' type='checkbox'><label for='data-76e16d37-2663-4fea-ac22-e8bee81983dc' 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>long_name :</span></dt><dd>Elevation difference from Topo DEM and SRTM DEM (viod version, v1)</dd><dt><span>units :</span></dt><dd>m</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,       nan, ...,       nan,       nan,       nan],\n",
       "       ...,\n",
       "       [      nan,       nan, -0.132459, ...,       nan,       nan,       nan],\n",
       "       [      nan,       nan,       nan, ...,       nan,       nan,       nan],\n",
       "       [      nan,       nan,       nan, ...,       nan,       nan,       nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>elevation</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d613cab3-5b06-4c03-bf72-45ec91832917' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d613cab3-5b06-4c03-bf72-45ec91832917' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0770ea76-4738-4145-bd01-2998f8b2c7c8' class='xr-var-data-in' type='checkbox'><label for='data-0770ea76-4738-4145-bd01-2998f8b2c7c8' 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>long_name :</span></dt><dd>elevation from viod filled SRTM DEM</dd><dt><span>units :</span></dt><dd>m a.s.l</dd></dl></div><div class='xr-var-data'><pre>[462 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>topo_dh_fill</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f31802f4-b50d-4e05-99ca-203673220d3b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f31802f4-b50d-4e05-99ca-203673220d3b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8c29e760-7cbc-4916-a092-4ec559fae5d4' class='xr-var-data-in' type='checkbox'><label for='data-8c29e760-7cbc-4916-a092-4ec559fae5d4' 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>long_name :</span></dt><dd>Filled elevation difference from Topo DEM and SRTM DEM</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>Fill method :</span></dt><dd>Normalized regional hypsometric interpolation</dd></dl></div><div class='xr-var-data'><pre>[462 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>dxt_id</span></div><div class='xr-var-dims'>(y, x)</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-79494d67-5786-4fa0-9b57-6d89002f6353' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-79494d67-5786-4fa0-9b57-6d89002f6353' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c92a512e-787a-4d83-8e98-0da0073955a5' class='xr-var-data-in' type='checkbox'><label for='data-c92a512e-787a-4d83-8e98-0da0073955a5' 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>long_name :</span></dt><dd>Photography time of aerial photos (Topo maps time)</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,   nan, ..., 1967.,   nan,   nan],\n",
       "       ...,\n",
       "       [  nan, 1967., 1967., ...,   nan,   nan,   nan],\n",
       "       [  nan,   nan,   nan, ...,   nan,   nan,   nan],\n",
       "       [  nan,   nan,   nan, ...,   nan,   nan,   nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>mb</span></div><div class='xr-var-dims'>(y, x)</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-c42c481d-a8bb-4619-98ad-644c7bbe0257' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c42c481d-a8bb-4619-98ad-644c7bbe0257' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-93db5d99-8465-4d46-99c7-391fecdf0a51' class='xr-var-data-in' type='checkbox'><label for='data-93db5d99-8465-4d46-99c7-391fecdf0a51' 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>long_name :</span></dt><dd>Mass balance (by using [topo_dh])</dd><dt><span>units :</span></dt><dd>m w.e. yr-1</dd></dl></div><div class='xr-var-data'><pre>array([[            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "        -2.70498599e-01, -1.10202002e-01,  1.25243658e-01,\n",
       "         1.73880409e-01,  8.43914947e-02,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan, -1.95761809e-01,\n",
       "        -2.19226340e-01, -1.17418400e-01, -2.69424137e-02,\n",
       "         4.99571075e-02,  1.27690200e-01,  1.03366719e-01,\n",
       "         3.19681163e-01,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "         1.49751761e-02, -1.18884849e-01, -1.81609660e-01,\n",
       "        -1.06236974e-01, -5.94100931e-02, -1.50055690e-01,\n",
       "...\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "        -1.07826075e-01, -3.06106724e-01,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan]])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-6b3a5871-8f80-410e-a72c-f3f402b26411' class='xr-section-summary-in' type='checkbox'  ><label for='section-6b3a5871-8f80-410e-a72c-f3f402b26411' class='xr-section-summary' >Indexes: <span>(2)</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-index-name'><div>x</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-bf197fbc-9a96-4242-a081-2505d6299724' class='xr-index-data-in' type='checkbox'/><label for='index-bf197fbc-9a96-4242-a081-2505d6299724' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([-357985.35090433015, -357955.35090433015, -357925.35090433015,\n",
       "       -357895.35090433015, -357865.35090433015, -357835.35090433015,\n",
       "       -357805.35090433015, -357775.35090433015, -357745.35090433015,\n",
       "       -357715.35090433015, -357685.35090433015, -357655.35090433015,\n",
       "       -357625.35090433015, -357595.35090433015, -357565.35090433015,\n",
       "       -357535.35090433015, -357505.35090433015, -357475.35090433015,\n",
       "       -357445.35090433015, -357415.35090433015, -357385.35090433015,\n",
       "       -357355.35090433015],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;x&#x27;))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>y</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-3ced869a-96c9-494d-b245-05357ff3357d' class='xr-index-data-in' type='checkbox'/><label for='index-3ced869a-96c9-494d-b245-05357ff3357d' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([3550055.936844331, 3550025.936844331, 3549995.936844331,\n",
       "       3549965.936844331, 3549935.936844331, 3549905.936844331,\n",
       "       3549875.936844331, 3549845.936844331, 3549815.936844331,\n",
       "       3549785.936844331, 3549755.936844331, 3549725.936844331,\n",
       "       3549695.936844331, 3549665.936844331, 3549635.936844331,\n",
       "       3549605.936844331, 3549575.936844331, 3549545.936844331,\n",
       "       3549515.936844331, 3549485.936844331, 3549455.936844331],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;y&#x27;))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-cee69176-3d54-4d91-ac06-99b99ca06f68' class='xr-section-summary-in' type='checkbox'  checked><label for='section-cee69176-3d54-4d91-ac06-99b99ca06f68' class='xr-section-summary' >Attributes: <span>(9)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>transform :</span></dt><dd>[ 3.00000000e+01  0.00000000e+00 -3.58000351e+05  0.00000000e+00\n",
       " -3.00000000e+01  3.55007094e+06]</dd><dt><span>res :</span></dt><dd>[30. 30.]</dd><dt><span>nodatavals :</span></dt><dd>nan</dd><dt><span>scales :</span></dt><dd>1.0</dd><dt><span>offsets :</span></dt><dd>0.0</dd><dt><span>topo process group :</span></dt><dd>g0018</dd><dt><span>effective samples :</span></dt><dd>10.5176</dd><dt><span>dh error :</span></dt><dd>3.6774</dd><dt><span>mb error :</span></dt><dd>0.0948</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset> Size: 15kB\n",
       "Dimensions:       (x: 22, y: 21)\n",
       "Coordinates:\n",
       "  * x             (x) float64 176B -3.58e+05 -3.58e+05 ... -3.574e+05 -3.574e+05\n",
       "  * y             (y) float64 168B 3.55e+06 3.55e+06 ... 3.549e+06 3.549e+06\n",
       "    spatial_ref   int32 4B 0\n",
       "Data variables:\n",
       "    topo_dh       (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    elevation     (y, x) float32 2kB ...\n",
       "    topo_dh_fill  (y, x) float64 4kB ...\n",
       "    dxt_id        (y, x) float32 2kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    mb            (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "Attributes:\n",
       "    transform:           [ 3.00000000e+01  0.00000000e+00 -3.58000351e+05  0....\n",
       "    res:                 [30. 30.]\n",
       "    nodatavals:          nan\n",
       "    scales:              1.0\n",
       "    offsets:             0.0\n",
       "    topo process group:  g0018\n",
       "    effective samples:   10.5176\n",
       "    dh error:            3.6774\n",
       "    mb error:            0.0948"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds_gla0"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 42,
   "id": "4fbf6557-a839-40f3-8507-22b39f2c5001",
   "metadata": {},
   "outputs": [],
   "source": [
    "dff_tst=pd.DataFrame(np.array(aaa),columns=['mb','mb_err'])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 47,
   "id": "5b198a74-ee64-4b15-805d-ab3a3d136aea",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(-0.31040319417826173, 0.1228385022605827)"
      ]
     },
     "execution_count": 47,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dff_tst['mb'].mean(),dff_tst['mb_err'].median()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 28,
   "id": "c2d41a72-bd47-442b-aea4-5fb52b97923c",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "array(0.09480903)"
      ]
     },
     "execution_count": 28,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "mb_err.values"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 31,
   "id": "8828503e-683a-4cc4-b94a-8a0542eed843",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": 37,
   "id": "cdb0ea33-2f94-455e-b098-907b13bad524",
   "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[data-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 !important;\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-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\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",
       ".xr-index-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",
       ".xr-index-data-in:checked ~ .xr-index-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-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-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",
       ".xr-no-icon {\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; Size: 15kB\n",
       "Dimensions:       (x: 22, y: 21)\n",
       "Coordinates:\n",
       "  * x             (x) float64 176B -3.58e+05 -3.58e+05 ... -3.574e+05 -3.574e+05\n",
       "  * y             (y) float64 168B 3.55e+06 3.55e+06 ... 3.549e+06 3.549e+06\n",
       "Data variables:\n",
       "    topo_dh       (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    elevation     (y, x) float32 2kB ...\n",
       "    topo_dh_fill  (y, x) float64 4kB ...\n",
       "    dxt_id        (y, x) float32 2kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    mb            (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "Attributes:\n",
       "    transform:           [ 3.00000000e+01  0.00000000e+00 -3.58000351e+05  0....\n",
       "    crs:                 +proj=aea +lat_0=0 +lon_0=105 +lat_1=25 +lat_2=47 +x...\n",
       "    res:                 [30. 30.]\n",
       "    nodatavals:          nan\n",
       "    scales:              1.0\n",
       "    offsets:             0.0\n",
       "    topo process group:  g0018\n",
       "    effective samples:   10.5176\n",
       "    dh error:            3.6774\n",
       "    mb error:            0.0948</pre><div class='xr-wrap' style='display:none'><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-c8802c0b-f56f-4605-bd41-10bf3a2f2d11' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-c8802c0b-f56f-4605-bd41-10bf3a2f2d11' 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'>x</span>: 22</li><li><span class='xr-has-index'>y</span>: 21</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-19688216-0932-4fcc-a5ec-148dad122df1' class='xr-section-summary-in' type='checkbox'  checked><label for='section-19688216-0932-4fcc-a5ec-148dad122df1' class='xr-section-summary' >Coordinates: <span>(2)</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'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-3.58e+05 -3.58e+05 ... -3.574e+05</div><input id='attrs-1b665e4a-2036-49dc-9c2c-10f217e89fb8' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1b665e4a-2036-49dc-9c2c-10f217e89fb8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6bf251fa-5cf5-44f2-85fc-abc069ada38e' class='xr-var-data-in' type='checkbox'><label for='data-6bf251fa-5cf5-44f2-85fc-abc069ada38e' 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([-357985.350904, -357955.350904, -357925.350904, -357895.350904,\n",
       "       -357865.350904, -357835.350904, -357805.350904, -357775.350904,\n",
       "       -357745.350904, -357715.350904, -357685.350904, -357655.350904,\n",
       "       -357625.350904, -357595.350904, -357565.350904, -357535.350904,\n",
       "       -357505.350904, -357475.350904, -357445.350904, -357415.350904,\n",
       "       -357385.350904, -357355.350904])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>3.55e+06 3.55e+06 ... 3.549e+06</div><input id='attrs-21b27d91-9ecf-47d8-924d-650f89433081' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-21b27d91-9ecf-47d8-924d-650f89433081' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ca26af5b-0d24-434a-aec6-b092c9236999' class='xr-var-data-in' type='checkbox'><label for='data-ca26af5b-0d24-434a-aec6-b092c9236999' 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([3550055.936844, 3550025.936844, 3549995.936844, 3549965.936844,\n",
       "       3549935.936844, 3549905.936844, 3549875.936844, 3549845.936844,\n",
       "       3549815.936844, 3549785.936844, 3549755.936844, 3549725.936844,\n",
       "       3549695.936844, 3549665.936844, 3549635.936844, 3549605.936844,\n",
       "       3549575.936844, 3549545.936844, 3549515.936844, 3549485.936844,\n",
       "       3549455.936844])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-929ad103-c398-432e-b9dc-ccb341a0ef39' class='xr-section-summary-in' type='checkbox'  checked><label for='section-929ad103-c398-432e-b9dc-ccb341a0ef39' class='xr-section-summary' >Data variables: <span>(5)</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>topo_dh</span></div><div class='xr-var-dims'>(y, x)</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-4f56fa35-1009-4cf5-894d-52426f5595d3' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4f56fa35-1009-4cf5-894d-52426f5595d3' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-085eaf79-79cc-4683-a372-4267311ed4b6' class='xr-var-data-in' type='checkbox'><label for='data-085eaf79-79cc-4683-a372-4267311ed4b6' 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>long_name :</span></dt><dd>Elevation difference from Topo DEM and SRTM DEM (viod version, v1)</dd><dt><span>units :</span></dt><dd>m</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,       nan, ...,       nan,       nan,       nan],\n",
       "       ...,\n",
       "       [      nan,       nan, -0.132459, ...,       nan,       nan,       nan],\n",
       "       [      nan,       nan,       nan, ...,       nan,       nan,       nan],\n",
       "       [      nan,       nan,       nan, ...,       nan,       nan,       nan]])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>elevation</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-09c48b72-9c8d-458d-8bd4-eb4266a5631b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-09c48b72-9c8d-458d-8bd4-eb4266a5631b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-aaf86127-7976-482d-91c6-e35f5a7e99e7' class='xr-var-data-in' type='checkbox'><label for='data-aaf86127-7976-482d-91c6-e35f5a7e99e7' 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>long_name :</span></dt><dd>elevation from viod filled SRTM DEM</dd><dt><span>units :</span></dt><dd>m a.s.l</dd></dl></div><div class='xr-var-data'><pre>[462 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>topo_dh_fill</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-39780930-457d-4dd1-92b4-09a2e1109e6b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-39780930-457d-4dd1-92b4-09a2e1109e6b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cba4c7be-d029-469e-8e62-9e926207636b' class='xr-var-data-in' type='checkbox'><label for='data-cba4c7be-d029-469e-8e62-9e926207636b' 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>long_name :</span></dt><dd>Filled elevation difference from Topo DEM and SRTM DEM</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>Fill method :</span></dt><dd>Normalized regional hypsometric interpolation</dd></dl></div><div class='xr-var-data'><pre>[462 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>dxt_id</span></div><div class='xr-var-dims'>(y, x)</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-d01c9d62-45a0-4936-a519-0e96cfb86e90' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d01c9d62-45a0-4936-a519-0e96cfb86e90' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-05b80e97-8f6a-454a-832d-41002aa3d8ac' class='xr-var-data-in' type='checkbox'><label for='data-05b80e97-8f6a-454a-832d-41002aa3d8ac' 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>long_name :</span></dt><dd>Photography time of aerial photos (Topo maps time)</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,   nan, ..., 1967.,   nan,   nan],\n",
       "       ...,\n",
       "       [  nan, 1967., 1967., ...,   nan,   nan,   nan],\n",
       "       [  nan,   nan,   nan, ...,   nan,   nan,   nan],\n",
       "       [  nan,   nan,   nan, ...,   nan,   nan,   nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>mb</span></div><div class='xr-var-dims'>(y, x)</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-c2faee10-a634-43a2-b694-9624f735628b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c2faee10-a634-43a2-b694-9624f735628b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e15a160c-f63f-4140-8570-bb1b5c95b078' class='xr-var-data-in' type='checkbox'><label for='data-e15a160c-f63f-4140-8570-bb1b5c95b078' 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>long_name :</span></dt><dd>Mass balance (by using [topo_dh])</dd><dt><span>units :</span></dt><dd>m w.e. yr-1</dd></dl></div><div class='xr-var-data'><pre>array([[            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "        -2.70498599e-01, -1.10202002e-01,  1.25243658e-01,\n",
       "         1.73880409e-01,  8.43914947e-02,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan, -1.95761809e-01,\n",
       "        -2.19226340e-01, -1.17418400e-01, -2.69424137e-02,\n",
       "         4.99571075e-02,  1.27690200e-01,  1.03366719e-01,\n",
       "         3.19681163e-01,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "         1.49751761e-02, -1.18884849e-01, -1.81609660e-01,\n",
       "        -1.06236974e-01, -5.94100931e-02, -1.50055690e-01,\n",
       "...\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "        -1.07826075e-01, -3.06106724e-01,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan],\n",
       "       [            nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan,             nan,             nan,\n",
       "                    nan]])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-c646aa99-79d7-425c-8356-a3770acd93d2' class='xr-section-summary-in' type='checkbox'  ><label for='section-c646aa99-79d7-425c-8356-a3770acd93d2' class='xr-section-summary' >Indexes: <span>(2)</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-index-name'><div>x</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-5e5750bd-6b6b-419a-b363-db6d65a052c2' class='xr-index-data-in' type='checkbox'/><label for='index-5e5750bd-6b6b-419a-b363-db6d65a052c2' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([-357985.35090433015, -357955.35090433015, -357925.35090433015,\n",
       "       -357895.35090433015, -357865.35090433015, -357835.35090433015,\n",
       "       -357805.35090433015, -357775.35090433015, -357745.35090433015,\n",
       "       -357715.35090433015, -357685.35090433015, -357655.35090433015,\n",
       "       -357625.35090433015, -357595.35090433015, -357565.35090433015,\n",
       "       -357535.35090433015, -357505.35090433015, -357475.35090433015,\n",
       "       -357445.35090433015, -357415.35090433015, -357385.35090433015,\n",
       "       -357355.35090433015],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;x&#x27;))</pre></div></li><li class='xr-var-item'><div class='xr-index-name'><div>y</div></div><div class='xr-index-preview'>PandasIndex</div><div></div><input id='index-8d7511e4-ef67-4af4-89ca-e91744f2f6d3' class='xr-index-data-in' type='checkbox'/><label for='index-8d7511e4-ef67-4af4-89ca-e91744f2f6d3' title='Show/Hide index repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-index-data'><pre>PandasIndex(Index([3550055.936844331, 3550025.936844331, 3549995.936844331,\n",
       "       3549965.936844331, 3549935.936844331, 3549905.936844331,\n",
       "       3549875.936844331, 3549845.936844331, 3549815.936844331,\n",
       "       3549785.936844331, 3549755.936844331, 3549725.936844331,\n",
       "       3549695.936844331, 3549665.936844331, 3549635.936844331,\n",
       "       3549605.936844331, 3549575.936844331, 3549545.936844331,\n",
       "       3549515.936844331, 3549485.936844331, 3549455.936844331],\n",
       "      dtype=&#x27;float64&#x27;, name=&#x27;y&#x27;))</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-f7ad61ba-d80b-4e3d-9865-177a8b118035' class='xr-section-summary-in' type='checkbox'  ><label for='section-f7ad61ba-d80b-4e3d-9865-177a8b118035' class='xr-section-summary' >Attributes: <span>(10)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>transform :</span></dt><dd>[ 3.00000000e+01  0.00000000e+00 -3.58000351e+05  0.00000000e+00\n",
       " -3.00000000e+01  3.55007094e+06]</dd><dt><span>crs :</span></dt><dd>+proj=aea +lat_0=0 +lon_0=105 +lat_1=25 +lat_2=47 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs=True</dd><dt><span>res :</span></dt><dd>[30. 30.]</dd><dt><span>nodatavals :</span></dt><dd>nan</dd><dt><span>scales :</span></dt><dd>1.0</dd><dt><span>offsets :</span></dt><dd>0.0</dd><dt><span>topo process group :</span></dt><dd>g0018</dd><dt><span>effective samples :</span></dt><dd>10.5176</dd><dt><span>dh error :</span></dt><dd>3.6774</dd><dt><span>mb error :</span></dt><dd>0.0948</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset> Size: 15kB\n",
       "Dimensions:       (x: 22, y: 21)\n",
       "Coordinates:\n",
       "  * x             (x) float64 176B -3.58e+05 -3.58e+05 ... -3.574e+05 -3.574e+05\n",
       "  * y             (y) float64 168B 3.55e+06 3.55e+06 ... 3.549e+06 3.549e+06\n",
       "Data variables:\n",
       "    topo_dh       (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    elevation     (y, x) float32 2kB ...\n",
       "    topo_dh_fill  (y, x) float64 4kB ...\n",
       "    dxt_id        (y, x) float32 2kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "    mb            (y, x) float64 4kB nan nan nan nan nan ... nan nan nan nan nan\n",
       "Attributes:\n",
       "    transform:           [ 3.00000000e+01  0.00000000e+00 -3.58000351e+05  0....\n",
       "    crs:                 +proj=aea +lat_0=0 +lon_0=105 +lat_1=25 +lat_2=47 +x...\n",
       "    res:                 [30. 30.]\n",
       "    nodatavals:          nan\n",
       "    scales:              1.0\n",
       "    offsets:             0.0\n",
       "    topo process group:  g0018\n",
       "    effective samples:   10.5176\n",
       "    dh error:            3.6774\n",
       "    mb error:            0.0948"
      ]
     },
     "execution_count": 37,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds_gla"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "941d549a-5832-4a9c-86d7-ef5274c36ad7",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "py39",
   "language": "python",
   "name": "py39"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.19"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}