{ "cells": [ { "cell_type": "markdown", "id": "be617771", "metadata": { "ExecuteTime": { "end_time": "2021-06-07T18:02:13.379282Z", "start_time": "2021-06-07T18:02:13.377450Z" }, "code_folding": [] }, "source": [ "# The Pyleoclim User Interface" ] }, { "cell_type": "code", "execution_count": null, "id": "2b991b26-5523-4cb0-8b8b-72de80903644", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "4748faa2", "metadata": {}, "source": [ "## Series" ] }, { "cell_type": "code", "execution_count": 22, "id": "2481a572-c427-4f50-8a94-3aa429e48dda", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\n", "[2021-06-15_10:29:04]: start //home/wenchang/wython/wystart.py\n", "[imported]: os.path, sys, os, datetime, glob\n", "[imported]: xarray(0.18.2) as xr, numpy(1.20.3) as np, pandas(1.2.4) as pd, matplotlib(3.4.2) as mpl\n", "[imported]: import matplotlib.pyplot as plt\n", "\n", "**wython plot settings (//home/wenchang/wython/wyconfig.py)**\n", "[config]: plt.rcParams['figure.dpi'] = 128.0\n", "[config]: plt.rcParams['figure.figsize'] = [6.4, 6.4*9/16]\n", "[config]: plt.rcParams['figure.constrained_layout.use'] = True\n", " plt.tight_layout() is NOT needed if plt.rcParams['figure.constrained_layout.use'] = True\n", "[registered colormaps]: tc and tc_r\n", "[registered colormaps]: parula and parula_r\n", "[imported]: import misc.colormaps\n", "[config]: xr.set_options(cmap_sequential=\"parula\")\n", "[shortcut functions]:\n", " constrained_layout_on(): plt.rcParams['figure.constrained_layout.use'] = True\n", " constrained_layout_off(): plt.rcParams['figure.constrained_layout.use'] = False\n", "[iPython config]: InlineBackend.figure_format ='retina'\n", "[imported]: from matplotlib.pyplot import plot, figure, close, show\n", "[executed]: plt.ion()\n", "\n", "[2021-06-15_10:29:04]: **done**; **0** seconds from \"start //home/wenchang/wython/wystart.py\"\n" ] } ], "source": [ "run -im wystart" ] }, { "cell_type": "code", "execution_count": 1, "id": "9a636e17", "metadata": { "ExecuteTime": { "end_time": "2021-06-07T18:03:59.950522Z", "start_time": "2021-06-07T18:03:27.069291Z" } }, "outputs": [], "source": [ "import pyleoclim as pyleo\n", "import pandas as pd" ] }, { "cell_type": "code", "execution_count": 2, "id": "27df6385", "metadata": { "ExecuteTime": { "end_time": "2021-06-07T18:04:01.880758Z", "start_time": "2021-06-07T18:04:01.723104Z" } }, "outputs": [], "source": [ "data=pd.read_csv(\n", " 'https://raw.githubusercontent.com/LinkedEarth/Pyleoclim_util/Development/example_data/soi_data.csv',skiprows=0, header=1)" ] }, { "cell_type": "code", "execution_count": 3, "id": "49463e75", "metadata": { "ExecuteTime": { "end_time": "2021-06-07T18:04:23.758526Z", "start_time": "2021-06-07T18:04:23.700479Z" } }, "outputs": [ { "data": { "text/html": [ "
| \n", " | Date | \n", "Year | \n", "Value | \n", "
|---|---|---|---|
| 0 | \n", "195101 | \n", "1951.000000 | \n", "1.5 | \n", "
| 1 | \n", "195102 | \n", "1951.083333 | \n", "0.9 | \n", "
| 2 | \n", "195103 | \n", "1951.166667 | \n", "-0.1 | \n", "
| 3 | \n", "195104 | \n", "1951.250000 | \n", "-0.3 | \n", "
| 4 | \n", "195105 | \n", "1951.333333 | \n", "-0.7 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "
| 823 | \n", "201908 | \n", "2019.583333 | \n", "-0.1 | \n", "
| 824 | \n", "201909 | \n", "2019.666667 | \n", "-1.2 | \n", "
| 825 | \n", "201910 | \n", "2019.750000 | \n", "-0.4 | \n", "
| 826 | \n", "201911 | \n", "2019.833333 | \n", "-0.8 | \n", "
| 827 | \n", "201912 | \n", "2019.916667 | \n", "-0.6 | \n", "
828 rows × 3 columns
\n", "| \n", " | t | \n", "air | \n", "nino | \n", "
|---|---|---|---|
| 0 | \n", "1871.00 | \n", "87.36090 | \n", "-0.358250 | \n", "
| 1 | \n", "1871.08 | \n", "-21.83460 | \n", "-0.292458 | \n", "
| 2 | \n", "1871.17 | \n", "-5.52632 | \n", "-0.143583 | \n", "
| 3 | \n", "1871.25 | \n", "75.73680 | \n", "-0.149625 | \n", "
| 4 | \n", "1871.33 | \n", "105.82000 | \n", "-0.274250 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "
| 1591 | \n", "2003.58 | \n", "-206.03800 | \n", "0.238497 | \n", "
| 1592 | \n", "2003.67 | \n", "103.90200 | \n", "0.411449 | \n", "
| 1593 | \n", "2003.75 | \n", "218.24100 | \n", "0.592756 | \n", "
| 1594 | \n", "2003.83 | \n", "-154.66200 | \n", "0.664131 | \n", "
| 1595 | \n", "2003.92 | \n", "15.04510 | \n", "0.604324 | \n", "
1596 rows × 3 columns
\n", "<xarray.Dataset>\n",
"Dimensions: (dim_0: 1596)\n",
"Coordinates:\n",
" * dim_0 (dim_0) int64 0 1 2 3 4 5 6 ... 1589 1590 1591 1592 1593 1594 1595\n",
"Data variables:\n",
" t (dim_0) float64 1.871e+03 1.871e+03 ... 2.004e+03 2.004e+03\n",
" air (dim_0) float64 87.36 -21.83 -5.526 75.74 ... 218.2 -154.7 15.05\n",
" nino (dim_0) float64 -0.3583 -0.2925 -0.1436 ... 0.5928 0.6641 0.6043array([ 0, 1, 2, ..., 1593, 1594, 1595])
array([1871. , 1871.08, 1871.17, ..., 2003.75, 2003.83, 2003.92])
array([ 87.3609 , -21.8346 , -5.52632, ..., 218.241 , -154.662 ,\n",
" 15.0451 ])array([-0.35825 , -0.2924584, -0.1435833, ..., 0.5927562, 0.6641306,\n",
" 0.6043244])<xarray.Dataset>\n",
"Dimensions: (dim_0: 1596)\n",
"Coordinates:\n",
" * dim_0 (dim_0) int64 0 1 2 3 4 5 6 ... 1590 1591 1592 1593 1594 1595\n",
"Data variables:\n",
" slope float64 -46.06\n",
" intercept float64 3.6\n",
" r float64 -0.1524\n",
" p float64 9.408e-10\n",
" stderr float64 7.482\n",
" predicted (dim_0) float64 20.1 17.07 10.21 10.49 ... -23.7 -26.99 -24.23\n",
" Nef int32 1217\n",
" pe float64 9.187e-08array([ 0, 1, 2, ..., 1593, 1594, 1595])
array(-46.06003419)
array(3.60029489)
array(-0.15239413)
array(9.40792838e-10)
array(7.481856)
array([ 20.10130214, 17.0709388 , 10.2137466 , ..., -23.70207595,\n",
" -26.98958325, -24.23490763])array(1217, dtype=int32)
array(9.18669389e-08)
<xarray.Dataset>\n",
"Dimensions: (year: 1596)\n",
"Dimensions without coordinates: year\n",
"Data variables:\n",
" t (year) float64 1.871e+03 1.871e+03 ... 2.004e+03 2.004e+03\n",
" air (year) float64 87.36 -21.83 -5.526 75.74 ... 218.2 -154.7 15.05\n",
" nino (year) float64 -0.3583 -0.2925 -0.1436 ... 0.5928 0.6641 0.6043array([1871. , 1871.08, 1871.17, ..., 2003.75, 2003.83, 2003.92])
array([ 87.3609 , -21.8346 , -5.52632, ..., 218.241 , -154.662 ,\n",
" 15.0451 ])array([-0.35825 , -0.2924584, -0.1435833, ..., 0.5927562, 0.6641306,\n",
" 0.6043244])<xarray.Dataset>\n",
"Dimensions: (year: 1516)\n",
"Coordinates:\n",
" * year (year) int64 500 501 502 503 504 ... 2012 2013 2014 2015\n",
"Data variables:\n",
" normTCcounts (year) float64 0.01712 0.01706 0.017 ... 0.01272 0.01105\n",
" nopost1700Vieq (year) float64 0.01712 0.01706 0.017 ... 0.01386 0.01202\n",
" smoothlower (year) float64 0.01361 0.01357 0.01352 ... 0.008644 0.006937\n",
" smoothupper (year) float64 0.02063 0.02056 0.02049 ... 0.0168 0.01517array([ 500, 501, 502, ..., 2013, 2014, 2015])
array([0.01711906, 0.01706153, 0.01700273, ..., 0.01439542, 0.01272127,\n",
" 0.01105489])array([0.01711906, 0.01706153, 0.01700273, ..., 0.01570766, 0.0138603 ,\n",
" 0.01201801])array([0.01360814, 0.01356591, 0.01351728, ..., 0.0102868 , 0.00864391,\n",
" 0.00693715])array([0.02062998, 0.02055714, 0.02048818, ..., 0.01850404, 0.01679863,\n",
" 0.01517264])<xarray.Dataset>\n",
"Dimensions: (mcrun: 20, year: 2001)\n",
"Coordinates:\n",
" * year (year) int32 0 1 2 3 4 5 6 7 ... 1994 1995 1996 1997 1998 1999 2000\n",
"Dimensions without coordinates: mcrun\n",
"Data variables:\n",
" HU (year, mcrun) float32 5.234 6.442 6.044 5.527 ... 5.552 5.026 5.276\n",
" MH (year, mcrun) float32 0.8484 1.135 0.9549 ... 0.9613 0.8497 0.8651\n",
" MHnoCO2 (year, mcrun) float32 0.7718 1.032 0.8687 ... 0.977 0.8636 0.8792\n",
" TS (year, mcrun) float32 7.701 8.858 8.373 7.981 ... 8.213 7.7 7.892\n",
" PDI (year, mcrun) float32 1.764 2.267 1.97 1.875 ... 2.13 1.912 1.951\n",
" MDR (year, mcrun) float32 -0.03904 -0.04728 -0.1791 ... 0.1807 0.1124\n",
" TROP (year, mcrun) float32 -0.05526 -0.04614 -0.146 ... 0.2754 0.1624\n",
" co2 (year) float32 277.5 277.1 277.2 277.2 ... 363.3 365.9 367.8 369.1array([ 0, 1, 2, ..., 1998, 1999, 2000], dtype=int32)
array([[5.233572, 6.442211, 6.044354, ..., 6.006451, 5.781193, 5.859242],\n",
" [5.150266, 6.517063, 6.253347, ..., 6.111382, 5.597956, 6.046769],\n",
" [5.273979, 6.42114 , 6.063237, ..., 6.156029, 5.663002, 5.857195],\n",
" ...,\n",
" [8.607759, 5.735194, 6.309206, ..., 5.426043, 7.792073, 6.349372],\n",
" [5.951543, 4.817696, 4.865606, ..., 5.416904, 5.246526, 5.486865],\n",
" [5.562082, 5.253688, 5.364417, ..., 5.551528, 5.025664, 5.275593]],\n",
" dtype=float32)array([[0.848376, 1.134567, 0.954869, ..., 1.034443, 0.954215, 1.024907],\n",
" [0.824386, 1.162915, 1.042043, ..., 1.056151, 0.915041, 1.063733],\n",
" [0.870265, 1.138879, 0.997377, ..., 1.070725, 0.929367, 1.000787],\n",
" ...,\n",
" [2.123141, 1.202116, 1.369071, ..., 1.021327, 1.726748, 1.324007],\n",
" [1.049104, 0.79278 , 0.799286, ..., 0.923721, 0.886139, 0.9061 ],\n",
" [0.991126, 0.916043, 0.935973, ..., 0.961338, 0.849738, 0.865122]],\n",
" dtype=float32)array([[0.771821, 1.032187, 0.868704, ..., 0.941097, 0.868109, 0.932423],\n",
" [0.749663, 1.057508, 0.947592, ..., 0.960421, 0.832101, 0.967316],\n",
" [0.79141 , 1.035684, 0.907003, ..., 0.973705, 0.845156, 0.910105],\n",
" ...,\n",
" [2.150444, 1.217575, 1.386677, ..., 1.034461, 1.748953, 1.341034],\n",
" [1.064745, 0.8046 , 0.811203, ..., 0.937493, 0.89935 , 0.919609],\n",
" [1.007259, 0.930954, 0.951208, ..., 0.976986, 0.863569, 0.879204]],\n",
" dtype=float32)array([[ 7.700802, 8.858367, 8.373077, ..., 8.45753 , 8.206414, 8.35135 ],\n",
" [ 7.609796, 8.939194, 8.62148 , ..., 8.551465, 8.03367 , 8.520521],\n",
" [ 7.758632, 8.850148, 8.443143, ..., 8.598428, 8.095691, 8.317599],\n",
" ...,\n",
" [11.35288 , 8.633804, 9.20207 , ..., 8.203881, 10.501794, 9.178846],\n",
" [ 8.591447, 7.469404, 7.512539, ..., 8.069045, 7.901129, 8.090614],\n",
" [ 8.2602 , 7.949969, 8.052602, ..., 8.212787, 7.699562, 7.89198 ]],\n",
" dtype=float32)array([[1.764139, 2.267454, 1.969975, ..., 2.092435, 1.956825, 2.070539],\n",
" [1.721461, 2.313551, 2.115247, ..., 2.130306, 1.886035, 2.13883 ],\n",
" [1.800178, 2.272291, 2.0369 , ..., 2.154965, 1.911669, 2.032453],\n",
" ...,\n",
" [4.136564, 2.531439, 2.833231, ..., 2.217835, 3.483728, 2.764105],\n",
" [2.296548, 1.800329, 1.814099, ..., 2.056627, 1.983548, 2.030182],\n",
" [2.180596, 2.036992, 2.076778, ..., 2.129763, 1.911505, 1.95099 ]],\n",
" dtype=float32)array([[-0.039044, -0.047281, -0.179137, ..., -0.067344, -0.104362, -0.003234],\n",
" [-0.057311, -0.028303, -0.095429, ..., -0.06747 , -0.112276, -0.003878],\n",
" [-0.010471, -0.034356, -0.108677, ..., -0.05825 , -0.108774, -0.045748],\n",
" ...,\n",
" [ 0.74415 , 0.491981, 0.565559, ..., 0.315151, 0.518547, 0.476624],\n",
" [ 0.265052, 0.118918, 0.154655, ..., 0.142988, 0.162764, 0.110104],\n",
" [ 0.308246, 0.197182, 0.233706, ..., 0.165041, 0.180703, 0.112441]],\n",
" dtype=float32)array([[-0.055265, -0.046137, -0.146011, ..., -0.048392, -0.076769, -0.012129],\n",
" [-0.061385, -0.036412, -0.091972, ..., -0.059893, -0.062815, -0.033429],\n",
" [-0.034247, -0.032187, -0.083764, ..., -0.056264, -0.067214, -0.050696],\n",
" ...,\n",
" [ 0.332312, 0.522402, 0.505372, ..., 0.367471, 0.295421, 0.372951],\n",
" [ 0.137719, 0.296573, 0.301219, ..., 0.211471, 0.230801, 0.134469],\n",
" [ 0.221633, 0.311034, 0.309191, ..., 0.215456, 0.275448, 0.162418]],\n",
" dtype=float32)array([277.454 , 277.137 , 277.16 , ..., 365.93295, 367.84497, 369.12497],\n",
" dtype=float32)<xarray.DataArray 'HU' (year: 1001)>\n",
"array([5.9251113, 5.7018557, 5.58632 , ..., 5.843549 , 5.6326075,\n",
" 5.0078006], dtype=float32)\n",
"Coordinates:\n",
" * year (year) int32 850 851 852 853 854 855 ... 1846 1847 1848 1849 1850array([5.9251113, 5.7018557, 5.58632 , ..., 5.843549 , 5.6326075,\n",
" 5.0078006], dtype=float32)array([ 850, 851, 852, ..., 1848, 1849, 1850], dtype=int32)
<xarray.DataArray 'normTCcounts' (year: 1001)>\n",
"array([0.0183048 , 0.01841119, 0.01853156, ..., 0.02311869, 0.02252923,\n",
" 0.02197403])\n",
"Coordinates:\n",
" * year (year) int64 850 851 852 853 854 855 ... 1846 1847 1848 1849 1850array([0.0183048 , 0.01841119, 0.01853156, ..., 0.02311869, 0.02252923,\n",
" 0.02197403])array([ 850, 851, 852, ..., 1848, 1849, 1850])
<xarray.Dataset>\n",
"Dimensions: (year: 1001)\n",
"Coordinates:\n",
" * year (year) int32 850 851 852 853 854 855 ... 1846 1847 1848 1849 1850\n",
"Data variables:\n",
" slope float64 0.004538\n",
" intercept float64 0.00153\n",
" r float64 0.31\n",
" p float64 9.793e-24\n",
" stderr float64 0.0004404\n",
" predicted (year) float64 0.02842 0.02741 0.02688 ... 0.02709 0.02426array([ 850, 851, 852, ..., 1848, 1849, 1850], dtype=int32)
array(0.0045384)
array(0.0015299)
array(0.30998597)
array(9.79348518e-24)
array(0.00044039)
array([0.02842043, 0.02740721, 0.02688286, ..., 0.02805027, 0.02709293,\n",
" 0.02425731])<xarray.Dataset>\n",
"Dimensions: (year: 1001)\n",
"Dimensions without coordinates: year\n",
"Data variables:\n",
" slope float64 0.004538\n",
" intercept float64 0.00153\n",
" r float64 0.31\n",
" p float64 9.793e-24\n",
" stderr float64 0.0004404\n",
" predicted (year) float64 0.02842 0.02741 0.02688 ... 0.02709 0.02426array(0.0045384)
array(0.0015299)
array(0.30998597)
array(9.79348518e-24)
array(0.00044039)
array([0.02842043, 0.02740721, 0.02688286, ..., 0.02805027, 0.02709293,\n",
" 0.02425731])<xarray.DataArray 'HU' (year: 1001)>\n",
"array([5.9251113, 5.7018557, 5.58632 , ..., 5.843549 , 5.6326075,\n",
" 5.0078006], dtype=float32)\n",
"Dimensions without coordinates: yeararray([5.9251113, 5.7018557, 5.58632 , ..., 5.843549 , 5.6326075,\n",
" 5.0078006], dtype=float32)<xarray.DataArray 'HU' (year: 1001)>\n",
"array([5.9251113, 5.7018557, 5.58632 , ..., 5.843549 , 5.6326075,\n",
" 5.0078006], dtype=float32)\n",
"Dimensions without coordinates: yeararray([5.9251113, 5.7018557, 5.58632 , ..., 5.843549 , 5.6326075,\n",
" 5.0078006], dtype=float32)<xarray.DataArray 'normTCcounts' (year: 1001)>\n",
"array([0.0183048 , 0.01841119, 0.01853156, ..., 0.02311869, 0.02252923,\n",
" 0.02197403])\n",
"Dimensions without coordinates: yeararray([0.0183048 , 0.01841119, 0.01853156, ..., 0.02311869, 0.02252923,\n",
" 0.02197403])<xarray.DataArray 'HU' (year: 1001)>\n",
"array([5.72471602, 5.7225005 , 5.72102525, ..., 5.45894145, 5.45703961,\n",
" 5.45519645])\n",
"Dimensions without coordinates: yeararray([5.72471602, 5.7225005 , 5.72102525, ..., 5.45894145, 5.45703961,\n",
" 5.45519645])<xarray.Dataset>\n",
"Dimensions: (year: 1001)\n",
"Coordinates:\n",
" * year (year) int64 850 851 852 853 854 855 ... 1846 1847 1848 1849 1850\n",
"Data variables:\n",
" slope float64 0.01175\n",
" intercept float64 -0.03948\n",
" r float64 0.4995\n",
" p float64 2.787e-64\n",
" stderr float64 0.000645\n",
" predicted (year) float64 0.02781 0.02778 0.02777 ... 0.02466 0.02464\n",
" Nef int32 1\n",
" pe float64 nanarray([ 850, 851, 852, ..., 1848, 1849, 1850])
array(0.01175475)
array(-0.03948217)
array(0.4994817)
array(2.78687362e-64)
array(0.00064505)
array([0.02781046, 0.02778441, 0.02776707, ..., 0.02468634, 0.02466399,\n",
" 0.02464232])array(1, dtype=int32)
array(nan)