generate_timeseries_data_chl_poc.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "5360977e-b3d9-48ba-8f1e-72e990b6d349",
   "metadata": {},
   "outputs": [],
   "source": [
    "import math, pylab, csv\n",
    "import xarray as xr\n",
    "import numpy as np\n",
    "from datetime import datetime\n",
    "from datetime import date\n",
    "from itertools import groupby\n",
    "from collections import Counter\n",
    "from matplotlib.path import Path\n",
    "import matplotlib.pyplot as plt\n",
    "import earthaccess\n",
    "import xarray as xr\n",
    "from xarray.backends.api import open_datatree\n",
    "import matplotlib.pyplot as plt\n",
    "import cartopy.crs as ccrs\n",
    "import numpy as np\n",
    "#%matplotlib widget\n",
    "from scipy.ndimage import generic_filter\n",
    "from scipy.ndimage import gaussian_filter\n",
    "from scipy.interpolate import griddata\n",
    "import seaborn as sns \n",
    "import pandas\n",
    "from xarray.backends.api import open_datatree\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "394e47f8-6837-451b-9523-41b84a46df1e",
   "metadata": {},
   "outputs": [],
   "source": [
    "# %matplotlib widget"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "a143ca43-b35b-4b30-851c-69667e4497b8",
   "metadata": {},
   "outputs": [],
   "source": [
    "fontsize = 20\n",
    "\n",
    "plt.rc('font', size=fontsize)          # controls default text sizes\n",
    "plt.rc('axes', titlesize=fontsize)     # fontsize of the axes title\n",
    "plt.rc('axes', labelsize=fontsize)    # fontsize of the x and y labels\n",
    "plt.rc('xtick', labelsize=fontsize)    # fontsize of the tick labels\n",
    "plt.rc('ytick', labelsize=fontsize)    # fontsize of the tick labels\n",
    "plt.rc('legend', fontsize=fontsize)    # legend fontsize\n",
    "plt.rc('figure', titlesize=fontsize)  # fontsize of the figure title"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "898ead40-d445-4fb7-a25e-3598c00a29bf",
   "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",
       "html[data-theme=dark],\n",
       "body[data-theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block !important;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-index-preview {\n",
       "  grid-column: 2 / 5;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data,\n",
       ".xr-index-data-in:checked ~ .xr-index-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-index-name div,\n",
       ".xr-index-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data,\n",
       ".xr-index-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2,\n",
       ".xr-no-icon {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.Dataset&gt; Size: 202kB\n",
       "Dimensions:                        (obs: 223, NbSample: 20)\n",
       "Dimensions without coordinates: obs, NbSample\n",
       "Data variables: (12/27)\n",
       "    amplitude                      (obs) float64 2kB ...\n",
       "    effective_area                 (obs) float32 892B ...\n",
       "    effective_contour_height       (obs) float32 892B ...\n",
       "    effective_contour_latitude     (obs, NbSample) float64 36kB ...\n",
       "    effective_contour_longitude    (obs, NbSample) float64 36kB ...\n",
       "    effective_contour_shape_error  (obs) float64 2kB ...\n",
       "    ...                             ...\n",
       "    speed_contour_longitude        (obs, NbSample) float64 36kB ...\n",
       "    speed_contour_shape_error      (obs) float64 2kB ...\n",
       "    speed_radius                   (obs) float64 2kB ...\n",
       "    time                           (obs) datetime64[ns] 2kB ...\n",
       "    track                          (obs) uint32 892B ...\n",
       "    uavg_profile                   (obs, NbSample) float64 36kB ...\n",
       "Attributes: (12/19)\n",
       "    Metadata_Conventions:      Unidata Dataset Discovery v1.0\n",
       "    comment:                   Surface product; mesoscale eddies\n",
       "    creator_email:             aviso@altimetry.fr\n",
       "    creator_url:               https://www.aviso.altimetry.fr\n",
       "    date_created:              2024-08-06T09:38:06Z\n",
       "    framework_used:            https://github.com/AntSimi/py-eddy-tracker\n",
       "    ...                        ...\n",
       "    standard_name_vocabulary:  NetCDF Climate and Forecast (CF) Metadata Conv...\n",
       "    summary:                   This dataset contains eddy atlas from all-sate...\n",
       "    time_coverage_duration:    P2396D\n",
       "    time_coverage_end:         2024-07-23T00:00:00Z\n",
       "    time_coverage_start:       2018-01-01T00:00:00Z\n",
       "    title:                     Mesoscale Cyclonic Eddies in Altimeter Observa...</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-98f5f99f-5a97-4216-ab60-5556438a3626' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-98f5f99f-5a97-4216-ab60-5556438a3626' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span>obs</span>: 223</li><li><span>NbSample</span>: 20</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-694ee8b0-a88f-4cc5-8e8b-ae2bb68ce765' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-694ee8b0-a88f-4cc5-8e8b-ae2bb68ce765' class='xr-section-summary'  title='Expand/collapse section'>Coordinates: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-5ef1e2bd-42e4-4f8d-9997-2059ad25aafc' class='xr-section-summary-in' type='checkbox'  ><label for='section-5ef1e2bd-42e4-4f8d-9997-2059ad25aafc' class='xr-section-summary' >Data variables: <span>(27)</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>amplitude</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3f43fe9f-439a-4557-8022-bd6c3917ba19' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3f43fe9f-439a-4557-8022-bd6c3917ba19' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-cbe66374-c1b9-4567-85e4-153922ce0921' class='xr-var-data-in' type='checkbox'><label for='data-cbe66374-c1b9-4567-85e4-153922ce0921' 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>comment :</span></dt><dd>Magnitude of the height difference between the extremum of SSH within the eddy and the SSH around the effective contour defining the eddy edge</dd><dt><span>long_name :</span></dt><dd>Amplitude</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>min :</span></dt><dd>0.004</dd><dt><span>max :</span></dt><dd>1.405</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>effective_area</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e74c7da1-d6f2-4705-b997-09dea8c9dc5e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e74c7da1-d6f2-4705-b997-09dea8c9dc5e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0dd80a89-209e-4fe4-8ac6-7d9592e8c37a' class='xr-var-data-in' type='checkbox'><label for='data-0dd80a89-209e-4fe4-8ac6-7d9592e8c37a' 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>comment :</span></dt><dd>Area enclosed by the effective contour in m^2</dd><dt><span>long_name :</span></dt><dd>Effective area</dd><dt><span>units :</span></dt><dd>m^2</dd><dt><span>min :</span></dt><dd>128509144.0</dd><dt><span>max :</span></dt><dd>325393650000.0</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>effective_contour_height</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-07cc3e66-7d89-41b0-8bfe-9cb009b9cc4c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-07cc3e66-7d89-41b0-8bfe-9cb009b9cc4c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8e3a4bd7-bf46-4393-856c-e720dadcbf07' class='xr-var-data-in' type='checkbox'><label for='data-8e3a4bd7-bf46-4393-856c-e720dadcbf07' 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>comment :</span></dt><dd>SSH filtered height for effective contour</dd><dt><span>long_name :</span></dt><dd>Effective Contour Height</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>min :</span></dt><dd>-0.918</dd><dt><span>max :</span></dt><dd>0.446</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>effective_contour_latitude</span></div><div class='xr-var-dims'>(obs, NbSample)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3bb2f457-4682-4231-81fc-279f788dae4a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3bb2f457-4682-4231-81fc-279f788dae4a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-50c8c1c2-a56e-4e45-92b0-acfb1be44c9f' class='xr-var-data-in' type='checkbox'><label for='data-50c8c1c2-a56e-4e45-92b0-acfb1be44c9f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>X</dd><dt><span>comment :</span></dt><dd>Latitudes of effective contour</dd><dt><span>long_name :</span></dt><dd>Effective Contour Latitudes</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>[4460 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>effective_contour_longitude</span></div><div class='xr-var-dims'>(obs, NbSample)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-3dbde689-4c4f-42bb-a4a1-b8b588b26f86' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3dbde689-4c4f-42bb-a4a1-b8b588b26f86' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3290aeba-db39-46f1-a0bc-f7dfea88491a' class='xr-var-data-in' type='checkbox'><label for='data-3290aeba-db39-46f1-a0bc-f7dfea88491a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>X</dd><dt><span>comment :</span></dt><dd>Longitudes of the effective contour</dd><dt><span>long_name :</span></dt><dd>Effective Contour Longitudes</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>[4460 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>effective_contour_shape_error</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-f447d5ce-2bef-4ff2-83c2-b83804308e75' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f447d5ce-2bef-4ff2-83c2-b83804308e75' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-458b63d9-704b-4e07-a5aa-27b2d8e4048f' class='xr-var-data-in' type='checkbox'><label for='data-458b63d9-704b-4e07-a5aa-27b2d8e4048f' 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>comment :</span></dt><dd>Error criterion between the effective contour and its best fit circle</dd><dt><span>long_name :</span></dt><dd>Effective Contour Shape Error</dd><dt><span>units :</span></dt><dd>%</dd><dt><span>min :</span></dt><dd>2.5</dd><dt><span>max :</span></dt><dd>75.5</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>effective_radius</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-fbb62839-0892-4a07-9c9c-e5c9b4176409' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-fbb62839-0892-4a07-9c9c-e5c9b4176409' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-5ffe2cb5-3a54-43ef-9343-0b6d91c81c72' class='xr-var-data-in' type='checkbox'><label for='data-5ffe2cb5-3a54-43ef-9343-0b6d91c81c72' 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>comment :</span></dt><dd>Radius of the best fit circle corresponding to the effective contour</dd><dt><span>long_name :</span></dt><dd>Effective Radius</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>min :</span></dt><dd>7050.0</dd><dt><span>max :</span></dt><dd>378400.0</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>inner_contour_height</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-577dae24-6701-424e-9d17-195663b74b20' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-577dae24-6701-424e-9d17-195663b74b20' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-db46127c-8fdd-4610-bccd-e9013c89af96' class='xr-var-data-in' type='checkbox'><label for='data-db46127c-8fdd-4610-bccd-e9013c89af96' 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>comment :</span></dt><dd>SSH filtered height for the smallest detected contour</dd><dt><span>long_name :</span></dt><dd>Inner Contour Height</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>min :</span></dt><dd>-1.462</dd><dt><span>max :</span></dt><dd>0.406</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>latitude</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4d97b07a-0d23-4a22-834f-aea8d0f81c4c' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4d97b07a-0d23-4a22-834f-aea8d0f81c4c' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-e1b77bff-5254-40d1-97a1-5b5669ffe9e3' class='xr-var-data-in' type='checkbox'><label for='data-e1b77bff-5254-40d1-97a1-5b5669ffe9e3' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>Y</dd><dt><span>comment :</span></dt><dd>Latitude center of the best fit circle</dd><dt><span>long_name :</span></dt><dd>Eddy Center Latitude</dd><dt><span>standard_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd><dt><span>min :</span></dt><dd>-78.05134</dd><dt><span>max :</span></dt><dd>84.03748</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>latitude_max</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-9ea71f0c-a9d9-44fa-830f-5c3799c81803' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9ea71f0c-a9d9-44fa-830f-5c3799c81803' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c8c9e7c3-974c-4f00-a2fc-a83200fb41b2' class='xr-var-data-in' type='checkbox'><label for='data-c8c9e7c3-974c-4f00-a2fc-a83200fb41b2' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>Y</dd><dt><span>comment :</span></dt><dd>Latitude of the inner contour</dd><dt><span>long_name :</span></dt><dd>Latitude of the SSH maximum</dd><dt><span>standard_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd><dt><span>min :</span></dt><dd>-90.0</dd><dt><span>max :</span></dt><dd>84.12464</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>longitude</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-39711bbc-71e6-4b0a-a86f-961bc9bffe93' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-39711bbc-71e6-4b0a-a86f-961bc9bffe93' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-902f9edf-f75a-4e10-9523-630ed6134ef0' class='xr-var-data-in' type='checkbox'><label for='data-902f9edf-f75a-4e10-9523-630ed6134ef0' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>X</dd><dt><span>comment :</span></dt><dd>Longitude center of the best fit circle</dd><dt><span>long_name :</span></dt><dd>Eddy Center Longitude</dd><dt><span>standard_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd><dt><span>min :</span></dt><dd>-28.260529</dd><dt><span>max :</span></dt><dd>367.2987</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>longitude_max</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d4523588-48ff-4a2c-b218-09c495d90754' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d4523588-48ff-4a2c-b218-09c495d90754' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-01691d0b-7721-44ec-8e5b-955709c4ff51' class='xr-var-data-in' type='checkbox'><label for='data-01691d0b-7721-44ec-8e5b-955709c4ff51' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>X</dd><dt><span>comment :</span></dt><dd>Longitude of the inner contour</dd><dt><span>long_name :</span></dt><dd>Longitude of the SSH maximum</dd><dt><span>standard_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd><dt><span>min :</span></dt><dd>-28.13443</dd><dt><span>max :</span></dt><dd>367.44983</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>num_contours</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e5ef0792-ee9b-4cd9-b18e-15f23cc92f34' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e5ef0792-ee9b-4cd9-b18e-15f23cc92f34' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0e99570c-f300-441a-9232-4d1cb67b519a' class='xr-var-data-in' type='checkbox'><label for='data-0e99570c-f300-441a-9232-4d1cb67b519a' 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>comment :</span></dt><dd>Number of contours selected for this eddy</dd><dt><span>long_name :</span></dt><dd>Number of contours</dd><dt><span>min :</span></dt><dd>1</dd><dt><span>max :</span></dt><dd>703</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>num_point_e</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-71d774c3-501c-492f-8341-abcfff3936ee' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-71d774c3-501c-492f-8341-abcfff3936ee' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-16280442-2969-409c-b6a5-eeb25e545c17' class='xr-var-data-in' type='checkbox'><label for='data-16280442-2969-409c-b6a5-eeb25e545c17' 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>description :</span></dt><dd>Number of points for effective contour before resampling</dd><dt><span>long_name :</span></dt><dd>number of points for effective contour</dd><dt><span>units :</span></dt><dd>ordinal</dd><dt><span>min :</span></dt><dd>10</dd><dt><span>max :</span></dt><dd>145</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>num_point_s</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-44aaed97-52ee-420e-b7e1-f27efdffa28f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-44aaed97-52ee-420e-b7e1-f27efdffa28f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d3555c60-3bc0-4ef5-9045-7c63e2ce30b1' class='xr-var-data-in' type='checkbox'><label for='data-d3555c60-3bc0-4ef5-9045-7c63e2ce30b1' 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>description :</span></dt><dd>Number of points for speed contour before resampling</dd><dt><span>long_name :</span></dt><dd>number of points for speed contour</dd><dt><span>units :</span></dt><dd>ordinal</dd><dt><span>min :</span></dt><dd>10</dd><dt><span>max :</span></dt><dd>135</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>observation_flag</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>int8</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-889488e9-47b4-4fcc-b5f9-3145ce5d9507' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-889488e9-47b4-4fcc-b5f9-3145ce5d9507' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fa95b9ee-e138-476b-9a5b-df1d1fd12d3d' class='xr-var-data-in' type='checkbox'><label for='data-fa95b9ee-e138-476b-9a5b-df1d1fd12d3d' 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>comment :</span></dt><dd>Flag indicating if the value is interpolated between two observations or not (0: observed eddy, 1: interpolated eddy)</dd><dt><span>long_name :</span></dt><dd>Virtual Eddy Position</dd><dt><span>min :</span></dt><dd>0</dd><dt><span>max :</span></dt><dd>1</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=int8]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>observation_number</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>uint16</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-171819c3-5e8b-4f59-b5c7-d5b72230d6c0' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-171819c3-5e8b-4f59-b5c7-d5b72230d6c0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8074171a-ed34-4862-ac0d-5ca72a6e1ece' class='xr-var-data-in' type='checkbox'><label for='data-8074171a-ed34-4862-ac0d-5ca72a6e1ece' 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>comment :</span></dt><dd>Observation sequence number, days starting at the eddy first detection</dd><dt><span>long_name :</span></dt><dd>Eddy temporal index in a trajectory</dd><dt><span>min :</span></dt><dd>0</dd><dt><span>max :</span></dt><dd>1961</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=uint16]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_area</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-32abfd96-3f12-40ab-8285-af7a1837a0ec' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-32abfd96-3f12-40ab-8285-af7a1837a0ec' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-2ab8c851-cc78-4769-bd61-c2b17d2597d7' class='xr-var-data-in' type='checkbox'><label for='data-2ab8c851-cc78-4769-bd61-c2b17d2597d7' 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>comment :</span></dt><dd>Area enclosed by the speed contour in m^2</dd><dt><span>long_name :</span></dt><dd>Speed area</dd><dt><span>units :</span></dt><dd>m^2</dd><dt><span>min :</span></dt><dd>128509144.0</dd><dt><span>max :</span></dt><dd>258590110000.0</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_average</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-ec37e62f-acfa-4e27-ab92-5f62259104ad' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-ec37e62f-acfa-4e27-ab92-5f62259104ad' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0de33144-0ac4-4472-bbef-dfc2917fc0c0' class='xr-var-data-in' type='checkbox'><label for='data-0de33144-0ac4-4472-bbef-dfc2917fc0c0' 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>comment :</span></dt><dd>Average speed of the contour defining the radius scale “speed_radius”</dd><dt><span>long_name :</span></dt><dd>Maximum circum-averaged Speed</dd><dt><span>units :</span></dt><dd>m/s</dd><dt><span>min :</span></dt><dd>0.0</dd><dt><span>max :</span></dt><dd>5.6524</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_contour_height</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-090f23e9-45bf-4423-8268-f61ad41816a6' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-090f23e9-45bf-4423-8268-f61ad41816a6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-8e5e2a2a-a7fb-4d28-a3fc-0a4931e1f20b' class='xr-var-data-in' type='checkbox'><label for='data-8e5e2a2a-a7fb-4d28-a3fc-0a4931e1f20b' 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>comment :</span></dt><dd>SSH filtered height for speed contour</dd><dt><span>long_name :</span></dt><dd>Speed Contour Height</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>min :</span></dt><dd>-0.918</dd><dt><span>max :</span></dt><dd>0.446</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_contour_latitude</span></div><div class='xr-var-dims'>(obs, NbSample)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-9636e434-70f2-47aa-a421-7c9ce4e7eb89' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-9636e434-70f2-47aa-a421-7c9ce4e7eb89' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-fd6aaf7e-4592-42e4-a223-e566faf9e174' class='xr-var-data-in' type='checkbox'><label for='data-fd6aaf7e-4592-42e4-a223-e566faf9e174' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>X</dd><dt><span>comment :</span></dt><dd>Latitudes of speed contour</dd><dt><span>long_name :</span></dt><dd>Speed Contour Latitudes</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>[4460 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_contour_longitude</span></div><div class='xr-var-dims'>(obs, NbSample)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-2c829909-7b54-4f7a-8daf-c4a71564a00f' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-2c829909-7b54-4f7a-8daf-c4a71564a00f' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b8af32a8-56e8-4e41-950a-2e3661395aa5' class='xr-var-data-in' type='checkbox'><label for='data-b8af32a8-56e8-4e41-950a-2e3661395aa5' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>X</dd><dt><span>comment :</span></dt><dd>Longitudes of speed contour</dd><dt><span>long_name :</span></dt><dd>Speed Contour Longitudes</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>[4460 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_contour_shape_error</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-e82ba1f0-0854-4954-8cf6-f1e518772451' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-e82ba1f0-0854-4954-8cf6-f1e518772451' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d007054b-849b-43ac-bdad-f46188e99cd9' class='xr-var-data-in' type='checkbox'><label for='data-d007054b-849b-43ac-bdad-f46188e99cd9' 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>comment :</span></dt><dd>Error criterion between the speed contour and its best fit circle</dd><dt><span>long_name :</span></dt><dd>Speed Contour Shape Error</dd><dt><span>units :</span></dt><dd>%</dd><dt><span>min :</span></dt><dd>1.5</dd><dt><span>max :</span></dt><dd>93.0</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>speed_radius</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-45fbf992-1bcd-4bdf-9090-34fbb3c73edd' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-45fbf992-1bcd-4bdf-9090-34fbb3c73edd' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-aa7a7181-706c-4f79-8fda-66da821d6537' class='xr-var-data-in' type='checkbox'><label for='data-aa7a7181-706c-4f79-8fda-66da821d6537' 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>comment :</span></dt><dd>Radius of the best fit circle corresponding to the contour of maximum circum-average speed</dd><dt><span>long_name :</span></dt><dd>Speed Radius</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>min :</span></dt><dd>7050.0</dd><dt><span>max :</span></dt><dd>342200.0</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>time</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4edd0fe4-847d-4a0e-af14-56f643c3d421' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4edd0fe4-847d-4a0e-af14-56f643c3d421' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-47bbe1be-6032-4b7c-8349-17c797636ff4' class='xr-var-data-in' type='checkbox'><label for='data-47bbe1be-6032-4b7c-8349-17c797636ff4' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>axis :</span></dt><dd>T</dd><dt><span>comment :</span></dt><dd>Date of this observation</dd><dt><span>long_name :</span></dt><dd>Time</dd><dt><span>standard_name :</span></dt><dd>time</dd><dt><span>min :</span></dt><dd>24837.0</dd><dt><span>max :</span></dt><dd>27232.0</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=datetime64[ns]]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>track</span></div><div class='xr-var-dims'>(obs)</div><div class='xr-var-dtype'>uint32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-02d1c6b4-e97b-45f6-916a-f84baf423f19' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-02d1c6b4-e97b-45f6-916a-f84baf423f19' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-9bbde46f-811b-4e7e-9e50-5a5596840332' class='xr-var-data-in' type='checkbox'><label for='data-9bbde46f-811b-4e7e-9e50-5a5596840332' 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>comment :</span></dt><dd>Trajectory identification number</dd><dt><span>long_name :</span></dt><dd>Trajectory number</dd><dt><span>min :</span></dt><dd>0</dd><dt><span>max :</span></dt><dd>191385</dd></dl></div><div class='xr-var-data'><pre>[223 values with dtype=uint32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>uavg_profile</span></div><div class='xr-var-dims'>(obs, NbSample)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-2ccfd07e-d0b6-4a94-8b4d-639eaf7f8e63' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-2ccfd07e-d0b6-4a94-8b4d-639eaf7f8e63' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-83d1361a-f0bc-4736-b5a5-910c44c8c271' class='xr-var-data-in' type='checkbox'><label for='data-83d1361a-f0bc-4736-b5a5-910c44c8c271' 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>comment :</span></dt><dd>Speed averaged values from the effective contour inwards to the smallest contour, evenly spaced points</dd><dt><span>long_name :</span></dt><dd>Radial Speed Profile</dd><dt><span>units :</span></dt><dd>m/s</dd></dl></div><div class='xr-var-data'><pre>[4460 values with dtype=float64]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-0bb98524-f968-4972-838f-297e5e6f1f91' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-0bb98524-f968-4972-838f-297e5e6f1f91' class='xr-section-summary'  title='Expand/collapse section'>Indexes: <span>(0)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'></ul></div></li><li class='xr-section-item'><input id='section-3140006e-b24f-4371-97c3-a8a797d57caa' class='xr-section-summary-in' type='checkbox'  ><label for='section-3140006e-b24f-4371-97c3-a8a797d57caa' class='xr-section-summary' >Attributes: <span>(19)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>Metadata_Conventions :</span></dt><dd>Unidata Dataset Discovery v1.0</dd><dt><span>comment :</span></dt><dd>Surface product; mesoscale eddies</dd><dt><span>creator_email :</span></dt><dd>aviso@altimetry.fr</dd><dt><span>creator_url :</span></dt><dd>https://www.aviso.altimetry.fr</dd><dt><span>date_created :</span></dt><dd>2024-08-06T09:38:06Z</dd><dt><span>framework_used :</span></dt><dd>https://github.com/AntSimi/py-eddy-tracker</dd><dt><span>framework_version :</span></dt><dd>v3.6.1</dd><dt><span>input_maps_processing :</span></dt><dd>DUACS DT2018 until 2021-12-07, DUACS DT2021 since 2021-12-08</dd><dt><span>input_product_reference :</span></dt><dd>https://www.aviso.altimetry.fr/en/data/products/value-added-products/global-mesoscale-eddy-trajectory-product.html</dd><dt><span>institution :</span></dt><dd>CLS/CNES in collaboration with IMEDEA</dd><dt><span>license :</span></dt><dd>https://www.aviso.altimetry.fr/fileadmin/documents/data/License_Aviso.pdf</dd><dt><span>product_version :</span></dt><dd>3.2exp</dd><dt><span>project :</span></dt><dd>SSALTO/DUACS</dd><dt><span>standard_name_vocabulary :</span></dt><dd>NetCDF Climate and Forecast (CF) Metadata Convention Standard Name Table</dd><dt><span>summary :</span></dt><dd>This dataset contains eddy atlas from all-satellite NRT CMEMS maps over global ocean (SEALEVEL_GLO_PHY_L4_NRT_OBSERVATIONS*)</dd><dt><span>time_coverage_duration :</span></dt><dd>P2396D</dd><dt><span>time_coverage_end :</span></dt><dd>2024-07-23T00:00:00Z</dd><dt><span>time_coverage_start :</span></dt><dd>2018-01-01T00:00:00Z</dd><dt><span>title :</span></dt><dd>Mesoscale Cyclonic Eddies in Altimeter Observations of ADT</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset> Size: 202kB\n",
       "Dimensions:                        (obs: 223, NbSample: 20)\n",
       "Dimensions without coordinates: obs, NbSample\n",
       "Data variables: (12/27)\n",
       "    amplitude                      (obs) float64 2kB ...\n",
       "    effective_area                 (obs) float32 892B ...\n",
       "    effective_contour_height       (obs) float32 892B ...\n",
       "    effective_contour_latitude     (obs, NbSample) float64 36kB ...\n",
       "    effective_contour_longitude    (obs, NbSample) float64 36kB ...\n",
       "    effective_contour_shape_error  (obs) float64 2kB ...\n",
       "    ...                             ...\n",
       "    speed_contour_longitude        (obs, NbSample) float64 36kB ...\n",
       "    speed_contour_shape_error      (obs) float64 2kB ...\n",
       "    speed_radius                   (obs) float64 2kB ...\n",
       "    time                           (obs) datetime64[ns] 2kB ...\n",
       "    track                          (obs) uint32 892B ...\n",
       "    uavg_profile                   (obs, NbSample) float64 36kB ...\n",
       "Attributes: (12/19)\n",
       "    Metadata_Conventions:      Unidata Dataset Discovery v1.0\n",
       "    comment:                   Surface product; mesoscale eddies\n",
       "    creator_email:             aviso@altimetry.fr\n",
       "    creator_url:               https://www.aviso.altimetry.fr\n",
       "    date_created:              2024-08-06T09:38:06Z\n",
       "    framework_used:            https://github.com/AntSimi/py-eddy-tracker\n",
       "    ...                        ...\n",
       "    standard_name_vocabulary:  NetCDF Climate and Forecast (CF) Metadata Conv...\n",
       "    summary:                   This dataset contains eddy atlas from all-sate...\n",
       "    time_coverage_duration:    P2396D\n",
       "    time_coverage_end:         2024-07-23T00:00:00Z\n",
       "    time_coverage_start:       2018-01-01T00:00:00Z\n",
       "    title:                     Mesoscale Cyclonic Eddies in Altimeter Observa..."
      ]
     },
     "execution_count": 4,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "ds = xr.open_dataset('Edward_Eddy_trajectory_nrt_3.2exp_cyclonic_20180101_20240723.nc')\n",
    "ds"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "d5ca866e-0650-4eaf-889d-5f583306aaef",
   "metadata": {},
   "outputs": [],
   "source": [
    "def in_eddy(ds,float_lat,float_lon,float_time):\n",
    "    \"\"\"\n",
    "    float_lat: degrees north\n",
    "    float_lon: degrees east\n",
    "    float_time: should be in format 'YYYY-MM-DD'\n",
    "    \"\"\"\n",
    "\n",
    "    def all_equal(iterable):\n",
    "        g = groupby(iterable)\n",
    "        return next(g, True) and not next(g, False)\n",
    "\n",
    "    in_eddy_flag = False \n",
    "\n",
    "    float_time = np.datetime64(float_time)\n",
    "    if float_time in ds.time: # some dates not in dateset        \n",
    "        for i in np.where(ds.time == float_time)[0]: # usually will only be 1 eddy, but sometimes there are 2 after a split\n",
    "            contour_lons = np.array(ds.effective_contour_longitude[i]) # eddy lons\n",
    "            contour_lats = np.array(ds.effective_contour_latitude[i]) # eddy lats\n",
    "            \n",
    "            if all_equal(contour_lons): # eddy break\n",
    "                pass\n",
    "            else:\n",
    "                poly = Path([(contour_lats[j],contour_lons[j]) for j in np.arange(0,len(contour_lats))]) # set up the polygon\n",
    "                if poly.contains_points([(float_lat,float_lon)]): #find if point is inside the polygon\n",
    "                    in_eddy_flag = True \n",
    "\n",
    "    return in_eddy_flag"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "fc64a12c-e504-4cdd-a475-12e8acec6c16",
   "metadata": {},
   "outputs": [],
   "source": [
    "# in_eddy(ds,35,289.5,'2024-05-08')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "a59ddcf3-60c9-4b10-81c4-2ce3045888b3",
   "metadata": {},
   "outputs": [],
   "source": [
    "satdate,N,chl_mean,chl_std,poc_mean,poc_std=[],[],[],[],[],[];"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "ddcd141d-a1b6-4c92-8d82-eb49e37ac174",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "7843e1695b4746228c76720876a7373a",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "QUEUEING TASKS | :   0%|          | 0/273 [00:00<?, ?it/s]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "208a76d2d9d9451fa5676e3eb58df156",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "PROCESSING TASKS | :   0%|          | 0/273 [00:00<?, ?it/s]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    },
    {
     "data": {
      "application/vnd.jupyter.widget-view+json": {
       "model_id": "4f8f4bd78c724b0fb8d8aa3c991478e7",
       "version_major": 2,
       "version_minor": 0
      },
      "text/plain": [
       "COLLECTING RESULTS | :   0%|          | 0/273 [00:00<?, ?it/s]"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# PACE CHL\n",
    "results = earthaccess.search_datasets(instrument=\"oci\")\n",
    "\n",
    "tspan = (\"2024-04-19\", \"2024-10-18\")\n",
    "bbox = (-73, 34, -70, 37.5) #edward\n",
    "clouds = (0, 100)\n",
    "results = earthaccess.search_data(\n",
    "    short_name=\"PACE_OCI_L2_BGC_NRT\",\n",
    "    temporal=tspan,\n",
    "    bounding_box=bbox,\n",
    "    cloud_cover=clouds,\n",
    ")\n",
    "\n",
    "paths = earthaccess.open(results) #'streaming' data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "ac2ca9d4-efaa-4361-b183-f66d1d5b4304",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "273"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "len(paths)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "d8c367a7-d8f2-4746-a768-cb1473540a78",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Mask via latitude only\n",
      "<xarray.DataArray 'longitude' ()> Size: 8B\n",
      "array(230957)\n"
     ]
    }
   ],
   "source": [
    "for IDX in range(261,len(paths)): #FOR EACH IMAGE. had to do in chunks \n",
    "#for IDX in range(0,1): #FOR EACH IMAGE\n",
    "    datatree = open_datatree(paths[IDX])\n",
    "    dataset = xr.merge(datatree.to_dict().values())\n",
    "    dataset = dataset.set_coords((\"longitude\", \"latitude\"))\n",
    "    h=str(paths[IDX])\n",
    "    sat_date=h[64:72]\n",
    "    mask = (dataset.longitude >= bbox[0]) & (dataset.longitude <= bbox[2]) & (dataset.latitude >= bbox[1]) & (dataset.latitude <= bbox[3])\n",
    "    if mask.sum() == 0:\n",
    "        mask = (dataset.longitude >= bbox[0]) & (dataset.longitude <= bbox[2])\n",
    "        if mask.sum() == 0:\n",
    "            print('Mask via longitude only')\n",
    "            mask = (dataset.latitude >= bbox[1]) & (dataset.latitude <= bbox[3])\n",
    "            print(mask.sum())\n",
    "        else:\n",
    "            print('Mask via latitude only')\n",
    "            print(mask.sum())\n",
    "    ds_masked = dataset.where(mask, drop=True)\n",
    "    lat=ds_masked.latitude.values.flatten()\n",
    "    lon=ds_masked.longitude.values.flatten()+360\n",
    "    chl=ds_masked.chlor_a.values.flatten()\n",
    "    poc=ds_masked.poc.values.flatten()\n",
    "    \n",
    "    chl_inside,poc_inside=[],[]\n",
    "    \n",
    "    for index in np.arange(0,len(lat)):\n",
    "        val=in_eddy(ds,lat[index],lon[index],'2024-05-08')\n",
    "        if val:\n",
    "            chl_inside.append(chl[index]) \n",
    "            poc_inside.append(poc[index])\n",
    "    \n",
    "    N1=np.size(np.where(np.isnan(chl_inside)==False))\n",
    "\n",
    "    if N1>5000:\n",
    "        chlmean=np.nanmean(chl_inside);\n",
    "        chlstd=np.nanstd(chl_inside);\n",
    "        pocmean=np.nanmean(poc_inside);\n",
    "        pocstd=np.nanstd(poc_inside);\n",
    "        \n",
    "        chl_mean.append(chlmean);\n",
    "        poc_mean.append(pocmean);\n",
    "        chl_std.append(chlstd);\n",
    "        poc_std.append(pocstd);\n",
    "        N.append(N1);\n",
    "        satdate.append(sat_date);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3a756b30-5677-4450-90f8-712cbf255129",
   "metadata": {},
   "outputs": [],
   "source": [
    "d = {'satdate': satdate, 'N':N, 'chlmean': chl_mean, 'pocmean': poc_mean, 'chlstd': chl_std, 'pocstd': poc_std}\n",
    "data = pandas.DataFrame(data=d)\n",
    "data"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "871dacbe-d4e7-47dc-bd46-569bed17841e",
   "metadata": {},
   "outputs": [],
   "source": [
    "data.to_csv('CHL_TIMESERIES_PACE.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "6ac4dfc1-7b45-40b6-bb47-f2f04cfcedb3",
   "metadata": {},
   "outputs": [],
   "source": [
    "IDX"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.14"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}