{ "cells": [ { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "# Get Specific Months Data and annual mean\n", "\n", "## spec_moth_dat\n", "\n", "get specific months' data \n", "\n", "### Parameter:\n", "+ DaArray (xr.DataArray): shape = (time, *number) original Dataarray\n", "+ months (list): get data from specific month\n", "\n", "### Raises:\n", "+ ValueError: xr.DataArray' input is required, not the %s\n", "\n", "### Returns:\n", "+ xr.DataArray: data in specific month" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "### Example" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Load Modules" ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import sacpy as scp" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "#### Load and Get Data" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ ":219: RuntimeWarning: numpy.ndarray size changed, may indicate binary incompatibility. Expected 16 from C header, got 96 from PyObject\n" ] }, { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'sst' (time: 62, lat: 13, lon: 57)>\n",
       "[45942 values with dtype=float64]\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1991-01-16T12:00:00 ... 2021-02-16T12:00:00\n",
       "  * lon      (lon) int64 20 25 30 35 40 45 50 55 ... 270 275 280 285 290 295 300\n",
       "  * lat      (lat) int64 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30
" ], "text/plain": [ "\n", "[45942 values with dtype=float64]\n", "Coordinates:\n", " * time (time) datetime64[ns] 1991-01-16T12:00:00 ... 2021-02-16T12:00:00\n", " * lon (lon) int64 20 25 30 35 40 45 50 55 ... 270 275 280 285 290 295 300\n", " * lat (lat) int64 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sst = scp.load_sst()['sst']\n", "scp.spec_moth_dat(sst,[1,2])" ] }, { "attachments": {}, "cell_type": "markdown", "metadata": {}, "source": [ "## spec_moth_yrmean\n", "\n", "get specific month data and average them in each year\n", "\n", "### Parameter:\n", "+ DaArray (xr.DataArray): shape = (time, *number) original Dataarray\n", "+ months (list): get data from specific month\n", "\n", "### Returns:\n", "+ xr.DataArray: data in specific month and average them in each year\n", "\n", "### Example" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'sst' (time: 30, lat: 13, lon: 57)>\n",
       "array([[[        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",
       "         25.82674297, 25.56359132],\n",
       "        [        nan,         nan,         nan, ..., 24.08079465,\n",
       "         23.95300039, 23.89301427],\n",
       "        [        nan,         nan,         nan, ..., 21.20026779,\n",
       "         21.11359294, 21.35712878]],\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",
       "         27.00849136, 26.57378038],\n",
       "        [        nan,         nan,         nan, ..., 25.40256437,\n",
       "         25.14014912, 24.91551399],\n",
       "        [        nan,         nan,         nan, ..., 22.52778212,\n",
       "         22.33078369, 22.41735236]],\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",
       "         26.71056239, 26.41532103],\n",
       "        [        nan,         nan,         nan, ..., 25.12457975,\n",
       "         24.86658351, 24.81672478],\n",
       "        [        nan,         nan,         nan, ..., 22.16162841,\n",
       "         22.10272376, 22.36174266]]])\n",
       "Coordinates:\n",
       "  * time     (time) int64 1991 1992 1993 1994 1995 ... 2016 2017 2018 2019 2020\n",
       "  * lon      (lon) int64 20 25 30 35 40 45 50 55 ... 270 275 280 285 290 295 300\n",
       "  * lat      (lat) int64 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30
" ], "text/plain": [ "\n", "array([[[ 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", " 25.82674297, 25.56359132],\n", " [ nan, nan, nan, ..., 24.08079465,\n", " 23.95300039, 23.89301427],\n", " [ nan, nan, nan, ..., 21.20026779,\n", " 21.11359294, 21.35712878]],\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", " 27.00849136, 26.57378038],\n", " [ nan, nan, nan, ..., 25.40256437,\n", " 25.14014912, 24.91551399],\n", " [ nan, nan, nan, ..., 22.52778212,\n", " 22.33078369, 22.41735236]],\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", " 26.71056239, 26.41532103],\n", " [ nan, nan, nan, ..., 25.12457975,\n", " 24.86658351, 24.81672478],\n", " [ nan, nan, nan, ..., 22.16162841,\n", " 22.10272376, 22.36174266]]])\n", "Coordinates:\n", " * time (time) int64 1991 1992 1993 1994 1995 ... 2016 2017 2018 2019 2020\n", " * lon (lon) int64 20 25 30 35 40 45 50 55 ... 270 275 280 285 290 295 300\n", " * lat (lat) int64 -30 -25 -20 -15 -10 -5 0 5 10 15 20 25 30" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "sst = scp.load_sst()['sst']\n", "scp.spec_moth_yrmean(sst,[12,1,2])" ] } ], "metadata": { "kernelspec": { "display_name": "sacpy", "language": "python", "name": "python3" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.8.16" }, "orig_nbformat": 4 }, "nbformat": 4, "nbformat_minor": 2 }