{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Nutrient Evaluation\n", "\n", "Compare nutrient fields to observations from the World Ocean Atlas. Fields included here are:\n", "* NO3\n", "* PO4\n", "* SiO3" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Imports" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Author: Mike Levy\n", "\n", "Compiler : GCC 9.3.0\n", "OS : Linux\n", "Release : 3.10.0-1127.18.2.el7.x86_64\n", "Machine : x86_64\n", "Processor : x86_64\n", "CPU cores : 72\n", "Architecture: 64bit\n", "\n", "Hostname: crhtc50\n", "\n", "Git hash: 135f85d57e8712d47b21027021f8409ad12c469c\n", "\n", "matplotlib: 3.4.2\n", "cmocean : 2.0\n", "cartopy : 0.19.0.post1\n", "xarray : 0.18.2\n", "numpy : 1.20.3\n", "\n" ] } ], "source": [ "%matplotlib inline\n", "import os\n", "\n", "import numpy as np\n", "import xarray as xr\n", "\n", "import utils\n", "\n", "import matplotlib.pyplot as plt\n", "import matplotlib.gridspec as gridspec\n", "import matplotlib.colors as colors\n", "\n", "import cartopy\n", "import cartopy.crs as ccrs\n", "import cmocean\n", "\n", "%load_ext watermark\n", "%watermark -a \"Mike Levy\" -d -iv -m -g -h" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Read in the Data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "nutrients = ['NO3', 'PO4', 'SiO3']" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Read in the surface data" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n", "Dimensions: (nlat: 384, nlon: 321)\n", "Dimensions without coordinates: nlat, nlon\n", "Data variables:\n", " NO3 (nlat, nlon) float32 dask.array<chunksize=(192, 321), meta=np.ndarray>\n", " NO3_obs (nlat, nlon) float32 dask.array<chunksize=(192, 321), meta=np.ndarray>\n", " PO4 (nlat, nlon) float32 dask.array<chunksize=(192, 321), meta=np.ndarray>\n", " PO4_obs (nlat, nlon) float32 dask.array<chunksize=(192, 321), meta=np.ndarray>\n", " SiO3 (nlat, nlon) float32 dask.array<chunksize=(192, 321), meta=np.ndarray>\n", " SiO3_obs (nlat, nlon) float32 dask.array<chunksize=(192, 321), meta=np.ndarray>\n", " TLAT (nlat, nlon) float64 dask.array<chunksize=(192, 161), meta=np.ndarray>\n", " TLONG (nlat, nlon) float64 dask.array<chunksize=(192, 161), meta=np.ndarray>
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
<xarray.Dataset>\n", "Dimensions: (basins: 4, lat_t: 394, z_t: 60)\n", "Coordinates:\n", " * lat_t (lat_t) float32 -79.22 -78.69 -78.15 -77.62 ... 88.84 89.37 89.9\n", " * basins (basins) <U14 'Global' 'Atlantic Ocean' ... 'Indian Ocean'\n", " * z_t (z_t) float64 5.0 15.0 25.0 ... 4.875e+03 5.125e+03 5.375e+03\n", " z_t_ (z_t) float64 500.0 1.5e+03 2.5e+03 ... 5.125e+05 5.375e+05\n", "Data variables:\n", " NO3 (basins, z_t, lat_t) float32 nan nan 19.02 19.07 ... nan nan nan\n", " NO3_obs (basins, z_t, lat_t) float32 nan nan 20.99 21.69 ... nan nan nan\n", " PO4 (basins, z_t, lat_t) float32 nan nan 1.423 1.428 ... nan nan nan\n", " PO4_obs (basins, z_t, lat_t) float32 nan nan 1.713 1.705 ... nan nan nan\n", " SiO3 (basins, z_t, lat_t) float32 nan nan 42.66 43.43 ... nan nan nan\n", " SiO3_obs (basins, z_t, lat_t) float32 nan nan 61.78 62.55 ... nan nan nan\n", " NO3_bias (basins, z_t, lat_t) float32 nan nan -1.967 ... nan nan nan\n", " PO4_bias (basins, z_t, lat_t) float32 nan nan -0.2902 ... nan nan nan\n", " SiO3_bias (basins, z_t, lat_t) float32 nan nan -19.12 ... nan nan nan
array([-79.22052 , -78.68631 , -78.15209 , ..., 88.836334, 89.370575,\n", " 89.904816], dtype=float32)
array(['Global', 'Atlantic Ocean', 'Pacific Ocean', 'Indian Ocean'],\n", " dtype='<U14')
array([5.000000e+00, 1.500000e+01, 2.500000e+01, 3.500000e+01, 4.500000e+01,\n", " 5.500000e+01, 6.500000e+01, 7.500000e+01, 8.500000e+01, 9.500000e+01,\n", " 1.050000e+02, 1.150000e+02, 1.250000e+02, 1.350000e+02, 1.450000e+02,\n", " 1.550000e+02, 1.650984e+02, 1.754791e+02, 1.862913e+02, 1.976603e+02,\n", " 2.097114e+02, 2.225783e+02, 2.364088e+02, 2.513702e+02, 2.676542e+02,\n", " 2.854837e+02, 3.051192e+02, 3.268680e+02, 3.510935e+02, 3.782276e+02,\n", " 4.087847e+02, 4.433777e+02, 4.827367e+02, 5.277280e+02, 5.793729e+02,\n", " 6.388626e+02, 7.075633e+02, 7.870025e+02, 8.788252e+02, 9.847059e+02,\n", " 1.106204e+03, 1.244567e+03, 1.400497e+03, 1.573946e+03, 1.764003e+03,\n", " 1.968944e+03, 2.186457e+03, 2.413972e+03, 2.649001e+03, 2.889385e+03,\n", " 3.133405e+03, 3.379794e+03, 3.627670e+03, 3.876452e+03, 4.125768e+03,\n", " 4.375393e+03, 4.625190e+03, 4.875083e+03, 5.125028e+03, 5.375000e+03])
array([5.00000000e+02, 1.50000000e+03, 2.50000000e+03, 3.50000000e+03,\n", " 4.50000000e+03, 5.50000000e+03, 6.50000000e+03, 7.50000000e+03,\n", " 8.50000000e+03, 9.50000000e+03, 1.05000000e+04, 1.15000000e+04,\n", " 1.25000000e+04, 1.35000000e+04, 1.45000000e+04, 1.55000000e+04,\n", " 1.65098398e+04, 1.75479043e+04, 1.86291270e+04, 1.97660273e+04,\n", " 2.09711387e+04, 2.22578281e+04, 2.36408828e+04, 2.51370156e+04,\n", " 2.67654199e+04, 2.85483652e+04, 3.05119219e+04, 3.26867988e+04,\n", " 3.51093477e+04, 3.78227617e+04, 4.08784648e+04, 4.43377695e+04,\n", " 4.82736719e+04, 5.27728008e+04, 5.79372891e+04, 6.38862617e+04,\n", " 7.07563281e+04, 7.87002500e+04, 8.78825234e+04, 9.84705859e+04,\n", " 1.10620422e+05, 1.24456688e+05, 1.40049719e+05, 1.57394641e+05,\n", " 1.76400328e+05, 1.96894422e+05, 2.18645656e+05, 2.41397156e+05,\n", " 2.64900125e+05, 2.88938469e+05, 3.13340469e+05, 3.37979344e+05,\n", " 3.62767031e+05, 3.87645188e+05, 4.12576812e+05, 4.37539250e+05,\n", " 4.62519031e+05, 4.87508344e+05, 5.12502812e+05, 5.37500000e+05])
array([[[ nan, nan, 19.019262 , ..., 3.127393 ,\n", " 3.1124809, 3.0936654],\n", " [ nan, nan, 19.192722 , ..., 3.2082498,\n", " 3.1962583, 3.17934 ],\n", " [ nan, nan, 19.747913 , ..., 3.527372 ,\n", " 3.5199492, 3.5046253],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, 21.196936 , ..., 3.127393 ,\n", " 3.1124809, 3.0936654],\n", " [ nan, nan, 21.252968 , ..., 3.2082498,\n", " 3.1962583, 3.17934 ],\n", " [ nan, nan, 21.355429 , ..., 3.527372 ,\n", " 3.5199492, 3.5046253],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
array([[[ nan, nan, 20.986591 , ..., 0.98718345,\n", " 0.9528614 , 0.956563 ],\n", " [ nan, nan, 23.002188 , ..., 1.2074015 ,\n", " 1.176309 , 1.1727661 ],\n", " [ nan, nan, 24.706089 , ..., 1.3017465 ,\n", " 1.2707111 , 1.2656828 ],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, 22.391838 , ..., 0.98718345,\n", " 0.9528614 , 0.956563 ],\n", " [ nan, nan, 24.235825 , ..., 1.2074015 ,\n", " 1.176309 , 1.1727661 ],\n", " [ nan, nan, 25.078253 , ..., 1.3017465 ,\n", " 1.2707111 , 1.2656828 ],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
array([[[ nan, nan, 1.4230961, ..., 0.9954035, 0.9968648,\n", " 0.9986534],\n", " [ nan, nan, 1.4344989, ..., 1.0013988, 1.0030755,\n", " 1.005006 ],\n", " [ nan, nan, 1.4741615, ..., 1.0225705, 1.0245651,\n", " 1.0266324],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]],\n", "\n", " [[ nan, nan, 1.5554473, ..., 0.9954035, 0.9968648,\n", " 0.9986534],\n", " [ nan, nan, 1.5587842, ..., 1.0013988, 1.0030755,\n", " 1.005006 ],\n", " [ nan, nan, 1.5648167, ..., 1.0225705, 1.0245651,\n", " 1.0266324],\n", "...\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]],\n", "\n", " [[ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]]], dtype=float32)
array([[[ nan, nan, 1.7132827 , ..., 0.7824321 ,\n", " 0.7851637 , 0.7851058 ],\n", " [ nan, nan, 1.8068694 , ..., 0.812785 ,\n", " 0.81371224, 0.813248 ],\n", " [ nan, nan, 1.8792542 , ..., 0.75705147,\n", " 0.7578078 , 0.7565225 ],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, 1.7572992 , ..., 0.7824321 ,\n", " 0.7851637 , 0.7851058 ],\n", " [ nan, nan, 1.8968039 , ..., 0.812785 ,\n", " 0.81371224, 0.813248 ],\n", " [ nan, nan, 1.9512222 , ..., 0.75705147,\n", " 0.7578078 , 0.7565225 ],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
array([[[ nan, nan, 42.659046 , ..., 7.488578 ,\n", " 7.4984303, 7.523059 ],\n", " [ nan, nan, 42.93306 , ..., 7.5811567,\n", " 7.5947866, 7.621899 ],\n", " [ nan, nan, 43.63968 , ..., 7.896808 ,\n", " 7.916356 , 7.9460306],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, 46.000652 , ..., 7.488578 ,\n", " 7.4984303, 7.523059 ],\n", " [ nan, nan, 46.0551 , ..., 7.5811567,\n", " 7.5947866, 7.621899 ],\n", " [ nan, nan, 46.14845 , ..., 7.896808 ,\n", " 7.916356 , 7.9460306],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
array([[[ nan, nan, 61.779434, ..., 10.102219, 10.192597,\n", " 10.225629],\n", " [ nan, nan, 63.352154, ..., 9.263759, 9.325229,\n", " 9.332197],\n", " [ nan, nan, 64.556366, ..., 9.001168, 9.07389 ,\n", " 9.077654],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]],\n", "\n", " [[ nan, nan, 58.09269 , ..., 10.102219, 10.192597,\n", " 10.225629],\n", " [ nan, nan, 58.541718, ..., 9.263759, 9.325229,\n", " 9.332197],\n", " [ nan, nan, 60.34995 , ..., 9.001168, 9.07389 ,\n", " 9.077654],\n", "...\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]],\n", "\n", " [[ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan],\n", " [ nan, nan, nan, ..., nan, nan,\n", " nan]]], dtype=float32)
array([[[ nan, nan, -1.967329 , ..., 2.1402097,\n", " 2.1596193, 2.1371024],\n", " [ nan, nan, -3.8094654, ..., 2.0008483,\n", " 2.0199494, 2.0065737],\n", " [ nan, nan, -4.9581757, ..., 2.2256255,\n", " 2.249238 , 2.2389426],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, -1.1949024, ..., 2.1402097,\n", " 2.1596193, 2.1371024],\n", " [ nan, nan, -2.9828568, ..., 2.0008483,\n", " 2.0199494, 2.0065737],\n", " [ nan, nan, -3.722824 , ..., 2.2256255,\n", " 2.249238 , 2.2389426],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
array([[[ nan, nan, -0.29018664, ..., 0.21297145,\n", " 0.2117011 , 0.21354759],\n", " [ nan, nan, -0.37237048, ..., 0.18861377,\n", " 0.18936324, 0.19175798],\n", " [ nan, nan, -0.40509272, ..., 0.26551902,\n", " 0.2667573 , 0.27010995],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, -0.20185184, ..., 0.21297145,\n", " 0.2117011 , 0.21354759],\n", " [ nan, nan, -0.3380196 , ..., 0.18861377,\n", " 0.18936324, 0.19175798],\n", " [ nan, nan, -0.38640547, ..., 0.26551902,\n", " 0.2667573 , 0.27010995],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
array([[[ nan, nan, -19.120388 , ..., -2.6136408,\n", " -2.6941671, -2.70257 ],\n", " [ nan, nan, -20.419094 , ..., -1.6826019,\n", " -1.730442 , -1.7102981],\n", " [ nan, nan, -20.916687 , ..., -1.1043601,\n", " -1.1575336, -1.1316233],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, -12.092037 , ..., -2.6136408,\n", " -2.6941671, -2.70257 ],\n", " [ nan, nan, -12.486618 , ..., -1.6826019,\n", " -1.730442 , -1.7102981],\n", " [ nan, nan, -14.2015 , ..., -1.1043601,\n", " -1.1575336, -1.1316233],\n", "...\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]],\n", "\n", " [[ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " ...,\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan],\n", " [ nan, nan, nan, ..., nan,\n", " nan, nan]]], dtype=float32)
<xarray.Dataset>\n", "Dimensions: (z_t: 60)\n", "Coordinates:\n", " * z_t (z_t) float64 5.0 15.0 25.0 35.0 ... 4.875e+03 5.125e+03 5.375e+03\n", "Data variables:\n", " NO3 (z_t) float64 dask.array<chunksize=(60,), meta=np.ndarray>\n", " NO3_obs (z_t) float64 dask.array<chunksize=(60,), meta=np.ndarray>\n", " PO4 (z_t) float64 dask.array<chunksize=(60,), meta=np.ndarray>\n", " PO4_obs (z_t) float64 dask.array<chunksize=(60,), meta=np.ndarray>\n", " SiO3 (z_t) float64 dask.array<chunksize=(60,), meta=np.ndarray>\n", " SiO3_obs (z_t) float64 dask.array<chunksize=(60,), meta=np.ndarray>
array([5.000000e+00, 1.500000e+01, 2.500000e+01, 3.500000e+01, 4.500000e+01,\n", " 5.500000e+01, 6.500000e+01, 7.500000e+01, 8.500000e+01, 9.500000e+01,\n", " 1.050000e+02, 1.150000e+02, 1.250000e+02, 1.350000e+02, 1.450000e+02,\n", " 1.550000e+02, 1.650984e+02, 1.754790e+02, 1.862913e+02, 1.976603e+02,\n", " 2.097114e+02, 2.225783e+02, 2.364088e+02, 2.513702e+02, 2.676542e+02,\n", " 2.854837e+02, 3.051192e+02, 3.268680e+02, 3.510935e+02, 3.782276e+02,\n", " 4.087846e+02, 4.433777e+02, 4.827367e+02, 5.277280e+02, 5.793729e+02,\n", " 6.388626e+02, 7.075633e+02, 7.870025e+02, 8.788252e+02, 9.847059e+02,\n", " 1.106204e+03, 1.244567e+03, 1.400497e+03, 1.573946e+03, 1.764003e+03,\n", " 1.968944e+03, 2.186457e+03, 2.413972e+03, 2.649001e+03, 2.889385e+03,\n", " 3.133405e+03, 3.379793e+03, 3.627670e+03, 3.876452e+03, 4.125768e+03,\n", " 4.375392e+03, 4.625190e+03, 4.875083e+03, 5.125028e+03, 5.375000e+03])
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "
\n",
"
| \n",
"\n", "\n", " | \n", "