Skip to content

Instantly share code, notes, and snippets.

@kuchaale

kuchaale/GW_activity_t_tide_analysis.ipynb

Created August 12, 2021 17:05
Show Gist options
  • Save kuchaale/ba28fb8348bbedf044a0c23ad5b9be5d to your computer and use it in GitHub Desktop.
Save kuchaale/ba28fb8348bbedf044a0c23ad5b9be5d to your computer and use it in GitHub Desktop.
Download ZIP
Display the source blob
Display the rendered blob
Raw
GW_activity_t_tide_analysis.ipynb
{
"cells": [
{
"cell_type": "code",
"execution_count": 22,
"metadata": {},
"outputs": [],
"source": [
"import xarray as xr\n",
"import matplotlib.pyplot as plt\n",
"#import matplotlib as mpl\n",
"#from matplotlib.ticker import MultipleLocator, FormatStrFormatter\n",
"from pathlib import Path\n",
"import numpy as np\n",
"import pandas as pd"
]
},
{
"cell_type": "code",
"execution_count": 21,
"metadata": {},
"outputs": [],
"source": [
"root_path = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/'\n",
"station = 'Davis'\n",
"alt = [\n",
" 70,\n",
" 72,\n",
" 74,\n",
" 76,\n",
" 78,\n",
" 80,\n",
" 82,\n",
" 84,\n",
" 86,\n",
" 88,\n",
" 90,\n",
" 92,\n",
" 94,\n",
" 96,\n",
" 98,\n",
" 100,\n",
" 102,\n",
" 104,\n",
" 106,\n",
" 108,\n",
" 110,\n",
" 112,\n",
" 114,\n",
" 116,\n",
" 118,\n",
" 120,\n",
"]\n",
"s_date = '2005-01-01T00'"
]
},
{
"cell_type": "code",
"execution_count": 23,
"metadata": {},
"outputs": [],
"source": [
"def preprocess(ds):\n",
" ds = ds.rename({\"phony_dim_8\": \"time\", \"phony_dim_9\": \"alt\"})\n",
" ds['alt'] = alt\n",
" ds['alt'].attrs['units'] = 'km'\n",
" ds['alt'].attrs['long_name'] = 'altitude'\n",
" ds['time'] = pd.date_range(s_date, periods=ds.time.shape[0], freq = '1D')\n",
" return ds\n"
]
},
{
"cell_type": "code",
"execution_count": 24,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
"<defs>\n",
"<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
"<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
"<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
"<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
"</symbol>\n",
"<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
"<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
"<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
"<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
"<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
"</symbol>\n",
"</defs>\n",
"</svg>\n",
"<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
" *\n",
" */\n",
"\n",
":root {\n",
" --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
" --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
" --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
" --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
" --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
" --xr-background-color: var(--jp-layout-color0, white);\n",
" --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
" --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
"}\n",
"\n",
"html[theme=dark],\n",
"body.vscode-dark {\n",
" --xr-font-color0: rgba(255, 255, 255, 1);\n",
" --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
" --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
" --xr-border-color: #1F1F1F;\n",
" --xr-disabled-color: #515151;\n",
" --xr-background-color: #111111;\n",
" --xr-background-color-row-even: #111111;\n",
" --xr-background-color-row-odd: #313131;\n",
"}\n",
"\n",
".xr-wrap {\n",
" display: block;\n",
" min-width: 300px;\n",
" max-width: 700px;\n",
"}\n",
"\n",
".xr-text-repr-fallback {\n",
" /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
" display: none;\n",
"}\n",
"\n",
".xr-header {\n",
" padding-top: 6px;\n",
" padding-bottom: 6px;\n",
" margin-bottom: 4px;\n",
" border-bottom: solid 1px var(--xr-border-color);\n",
"}\n",
"\n",
".xr-header > div,\n",
".xr-header > ul {\n",
" display: inline;\n",
" margin-top: 0;\n",
" margin-bottom: 0;\n",
"}\n",
"\n",
".xr-obj-type,\n",
".xr-array-name {\n",
" margin-left: 2px;\n",
" margin-right: 10px;\n",
"}\n",
"\n",
".xr-obj-type {\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-sections {\n",
" padding-left: 0 !important;\n",
" display: grid;\n",
" grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
"}\n",
"\n",
".xr-section-item {\n",
" display: contents;\n",
"}\n",
"\n",
".xr-section-item input {\n",
" display: none;\n",
"}\n",
"\n",
".xr-section-item input + label {\n",
" color: var(--xr-disabled-color);\n",
"}\n",
"\n",
".xr-section-item input:enabled + label {\n",
" cursor: pointer;\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-section-item input:enabled + label:hover {\n",
" color: var(--xr-font-color0);\n",
"}\n",
"\n",
".xr-section-summary {\n",
" grid-column: 1;\n",
" color: var(--xr-font-color2);\n",
" font-weight: 500;\n",
"}\n",
"\n",
".xr-section-summary > span {\n",
" display: inline-block;\n",
" padding-left: 0.5em;\n",
"}\n",
"\n",
".xr-section-summary-in:disabled + label {\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-section-summary-in + label:before {\n",
" display: inline-block;\n",
" content: '►';\n",
" font-size: 11px;\n",
" width: 15px;\n",
" text-align: center;\n",
"}\n",
"\n",
".xr-section-summary-in:disabled + label:before {\n",
" color: var(--xr-disabled-color);\n",
"}\n",
"\n",
".xr-section-summary-in:checked + label:before {\n",
" content: '▼';\n",
"}\n",
"\n",
".xr-section-summary-in:checked + label > span {\n",
" display: none;\n",
"}\n",
"\n",
".xr-section-summary,\n",
".xr-section-inline-details {\n",
" padding-top: 4px;\n",
" padding-bottom: 4px;\n",
"}\n",
"\n",
".xr-section-inline-details {\n",
" grid-column: 2 / -1;\n",
"}\n",
"\n",
".xr-section-details {\n",
" display: none;\n",
" grid-column: 1 / -1;\n",
" margin-bottom: 5px;\n",
"}\n",
"\n",
".xr-section-summary-in:checked ~ .xr-section-details {\n",
" display: contents;\n",
"}\n",
"\n",
".xr-array-wrap {\n",
" grid-column: 1 / -1;\n",
" display: grid;\n",
" grid-template-columns: 20px auto;\n",
"}\n",
"\n",
".xr-array-wrap > label {\n",
" grid-column: 1;\n",
" vertical-align: top;\n",
"}\n",
"\n",
".xr-preview {\n",
" color: var(--xr-font-color3);\n",
"}\n",
"\n",
".xr-array-preview,\n",
".xr-array-data {\n",
" padding: 0 5px !important;\n",
" grid-column: 2;\n",
"}\n",
"\n",
".xr-array-data,\n",
".xr-array-in:checked ~ .xr-array-preview {\n",
" display: none;\n",
"}\n",
"\n",
".xr-array-in:checked ~ .xr-array-data,\n",
".xr-array-preview {\n",
" display: inline-block;\n",
"}\n",
"\n",
".xr-dim-list {\n",
" display: inline-block !important;\n",
" list-style: none;\n",
" padding: 0 !important;\n",
" margin: 0;\n",
"}\n",
"\n",
".xr-dim-list li {\n",
" display: inline-block;\n",
" padding: 0;\n",
" margin: 0;\n",
"}\n",
"\n",
".xr-dim-list:before {\n",
" content: '(';\n",
"}\n",
"\n",
".xr-dim-list:after {\n",
" content: ')';\n",
"}\n",
"\n",
".xr-dim-list li:not(:last-child):after {\n",
" content: ',';\n",
" padding-right: 5px;\n",
"}\n",
"\n",
".xr-has-index {\n",
" font-weight: bold;\n",
"}\n",
"\n",
".xr-var-list,\n",
".xr-var-item {\n",
" display: contents;\n",
"}\n",
"\n",
".xr-var-item > div,\n",
".xr-var-item label,\n",
".xr-var-item > .xr-var-name span {\n",
" background-color: var(--xr-background-color-row-even);\n",
" margin-bottom: 0;\n",
"}\n",
"\n",
".xr-var-item > .xr-var-name:hover span {\n",
" padding-right: 5px;\n",
"}\n",
"\n",
".xr-var-list > li:nth-child(odd) > div,\n",
".xr-var-list > li:nth-child(odd) > label,\n",
".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
" background-color: var(--xr-background-color-row-odd);\n",
"}\n",
"\n",
".xr-var-name {\n",
" grid-column: 1;\n",
"}\n",
"\n",
".xr-var-dims {\n",
" grid-column: 2;\n",
"}\n",
"\n",
".xr-var-dtype {\n",
" grid-column: 3;\n",
" text-align: right;\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-var-preview {\n",
" grid-column: 4;\n",
"}\n",
"\n",
".xr-var-name,\n",
".xr-var-dims,\n",
".xr-var-dtype,\n",
".xr-preview,\n",
".xr-attrs dt {\n",
" white-space: nowrap;\n",
" overflow: hidden;\n",
" text-overflow: ellipsis;\n",
" padding-right: 10px;\n",
"}\n",
"\n",
".xr-var-name:hover,\n",
".xr-var-dims:hover,\n",
".xr-var-dtype:hover,\n",
".xr-attrs dt:hover {\n",
" overflow: visible;\n",
" width: auto;\n",
" z-index: 1;\n",
"}\n",
"\n",
".xr-var-attrs,\n",
".xr-var-data {\n",
" display: none;\n",
" background-color: var(--xr-background-color) !important;\n",
" padding-bottom: 5px !important;\n",
"}\n",
"\n",
".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
".xr-var-data-in:checked ~ .xr-var-data {\n",
" display: block;\n",
"}\n",
"\n",
".xr-var-data > table {\n",
" float: right;\n",
"}\n",
"\n",
".xr-var-name span,\n",
".xr-var-data,\n",
".xr-attrs {\n",
" padding-left: 25px !important;\n",
"}\n",
"\n",
".xr-attrs,\n",
".xr-var-attrs,\n",
".xr-var-data {\n",
" grid-column: 1 / -1;\n",
"}\n",
"\n",
"dl.xr-attrs {\n",
" padding: 0;\n",
" margin: 0;\n",
" display: grid;\n",
" grid-template-columns: 125px auto;\n",
"}\n",
"\n",
".xr-attrs dt,\n",
".xr-attrs dd {\n",
" padding: 0;\n",
" margin: 0;\n",
" float: left;\n",
" padding-right: 10px;\n",
" width: auto;\n",
"}\n",
"\n",
".xr-attrs dt {\n",
" font-weight: normal;\n",
" grid-column: 1;\n",
"}\n",
"\n",
".xr-attrs dt:hover span {\n",
" display: inline-block;\n",
" background: var(--xr-background-color);\n",
" padding-right: 10px;\n",
"}\n",
"\n",
".xr-attrs dd {\n",
" grid-column: 2;\n",
" white-space: pre-wrap;\n",
" word-break: break-all;\n",
"}\n",
"\n",
".xr-icon-database,\n",
".xr-icon-file-text2 {\n",
" display: inline-block;\n",
" vertical-align: middle;\n",
" width: 1em;\n",
" height: 1.5em !important;\n",
" stroke-width: 0;\n",
" stroke: currentColor;\n",
" fill: currentColor;\n",
"}\n",
"</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
"Dimensions: (alt: 26, time: 5844)\n",
"Coordinates:\n",
" * alt (alt) int64 70 72 74 76 78 80 82 84 ... 108 110 112 114 116 118 120\n",
" * time (time) datetime64[ns] 2005-01-01 2005-01-02 ... 2020-12-31\n",
"Data variables:\n",
" GW (time, alt) float64 ...\n",
" GW_err (time, alt) float64 ...\n",
" PW (time, alt) float64 ...\n",
" PW_err (time, alt) float64 ...</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-d4d01dc5-8c82-4851-b2c1-7689c3a72a89' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-d4d01dc5-8c82-4851-b2c1-7689c3a72a89' class='xr-section-summary' title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>alt</span>: 26</li><li><span class='xr-has-index'>time</span>: 5844</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-d1568142-1436-47da-b3db-9e7d006b231e' class='xr-section-summary-in' type='checkbox' checked><label for='section-d1568142-1436-47da-b3db-9e7d006b231e' class='xr-section-summary' >Coordinates: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>alt</span></div><div class='xr-var-dims'>(alt)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>70 72 74 76 78 ... 114 116 118 120</div><input id='attrs-f7c70ce1-6f9d-4f0d-82c6-b12711d8042f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f7c70ce1-6f9d-4f0d-82c6-b12711d8042f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d6d8f35f-7527-4d67-b9fa-1cad8c90c37f' class='xr-var-data-in' type='checkbox'><label for='data-d6d8f35f-7527-4d67-b9fa-1cad8c90c37f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>km</dd><dt><span>long_name :</span></dt><dd>altitude</dd></dl></div><div class='xr-var-data'><pre>array([ 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96,\n",
" 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2005-01-01 ... 2020-12-31</div><input id='attrs-48cc6088-bff4-43e1-b180-6264cd98586c' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-48cc6088-bff4-43e1-b180-6264cd98586c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fdbcae2d-3ea4-4464-bfab-792fe93ea950' class='xr-var-data-in' type='checkbox'><label for='data-fdbcae2d-3ea4-4464-bfab-792fe93ea950' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2005-01-01T00:00:00.000000000&#x27;, &#x27;2005-01-02T00:00:00.000000000&#x27;,\n",
" &#x27;2005-01-03T00:00:00.000000000&#x27;, ..., &#x27;2020-12-29T00:00:00.000000000&#x27;,\n",
" &#x27;2020-12-30T00:00:00.000000000&#x27;, &#x27;2020-12-31T00:00:00.000000000&#x27;],\n",
" dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-21895a93-383b-4d0e-9d1d-411d8f183cf8' class='xr-section-summary-in' type='checkbox' checked><label for='section-21895a93-383b-4d0e-9d1d-411d8f183cf8' class='xr-section-summary' >Data variables: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>GW</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-cd3aa82b-62c6-47ba-81af-5a11c50e1c8f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-cd3aa82b-62c6-47ba-81af-5a11c50e1c8f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c7613533-e9db-490d-9770-af68f59f7393' class='xr-var-data-in' type='checkbox'><label for='data-c7613533-e9db-490d-9770-af68f59f7393' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>GW energy :</span></dt><dd>GW energy / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>GW_err</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-df35135d-aaec-4bdb-a265-197713511afe' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-df35135d-aaec-4bdb-a265-197713511afe' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0512ca8b-aaf0-43ce-8889-79e175bbfeb1' class='xr-var-data-in' type='checkbox'><label for='data-0512ca8b-aaf0-43ce-8889-79e175bbfeb1' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>PW energy error :</span></dt><dd>GW energy error / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>PW</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-0f3d108e-9c48-45a7-ab5f-d2051d02a3cd' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-0f3d108e-9c48-45a7-ab5f-d2051d02a3cd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2a2e3352-05f7-4f2f-a86e-3d4824bbc671' class='xr-var-data-in' type='checkbox'><label for='data-2a2e3352-05f7-4f2f-a86e-3d4824bbc671' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>PW energy :</span></dt><dd>PW energyl / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>PW_err</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-5095204f-90a1-45a0-92d5-fc5d7eb06ecf' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-5095204f-90a1-45a0-92d5-fc5d7eb06ecf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8803e5a6-66c2-478d-8bfa-228ac4b5dca2' class='xr-var-data-in' type='checkbox'><label for='data-8803e5a6-66c2-478d-8bfa-228ac4b5dca2' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>PW energy error :</span></dt><dd>PW energy error / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-a3884559-3fc5-48e4-ab5b-1d85d42495f4' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-a3884559-3fc5-48e4-ab5b-1d85d42495f4' class='xr-section-summary' title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
],
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (alt: 26, time: 5844)\n",
"Coordinates:\n",
" * alt (alt) int64 70 72 74 76 78 80 82 84 ... 108 110 112 114 116 118 120\n",
" * time (time) datetime64[ns] 2005-01-01 2005-01-02 ... 2020-12-31\n",
"Data variables:\n",
" GW (time, alt) float64 ...\n",
" GW_err (time, alt) float64 ...\n",
" PW (time, alt) float64 ...\n",
" PW_err (time, alt) float64 ..."
]
},
"execution_count": 24,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds = xr.open_dataset(f'{root_path}{station}_solar_cycle.h5', group = 'waves')\n",
"ds = preprocess(ds)\n",
"ds"
]
},
{
"cell_type": "code",
"execution_count": 25,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
"<defs>\n",
"<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
"<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
"<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
"<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
"</symbol>\n",
"<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
"<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
"<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
"<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
"<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
"</symbol>\n",
"</defs>\n",
"</svg>\n",
"<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
" *\n",
" */\n",
"\n",
":root {\n",
" --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
" --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
" --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
" --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
" --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
" --xr-background-color: var(--jp-layout-color0, white);\n",
" --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
" --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
"}\n",
"\n",
"html[theme=dark],\n",
"body.vscode-dark {\n",
" --xr-font-color0: rgba(255, 255, 255, 1);\n",
" --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
" --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
" --xr-border-color: #1F1F1F;\n",
" --xr-disabled-color: #515151;\n",
" --xr-background-color: #111111;\n",
" --xr-background-color-row-even: #111111;\n",
" --xr-background-color-row-odd: #313131;\n",
"}\n",
"\n",
".xr-wrap {\n",
" display: block;\n",
" min-width: 300px;\n",
" max-width: 700px;\n",
"}\n",
"\n",
".xr-text-repr-fallback {\n",
" /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
" display: none;\n",
"}\n",
"\n",
".xr-header {\n",
" padding-top: 6px;\n",
" padding-bottom: 6px;\n",
" margin-bottom: 4px;\n",
" border-bottom: solid 1px var(--xr-border-color);\n",
"}\n",
"\n",
".xr-header > div,\n",
".xr-header > ul {\n",
" display: inline;\n",
" margin-top: 0;\n",
" margin-bottom: 0;\n",
"}\n",
"\n",
".xr-obj-type,\n",
".xr-array-name {\n",
" margin-left: 2px;\n",
" margin-right: 10px;\n",
"}\n",
"\n",
".xr-obj-type {\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-sections {\n",
" padding-left: 0 !important;\n",
" display: grid;\n",
" grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
"}\n",
"\n",
".xr-section-item {\n",
" display: contents;\n",
"}\n",
"\n",
".xr-section-item input {\n",
" display: none;\n",
"}\n",
"\n",
".xr-section-item input + label {\n",
" color: var(--xr-disabled-color);\n",
"}\n",
"\n",
".xr-section-item input:enabled + label {\n",
" cursor: pointer;\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-section-item input:enabled + label:hover {\n",
" color: var(--xr-font-color0);\n",
"}\n",
"\n",
".xr-section-summary {\n",
" grid-column: 1;\n",
" color: var(--xr-font-color2);\n",
" font-weight: 500;\n",
"}\n",
"\n",
".xr-section-summary > span {\n",
" display: inline-block;\n",
" padding-left: 0.5em;\n",
"}\n",
"\n",
".xr-section-summary-in:disabled + label {\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-section-summary-in + label:before {\n",
" display: inline-block;\n",
" content: '►';\n",
" font-size: 11px;\n",
" width: 15px;\n",
" text-align: center;\n",
"}\n",
"\n",
".xr-section-summary-in:disabled + label:before {\n",
" color: var(--xr-disabled-color);\n",
"}\n",
"\n",
".xr-section-summary-in:checked + label:before {\n",
" content: '▼';\n",
"}\n",
"\n",
".xr-section-summary-in:checked + label > span {\n",
" display: none;\n",
"}\n",
"\n",
".xr-section-summary,\n",
".xr-section-inline-details {\n",
" padding-top: 4px;\n",
" padding-bottom: 4px;\n",
"}\n",
"\n",
".xr-section-inline-details {\n",
" grid-column: 2 / -1;\n",
"}\n",
"\n",
".xr-section-details {\n",
" display: none;\n",
" grid-column: 1 / -1;\n",
" margin-bottom: 5px;\n",
"}\n",
"\n",
".xr-section-summary-in:checked ~ .xr-section-details {\n",
" display: contents;\n",
"}\n",
"\n",
".xr-array-wrap {\n",
" grid-column: 1 / -1;\n",
" display: grid;\n",
" grid-template-columns: 20px auto;\n",
"}\n",
"\n",
".xr-array-wrap > label {\n",
" grid-column: 1;\n",
" vertical-align: top;\n",
"}\n",
"\n",
".xr-preview {\n",
" color: var(--xr-font-color3);\n",
"}\n",
"\n",
".xr-array-preview,\n",
".xr-array-data {\n",
" padding: 0 5px !important;\n",
" grid-column: 2;\n",
"}\n",
"\n",
".xr-array-data,\n",
".xr-array-in:checked ~ .xr-array-preview {\n",
" display: none;\n",
"}\n",
"\n",
".xr-array-in:checked ~ .xr-array-data,\n",
".xr-array-preview {\n",
" display: inline-block;\n",
"}\n",
"\n",
".xr-dim-list {\n",
" display: inline-block !important;\n",
" list-style: none;\n",
" padding: 0 !important;\n",
" margin: 0;\n",
"}\n",
"\n",
".xr-dim-list li {\n",
" display: inline-block;\n",
" padding: 0;\n",
" margin: 0;\n",
"}\n",
"\n",
".xr-dim-list:before {\n",
" content: '(';\n",
"}\n",
"\n",
".xr-dim-list:after {\n",
" content: ')';\n",
"}\n",
"\n",
".xr-dim-list li:not(:last-child):after {\n",
" content: ',';\n",
" padding-right: 5px;\n",
"}\n",
"\n",
".xr-has-index {\n",
" font-weight: bold;\n",
"}\n",
"\n",
".xr-var-list,\n",
".xr-var-item {\n",
" display: contents;\n",
"}\n",
"\n",
".xr-var-item > div,\n",
".xr-var-item label,\n",
".xr-var-item > .xr-var-name span {\n",
" background-color: var(--xr-background-color-row-even);\n",
" margin-bottom: 0;\n",
"}\n",
"\n",
".xr-var-item > .xr-var-name:hover span {\n",
" padding-right: 5px;\n",
"}\n",
"\n",
".xr-var-list > li:nth-child(odd) > div,\n",
".xr-var-list > li:nth-child(odd) > label,\n",
".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
" background-color: var(--xr-background-color-row-odd);\n",
"}\n",
"\n",
".xr-var-name {\n",
" grid-column: 1;\n",
"}\n",
"\n",
".xr-var-dims {\n",
" grid-column: 2;\n",
"}\n",
"\n",
".xr-var-dtype {\n",
" grid-column: 3;\n",
" text-align: right;\n",
" color: var(--xr-font-color2);\n",
"}\n",
"\n",
".xr-var-preview {\n",
" grid-column: 4;\n",
"}\n",
"\n",
".xr-var-name,\n",
".xr-var-dims,\n",
".xr-var-dtype,\n",
".xr-preview,\n",
".xr-attrs dt {\n",
" white-space: nowrap;\n",
" overflow: hidden;\n",
" text-overflow: ellipsis;\n",
" padding-right: 10px;\n",
"}\n",
"\n",
".xr-var-name:hover,\n",
".xr-var-dims:hover,\n",
".xr-var-dtype:hover,\n",
".xr-attrs dt:hover {\n",
" overflow: visible;\n",
" width: auto;\n",
" z-index: 1;\n",
"}\n",
"\n",
".xr-var-attrs,\n",
".xr-var-data {\n",
" display: none;\n",
" background-color: var(--xr-background-color) !important;\n",
" padding-bottom: 5px !important;\n",
"}\n",
"\n",
".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
".xr-var-data-in:checked ~ .xr-var-data {\n",
" display: block;\n",
"}\n",
"\n",
".xr-var-data > table {\n",
" float: right;\n",
"}\n",
"\n",
".xr-var-name span,\n",
".xr-var-data,\n",
".xr-attrs {\n",
" padding-left: 25px !important;\n",
"}\n",
"\n",
".xr-attrs,\n",
".xr-var-attrs,\n",
".xr-var-data {\n",
" grid-column: 1 / -1;\n",
"}\n",
"\n",
"dl.xr-attrs {\n",
" padding: 0;\n",
" margin: 0;\n",
" display: grid;\n",
" grid-template-columns: 125px auto;\n",
"}\n",
"\n",
".xr-attrs dt,\n",
".xr-attrs dd {\n",
" padding: 0;\n",
" margin: 0;\n",
" float: left;\n",
" padding-right: 10px;\n",
" width: auto;\n",
"}\n",
"\n",
".xr-attrs dt {\n",
" font-weight: normal;\n",
" grid-column: 1;\n",
"}\n",
"\n",
".xr-attrs dt:hover span {\n",
" display: inline-block;\n",
" background: var(--xr-background-color);\n",
" padding-right: 10px;\n",
"}\n",
"\n",
".xr-attrs dd {\n",
" grid-column: 2;\n",
" white-space: pre-wrap;\n",
" word-break: break-all;\n",
"}\n",
"\n",
".xr-icon-database,\n",
".xr-icon-file-text2 {\n",
" display: inline-block;\n",
" vertical-align: middle;\n",
" width: 1em;\n",
" height: 1.5em !important;\n",
" stroke-width: 0;\n",
" stroke: currentColor;\n",
" fill: currentColor;\n",
"}\n",
"</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt;\n",
"Dimensions: (alt: 26, time: 5844)\n",
"Coordinates:\n",
" * alt (alt) int64 70 72 74 76 78 80 82 84 ... 108 110 112 114 116 118 120\n",
" * time (time) datetime64[ns] 2005-01-01 2005-01-02 ... 2020-12-31\n",
"Data variables:\n",
" GW (time, alt) float64 ...\n",
" GW_err (time, alt) float64 ...\n",
" PW (time, alt) float64 ...\n",
" PW_err (time, alt) float64 ...</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-45a12275-847c-4548-a531-28c66fbba34a' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-45a12275-847c-4548-a531-28c66fbba34a' class='xr-section-summary' title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>alt</span>: 26</li><li><span class='xr-has-index'>time</span>: 5844</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-2f27386f-099c-4be9-9e26-b66aa3ef3c57' class='xr-section-summary-in' type='checkbox' checked><label for='section-2f27386f-099c-4be9-9e26-b66aa3ef3c57' class='xr-section-summary' >Coordinates: <span>(2)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>alt</span></div><div class='xr-var-dims'>(alt)</div><div class='xr-var-dtype'>int64</div><div class='xr-var-preview xr-preview'>70 72 74 76 78 ... 114 116 118 120</div><input id='attrs-9391fcdd-459d-4560-829b-1f9b335f9621' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9391fcdd-459d-4560-829b-1f9b335f9621' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7bddee88-3fb0-4f9d-9819-2c8a496aeb57' class='xr-var-data-in' type='checkbox'><label for='data-7bddee88-3fb0-4f9d-9819-2c8a496aeb57' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>units :</span></dt><dd>km</dd><dt><span>long_name :</span></dt><dd>altitude</dd></dl></div><div class='xr-var-data'><pre>array([ 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96,\n",
" 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>2005-01-01 ... 2020-12-31</div><input id='attrs-1e483add-0b26-4ebc-858f-a80bb805e6fc' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-1e483add-0b26-4ebc-858f-a80bb805e6fc' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-600d4473-581f-4311-9e3a-1ce3104111c6' class='xr-var-data-in' type='checkbox'><label for='data-600d4473-581f-4311-9e3a-1ce3104111c6' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'></dl></div><div class='xr-var-data'><pre>array([&#x27;2005-01-01T00:00:00.000000000&#x27;, &#x27;2005-01-02T00:00:00.000000000&#x27;,\n",
" &#x27;2005-01-03T00:00:00.000000000&#x27;, ..., &#x27;2020-12-29T00:00:00.000000000&#x27;,\n",
" &#x27;2020-12-30T00:00:00.000000000&#x27;, &#x27;2020-12-31T00:00:00.000000000&#x27;],\n",
" dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-4d17225e-5a88-4d14-989d-b69374d0ecaa' class='xr-section-summary-in' type='checkbox' checked><label for='section-4d17225e-5a88-4d14-989d-b69374d0ecaa' class='xr-section-summary' >Data variables: <span>(4)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span>GW</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-aa875e99-f2f2-4d8f-af33-bf9213d78503' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-aa875e99-f2f2-4d8f-af33-bf9213d78503' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-156effc9-c07e-4a75-b21f-a5e245df0d2c' class='xr-var-data-in' type='checkbox'><label for='data-156effc9-c07e-4a75-b21f-a5e245df0d2c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>GW energy :</span></dt><dd>GW energy without terdiunal tide / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>GW_err</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f99172ed-0778-4ce1-a25b-ca4143b35ed5' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f99172ed-0778-4ce1-a25b-ca4143b35ed5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-889bfa5e-1e6a-42df-84fb-c65399857e4b' class='xr-var-data-in' type='checkbox'><label for='data-889bfa5e-1e6a-42df-84fb-c65399857e4b' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>PW energy error :</span></dt><dd>GW energy error without terdiunal tide / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>PW</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-2d7e81a5-6a24-407b-9838-7c32bb5ada0e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-2d7e81a5-6a24-407b-9838-7c32bb5ada0e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9791d6ef-95e6-4b42-aff3-27df8980358c' class='xr-var-data-in' type='checkbox'><label for='data-9791d6ef-95e6-4b42-aff3-27df8980358c' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>PW energy :</span></dt><dd>PW energy / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>PW_err</span></div><div class='xr-var-dims'>(time, alt)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-7080ce62-3a2a-49d0-b58b-95260053c83e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7080ce62-3a2a-49d0-b58b-95260053c83e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-347660ce-47a8-4157-80dc-e113bc5d6a7d' class='xr-var-data-in' type='checkbox'><label for='data-347660ce-47a8-4157-80dc-e113bc5d6a7d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>PW energy error :</span></dt><dd>PW energy error / m^2/s^2</dd></dl></div><div class='xr-var-data'><pre>[151944 values with dtype=float64]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-6a9ee300-4bd2-4ad0-9e72-548949c399bc' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-6a9ee300-4bd2-4ad0-9e72-548949c399bc' class='xr-section-summary' title='Expand/collapse section'>Attributes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'></dl></div></li></ul></div></div>"
],
"text/plain": [
"<xarray.Dataset>\n",
"Dimensions: (alt: 26, time: 5844)\n",
"Coordinates:\n",
" * alt (alt) int64 70 72 74 76 78 80 82 84 ... 108 110 112 114 116 118 120\n",
" * time (time) datetime64[ns] 2005-01-01 2005-01-02 ... 2020-12-31\n",
"Data variables:\n",
" GW (time, alt) float64 ...\n",
" GW_err (time, alt) float64 ...\n",
" PW (time, alt) float64 ...\n",
" PW_err (time, alt) float64 ..."
]
},
"execution_count": 25,
"metadata": {},
"output_type": "execute_result"
}
],
"source": [
"ds_wo = xr.open_dataset(f'{root_path}t_tide_sa/{station}_solar_cycle.h5', group = 'waves')\n",
"ds_wo = preprocess(ds_wo)\n",
"ds_wo"
]
},
{
"cell_type": "code",
"execution_count": 30,
"metadata": {},
"outputs": [
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3df5xcdX3v8dd784MkRISQCCFBgxhRoBht5CLYSgVqQArILd5Qa9OWNrRixR+3KvVWsS1XemvRai+2Ab2kLUIRg+SBRU1RmoetgEkgmBCRaIAsiUn4HYkkuzuf+8c5s5msu3POzs53Z2b3/Xw8zmNnzpz5nO8sYT/z/a2IwMzMDKCr1QUwM7P24aRgZmb9nBTMzKyfk4KZmfVzUjAzs34TW12AkZg5c2bMmzev1cUwsw6wdu3aJyNiVjNj/tGSQ+P006ayeOl2NTNuKyVLCpK+BJwL7IyIE/NzfwP8BrAP+DHwexHxbP7aFcAlQB/wvoj4ZtE95s2bx5o1axJ9AjMbSyQ91uR4J532xil8//4XufhSTYyI3mbGb5WUzUc3AIsGnFsFnBgRJwE/Aq4AkHQ8sBg4IX/PtZImJCybmdmIXPj26euv/vOZnPmWaSz//BE9rS5PsyRLChGxGnh6wLlv1WTTe4C5+ePzgZsjYm9EbAE2AyenKpuZ2UhIOulnL1Q49Y1T+eClh/J3y55FUkc3x1e1sqP594E788dzgK01r3Xn58zM2s6Fb5++/s8/OAOAWTMnjqnaQkuSgqSPAb3AjdVTg1w26PobkpZKWiNpza5du1IV0cxsULW1hKoPXnoon/3HsVFbGPWkIGkJWQf0u2L/wkvdwNE1l80Ftg32/ohYFhELI2LhrFlNHUhgZlaotpZQNZZqC6OaFCQtAj4CnBcRe2peWgkslnSQpGOA+cB9o1k2M7Mig9USqj70R2OjtpBySOpNwOnATEndwCfIRhsdBKySBHBPRPxRRGyUdAvwEFmz0mUR0ZeqbGZmjbjw7dPXf+DSQwd9rVpbeP+lh/YweJN4R1AnL529cOHC8DwFMytD0tqIWDiC95/066dPW3/nTUOPgdn1ZC9nX7yN+zfsnVQ0byEfdr8GeCIizpU0A/hXYB7wKPDOiHgmv3bY87ga5WUuzMzKueD3Lj6k7gWzZk7kNfMnQzbnqsjlwKaa5x8F7oqI+cBd+fNRn8flpGBmVo6mHFTcKjRliqCg+UjSXODtwPU1p88HluePlwMX1JwftXlcTgpmZiX1RR890Vv3qEQF4MLq0Pn8WDog1GeBDwOVmnNHRMR2gPzny/LzozqPq6N7yc3MRlMFqAw+hapf/uqKiPj4YK9Lqq4Jt1bS6SVuW3oeVzM4KZiZlRRUqBzw5X7wqwqcBpwn6RxgCnCIpH8BdkiaHRHbJc0GdubXl57H1QxuPjIzK6mXoCcqdY/CmkTEFRExNyLmkXUgfzsifptsvtaS/LIlwO3541Gdx+WagplZSRWCvnLNR424GrhF0iXA48BFAKM9j8tJwcyspKB0n0K5eBF3A3fnj58CzhjiuquAq4YRumFOCmZmJfVF0FMw4bfSwROCwUnBzKy0vrTNR23BScHMrKRKQF/BX30nBTOzcSLrUyi+ppM5KZiZldQXoifqL3VRlDTanZOCmVlJfYi+glWxo3NXzQacFMzMSqtAiaTQ2ZwUzMxK6qOLnqi/EESloHmp3TkpmJmVVEH0FawO5D4FM7NxohIqURNwTcHMbFxwn4KZmfXriwn0RP0/m5UOX3zaScHMrKRKiSGp7lMwMxsnKoi+gtFH4dFHZmbjQy9d9DCh7jUVdzSbmY0Plegqrik4KZiZjQ+BCjuSnRTMzMaJ3uhiXxQ0H7lPwcxsfKjQRcXNR2ZmBuWWufDkNTOzcaISoq+gecg1BTOzcaIvuopnNLtPwcxsfKigwnkInV5T6OxFOszMRlGFbJ5CvaNMn4KkKZLuk7Re0kZJn8zPXynpCUkP5Mc5Ne+5QtJmSQ9Leluqz+iawih6yzn/J2n8//i3DyeNbzbeZc1HRUNSS33X3gu8NSJ+JmkS8F1Jd+avfSYiPl17saTjgcXACcBRwL9LenVE9A33MxRJVlOQ9CVJOyVtqDl3UZ4VK5IWDrh+VLKgmVmjqqOP6h1lagqR+Vn+dFJ+1Hvr+cDNEbE3IrYAm4GTR/ZpBpey+egGYNGAcxuAC4HVtScHZMFFwLWS6qdjM7NRlm2y01X3yPsULpS0puZYOjCWpAmSHgB2Aqsi4t78pfdKejD/Yn1Yfm4OsLXm7d35uaZLlhQiYjXw9IBzmyLi4UEuH7UsaGbWqGw/hfpH3hG9IiIW1hzLBsaKiL6IWADMBU6WdCLwBeBYYAGwHfjb/PLBeq+TTIlol47m0llQ0tJq9t21a9eoFM7MDPbvvFbvGO7oo4h4FrgbWBQRO/JkUQGuY/+X427g6Jq3zQW2jfgDDaJdkkLpLBgRy6rZd9asWYmLZWa2X3WZi7rNRyW+v0uaJenQ/PFU4Ezgh5Jm11z2DrImd4CVwGJJB0k6BpgP3NfMz1bVLqOPRi0Lmpk1KpvR3JS1j2YDy/O+0y7gloi4Q9I/S1pA9qX4UeBSgIjYKOkW4CGgF7gsxcgjaJ+ksBL4sqRryIZbJcuCZmaNKjUktURSiIgHgdcPcv7ddd5zFXBVcSlHJllSkHQTcDowU1I38AmyjufPA7OAr0t6ICLeNppZsIyzui5KE3jRG9PEzVV++uqmx3zVt3+36TEBjrjjoCRxq7735Q8ljT/eJft/BFhV+Uqy2CNVQcXLWHiZi8FFxMVDvHTbENePShZspUs/99VksRdPfyZZbDPLlFkl1dtxmpmNE33RRW9B81Gnr33kpGBmVlKppbM7fEMFJwUzs5KiRJ9Ch+cEJwUzs7IiX+ai7jVuPjIzGx96Q/QUJAVvsmNmNk5EPmu57jWuKbTOj9b+hLMm/I9WF8Pa0LUP/1qSuO857jtJ4qaSYu5K5nWJ4sK8L3y6+KIWqdD5Q06LdHRSSKXn1xcWX9SATz10fJK4AItPvilJ3D95/d1J4t56W7otM971V3cki22Zb25bT9eRP2p63HZOCJCvklop2mSns5OGk4KZWUmVKP6j7+YjM7NxIlBh85GHpJqZjRN90UVvQfNRuPnIzGx8yOYpuPnIzMwoN/rIzUdmZuNEqWUuOjwrOCmYmZXUV+mit1Iwec19CmY2Hi2aubT5Qf8y1WS75gg8JLWt9RxxMNuWnNL0uC+8vNL0mAAvSRI186GfviFJ3DkHpdm8Z+XffSZJXICPbDszWexOc86ixWkCd/80Sdhjv9KTJC7AY02I4SGpZmbWr68iNx+ZmVmmUmpIamdzUjAzK8mjj8zMrF9Emeahzm4+qt84ZmZm/fqii96Co8zS2pKmSLpP0npJGyV9Mj8/Q9IqSY/kPw+rec8VkjZLelhSsmWGnRTMzIYhQnWPkp0Ke4G3RsTrgAXAIkmnAB8F7oqI+cBd+XMkHQ8sBk4AFgHXSqq/CFOD3HxkNgy/tPITTY+5b/2hTY9ZdUyyyONTszqaIyKAn+VPJ+VHAOcDp+fnlwN3Ax/Jz98cEXuBLZI2AycD3xveJyjW0Ulh0kv3cdTZjzc97mNPzmh6TICurjTzHwBOmPZEkrivnLwzSdzNPZOSxAV4au/BSeJufmpmkrgpbT37sOKLGtA7NU3cSPkXqQmb5lWii75yQ1IvlHR9zellEbGs9rr8m/5a4FXA/42IeyUdERHbszixXdLL8svnAPfUvL07P9d0HZ0UzMxGVZQeXbQiIj5eN1REH7BA0qHAbZJOrHP5YNWTJOOc6iYFSWW+Mlci4tkmlcfMrG2VWyV1eKOPIuJZSXeT9RXskDQ7ryXMBqpV9W7g6Jq3zQW2DetGJRXVFLblR71POQF4edNKZGbWpso1HxXHkTQL6MkTwlTgTOCvgZXAEuDq/Oft+VtWAl+WdA1wFDAfuK+xT1FfUVLYFBGvr3eBpPubWB4zs7YVJZqPSrbpzAaW5/0KXcAtEXGHpO8Bt0i6BHgcuCi7b2yUdAvwENALXJY3PzVdUVJ4U4kYZa4xM+t8ZSavlVj7KCIeBH7hC3dEPAWcMcR7rgKuKlPMkaibFCLixerjfBLF0bXviYh1tdeYmY1l/XMR6l0zSmVJpdToI0l/Cfwu8GP2f+YA3pqmWGZm7acSoq/itY8A3gkcGxH7UhZmuCZ1VThy6u6mx/1J76ymxwR4dke6HRW2vDxNmX//kB1J4mZzddK4/pjbiy9qwLLDTkoSd3lv8/cEqXph57QkcQ+Z83ySuKcetSVJXIBHPjbyGF77aL8NwLCmXUr6kqSdkjbUnGv5uh5mZo2rv8TFMJa5aFtlk8KngPslfVPSyupR8J4byMbd1mr5uh5mZo2qLnNR7xgv23EuJxtD+wOy+RuFImK1pHkDTrd8XQ8zs4YFxTWBDq8plE0KT0bE55pwvxGv6yFpKbAU4OAj06xxY2Y2FI8+yqyV9CmyWXV7qycjYl2TylF6XY98UallAIe/dlan//7NrIOUmbzW6comheoki9phEo0MSW35uh5mZo2KShdRsMxFmclr7axsUjh74CQ1SYc3cL+Wr+thZtaooGnLXLStsqOPviqpP4FIOhL4Vr03SLqJrKP4OEnd+VoeVwNnSXoEOCt/TkRsBKrrenyDhOt6mJk1LEoeHaxsTeFrwK2S/jtZM89K4H/We0NEXDzES01b12NPzyQe2NH8fSb0+NSmx4SU07XgnqfmJYl749Q0rXivnpxqUhy8YmKazYymqDdJ3JR/Q6Y+kWZk9/OVQ5LEXd13bJK4TRMiCmY0j4ukEBHXSZpMlhzmAZdGxH+lLJiZWbspM6N5TM9TkPTB2qdktYQHgFMknRIR16QsnJlZWxkDzUNFimoKAxfruW2I82Zm44CKRxd1eNIoWjr7k6NVEDOztjcOagp1Rx9JurIoQJlrzMzGhFCJo9WFHJmi5qM/kFRvjVyRLWR3ZdNKZGbWpsrMU+h0RUnhOor7D65rUlnMzNpbBSgckjqGRx+1e59CZd8EXtjy0qbHnfp000MC0Jtw/b4upfn60t0zI0ncp/umJ4kLB66s2Ey3bXtdkri9P0n3u5i2J1Hc7jTzH17sbe8xLIrsGMvKTl4zMzPo+D6DImWXuTAzs4qKjxYlDUkLJI247co1BTOz4WjfTXauB46RtA74T+C/gHsiYlgbapeqKUh6taS7qvstSzpJ0v8abonNzDpakxbDk3S0pO9I2iRpo6TL8/NXSnpC0gP5cU7Ne+ruYx8RC8lWnbgK2Ae8D3hE0npJ15b9iGVrCtcBfwr8Y37zByV9GfirsjcyM+t0CqGC0UcqN/qoF/hQRKyT9BKyjcxW5a99JiI+fUDMA/exPwr4d0mvHriadETsAe6W9H3gXuA04HeARWUKBeWTwrSIuG9Ac1WaJSPNzNpVk2Y059sSV7cm3i1pE0NsQZwr3Mde0m8BpwILyHbIrCaGN0fET8uWrfQezZKOJf91SPrN6gcyM7PGSZpHtrtl9Zv9eyX9DrCGrDbxDOX2sV8G/BD4B2B1RPyokfKUHX10GVnT0WskPQG8H/jjRm5oZtaxYv9chaGOvCZxoaQ1NcfSwcJJmg58FXh/3iH8BeBYsm/724G/rV46eGkO8FJgKTAFuFLSWkl3SPqYpNJbJ5fdT+EnwJmSDga6ImJ32Ruk1NUDU3Y1f/bg5ESfrjI5TVyAnr40k4me652WJO6EiemGaDydaJbgMy+k+V2oL90M2L40+0WxL80eO+0/GTgontGcWRERH693gaRJZAnhxohYARARO2pevw64I39auI993r+wLj/+XtIRwG8CHwD+Aij1R2I4+ynUnq8WwvspmNn4UaZPocT3nXw+wReBTbV/RyXNzvsbAN4BbMgfF+5jL+kksj6F6jGZrMnp82RDVEspu5/CccAb84IB/AawuuxNzMzGgjLLXJSs7JwGvBv4gaQH8nN/BlwsaQFZankUuBSyfewlVfex72XwfexvIPvjfyfw5xHxWLmiHKjU2keSvgW8odpslC+X/ZVGbmhm1rGCbFG8omuKLon4LoPnj3+r856ifez/ApgbETcBSLoXmJW/9uGIuLW4ZOU7ml9ONhmiah/ZXs1mZuNGUSdzixfL+1Pg9prnB5G18JzOMAYGlR2S+s/AfZJuI8uD7wD+qexNzMzGhrbejnNyRGytef7diHgKeCofJFRK2dFHV0m6E/iV/NTvRcT95ctqZjYGtPd2nIfVPomI99Y8nUVJZdc+ejnwJHBbfjyVnzMzGzcUoEr9o4VJ415JfzjwpKRLGTBSqZ6yzUdfZ/9HnQocAzxMtg5Hy6gCExNsItK1r/iaRqQcg71zd5qNWvYenmYh3VRzCQC2v9j8jZcA+sqNTx+2ypR0f0UqP09T5t4jE/1PknDORlM0b/RRCh8AvpYvd7EuP/fLZH0LF5QNUrb56Jdqn0t6A/lQKTOzcaNJ8xRSiIidwKn57OXqF/avR8S3hxOnoa+B+cp+b2zkvWZmnarafNTO8iQwrERQq1RSGDCzuQt4A7Cr0ZuamVl7KltTqN1Nu5esj+GrzS+OmVkba+Pmo2YpmxQeiogDZjBLugjPajazcaRM81GLJ7CNWNkZzVeUPGdmNnY1aTvOdla0SurZwDnAHEmfq3npELzzmpmNN61fyiK5ouajbWS7/5wHrK05v5tsTGxrBUz8efPDdvWm+a9etLfrSOx5Ns3C+Zt3l54IOSzTJ+1NEhdgXyXN3IpKpWzFenjUkyQskO05ksKU6Wn+++3dlm7+StOM5z6FiFgPrJd0Y0Q0rWYg6XLgD8nmeVwXEZ+VNAP4V7KF9h4F3plvQ2dm1hb6Zy2PYXW/+uTrdwPcL+nBgUcjN5R0IllCOBl4HXCupPnAR4G7ImI+cFf+3MysvYzxfoWievbl+c9zm3jP1wL3RMQeAEn/Qbbq6vlkS7wCLAfuBj7SxPuamY1IEzfZaVt1awo128K9JyIeqz2A9zR4zw3Ar0o6XNI0so7so4EjqvfLf75ssDdLWlrdDLv35y80WAQzswZUShwdXlMo23N21iDnzm7khhGxCfhrYBXwDWA9wxjJFBHLImJhRCycOLUDOqXMbExp4012mqKoT+GPJf0AOG5Af8IWoKE+BYCI+GJEvCEifhV4GngE2CFpdn7f2cDORuObmSVRpj+hwxNDUZ/Cl8k2gf4UB3b87o6Ipxu9qaSXRcTOfE+GC4E3kS3HvQS4Ov95e50QZmajbqzUBuopSgoTgOeBywa+IGnGCBLDVyUdDvQAl0XEM5KuBm6RdAnwOHBRg7HNzNIIsn6Doms6WFFSWMv+jziwUz2AVzZy04j4lUHOPQWcMZw4Xb0w9cnmDxree0ia8QMTE/aLT9wxOUncjZNnJ4k7c8bPksQF2Ns7IU3cPZOSxJ2caCMcSLdh1J4n00yWPOTHaSYINosoHl3U6aOPiiavHTNaBTEza3tjoM+gSOn1ACQdBswHplTPRcTqFIUyM2tLZTbZ6fCkUaquJukPgNXAN4FP5j+vTFcsM7M21KSRR5KOlvQdSZskbcyX/kHSDEmrJD2S/zys5j1XSNos6WFJb2vuB9uvbAPe5cAbgcci4teA1+Od18xsnBHF8xRKjk7qBT4UEa8FTgEuk3Q8Qyz3k7+2mGzv5UXAtZKSdJ6VTQovRsSLeeEOiogfAselKJCZWbuqLohX7yhTW4iI7RGxLn+8G9gEzCFb7md5ftly4IL88fnAzRGxNyK2AJvJ1o9rurJJoVvSocDXgFWSbidbVtvMbPwo33x0YXU5nvxYOlRISfPIWl/uZejlfuYAW2ve1p2fa7pSHc0R8Y784ZWSvgO8lGyJCjOzcaPafFTCioj4eGE8aTrZfvfvj4jnpSEHtA72QpIu7WHvRhIR/5GiII3o6gmmPJlgIHakGY8+8cV0I5gnvphmfPdupRmPvjPRXAKA6E00z+SZNJv3THo+SVgAJieaDnLY+jS/i0MeTbjjUDOU6Uwu39k8iSwh3BgRK/LTOyTNjojtA5b76SZbOLRqLolaa9p7poiZWTuJEv0KJZKCsirBF4FNEXFNzUsryZb5gQOX+1kJLJZ0kKRjyKYH3Nesj1UrTbo3MxuLmjd57TTg3cAPJD2Qn/szsrXffmG5n4jYmG969hDZyKXLIqKvKSUZwEnBzKwkRaAYeVaIiO8y9IoYgy73ExFXAVeN+OYFnBTMzEpSiRnNnb6KqpOCmVlZXvvIzMz6lZmx3OFJw0nBzGw4OvyPfpGOTgqVSWLPkc3fR2D61hebHhOga2/praiH7cVZaeYTRFeaORt7etLs/wAlVrFs0ORn08R9SXeiAgPTH0uziceEZ/YkifviKw4rvqiF+peyqKfDk0ZHJwUzs9FUZkbzmN5kx8zMagRQNCTVNQUzs3FiHGyy46RgZlZSqT6FDuekYGY2HB6SamZmAKoEqtT/q68OzwpOCmZmZXnympmZHaDD/+gX6eik0DMddryp+XG7eg5qflBg+h0PJYkLMO3IlxVf1ICuviOTxJ348zST4gC6Es0RPHh7mg1gpjz2TJK4AL2bt6QJfFSafxfPHXtEkrgArBp5CJWoKXhBPDOzcUJR3KfQ6ZwUzMzKauJ2nO3KScHMrKQyzUedzknBzKysCtBXtMxFZ2cNJwUzs2FwR7OZmWUiOr4mUMRJwcysJO+n0Oa6Jvcx5eW7mx73qRNe2vSYANNW7EsSF0DPNf/3ADBla5rNe7r2vSRJXICuvX1J4k7cuitJ3L4daeICEGlWb+s9elaSuC8clSRs04hABTWFTt9PoasVN5X0AUkbJW2QdJOkKZJmSFol6ZH8Z3tvwWRm40+lxNHhNYVRTwqS5gDvAxZGxInABGAx8FHgroiYD9yVPzczayuKqHuUyQqSviRpp6QNNeeulPSEpAfy45ya166QtFnSw5LeluaTZVpSUyBrtpoqaSIwDdgGnA8sz19fDlzQorKZmQ2uEsVHuZrCDcCiQc5/JiIW5Me/AUg6nuyL8wn5e66VNKE5H+gXjXpSiIgngE8DjwPbgeci4lvAERGxPb9mOzDoYj6SlkpaI2lN7/NpNg83MxtMdY/mekeZpBARq4GnS972fODmiNgbEVuAzcDJDX6EQq1oPjqM7EMeAxwFHCzpt8u+PyKWRcTCiFg48ZBpqYppZvaLqkNS6x1ZVriw+uU1P5aWvMN7JT2YNy9V+1XnAFtrrunOzyXRiuajM4EtEbErInqAFcCpwA5JswHynztbUDYzsyGpEqiv4MhqCiuqX17zY1mJ8F8AjgUWkLWi/G31toNcm6w7uxVJ4XHgFEnTJAk4A9gErASW5NcsAW5vQdnMzIYWJY9GQkfsiIi+iKgA17G/iagbOLrm0rlk/bBJjPo8hYi4V9KtwDqgF7gfWAZMB26RdAlZ4rioKFalr4s9zzV/HP2EQ9Ik4YlHzU4SF4CD0zSlxaQ0/Vld+9LMJQDo+nmafQ9SzWTtOmR6krgAmpRm34o9h01OEjeSdZ82SVA4T6FRkmZX+1WBdwDVkUkrgS9LuoasyX0+cF+SQtCiyWsR8QngEwNO7yWrNZiZtakyy1yUGpJ6E3A6MFNSN9nfw9MlLcgDPApcChARGyXdAjxE9kX6sohI9q2qo2c0m5mNpmqfQt1ryo0+uniQ01+sc/1VwFXFkUfOScHMrKyguKbQ4QvmOSmYmZXlVVLNzKyqTPNRp6995KRgZlaWm4/MzGw/Nx+ZmVlVJUrs0Tw6RUmlo5PCLx1+JGuWfKTpcV/3/s80PSZAZfbhSeIC7JsxJUncyoQ0W4bEpHRbkUz4eZoJW5OSTeTrTRIXIHrSDGef9tBPk8Q99iG4c8s1SWLrYx8ceYwyk9c6vCbR0UnBzGxUefSRmZn1i3zPhKJrOpiTgplZWRFQKdj3urNzgpOCmVlppZqPOjsrOCmYmZVVKdN8NDpFScVJwcysrKhApWBEl/sUzMzGiaC4ptDhVQUnhUGs/+wHksQ9683pVr59YXaasfkq6FNrVCXhv7yu3jTzCab3pvllVBKVF2DSj9PMJwC484nPJ4vdtsr0KXR2TnBSMDMrzfMUzMysX6UCfUV9Comq16PEScHMrCw3H5mZ2QHcfGRmZgBRqRBuPjIzM8CT18zMrJZHH5mZWVWp0UednTScFEbTfz2QLPQLbz41SdyufUnCsndGmrgAfdPS/E/ZN/mgJHEB1lw/8g1gLL2oBFG4Smrxvz9JXwLOBXZGxIn5uRnAvwLzgEeBd0bEM/lrVwCXAH3A+yLim41+hiJdqQKbmY051SGpRUexG4BFA859FLgrIuYDd+XPkXQ8sBg4IX/PtZKSTYN3UjAzK6u6yU69o1SYWA08PeD0+cDy/PFy4IKa8zdHxN6I2AJsBk4e+YcZnJOCmVlJ1SGpdY8sMVwoaU3NsbRE+CMiYjtA/vNl+fk5wNaa67rzc0m4T8HMrKyIEvMQAmBFRHy8SXfVUDdJwTUFM7OSIip5Z3Odo/HRRzskzQbIf+7Mz3cDR9dcNxfY1vCHKOCkYGZWzuqdsTWrKQxx9FX2sZtnAH7UQPyVwJL88RLg9prziyUdJOkYYD5w30g/zFCcFMzMyrl7Ly/yQjw/5AXd/JjZvIKI2FMvkKSbgO8Bx0nqlnQJcDVwlqRHgLPy50TERuAW4CHgG8BlEVEwWaJxGkFVp7EbSseRjcWteiXwceCfGGKM7lAWLlwYa9asSVNQMxtTJK2NiIUjifHLektsYwsn6r/9wmt90cv3+Q4/47mDi5JCOxv1mkJEPBwRCyJiAfDLwB7gNoYYo2tm1i7WsbprqNpC2VpCu2t189EZwI8j4jGGHqNrZtYWIiKO4bVsYdMB5/uil+08ziM8eHCLitY0rU4Ki4Gb8sdDjdE1M2sbg9UWxkotAVqYFCRNBs4DvjLM9y2tTgjZtWtXmsKZmQ1hYG1hLNUSoLU1hbOBdRGxI38+1BjdA0TEsohYGBELZ82aNUpFNTPbr7a2MJZqCdDapHAx+5uOYOgxumZmbaVaW9jMhjFVS4AWJQVJ08jG4a6oOT3oGF0zs5s6I8UAAASmSURBVHa0jtVdvfSMqVoCtGCeQjNJ2gU81upyADOBJ1tdiGHotPKCyzwaOq28MLwyvyIi3OZcoKOTQruQtGakk2JGU6eVF1zm0dBp5YXOLHO7a/WQVDMzayNOCmZm1s9JoTmWtboAw9Rp5QWXeTR0WnmhM8vc1tynYGZm/VxTMDOzfk4KZmbWz0mhQZKOlvQdSZskbZR0eavLVJakCZLul3RHq8tShqRDJd0q6Yf57/tNrS5TPZI+kP+b2CDpJklTWl2mgSR9SdJOSRtqzs2QtErSI/nPw1pZxoGGKPPf5P8uHpR0m6RDW1nGscBJoXG9wIci4rXAKcBlko5vcZnKuhwGrP3b3v4O+EZEvAZ4HW1cdklzgPcBCyPiRGAC2WrA7eYGYNGAc+2+p8kN/GKZVwEnRsRJZFtgXjHahRprnBQaFBHbI2Jd/ng32R+qOa0tVTFJc4G3A9e3uixlSDoE+FXgiwARsS8inm1tqQpNBKZKmghMI+Em642KiNXA0wNOt/WeJoOVOSK+FRG9+dN7yDa1txFwUmgCSfOA1wP3trYkpXwW+DBQaXVBSnolsAv4f3mT1/WS2nbxsYh4Avg08DiwHXguIr7V2lKV1ul7mvw+cGerC9HpnBRGSNJ04KvA+yPq7OjdBiSdC+yMiLWtLsswTATeAHwhIl4PvED7NWv0y9vhzweOAY4CDpb0260t1dgn6WNkTbo3trosnc5JYQQkTSJLCDdGxIqi69vAacB5kh4FbgbeKulfWlukQt1Ad0RUa2G3kiWJdnUmsCUidkVED9lKwKe2uExlldrTpN1IWgKcC7wrPPFqxJwUGiRJZO3cmyLimlaXp4yIuCIi5kbEPLLOz29HRFt/i42InwJbJR2XnzoDeKiFRSryOHCKpGn5v5EzaOOO8QE6bk8TSYuAjwDnjaXlq1vJSaFxpwHvJvu2/UB+nNPqQo1RfwLcKOlBYAHwv1tcniHlNZpbgXXAD8j+H2u7pRgk3QR8DzhOUrekS2jzPU2GKPPfAy8BVuX/D/5DSws5BniZCzMz6+eagpmZ9XNSMDOzfk4KZmbWz0nBzMz6OSmYmVk/JwUbN/LVVt9T8/z0Tlkp1my0OCnYeHIo8J7Cq8zGMScFa0uS5uXr5F+f70two6QzJf1nvt7/yfn6/1/L19K/R9JJ+XuvzNfev1vSTyS9Lw97NXBsPsnpb/Jz02v2argxn4VsNm5NbHUBzOp4FXARsBT4PvBbwJuB84A/A7YC90fEBZLeCvwT2YxngNcAv0Y22/VhSV8gW0jvxIhYAFnzEdnqtieQLW/9n2Qz1b87Gh/OrB25pmDtbEtE/CAiKsBGsg1ggmz5iHlkCeKfASLi28Dhkl6av/frEbE3Ip4kW9jtiCHucV9EdOf3eCCPazZuOSlYO9tb87hS87xCVssdrKmnum5L7Xv7GLpWXPY6s3HBScE62WrgXdDfFPRkwZ4Wu8mak8xsCP5WZJ3sSrId2R4E9rB/2edBRcRTeUf1BrIdur6evohmncWrpJqZWT83H5mZWT8nBTMz6+ekYGZm/ZwUzMysn5OCmZn1c1IwM7N+TgpmZtbv/wN+iO83LNU/OgAAAABJRU5ErkJggg==\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"p = ds['GW'].groupby('time.month').mean('time').plot(x = 'month',\n",
" robust = True)"
]
},
{
"cell_type": "code",
"execution_count": 36,
"metadata": {},
"outputs": [],
"source": [
"vmax = p.colorbar.vmax\n",
"vmin = p.colorbar.vmin"
]
},
{
"cell_type": "code",
"execution_count": 37,
"metadata": {},
"outputs": [
{
"data": {
"text/plain": [
"<matplotlib.collections.QuadMesh at 0x7f9635f15b00>"
]
},
"execution_count": 37,
"metadata": {},
"output_type": "execute_result"
},
{
"data": {
"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEGCAYAAACKB4k+AAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3de5xcZZ3n8c+3cw9RAiZgIGgQIwosBo0MwlxQYAzIAGYHN4wX1mEmOOKCyI6KroizmxVXRWecxZmALNHBIGKQrA5KRDAvHS4m4ZYQkUu4BDJJuCkXSdJdv/njnOpUmu46T1fX6arq/r5fr/PqqlOnnvNUaPpXz+33KCIwMzMD6Gp1BczMrH04KJiZWS8HBTMz6+WgYGZmvRwUzMys19hWV2Aopk2bFrNmzWp1NcysA6xevfrJiJjezDI/fPrUOPqoSSxYuEnNLLeVSgsKki4HTgS2RMQh+bkvAX8GbAceBD4UEc/mr50PnAH0AGdHxE+K7jFr1ixWrVpV0icws5FE0iNNLu/Qo942kV/d8RKnnamxEdHdzPJbpczuoyuAeX3OrQAOiYhDgd8A5wNIOghYABycv+cSSWNKrJuZ2ZDMf/eUuy767DSO/ZPJLPn63jtaXZ9mKS0oRMRK4Ok+526oiaa3AjPzxycDV0XEtojYADwAHF5W3czMhkLSoc+/UOHIt03i42dO5e8XP4ukju6Or2rlQPNfAtfnj/cFHqt5bWN+zsys7cx/95S7PvvxPQGYPm3siGottCQoSPoM0A1cWT3Vz2X95t+QtFDSKkmrtm7dWlYVzcz6VdtKqPr4mVP52j+PjNbCsAcFSaeTDUC/L3YmXtoI7Fdz2Uzgif7eHxGLI2JuRMydPr2pEwnMzArVthKqRlJrYViDgqR5wCeBkyLixZqXlgMLJE2QtD8wG7h9OOtmZlakv1ZC1XkfHhmthTKnpC4FjgamSdoIfI5sttEEYIUkgFsj4sMRsU7S1cC9ZN1KZ0VET1l1MzNrxPx3T7nr3DOn9vtatbXwsTOn7qD/LvGOoE5OnT137tzwOgUzSyFpdUTMHcL7D/3Toyffdf3SgefAbH2ym+NPe4I71m4bV7RuIZ92vwp4PCJOlLQn8F1gFvAw8N6IeCa/dtDruBrlNBdmZmlO+dBpr6x7wfRpY3nj7PGQrbkqcg6wvub5p4AbI2I2cGP+fNjXcTkomJml0cQJxb1CEycKCrqPJM0E3g1cVnP6ZGBJ/ngJcErN+WFbx+WgYGaWqCd62BHddY9KVADmV6fO58fCPkV9DfgEUKk5t3dEbALIf+6Vnx/WdVwdPUpuZjacKkCl/yVUvfJXl0XEBf29LqmaE261pKMTbpu8jqsZHBTMzBIFFSq7fLnv/6oCRwEnSToBmAi8UtK/AJslzYiITZJmAFvy65PXcTWDu4/MzBJ1E+yISt2jsCURcX5EzIyIWWQDyD+LiPeTrdc6Pb/sdOC6/PGwruNyS8HMLFGFoCet+6gRFwFXSzoDeBQ4FWC413E5KJiZJQqSxxTSyou4Gbg5f/wUcMwA1y0CFg2i6IY5KJiZJeqJYEfBgt9KBy8IBgcFM7NkPeV2H7UFBwUzs0SVgJ6Cv/oOCmZmo0Q2plB8TSdzUDAzS9QTYkfUT3VRFDTanYOCmVmiHkRPQVbs6Nys2YCDgplZsgokBIXO5qBgZpaohy52RP1EEJWC7qV256BgZpaogugpyA7kMQUzs1GiEkpoCbilYGY2KnhMwczMevXEGHZE/T+blQ5PPu2gYGaWqJIwJdVjCmZmo0QF0VMw+yg8+8jMbHToposdjKl7TcUDzWZmo0MluopbCg4KZmajQ6DCgWQHBTOzUaI7utgeBd1HHlMwMxsdKnRRcfeRmZlBWpoLL14zMxslKiF6CrqH3FIwMxsleqKreEWzxxTMzEaHCipch9DpLYXOTtJhZjaMKmTrFOodKWMKkiZKul3SXZLWSfp8fv5CSY9LujM/Tqh5z/mSHpB0n6R3lfUZ3VIYRsd1nVpq+Ssq3yu1fLPRLus+KpqSmvRdexvwzoh4XtI44BeSrs9f+2pEfLn2YkkHAQuAg4F9gJ9KekNE9Az2MxQpraUg6XJJWyStrTl3ah4VK5Lm9rl+WKKgmVmjqrOP6h0pLYXIPJ8/HZcf9d56MnBVRGyLiA3AA8DhQ/s0/Suz++gKYF6fc2uB+cDK2pN9ouA84BJJ9cOxmdkwyzbZ6ap75GMK8yWtqjkW9i1L0hhJdwJbgBURcVv+0kcl3Z1/sd4jP7cv8FjN2zfm55qutKAQESuBp/ucWx8R9/Vz+bBFQTOzRmX7KdQ/8oHoZRExt+ZY3LesiOiJiDnATOBwSYcA3wAOAOYAm4Cv5Jf3N3pdypKIdhloTo6CkhZWo+/WrVuHpXJmZrBz57V6x2BnH0XEs8DNwLyI2JwHiwpwKTu/HG8E9qt520zgiSF/oH60S1BIjoIRsbgafadPn15ytczMdqqmuajbfZTw/V3SdElT88eTgGOBX0uaUXPZe8i63AGWAwskTZC0PzAbuL2Zn62qXWYfDVsUNDNrVLaiuSm5j2YAS/Kx0y7g6oj4oaRvS5pD9qX4YeBMgIhYJ+lq4F6gGzirjJlH0D5BYTnwHUkXk023Ki0Kmpk1KmlKakJQiIi7gcP6Of+BOu9ZBCwqruXQlBYUJC0FjgamSdoIfI5s4PnrwHTgR5LujIh3DWcUTPGOP/1iKeWWHYH/4IMXN73MPf7/vU0vE+A3FxxUSrlVD5378VLLH+3KXHPTzuttKqg4jYXTXPQvIk4b4KVrB7h+WKJgK4193azSyn76yBnFF5nZkKRkSfV2nGZmo0RPdNFd0H3U6bmPHBTMzBIlpc7u8A0VHBTMzBJFwphCh8cEBwUzs1SRp7moe427j8zMRofuEDsKgoI32TEzGyUiX7Vc9xq3FFrnN6sfKme+9LFzi6+xtlbWPPp2nkPfn7ef9pXiixowpZRSM2XvOzIUFTp/ymmRjg4KZRn701WllNtdSqmZyQdMK6VcjS3nV+SA824ppVwbHs//lyO4Zel5TS+3nQMC5FlSK0Wb7HR20HBQMDNLVIniP/ruPjIzGyUCFXYfeUqqmdko0RNddBd0H4W7j8zMRodsnYK7j8zMjLTZR+4+MjMbJZLSXHR4VHBQMDNL1FPportSsHjNYwpmNhq1+5qCMgSekmpN1DVpUmllT7jl1+UUvPsrSym2662HlFIuQM8k/1pXTfnura2uwojiKalmZtarpyJ3H5mZWaaSNCW1szkomJkl8uwjMzPrFZHSPdTZ3Uf1O8fMzKxXT3TRXXCkpNaWNFHS7ZLukrRO0ufz83tKWiHp/vznHjXvOV/SA5Luk/Susj6jg4KZ2SBEqO6ROKiwDXhnRLwZmAPMk3QE8CngxoiYDdyYP0fSQcAC4GBgHnCJpPpJmBrk7iOzQTh+1rlNL7P70Y1NL9PK0ayB5ogI4Pn86bj8COBk4Oj8/BLgZuCT+fmrImIbsEHSA8DhQNM3JunsoCChseOaXmx072h6mQCV3/++lHIBug47qJRy9fty/i1Y/1A55QJd28up85h99i6l3E7UNXlyKeVq4oRSygXgqaEXUYkuetKmpM6XdFnN6cURsbj2uvyb/mrg9cD/jYjbJO0dEZuycmKTpL3yy/cFahedbMzPNV1nBwUzs+EUybOLlkXEBXWLiugB5kiaClwrqd6Kzv6aJ6XMc6obFCTtmVBGJSKebVJ9zMzaVlqW1MHNPoqIZyXdTDZWsFnSjLyVMAPYkl+2Ediv5m0zgScGdaNERS2FJ/Kj3qccA7ymaTUyM2tTad1HxeVImg7syAPCJOBY4IvAcuB04KL853X5W5YD35F0MbAPMBu4vbFPUV9RUFgfEYfVu0DSHU2sj5lZ24qE7qPEPp0ZwJJ8XKELuDoifijpFuBqSWcAjwKnZveNdZKuBu4FuoGz8u6npisKCm9PKCPlGjOzzpeyeC0h91FE3A287At3RDwFHDPAexYBi1KqORR1g0JEvFR9nC+i2K/2PRGxpvYaM7ORrHctQr1rhqkuZUmafSTpfwL/FXiQnZ85gHeWUy0zs/ZTCdFTce4jgPcCB0TE9jIrM1iS6Jo0senl9jxX0tz8EnW9WNJ/mh3d5ZRb4v85Y169V/FFDShrkZnGlLIwFYDoKaXbma7pryql3J5X71F8UaOasE7BuY92WgtMHUzBki6XtEXS2ppzLc/rYWbWuPopLgaR5qJtpQaFLwB3SPqJpOXVo+A9V5DNu63V8rweZmaNqqa5qHeMlu04l5DNob2HbP1GoYhYKWlWn9Mtz+thZtawoLgl0OEthdSg8GRE/EMT7jfkvB6SFgILASZqtyZUycwsnWcfZVZL+gLZqrpt1ZMRsaZJ9UjO65EnlVoMsPuYaZ3+729mHSRl8VqnSw0K1UUWR9Sca2RKasvzepiZNSoqXURBmouUxWvtLDUoHN93kZqkRuaktTyvh5lZo4KmpbloW6mzj74vqTeASHo1cEO9N0haSjZQfKCkjXkuj4uA4yTdDxyXPyci1gHVvB4/psS8HmZmDYvEo4OlthR+AFwj6T+TdfMsB/57vTdExGkDvNS8vB4RxLZtxde1ibI2JgFAndVkLXPDodj076WVXYayFpiVqfvRx0spd+yL5f1eNEWIKFjRPCqCQkRcKmk8WXCYBZwZEf9WZsXMzNpNyormEb1OQdLHa5+StRLuBI6QdEREXFxm5czM2soI6B4qUtRSeEWf59cOcN7MbBRQ8eyiDg8aRamzPz9cFTEza3ujoKVQd/aRpAuLCki5xsxsRAglHK2u5NAUdR/9laTf1XldZInsLmxajczM2lTKOoVOVxQULqV4/ODSJtXFzKy9VYDCKakjePZRu48pRASV7W21709dXbu/srSyX5q5eynljn+mnN1WuyY0f3OkXpGUyHfwxXbgeoKy/p0r28r5veje+mQp5TaLIjtGstTFa2ZmBh0/ZlAkNc2FmZlVVHy0KGhImiMNPbWBWwpmZoPRvpvsXAbsL2kN8Evg34BbI6LeZKGXSWopSHqDpBur+y1LOlTS/xhsjc3MOlqTkuFJ2k/STZLWS1on6Zz8/IWSHpd0Z36cUPOeuvvYR8RcsqwTi4DtwNnA/ZLuknRJ6kdMbSlcCvwt8M/5ze+W9B3gf6XeyMys0ymECmYfKW32UTdwXkSskfQKso3MVuSvfTUivrxLmbvuY78P8FNJb+ibTToiXgRulvQr4DbgKOCDwLyUSkF6UJgcEbf36a7qTr2JmdmI0KQVzfm2xNWtiZ+TtJ4BtiDOFe5jL+kvgCOBOWQ7ZFYDwx9GRHLq4OQ9miUdQP7PIenPqx/IzMwaJ2kW2e6W1W/2H5X0QWAVWWviGdL2sV8M/Br4J2BlRPymkfqkzj46i6zr6I2SHgc+BvxNIzc0M+tYsXOtwkBH3pKYL2lVzbGwv+IkTQG+D3wsHxD+BnAA2bf9TcBXqpf2X5td7A4sBCYCF0paLemHkj4jKXnr5NT9FB4CjpW0G9AVEc+l3sBqTCpvwdaOV4wppdzoKqfOk2bOKKVcALpL6tnc+lQpxcb2HaWUC6DZry2l3LE7SlrI9+Qz5ZQL0Ix1cUHxiubMsoi4oN4FksaRBYQrI2IZQERsrnn9UuCH+dPCfezz8YU1+fGPkvYG/hw4F/g7IOmPxGD2U6g9X62E91Mws9EjZUwhYcwhX0/wTWB97d9RSTPy8QaA9wBr88eF+9hLOpRsTKF6jCfrcvo62RTVJKn7KRwIvC2vGMCfAStTb2JmNhKkpLlIXD12FPAB4B5Jd+bnPg2cJmkOWWh5GDgTsn3sJVX3se+m/33sryD743898NmIeCStKrtKyn0k6QbgLdVuozxd9vcauaGZWccKsqR4RdcUXRLxC/qPH/9a5z1F+9j/HTAzIpYCSLoNmJ6/9omIuKa4ZukDza8hWwxRtZ1sr2Yzs1GjaJC5xcny/ha4rub5BLIenqMZxMSg1Cmp3wZul3QtWRx8D/Ct1JuYmY0Mbb0d5/iIeKzm+S8i4ingqXySUJLU2UeLJF0P/FF+6kMRcUd6Xc3MRoD23o5zj9onEfHRmqfTSZSa++g1ZBO6rs2Pp/JzZmajhgJUqX+0MGjcJumv+56UdCZ9ZirVk9p99CN2ftRJwP7AfWR5OFpGUimbiGjihKaXCVCZmtyCaxvdu5Wz/qFnjymllAugSjmb7GjKpHLKHXq24wFtn17Ov/ML+4wvpdwpjxVt9DgEP29CGc2bfVSGc4Ef5Oku1uTn3ko2tnBKaiGp3Uf/qfa5pLeQT5UyMxs1mrROoQwRsQU4Ml+9XP3C/qOI+NlgymloP4U8s9/bGnmvmVmnqnYftbM8CAwqENRKCgp9VjZ3AW8BtjZ6UzMza0+pLYXajr5usjGG7ze/OmZmbayNu4+aJTUo3BsRu6xglnQqXtVsZqNISvdRixewDVnqiubzE8+ZmY1cTdqOs50VZUk9HjgB2FfSP9S89Eq885qZjTatT2VRuqLuoyfIdv85CVhdc/45sjmxrdXVVcqaAk3bs+llAux4RTnrH4CO+4YS41IbqYNXGVPS3hIl/ffr2lbS3gRAz8Ry/p2f36ec2fiTN5fz366pRvOYQkTcBdwl6cqIaFrLQNI5wF+TrfO4NCK+JmlP4LtkifYeBt6bb0NnZtYWelctj2B1v0bk+bsB7pB0d9+jkRtKOoQsIBwOvBk4UdJs4FPAjRExG7gxf25m1l5G+LhCUffROfnPE5t4zzcBt0bEiwCSfk6WdfVkshSvAEuAm4FPNvG+ZmZD0sRNdtpW3ZZCzbZwH4mIR2oP4CMN3nMt8MeSXiVpMtlA9n7A3tX75T/36u/NkhZWN8PeXvl9g1UwM2tAJeHo8JZC6ijUcf2cO76RG0bEeuCLwArgx8BdDGImU0Qsjoi5ETF3fFc5CcrMzAbSxpvsNEXRmMLfSLoHOLDPeMIGoKExBYCI+GZEvCUi/hh4Grgf2CxpRn7fGcCWRss3MytFynhChweGojGF75BtAv0Fdh34fS4inm70ppL2iogt+Z4M84G3k6XjPh24KP95XZ0izMyG3UhpDdRTFBTGAL8Dzur7gqQ9hxAYvi/pVcAO4KyIeEbSRcDVks4AHgVObbBsM7NyBNm4QdE1HawoKKxm50fsO6gewOsauWlE/FE/554CjhlUOT099Pz2t41Uoa4xr92n6WUCxNjy5iWMe6GcBVCqlPMb3rW9vAVblXElbQw0qZyFYGX9G0N5v3PjXiylWCY+9GQ5BTeJKJ5d1Omzj4oWr+0/XBUxM2t7I2DMoEjyJjuS9gBmA737X0bEyjIqZWbWllI22enwoJHUHpb0V8BK4CfA5/OfF5ZXLTOzNtSkmUeS9pN0k6T1ktblqX+QtKekFZLuz3/uUfOe8yU9IOk+Se9q7gfbKbWT9BzgbcAjEfEO4DC885qZjTKieJ1C4uykbuC8iHgTcARwlqSDGCDdT/7aArK9l+cBl0gqZfAsNSi8FBEv5ZWbEBG/Bg4so0JmZu2qmhCv3pHSWoiITRGxJn/8HLAe2Jcs3c+S/LIlwCn545OBqyJiW0RsAB4gyx/XdKlBYaOkqcAPgBWSriNLq21mNnqkdx/Nr6bjyY+FAxUpaRZZ78ttDJzuZ1/gsZq3bczPNV3SQHNEvCd/eKGkm4DdyVJUmJmNGtXuowTLIuKCwvKkKWT73X8sIn4nDTihtb8XShnSTp591FuLiJ+XUZF2oqeav/YBYNyY8mYwq6ecKQ+V8YP+FUmi7eVt3NfVXU7C+zElzfkvMz9/1/Zyfi+m3VXOQoXuhx4updymSRlMTh9sHkcWEK6MiGX56c2SZkTEpj7pfjaSJQ6tmklJvTXlbX9lZjbSRMK4QkJQUNYk+CawPiIurnlpOVmaH9g13c9yYIGkCZL2J1secHuzPlatcr4GmpmNRM1bvHYU8AHgHkl35uc+TZb77WXpfiJiXb7p2b1kM5fOiohS0gI4KJiZJVIEiqFHhYj4BQNnxOg33U9ELAIWDfnmBRwUzMwSKWFFc6dnUXVQMDNL5dxHZmbWK2XFcocHDQcFM7PB6PA/+kUcFPrR/Xg5i7XHvFBSEnpAE8aXUu6YceNKKZex5f3qqauk9QRNGGDst9wXt5VSLsDYB58vpdzuzaNzt9zeVBb1dHjQcFAwM0uUsqJ5RG+yY2ZmNQIoajG6pWBmNkqMgk12HBTMzBIljSl0OAcFM7PB8JRUMzMDUCVQpf5ffXV4VHBQMDNL5cVrZma2iw7/o1+ko4OCxo1j7F4zml5u9xObml4mQM+zz5ZSbifqmjCxtLKjp5SMwqWJ7h2troIlUkJLwQnxzMxGCUXxmEKnc1AwM0vVxO0425WDgplZopTuo07noGBmlqoC9BSluejsqOGgYGY2CB5oNjOzTETHtwSKOCiYmSXyfgptLiaMZfvrX930csfsM63pZQLEqntKKbcTVba91OoqWBsa+7pZ5RX+4NCLEFG42VKn76fQ1YqbSjpX0jpJayUtlTRR0p6SVki6P/+5RyvqZmY2oErC0eEthWEPCpL2Bc4G5kbEIcAYYAHwKeDGiJgN3Jg/NzNrK4qoe6REBUmXS9oiaW3NuQslPS7pzvw4oea18yU9IOk+Se8q55NlWtJSIOu2miRpLDAZeAI4GViSv74EOKVFdTMz618lio+0lsIVwLx+zn81Iubkx78CSDqI7Ivzwfl7LpE0pjkf6OWGPShExOPAl4FHgU3AbyPiBmDviNiUX7MJ2Ku/90taKGmVpFXbd7wwXNU2M+vdo7nekRIUImIl8HTibU8GroqIbRGxAXgAOLzBj1CoFd1He5B9yP2BfYDdJL0/9f0RsTgi5kbE3PHjdiurmmZmL1edklrvyKLC/OqX1/xYmHiHj0q6O+9eqo6r7gs8VnPNxvxcKVrRfXQssCEitkbEDmAZcCSwWdIMgPznlhbUzcxsQKoE6ik4spbCsuqX1/xYnFD8N4ADgDlkvShfqd62n2tLG85uRVB4FDhC0mRJAo4B1gPLgdPza04HrmtB3czMBhaJRyNFR2yOiJ6IqACXsrOLaCOwX82lM8nGYUsx7OsUIuI2SdcAa4Bu4A5gMTAFuFrSGWSB49TCsrpEZULzx1sq48sZw+noRSFmNbrGjy+l3MrUNu8SDgrXKTRK0ozquCrwHqA6M2k58B1JF5N1uc8Gbi+lErTo71REfA74XJ/T28haDWZmbSolzUXSlNSlwNHANEkbyf4eHi1pTl7Aw8CZABGxTtLVwL1kX6TPiojSdpLyl1czs0TVMYW616TNPjqtn9PfrHP9ImBRcclD56BgZpYqKG4pdHjCPAcFM7NUzpJqZmZVKd1HnZ77yEHBzCyVu4/MzGwndx+ZmVlVJRL2aB6eqpSlo4PCgbNfzU0/+WTTyz36+P/T9DLNRpLK9u3lFLxmHSsq3yulaOnzQy8jZfFah7ckOjoomJkNK88+MjOzXpHvmVB0TQdzUDAzSxUBlUrBNcNTlbI4KJiZpUrqPursqOCgYGaWqpLSfTQ8VSmLg4KZWaqoQKUgQanHFMzMRomguKXQ4U0FB4V+3Hz9J0op97iuwn2DzDpGWesJ2lrKmEJnxwQHBTOzZF6nYGZmvSoV6CkaUyiYstrmHBTMzFK5+8jMzHbh7iMzMwOISoVw95GZmQFevGZmZrU8+8jMzKqSZh91dtBwUDBrA6NyIVgHikoQhVlSi4OCpMuBE4EtEXFIfm5P4LvALOBh4L0R8Uz+2vnAGUAPcHZE/KTRz1Ckq6yCzcxGnOqU1KKj2BXAvD7nPgXcGBGzgRvz50g6CFgAHJy/5xJJY5r1kfpyUDAzS1XdZKfekVRMrASe7nP6ZGBJ/ngJcErN+asiYltEbAAeAA4f+ofpn4OCmVmi6pTUukcWGOZLWlVzLEwofu+I2ASQ/9wrP78v8FjNdRvzc6XwmIKZWaqIhHUIAbAsIi5o0l010E3K4JaCmVmiiEo+2FznaHz20WZJMwDyn1vy8xuB/Wqumwk80fCHKOCgYGaWZuWWeCxrKQxw9FS28xzPAPymgfKXA6fnj08Hrqs5v0DSBEn7A7OB24f6YQbioGBmlubmbbzEC/G7AS/YyIPM4LVExIv1CpK0FLgFOFDSRklnABcBx0m6Hzguf05ErAOuBu4FfgycFREFiyUapyE0dRq7oXQg2VzcqtcBFwDfYoA5ugOZO3durFq1qpyKmtmIIml1RMwdShlv1Z/EE2zgEP3By17riW5+xU08z293KwoK7WzYWwoRcV9EzImIOcBbgReBaxlgjq6ZWbtYw8qugVoLqa2Edtfq7qNjgAcj4hEGnqNrZtYWIiL2501sYP0u53uim008yv3cvVuLqtY0rQ4KC4Cl+eOB5uiambWN/loLI6WVAC0MCpLGAycBg0r6ImlhdUHI1q1by6mcmdkA+rYWRlIrAVrbUjgeWBMRm/PnA83R3UVELI6IuRExd/r06cNUVTOznWpbCyOplQCtDQqnsbPrCAaeo2tm1laqrYUHWDuiWgnQoqAgaTLZPNxlNaf7naNrZtaO1rCyq5sdI6qVAC1Yp9BMkrYCj7S6HsA04MlWV2IQOq2+4DoPh06rLwyuzq+NCPc5F+jooNAuJK0a6qKY4dRp9QXXeTh0Wn2hM+vc7lo9JdXMzNqIg4KZmfVyUGiOxa2uwCB1Wn3BdR4OnVZf6Mw6tzWPKZiZWS+3FMzMrJeDgpmZ9XJQaJCk/STdJGm9pHWSzml1nVJJGiPpDkk/bHVdUkiaKukaSb/O/73f3uo61SPp3Px3Yq2kpZImtrpOfUm6XNIWSWtrzu0paYWk+/Ofe7Syjn0NUOcv5b8Xd0u6VtLUVtZxJHBQaFw3cF5EvAk4AjhL0kEtrlOqc6BP7t/29vfAjyPijcCbaeO6S9oXOBuYGxGHAGPIsgG3myuAeX3OtfueJlfw8jqvAA6JiEPJtsA8f7grNdI4KDQoIjZFxJr88XNkf6j2bW2tikmaCbwbuKzVdUkh6ZXAHwPfBIiI7RHxbGtrVWgsMEnSWGAyJW6y3qiIWAk83ed0W+9p0l+dI+KGiOjOn/XmGoYAAAMYSURBVN5Ktqm9DYGDQhNImgUcBtzW2pok+RrwCaDS6ookeh2wFfh/eZfXZZLaNvlYRDwOfBl4FNgE/DYibmhtrZJ1+p4mfwlc3+pKdDoHhSGSNAX4PvCxiDo7ercBSScCWyJidavrMghjgbcA34iIw4AXaL9ujV55P/zJwP7APsBukt7f2lqNfJI+Q9ale2Wr69LpHBSGQNI4soBwZUQsK7q+DRwFnCTpYeAq4J2S/qW1VSq0EdgYEdVW2DVkQaJdHQtsiIitEbGDLBPwkS2uU6qkPU3ajaTTgROB94UXXg2Zg0KDJImsn3t9RFzc6vqkiIjzI2JmRMwiG/z8WUS09bfYiPh34DFJB+anjgHubWGVijwKHCFpcv47cgxtPDDeR8ftaSJpHvBJ4KSRlL66lRwUGncU8AGyb9t35scJra7UCPXfgCsl3Q3MAf53i+szoLxFcw2wBriH7P+xtkvFIGkpcAtwoKSNks6gzfc0GaDO/wi8AliR/z/4Ty2t5AjgNBdmZtbLLQUzM+vloGBmZr0cFMzMrJeDgpmZ9XJQMDOzXg4KNmrk2VY/UvP86E7JFGs2XBwUbDSZCnyk8CqzUcxBwdqSpFl5nvzL8n0JrpR0rKRf5vn+D8/z//8gz6V/q6RD8/demOfev1nSQ5LOzou9CDggX+T0pfzclJq9Gq7MVyGbjVpjW10BszpeD5wKLAR+BfwF8IfAScCngceAOyLiFEnvBL5FtuIZ4I3AO8hWu94n6RtkifQOiYg5kHUfkWW3PZgsvfUvyVaq/2I4PpxZO3JLwdrZhoi4JyIqwDqyDWCCLH3ELLIA8W2AiPgZ8CpJu+fv/VFEbIuIJ8kSu+09wD1uj4iN+T3uzMs1G7UcFKydbat5XKl5XiFr5fbX1VPN21L73h4GbhWnXmc2KjgoWCdbCbwPeruCnizY0+I5su4kMxuAvxVZJ7uQbEe2u4EX2Zn2uV8R8VQ+UL2WbIeuH5VfRbPO4iypZmbWy91HZmbWy0HBzMx6OSiYmVkvBwUzM+vloGBmZr0cFMzMrJeDgpmZ9foP9EONHUhjtiEAAAAASUVORK5CYII=\n",
"text/plain": [
"<Figure size 432x288 with 2 Axes>"
]
},
"metadata": {
"needs_background": "light"
},
"output_type": "display_data"
}
],
"source": [
"ds_wo['GW'].groupby('time.month').mean('time').plot(x = 'month', vmin = vmin,\n",
" vmax = vmax, extend = 'both')"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.6.10"
}
},
"nbformat": 4,
"nbformat_minor": 4
}
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment