{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# MIKE IO 1D Pandas DataFrame examples\n",
    "\n",
    "Results in MIKE IO 1D are fundamentally extracted as Pandas DataFrames. This notebook shows\n",
    "specific methods for working with those DataFrames."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from mikeio1d import Res1D\n",
    "\n",
    "res = Res1D(\"../tests/testdata/network_river.res1d\")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### Column modes\n",
    "MIKE IO 1D supports different types of column headers. Depending on your use case, one or the other may be preferable. There are four main modes:\n",
    "- string (default)\n",
    "- all\n",
    "- compact\n",
    "- timeseries\n",
    "\n",
    "These are illustrated in the following cells."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "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>Discharge:basin_left1:1.002</th>\n",
       "      <th>Discharge:basin_left1:17.889</th>\n",
       "      <th>Discharge:basin_left1:57.3595</th>\n",
       "      <th>Discharge:basin_left1:101.493</th>\n",
       "      <th>Discharge:basin_left1:144.075</th>\n",
       "      <th>Discharge:basin_left1:175.496</th>\n",
       "      <th>Discharge:basin_left1:197.549</th>\n",
       "      <th>Discharge:basin_left1:228.411</th>\n",
       "      <th>Discharge:basin_left1:255.612</th>\n",
       "      <th>Discharge:basin_left1:274.626</th>\n",
       "      <th>...</th>\n",
       "      <th>Discharge:tributary:95</th>\n",
       "      <th>Discharge:tributary:110</th>\n",
       "      <th>Discharge:tributary:125</th>\n",
       "      <th>Discharge:tributary:140</th>\n",
       "      <th>Discharge:tributary:175</th>\n",
       "      <th>Discharge:tributary:225</th>\n",
       "      <th>Discharge:tributary:275</th>\n",
       "      <th>Discharge:tributary:350</th>\n",
       "      <th>Discharge:tributary:425</th>\n",
       "      <th>Discharge:tributary:475</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.013857</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>0.050522</td>\n",
       "      <td>0.053607</td>\n",
       "      <td>0.090468</td>\n",
       "      <td>0.107519</td>\n",
       "      <td>0.144026</td>\n",
       "      <td>0.177496</td>\n",
       "      <td>0.191049</td>\n",
       "      <td>0.197755</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 101 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                     Discharge:basin_left1:1.002  \\\n",
       "2000-02-18 00:06:00                          0.0   \n",
       "2000-02-18 00:16:00                         -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:17.889  \\\n",
       "2000-02-18 00:06:00                           0.0   \n",
       "2000-02-18 00:16:00                          -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:57.3595  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:101.493  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:144.075  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:175.496  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:197.549  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:228.411  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:255.612  \\\n",
       "2000-02-18 00:06:00                            0.0   \n",
       "2000-02-18 00:16:00                           -0.0   \n",
       "\n",
       "                     Discharge:basin_left1:274.626  ...  \\\n",
       "2000-02-18 00:06:00                            0.0  ...   \n",
       "2000-02-18 00:16:00                           -0.0  ...   \n",
       "\n",
       "                     Discharge:tributary:95  Discharge:tributary:110  \\\n",
       "2000-02-18 00:06:00                0.000000                      0.0   \n",
       "2000-02-18 00:16:00                0.013857                     -0.0   \n",
       "\n",
       "                     Discharge:tributary:125  Discharge:tributary:140  \\\n",
       "2000-02-18 00:06:00                 0.000000                 0.000000   \n",
       "2000-02-18 00:16:00                 0.050522                 0.053607   \n",
       "\n",
       "                     Discharge:tributary:175  Discharge:tributary:225  \\\n",
       "2000-02-18 00:06:00                 0.000000                 0.000000   \n",
       "2000-02-18 00:16:00                 0.090468                 0.107519   \n",
       "\n",
       "                     Discharge:tributary:275  Discharge:tributary:350  \\\n",
       "2000-02-18 00:06:00                 0.000000                 0.000000   \n",
       "2000-02-18 00:16:00                 0.144026                 0.177496   \n",
       "\n",
       "                     Discharge:tributary:425  Discharge:tributary:475  \n",
       "2000-02-18 00:06:00                 0.000000                 0.000000  \n",
       "2000-02-18 00:16:00                 0.191049                 0.197755  \n",
       "\n",
       "[2 rows x 101 columns]"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# The default column mode concatenates information about the column with a ':' delimiter.\n",
    "df = res.reaches.Discharge.read()\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"21\" halign=\"left\">Discharge</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"21\" halign=\"left\">Reach</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th colspan=\"10\" halign=\"left\">basin_left1</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"10\" halign=\"left\">tributary</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>1.0020</th>\n",
       "      <th>17.8890</th>\n",
       "      <th>57.3595</th>\n",
       "      <th>101.4935</th>\n",
       "      <th>144.0745</th>\n",
       "      <th>175.4965</th>\n",
       "      <th>197.5490</th>\n",
       "      <th>228.4110</th>\n",
       "      <th>255.6120</th>\n",
       "      <th>274.6260</th>\n",
       "      <th>...</th>\n",
       "      <th>95.0000</th>\n",
       "      <th>110.0000</th>\n",
       "      <th>125.0000</th>\n",
       "      <th>140.0000</th>\n",
       "      <th>175.0000</th>\n",
       "      <th>225.0000</th>\n",
       "      <th>275.0000</th>\n",
       "      <th>350.0000</th>\n",
       "      <th>425.0000</th>\n",
       "      <th>475.0000</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>...</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.013857</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>0.050522</td>\n",
       "      <td>0.053607</td>\n",
       "      <td>0.090468</td>\n",
       "      <td>0.107519</td>\n",
       "      <td>0.144026</td>\n",
       "      <td>0.177496</td>\n",
       "      <td>0.191049</td>\n",
       "      <td>0.197755</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 101 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity              Discharge                                          \\\n",
       "group                     Reach                                           \n",
       "name                basin_left1                                           \n",
       "chainage               1.0020    17.8890   57.3595   101.4935  144.0745   \n",
       "tag                   0.0-246.6 0.0-246.6 0.0-246.6 0.0-246.6 0.0-246.6   \n",
       "2000-02-18 00:06:00         0.0       0.0       0.0       0.0       0.0   \n",
       "2000-02-18 00:16:00        -0.0      -0.0      -0.0      -0.0      -0.0   \n",
       "\n",
       "quantity                                                                   \\\n",
       "group                                                                       \n",
       "name                                                                        \n",
       "chainage             175.4965  197.5490  228.4110    255.6120    274.6260   \n",
       "tag                 0.0-246.6 0.0-246.6 0.0-246.6 246.6-413.6 246.6-413.6   \n",
       "2000-02-18 00:06:00       0.0       0.0       0.0         0.0         0.0   \n",
       "2000-02-18 00:16:00      -0.0      -0.0      -0.0        -0.0        -0.0   \n",
       "\n",
       "quantity             ...                                              \\\n",
       "group                ...                                               \n",
       "name                 ...  tributary                                    \n",
       "chainage             ...   95.0000    110.0000   125.0000   140.0000   \n",
       "tag                  ... 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "2000-02-18 00:06:00  ...   0.000000        0.0   0.000000   0.000000   \n",
       "2000-02-18 00:16:00  ...   0.013857       -0.0   0.050522   0.053607   \n",
       "\n",
       "quantity                                                                    \\\n",
       "group                                                                        \n",
       "name                                                                         \n",
       "chainage              175.0000   225.0000   275.0000   350.0000   425.0000   \n",
       "tag                 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "2000-02-18 00:06:00   0.000000   0.000000   0.000000   0.000000   0.000000   \n",
       "2000-02-18 00:16:00   0.090468   0.107519   0.144026   0.177496   0.191049   \n",
       "\n",
       "quantity                        \n",
       "group                           \n",
       "name                            \n",
       "chainage              475.0000  \n",
       "tag                 50.0-500.0  \n",
       "2000-02-18 00:06:00   0.000000  \n",
       "2000-02-18 00:16:00   0.197755  \n",
       "\n",
       "[2 rows x 101 columns]"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# A hierarchical column mode is also supported. Only relevant levels are included with 'compact'. For full hierarchy use 'all'.\n",
    "df = res.reaches.Discharge.read(column_mode=\"compact\")\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "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>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=1.002, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=17.889, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=57.3595, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=101.4935, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=144.0745, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=175.4965, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=197.54899999999998, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=228.411, tag='0.0-246.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=255.612, tag='246.6-413.6', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=274.626, tag='246.6-413.6', duplicate=0, derived=False)</th>\n",
       "      <th>...</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=95.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=110.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=125.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=140.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=175.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=225.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=275.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=350.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=425.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "      <th>TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=475.0, tag='50.0-500.0', duplicate=0, derived=False)</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.013857</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>0.050522</td>\n",
       "      <td>0.053607</td>\n",
       "      <td>0.090468</td>\n",
       "      <td>0.107519</td>\n",
       "      <td>0.144026</td>\n",
       "      <td>0.177496</td>\n",
       "      <td>0.191049</td>\n",
       "      <td>0.197755</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 101 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=1.002, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                    \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                    \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=17.889, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                     \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                     \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=57.3595, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                      \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                      \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=101.4935, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                       \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                       \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=144.0745, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                       \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                       \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=175.4965, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                       \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                       \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=197.54899999999998, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                                 \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                                 \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=228.411, tag='0.0-246.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                      \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                      \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=255.612, tag='246.6-413.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                        \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                        \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='basin_left1', chainage=274.626, tag='246.6-413.6', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                        \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                        \n",
       "\n",
       "                     ...  \\\n",
       "2000-02-18 00:06:00  ...   \n",
       "2000-02-18 00:16:00  ...   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=95.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                  \n",
       "2000-02-18 00:16:00                                           0.013857                                                                                  \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=110.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                                0.0                                                                                   \n",
       "2000-02-18 00:16:00                                               -0.0                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=125.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.050522                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=140.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.053607                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=175.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.090468                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=225.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.107519                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=275.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.144026                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=350.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.177496                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=425.0, tag='50.0-500.0', duplicate=0, derived=False)  \\\n",
       "2000-02-18 00:06:00                                           0.000000                                                                                   \n",
       "2000-02-18 00:16:00                                           0.191049                                                                                   \n",
       "\n",
       "                     TimeSeriesId(quantity='Discharge', group='Reach', name='tributary', chainage=475.0, tag='50.0-500.0', duplicate=0, derived=False)  \n",
       "2000-02-18 00:06:00                                           0.000000                                                                                  \n",
       "2000-02-18 00:16:00                                           0.197755                                                                                  \n",
       "\n",
       "[2 rows x 101 columns]"
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# It's also possible to have a TimeSeriesId column index.\n",
    "# TimeSeriesId is an object that uniquely identifies each time series.\n",
    "df = res.reaches.Discharge.read(column_mode=\"timeseries\")\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "### MIKE IO 1D extends Pandas by providing a '.m1d' accessor on all DataFrames.\n",
    "The accessor provides several useful methods for working with DataFrames within MIKE IO 1D:\n",
    "- .m1d.query()\n",
    "- .m1d.compact()\n",
    "- .m1d.groupby()\n",
    "- .m1d.groupby_chainage()\n",
    "- .m1d.agg_chainage()\n",
    "\n",
    "These methods are illustrated below."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<mikeio1d.pandas_extension.mikeio1d_accessor.Mikeio1dAccessor at 0x7ebdf5465d30>"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# To use the .m1d accessor, the DataFrame must have a MultiIndex column (e.g. column_mode='all' or 'compact').\n",
    "# The .m1d accessor exists on the DataFrame itself.\n",
    "df = res.read(column_mode=\"all\")\n",
    "df.m1d"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### .m1d.compact()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"10\" halign=\"left\">WaterLevel</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"6\" halign=\"left\">FlowVelocity</th>\n",
       "      <th>Variable:TwoTimeSensorGateLevel</th>\n",
       "      <th>Water level:Sensor:s.h.river53745.34</th>\n",
       "      <th>Gate level:Sensor:SensorGateLevel</th>\n",
       "      <th>Discharge:Sensor:SensorGauge1</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"10\" halign=\"left\">Node</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"6\" halign=\"left\">Reach</th>\n",
       "      <th>Global</th>\n",
       "      <th>Global</th>\n",
       "      <th>Global</th>\n",
       "      <th>Global</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th>'basin_left1', 0</th>\n",
       "      <th>'basin_left1', 483.74285236705026</th>\n",
       "      <th>'basin_left2', 0</th>\n",
       "      <th>'basin_left2', 441.71489580142429</th>\n",
       "      <th>'basin_right', -10</th>\n",
       "      <th>'basin_right', 720</th>\n",
       "      <th>'river', 53100</th>\n",
       "      <th>'river', 55124.2276598819</th>\n",
       "      <th>'tributary', 500</th>\n",
       "      <th>'basin_left1', 246.61</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"6\" halign=\"left\">tributary</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>...</th>\n",
       "      <th>175.00000</th>\n",
       "      <th>225.00000</th>\n",
       "      <th>275.00000</th>\n",
       "      <th>350.00000</th>\n",
       "      <th>425.00000</th>\n",
       "      <th>475.00000</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>...</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>duplicate</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>...</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "      <th>0</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>derived</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>...</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "      <th>False</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>57.426998</td>\n",
       "      <td>53.495998</td>\n",
       "      <td>54.002998</td>\n",
       "      <td>59.200001</td>\n",
       "      <td>53.652</td>\n",
       "      <td>58.342999</td>\n",
       "      <td>51.835999</td>\n",
       "      <td>53.194</td>\n",
       "      <td>55.480000</td>\n",
       "      <td>53.065304</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>51.905998</td>\n",
       "      <td>53.77</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>57.426998</td>\n",
       "      <td>53.495998</td>\n",
       "      <td>54.002998</td>\n",
       "      <td>59.200001</td>\n",
       "      <td>53.652</td>\n",
       "      <td>58.342999</td>\n",
       "      <td>51.836006</td>\n",
       "      <td>53.194</td>\n",
       "      <td>55.946224</td>\n",
       "      <td>53.065304</td>\n",
       "      <td>...</td>\n",
       "      <td>0.608204</td>\n",
       "      <td>0.517946</td>\n",
       "      <td>0.449898</td>\n",
       "      <td>0.446345</td>\n",
       "      <td>0.419857</td>\n",
       "      <td>0.417569</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>51.905998</td>\n",
       "      <td>53.77</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 498 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity                  WaterLevel                                    \\\n",
       "group                           Node                                     \n",
       "name                'basin_left1', 0 'basin_left1', 483.74285236705026   \n",
       "chainage                         NaN                               NaN   \n",
       "tag                                                                      \n",
       "duplicate                          0                                 0   \n",
       "derived                        False                             False   \n",
       "2000-02-18 00:06:00        57.426998                         53.495998   \n",
       "2000-02-18 00:16:00        57.426998                         53.495998   \n",
       "\n",
       "quantity                                                                \\\n",
       "group                                                                    \n",
       "name                'basin_left2', 0 'basin_left2', 441.71489580142429   \n",
       "chainage                         NaN                               NaN   \n",
       "tag                                                                      \n",
       "duplicate                          0                                 0   \n",
       "derived                        False                             False   \n",
       "2000-02-18 00:06:00        54.002998                         59.200001   \n",
       "2000-02-18 00:16:00        54.002998                         59.200001   \n",
       "\n",
       "quantity                                                                  \\\n",
       "group                                                                      \n",
       "name                'basin_right', -10 'basin_right', 720 'river', 53100   \n",
       "chainage                           NaN                NaN            NaN   \n",
       "tag                                                                        \n",
       "duplicate                            0                  0              0   \n",
       "derived                          False              False          False   \n",
       "2000-02-18 00:06:00             53.652          58.342999      51.835999   \n",
       "2000-02-18 00:16:00             53.652          58.342999      51.836006   \n",
       "\n",
       "quantity                                                        \\\n",
       "group                                                            \n",
       "name                'river', 55124.2276598819 'tributary', 500   \n",
       "chainage                                  NaN              NaN   \n",
       "tag                                                              \n",
       "duplicate                                   0                0   \n",
       "derived                                 False            False   \n",
       "2000-02-18 00:06:00                    53.194        55.480000   \n",
       "2000-02-18 00:16:00                    53.194        55.946224   \n",
       "\n",
       "quantity                                   ... FlowVelocity              \\\n",
       "group                                      ...        Reach               \n",
       "name                'basin_left1', 246.61  ...    tributary               \n",
       "chainage                              NaN  ...  175.00000   225.00000     \n",
       "tag                                        ...   50.0-500.0  50.0-500.0   \n",
       "duplicate                               0  ...            0           0   \n",
       "derived                             False  ...        False       False   \n",
       "2000-02-18 00:06:00             53.065304  ...     0.000000    0.000000   \n",
       "2000-02-18 00:16:00             53.065304  ...     0.608204    0.517946   \n",
       "\n",
       "quantity                                                             \\\n",
       "group                                                                 \n",
       "name                                                                  \n",
       "chainage            275.00000   350.00000   425.00000   475.00000     \n",
       "tag                  50.0-500.0  50.0-500.0  50.0-500.0  50.0-500.0   \n",
       "duplicate                     0           0           0           0   \n",
       "derived                   False       False       False       False   \n",
       "2000-02-18 00:06:00    0.000000    0.000000    0.000000    0.000000   \n",
       "2000-02-18 00:16:00    0.449898    0.446345    0.419857    0.417569   \n",
       "\n",
       "quantity            Variable:TwoTimeSensorGateLevel  \\\n",
       "group                                        Global   \n",
       "name                                                  \n",
       "chainage                                        NaN   \n",
       "tag                                                   \n",
       "duplicate                                         0   \n",
       "derived                                       False   \n",
       "2000-02-18 00:06:00                      107.540001   \n",
       "2000-02-18 00:16:00                      107.540001   \n",
       "\n",
       "quantity            Water level:Sensor:s.h.river53745.34  \\\n",
       "group                                             Global   \n",
       "name                                                       \n",
       "chainage                                             NaN   \n",
       "tag                                                        \n",
       "duplicate                                              0   \n",
       "derived                                            False   \n",
       "2000-02-18 00:06:00                            51.905998   \n",
       "2000-02-18 00:16:00                            51.905998   \n",
       "\n",
       "quantity            Gate level:Sensor:SensorGateLevel  \\\n",
       "group                                          Global   \n",
       "name                                                    \n",
       "chainage                                          NaN   \n",
       "tag                                                     \n",
       "duplicate                                           0   \n",
       "derived                                         False   \n",
       "2000-02-18 00:06:00                             53.77   \n",
       "2000-02-18 00:16:00                             53.77   \n",
       "\n",
       "quantity            Discharge:Sensor:SensorGauge1  \n",
       "group                                      Global  \n",
       "name                                               \n",
       "chainage                                      NaN  \n",
       "tag                                                \n",
       "duplicate                                       0  \n",
       "derived                                     False  \n",
       "2000-02-18 00:06:00                           0.0  \n",
       "2000-02-18 00:16:00                          -0.0  \n",
       "\n",
       "[2 rows x 498 columns]"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Let's make some queries on the DataFrame itself with help from the .m1d accessor.\n",
    "# Let's read the entire file into a DataFrame with column_mode='all'.\n",
    "df = res.read(column_mode=\"all\")\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"10\" halign=\"left\">WaterLevel</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"6\" halign=\"left\">FlowVelocity</th>\n",
       "      <th>Variable:TwoTimeSensorGateLevel</th>\n",
       "      <th>Water level:Sensor:s.h.river53745.34</th>\n",
       "      <th>Gate level:Sensor:SensorGateLevel</th>\n",
       "      <th>Discharge:Sensor:SensorGauge1</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"10\" halign=\"left\">Node</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"6\" halign=\"left\">Reach</th>\n",
       "      <th>Global</th>\n",
       "      <th>Global</th>\n",
       "      <th>Global</th>\n",
       "      <th>Global</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th>'basin_left1', 0</th>\n",
       "      <th>'basin_left1', 483.74285236705026</th>\n",
       "      <th>'basin_left2', 0</th>\n",
       "      <th>'basin_left2', 441.71489580142429</th>\n",
       "      <th>'basin_right', -10</th>\n",
       "      <th>'basin_right', 720</th>\n",
       "      <th>'river', 53100</th>\n",
       "      <th>'river', 55124.2276598819</th>\n",
       "      <th>'tributary', 500</th>\n",
       "      <th>'basin_left1', 246.61</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"6\" halign=\"left\">tributary</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>...</th>\n",
       "      <th>175.00000</th>\n",
       "      <th>225.00000</th>\n",
       "      <th>275.00000</th>\n",
       "      <th>350.00000</th>\n",
       "      <th>425.00000</th>\n",
       "      <th>475.00000</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "      <th>NaN</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>...</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>57.426998</td>\n",
       "      <td>53.495998</td>\n",
       "      <td>54.002998</td>\n",
       "      <td>59.200001</td>\n",
       "      <td>53.652</td>\n",
       "      <td>58.342999</td>\n",
       "      <td>51.835999</td>\n",
       "      <td>53.194</td>\n",
       "      <td>55.480000</td>\n",
       "      <td>53.065304</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>51.905998</td>\n",
       "      <td>53.77</td>\n",
       "      <td>0.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>57.426998</td>\n",
       "      <td>53.495998</td>\n",
       "      <td>54.002998</td>\n",
       "      <td>59.200001</td>\n",
       "      <td>53.652</td>\n",
       "      <td>58.342999</td>\n",
       "      <td>51.836006</td>\n",
       "      <td>53.194</td>\n",
       "      <td>55.946224</td>\n",
       "      <td>53.065304</td>\n",
       "      <td>...</td>\n",
       "      <td>0.608204</td>\n",
       "      <td>0.517946</td>\n",
       "      <td>0.449898</td>\n",
       "      <td>0.446345</td>\n",
       "      <td>0.419857</td>\n",
       "      <td>0.417569</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>51.905998</td>\n",
       "      <td>53.77</td>\n",
       "      <td>-0.0</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 498 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity                  WaterLevel                                    \\\n",
       "group                           Node                                     \n",
       "name                'basin_left1', 0 'basin_left1', 483.74285236705026   \n",
       "chainage                         NaN                               NaN   \n",
       "tag                                                                      \n",
       "2000-02-18 00:06:00        57.426998                         53.495998   \n",
       "2000-02-18 00:16:00        57.426998                         53.495998   \n",
       "\n",
       "quantity                                                                \\\n",
       "group                                                                    \n",
       "name                'basin_left2', 0 'basin_left2', 441.71489580142429   \n",
       "chainage                         NaN                               NaN   \n",
       "tag                                                                      \n",
       "2000-02-18 00:06:00        54.002998                         59.200001   \n",
       "2000-02-18 00:16:00        54.002998                         59.200001   \n",
       "\n",
       "quantity                                                                  \\\n",
       "group                                                                      \n",
       "name                'basin_right', -10 'basin_right', 720 'river', 53100   \n",
       "chainage                           NaN                NaN            NaN   \n",
       "tag                                                                        \n",
       "2000-02-18 00:06:00             53.652          58.342999      51.835999   \n",
       "2000-02-18 00:16:00             53.652          58.342999      51.836006   \n",
       "\n",
       "quantity                                                        \\\n",
       "group                                                            \n",
       "name                'river', 55124.2276598819 'tributary', 500   \n",
       "chainage                                  NaN              NaN   \n",
       "tag                                                              \n",
       "2000-02-18 00:06:00                    53.194        55.480000   \n",
       "2000-02-18 00:16:00                    53.194        55.946224   \n",
       "\n",
       "quantity                                   ... FlowVelocity              \\\n",
       "group                                      ...        Reach               \n",
       "name                'basin_left1', 246.61  ...    tributary               \n",
       "chainage                              NaN  ...  175.00000   225.00000     \n",
       "tag                                        ...   50.0-500.0  50.0-500.0   \n",
       "2000-02-18 00:06:00             53.065304  ...     0.000000    0.000000   \n",
       "2000-02-18 00:16:00             53.065304  ...     0.608204    0.517946   \n",
       "\n",
       "quantity                                                             \\\n",
       "group                                                                 \n",
       "name                                                                  \n",
       "chainage            275.00000   350.00000   425.00000   475.00000     \n",
       "tag                  50.0-500.0  50.0-500.0  50.0-500.0  50.0-500.0   \n",
       "2000-02-18 00:06:00    0.000000    0.000000    0.000000    0.000000   \n",
       "2000-02-18 00:16:00    0.449898    0.446345    0.419857    0.417569   \n",
       "\n",
       "quantity            Variable:TwoTimeSensorGateLevel  \\\n",
       "group                                        Global   \n",
       "name                                                  \n",
       "chainage                                        NaN   \n",
       "tag                                                   \n",
       "2000-02-18 00:06:00                      107.540001   \n",
       "2000-02-18 00:16:00                      107.540001   \n",
       "\n",
       "quantity            Water level:Sensor:s.h.river53745.34  \\\n",
       "group                                             Global   \n",
       "name                                                       \n",
       "chainage                                             NaN   \n",
       "tag                                                        \n",
       "2000-02-18 00:06:00                            51.905998   \n",
       "2000-02-18 00:16:00                            51.905998   \n",
       "\n",
       "quantity            Gate level:Sensor:SensorGateLevel  \\\n",
       "group                                          Global   \n",
       "name                                                    \n",
       "chainage                                          NaN   \n",
       "tag                                                     \n",
       "2000-02-18 00:06:00                             53.77   \n",
       "2000-02-18 00:16:00                             53.77   \n",
       "\n",
       "quantity            Discharge:Sensor:SensorGauge1  \n",
       "group                                      Global  \n",
       "name                                               \n",
       "chainage                                      NaN  \n",
       "tag                                                \n",
       "2000-02-18 00:06:00                           0.0  \n",
       "2000-02-18 00:16:00                          -0.0  \n",
       "\n",
       "[2 rows x 498 columns]"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# That's a lot of detail included and is a bit busy. We can use .m1d.compact() to remove redundant levels.\n",
    "df = df.m1d.compact()\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### .m1d.query()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"21\" halign=\"left\">Discharge</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"21\" halign=\"left\">Reach</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th colspan=\"10\" halign=\"left\">basin_left1</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"10\" halign=\"left\">tributary</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>1.0020</th>\n",
       "      <th>17.8890</th>\n",
       "      <th>57.3595</th>\n",
       "      <th>101.4935</th>\n",
       "      <th>144.0745</th>\n",
       "      <th>175.4965</th>\n",
       "      <th>197.5490</th>\n",
       "      <th>228.4110</th>\n",
       "      <th>255.6120</th>\n",
       "      <th>274.6260</th>\n",
       "      <th>...</th>\n",
       "      <th>95.0000</th>\n",
       "      <th>110.0000</th>\n",
       "      <th>125.0000</th>\n",
       "      <th>140.0000</th>\n",
       "      <th>175.0000</th>\n",
       "      <th>225.0000</th>\n",
       "      <th>275.0000</th>\n",
       "      <th>350.0000</th>\n",
       "      <th>425.0000</th>\n",
       "      <th>475.0000</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>...</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.013857</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>0.050522</td>\n",
       "      <td>0.053607</td>\n",
       "      <td>0.090468</td>\n",
       "      <td>0.107519</td>\n",
       "      <td>0.144026</td>\n",
       "      <td>0.177496</td>\n",
       "      <td>0.191049</td>\n",
       "      <td>0.197755</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>2 rows × 101 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity              Discharge                                          \\\n",
       "group                     Reach                                           \n",
       "name                basin_left1                                           \n",
       "chainage               1.0020    17.8890   57.3595   101.4935  144.0745   \n",
       "tag                   0.0-246.6 0.0-246.6 0.0-246.6 0.0-246.6 0.0-246.6   \n",
       "2000-02-18 00:06:00         0.0       0.0       0.0       0.0       0.0   \n",
       "2000-02-18 00:16:00        -0.0      -0.0      -0.0      -0.0      -0.0   \n",
       "\n",
       "quantity                                                                   \\\n",
       "group                                                                       \n",
       "name                                                                        \n",
       "chainage             175.4965  197.5490  228.4110    255.6120    274.6260   \n",
       "tag                 0.0-246.6 0.0-246.6 0.0-246.6 246.6-413.6 246.6-413.6   \n",
       "2000-02-18 00:06:00       0.0       0.0       0.0         0.0         0.0   \n",
       "2000-02-18 00:16:00      -0.0      -0.0      -0.0        -0.0        -0.0   \n",
       "\n",
       "quantity             ...                                              \\\n",
       "group                ...                                               \n",
       "name                 ...  tributary                                    \n",
       "chainage             ...   95.0000    110.0000   125.0000   140.0000   \n",
       "tag                  ... 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "2000-02-18 00:06:00  ...   0.000000        0.0   0.000000   0.000000   \n",
       "2000-02-18 00:16:00  ...   0.013857       -0.0   0.050522   0.053607   \n",
       "\n",
       "quantity                                                                    \\\n",
       "group                                                                        \n",
       "name                                                                         \n",
       "chainage              175.0000   225.0000   275.0000   350.0000   425.0000   \n",
       "tag                 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "2000-02-18 00:06:00   0.000000   0.000000   0.000000   0.000000   0.000000   \n",
       "2000-02-18 00:16:00   0.090468   0.107519   0.144026   0.177496   0.191049   \n",
       "\n",
       "quantity                        \n",
       "group                           \n",
       "name                            \n",
       "chainage              475.0000  \n",
       "tag                 50.0-500.0  \n",
       "2000-02-18 00:06:00   0.000000  \n",
       "2000-02-18 00:16:00   0.197755  \n",
       "\n",
       "[2 rows x 101 columns]"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Let's get Discharge for the reaches.\n",
    "df = df.m1d.query(\"group=='Reach' and quantity=='Discharge'\")\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"11\" halign=\"left\">Discharge</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"11\" halign=\"left\">Reach</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th colspan=\"11\" halign=\"left\">tributary</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>70.0</th>\n",
       "      <th>95.0</th>\n",
       "      <th>110.0</th>\n",
       "      <th>125.0</th>\n",
       "      <th>140.0</th>\n",
       "      <th>175.0</th>\n",
       "      <th>225.0</th>\n",
       "      <th>275.0</th>\n",
       "      <th>350.0</th>\n",
       "      <th>425.0</th>\n",
       "      <th>475.0</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:06:00</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2000-02-18 00:16:00</th>\n",
       "      <td>0.020551</td>\n",
       "      <td>0.013857</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>0.050522</td>\n",
       "      <td>0.053607</td>\n",
       "      <td>0.090468</td>\n",
       "      <td>0.107519</td>\n",
       "      <td>0.144026</td>\n",
       "      <td>0.177496</td>\n",
       "      <td>0.191049</td>\n",
       "      <td>0.197755</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity             Discharge                                              \\\n",
       "group                    Reach                                               \n",
       "name                 tributary                                               \n",
       "chainage                 70.0       95.0       110.0      125.0      140.0   \n",
       "tag                 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "2000-02-18 00:06:00   0.000000   0.000000        0.0   0.000000   0.000000   \n",
       "2000-02-18 00:16:00   0.020551   0.013857       -0.0   0.050522   0.053607   \n",
       "\n",
       "quantity                                                                    \\\n",
       "group                                                                        \n",
       "name                                                                         \n",
       "chainage                 175.0      225.0      275.0      350.0      425.0   \n",
       "tag                 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "2000-02-18 00:06:00   0.000000   0.000000   0.000000   0.000000   0.000000   \n",
       "2000-02-18 00:16:00   0.090468   0.107519   0.144026   0.177496   0.191049   \n",
       "\n",
       "quantity                        \n",
       "group                           \n",
       "name                            \n",
       "chainage                 475.0  \n",
       "tag                 50.0-500.0  \n",
       "2000-02-18 00:06:00   0.000000  \n",
       "2000-02-18 00:16:00   0.197755  "
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Let's look at all the reaches with 'trib' in their name.\n",
    "df = df.m1d.query(\"name.str.contains('trib')\")\n",
    "df.head(2)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"11\" halign=\"left\">Discharge</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"11\" halign=\"left\">Reach</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th colspan=\"11\" halign=\"left\">tributary</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>70.0</th>\n",
       "      <th>95.0</th>\n",
       "      <th>110.0</th>\n",
       "      <th>125.0</th>\n",
       "      <th>140.0</th>\n",
       "      <th>175.0</th>\n",
       "      <th>225.0</th>\n",
       "      <th>275.0</th>\n",
       "      <th>350.0</th>\n",
       "      <th>425.0</th>\n",
       "      <th>475.0</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>2.027786</td>\n",
       "      <td>1.698108</td>\n",
       "      <td>1.452274</td>\n",
       "      <td>1.198294</td>\n",
       "      <td>0.945816</td>\n",
       "      <td>0.696441</td>\n",
       "      <td>0.665936</td>\n",
       "      <td>0.59445</td>\n",
       "      <td>0.589273</td>\n",
       "      <td>0.495095</td>\n",
       "      <td>0.296315</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity  Discharge                                                         \\\n",
       "group         Reach                                                          \n",
       "name      tributary                                                          \n",
       "chainage      70.0       95.0       110.0      125.0      140.0      175.0   \n",
       "tag      50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0   \n",
       "max        2.027786   1.698108   1.452274   1.198294   0.945816   0.696441   \n",
       "\n",
       "quantity                                                         \n",
       "group                                                            \n",
       "name                                                             \n",
       "chainage      225.0      275.0      350.0      425.0      475.0  \n",
       "tag      50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0 50.0-500.0  \n",
       "max        0.665936    0.59445   0.589273   0.495095   0.296315  "
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Let's look at the max discharge for each reach.\n",
    "df.agg([\"max\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe thead tr:last-of-type th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>max</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>tag</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th>name</th>\n",
       "      <th>chainage</th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"11\" valign=\"top\">Discharge</th>\n",
       "      <th rowspan=\"11\" valign=\"top\">tributary</th>\n",
       "      <th>70.0</th>\n",
       "      <td>2.027786</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95.0</th>\n",
       "      <td>1.698108</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>110.0</th>\n",
       "      <td>1.452274</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>125.0</th>\n",
       "      <td>1.198294</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>140.0</th>\n",
       "      <td>0.945816</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>175.0</th>\n",
       "      <td>0.696441</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>225.0</th>\n",
       "      <td>0.665936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>275.0</th>\n",
       "      <td>0.594450</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>350.0</th>\n",
       "      <td>0.589273</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>425.0</th>\n",
       "      <td>0.495095</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>475.0</th>\n",
       "      <td>0.296315</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                    max\n",
       "tag                          50.0-500.0\n",
       "quantity  name      chainage           \n",
       "Discharge tributary 70.0       2.027786\n",
       "                    95.0       1.698108\n",
       "                    110.0      1.452274\n",
       "                    125.0      1.198294\n",
       "                    140.0      0.945816\n",
       "                    175.0      0.696441\n",
       "                    225.0      0.665936\n",
       "                    275.0      0.594450\n",
       "                    350.0      0.589273\n",
       "                    425.0      0.495095\n",
       "                    475.0      0.296315"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Use some standard pandas methods to format the table a different way\n",
    "# Tip: Chaining methods in brackets is a great way to explore the data. Comment out lines from bottom up to see the effect.\n",
    "(\n",
    "    df.agg([\"max\"]).T.droplevel([\"group\"]).unstack()\n",
    "    # .pipe(lambda df: df * 2) # Uncomment this line to see the effect of the .pipe() method.\n",
    ")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "\n",
       "    .dataframe thead tr:last-of-type th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>max</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th></th>\n",
       "      <th></th>\n",
       "      <th>tag</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th>name</th>\n",
       "      <th>chainage</th>\n",
       "      <th></th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th rowspan=\"11\" valign=\"top\">Discharge</th>\n",
       "      <th rowspan=\"11\" valign=\"top\">tributary</th>\n",
       "      <th>70.0</th>\n",
       "      <td>2.027786</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>95.0</th>\n",
       "      <td>1.698108</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>110.0</th>\n",
       "      <td>1.452274</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>125.0</th>\n",
       "      <td>1.198294</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>140.0</th>\n",
       "      <td>0.945816</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>175.0</th>\n",
       "      <td>0.696441</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>225.0</th>\n",
       "      <td>0.665936</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>275.0</th>\n",
       "      <td>0.594450</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>350.0</th>\n",
       "      <td>0.589273</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>425.0</th>\n",
       "      <td>0.495095</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>475.0</th>\n",
       "      <td>0.296315</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "                                    max\n",
       "tag                          50.0-500.0\n",
       "quantity  name      chainage           \n",
       "Discharge tributary 70.0       2.027786\n",
       "                    95.0       1.698108\n",
       "                    110.0      1.452274\n",
       "                    125.0      1.198294\n",
       "                    140.0      0.945816\n",
       "                    175.0      0.696441\n",
       "                    225.0      0.665936\n",
       "                    275.0      0.594450\n",
       "                    350.0      0.589273\n",
       "                    425.0      0.495095\n",
       "                    475.0      0.296315"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Let's start from scratch and use bracket chaining to create the same table\n",
    "df = (\n",
    "    res.read(column_mode=\"all\")\n",
    "    .m1d.query(\"name.str.contains('trib')\")\n",
    "    .m1d.compact()\n",
    "    .m1d.query(\"quantity=='Discharge'\")\n",
    "    .agg([\"max\"])\n",
    "    .T.droplevel(\"group\")\n",
    "    .unstack()\n",
    ")\n",
    "df"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### .m1d.agg_chainage()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>quantity</th>\n",
       "      <th colspan=\"18\" halign=\"left\">Discharge</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>group</th>\n",
       "      <th colspan=\"18\" halign=\"left\">Reach</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th colspan=\"3\" halign=\"left\">basin_left1</th>\n",
       "      <th colspan=\"2\" halign=\"left\">basin_left2</th>\n",
       "      <th colspan=\"2\" halign=\"left\">basin_right</th>\n",
       "      <th>link_basin_left</th>\n",
       "      <th>link_basin_left1_2</th>\n",
       "      <th>link_basin_left2_2</th>\n",
       "      <th>link_basin_right</th>\n",
       "      <th colspan=\"6\" halign=\"left\">river</th>\n",
       "      <th>tributary</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>413.6-483.7</th>\n",
       "      <th>0.0-159.1</th>\n",
       "      <th>159.1-441.7</th>\n",
       "      <th>-10.0-636.4</th>\n",
       "      <th>636.4-720.0</th>\n",
       "      <th>0.0-64.0</th>\n",
       "      <th>0.0-58.6</th>\n",
       "      <th>0.0-69.0</th>\n",
       "      <th>0.0-80.4</th>\n",
       "      <th>53100.0-53950.0</th>\n",
       "      <th>53950.0-54000.0</th>\n",
       "      <th>54000.0-54050.0</th>\n",
       "      <th>54050.0-54150.0</th>\n",
       "      <th>54150.0-54350.0</th>\n",
       "      <th>54350.0-55124.2</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>2.649067</td>\n",
       "      <td>-0.360366</td>\n",
       "      <td>0.025040</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>4.054045</td>\n",
       "      <td>0.008020</td>\n",
       "      <td>4.492208</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>5.737202</td>\n",
       "      <td>72.182541</td>\n",
       "      <td>70.029579</td>\n",
       "      <td>74.721962</td>\n",
       "      <td>74.971359</td>\n",
       "      <td>75.507965</td>\n",
       "      <td>76.910072</td>\n",
       "      <td>0.548908</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>2.888455</td>\n",
       "      <td>0.820759</td>\n",
       "      <td>0.059312</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>2.918259</td>\n",
       "      <td>0.044539</td>\n",
       "      <td>5.646258</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>3.961021</td>\n",
       "      <td>29.658510</td>\n",
       "      <td>27.060493</td>\n",
       "      <td>29.563745</td>\n",
       "      <td>29.512341</td>\n",
       "      <td>28.967613</td>\n",
       "      <td>27.771280</td>\n",
       "      <td>0.489293</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>-0.000545</td>\n",
       "      <td>-2.174818</td>\n",
       "      <td>-0.009495</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.000000</td>\n",
       "      <td>-0.000000</td>\n",
       "      <td>-2.671227</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-0.049854</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>-1.117922</td>\n",
       "      <td>-0.003883</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.036112</td>\n",
       "      <td>75.885757</td>\n",
       "      <td>76.290787</td>\n",
       "      <td>76.206879</td>\n",
       "      <td>75.494148</td>\n",
       "      <td>75.799934</td>\n",
       "      <td>76.689796</td>\n",
       "      <td>0.200483</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>1.183094</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-0.000986</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>4.485793</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.286976</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>6.579083</td>\n",
       "      <td>78.726158</td>\n",
       "      <td>81.531731</td>\n",
       "      <td>85.145889</td>\n",
       "      <td>84.436142</td>\n",
       "      <td>84.838776</td>\n",
       "      <td>86.002777</td>\n",
       "      <td>0.393951</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>5.034901</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>5.796013</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>9.329815</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>8.944139</td>\n",
       "      <td>86.603371</td>\n",
       "      <td>84.478760</td>\n",
       "      <td>93.835732</td>\n",
       "      <td>94.904182</td>\n",
       "      <td>93.394226</td>\n",
       "      <td>93.584366</td>\n",
       "      <td>0.627040</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>7.773971</td>\n",
       "      <td>1.041439</td>\n",
       "      <td>0.202153</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>9.116076</td>\n",
       "      <td>0.366389</td>\n",
       "      <td>13.542682</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>11.380744</td>\n",
       "      <td>102.722717</td>\n",
       "      <td>90.085312</td>\n",
       "      <td>99.735207</td>\n",
       "      <td>100.169701</td>\n",
       "      <td>100.301697</td>\n",
       "      <td>100.134666</td>\n",
       "      <td>2.027786</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity   Discharge                                                  \\\n",
       "group          Reach                                                   \n",
       "name     basin_left1                         basin_left2               \n",
       "tag        0.0-246.6 246.6-413.6 413.6-483.7   0.0-159.1 159.1-441.7   \n",
       "count      73.000000   73.000000   73.000000        73.0        73.0   \n",
       "mean        2.649067   -0.360366    0.025040         0.0         0.0   \n",
       "std         2.888455    0.820759    0.059312         0.0         0.0   \n",
       "min        -0.000545   -2.174818   -0.009495        -0.0        -0.0   \n",
       "25%         0.000000   -1.117922   -0.003883         0.0         0.0   \n",
       "50%         1.183094    0.000000   -0.000986         0.0         0.0   \n",
       "75%         5.034901    0.000000    0.000000         0.0         0.0   \n",
       "max         7.773971    1.041439    0.202153        -0.0        -0.0   \n",
       "\n",
       "quantity                                                             \\\n",
       "group                                                                 \n",
       "name     basin_right             link_basin_left link_basin_left1_2   \n",
       "tag      -10.0-636.4 636.4-720.0        0.0-64.0           0.0-58.6   \n",
       "count      73.000000   73.000000       73.000000               73.0   \n",
       "mean        4.054045    0.008020        4.492208                0.0   \n",
       "std         2.918259    0.044539        5.646258                0.0   \n",
       "min        -0.000000   -0.000000       -2.671227               -0.0   \n",
       "25%         0.000000    0.000000        0.000000                0.0   \n",
       "50%         4.485793    0.000000        1.286976                0.0   \n",
       "75%         5.796013    0.000000        9.329815                0.0   \n",
       "max         9.116076    0.366389       13.542682               -0.0   \n",
       "\n",
       "quantity                                                                      \\\n",
       "group                                                                          \n",
       "name     link_basin_left2_2 link_basin_right           river                   \n",
       "tag                0.0-69.0         0.0-80.4 53100.0-53950.0 53950.0-54000.0   \n",
       "count                  73.0        73.000000       73.000000       73.000000   \n",
       "mean                    0.0         5.737202       72.182541       70.029579   \n",
       "std                     0.0         3.961021       29.658510       27.060493   \n",
       "min                    -0.0         0.000000        0.000000        0.000000   \n",
       "25%                     0.0         0.036112       75.885757       76.290787   \n",
       "50%                     0.0         6.579083       78.726158       81.531731   \n",
       "75%                     0.0         8.944139       86.603371       84.478760   \n",
       "max                    -0.0        11.380744      102.722717       90.085312   \n",
       "\n",
       "quantity                                                                  \\\n",
       "group                                                                      \n",
       "name                                                                       \n",
       "tag      54000.0-54050.0 54050.0-54150.0 54150.0-54350.0 54350.0-55124.2   \n",
       "count          73.000000       73.000000       73.000000       73.000000   \n",
       "mean           74.721962       74.971359       75.507965       76.910072   \n",
       "std            29.563745       29.512341       28.967613       27.771280   \n",
       "min             0.000000        0.000000       -0.049854        0.000000   \n",
       "25%            76.206879       75.494148       75.799934       76.689796   \n",
       "50%            85.145889       84.436142       84.838776       86.002777   \n",
       "75%            93.835732       94.904182       93.394226       93.584366   \n",
       "max            99.735207      100.169701      100.301697      100.134666   \n",
       "\n",
       "quantity             \n",
       "group                \n",
       "name      tributary  \n",
       "tag      50.0-500.0  \n",
       "count     73.000000  \n",
       "mean       0.548908  \n",
       "std        0.489293  \n",
       "min        0.000000  \n",
       "25%        0.200483  \n",
       "50%        0.393951  \n",
       "75%        0.627040  \n",
       "max        2.027786  "
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Now let's try something different. We aggregate the max discharge for each reach, then look at descriptive staistics along the time axis.\n",
    "# Here, 'count' is the number of time steps and 'mean' is the mean of the max discharges of all Q-points along a reach.\n",
    "(\n",
    "    res.read(column_mode=\"all\")\n",
    "    .m1d.query(\"quantity=='Discharge'\")\n",
    "    .m1d.compact()\n",
    "    .m1d.agg_chainage(\"max\")\n",
    "    .describe()\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### .m1d.groupby_chainage()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "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 tr th {\n",
       "        text-align: left;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr>\n",
       "      <th>name</th>\n",
       "      <th colspan=\"6\" halign=\"left\">basin_left1</th>\n",
       "      <th colspan=\"4\" halign=\"left\">basin_left2</th>\n",
       "      <th>...</th>\n",
       "      <th colspan=\"8\" halign=\"left\">river</th>\n",
       "      <th colspan=\"2\" halign=\"left\">tributary</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>chainage</th>\n",
       "      <th>1.002000</th>\n",
       "      <th>228.411000</th>\n",
       "      <th>255.612000</th>\n",
       "      <th>377.384500</th>\n",
       "      <th>447.534000</th>\n",
       "      <th>482.596926</th>\n",
       "      <th>14.597000</th>\n",
       "      <th>126.599500</th>\n",
       "      <th>186.746500</th>\n",
       "      <th>437.951448</th>\n",
       "      <th>...</th>\n",
       "      <th>53975.000000</th>\n",
       "      <th>54025.000000</th>\n",
       "      <th>54075.000000</th>\n",
       "      <th>54125.000000</th>\n",
       "      <th>54175.000000</th>\n",
       "      <th>54324.998000</th>\n",
       "      <th>54375.000000</th>\n",
       "      <th>55112.113830</th>\n",
       "      <th>70.000000</th>\n",
       "      <th>475.000000</th>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>tag</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>0.0-246.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>246.6-413.6</th>\n",
       "      <th>413.6-483.7</th>\n",
       "      <th>413.6-483.7</th>\n",
       "      <th>0.0-159.1</th>\n",
       "      <th>0.0-159.1</th>\n",
       "      <th>159.1-441.7</th>\n",
       "      <th>159.1-441.7</th>\n",
       "      <th>...</th>\n",
       "      <th>53950.0-54000.0</th>\n",
       "      <th>54000.0-54050.0</th>\n",
       "      <th>54050.0-54150.0</th>\n",
       "      <th>54050.0-54150.0</th>\n",
       "      <th>54150.0-54350.0</th>\n",
       "      <th>54150.0-54350.0</th>\n",
       "      <th>54350.0-55124.2</th>\n",
       "      <th>54350.0-55124.2</th>\n",
       "      <th>50.0-500.0</th>\n",
       "      <th>50.0-500.0</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>73.0</td>\n",
       "      <td>...</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>-0.000061</td>\n",
       "      <td>2.306082</td>\n",
       "      <td>-1.521666</td>\n",
       "      <td>-0.416949</td>\n",
       "      <td>-0.049141</td>\n",
       "      <td>-0.001555</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>70.029579</td>\n",
       "      <td>74.721962</td>\n",
       "      <td>74.744926</td>\n",
       "      <td>74.757668</td>\n",
       "      <td>75.087944</td>\n",
       "      <td>75.408974</td>\n",
       "      <td>75.508713</td>\n",
       "      <td>76.717072</td>\n",
       "      <td>-0.211621</td>\n",
       "      <td>0.188293</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>0.000171</td>\n",
       "      <td>3.264935</td>\n",
       "      <td>1.941555</td>\n",
       "      <td>0.753996</td>\n",
       "      <td>0.129560</td>\n",
       "      <td>0.004101</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>27.060493</td>\n",
       "      <td>29.563745</td>\n",
       "      <td>29.496099</td>\n",
       "      <td>29.449484</td>\n",
       "      <td>29.105959</td>\n",
       "      <td>28.908295</td>\n",
       "      <td>28.846025</td>\n",
       "      <td>27.701380</td>\n",
       "      <td>1.107742</td>\n",
       "      <td>0.049437</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>-0.000545</td>\n",
       "      <td>-2.167155</td>\n",
       "      <td>-4.498637</td>\n",
       "      <td>-2.174818</td>\n",
       "      <td>-0.299972</td>\n",
       "      <td>-0.009495</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-0.120941</td>\n",
       "      <td>-0.049854</td>\n",
       "      <td>-0.032419</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-2.884536</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>-3.195148</td>\n",
       "      <td>-1.117922</td>\n",
       "      <td>-0.122781</td>\n",
       "      <td>-0.003883</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>76.290787</td>\n",
       "      <td>76.206879</td>\n",
       "      <td>75.494148</td>\n",
       "      <td>74.803177</td>\n",
       "      <td>75.529945</td>\n",
       "      <td>75.799934</td>\n",
       "      <td>75.881004</td>\n",
       "      <td>76.689796</td>\n",
       "      <td>-0.840475</td>\n",
       "      <td>0.164611</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.138214</td>\n",
       "      <td>-0.351519</td>\n",
       "      <td>-0.255342</td>\n",
       "      <td>-0.031151</td>\n",
       "      <td>-0.000986</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>81.531731</td>\n",
       "      <td>85.145889</td>\n",
       "      <td>84.436142</td>\n",
       "      <td>83.805542</td>\n",
       "      <td>84.253838</td>\n",
       "      <td>84.838776</td>\n",
       "      <td>85.320358</td>\n",
       "      <td>85.784035</td>\n",
       "      <td>-0.223872</td>\n",
       "      <td>0.186661</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>5.034901</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>84.478760</td>\n",
       "      <td>93.835732</td>\n",
       "      <td>94.714996</td>\n",
       "      <td>94.904182</td>\n",
       "      <td>93.394226</td>\n",
       "      <td>93.146133</td>\n",
       "      <td>93.251762</td>\n",
       "      <td>93.584366</td>\n",
       "      <td>0.340223</td>\n",
       "      <td>0.211712</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>0.000301</td>\n",
       "      <td>7.773971</td>\n",
       "      <td>0.993395</td>\n",
       "      <td>0.848020</td>\n",
       "      <td>0.202153</td>\n",
       "      <td>0.006400</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>-0.0</td>\n",
       "      <td>...</td>\n",
       "      <td>90.085312</td>\n",
       "      <td>99.735207</td>\n",
       "      <td>99.707542</td>\n",
       "      <td>100.169701</td>\n",
       "      <td>100.264687</td>\n",
       "      <td>100.144592</td>\n",
       "      <td>100.042068</td>\n",
       "      <td>99.999596</td>\n",
       "      <td>2.027786</td>\n",
       "      <td>0.296315</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "<p>8 rows × 32 columns</p>\n",
       "</div>"
      ],
      "text/plain": [
       "name      basin_left1                                                      \\\n",
       "chainage 1.002000     228.411000   255.612000   377.384500   447.534000     \n",
       "tag         0.0-246.6    0.0-246.6  246.6-413.6  246.6-413.6  413.6-483.7   \n",
       "count       73.000000    73.000000    73.000000    73.000000    73.000000   \n",
       "mean        -0.000061     2.306082    -1.521666    -0.416949    -0.049141   \n",
       "std          0.000171     3.264935     1.941555     0.753996     0.129560   \n",
       "min         -0.000545    -2.167155    -4.498637    -2.174818    -0.299972   \n",
       "25%          0.000000     0.000000    -3.195148    -1.117922    -0.122781   \n",
       "50%          0.000000     1.138214    -0.351519    -0.255342    -0.031151   \n",
       "75%          0.000000     5.034901     0.000000     0.000000     0.000000   \n",
       "max          0.000301     7.773971     0.993395     0.848020     0.202153   \n",
       "\n",
       "name                   basin_left2                                         \\\n",
       "chainage 482.596926   14.597000    126.599500   186.746500   437.951448     \n",
       "tag       413.6-483.7    0.0-159.1    0.0-159.1  159.1-441.7  159.1-441.7   \n",
       "count       73.000000         73.0         73.0         73.0         73.0   \n",
       "mean        -0.001555          0.0          0.0          0.0          0.0   \n",
       "std          0.004101          0.0          0.0          0.0          0.0   \n",
       "min         -0.009495         -0.0         -0.0         -0.0         -0.0   \n",
       "25%         -0.003883          0.0          0.0          0.0          0.0   \n",
       "50%         -0.000986          0.0          0.0          0.0          0.0   \n",
       "75%          0.000000          0.0          0.0          0.0          0.0   \n",
       "max          0.006400         -0.0         -0.0         -0.0         -0.0   \n",
       "\n",
       "name      ...           river                                                  \\\n",
       "chainage  ...    53975.000000    54025.000000    54075.000000    54125.000000   \n",
       "tag       ... 53950.0-54000.0 54000.0-54050.0 54050.0-54150.0 54050.0-54150.0   \n",
       "count     ...       73.000000       73.000000       73.000000       73.000000   \n",
       "mean      ...       70.029579       74.721962       74.744926       74.757668   \n",
       "std       ...       27.060493       29.563745       29.496099       29.449484   \n",
       "min       ...        0.000000        0.000000        0.000000        0.000000   \n",
       "25%       ...       76.290787       76.206879       75.494148       74.803177   \n",
       "50%       ...       81.531731       85.145889       84.436142       83.805542   \n",
       "75%       ...       84.478760       93.835732       94.714996       94.904182   \n",
       "max       ...       90.085312       99.735207       99.707542      100.169701   \n",
       "\n",
       "name                                                                      \\\n",
       "chainage    54175.000000    54324.998000    54375.000000    55112.113830   \n",
       "tag      54150.0-54350.0 54150.0-54350.0 54350.0-55124.2 54350.0-55124.2   \n",
       "count          73.000000       73.000000       73.000000       73.000000   \n",
       "mean           75.087944       75.408974       75.508713       76.717072   \n",
       "std            29.105959       28.908295       28.846025       27.701380   \n",
       "min            -0.120941       -0.049854       -0.032419        0.000000   \n",
       "25%            75.529945       75.799934       75.881004       76.689796   \n",
       "50%            84.253838       84.838776       85.320358       85.784035   \n",
       "75%            93.394226       93.146133       93.251762       93.584366   \n",
       "max           100.264687      100.144592      100.042068       99.999596   \n",
       "\n",
       "name        tributary               \n",
       "chainage 70.000000    475.000000    \n",
       "tag        50.0-500.0   50.0-500.0  \n",
       "count       73.000000    73.000000  \n",
       "mean        -0.211621     0.188293  \n",
       "std          1.107742     0.049437  \n",
       "min         -2.884536     0.000000  \n",
       "25%         -0.840475     0.164611  \n",
       "50%         -0.223872     0.186661  \n",
       "75%          0.340223     0.211712  \n",
       "max          2.027786     0.296315  \n",
       "\n",
       "[8 rows x 32 columns]"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# Similarly, let's look at the mean of the first and last Q points.\n",
    "(\n",
    "    res.read(column_mode=\"all\")\n",
    "    .m1d.query(\"quantity=='Discharge'\")\n",
    "    .m1d.compact()\n",
    "    .m1d.groupby_chainage()\n",
    "    .nth([0, -1])  # First we groupby, then select the first and last chaianges.\n",
    "    .describe()\n",
    "    .droplevel([\"quantity\", \"group\"], axis=1)\n",
    ")"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "#### .m1d.groupby()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "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>quantity</th>\n",
       "      <th>ControlStrategyId</th>\n",
       "      <th>Discharge</th>\n",
       "      <th>Discharge:Sensor:SensorGauge1</th>\n",
       "      <th>DischargeInStructure</th>\n",
       "      <th>FlowAreaInStructure</th>\n",
       "      <th>FlowVelocity</th>\n",
       "      <th>FlowVelocityInStructure</th>\n",
       "      <th>Gate level:Sensor:SensorGateLevel</th>\n",
       "      <th>GateLevel</th>\n",
       "      <th>ManningResistanceNumber</th>\n",
       "      <th>Variable:TwoTimeSensorGateLevel</th>\n",
       "      <th>Water level:Sensor:s.h.river53745.34</th>\n",
       "      <th>WaterLevel</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>count</th>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "      <td>73.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>mean</th>\n",
       "      <td>5.917808</td>\n",
       "      <td>79.343079</td>\n",
       "      <td>61.396931</td>\n",
       "      <td>77.747810</td>\n",
       "      <td>47.783684</td>\n",
       "      <td>3.353872</td>\n",
       "      <td>11.329100</td>\n",
       "      <td>53.769997</td>\n",
       "      <td>53.769997</td>\n",
       "      <td>32.665512</td>\n",
       "      <td>107.539993</td>\n",
       "      <td>56.858250</td>\n",
       "      <td>59.200008</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>std</th>\n",
       "      <td>0.702247</td>\n",
       "      <td>28.940966</td>\n",
       "      <td>26.605530</td>\n",
       "      <td>29.185284</td>\n",
       "      <td>18.166550</td>\n",
       "      <td>0.818552</td>\n",
       "      <td>5.346817</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>2.485120</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>1.846656</td>\n",
       "      <td>0.000000</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>min</th>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>0.000000</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>25.000000</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>51.905998</td>\n",
       "      <td>59.200001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25%</th>\n",
       "      <td>6.000000</td>\n",
       "      <td>80.783112</td>\n",
       "      <td>69.281067</td>\n",
       "      <td>76.470184</td>\n",
       "      <td>53.639187</td>\n",
       "      <td>3.652213</td>\n",
       "      <td>11.170768</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>33.597073</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>57.567822</td>\n",
       "      <td>59.200001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>50%</th>\n",
       "      <td>6.000000</td>\n",
       "      <td>86.655334</td>\n",
       "      <td>74.483246</td>\n",
       "      <td>86.155029</td>\n",
       "      <td>56.015728</td>\n",
       "      <td>3.719032</td>\n",
       "      <td>11.416133</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>33.614704</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>57.708385</td>\n",
       "      <td>59.200001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>75%</th>\n",
       "      <td>6.000000</td>\n",
       "      <td>99.554710</td>\n",
       "      <td>75.963501</td>\n",
       "      <td>99.477150</td>\n",
       "      <td>57.551910</td>\n",
       "      <td>3.736953</td>\n",
       "      <td>11.514374</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>33.618801</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>57.743835</td>\n",
       "      <td>59.200001</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>max</th>\n",
       "      <td>6.000000</td>\n",
       "      <td>102.722717</td>\n",
       "      <td>77.521797</td>\n",
       "      <td>100.247269</td>\n",
       "      <td>59.036644</td>\n",
       "      <td>3.754651</td>\n",
       "      <td>19.363970</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>53.770000</td>\n",
       "      <td>33.622650</td>\n",
       "      <td>107.540001</td>\n",
       "      <td>57.777584</td>\n",
       "      <td>59.200001</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "quantity  ControlStrategyId   Discharge  Discharge:Sensor:SensorGauge1  \\\n",
       "count             73.000000   73.000000                      73.000000   \n",
       "mean               5.917808   79.343079                      61.396931   \n",
       "std                0.702247   28.940966                      26.605530   \n",
       "min                0.000000    0.000000                       0.000000   \n",
       "25%                6.000000   80.783112                      69.281067   \n",
       "50%                6.000000   86.655334                      74.483246   \n",
       "75%                6.000000   99.554710                      75.963501   \n",
       "max                6.000000  102.722717                      77.521797   \n",
       "\n",
       "quantity  DischargeInStructure  FlowAreaInStructure  FlowVelocity  \\\n",
       "count                73.000000            73.000000     73.000000   \n",
       "mean                 77.747810            47.783684      3.353872   \n",
       "std                  29.185284            18.166550      0.818552   \n",
       "min                   0.000000             0.000000      0.000000   \n",
       "25%                  76.470184            53.639187      3.652213   \n",
       "50%                  86.155029            56.015728      3.719032   \n",
       "75%                  99.477150            57.551910      3.736953   \n",
       "max                 100.247269            59.036644      3.754651   \n",
       "\n",
       "quantity  FlowVelocityInStructure  Gate level:Sensor:SensorGateLevel  \\\n",
       "count                   73.000000                          73.000000   \n",
       "mean                    11.329100                          53.769997   \n",
       "std                      5.346817                           0.000000   \n",
       "min                      0.000000                          53.770000   \n",
       "25%                     11.170768                          53.770000   \n",
       "50%                     11.416133                          53.770000   \n",
       "75%                     11.514374                          53.770000   \n",
       "max                     19.363970                          53.770000   \n",
       "\n",
       "quantity  GateLevel  ManningResistanceNumber  Variable:TwoTimeSensorGateLevel  \\\n",
       "count     73.000000                73.000000                        73.000000   \n",
       "mean      53.769997                32.665512                       107.539993   \n",
       "std        0.000000                 2.485120                         0.000000   \n",
       "min       53.770000                25.000000                       107.540001   \n",
       "25%       53.770000                33.597073                       107.540001   \n",
       "50%       53.770000                33.614704                       107.540001   \n",
       "75%       53.770000                33.618801                       107.540001   \n",
       "max       53.770000                33.622650                       107.540001   \n",
       "\n",
       "quantity  Water level:Sensor:s.h.river53745.34  WaterLevel  \n",
       "count                                73.000000   73.000000  \n",
       "mean                                 56.858250   59.200008  \n",
       "std                                   1.846656    0.000000  \n",
       "min                                  51.905998   59.200001  \n",
       "25%                                  57.567822   59.200001  \n",
       "50%                                  57.708385   59.200001  \n",
       "75%                                  57.743835   59.200001  \n",
       "max                                  57.777584   59.200001  "
      ]
     },
     "execution_count": 15,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# We can similarly use the .m1d accessor to group by dimensions other than chainage.\n",
    "# Below we describe how the global maximum of each quantity in the system varies with time.\n",
    "(res.read(column_mode=\"all\").m1d.groupby(\"quantity\").agg(\"max\").describe())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": ".venv",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.12.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}