{ "cells": [ { "cell_type": "markdown", "id": "15c85d44", "metadata": {}, "source": [ "## Import packages" ] }, { "cell_type": "code", "execution_count": 1, "id": "9a723e0b", "metadata": { "ExecuteTime": { "end_time": "2022-07-14T14:21:53.039629Z", "start_time": "2022-07-14T14:21:49.985325Z" } }, "outputs": [], "source": [ "import random\n", "import xarray as xr\n", "import pandas as pd\n", "import seaborn as sns\n", "import matplotlib.pyplot as plt\n", "import binsreg\n", "from ipynb.fs.full.asm import binscatter\n", "sns.set_theme(style='darkgrid')" ] }, { "cell_type": "markdown", "id": "cda6be89", "metadata": {}, "source": [ "## Load tdiff, winds, and pblh (not the anomaly)" ] }, { "cell_type": "code", "execution_count": 2, "id": "1d312ccb", "metadata": { "ExecuteTime": { "end_time": "2022-07-14T14:21:53.075827Z", "start_time": "2022-07-14T14:21:53.041270Z" } }, "outputs": [], "source": [ "tdif = xr.open_dataset('/bog/amuttaqin/Datasets/Derived/tdiff_MalayPenin_Sumatra.nc')\n", "wind_anom = xr.open_dataset('/bog/amuttaqin/Datasets/Derived/windslm_anom_MalayPenin_Sumatra.nc')" ] }, { "cell_type": "code", "execution_count": 3, "id": "22ae37b8", "metadata": { "ExecuteTime": { "end_time": "2022-07-14T14:21:53.082088Z", "start_time": "2022-07-14T14:21:53.077126Z" } }, "outputs": [ { "data": { "text/plain": [ "Frozen({'time': 262992, 'segment': 153})" ] }, "execution_count": 3, "metadata": {}, "output_type": "execute_result" } ], "source": [ "tdif.dims" ] }, { "cell_type": "code", "execution_count": 4, "id": "2d87b349", "metadata": { "ExecuteTime": { "end_time": "2022-07-14T14:21:53.094875Z", "start_time": "2022-07-14T14:21:53.083574Z" } }, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n",
"Dimensions: (time: 262992, segment: 152)\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 1991-01-01 ... 2020-12-31T23:00:00\n",
" * segment (segment) int64 1 1 1 1 2 2 2 2 2 2 ... 31 31 31 32 32 32 33 33 33\n",
" hour (time) int64 ...\n",
"Data variables:\n",
" spd (time, segment) float32 ...<xarray.Dataset>\n",
"Dimensions: (time: 43824)\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 1991-01-01 ... 1995-12-31T23:00:00\n",
" hour (time) int64 0 1 2 3 4 5 6 7 8 9 ... 14 15 16 17 18 19 20 21 22 23\n",
"Data variables:\n",
" skt (time) float32 -4.793 -3.941 -2.912 -1.208 ... -5.293 -5.528 -5.692\n",
" spd (time) float32 0.9362 0.7959 0.6947 1.252 ... 2.5 2.622 3.429 3.385| \n", " | skt | \n", "hour | \n", "spd | \n", "
|---|---|---|---|
| time | \n", "\n", " | \n", " | \n", " |
| 1991-01-01 00:00:00 | \n", "-4.792969 | \n", "0 | \n", "0.936177 | \n", "
| 1991-01-01 01:00:00 | \n", "-3.941162 | \n", "1 | \n", "0.795923 | \n", "
| 1991-01-01 02:00:00 | \n", "-2.912354 | \n", "2 | \n", "0.694712 | \n", "
| 1991-01-01 03:00:00 | \n", "-1.208252 | \n", "3 | \n", "1.252381 | \n", "
| 1991-01-01 04:00:00 | \n", "0.089355 | \n", "4 | \n", "1.056002 | \n", "
| ... | \n", "... | \n", "... | \n", "... | \n", "
| 1995-12-31 19:00:00 | \n", "-4.867676 | \n", "19 | \n", "2.540793 | \n", "
| 1995-12-31 20:00:00 | \n", "-5.063721 | \n", "20 | \n", "2.500491 | \n", "
| 1995-12-31 21:00:00 | \n", "-5.293213 | \n", "21 | \n", "2.622429 | \n", "
| 1995-12-31 22:00:00 | \n", "-5.528076 | \n", "22 | \n", "3.429289 | \n", "
| 1995-12-31 23:00:00 | \n", "-5.692139 | \n", "23 | \n", "3.385187 | \n", "
43824 rows × 3 columns
\n", "| \n", " | group | \n", "skt | \n", "bin | \n", "isknot | \n", "mid | \n", "spd | \n", "ci_l | \n", "ci_r | \n", "ci | \n", "
|---|---|---|---|---|---|---|---|---|---|
| 0 | \n", "Full Sample | \n", "-7.120469 | \n", "0 | \n", "0 | \n", "0 | \n", "0.031679 | \n", "-0.091781 | \n", "0.155061 | \n", "0.246841 | \n", "
| 1 | \n", "Full Sample | \n", "-6.599748 | \n", "1 | \n", "0 | \n", "0 | \n", "-0.232770 | \n", "-0.316538 | \n", "-0.123449 | \n", "0.193089 | \n", "
| 2 | \n", "Full Sample | \n", "-6.388688 | \n", "2 | \n", "0 | \n", "0 | \n", "-0.175899 | \n", "-0.324673 | \n", "-0.128822 | \n", "0.195851 | \n", "
| 3 | \n", "Full Sample | \n", "-6.234826 | \n", "3 | \n", "0 | \n", "0 | \n", "-0.379452 | \n", "-0.434728 | \n", "-0.237376 | \n", "0.197352 | \n", "
| 4 | \n", "Full Sample | \n", "-6.104795 | \n", "4 | \n", "0 | \n", "0 | \n", "-0.337174 | \n", "-0.436070 | \n", "-0.239627 | \n", "0.196443 | \n", "