Skip to content
Snippets Groups Projects
MatfiletoNetCDF_CL_ARTofMELT_v1.ipynb 485 KiB
Newer Older
       ".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:                       (cloud_layer: 3, range_levels: 770, sky_condition_layer: 5, time: 2713)\n",
       "Coordinates:\n",
       "  * time                          (time) datetime64[ns] 2023-05-07T00:10:22.2...\n",
       "  * range_levels                  (range_levels) int32 1 2 3 4 ... 768 769 770\n",
       "  * cloud_layer                   (cloud_layer) int32 1 2 3\n",
       "  * sky_condition_layer           (sky_condition_layer) int32 1 2 3 4 5\n",
       "Data variables: (12/15)\n",
       "    day_of_year                   (time) float32 127.0 127.0 ... 164.7 164.7\n",
       "    year                          (time) int32 2023 2023 2023 ... 2023 2023 2023\n",
       "    month                         (time) int32 5 5 5 5 5 5 5 5 ... 6 6 6 6 6 6 6\n",
       "    day                           (time) int32 7 7 7 7 7 7 ... 13 13 13 13 13 13\n",
       "    hour                          (time) int32 0 0 0 1 1 1 ... 14 14 15 15 15 16\n",
       "    minute                        (time) int32 10 30 50 10 30 ... 50 10 30 50 10\n",
       "    ...                            ...\n",
       "    vertical_visibility           (time) float32 nan nan nan nan ... nan nan nan\n",
       "    highest_detected_signal       (time) float32 nan nan nan nan ... nan nan nan\n",
       "    cloud_base_altitude           (time, cloud_layer) float32 783.5 ... nan\n",
       "    sky_condition_cloud_fraction  (time, sky_condition_layer) float32 7.0 ......\n",
       "    sky_condition_cloud_altitude  (time, sky_condition_layer) float32 725.0 ....\n",
       "    backscatter_profile           (time, range_levels) float32 0.0005729 ... ...\n",
       "Attributes: (12/27)\n",
       "    Conventions:                           CF-1.8\n",
       "    source:                                Ceilometer\n",
       "    instrument_model:                      Vaisala Ceilometer CL31\n",
       "    creator_name:                          Sonja Murto\n",
       "    creator_email:                         sonja.murto@misu.su.se\n",
       "    creator_url:                           https://orcid.org/0000-0002-4966-9077\n",
       "    ...                                    ...\n",
       "    time_coverage_start:                   2023-05-07T00:10:22\n",
       "    time_coverage_end:                     2023-06-13T16:10:07\n",
       "    geospatial_bounds:                     80.52392166666667N, -3.87377499999...\n",
       "    platform_altitude:                     Located at approximately 25 m a.s.l\n",
       "    location_keywords:                     Oden, Arctic Ocean, Fram Strait, a...\n",
Sonja Murto's avatar
Sonja Murto committed
       "    comments:                              This file consists of 20 min avera...</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-cc4dbf4e-95da-4b51-b531-18898ed59713' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-cc4dbf4e-95da-4b51-b531-18898ed59713' 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'>cloud_layer</span>: 3</li><li><span class='xr-has-index'>range_levels</span>: 770</li><li><span class='xr-has-index'>sky_condition_layer</span>: 5</li><li><span class='xr-has-index'>time</span>: 2713</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-ff9271e5-d385-422b-8a6d-2cf645e8d253' class='xr-section-summary-in' type='checkbox'  checked><label for='section-ff9271e5-d385-422b-8a6d-2cf645e8d253' class='xr-section-summary' >Coordinates: <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 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'>2023-05-07T00:10:22.255854 ... 2...</div><input id='attrs-ca7af933-abd2-4b70-b5f7-ac1d78f068fc' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-ca7af933-abd2-4b70-b5f7-ac1d78f068fc' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1cb44d70-1265-41c7-8520-d322246eb4e5' class='xr-var-data-in' type='checkbox'><label for='data-1cb44d70-1265-41c7-8520-d322246eb4e5' 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;2023-05-07T00:10:22.255854000&#x27;, &#x27;2023-05-07T00:30:37.234550000&#x27;,\n",
       "       &#x27;2023-05-07T00:50:37.210550000&#x27;, ..., &#x27;2023-06-13T15:30:07.349600000&#x27;,\n",
       "       &#x27;2023-06-13T15:50:07.329201000&#x27;, &#x27;2023-06-13T16:10:07.312075000&#x27;],\n",
Sonja Murto's avatar
Sonja Murto committed
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>range_levels</span></div><div class='xr-var-dims'>(range_levels)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>1 2 3 4 5 6 ... 766 767 768 769 770</div><input id='attrs-bb65f784-f4a0-4b2a-84bb-f424d235c11d' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-bb65f784-f4a0-4b2a-84bb-f424d235c11d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-63e2f4f3-7a71-4f2b-a82e-6fb519f01844' class='xr-var-data-in' type='checkbox'><label for='data-63e2f4f3-7a71-4f2b-a82e-6fb519f01844' 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([  1,   2,   3, ..., 768, 769, 770], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>cloud_layer</span></div><div class='xr-var-dims'>(cloud_layer)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>1 2 3</div><input id='attrs-9a294552-fbbc-4671-9184-dad746b931dc' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-9a294552-fbbc-4671-9184-dad746b931dc' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7551b6b7-c2d4-4512-b620-a15abc2cba21' class='xr-var-data-in' type='checkbox'><label for='data-7551b6b7-c2d4-4512-b620-a15abc2cba21' 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([1, 2, 3], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>sky_condition_layer</span></div><div class='xr-var-dims'>(sky_condition_layer)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>1 2 3 4 5</div><input id='attrs-f292efa9-794d-49d0-944b-b99b190bc530' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-f292efa9-794d-49d0-944b-b99b190bc530' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-105db12f-76e9-4ec2-a708-5aa20607dd57' class='xr-var-data-in' type='checkbox'><label for='data-105db12f-76e9-4ec2-a708-5aa20607dd57' 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([1, 2, 3, 4, 5], dtype=int32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-ec238ad4-275a-4509-9d1e-73c30a3d32ed' class='xr-section-summary-in' type='checkbox'  ><label for='section-ec238ad4-275a-4509-9d1e-73c30a3d32ed' class='xr-section-summary' >Data variables: <span>(15)</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>day_of_year</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>127.0 127.0 127.0 ... 164.7 164.7</div><input id='attrs-9a1c3831-4179-47c5-a773-2ce2fc90b17c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9a1c3831-4179-47c5-a773-2ce2fc90b17c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b218596e-b849-4ba4-ae28-54aff9a9a2f2' class='xr-var-data-in' type='checkbox'><label for='data-b218596e-b849-4ba4-ae28-54aff9a9a2f2' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Day of Year</dd><dt><span>description :</span></dt><dd>time as decimal day of year</dd></dl></div><div class='xr-var-data'><pre>array([127.0072 , 127.02126, 127.03516, ..., 164.64592, 164.6598 ,\n",
       "       164.67369], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>year</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>2023 2023 2023 ... 2023 2023 2023</div><input id='attrs-22746881-afbc-444d-ad91-f9a85b373dce' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-22746881-afbc-444d-ad91-f9a85b373dce' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-23758d8c-63e0-438f-9fcb-56059b539843' class='xr-var-data-in' type='checkbox'><label for='data-23758d8c-63e0-438f-9fcb-56059b539843' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Year</dd></dl></div><div class='xr-var-data'><pre>array([2023, 2023, 2023, ..., 2023, 2023, 2023], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>month</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>5 5 5 5 5 5 5 5 ... 6 6 6 6 6 6 6 6</div><input id='attrs-81bcfd51-69a6-41e5-84ae-39d6549e2cf9' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-81bcfd51-69a6-41e5-84ae-39d6549e2cf9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5d24b69e-b054-492d-b7d6-f95dfd6b8f3c' class='xr-var-data-in' type='checkbox'><label for='data-5d24b69e-b054-492d-b7d6-f95dfd6b8f3c' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Month</dd></dl></div><div class='xr-var-data'><pre>array([5, 5, 5, ..., 6, 6, 6], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>day</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>7 7 7 7 7 7 7 ... 13 13 13 13 13 13</div><input id='attrs-8c916d9a-62a8-4965-8cd8-5ceb8a51e513' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-8c916d9a-62a8-4965-8cd8-5ceb8a51e513' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a74f4506-639b-4ee8-b200-0fb7804d61e0' class='xr-var-data-in' type='checkbox'><label for='data-a74f4506-639b-4ee8-b200-0fb7804d61e0' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Day</dd></dl></div><div class='xr-var-data'><pre>array([ 7,  7,  7, ..., 13, 13, 13], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>hour</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>0 0 0 1 1 1 2 ... 14 14 15 15 15 16</div><input id='attrs-8e752ec1-d9e7-4557-a74e-19536fdddf9a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-8e752ec1-d9e7-4557-a74e-19536fdddf9a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ef506819-1c9b-46d7-9a6a-c5a08666028c' class='xr-var-data-in' type='checkbox'><label for='data-ef506819-1c9b-46d7-9a6a-c5a08666028c' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Hour</dd></dl></div><div class='xr-var-data'><pre>array([ 0,  0,  0, ..., 15, 15, 16], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>minute</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>10 30 50 10 30 ... 50 10 30 50 10</div><input id='attrs-3e519d28-6fed-4ba8-8745-2183f24ae734' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3e519d28-6fed-4ba8-8745-2183f24ae734' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-10a487fd-a0e2-4495-aab1-5ecf13f402d4' class='xr-var-data-in' type='checkbox'><label for='data-10a487fd-a0e2-4495-aab1-5ecf13f402d4' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Minute</dd></dl></div><div class='xr-var-data'><pre>array([10, 30, 50, ..., 30, 50, 10], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>second</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>22 37 37 37 37 37 ... 7 7 7 7 7 7</div><input id='attrs-5b9130c2-ef3f-4a19-9823-06ab75a6d025' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-5b9130c2-ef3f-4a19-9823-06ab75a6d025' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-18ebf488-bc4c-4bc5-8eae-37584ae45c17' class='xr-var-data-in' type='checkbox'><label for='data-18ebf488-bc4c-4bc5-8eae-37584ae45c17' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Second</dd><dt><span>description :</span></dt><dd>Time averaged to closest second</dd></dl></div><div class='xr-var-data'><pre>array([22, 37, 37, ...,  7,  7,  7], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>ceilometer_range</span></div><div class='xr-var-dims'>(range_levels)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>35 45 55 65 ... 7695 7705 7715 7725</div><input id='attrs-3078d9d9-24fb-402f-8e45-dac81614bd97' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3078d9d9-24fb-402f-8e45-dac81614bd97' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-321231a1-a422-43ee-ba5a-cfa1ce4b0849' class='xr-var-data-in' type='checkbox'><label for='data-321231a1-a422-43ee-ba5a-cfa1ce4b0849' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>range_levels</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Ceilometer Range</dd><dt><span>description :</span></dt><dd>Ranges for the ceilometer backscatter profile, including the instrument height</dd></dl></div><div class='xr-var-data'><pre>array([  35,   45,   55,   65,   75,   85,   95,  105,  115,  125,  135,\n",
       "        145,  155,  165,  175,  185,  195,  205,  215,  225,  235,  245,\n",
       "        255,  265,  275,  285,  295,  305,  315,  325,  335,  345,  355,\n",
       "        365,  375,  385,  395,  405,  415,  425,  435,  445,  455,  465,\n",
       "        475,  485,  495,  505,  515,  525,  535,  545,  555,  565,  575,\n",
       "        585,  595,  605,  615,  625,  635,  645,  655,  665,  675,  685,\n",
       "        695,  705,  715,  725,  735,  745,  755,  765,  775,  785,  795,\n",
       "        805,  815,  825,  835,  845,  855,  865,  875,  885,  895,  905,\n",
       "        915,  925,  935,  945,  955,  965,  975,  985,  995, 1005, 1015,\n",
       "       1025, 1035, 1045, 1055, 1065, 1075, 1085, 1095, 1105, 1115, 1125,\n",
       "       1135, 1145, 1155, 1165, 1175, 1185, 1195, 1205, 1215, 1225, 1235,\n",
       "       1245, 1255, 1265, 1275, 1285, 1295, 1305, 1315, 1325, 1335, 1345,\n",
       "       1355, 1365, 1375, 1385, 1395, 1405, 1415, 1425, 1435, 1445, 1455,\n",
       "       1465, 1475, 1485, 1495, 1505, 1515, 1525, 1535, 1545, 1555, 1565,\n",
       "       1575, 1585, 1595, 1605, 1615, 1625, 1635, 1645, 1655, 1665, 1675,\n",
       "       1685, 1695, 1705, 1715, 1725, 1735, 1745, 1755, 1765, 1775, 1785,\n",
       "       1795, 1805, 1815, 1825, 1835, 1845, 1855, 1865, 1875, 1885, 1895,\n",
       "       1905, 1915, 1925, 1935, 1945, 1955, 1965, 1975, 1985, 1995, 2005,\n",
       "       2015, 2025, 2035, 2045, 2055, 2065, 2075, 2085, 2095, 2105, 2115,\n",
       "       2125, 2135, 2145, 2155, 2165, 2175, 2185, 2195, 2205, 2215, 2225,\n",
       "...\n",
       "       5645, 5655, 5665, 5675, 5685, 5695, 5705, 5715, 5725, 5735, 5745,\n",
       "       5755, 5765, 5775, 5785, 5795, 5805, 5815, 5825, 5835, 5845, 5855,\n",
       "       5865, 5875, 5885, 5895, 5905, 5915, 5925, 5935, 5945, 5955, 5965,\n",
       "       5975, 5985, 5995, 6005, 6015, 6025, 6035, 6045, 6055, 6065, 6075,\n",
       "       6085, 6095, 6105, 6115, 6125, 6135, 6145, 6155, 6165, 6175, 6185,\n",
       "       6195, 6205, 6215, 6225, 6235, 6245, 6255, 6265, 6275, 6285, 6295,\n",
       "       6305, 6315, 6325, 6335, 6345, 6355, 6365, 6375, 6385, 6395, 6405,\n",
       "       6415, 6425, 6435, 6445, 6455, 6465, 6475, 6485, 6495, 6505, 6515,\n",
       "       6525, 6535, 6545, 6555, 6565, 6575, 6585, 6595, 6605, 6615, 6625,\n",
       "       6635, 6645, 6655, 6665, 6675, 6685, 6695, 6705, 6715, 6725, 6735,\n",
       "       6745, 6755, 6765, 6775, 6785, 6795, 6805, 6815, 6825, 6835, 6845,\n",
       "       6855, 6865, 6875, 6885, 6895, 6905, 6915, 6925, 6935, 6945, 6955,\n",
       "       6965, 6975, 6985, 6995, 7005, 7015, 7025, 7035, 7045, 7055, 7065,\n",
       "       7075, 7085, 7095, 7105, 7115, 7125, 7135, 7145, 7155, 7165, 7175,\n",
       "       7185, 7195, 7205, 7215, 7225, 7235, 7245, 7255, 7265, 7275, 7285,\n",
       "       7295, 7305, 7315, 7325, 7335, 7345, 7355, 7365, 7375, 7385, 7395,\n",
       "       7405, 7415, 7425, 7435, 7445, 7455, 7465, 7475, 7485, 7495, 7505,\n",
       "       7515, 7525, 7535, 7545, 7555, 7565, 7575, 7585, 7595, 7605, 7615,\n",
       "       7625, 7635, 7645, 7655, 7665, 7675, 7685, 7695, 7705, 7715, 7725],\n",
Sonja Murto's avatar
Sonja Murto committed
       "      dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>flag_cloudcode</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>int32</div><div class='xr-var-preview xr-preview'>1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1 1</div><input id='attrs-6a2fb53a-b3e1-477d-9157-7f7cedd693f6' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6a2fb53a-b3e1-477d-9157-7f7cedd693f6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f1b13da4-9e2f-4789-a929-8afbbb87bc5a' class='xr-var-data-in' type='checkbox'><label for='data-f1b13da4-9e2f-4789-a929-8afbbb87bc5a' 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>type :</span></dt><dd>int32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>1</dd><dt><span>long_name :</span></dt><dd>Data flag: Cloudcode</dd><dt><span>flag_values :</span></dt><dd>-1,0,1,2,3,4</dd><dt><span>flag_meanings :</span></dt><dd>missing_data\n",
       "no_significant_backscatter\n",
       "one_cloud_base_detected\n",
       "two_cloud_bases_detected\n",
       "three_cloud_bases_detected\n",
Sonja Murto's avatar
Sonja Murto committed
       "full_obscuration</dd><dt><span>description :</span></dt><dd>Code for number of cloud bases detected; see Readme document for more information</dd></dl></div><div class='xr-var-data'><pre>array([1, 1, 1, ..., 1, 1, 1], dtype=int32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>vertical_visibility</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>nan nan nan nan ... nan nan nan nan</div><input id='attrs-47a93c8e-1c1d-47c8-be30-6d490780e6f1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-47a93c8e-1c1d-47c8-be30-6d490780e6f1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5489a621-0013-43e8-a8b5-c293299cd1ab' class='xr-var-data-in' type='checkbox'><label for='data-5489a621-0013-43e8-a8b5-c293299cd1ab' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Vertical Visibilility</dd><dt><span>description :</span></dt><dd>Vertical visibility given in case of obscured cloud base (at flag_cloudcode 4), else NaN</dd></dl></div><div class='xr-var-data'><pre>array([nan, nan, nan, ..., nan, nan, nan], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>highest_detected_signal</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>nan nan nan nan ... nan nan nan nan</div><input id='attrs-37918ffb-2c4b-406d-a1a8-095eb3773454' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-37918ffb-2c4b-406d-a1a8-095eb3773454' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-bd472c4e-3a0b-47e1-8bc1-d7ce1ed64967' class='xr-var-data-in' type='checkbox'><label for='data-bd472c4e-3a0b-47e1-8bc1-d7ce1ed64967' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Highest Signal Detected</dd><dt><span>description :</span></dt><dd>Highest signal detected given in case of obscured cloud base (at flag_cloudcode 4), else NaN</dd></dl></div><div class='xr-var-data'><pre>array([nan, nan, nan, ..., nan, nan, nan], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>cloud_base_altitude</span></div><div class='xr-var-dims'>(time, cloud_layer)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>783.5 938.3 nan ... 411.8 nan nan</div><input id='attrs-72d149b8-8aa0-4766-9e43-6321d263ef73' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-72d149b8-8aa0-4766-9e43-6321d263ef73' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-adda3c4f-5163-4f74-8b79-4dc7426a5a42' class='xr-var-data-in' type='checkbox'><label for='data-adda3c4f-5163-4f74-8b79-4dc7426a5a42' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time, cloud_layer</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Cloud Base Altitude</dd><dt><span>description :</span></dt><dd>cloud base height of 1-3 cloud layers; NaN if no layer detected. Instrument height incorporated.</dd></dl></div><div class='xr-var-data'><pre>array([[ 783.53656,  938.3333 ,        nan],\n",
       "       [ 709.75   ,  971.9231 , 1180.     ],\n",
       "       [ 677.25   ,  851.6667 , 1135.     ],\n",
       "       ...,\n",
       "       [ 255.25   ,  470.45456,        nan],\n",
       "       [ 596.5    ,  755.     ,        nan],\n",
Sonja Murto's avatar
Sonja Murto committed
       "       [ 411.75   ,        nan,        nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sky_condition_cloud_fraction</span></div><div class='xr-var-dims'>(time, sky_condition_layer)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>7.0 8.0 0.0 0.0 ... 0.0 0.0 0.0 0.0</div><input id='attrs-721b08f2-de6c-41bd-b5eb-571d309bf65a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-721b08f2-de6c-41bd-b5eb-571d309bf65a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-581def9f-3516-4d1b-8284-f0e4f70d8785' class='xr-var-data-in' type='checkbox'><label for='data-581def9f-3516-4d1b-8284-f0e4f70d8785' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time, sky_condition_layer</dd><dt><span>units :</span></dt><dd>octal</dd><dt><span>long_name :</span></dt><dd>Sky Condition Cloud Fraction</dd><dt><span>description :</span></dt><dd>Cloud fraction calculated with the sky condition algorithm. 0-8 = cloud coverage of up to 5 levels; 9 = obscuration. NaN = missing data or no detected layer.</dd></dl></div><div class='xr-var-data'><pre>array([[7., 8., 0., 0., 0.],\n",
       "       [8., 0., 0., 0., 0.],\n",
       "       [8., 0., 0., 0., 0.],\n",
       "       ...,\n",
       "       [8., 0., 0., 0., 0.],\n",
       "       [8., 0., 0., 0., 0.],\n",
Sonja Murto's avatar
Sonja Murto committed
       "       [8., 0., 0., 0., 0.]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>sky_condition_cloud_altitude</span></div><div class='xr-var-dims'>(time, sky_condition_layer)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>725.0 1.525e+03 nan ... nan nan nan</div><input id='attrs-58593bd7-06c4-4a1b-abfd-5c8bde7ee87c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-58593bd7-06c4-4a1b-abfd-5c8bde7ee87c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-58f9031a-a32a-43e5-a9b5-cb4d2d3326a4' class='xr-var-data-in' type='checkbox'><label for='data-58f9031a-a32a-43e5-a9b5-cb4d2d3326a4' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time, sky_condition_layer</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>long_name :</span></dt><dd>Sky Condition Cloud Altitude</dd><dt><span>description :</span></dt><dd>Cloud layer height calculated with the sky condition algorithm. Cloud layer height given for 1-5 sky condition layers; NaN if no layer detected. Vertical visibility is reported as height if obscuration (at sky_condition_cloud_fraction 9). Instrument height incorporated.</dd></dl></div><div class='xr-var-data'><pre>array([[ 725., 1525.,   nan,   nan,   nan],\n",
       "       [ 655.,   nan,   nan,   nan,   nan],\n",
       "       [ 585.,   nan,   nan,   nan,   nan],\n",
       "       ...,\n",
       "       [ 155.,   nan,   nan,   nan,   nan],\n",
       "       [ 235.,   nan,   nan,   nan,   nan],\n",
Sonja Murto's avatar
Sonja Murto committed
       "       [ 365.,   nan,   nan,   nan,   nan]], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>backscatter_profile</span></div><div class='xr-var-dims'>(time, range_levels)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>0.0005729 0.0004571 ... -0.0003775</div><input id='attrs-96012a65-865f-494c-9765-81b9a1a5c4ca' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-96012a65-865f-494c-9765-81b9a1a5c4ca' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6dc18bb6-25fc-4353-b06d-864dd60e4c95' class='xr-var-data-in' type='checkbox'><label for='data-6dc18bb6-25fc-4353-b06d-864dd60e4c95' 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>type :</span></dt><dd>float32</dd><dt><span>dimension :</span></dt><dd>time, range_levels</dd><dt><span>units :</span></dt><dd>1 km-1 steradians-1</dd><dt><span>long_name :</span></dt><dd>Backscatter Profile</dd><dt><span>description :</span></dt><dd>backscatter coefficient profile</dd></dl></div><div class='xr-var-data'><pre>array([[ 5.7292683e-04,  4.5707318e-04,  4.5439025e-04, ...,\n",
       "         9.9268291e-05,  3.9487804e-04,  4.8731707e-04],\n",
       "       [ 1.5730000e-03,  1.4250000e-03,  1.4587500e-03, ...,\n",
       "        -4.4500001e-05, -7.3224999e-04, -5.9900002e-04],\n",
       "       [ 1.4905750e-02,  8.1692496e-03,  9.2337504e-03, ...,\n",
       "        -1.0550000e-04, -2.7724999e-04,  6.9250003e-05],\n",
       "       ...,\n",
       "       [ 2.8947501e-03,  3.6892500e-03,  4.4765002e-03, ...,\n",
       "         2.6700000e-04,  7.8499998e-04,  3.2774999e-04],\n",
       "       [ 3.1457499e-03,  2.5472499e-03,  3.0722499e-03, ...,\n",
       "        -6.9825002e-04, -1.0052500e-03, -8.0824998e-04],\n",
       "       [ 1.5090750e-02,  9.6795000e-03,  1.4006000e-02, ...,\n",
Sonja Murto's avatar
Sonja Murto committed
       "        -2.8924999e-04,  1.2325001e-04, -3.7749999e-04]], dtype=float32)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-86683f8c-5302-4e41-bdf2-991b861c6427' class='xr-section-summary-in' type='checkbox'  ><label for='section-86683f8c-5302-4e41-bdf2-991b861c6427' class='xr-section-summary' >Attributes: <span>(27)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>Conventions :</span></dt><dd>CF-1.8</dd><dt><span>source :</span></dt><dd>Ceilometer</dd><dt><span>instrument_model :</span></dt><dd>Vaisala Ceilometer CL31</dd><dt><span>creator_name :</span></dt><dd>Sonja Murto</dd><dt><span>creator_email :</span></dt><dd>sonja.murto@misu.su.se</dd><dt><span>creator_url :</span></dt><dd>https://orcid.org/0000-0002-4966-9077</dd><dt><span>institution :</span></dt><dd>Stockholm University</dd><dt><span>processing_software :</span></dt><dd>Matlab (for creating the matlab file) and a jupyter notebook script (netCDF)</dd><dt><span>sampling_interval :</span></dt><dd>original 30s; here 20min averages</dd><dt><span>product_version :</span></dt><dd>v01</dd><dt><span>last_revised_date :</span></dt><dd>2024-05-31T15:00:00</dd><dt><span>project :</span></dt><dd>ARTofMELT</dd><dt><span>project_principal_investigator :</span></dt><dd>Michael Tjernström</dd><dt><span>project_principal_investigator_email :</span></dt><dd>michaelt@misu.su.se</dd><dt><span>project_principal_investigator_url :</span></dt><dd>https://orcid.org/0000-0002-6908-7410</dd><dt><span>acknowledgement :</span></dt><dd> Knut och Alice Wallenbergs Stiftelse, Grant 2016-0024</dd><dt><span>platform :</span></dt><dd>Swedish Icebreaker Oden</dd><dt><span>platform_type :</span></dt><dd>On Oden&#x27;s 7th deck above the bridge</dd><dt><span>deployment_mode :</span></dt><dd>ship</dd><dt><span>title :</span></dt><dd>Ceilometer cloud base height, vertical visibility and backscatter profiles</dd><dt><span>feature_type :</span></dt><dd>time series</dd><dt><span>time_coverage_start :</span></dt><dd>2023-05-07T00:10:22</dd><dt><span>time_coverage_end :</span></dt><dd>2023-06-13T16:10:07</dd><dt><span>geospatial_bounds :</span></dt><dd>80.52392166666667N, -3.8737749999999997E, 78.04355166666666N, 15.660881666666667E</dd><dt><span>platform_altitude :</span></dt><dd>Located at approximately 25 m a.s.l</dd><dt><span>location_keywords :</span></dt><dd>Oden, Arctic Ocean, Fram Strait, atmosphere, on the ship</dd><dt><span>comments :</span></dt><dd>This file consists of 20 min averages of ceilometer data measured with the Vaisala Ceilometer CL31 that was located on the 7th deck, above the bridge (at approximately 25m).The sky condition measurements are time averages to represent an area average. The vertical resolution is 10m * 770, but the measurement height (25m) is included in the backscatter profile ranges, as well as in the cloud base heights (cloud_base_altitude and sky_condition_cloud_altitude). Geospatial bounds are taken from the gps location of the weather station dataset located on Oden. Time variables month, day, hour, minute and second are approximated to the nearest second. Data produced by Sonja Murto. See the document - Readme_CL.txt - for more details.</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:                       (cloud_layer: 3, range_levels: 770, sky_condition_layer: 5, time: 2713)\n",
       "Coordinates:\n",
       "  * time                          (time) datetime64[ns] 2023-05-07T00:10:22.2...\n",
       "  * range_levels                  (range_levels) int32 1 2 3 4 ... 768 769 770\n",
       "  * cloud_layer                   (cloud_layer) int32 1 2 3\n",
       "  * sky_condition_layer           (sky_condition_layer) int32 1 2 3 4 5\n",
       "Data variables: (12/15)\n",
       "    day_of_year                   (time) float32 127.0 127.0 ... 164.7 164.7\n",
       "    year                          (time) int32 2023 2023 2023 ... 2023 2023 2023\n",
       "    month                         (time) int32 5 5 5 5 5 5 5 5 ... 6 6 6 6 6 6 6\n",
       "    day                           (time) int32 7 7 7 7 7 7 ... 13 13 13 13 13 13\n",
       "    hour                          (time) int32 0 0 0 1 1 1 ... 14 14 15 15 15 16\n",
       "    minute                        (time) int32 10 30 50 10 30 ... 50 10 30 50 10\n",
       "    ...                            ...\n",
       "    vertical_visibility           (time) float32 nan nan nan nan ... nan nan nan\n",
       "    highest_detected_signal       (time) float32 nan nan nan nan ... nan nan nan\n",
       "    cloud_base_altitude           (time, cloud_layer) float32 783.5 ... nan\n",
       "    sky_condition_cloud_fraction  (time, sky_condition_layer) float32 7.0 ......\n",
       "    sky_condition_cloud_altitude  (time, sky_condition_layer) float32 725.0 ....\n",
       "    backscatter_profile           (time, range_levels) float32 0.0005729 ... ...\n",
       "Attributes: (12/27)\n",
       "    Conventions:                           CF-1.8\n",
       "    source:                                Ceilometer\n",
       "    instrument_model:                      Vaisala Ceilometer CL31\n",
       "    creator_name:                          Sonja Murto\n",
       "    creator_email:                         sonja.murto@misu.su.se\n",
       "    creator_url:                           https://orcid.org/0000-0002-4966-9077\n",
       "    ...                                    ...\n",
       "    time_coverage_start:                   2023-05-07T00:10:22\n",
       "    time_coverage_end:                     2023-06-13T16:10:07\n",
       "    geospatial_bounds:                     80.52392166666667N, -3.87377499999...\n",
       "    platform_altitude:                     Located at approximately 25 m a.s.l\n",
       "    location_keywords:                     Oden, Arctic Ocean, Fram Strait, a...\n",
       "    comments:                              This file consists of 20 min avera..."
      ]
     },
Sonja Murto's avatar
Sonja Murto committed
     "execution_count": 77,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds_all"
   ]
  },
  {
   "cell_type": "code",
Sonja Murto's avatar
Sonja Murto committed
   "execution_count": 78,
   "metadata": {},
   "outputs": [],
   "source": [
    "#save to netCDF\n",
    "ds_all.to_netcdf(load_data + 'CL31_ceilometer_ARTofMELT_20230507_20230613_20min_v01.nc')"
   ]
  },
  {
   "cell_type": "code",
Sonja Murto's avatar
Sonja Murto committed
   "execution_count": 79,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYsAAAElCAYAAAAV9s4VAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8vihELAAAACXBIWXMAAAsTAAALEwEAmpwYAABroklEQVR4nO2ddbwc1fXAv2f3SdwTiBASIMETJIRQJGjxAkUKLdoCRYqUXylWKKVYkQLBnQDFKcUdgoaEBCLEE4gRd3m+e35/3JndWZe39jb3+/nse7N37N7ZmTn3yD1XVBWLxWKxWJLhK3YFLBaLxVL6WGFhsVgslpRYYWGxWCyWlFhhYbFYLJaUWGFhsVgslpRYYWGxWCyWlORVWIjIn0Vkioj8ICLPi0grEekiIh+KyCznf2fP9leJyGwRmSEih3rKdxeRyc66ESIi+ay3xWKxWCLJm7AQkd7AxcAQVd0J8AMnA1cCH6vqAOBj5zsisoOzfkfgMOABEfE7h3sQOBcY4HwOy1e9LRaLxRJLvs1QFUBrEakA2gCLgGOAkc76kcCxzvIxwAuqWq+qPwGzgaEi0hPooKqj1YwgfNqzj8VisVgKQEW+DqyqP4vIHcB8oBb4QFU/EJHNVHWxs81iEenh7NIb+MZziIVOWaOzHF0eg4ici9FAaNu27e7bbbddLptksVgsZc/48eNXqGr36PK8CQvHF3EM0B9YA7wsIqcm2yVOmSYpjy1UfQR4BGDIkCE6bty4TKpssVgsmzwiMi9eeT7NUAcDP6nqclVtBP4L/AJY6piWcP4vc7ZfCGzh2b8Pxmy10FmOLrdYLBZLgcinsJgPDBORNk700kHANOAN4AxnmzOA153lN4CTRaRaRPpjHNljHZPVehEZ5hzndM8+FovFYikA+fRZjBGRV4DvgCbge4yJqB3wkoj8ASNQTnS2nyIiLwFTne0vVNWAc7jzgaeA1sC7zsdisVgsBULKNUW59VlYLPmhsbGRhQsXUldXV+yqWJpBq1at6NOnD5WVlRHlIjJeVYdEb583zcJisZQnCxcupH379vTr1w87PrZloqqsXLmShQsX0r9//7T2sek+LBZLRtTV1dG1a1crKFowIkLXrl0z0g6tsLBYLBljBUXLJ9Pf0AoLi8VisaTECguLxVIWXH/99dxxxx05Odb+++9PsgCZfv36sWLFipycKx8cdthhdOrUiaOOOipnx7TCwmKxWEqMQCCQeqMkXH755TzzzDM5qo3BRkNZLJbsufRSmDAht8fcZRe4++6kmzz99NPccccdiAiDBg2KeTFOmDCB8847j5qaGrbeemueeOIJOnfuzP77788dd9zBkCFDWLFiBUOGDGHu3LnU1tZy1llnMXXqVLbffntqa2vTru6xxx7LggULqKur45JLLuHcc8/l8ccf54cffuCuu+4C4NFHH2XatGn8+9//5tlnn2XEiBE0NDSw55578sADD+D3+2nXrh2XXXYZ77//PnfeeSf77LNPplcuxEEHHcSoUaOy3j8eVrOwWCwtiilTpnDTTTfxySefMHHiRO65556YbU4//XT+9a9/MWnSJHbeeWf+8Y9/JD3mgw8+SJs2bZg0aRLXXHMN48ePT7s+TzzxBOPHj2fcuHGMGDGClStXcvLJJ/PGG2/Q2NgIwJNPPslZZ53FtGnTePHFF/nqq6+YMGECfr+f//znPwBs3LiRnXbaiTFjxsQIittvv51ddtkl5nPxxRenXc/mYjULi8WSPSk0gHzwySefcMIJJ9CtWzcAunTpErF+7dq1rFmzhuHDhwNwxhlncOKJJyY95ueffx568Q4aNIhBgwalXZ8RI0bw2muvAbBgwQJmzZrFsGHDOPDAA3nrrbfYfvvtaWxsZOedd+a+++5j/Pjx7LHHHgDU1tbSo4dJvO33+zn++OPjnuPyyy/n8ssvT7tO+cAKC4vF0qJQ1axDdysqKggGgwAxYwyyOeaoUaP46KOPGD16NG3atGH//fcPHffss8/m5ptvZrvttuOss84K1f2MM87glltuiTlWq1at8Pv9MeVgNAtXA/Gy3377MWLEiIzrnQ3WDGWxWFoUBx10EC+99BIrV64EYNWqVRHrO3bsSOfOnfniiy8AeOaZZ0JaRr9+/UImpldeeSW0z3777Rd6Gf/www9MmjQprbqsXbuWzp0706ZNG6ZPn84334Sn5Nlzzz1ZsGABzz33HKecckqo7q+88grLli0L1X3evLgZwSO4/PLLmTBhQsynUIICrGZhsVhaGDvuuCPXXHMNw4cPx+/3s+uuu/LUU09FbDNy5MiQg3urrbbiySefBOAvf/kLJ510Es888wwHHnhgaPvzzz+fs846i0GDBrHLLrswdOjQtOpy2GGH8dBDDzFo0CC23XZbhg0bFrH+pJNOYsKECXTu3BmAHXbYgRtvvJFf/vKXBINBKisruf/++9lyyy2bcUVi2XfffZk+fTobNmygT58+PP744xx66KHNOqZNJGixWDJi2rRpbL/99sWuRovgqKOO4s9//jMHHXRQsasSl3i/ZaJEgtYMZbFYLDlmzZo1DBw4kNatW5esoMgUa4ayWCyWBOy5557U19dHlD3zzDPsvPPOSffr1KkTM2fOzGfVCo4VFhaLxZKAMWPGFLsKJYM1Q1ksFoslJXkTFiKyrYhM8HzWicilItJFRD4UkVnO/86efa4SkdkiMkNEDvWU7y4ik511I8TmR7ZYLJaCkjdhoaozVHUXVd0F2B2oAV4DrgQ+VtUBwMfOd0RkB+BkYEfgMOABEXFHqDwInAsMcD6H5aveFovFYomlUGaog4A5qjoPOAYY6ZSPBI51lo8BXlDVelX9CZgNDBWRnkAHVR2tJs73ac8+FovFYikAhRIWJwPPO8ubqepiAOd/D6e8N7DAs89Cp6y3sxxdHoOInCsi40Rk3PLly3NYfYvFUurY+SwMEyZMYK+99mLHHXdk0KBBvPjiizk5bt6FhYhUAb8CXk61aZwyTVIeW6j6iKoOUdUh3bt3z6yiFovFUiI0Zz6LNm3a8PTTTzNlyhTee+89Lr30UtasWdPsOhUidPZw4DtVXep8XyoiPVV1sWNiWuaULwS28OzXB1jklPeJU26xWIpMkaazsPNZJGHgwIGh5V69etGjRw+WL19Op06dsjqeSyHMUKcQNkEBvAGc4SyfAbzuKT9ZRKpFpD/GkT3WMVWtF5FhThTU6Z59LBbLJoadzyL9+SzGjh1LQ0MDW2+9ddrtSUReNQsRaQMcAvzRU3wr8JKI/AGYD5wIoKpTROQlYCrQBFyoqq4udj7wFNAaeNf5WCyWIlOE6SzsfBZpsnjxYk477TRGjhyJz9d8vSCvwkJVa4CuUWUrMdFR8ba/CbgpTvk4YKd81NFisbQs7HwWYRLNZ7Fu3TqOPPJIbrzxxphMuNliR3BbLJYWhZ3PIvl8Fg0NDRx33HGcfvrpKTWqTLC5oSwWS4vCzmeRnJdeeonPP/+clStXhq7LU089xS677NKs49r5LCwWS0bY+SzSx85nYbFYLJaE2PksLBaLZRPCzmcRxgoLi8ViSYCdzyKMNUNZLOVIXR1ceCF8/HGxa2IpE6ywsFjKkcZGeOCB3OfisGyyWGFhsVgslpRYYWGxWCyWlFhhYbFYygI7n4Vh3rx57L777uyyyy7suOOOPPTQQzk5ro2GslgslhIjEAgkzBOVip49e/L1119TXV3Nhg0b2GmnnfjVr35Fr169mlUnKywslnImzxkaZs26lA0bJuT0mO3a7cKAAXcn3cbOZ5GYqqqq0HJ9fX0ocWJzsWYoi6UcyTIra0vAzmeRej6LBQsWMGjQILbYYguuuOKKZmsVYDULi8XSDFJpAPnAzmeRmi222IJJkyaxaNEijj32WE444QQ222yztPePhxUWFoulRWHnswiTaD4Ll169erHjjjvyxRdfcMIJJ2TcPi/WDGWxWFoUdj6L5PNZLFy4MORzWb16NV999RXbbrttWu1JRl6FhYh0EpFXRGS6iEwTkb1EpIuIfCgis5z/nT3bXyUis0Vkhogc6infXUQmO+tGSLbdCotlU6MMpyDwzmcxePBgLrvssphtRo4cyeWXX86gQYOYMGEC1113HWDms3jwwQf5xS9+ERH6ev7557NhwwYGDRrEbbfdltF8Fk1NTQwaNIhrr7027nwWe++9d9z5LAYNGsQhhxzC4sWLs70UcZk2bRp77rkngwcPZvjw4fzlL39JmfgwLVQ1bx9gJHC2s1wFdAJuA650yq4E/uUs7wBMBKqB/sAcwO+sGwvsBQhm/u3DU5179913V4tlk2XDBlVQve22nB966tSpOT9muXLkkUfqRx99VOxqJCTebwmM0zjv1KQ+CxHZLQ1506iqk+Ps2wHYDzjTEUoNQIOIHAPs72w2EhgFXAEcA7ygqvXATyIyGxgqInOBDqo62jnu08CxjtCwWFomL74Iu+8O22xT7JpY8sCaNWsYOnQogwcP3mTms/gM+BbTo09Ef6BfnPKtgOXAkyIyGBgPXAJspqqLAVR1sYj0cLbvDXzj2X+hU9boLEeXxyAi5wLnAvTt2zdZuyyW4nLyydCqFWQQz28pPHY+izCphMW3qnpgsg1E5JMkx94NuEhVx4jIPRizU8JDxSnTJOWxhaqPAI+AmVY1ybkslqJzV9151N4MV19d7JpkjjYjIqklUc7zWWiG/qykDu5UgiLFNguBharqXu1XMMJjqYj0BHD+L/Nsv4Vn/z7AIqe8T5xySyEZPx6amopdi7LiMu7immuKXYvMadWqFStXrsz4ZWMpHVSVlStX0qpVq7T3SXuchYgMwpibQvuo6n+TVGaJiCwQkW1VdQZwEDDV+ZwB3Or8f93Z5Q3gORH5N9ALGACMVdWAiKwXkWHAGOB04N60W2hpPpMmwZAhcOWVECc+3FLC5OGF3qdPHxYuXMjy5ctzfmxL4WjVqhV9+vRJvaFDWsJCRJ4ABgFTADfRiAIJhYXDRcB/RKQK+BE4C6PNvCQifwDmAycCqOoUEXkJI0yagAtVNeAc53zgKaA1xrFtnduFxA3t++674tbDkj55NBFVVlbSv3//vB3fUpqkq1kMU9UdMj24qk4AhsRZFTc8QFVvAm6KUz4O2CnT81tyhJuq2RkEZbFYNj3SHZQ3WkQyFhaWMmHqVPN/7dri1sNisRSNdDWLkRiBsQSox0Qoqaqmn23L0nJxTRqzZxe3HhaLpWikKyyeAE4DJhP2WVg2FTaBEMlC09AZ+nWcwty5Oxa7KhZLWqQrLOar6ht5rYmldLHCIucsPA4e++1gDj44z+HINrzVkiPSFRbTReQ54E2MGQpIHjprsViSkG/5awW8JcekKyxaY4TELz1l6YTOWsoBn81kn3MEVO0L3dJySEtYqOpZ+a6IpYSxvdS8YIWFpSWRtMvoJOZLSjrbWFo4VlhYLJs8qTSLK0VkRZL1gskk+0juqmQpOaKFxZ13wiuvwOjRxalPGaCFMkNZB7clR6STovzoFNt8mKO6WEqVaGGxZAlMnFicupQLEvqTp+NbbdCSW5IKC+ursACxDm6/HwKB+NtaUuP09q3PwtKSsGEultRE91KtsGgeqvkPnbVYcowVFpbMcYWFtYc3C6tZWFoSVlhYUhOtWTz0kPn/zTex21pSo1owB3dDk481a/J+GssmQFrCQkQ2E5HHReRd5/sOznwUlk0Br7BQBXfSm1mzilOfssFc12Aes62Nmd+Tzp3hQxuGYmkm6WoWTwHvY2awA5gJXJqH+lhKEVeTgEjTkzVDZUeUzyIv7h9HwM9Y0RWArbfOwzksmxTpCotuqvoSTsZZVW0CrIdzUySf3eBNCY8ZKp9Tm/9m0DS++Qa23DJ/57BsGqQrLDaKSFdMPiic+bBTzoQjInNFZLKITBCRcU5ZFxH5UERmOf87e7a/SkRmi8gMETnUU767c5zZIjJCxAaRF5p6qlhLBzRghUWziQqdzaewaF/dwJ57mpgEi6U5pCssLgPeALYWka+ApzHza6fDAaq6i6q606teCXysqgOAj53vODPxnQzsCBwGPCAi7i3+IHAuMMD5HJbmuS054i7+TCfWUler8M9/msLevYtbqZZKlPnORiFbWgJpCQtV/Q4YDvwC+COwo6pOyvKcx2Bm3sP5f6yn/AVVrVfVn4DZwFAR6Ql0UNXRqqoYQXUsloLic+a80kAQhg83hVbByxotkBmKV14xv1NVVR5PYtkUSDqCW0R+nWDVQBFJZz4LBT4QEQUeVtVHgM1UdTGAqi4WkR7Otr0BbyzmQqes0VmOLo9X33MxGgh9+/ZNUTVLJoixQBJsCkKFc9vk9S1XxoQc3AUQFt99Z/43NubxJJZNgVS5ody8UD0wWsUnzvcDgFGkns9ib1Vd5AiED0VkepJt43VTNUl5bKERRo8ADBkyxIbq5IqTTsL3kqNZBDVsALfColnkVbOwWp8lx6SVG0pE3gJ2cDUCxzR0f6qDq+oi5/8yEXkNGAosFZGejlbRE1jmbL4Q2MKzex9gkVPeJ065pVBUVITMUMGmIOy6qxlr0aFDkSvWMqnZqIyWYXRmOgfyMb16HZT/k261Vf7PYSlr0nVw93MFhcNSYGCyHUSkrYi0d5cxs+z9gHGUn+FsdgbwurP8BnCyiFSLSH+MI3usc971IjLMiYI63bOPpUCEzFBBhcpK6NbN2sGzpKYGxjEEVaEdG/B9+nH+T1pZmf9z5Iv77oNPPkm9nSWvpDut6igReR94HmMCOhn4NMU+mwGvOVGuFcBzqvqeiHwLvOSMAJ8PnAigqlNE5CVgKtAEXKiqbpzI+ZiBga2Bd52PpVCIeBzc1rrXXIKBcOisoPDWW3BQnrWLliwsLnICL+0g0KKS7rSqf3Kc3fs6RY+o6msp9vkRGBynfCUQ98lQ1ZuAm+KUjwN2SqeulvwQ4eC2NItgEEzMh3NdGxryf1KrBVqaSbqahRv5lMqhbSlHvJpF0PbumktYWDiaRX19/k/akjULS0mQlrAQkfWEI5CqgEpgo6paD+emgEdYWM2i+QQDioiGs87mQ1hER0NZYWFpJukOymuvqh2cTyvgeOC+/FbNUjKIhM1Q1mfRbNz0WiGfRSE0iyKaoV59OcBJJ2lWIcKq8Dn7xo+VtxSUrOazUNX/AQfmtiqWksXns2aoHGKuocdnUeZmqKknXc/LL0tWQz9eew1O3eIRHu1wRuqNLXklXTOUdyS3DxhCgoFxljxy/PEmhHD16sKet7KypMxQTz1lsqgecECxa5Idrs+ioJpFmzb5P0cCgk6fNHoq93T48Ufl6ae3Z8bIQ3JcK0umpOvgPtqz3ATMxeRyshSS/xYvviA8zqJoVQhx9dVw5JEtWFiETHliXqSFEBY335z/cyQggBnxn41mIWJsV9pkJ/UsNukKi8dU9StvgYjsTXj0taUQdOoEa9bAL38JH3xQ0FNHJBIsMq1aQV1dsWuRPd7Q2f9yPNy9Te5P4n0z9+wJ222X+3OkSRAffprIIPgyhIjJaaVNNsd6sUlXXN+bZpkln7gDt6KTwn36qbHN5AvVknJwl4uwCEVD9eiRfIeWTDBIkLDPK3PMGJSgFRZFJ1XW2b0wCQS7i8hlnlUdAPvrFRp3BGu0sHjuOXjnHTjzzLydupQc3K1aFcZyky9cB3dIWOR7ZqJiJhVsamqWsAhrFtYMVWxS/QJVQDuMUGnv+awDTshv1SwJKUK66VJycJeDZgEmdPY4/lvemkVjIwH85v7JZpYnNb0Cq1kUn1RZZz8DPhORp1R1XoHqZElEIs2iADlzrBkqd7iD8gAe4jwg0bQxZYCjWfgJmB+tbduMdtfJM/lw5e9ovyCz/Sy5J6lmISJ3O4v3icgb0Z/8V88SlyhhEfQFWDaslo0bp+bnfKocydv8MOi3bNm/+OaA6uoWLiw86T6Kwfr141m9OlUe0BzR2Bg2Q2WRA6tpeQU33/wsqyZunYfKWTIhVXjCM87/O/JdEUsauBpE1FBYFWXqRavpv+IN2rbdIS+n7tS7HZ0mPpeXY2dKi9csao1pRVXowfL8nMTrp4jyWcyffxsbNkxkzz2TzUWWI7w+iywcTVokgWqJJZUZaryzuIuq3uNdJyKXAJ/lq2KWOCQwQ81d3ZOGulbMnLmCLbcsQr0KzJ/+BOvXF7sW2RN87wOkorjmPCmU09vrs2iGhBc7BrjopGtTiDfW/swc1sOSAQtru7Lc2yFVpba+DXUtubudAQcdBMceW+xaZI/6fEREQxW+BoU7lddnkcWITqtZlA6pQmdPAX4L9I/yUbQHVuazYpY4OJrFvoteYN9Lm3j6P+bna+VvJBj0E54rKj/nteSGQEUdhx/6VBFrUEB/iddnYe+jFk0qn8XXwGKgG3Cnp3w9MClflbIkx0cQra3H/flcYSGSJ2EBxY3VLzOC/izSrzaHqN9OtYDCwuuzaIawsGao4pPUDKWq81R1lKrupaqfeT7fqWpad7yI+EXkexF5y/neRUQ+FJFZzv/Onm2vEpHZIjJDRA71lO8uIpOddSOkYAbXEkPDmUqDEo47b+VrIBj0E8gmjt1ScII+jznmyiuLUIMCCovttw/7LLIQFtYMVTqkCp1dLyLr4nzWi8i6NM9xCTDN8/1K4GNVHQB87HxHRHbAzO29I3AY8IBI6I34IHAuMMD5HJbmucsSH8GIfla1v4lAwE8waIVFSyDo9/xOnTrl5yRJ+1NaOAe3SLN8Fi4v8pscVsqSDak0C++kR95P+3RmyRORPsCRwGOe4mOAkc7ySOBYT/kLqlqvqj8Bs4GhItIT6KCqo9Xoz0979tm08GoWnvdNhS+IWs2ixRBs75lbIt+pPhJSuB57c8xQrmYxl/4EAjUEgwWYr9wSl4xHWIlIWxH5nYi8ncbmdwN/hYjEMJup6mIA57+b66A3sMCz3UKnrLezHF0er27nisg4ERm3fHme4tdLAB9Bx+4cRoO+/GkW1jGZU4K7nRv+UhRhUdjfsznCIuCJGPvii7Z8992euayaJQPSEhYiUiUix4rISxiH98HAQyn2OQpY5hmrkfI0ccoSGVfj3nWq+oiqDlHVId27d0/ztC0Ir2YRjCzXoC9/0VBgHdx54LnnrshuRqBMKaaDG5rlswhq5PXZsGFCjmplyZRUobOHAKcAhwKfYkZ0D1XVs9I49t7Ar0TkCKAV0EFEngWWikhPVV3smJjcOTEWAlt49u8DLHLK+8Qp32TxEUSjzL/BoA/81gzVkvj220NhcJ5StCSlsMIi5LPIRrMIClVVtewx5P081MySCam6Ne8DWwP7qOqpqvompJdrWFWvUtU+qtoP47j+RFVPBd4gPMjvDOB1Z/kN4GQRqRaR/hhH9ljHVLVeRIY5UVCne/bZtEikWdTWogEfrLFDX1oSjY3VcPTRqTfMhlJxcNN8M1T79qu58abj8lAzSyakEha7A98AHzlhrn+g+fNY3AocIiKzgEOc76jqFOAlYCrwHnChhu0q52Oc5LOBOcC7zaxDiybGZ1FTg6qg+cqBYX0WSWlshPHjYVmG80be3Xgp9O2blzolp/CahY9gVtFQQfURDNr05KVAqtxQ3wPfA1c406ieAlSJyLvAa6r6SDonUdVRwChneSVwUILtbgJuilM+DtgpnXOVNap8zV5MYjDM90SF+P0QFJA8vtStzyIhK1fCkCHwwANw/vnp71ddqHlJ4v52LcNnEVAhELDCohRI27umql+p6p8wkUh3A3vlq1KWBKiyzAkeW73WF7Uqz8LCkpBKJxI243d/U3EEcHQkXb5pjs8iGBQCgczn7rbknoxDMVQ1qKrvp+nktuQYP8YyF/FiUs2/ZmFJSGUl/Pa3t1BdndkcEdpQrB5zkcxQ2QgLxJqhSoTiz2RjSR9VKjBZVhq9vVJVNCjgs8KiGFRWwpln/p127T7IaD9tLJ6wKKSDex0d+I7dmfZdbcb7BoLWZ1EqWGHRwnA1i4YGz8MeDBoHd76ef+vgTkplJQSDfoIZOnD99QXyMcUIhsJqFt1YAcCVl2QjLKxmUSqkLSxEZB8ROctZ7u6Et1oKSYRmEVlOUBBf9rl3UmId3Anx+804l3QGRXozslTe9Wgea5WKwv2eJ/MCAJWBzOdbCajP+ixKhHRHcP8duAK4yimqBJ7NV6UsiUlohlJfsaZ03uQRIe35RGbOhFGjTmTevO2oLNJseYV2cDc5QZfuvZsJQRUz4NRSdNL9FY4DfgVsBFDVRZgJkCyFRDVkhgoEIs1QBEGsz6JoGFNJamEhAk1Nlfj9TVT686gJJqWwZihXWFRK5sJiq27rAEEDtidUbNIVFg1OxlcFk0wwf1WyJCNu70w1v6Gz1meRknQnn1q0CD766Hd8/PyZ7NhvYwFqFo/COrhDmkUWwqJbe2O6koC9B4tNusLiJRF5GOgkIucAHwHFNLhumqjShVUAtG0djCjPezSU9VkkxWgWqTWF1athzJgjuOidV/FVFshxW2QH9x58C8BOlTMz3jcYNPWUoL3/ik1awkJV7wBeAV4FtgWuU9V781kxS3wGMJsPOITX750fLgwGQQWx4yyKhqqPdMxQTU7nujW1RZvLotBZZzuzGoD2siHjfUMz5aVjsfv8c7jiCqjL3JFuSU26Du62mESAl2M0itYiUpliN0uuccxBh/ARB+25IbI8KNZnUUTSNUO5wqKCJqgoZpRP4YSFz3nT+7ttoKZmdkb7uhZQSUdYfPst3HZbFkPpLemQrhnqc6BaRHpjTFBnAU/lq1KWNPD6EfLts7CkJF0Hd9PGeqDYwqKw94nvwgsA6HHn84wdO4Dly19Ne9+gM3goDTlsyTPpCgtR1Rrg18C9qnocsEP+qmWJS5SAiFjOp8/COrhTkrZmsdY4tf0ECpdxNo7PIhsH98yZszn99OWMGZPZfr4/ngNA+w7GHPXzzw+mva9762kgjVeVvU/zStrCQkT2An4HuNOp2pEyhcb7MAQjHdx591lYB3dSNM3h8yHN4v4R0KlTHmuUjOx8FosW7UTr1nfw00+Z7eerNhbr0aOPNGfX9KOiXM1iWbBHii092Hs1L6T7wr8EMyDvNVWdIiJbYWbOsxSLaDOU9VmUAKmvf7B+HW3atKWifesC1Cc+2Tq4RSqpqGgk02lTfBU+Z3/TwQkG03dAu1fUjuIuPmn9Aqr6OcZv4X7/Ebg4X5WyJCCZGUqxPosikq5m0XOH/+PBB2dR0b4AkeciCUwz2QkLn68Sv7+RDRkGNYlPnP0dYVG3Nu193euaVn4oa4bKK2kJCxHpDvwV2BEznzYAqnpgnuplSYXXDBUM5jc3lH0IU2IsH6mvU6DJh7+qiYoObfJep2Rk47Pw1Qc4ufUrbHvo34DN0t/P7zipHc1CgxmYoYLu/wzCjK0ZKi+k67P4DzAd6A/8A5gLzkgbS+FIqlkk91msWzeWiRMPo6ZmRnbntg9gUlQlLZkabPTh9zdR0bY6/5VyiTsoL4vD+CroElxC1/NOzGg/1wzlC3Vm0j+/O87CzpZXfNIVFl1V9XGgUVU/U9XfA8OS7SAirURkrIhMFJEpIvIPp7yLM5/3LOd/Z88+V4nIbBGZISKHesp3F5HJzroRUshcBaVKTDQUSYVFY+MKVq9+n6am9E0AltwzZmZ3/P4m/G0KKCxiyNJnUdkKrQAWLMhov7DPwtyfK9akf273Ng8G0zCCWA04r6QrLNxRLotF5EgR2RXok2KfeuBAVR0M7AIcJiLDgCuBj1V1APCx8x0R2QE4GWPqOgx4QETc7sSDwLnAAOdzWJr1Li+SaBaqqUJnxdm0WMnrNgXSMEMFKoxm0bp4Y1qzdXD7pNIIi0z3izJDrcqgv+JGQ2WkWdi+ZF5IV1jcKCIdgf8D/gI8Bvw52Q5qcF1hlc5HgWOAkU75SOBYZ/kY4AVVrVfVn4DZwFAR6Ql0UNXRTjLDpz37bLrE+CzCD2N83AfI9r7yQ3ovqJ12+oouXZYRrMx8IqCMSfjSzE5YbL3FLfR+DfjtbzOrhj/SDJX8Po0kmImD25JX0o2GestZXAsckO7BHc1gPLANcL+qjhGRzVR1sXPcxSLiBlD3Br7x7L7QKWt0lqPL453vXIwGQt9CDXgqJKl8Fkk0i7DlLgthYdX7NEl9nbbeejIAUt0qxZY5JEeD8tptdgq+T4+ArpnNThAdOuvLIBAj6FzStEJn7X2aV5JqFo7f4QwR+ZUYrhCRt0TkHhHplurgqhpQ1V0wJquhIrJTstPFO0SS8njne0RVh6jqkO7du6eqXssmzgju5IPymqlZWNU+KemGzi5e3A8Af6sueaxNOmT+e3buDNfe0RF8mU1GFO3g7tlzbto5otysszYaqvik+tWfBn4J/B4YBfQF7gPWk0FuKFVd4+x/GLDUMS3h/F/mbLYQ2MKzWx9gkVPeJ075pkcizSIYBCVF6Kzrs7C9r/yR+tpefPGXPPR/9yPVxXZwZ05dHTz/fOb7uT4L73nHjh2Q1r6uZmHHWRSfVMJiB1X9HXACsK2qXqiq76nq34h8scfgzNPdyVluDRyMCb99AzjD2ewM4HVn+Q3gZBGpdub3HgCMdUxW60VkmBMFdbpnn02LJOk+JEU0lIj7U9sHKj+kFzpbv6I9Pb+rh6qq/FcpAc1JUZ5hIBQQ67PIhLDPIgNtxmoWeSHVL9AAoCaZS3RvPlXWtJ7ApyIyCTMm40PH93ErcIiIzAIOcb6jqlOAl4CpwHvAhRqe1Ph8jFN9NjAHeDd108oU90H47LNwmeOz8PmT/STuA5RFNJTtsaVFOrm5AvhNxtlCzmVR5MmPEMFHIMaxvWrV+yl3dc1Q6Zr5LPkjldeoj4iMwNxZ7jLO97hOZhdVnQTsGqd8JXBQgn1uAm6KUz4OSObv2DRQNfbiQACuuw6uvdaUB4MmN1QaPouszVC2t5YUDYKmMY9CExWhedTzTpJoqEIPVRI0RrOYNOkw9t8/+f0YzERI2E5NXkklLC73LI+LWhf93VIIXGHhxR2Ul4bPwpqh8kQwCAtT22gC+PG3KWAkVFwKrFlgJkCKFzIbCNTg9ydOffL2D1s6S6nr+9T4nWnLCZxoOzZ5IamwUNWRydZbCoyrWcQrT5Huo1mhs5bUKJF+pAQEqMDfo2v+65OENWsKryj6CMb1WXTtupH77mvDqafG3++7BSaqMR0z1FmvHgUchVLfnKpaEpBZDNwmQm3tXEaP7ktd3fzUGxeaeE+5KhJUfL50fBZWWBQLV5b4K4rb8509Wxk7Nrs6tM4ys7rRLGLvvbVru3Paaan3TyUsamqyq5clfaywiMPixY9RX7+AJUueLnZVIkmkWaSVddbnHMIOyssPqV++69ZNZY893qOissDCIqqDIaJZOYznzcsuGgrim6Hmz/9naHnu3Bv54ouOCfdPVd8lSzxfrBkqL1hhEYdwmGkJTvybyAyVwsEdNkNlmRvKPoDNZsmSJ/jHP04oumYB2QmLvn2ha5YWtHhmqKOP3ovddoMDDviKBQv+RSCwLskRktd3q63ghJ2msQNTsqugJSVJfRYici9J7BaqWqYTIJmwxmAaNuiCktRnQZqhs1ZLyDmu5pXi/RsMNppEgoUSFgkEfLaaRbOqgsZoFp9+ehD77AMbNrxLIGDSyKlmH6n10K/eRX+4mU11zG6+SaVZjMPkdmoF7AbMcj67UJLd7tzw7bfmsnz/fYkJC4j/AnDNUGkkErQjuPOAk/U31ZUNBptoaqrEX1kKCn2ho6FihcXxx8Nm3QM0NnpHs0dexWu5wZSmIdy6tqmlGyutFpwn0oqGEpEzgQNUtdH5/hDwQd5rVyTWr/fRowds2FBiwiKJZiHBVCNkbSLBvKGKGbuQfLOQZlFoYRHHZ1FofBLp4H700ZsBOGrvNfww1ZuuPbJuN/B3ZrONHZRXAqR71/YCvKkm2zllZYnrs0h37of167+jtnZuyu0aG9c0o1YOCYRFq87r8PsDBALxw0KaHTpre2uJcUKX09UsKqqKey2LYYaqoiGiM+PmeurSoYnOVStC5fGeubN5LL36tvBOzYYNE9mwYVKxq5GQdIXFrcD3IvKUiDwFfAfcnLdaFZlMhcXEiQezcOGdCdfX1c1jwoSD+OqrzixceE/2FUuiWXTbfRYADQ2LE+xsc0PljTRfUsuWuT6LYpuhCi8sqqUhwgylKlxxBfTZvIlanzceN/Za9mIRGZnNWmjHZty4XRg3bnCxq5GQdOezeFJE3gX2dIquVNUlyfZpyfhCL+R0zVCVBIOJUz00Nm5gzZpPzJaVKTO7JyZJ6GxFmzoA/P4OCXa2M+XlDU3v5fv99yvo2bOSDm2Lm+6jGJpFtTTEmEnbtweCQVbM7B8qi3d/1tLamqFKgLS6OCKyHzAQWO18BjplZUlYs0j9UF944WU0NS1j4cKmhNs8/3z4Mn/5ZeLUBmlWLrZMlTlPHAxARUXn2PXYyY/yiqoRxSneZ+3arWHt2u706NRQiFqFieuzKL5mUVUFBAKct+gpT3nsvVZHmulR7H2aV9LVhy/3fK4F3gSuz1Odio479XeqXvhPP0G/fiaue8qUGQm1i3nzwg/me+9lMYmxSxIz1KI39uDAAwL4fPGPP2WKqcPPP1ufRc5RTUsE+3wBamra07lD8QMJC91Tr6QpwsGtKlRXA8Eg/jU+1q/v5JTHPnPD+IYOmmwMRhT2Xs0LaQkLVT3a8zkEkwF2aX6rVjzCg/KSC4uVK8HvNwJi0KAvqa9PNLw1fJlraupYseIN6uoWJtg2ZeViy4JBfARRfAk7V++9Z/b79lvb+8o5oXEWya+t3x8gGPTTo3Pq7LT5pBjRUBUSmaJcREOaRWB1G5599hogvmYhgM+aT4tOtp62hZRxynDXZ5FKs6iogIqK8IO/fPlHiY4YWqqqWs4PPxzDmjWjMq9YEs3CTXudSFi0aWOERV2dFRY5x4mGSmXaad2qCQ362GKzApuhYii8z6JCmiJ8Fn5/U0izaKIiVJ9AINEzV/7RUKVOWjaRqJHcPsygvIl5qlPRSdfBHQhECouVK9+ib99zI7ZpbFzFIYeEp5CsrzeRH86QlcxJ4LNwhUVDA7SKY+KtqvI59bEPVM5JV7PwNdEtsAI6JApCKAzF8FlEaxZ+f1NIswjiQzVB7jKfLzRtMEDnzoekPpk1Q+WFdDULdyT3eGA0cIWqJkgqXAYs+BkAXb826WbRwqK+PtYyJ1IZ8b2+1oxWnTHj96xfPz6zeiWJhurlpDhYlDDTgXmAspna0vbYUuCkW0mFr6nObNc5fhBCznFfmjlKJMioUfDyy1lVpUICEfeezxcMaRYB/KH6xKTYce73gXt94xS0rHtx9WrYaSd47rli16T5pOuzGAk8jxEWE4GxqfYRkS1E5FMRmSYiU0TkEqe8i4h8KCKznP+dPftcJSKzRWSGiBzqKd9dRCY760ZInqf58q1ebRbqamloWMH06WezZs2XMdsFAmGfBYDP1zf2WL7InM51q8PL9fWLmTfvlszCWRNoFtVODv+mhEFZdlBe3khTmIovCEFf4YRFQrIUFg88AFddldUZ/RLk1FNn88ILZj41ny9gNIuQGSqBZhF1302enOQkJdipCQRgyhTj32zppBs6uz8mJ9T9wAPAzDRCZ5uA/1PV7YFhwIUisgNwJfCxqg4APna+46w7GdgROAx4QNywJHgQOBcY4HwOS7N92eFzoqEIEAzWsmTJ49TWzojZLNCkfPHFr0PfVWNfAtHRSQ314Tw4c+dex08/Xc3KlW+mV68ffzQhWNGomnmdSSwsevR4zN04vXNZ0ifNcRYiigYF2rYtQKWSk5WwWL4cevbM6nwVEmDjxo7U1JhEED5fwGgWITNUAp9FlLBYks7orhLq2LhVKUE5ljHpmqHuBH6pqsNVdT/gUOCuZDuo6mJV/c5ZXg9Mw8zbfQzgzsA3EjjWWT4GeEFV61X1J2A2MFREegIdVHW0mm7H05598kIQ13SkiJiXfbyw2Jp5yxk58nomTBgOQCBxtz7E0o09QssbNnxv9nMybqZkbQKzmMdnkagKvXq5I8fL4K4tNUJvAvNm2LhxKsuX/zeml2w0C6Bdu8LWL4qso6GCQfD7U28XhwoJOIcwrxyfL8Bee5ljdmStxwyVXLMoRiRXcyghudVs0hUWlaoa6lqr6kygMsn2EYhIP2BXYAywmaoudo6zGHDfnr0Bb+zpQqest7McXR7vPOeKyDgRGbd8+fJ0qxeDehzcNTWmmbNnx+q/h55mqh4IGIGyfk2sQJk5M/L7WWddF7NNMNjM6JhgMKxZxHFgR7607KC8nBMalGeu07JlLzJlyvExm/n8AQiKmXyhkORqUF4in1kazG3sA8DMmbsDMGvWbsbPHwhwOO/SRxc6p0guLFLWr0Qp4aqlTdoObhF5XET2dz6PYvwXKRGRdsCrwKWqSUfWxLsrEt3VcS+9qj6iqkNUdUj37t3TqV5c1O+qwkEmTapwjv1Awu2bmoxAqW77LBs3Jp98ZdddR8WeL9vIqPABwsKiIdb/oRoWRm3bTs/uHOXURco1UYPyzMh/X8y8DOILgmbXM88lWTu4g8GshcXkhm0B+O7bg/ntb+fw2Wcnho4pwFAd63xNboZKS7MooXt1UzRDnQ9MAS4GLgGmAuel2klMKNCrwH9U9b9O8VLHtITzf5lTvhDYwrN7H8wsJgud5ejyvNFruGthC6KaOrp4w4ZOoeUFC+6IWNelS+S2tbVt+fjjkyPKmq1ZeIRFYM36mNXBYHgC+222+UfzzmWJJWbyowBhd1sY8SlUNzPdSyYkfGmGhUVd3QIWL36choZlCbb10Axh4dKBdSxe7NGsAgG3SuZfczSLEmSTExaOH+HfqvprVT1OVe9S1fpk+zgRS48D01T1355VbwBnOMtnAK97yk8WkWoR6Y9xZI91TFXrRWSYc8zTPfvklYF7vU86Q1FuvPH50HJDQ2T4bOvWMH78QaHvgYCfW28dycknhx3VM2bk0Ay1Jtb/0WxhZElOVIpy1fjConO3pc7gveLi1Sw2bpzMjBlnU1c3N/WOORAWq4nqPTmahKswpJ6dMslbtwTfyC1c1kWQalrVyST5dVR1UJLd9wZOAyaLyASn7GpMuvOXROQPwHzgROdYU0TkJYzW0gRcqOFMfucDTwGtgXedT0FQrczoB+/Z85yI79XVMH36Huy++8fO91qamqpYurRfaJva2qRyN51Khh3ccc1QzTy+JTnOS8q9T2prA6H5GlxWrx4FwJZDYkOw807cGzg6C3EaQiAHwqKKehrwzIznCIe1Hw3g0tHPMm3aZlHVjKz7brt92qzzF4sSlGMZk6rbfFS2B1bVL0nsRTsoXqGq3gTcFKd8HAVKLxLdC1f1ZyQs2rX7RcT3igrYsCEcUrt4cf/oXeLmw8kIr8+iPjZJndcMle3xLUmIGpQ3evQEBgzYGPFudSPfWnVIPtCzEETa/Z2efTo3eQ6ExaXczW1cES5wzVB1Fays6xV7q0XVa/bswey/f4KDl+B9mq4ZyqtQzZ79Z3r2/CNt226Xv4plQapfvhLoo6rzvB+gL2mmCmlprF79ScR3Vbjppme47rpXE+7zMieElu+9N/ayNDSY/BsTJgynb9+Z/P73f4tYHxMumClen0VjrGbRbGEB5aVP55qoN8GAAeYeCsTK7ZLAa4YqlGbRr5WZlKsNUTM5Om9JnyO0UgmLc86ZkPxEJXafpissGhrCz+jChXfzww/HsGQJtG69kQ8/fIf6+p/zWMv0SPXL3w3Eekyh1llXdvj9kSOuVeGjj06NGHwXTR9PZG99fewldU0Sbm4cN625S7NfKimioTZujNaWSq8HlgrVAOvWjaO2dk6xqxKLMyhPMVOnujQ2Blm3bgzffbcXwWD86W6LgzcayjWh5VdYfNf3OGYygFbURa6YbqLzXI9PMAisW+cJGiitl3+mpFv99eujQ/OFb7+FLl0WU1l5ZEwnthik+uX7qWrMpLCOWahfXmpUZHy+yGiVdN6rirBg5vYA1NXFOjbdcRhubhw3tcHEifs566uyrq9bSe84i59/jlRr166NtAMvWHAHNTWz05rcqVQIBuv57rs9GDNmm2JXJZaQz0IJBMJ9q5qaWXz33TDWrfuG1as/K3y93DfV1KkRxXPmDPaYQwujWXSeOYYBzOZQ3o9ccdFFztmdZ2Phz7DjjvDQQ2b9JjLOIhCIHE0bHgzszq1T/Gc11S+fbIqq1knWtViiczmlKyweu+IerrvuVTZujM0o+uOPO/Pf/15EXV0b55huagNzQ3jzS2WFx8E9ZXY1ffrAv/7lrlLq61+Iqs9fGTt2ALNmXZz28YuN3x8W4o2Nq5NsWQSc66MIfn97pk0bSmNjJcFgOK3HkiWpR/cXiieffINA4O9A2AzlahZ1dQsZNUpYsuTp0ParVsHSpeTEZ7ELEzl892UMGRJZ/ks+4DWOpfMOPWHoUBMZkg0lpomka4YKRJkXRCoQCb8jIMD06b9n/Pg9cl/JNEn1y38rIudEFzqRTBmmTG0Z1NZGDubz/siJfnC97wHq17Tniy9+TX0c98AT+87g3ntHsH59F+c45rIHAqbXUFXVzBeJJ3T2b49uCcCNN5pVIsLKlfGTvy1d+kz65yjSQ7h+/XeMGiUsWHB3qKyxscSysnk0C5+vgssv/4Df/GYBqmFh0dBQOsIi8j52NQvz+9bUTAOIEBZ/+hPssw/NExannx5aND4T54vjre7HPI7ldVq39cH558OLL8LixfFqn935i4Q0OR3B9fGs+WFSaxZNLFnyJOvXj8t9JdMk1S9/KXCWiIwSkTudz2fA2ZjBeWXH8uWRJhuvff/BB6Emjul56x2qqXBMTL/7Xez6oVssdo5lLrebH8d9QIMr0xgQlQxVKonUTrz1FOkYsS4QcIXVehoaVjTv3HnGjSJauzYcctrsEe+5JhgEFfzVDYwaJbz+eldWr96Mn38OC4sOHb4oYgUjUfX2eM39/eOPVztrY7MTh7ZvjrAYORLeeCO2/IgjYst+/hk++ABqa+N2UhIOxSgBDTgGJ5+bTkqWLheCjeaenjVrF6ekgqOPDmsWJW+GUtWlqvoL4B/AXOfzD1XdS1XTyf/Y4oi+Eb3pwx9++HMee2xfFi0yTtYjdl3EEL5l882hUprYt/cc0wNLQDgCxfx3Hd6BdWuaV2nV2CgTD8Hgqojvfn+4TT//fF/zzp13HB9Qbbh9JScsnJdUq/ZmQKTfbx7s77+vYvToIyM2ffHF/yts3eLgFRbLlv0HgJUrzTjXcAit0ti4hp9++jtmRDrNN0P16mXOgUeziBfd4XVuO/WZ9PqhTpHSkGyMaamZoXzOM58i4rGx0VwHNw1KU5O5713Noq5uXr6qmDbpjuD+VFXvdT7Fd8vnkUAAmprC4a/DhoXVwzZt1jNo0JfMn22in4KzfzKOOb+fSgnQGEhwOcXNqOnm7Pc5xU4ESKCZL79gMEZYbOFNnEKy45f23MahkdDvhsdhamNdgq3N/Dx//3vzIsymTTudBQtikyr/8MMJLFr0SOwOzsutdcdIU0MwCFtuGelcfuqpAqZbSfDi9AqLlSvfit4ptDRnzmXMm3cDW275upERzRUWzkkFkgsL7/bOPrWrOjpFSYTFunUlF68sfveZT75dbY15Rs8++xqnxLyDXFP1ggW35aV+mdA8b1UZEgjA1Ve/xTffHM5xB6+m2okkadVqA5tvblJ01Kwzg4oa19caX0FFBZW+psTCQpUTeNkzT7arWTjCorlpl+NoFs8+G17+4Yck6k46NuAiqvcbNxr1fcV+sHZtV1OdX/8q4fajRsENN2T/TmtqWsvSpc/F9YusWvU2tbWzY3dyBuW17bo6prhXr8j5RxoasnTc5hCvsGjdekDEOjdHlGoglDq/qmpVWLNoTs9dwvd9WpqFZ5/qdhsBqKhoSCws7r03+7rlCVdYpHqG6mMG07rConSGs1lhEUUgAN9+eyhXXfUOawKdmDfHaBZDh77HJZeYMD+pNL79dXSgA+vA76dCgjQlEhbAy5xE/6B5cbiahfui1mymOvUSx2exp+/b0PJFF33BjTf+J2L9U0/93VlKMwtqkdT72tofQ8sdO5oXeMOKxH6WpiYz5UK21a2o6Miee86gd+/zI8qXL3+VYLAOvz/OxEUJXgTxbOsaLP4j1707dHTcWN4cVqoBZs78IwBr137O8uVmCtXhwy/MjRkqXc0ijhlqlxPfBuCQQ56NG0RSsohrek5hhqqPdHCvXh3p4C4Fin/nlhjR9+6jb24OQENDOKR2ypi9UYU5bE1raqFjRyp9ARoTvQicG2aIGiHjmqNCmkVlM38Gdz4FD5XX/DXi+/jxB4eWP/74FHy+gFOH0r0FampgwYLYh2WRP/FsbY2NUJn2TCvxWbJka6ZOjZwyxfXtzJv3Y+wOCWbKiycsTuxf/CDCcePgHncuLM8rQMQfd7S/39+QG2EROo9Hs0iWONAjLFzatFmX3GdRYoTMUCn6g9HRUA0NfqfcahYlS7SwWLfRXKL6eiMs5s8fyEV/H8Ubb8AquvJHHoYOHZgd6MfUNb2TaptVjmNWVfgLt9NKas13YsMqL70Ubr89zUrHOalv2ZKoTcIP3ZQpe4WcsPGyo5YK06fD6NHhen/++XEANHZMtIcRFhVJnq9ly16htvanxBtg5ibabbfIsqamNQB8+OGBsTvEEdYAwUCQyZP3jih7bL+n42xZYC67zDOXduRbrKkp3BLXRDVlyh/MO7sZkx8BiTWLigrYe2/YwxlDkOQhat16Y4vSLNL1WUQLi6lTrWZR8kQLi5q6SGHhmpD6OwNgV2wzDHw+JjbsADiDlxLgBD+h6uNiRqBur94Xe/d/8AGMHZtZ3evxjAT/OTKXTNj0ZUaSv/nmuU55GiawIvksoq/lf/9rBhGu/8fCOFsbmprgmGPuZdmyl2LWqQaYOvVEvv8+mQ8nPptv/nsAVqyIkwMzavKj229/FIBgIDwxFpjcYK0qS8ABe9ddcOutANTWGh/K3XffD0BjY1hYtG8/FIDevT+jR4/ZOTRDRfks/H7o0CF8n8UxQy2ab+bA+P77A1qkZpHqGQpGvXgaG92Q2dJ5RZdOTUqEaGHx/KfG5LFkST/ApEoA6Nm9iV5Vy6n1Rc6nvCoySjUu/XQunzGcRseJpVJPU9NaAoFwlI/b4cqEKq/fYuPGiHVezWLmzN1YurQfgYCf+vp6VqxIQx4UwWdxeLdvI7672lAy7r0Xzj77YqZO/U3MOlc7aGxMb2zJMs/wl969/8Tpp68nch4uB1X67xTOivPOO2cDEAwokyfvGyoX0Yx/02aRxm/2yivDARg92k0wHd5n9myz3KnTbK69dkDuhEW0g9t1MiURFu0bTABHTU37bKcBLw6uzyI6dHbpUjjoINMrJHYEdzBYOuYnFyssolC9kKuuOp3uRA6UW7Vqcy699FMef9wMjZ713Tp+3v0Yzt7yw4jttkuSVbiP09kNqJ8qGkI+i5/X+Pnyy058/304vbnrqM2U2WzNOxwO50VOZOgKi2eeuSb0AmtoaMXIkXV0727GP5Ucf/oT3erSS+1RUwN/+UvybczEjdCx497JN3R41ZNoeONGYcGCdvHfv4kc3E1Bnnrqeu+G4Z5mifDgg3dwxhlTWbasLxDZqVi6NCrkOl8Obr/fHDfeLHkeAQNw5pkmdVRLIaEZat06+OSTUI8kOny+lHwVLqV155YESxk4cBxP8PtQSX9+5BHOZeLE/Vm0yCSyW7cqAPX1UGVMPxd0exFI8Cw5N3zbBbBuVVcIgp8AIkpjYyVPPmli793RyhB+hjJla37kcN6DtpFRO/X1bbjvvrv5/PPjQ2UNDa2orq53zp35ufJGIGAepOXL+fyBM1m2rA8nnrgg6S4NDXDnnYnXP/443HRTB/z+DrRtm2zOrjBjxoSX73KGXTwTL0NKAmHRvl0w1oyQAwdxLgkEKpk/f/vQd6+wWLx438iNcyQs9p14H0f8Yo1bgbBmEQzZaRMewtV2WgoJhYVrS3PeH4GmSHNwu3alpz6V1p1bAgQCA9mwoRMdCU9S04YazuEx3uFwPnLmbbr/yTbmB3d+7Pv73YEeHk5d0NAAX38dm95GRBENC4uJE4cze/YuceqRuRkqghEjQiouQGNjNa++egmzZ+8aKnv//TNYsGAYUGLCYtw4o6I3NhLY2Irf/GYBK1b0CY14j0e8a7VqVdi4/eGH8NxzrkM/Pb/BJZ6ENuvWmf9xnatxXm7L6M5WW0bWd9q0PUtOWETjNX9MnnxhaPmLL47NmbC4pO5f3D7VGdmephmqwdEKv/iiZQqLmHvEFRZOwsSAk9r+nntM1F379slyuBaHvN25IvKEiCwTkR88ZV1E5EMRmeX87+xZd5WIzBaRGSJyqKd8dxGZ7KwbIWlN6ZU99fU3c9FFX9OKOlbTCYApziR9h/MeA5kJwFuftjU/uJsds6LC2I4cVq82AR7/+x8R5UQIiyCVlQ384df/iqlHtmYoT0Pg5ZeTbvLgg3cyZ87v2HXXFD6LQju4m5pMgrlzzgmlrgaYPt04XGsnbhazi1dYBIPCzJm78f33/oj1gYCJLkk9z7NhV0euqsK3jvvk44/jbBjn+nRnBf93telIrF3blUmT9mGPpzeUvLAQj6v++7CiS+vWG3ImLICwFpFKWDi48WbDh7csYZHQZ+H2OpzOZnX7JwH48stjqa1tQ0NDx6gZDYtPPu/cp4DDosquBD5W1QHAx853RGQH4GRgR2efByQc0/kgcC4wwPlEHzOnBBYbG6KfAJ1Yy7/5M6MYHlq/hTPR0QUnrYjQLKiogCnhSY3cWP+GBuDaa0PlIhpKKS6itG+/it+cd31sPbI0Q0WQ4AC7VIdTUIwaBSefDFtvneJYhXRw7703fPopdO0aISw2buxIfX0r/jf1jJhdvMJi9erNmDFjCKqRwsLM87GC6dM/TzlICoBXXgHM1AqfOdNR7LdfnO0SHGvObFPe0FDNwoUD+YPvmZITFn/72ymceaabrhzEM0B0r73gj3/8lg8//B1vv31OboWFe82CwVhhEWd7dZZ/nXgOspIg+r5KmKI8ygxV1eo7wKT3OO+8cfzww18SPnJpRTDmgbzduar6ORAdG3QMMNJZHgkc6yl/QVXrVfUnYDYwVER6Ah1UdbSaX+Fpzz55Yc77Jp2Dm5rjz9zNcD6P2EYR7n+pO8yda1QIgC+/hEWL4HWTkG3NGlMcHeY3feKerJ+/eUhYVFTEjwNsamqmGQrA5+Oaa+AGro0obl2/JuL77DgZLEoCVbqzPKIoGPTFNUd55aLf3xTjIPT7oc4JNuvY8QfS4kST1O2kk1LX06WxMRwqu2U3E5HW0NCKysp6E6FWSKGbxrkGD/6Mbt3CYdauZjH/k71p3x5mzhzCzTc/y6hRJ+VWWLjRP4GAOabXwe19szpRC0HX79e2NDWLQKCWzz6rZP78WCsBJBEWUfN2BIN+5s/fnqamzUO/xQEHRO5crAy0he7mbKaqiwGc/z2c8t6A14O50Cnr7SxHl8dFRM4VkXEiMm758uWJNkvK/71qIpL8adq1efPNyO/fmR7CaaeZrzNnRq6+5++PMvN/B+AnwLffHsqUKb+IWO/2THKiWdTWcvPNcB3/5BCM/+Jxfs/tXM6WzA1tVmK518KocjeX0hpvxllf3BBa73uobdu17LXXm4iEb51owZuJNbNLF2jdGv4vUcJYz5tAVajEvAj+cJDJFNrYWG2EBaSc16CQBINKhw4rWbPGPIYiUCEmKuf7FYOJeYTybYaKdnCLmAGERx1F0HlVtWlTmsLirruqUW3i+utrYzKxC8FYaRFlhnJxx+VUVERe6iVL+oaWVYszN0qp6MRxAxKTlMdFVR9R1SGqOqR79+6JNksLX7rZWKO7nY7wcAeURb8bAvjxE8BPgCefvIFXXw17US+8cEnoJdZsBzfAU0+FFp/nFB7hHH7Pk+zN18ylf7hOqYRFsRIJqlJFIzWEI7uCQT9bSfJ5uK+88l0222wBeDTC5lxLEehVtZxFr37NHnvA3XfH1tPF5wuafGFAk2MGixAWDz+cfUVyjKpwxBEbeOaZv4XKaiabcUVjJx7qWuHC5CjdBwDTzCRLjBwJ8+cn91mIUOkIsY4dS1NYBIM+6utbUV1dGzcqL8ZnEWWGcqmrM/d6RUWk/6iqKjwG64svNg3NYqljWsL57w5mWAh4k2r3ARY55X3ilOcdPwG4+ebkG339dXiuYJcdd2TUKJjjvM/WL49Mpx0SFr1MzqmqKvMS+dvfXmPq1LDjttkO7ii6sopzeCyirC0mBCotzaIYiQTjCKl2wY0c5XszzsZh3Oy03pl/0xUWieTiFWuv5qS5tzFuHPz5z7E7PfrozXzzzRHceONzIWdso5M647PPTmDcuF+mV4F842lgU5MJna2vbxMau7Bu4pYAfPttnPrmUrOIHtiTwmfRy3nst9iiNIVFVZXJ8lBVVRs95bkZsR69Q5QZauPG9oD5PcAEBDY0hV8Ab7xhEluuXduVK66IvUmfeALuy/PUNIUWFm8ArnfyDOB1T/nJIlItIv0xjuyxjqlqvYgMc6KgTvfskxeqKsyb00/A5M+JnizY5eSTjQewsxPQ5YbLrFvHAQeYxTZsZPsZr0XsFhIWW/cDCPU4a2oi5+5Oywy1IGrsgWfqynQ4ptvXoXOVJHHe3BXahPg0bgzrGc6d1apVjbN7eN5u91r+9a/vcdFFiWetE4FruSGm/Bwe4wA+TVjP5567iquuepvPPjsxJCyO3G0xj/N7/vOfa/jf//6U8JwFpS7cefEG6f3JqZ7060og4KOpqYoYciksoknks3BXJwmbLgWqquDTT09m5szdiU6KLGhsIsE2bWDgQPMf+Oc/X+DBB8PJ4MZH5ZscOfJ6DjhAOfbYFaxc2T7m/P/7Hzz2WExxTsln6OzzwGhgWxFZ6MzbfStwiIjMAg5xvqOqU4CXgKnAe8CFGvbinA88hnF6zwHeJY9UBc3DFPJZfPIJ9O1LRHfhttvggQcidxwyxITKeLzF1dTTUBc12IYNtKEGfyvTg3CFRWNjpKMrLQf3889Hfs9QFXGngm1JwkICGONk9NME3HST+V9dbXqtCxbEahbffntoivk9wr99MMoK2pF19Gm/hl/uuoymJhg9Gn78MbaerrAYvMUqfs+TkQfv1i3puXNKvJdz6/A1cX/3P3EvfzSZyfF1bmLduq7xcxLlaD6LhOvi+Sxi9itdzeLuux/g/ffPjFknaOy9fPTRMGMGbGk0uTFjjuCll1KkIHBwTdtTp5pL86tfmdsqSeb+nJDPaKhTVLWnqlaqah9VfVxVV6rqQao6wPm/yrP9Taq6tapuq6rvesrHqepOzro/aVoxj9njDxrbaGh+iPbtYd482H572GUX2GwzuPzysEbhpUOHCKHiI0jdushh/HPpz+38FV+VeXu5wsI7KY77zKR890clC8z0rV/htDHlsIMi+ixcFtKbFXRl0F+h7/MkzU8yffoe/PGP33LhheEBiJn4LFxhEYgz10fH9Qv54Pse1NfDvvvCk09G1nPIEELOWBrjzFBYSGGRiN/8BrbdNqRZDGBW6H085/MzuOKK90KbDtxsTXi/5moWyfZNMc6CkNwoXWHhcvyvY5+XVI/Q/529NvkGHsyEgDXMmjUKMANOf/97+Pe/TW60E05I+1AZUSoO7pJhCycoqyuxM6UxdqxxxiUi6o5YSTcebTwz7qauGcr1WXg1C/edn/IFF+1cTzYBeBwq2rWOOF/J4bmevVlEV1bRfg60WgqsjP192jeavkdNTQdmzhxCbW1YXc9E6UomLE7gFc7kSaqq4MgjoV27yHrus09Ys3CFxdfsxffsYsqq4ph3Co3PB8Eg1dVwHf9gKGNDN8H61Zsza1Y4P7vfp7H7ZkuyGzqZz0LEo1CUvrBoNXd6xLqILLsJuOPKFfQgccrqs3gitKxq5lfp2PEAunX7mXPOgb/+1bwOpk6FLxJbWZtF6WWrKjLvcRjfMIzWxJnnOdWsOkccAW+/ndZ5/NuaHFOTJ+/NlVe+Hcpqa+YAd7ZJ9YKLjvhynSWJGDLEeM4cKn6cARzcohzcIeIIiw5jP2IJF7G589Cdckp4nfc9lSrhYIUzv0g8YXE9Jo8XlWe5Q2pgfJC/8i/qqcYXuDAsLBwn5l58Ez5AKQzKc4RF27bwD643ZY2N4PfTGIhsc0zew+bUP9nzE0+z8NDhB6jtA35/h5h1pYBXWFR/9zUQzrdlfBbhNi1dagbDHnig5xFWpYrEudcHMpNzeIRHOZd77oGuXY9i7dovmTVrHZ07m9EEP/yQo/FZCSiBO7e06M0ijue/2e3cKjKfyxG8zZDuc2O3O+YY/G3NtqtXb86YMUeEesG1tWHbY4yw2GILOPXU8PdttonMiZ7khT6HrWI8YDtgTGYtQbOIwWn38OGe8Q/BIJt5sgXvv394c+8DdM45yU/rahZNyfpSNZ45z599ln9xJXfzZ6ShPmyGijfxQq6E7n//Cw8+mN2+jrCIuL6OFtQQJSwq/EF8BGjjrwvvmy3J3mLxHNyeazXw3zDkD1BdvXn2588jVV9+Elp+iPMikq2ZENjwtZ461cTH/OAdGxoM0ipeB9VdjS8Uzh8Mwtdf78ANN7yBSFgoPfigFRbFYdKk1NtEEyUs6iraMWFFnPkPli9PqDbU1pogK4hJHBs7U5nPxy3d7mAkThRUkhfRDLY16y+6KFR2IQ8wbFiBlYa1a+NOlxmXNDSLZcs8QWGO86WLY0L0ylXv5W490dPTj0MyM1SIJZ6ZCAcONP8rKhCfxJihIsjFxV66FI4/Hi64IDt/kissvPVzBNtpgyPv+wnzu7KBdqw87pzwvtmSSrNI5OAG/I3QLs6MtqXAyJHw/ifmDX0o71FJU8zENt6fyY0xiHC7LV/OA1zAQGaEirwRYD6C4cCLoLn9XnnFuC1/+gmGDjWmqKblq/HX5icrqBUWiXCiFDIiSlh80jScJq2gDscf4Xqevv6axgTTJTY1hf3WXbpErVy4MDJkdNEi7tcL+AAnJj7Oi0ic3khvnIMecUTE+tGjYwehx5BLB/fTnmlF471M0z2vIyzaywbWr418yXzJPtxy7UY3KhGI7G21OunopKc9jWeYzrYRmYe9nM5IDjolbAKctbaHyR/Wrx+iweTCIhdmqK+/Di+vTeIYTSSYXGHh1Xycug7rFeuTa00drbQ2vG+2NMdnUcKceSbcN9EkDTsOJ1TezfcDHFv9HjtN+A8MGwaqIWFR51Uk9tmHQ/iIGYQnxDlG/xdavpgRVGOe/VatwnES69ZBv34mnf6WW0LT2+9RsTw/Q9GssEhEKv9EPBxh8SEHM4rhdHBeNhMZDNdcE5EHpsGTCvpVfs3RmBwBgVXhhz8i4MrtqXge8Pm9h/EzfajBeSvGeahqac1M33Zsx3To2hUOO8yoyGedBX3iaD2JyMUDqwoXXxz+nmrQYyph8f77DJr2An1rHYeiM3R2e6Zz5cWR0VIRwsKfPF1CF1azLTPxJxjFv4webJgbjlN89JOtOZx3oaoKH3GExfHhOURych1dYbF6NXTqlPn+8YSFu5zIJuk60nKpWXj9Tm3bhqcmTDI4j3vvzf78zWD16tSzYC6kN+fyiPniERaXtn2ULfUnMxZLJL5mEYfb+GtouQ21XMmtXMpdnHpquF/q7TuuXAmfs1/YDJpjrLBIRDaGPyep4MF8zHA+Zy79GMVwdmayGavRr19oU39l+AH4da8x/NrxkzQtXEKPDqbLsduunofGzWh79tmhos1Yyrk8zAU4Yz7ivIiqaWBA3WSqZ/4AvZ20Wm3bGsHl3GmBQPIOas447rjI79dfn3z7ZMJizRr4+mse4xwe3udZU+bNqf3aaxFdtx49wqtaByOnnE3KccdFPPgANbShzbK5JjHgrFm8Nn4L6mgN1dWIBmNDZ//2N/MB02loLsOGmWHk2QgKSKpZRAsLN89VToRF9DN1gzP4cfBg8wPV1UX6UuIJVm9no4D8858wYKsm1q2NvCe9Yec9WRyO1fLcM1dvuIqruMVsXF/PhO/MTomExSk8B0BrIjfowXLu4jKqq8P9Tq928uabsIjezGGbTJuXFlZYROO+jLPJtRGVv7qean7Di7zAyaZXde21oZfPkEGNtKaGRzgHKitDEThNDUGOXP88fVhA9zbOS2358vCxd9opdPxqGniY8zhkpyVm8oVEYZmVlTBgQGRZVVXoZfHnPztybM0a86B+8EHue3AbN4Yy8ob8JgccYO7wRE9NMmGxenU4p4rGGfR07rnm49CpVfipqtAU5i+HlXRh9v8mE3g4MjCghja0ZaO5pgMHMnuVk16kqgrRQKxm4febt40qHHNMWudOyvHHm6D6bPH70xYWIaer++LO5TiLQMB0XA4+OKx1eE13pWJ+amxkm7su4A9r/02Hlx+PWOVOigWwDk+klpuNGtjMv4KlOKl8WrXiz2cYFSV027t5shz+xRW8vcV5dCPxKDtXWHg1i4jgyDyMjbLCIpqHHjLd7Gweii22MD9SmzbQuzedWMNSNmcxPc3LubraRLF07IhU+KmhrcnXVFERISz+r9UDPMupoRsueMe/I8/h4moqp5xist327Jl+XT2aRafR77B2TZBg5y6w1VZw6KGmB+feic298Vatgt13N8vPPmtm8dtzT5g82Qw/9Zpp4rF0Kfzud5Flq1eHU6yMGRP/9/okHKHS6k9hjYy905uD+2VOZACzWXJ75FyqFzOC3/NEaBrEbZhlVlRVIY88HBM6m9MkX7kgAzNUqHfr3gu5DP1tagqPCvcKi2INAk3Eo49yAQ9yG1eYEXAeOnYML3fy+ric7NMAp7V9jSsIpy4/sMGMOa5Z7Gz/XngQJJg5c45oeoPqJKG08YTFln2V/fiMr/hFZC6XHGGFRTR+vxmJ3RzWroW5c2lFPVsxx4TORff6vS8Qr7BoVHb0TTNzaHz1FY/5/8hOt51GIxUmAYy3t+WmvPB6cr2MGGEGEsajujr0gug87kMUH2vpaObocDnvvPByc3p5f/mLSW0A4cEPFRXG/HD//fDLBEn23JdGp05w662RMa+LF4fzv3+aIGeT+8DMnUvrxUYLeY9D4auvIrdbscJcp6hQ19lsQ2tq6LlyckT56TzDr11HJjCJQayhI/z4Y6TPwj1evmIZ166NjL+cMyeiR5uQFGaoLp4BqaHerWvHy7WwcCP84gmLUnFwe6/xSy9FjKUSgSN5i8FMiNxnxIjQvLxHtBkVTuI5bBi/cvyTvae8b8ri+Ymi52OOIkJYrF0L8+ez0w8v8Bn78wtG5yUe3gqLfFBREXpBzGEbruaWiJw8QORDt2ZNhGbBRsf8dMop9AouYBk9mHbGv2JNGB06mBs5utftsvXWsMce8detXGleGGeeyek8zWy2pkO18/I4+GDzP1dDQa+80njk3n473O7GRvNAXHABXHpp/P28L40+feCRR4z5xU2R6pLIdr90qXlo+vdnM5ZyPK/EV+132sloOtXV8G449dh13MA4huCLlxXfM7lFa+royDrYuJEtmccvcBzQbkxvNsES6XDwwbDzzuZ3nDLFjLvxhtCliobydkvdcSMzZnBGaH4yOB9nLIcbKhx9HzcHM89tamERTb6GKCfCG30GcNRR4eXnn+ctjmYCuxLDiBFwyCGRWR+qqzmRV1CEU9/4jfHbJNECfqQ/Uz0D/DyHAZyfcOhQEwrl7ehFRWbmBFUty8/uu++uRWf0aNearjpqVOS6uXPD60BHs6cezts68z9jI8oDiNb8/dbIfUH1j3+Mf07PvvrWW4nrduedkduCal2d6tSpZv1hh5mypUtVjzpKdbfdsr8OyeqZjBtuMNs0NkaWz5unet55pnz4cNV9941tS6qPS0ND5vvedFN437feCpeLxN++oSGnly7EL35hjj92rOpf/hI+34wZZn2HDvGv86WXqnbsqPrNN+H1b76pOm6cKmgQ9Cpu0uXHnh1eX1Fh/o8Z07w6e6/Lqaeq+v2qV1+t+tBDpuznn1X/9S+zvHGj2ef44yP3e/TR5MfOJStXxv9d166NbY/3+oPqFlukvpe6dVO98caM790lS8ziffclqEMzAMZpnHeq1SzyyfaeHkH7qLTCUeM4hjGGdziSARo5tZ4PpfXWvSL3bWpKPHr33XfDPe1kGQK9uTDAOLWrq8N1dnuam22W2JSVC2prTWK7Sy4Jm0Xefdc4j6+7znyP7mH27WvaX1FhjMbr1kWazNLBDf/yOvKHD09vX1c7qqw0CaJcNI4G0qVL/jSLhx82I7GmTIE77jBlP/xgBgnOnRvpfXU58UQze1MgEGmGOvroUDp+AW7mGrpVefZvajLO6F12yV39Fy0y9aiszMxncf/9ESOkgfz4OYJBcx+qhs3Gbrj3gQfGzkd8xhnGDDhwoNln/nxzbyejffuwZvH992lk9QRuvJGOHeH+y2YzfN8Cpm6PJ0HK4VMSmkVTU1jSz5wZux5UN988skfw+OOxvQS3J5suRx1l9nv99fTqduihsesPPjiyDkOHZlaHVFx/fbgn6Z5jiy1U162LbX8gkPg4p56q2r+/6vz5iXtil12mumhRZNn335v93e99+6pOmRJ//9NOi6zneedF1iHRec88UzUYzO11i8df/2rOd/HF4bJ33ontaXqvbXV14nrvv7/q4MHmN/eWH3RQbupbWRl53BtvVH36abM8a5bqLbeY5Zoas320ZgGqt99uNMsVK1RPPtn8nrnQLNatM1r3zz+rtmsXPub69Uajra9PfN0efjj+MZNpCYMHh5fnzTPbB4OptY3//c/8j76WoHrHHc26BFjNogh4ndgxw7ExDtqIBDGEe7xeJ3umkQ1uT1yT9Lb8fujlaCzxopGiI3i86S1ywfDhJk3mW2+Fy3bYASZONMve3ngyB2fHjuaaRY+W9vvhnXdM8v8774TNo3IKXX+9GYvh8sknRqty/SGDB4fXPf20uVbnm9nKQuMDXD77LLZeb75p8pcXwjl7yy1m2P8994TLIoYHO7Rvb3rk99wTd/IowFzHxkbzO7gapTuPrCdVTLP4JirdSkVFZj4LMNME7LyzucYvvJC5ZhnN+vUmz06HDsYf1bt3WHv58UeTXrhvXxOoEp3t+e67zTU988z4x3bD3uOlp/f6gCZMMP9FzHicZLnG3WyY3vv+lVfMtUs4WXwziSdByuFTEpqFqurChcl9B6rhniGo/u1v4V6pW/bZZ5md81e/Mvu99lry7ZYuVT3hBNVly2LXHXFEZG+ldevM6pAurl/nyCPDGsR776k+/3x6PcWrrzb29B9+iKzvOefEbpuolzZ+fHibyy4zZa4fpF279Npx1lmRx/z88/T2yxfPPRdZHy9Ll6p26RJet3ix6rRpqrW1Zv0ll4TXebWVTJk82fhT4vHgg+Fz3Hab6iuvmOVJk1Rvvjm8buJE1QMPTN7L9n7uuSezOn70kepvfqO61Vaxx2rXzjgHopk1K3K7H39Mfo6bbjLbtW+fvO5PPBG53+LFsdt07Jh4/xxBAs2i6C/1fH1KRlikg9ck5H4++kj1kUdUly/P/HgnnKAhVTVbXNX/3/9Wvftu1TlzYreZM8e8eJrLiy+qbtgQWTZqVHoPwa23aoRwffhh1W23Vf3229ht3eMNGxZ5rb1mrkmTVLffXvWpp8y6ww9Prw319WGhO2tWevvki6++Uu3aNdy+Ll1itwkGjekunqlywwbV7bZTffnl5tXjoINU99or/rqJE8P1u+CCsFnl9dcjTTAiqgMHmuXHHze/c7duiV+Yq1dnVkev2bd7dxN4kg7vvad63XWq06en3vaeexLX9667wsu33Ra77847m3X9+5v/s2fHHuOii8y7Ike0eGEBHAbMwEyvemWq7VuUsFA1D2ifPuEboDlRJ0uWmBuoOVE4n32mesABxiaciN/8xvhckvkUsmXt2vSExWefmYiat982giI6csrLTTepnn++qe+TT6oefbTpUccjGDShJosXZ92EohGt5Xz6aXHqceqpqv36xV/X2Kg6aJCp3zvvqL7/fuIX6iefqO66q4nWUzVCLN52F1yQeR29vqjZs7NvazJef934ga67ztTRPd9RR0Xe51dckfgYGzcaAasa2+4ck0hYiFlX2oiIH5iJmbd7IfAtcIqqTk20z5AhQ3ScZ6KfFkFdHbz6qrHFDhyYn1jpXPIPZyKgv/89P8dvbDSpTnr1Sr7dkiWxPolNmdWrjQ1+1ixje99uu9T75IMvvzTJAX/968TbqBobfX29scPfd58pf+IJMxDzuOPM/+h9Zs+OTGFTW5v9GJCGBhOFVKznbeVK45fs0SM9H1dNjfHxrF5tfCjZ5ghLgIiMV9UhMeUtRFjsBVyvqoc6368CUNVbEu3TIoWFxbKpEwyaT75GvVtSkkhYtJRoqN7AAs/3hU5ZBCJyroiME5Fxy5cvL1jlLBZLjvD5rKAoUVqKsIinm8WoRKr6iKoOUdUh3aPnp7ZYLBZL1rQUYbEQ8KRbpQ+Qn+mgLBaLxRJDSxEW3wIDRKS/iFQBJ4OTutFisVgseadFGAdVtUlE/gS8D/iBJ1R1SpGrZbFYLJsMLUJYAKjqO8A7xa6HxWKxbIq0FDOUxWKxWIqIFRYWi8ViSUmLGJSXDSKyHJiXp8N3gySzqZcH5d7Gcm9fISj3a1ju7YP4bdxSVWPGHpStsMgnIjIu3gjHcqLc21ju7SsE5X4Ny719kFkbrRnKYrFYLCmxwsJisVgsKbHCIjseKXYFCkC5t7Hc21cIyv0alnv7IIM2Wp+FxWKxWFJiNQuLxWKxpMQKC4vFYrGkxAqLTRyRdKbmslgsLZlcPOdWWEQhhiNFpEex65IvnDb+WUT6aJk6rURkaxHJcp5NC4SmMy7bDkW5P+u5fs6tsPAgIscC04A/AY+LyI7FrVHuEZHTgU+BXYF15fYiEJHficgU4HbgVRFpMckySwUROVNEvgcuKXZd8kW5P+v5eM6tsHAQka7A74GzVPVwoBLYwVlXFi9UEdkbeAr4i6qerqrr3B5HObRRRI4B/gj8QVV/DbQBznHW2Xs9DURkO+AC4C1gPxHZSlW1nK5fuT/r+XrOy+YGyAYRae/52gC0AzZ3vgeBniLSs+AVyyGuKQFAVb8CxgDbO+uuFJGjRaRdSzVHiUg7z9eJwBmq+o3zfQRwLICqBgtctRaD9zlQ1enA6cBdwFRMz7vFX79yf9YL8ZxvssJCRK4EvheR20TkVFVdDzwL/FZElgHzgUHAP4GDi1jVrBGRG4DrRMSbFOw8YKSITAA6ARcBtzs9yhaFiPwVGCUit4vIKao6l8jkkVsBXxelci0Ez3PwLxE50ymeoaqrgNeArUVkP2fbFvm+KPdnvVDP+SZpzxWRA4EjgEOAfsBzIjJRVZ8QkdXAQlX9s9MbuQ7YGviwaBXOEBGpBi4DzgC+w9gtPwBQ1YnOrINjVHW8c4M9CmwJTC9SlTPCMSPciTEznQ1sA1wuIp+r6s8iUqmqjZie48QiVrWkSfAcfKeqk5xNpmHs3n8EPlfVoIj4VTVQlApnQTk/64V+zltkTyEHVALfq+pPqvopcA9wm7NuFRAQkc5OD0QwL6WWRCPG5rwD8A1wgIj0d1eq6gOqOt5ZXo5pc5diVDRLNgIfqOpJqjoB+BiYDPQBcAQFwM7AFyIyUESud4SMJUy85+AWz/qNwCvABhH5p4jchnnhtiTK+Vkv6HO+qQqLNkBXEWkFoKq3ApuLyNHAEmAzjMp2J6ZHMrZoNc0Cx748U1U3Ai9iXqJDnZ5IyMklIl2cNg4Cvi1WfTNFVeuANz1FTcAuwGK3QET6Yh6MfwAvACtUdWUBq9kSiPcc9BSRE53vCtRhhO75wHJVnVOsymZJ2T7rhX7Oy1pYiMj+UXY8AFT1NYy6eZSn+Hbgr6o6A7ge+BnYAPxCVb8sQHWzIkkb653/c4EvgeHAdk6ZOj2QFzE9r+GqOrtglc6AJO1b7/naFVimqvM9ZdWY33gZsLeq3pffmpYuInKaiOwcXZ7gObgNY9pwuQWYAvRV1dvzWtFmkGEbW9yznqR9hXvOVbXsPsBhwOcYG10HT7kA1c7yyc42/ZzvfYEHgLbutsVuRzPa6HOW/c7/DsC9wG+B04CjnfKuxW5Hjtq3D/AfZ/lQYD/n4ehZ7HYU+RoOxvhs/gcMjrqGyZ6D+4H2zvdWxW5HntrYIp71FO0r6HNeNg5uR+US4DfAw5hY+5e969VcuXoR2QojbXcA/iZmANLRwFw1Kh3OtiVFBm1Up42rgDWquk5EZmFeAiuBiwG0xMwy2bYPIxyqRORBjKp9pRq/xWI2bY4A7lfViDTUaT4H651t6wpc50xpThtL9ln3kKx9BX3Oy8IM5b5E1NjwFgFPA7OddSeKSB9MTxMRuQYTg7wPJqLmSWAo8LGqnleM+qdDFm38CthLDNsBlwK3q+o2qvpOURqRhCzbt7ez+47AnsB0Vd1bVb8oeANKANdG7WE7jF0eMWkfDhORTs73KzD2+Rb3HEQVlVUbs2jf1xToOW/x81k44WEHAV9gXjCrgT8AZwGdMVEyFcBSVT3P2f4/qrrac4ySDgdsbhsd557f7UmVGjlo31HAl6q6pgjVLwk81/Az4GU1IcS3AHOAIzEvnM6AH7gR2At4sYU+B2XZxua2L+/PeTHtcc39AMdhvPsHYHoN9wMDgV7ArcCuznZdMWrZYM++fkrYVpmjNlYUu/55bl9lsetfCp8E17Avxlb/Caan6d7zn+DYsj1lLfU5KJs2NrN9BXnOW7oZak/gQTXx09cDczFRDouA61X1ewjZ7P6LkcquySOgzpUucZrTxqai1DgzmtO+xrhH3PSIvobzgKtU9QVMb7RSRDZT06P+GhMd1NKfg3JrY3PaV5DnvEUIi2g7nuf7jxjPP6o6DxN7315EfqUex5yIXIuxa09zti25G6fc21ju7SsEGVzDN4DuIrIPJky0EbjKuYYnYMwcJXkNy72NLbl9LUJY4Dg2XTwX6BWgRky2UTDRL6MIZ5DcV0Q+xZg1jlfVpYWpblaUexvLvX2FIJNr+Alm3MD3mLES0zED1A52tbUSpdzb2GLbV9KhsyIyBLgCWCQiLwOjVTUgIj41UTOrMcnOzheRN1R1rZgspO6kN3OBC1V1ajHqnw7l3sZyb18hyPIatgXaAqhJCvhQseqfDuXexnJoX0lqFk4Y2K2Yi/MWsBSTKrkvRKRLbg28j5HCj4hIL0wyrQZnuwWl+pIp9zaWe/sKQa6uYSlT7m0sq/blwkuejw8mVKyzs9wTeB5o51l/A+bi7orJAXQjxnzxAM6IxlL/lHsby7199hraNm5K7SuZcRZipjncAZikqm95yvfF5J5fhBlM9zomxOwR4Dr15DoRkTaqWlPIemdCubex3NtXCDaFa1jubSzb9hVbWgHdMXlPPsdM2LEMOM6zfkfgAGf5LMygrf6e9b5it2FTb2O5t89eQ9tG2z4tCWExDLjc8/004OsE226FcQL1aQkXd1NpY7m3z15D20bbviIlEhSR0zFTGY4FxgM/OeV+zLy/kxPs+kuMU95Nclay8wKXexvLvX2FYFO4huXexnJvn5eCCQtn8MnmwHOYCdLnAOcAl6jqUnFytojI9jijdJ39qjBZRf+FyTv/V1VdW6h6Z0K5t7Hc21cINoVrWO5tLPf2JaIgobPOxVOgPfCzqh4EXIBJrftI1Oa/xAxQQUS6qWoDxvb3d1X9laqW5DzR5d7Gcm9fIdgUrmG5t7Hc25eMvGoWIlKBCQvzi8g7mMk5AgCq2iQiF2MGqQxX1c+c3TYAP4nIDcCvReRwNRPIT4pziqJT7m0s9/YVgk3hGpZ7G8u9femQN81CRIZjbHidMfMS/BOT3+QAERkKoaHuN2ASZ7l2vt9jpHEHTOTAgnzVsbmUexvLvX2FYFO4huXexnJvX9pk6xlP9QH2BU7zfH8AM+n7mcB4p8yHsf29BGyJyaR4N7Bbvupl22jbZ6+hbaNtXxbXIY8XuA1QTXh+2N8BtzjLE4CLnOUhwAvFvhC2jZte++w1tG207Uv/kzczlKrWqGq9hmelOgRY7iyfBWwvIm9hhr6Ph7hTCpY05d7Gcm9fIdgUrmG5t7Hc25cueQ+ddWx3CmyGydEOJrb4amAn4CdV/RlKL/d8upR7G8u9fYVgU7iG5d7Gcm9fKgoROhvE5HBfAQxyJPC1QFBVv3Qvbgun3NtY7u0rBJvCNSz3NpZ7+5JSkESCIjIMMxXg18CTqvp43k9aYMq9jeXevkKwKVzDcm9jubcvGYUSFn0weVL+rar1eT9hESj3NpZ7+wrBpnANy72N5d6+ZJRMinKLxWKxlC4lOVOexWKxWEoLKywsFovFkhIrLCwWi8WSEissLBaLxZISKywsFovFkhIrLCyWHCAinUTkAme5l4i8Uuw6WSy5xIbOWiw5QET6AW+p6k7FrovFkg+KMge3xVKG3ApsLSITgFnA9qq6k4icCRwL+DH5g+4EqjADu+qBI1R1lYhsDdwPdAdqgHO0hc2kZilvrBnKYskNVwJzVHUX4PKodTsBvwWGAjcBNaq6KzAaON3Z5hFMquvdgb9g5kywWEoGq1lYLPnnU1VdD6wXkbXAm075ZExCunbAL4CXPZmtqwtfTYslMVZYWCz5x5tDKOj5HsQ8gz5gjaOVWCwliTVDWSy5YT3QPpsdVXUd8JOInAhm4hwRGZzLylkszcUKC4slB6jqSuArEfkBuD2LQ/wO+IOITASmAMfksn4WS3OxobMWi8ViSYnVLCwWi8WSEissLBaLxZISKywsFovFkhIrLCwWi8WSEissLBaLxZISKywsFovFkhIrLCwWi8WSkv8HWsD81j4sOcYAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "#test plot\n",
    "ds_all.sel(cloud_layer=1).cloud_base_altitude.plot(color='r',label='cloud_layer = 1')\n",
    "ds_all.sel(cloud_layer=2).cloud_base_altitude.plot(color='b',label='cloud_layer = 2')\n",
    "ds_all.sel(cloud_layer=3).cloud_base_altitude.plot(color='y',label='cloud_layer = 3')\n",
    "\n",
    "plt.title('')\n",
    "plt.legend()\n",
    "plt.show()"
   ]
Sonja Murto's avatar
Sonja Murto committed
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "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.8.5"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}