{ "cells": [ { "cell_type": "markdown", "id": "eca44042", "metadata": { "papermill": { "duration": 0.023375, "end_time": "2021-06-07T20:43:04.877986", "exception": false, "start_time": "2021-06-07T20:43:04.854611", "status": "completed" }, "tags": [] }, "source": [ "# Analysis of the Gravity wave drag" ] }, { "cell_type": "code", "execution_count": 1, "id": "d21abe05", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:04.957934Z", "iopub.status.busy": "2021-06-07T20:43:04.957239Z", "iopub.status.idle": "2021-06-07T20:43:05.847823Z", "shell.execute_reply": "2021-06-07T20:43:05.846720Z" }, "papermill": { "duration": 0.931382, "end_time": "2021-06-07T20:43:05.848135", "exception": false, "start_time": "2021-06-07T20:43:04.916753", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "from netCDF4 import Dataset, num2date\n", "import xarray as xr\n", "import xarray.ufuncs as xrf\n", "import numpy as np\n", "import matplotlib.pyplot as plt\n", "import pandas as pd\n", "import glob" ] }, { "cell_type": "markdown", "id": "bc349ecf", "metadata": { "papermill": { "duration": 0.006214, "end_time": "2021-06-07T20:43:05.866740", "exception": false, "start_time": "2021-06-07T20:43:05.860526", "status": "completed" }, "tags": [] }, "source": [ "## Definitions" ] }, { "cell_type": "code", "execution_count": 2, "id": "cb9688c1", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:05.886105Z", "iopub.status.busy": "2021-06-07T20:43:05.885689Z", "iopub.status.idle": "2021-06-07T20:43:05.888634Z", "shell.execute_reply": "2021-06-07T20:43:05.888952Z" }, "papermill": { "duration": 0.015196, "end_time": "2021-06-07T20:43:05.889087", "exception": false, "start_time": "2021-06-07T20:43:05.873891", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "def anomalie(step,var):\n", " string = \"time.\"+step \n", " climatology = var.groupby(string).mean(\"time\")\n", " anomalies = var.groupby(string) - climatology\n", " \n", " return anomalies" ] }, { "cell_type": "code", "execution_count": 3, "id": "374d5895", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:05.913535Z", "iopub.status.busy": "2021-06-07T20:43:05.913116Z", "iopub.status.idle": "2021-06-07T20:43:05.915996Z", "shell.execute_reply": "2021-06-07T20:43:05.915668Z" }, "papermill": { "duration": 0.01975, "end_time": "2021-06-07T20:43:05.916075", "exception": false, "start_time": "2021-06-07T20:43:05.896325", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "def sea(days_period,station,var):\n", " df_dates = pd.read_csv('/home/hochatmstud/bene/dates/without_final_warmings/ssw_dates_without_final_warmings_'+station+'.csv') # you can load SSWs from a csv file like attached\n", " dates = df_dates.set_index('BeginDate')\n", "\n", " xa_ls = []\n", " max_lag = days_period\n", " for il, days in enumerate(range(-max_lag,max_lag+1)): \n", " sel_dates = pd.to_datetime(dates.index) + pd.Timedelta(str(days)+' days')\n", " mask = np.in1d(var.time.dt.floor('1D'),sel_dates)\n", " comp_m = var.sel(time = mask).mean('time')\n", " xa_ls.append(comp_m)\n", " \n", " xa_comp = xr.concat(xa_ls, dim = 'days')\n", " xa_comp['days'] = range(-max_lag, max_lag+1)\n", "\n", " # print(xa_comp)\n", " return xa_comp" ] }, { "cell_type": "markdown", "id": "9afbcb68", "metadata": { "papermill": { "duration": 0.007335, "end_time": "2021-06-07T20:43:05.933613", "exception": false, "start_time": "2021-06-07T20:43:05.926278", "status": "completed" }, "tags": [] }, "source": [ "## Analysis" ] }, { "cell_type": "code", "execution_count": 4, "id": "7d8ad0fe", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:05.950674Z", "iopub.status.busy": "2021-06-07T20:43:05.950250Z", "iopub.status.idle": "2021-06-07T20:43:05.953297Z", "shell.execute_reply": "2021-06-07T20:43:05.952938Z" }, "papermill": { "duration": 0.013748, "end_time": "2021-06-07T20:43:05.953394", "exception": false, "start_time": "2021-06-07T20:43:05.939646", "status": "completed" }, "tags": [ "parameters" ] }, "outputs": [], "source": [ "station = 'Leipzig'" ] }, { "cell_type": "code", "execution_count": 5, "id": "facb7378", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:05.971127Z", "iopub.status.busy": "2021-06-07T20:43:05.970711Z", "iopub.status.idle": "2021-06-07T20:43:05.973449Z", "shell.execute_reply": "2021-06-07T20:43:05.973776Z" }, "papermill": { "duration": 0.013757, "end_time": "2021-06-07T20:43:05.973875", "exception": false, "start_time": "2021-06-07T20:43:05.960118", "status": "completed" }, "tags": [ "injected-parameters" ] }, "outputs": [], "source": [ "# Parameters\n", "station = \"Andenes\"\n" ] }, { "cell_type": "code", "execution_count": 6, "id": "2a72c17a", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:05.995091Z", "iopub.status.busy": "2021-06-07T20:43:05.994662Z", "iopub.status.idle": "2021-06-07T20:43:05.997256Z", "shell.execute_reply": "2021-06-07T20:43:05.996911Z" }, "papermill": { "duration": 0.016401, "end_time": "2021-06-07T20:43:05.997350", "exception": false, "start_time": "2021-06-07T20:43:05.980949", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "if station == 'Leipzig':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_2004-08-01_2020-12-01_Col.nc'\n", "elif station == 'CMOR':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_2002-01-01_2020-11-01_CMA.nc'\n", "elif station == 'Esrange':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_1999-08-01_2021-03-01_Kir.nc'\n", "elif station == 'Sodankyla':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_2008-10-01_2021-03-01_Sod.nc'\n", "elif station == 'Davis':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_2005-01-01_2020-12-01_Dav.nc'\n", "elif station == 'RioGrande':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_2008-02-01_2021-01-01_Rio.nc'\n", "elif station == 'Andenes':\n", " dir = '/home/gemeinsam_tmp/UA_students/data/PW_GW_analysis/GWD_1999-08-01_2019-12-01_SES.nc'" ] }, { "cell_type": "code", "execution_count": 7, "id": "42caf289", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:06.021675Z", "iopub.status.busy": "2021-06-07T20:43:06.021275Z", "iopub.status.idle": "2021-06-07T20:43:06.040922Z", "shell.execute_reply": "2021-06-07T20:43:06.040139Z" }, "papermill": { "duration": 0.036342, "end_time": "2021-06-07T20:43:06.041131", "exception": false, "start_time": "2021-06-07T20:43:06.004789", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "ds=xr.open_dataset(dir)" ] }, { "cell_type": "code", "execution_count": 8, "id": "7b3a0cbe", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:06.073647Z", "iopub.status.busy": "2021-06-07T20:43:06.070562Z", "iopub.status.idle": "2021-06-07T20:43:06.076989Z", "shell.execute_reply": "2021-06-07T20:43:06.076448Z" }, "papermill": { "duration": 0.024815, "end_time": "2021-06-07T20:43:06.077133", "exception": false, "start_time": "2021-06-07T20:43:06.052318", "status": "completed" }, "tags": [] }, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "Show/Hide data repr\n", "\n", "\n", "\n", "\n", "\n", "Show/Hide attributes\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
xarray.Dataset
" ], "text/plain": [ "\n", "Dimensions: (alt: 26, time: 178992)\n", "Coordinates:\n", " month (time) int64 ...\n", " * alt (alt) int64 70 72 74 76 78 80 82 84 ... 108 110 112 114 116 118 120\n", " * time (time) datetime64[ns] 1999-08-01 ... 2019-12-31T23:00:00\n", "Data variables:\n", " GWD (time, alt) float64 ...\n", "Attributes:\n", " history: Tue Jun 1 17:49:59 2021: ncrename -v __xarray_dataarray_variab...\n", " NCO: netCDF Operators version 4.8.1 (Homepage = http://nco.sf.net, C..." ] }, "execution_count": 8, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ds" ] }, { "cell_type": "code", "execution_count": 9, "id": "0deb369b", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:06.110059Z", "iopub.status.busy": "2021-06-07T20:43:06.109555Z", "iopub.status.idle": "2021-06-07T20:43:22.773541Z", "shell.execute_reply": "2021-06-07T20:43:22.772388Z" }, "papermill": { "duration": 16.681522, "end_time": "2021-06-07T20:43:22.773917", "exception": false, "start_time": "2021-06-07T20:43:06.092395", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "gwd = ds['GWD'] \n", "\n", "gwd_anomalie = anomalie('day',gwd)\n", "\n", "sea_gwd = sea(40,station,gwd)\n", "sea_gwd_anomalie = sea(40,'Leipzig',gwd_anomalie)" ] }, { "cell_type": "code", "execution_count": 10, "id": "ff7a6abe", "metadata": { "execution": { "iopub.execute_input": "2021-06-07T20:43:22.808306Z", "iopub.status.busy": "2021-06-07T20:43:22.807690Z", "iopub.status.idle": "2021-06-07T20:43:22.990431Z", "shell.execute_reply": "2021-06-07T20:43:22.991473Z" }, "papermill": { "duration": 0.200648, "end_time": "2021-06-07T20:43:22.991774", "exception": false, "start_time": "2021-06-07T20:43:22.791126", "status": "completed" }, "tags": [] }, "outputs": [ { "data": { "text/plain": [ "" ] }, "execution_count": 10, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAaEAAAEGCAYAAADMsSqUAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4yLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+j8jraAAAgAElEQVR4nO3de5ycVZ3n8c+3qnMFAkm4iASHoBkUGAdJBlFHxxWB6DgCvmQIO2JGWaMOvkTHmRF0dnFVdmHHGWfnojNRWNBxQBYvsKOIEVF2ZrlFRLkNErkZiOGScI1Jurt++8dzOl3dqa7nVHdXVVf39/16Pa88dZ5LnXq6U78+z3PO7ygiMDMz64ZKtytgZmYzl4OQmZl1jYOQmZl1jYOQmZl1jYOQmZl1TV+3KzARi+bPjSX77NXtarRPs56LUt5+1pr66zrFTN2adUDuz6XJbj995PEnImK/yalQ4XVaHM8ywI/j6Rn945mIng5CS/bZi2+/5+RuV6NtYrA25jZVK1n7WWvqr+tUo8rUrVu75f5cmu33wo997qHJqg+ApHmLmEUfQtKvRcSknn+maNtvtaSLJT0m6c66slMl3SWpJmnFqP3PlbRB0r2STmxXvczMJsOrWLhtGXtwNHvzMvZ8sNv16VXt/NPqEmDlqLI7gbcBN9QXSjocWAUckY75nKRqG+tmZjZukubdy3McwV68iHlsoR9Jv9btevWitt2Oi4gbJB0yquweAO1+f/ck4PKI2AE8IGkDcAxwY7vqZ2Y2Xq9i4bYawaz0d/wrWMBD/OpBZviju/GYKjeZDwJ+Ufd6YyozM5tS6ltBQ9waGr+p0jGh0V8PDbt8SVoDrAE4aO8921knM7PdjG4FAQi5NTROU6UltBE4uO71EuDRRjtGxNqIWBERKxbNn9uRypmZQeNW0BC3hsZnqgShq4FVkuZIWgosA27pcp3MzEYY6hE3q8FX51BryD3lWtPOLtqXUXQsOEzSRklnSjpF0kbgVcC3JF0LEBF3AVcAdwPfAc6KiMF21c3MrFXNWkFD3BpqXTt7x50+xqZvjLH/+cD57aqPmdkEnfGSMVpBQ4Q4igVcy+MfBj7Uuar1rqnSMWFcqrNnseCQA8d9/FQeHd+KqZ5ZIWqN33d0fZ575PFd6wPbd7a1TmPqz9utkpG9oDbG5x6vnPecriYjY8IkOHAhs0p3SvuM/4tphunpIGRm1kkVoFrS963qVI4tcRAyM8tUkaiWJFOtoDEGmFgjDkJmZpn6BLMrzYPQrBC4W1U2ByEzs0wVZdyO60xVpg0HITOzTNXc23GWzUHIzCxTVseEjtRk+nAQMjPLVJWYVdIS6gu3hFrhIGRmlimnd1xZS8lGchAyM8uUczuupPOcjeIgZGaWqegdV9YSchRqRU8Hoepee7Pna99cvKgNd8yvzd9nxH4xZ4yEgzEyrUrl+S271hWdSXEz+OTIGSti+7as42JwCgxEqE1uHea95KWTVofdUgUNZKYBqgw/Vt65aePw2/QPtFy1yTCeVEujU9dMNF3T9qee3bW+7ZdPTuhc46UpkrIoq4u2Y1BLejoImZl10iypfLCqu2i3xEHIzCxTdtoey+YgZGaWybfjJp+DkJlZJveOm3wOQmZmmfoq5QlMZ3uwakvaOb33xZIek3RnXdkiSesk3Zf+XZjKD5H0K0m3p+Uf2lUvM7PxqqBd+ePGWiruot2SdvZ7vARYOarsHOC6iFgGXJdeD/l5RByVlve1sV5mZuMy9EyobLF8bQtCEXEDsGVU8UnApWn9UuDkdr2/mdlkq5S0gqqSe8e1qNMjwA6IiE0A6d/967YtlfRjST+U9NqxTiBpjaT1ktY/vuWpdtfXzGyXoUntmi/drmVvmSodEzYBL4qIJyUtB74p6YiIeGb0jhGxFlgLcPTy5bFjyVGTV4sFL5i8czWhqJv798Ajxt6v/1e71qvPjxqp3qGMDpOirq6qy1wQT20esVtljwXDL9T9/8mzl03wd6uXfkZN7D3JmTFy7Zb1Ysjo+tS9Hnjil8PH79w+cr+/vmzCdapkPPORnwm1pNP/0zdLOhAg/fsYQETsiIgn0/qPgJ8Dv97hupmZNaWKULV8sXydDkJXA6vT+mrgKgBJ+0mqpvVDgWXA/R2um5lZU6pAparmiwcKtaSdXbQvA24EDpO0UdKZwAXA8ZLuA45PrwFeB/xU0k+AK4H3RcToTg1mZl2lviqV2SXLrPK5VSUdLOl6SfdIukvS2am84TCWtO1cSRsk3SvpxLry5ZLuSNv+Rul+oKQ5kr6aym+WdEjdMavTe9wnaTVd1LZnQhFx+hibjmuw79eAr7WrLmZmk6GSWkLNKG9+7wHgIxFxm6S9gB9JWgf8IcUwlgsknUMxjOWjkg4HVgFHAC8Evifp1yNiEPg8sAa4Cfg2xdCYa4Azga0R8RJJq4ALgdMkLQLOA1YAkd776ojY2sKlmDTdf/prZtYrKhWUsZSJiE0RcVtafxa4BziIsYexnARcnp6fPwBsAI5Jz9YXRMSNERHAl0YdM3SuK4HjUivpRGBdRGxJgWcdu4/p7Jip0jvOzGzKy2kJpWdCh0paX1e8NvXs3U26TfYK4GZGDWORNDSM5SCKls6QjamsP62PLh865hfpXAOSngYW15c3OKbjHITMzHJVynu/pe33R8RpZaeTtCfFo4gPRcQzTbp3N9oQTcrHe0zH+XacmVmmal+F6uxq06UyK+9rVdIsigD0lYj4eipuOIyForVycN3hS4BHU/mSBuUjjpHUB+xNkcVmrHN1hYOQmVkuqRgrVLKUn0YCLgLuiYi/qtvUcBhLKl+VerwtpRjGcku6dfespGPTOd856pihc70d+H56bnQtcIKkhan33QmprCt6+nbctv5BfvzL5wHYf4/Zu8rnjmouzx6j+Tx6PPaOgeEWaS060zrdMTjyfQZHVKr+x3PAiP2mQEKB8anvOXTA0jF3yx1qUX/3Yr95w9er2U9vTt3vw85R17/+fedWhrdpdPaDutfVp8f4I7KXMybU/4JV6q5r35wuVCZf5YBl7T1/RVSqzf/zVTI6JgCvAc4A7pB0eyr7GMWwlSvSkJaHgVMBIuIuSVcAd1P0rDsr9YwDeD9Fwuh5FL3irknlFwFflrSBogW0Kp1ri6RPAbem/T7ZzSExPR2EzMw6KiMjQk7GhIj4Vxo/m4EGw1jSMecD5zcoXw8c2aB8OymINdh2MXBxaUU7wEHIzCxTJT0TarrPQK/epugOByEzs0zK6R3ntD0tcRAyM8skUZobzrnjWuMgZGaWSZUKKumYULbdRnIQMjPLVOkT1ZJZ6yo7HYRa4SBkZpZJ1ZyWkG/HtcJByMwsk/Jzx1kmByEzs0xSeUYE945rjYOQmVmmrHFCGZPa2bDeDkIxnOZm07M7dxV3KuVOu/XXhj/Htv7BEdtevHB+p6szbrW6JDr1P5qntw9M6vs89lz/hM9R/7vz863bMo/aY8Lva/lqtbH/fz9f9/9k58Dkp00qxgn5mdBkauf03hdLekzSnXVlLU9da2Y2VeRMaKeeTezYHe28Wpew+2x951BMXbsMuC69ZtTUtSuBz0mZk+SamXWIUgLTssXyte1qRcQNFJlb67U0dW276mZmNh5Dg1WbLnlZtC3p9DOhVqeu3Y2kNcAagANeuKTRLmZmbaFqhcqs5l+bcseElkyVkJ093WxErI2IFRGxYp9Fi9tcLTOzYUVLqNp8qTgItaLTLaHNkg5MraCcqWvNzKaOrNxx7h3Xik63hFqaurbDdTMza6pSEZVKpXSxfG1rCUm6DHg9sK+kjcB5jG/qWjOzKUF9VSqz/UxoMrUtCEXE6WNsamnqWjOzqaJI21NyO87jhFrS0xkT5vZVOGzxXAB+NTDcj2G6ZEyoKO/e8t5zhn/p5/bU/ejZLR8RTa7Jtv7hEfLP7Rw5Wr7+dX9teL3aJM/XrOpwJoTBJoPvn9i2c+yNGZrVwaYWzyc0+Xo6CJmZdVQ1J22Pg1ArHITMzDJVqlWqJeOEqn1+JtQKByEzs0w5CUxx77iWOAiZmWXKeybkZ3ytcBAyM8tVyegd55ZQSxyEzMwyVfr6qM6e1XSfsmdGNpKvlplZpqxJ7dwSaomDkJlZppzBqu6Y0BoHITOzTB6sOvl6OghVBQtmFz/wBa0Pvp/ycrMDzNNwmj1tf76tdZrK9pw9nOGA2dVR2/K+GOqv+LyFczLfeY/yXWxaUF/5fEKVjHFCki4G3gI8FhFHprJFwFeBQ4AHgd+PiK1p27nAmcAg8MGIuDaVL6eYxXoe8G3g7IgISXOALwHLgSeB0yLiwXTMauDPU1U+HRFDE412hUO2mVkmqZgvqGzJcAmwclTZOcB1EbEMuC69RtLhwCrgiHTM5yQNvcnnKSb5XJaWoXOeCWyNiJcAnwUuTOdaRJFM+pUUs1efJ2lhyxdiEjkImZnlqlSgUi1Zyr9WI+IGYMuo4pOAoVbJpcDJdeWXR8SOiHgA2AAck+ZkWxARN0ZEULR8Tm5wriuB4yQJOBFYFxFbUitrHbsHw47q6dtxZmYdJZUHmSKL9qGS1teVro2ItSVnPyAiNgGkiT/3T+UHATfV7bcxlfWn9dHlQ8f8Ip1rQNLTwOL68gbHdIWDkJlZrjS9dzOpY8L9EXHaJL1ro4fD0aR8vMd0hW/HmZnlqvZB3+zmS3XcvaQ2p1tspH8fS+UbgYPr9lsCPJrKlzQoH3GMpD5gb4rbf2Odq2schMzMMkmVopt202XcueOuBlan9dXAVXXlqyTNkbSUogPCLenW3bOSjk3Pe9456pihc70d+H56bnQtcIKkhalDwgmprGu6EoQknS3pTkl3SfpQKvuEpEck3Z6WN3ejbmZmY8rqmJDVRfsy4EbgMEkbJZ0JXAAcL+k+4Pj0moi4C7gCuBv4DnBWRAyNy3g/8EWKzgo/B65J5RcBiyVtAP6Y1NMuIrYAnwJuTcsnU1nXdPyZkKQjgfdQdA/cCXxH0rfS5s9GxGc6XSczsyyqlAeZjOm9I+L0MTYdN8b+5wPnNyhfDxzZoHw7cOoY57oYuLi0kh3SjY4JLwNuiohtAJJ+CJzShXqYmbVE1Srqa57AtGy7jdSN23F3Aq+TtFjSfODNDD8o+4Ckn0q6eKwBVJLWSFovaf3jTzzRqTqbmQ23hCY4TsiGdbwlFBH3SLqQYpDUc8BPgAGKkb+fougu+CngL4F3Nzh+LbAWYPnRR7eta2GzlDm5FOXVa/Y+zY7fo6/uuP7+uoNm7n+AqPvsfaMeDs/2RGM2GTK6aOc8E7JhXfnGioiLIuLoiHgdRbfB+yJic0QMRkQN+ALFMyMzs6mjUslbLFtXBqtK2j8iHpP0IuBtwKskHTg0WpjiGdGd3aibmdlYVOlDfc3HAanqZ0Kt6FbGhK9JWkyRduKsiNgq6cuSjqK4Hfcg8N4u1c3MrLFKRtoet4Ra0pUgFBGvbVB2RjfqYmaWayiLdtN9HIRa4txxZma5KjnjhNwxoRUOQmZmmVStolkl44RKJr2zkXy1zMxyZWVMcEuoFU2DkKQ/zjjH8xHxj5NUHzOzqSvndpyfCbWk7Gr9KbAnsFeT5SPtrKCZ2ZSRlUXbQagVZbfjvhwRn2y2g6Q9JrE+U0ZOtoN2y83aoBmcJaGeorZrfc6oL4Lu/zRtWnBLaNI1DUIR8WdlJ8jZx8xsWlAFKiV/u8uP2luRdbUk7UMxYdIh9cdExAfbUy0zsylIlRE5Chsp224j5YbsbwM3AXcAtZJ9zcymKZUnCZ6E5MczSW4QmhsROT3lzMymL6k8yDgItSQ3CH1Z0nuAfwF2DBV2e1pYM7NOikqVqDb/2oyyZ0Y2Qu7V2gn8BfBxhjsaBXBoOyplZjYlybfj6kk6DFgDvDQV3QN8ISLuzT1HbhD6Y+AlEeGpTM1s5lIlIwjNjI4Jkl4FfB34R4qJRgW8Arhe0tsi4qac8+QGobuAbeOpqJnZ9JHTEpoZQQj4L8DpEfGDurJvSvo+cB7wppyT5AahQeB2Sdcz8pmQu2ib2YwRlWrpM5+YOdN7v3hUAAIgIn4oaW3uSXKD0DfTMqWElJ1VYFrzNSjE2KMHpkIGDJsm3DtuyLNNtj2fe5LcIHRnRPyovkDS7+W+iZnZtJDzTKg0Jee0cbCkv2lQLuCg3JPkBqEvSFodEXcASDod+BDwf3LfqJ6ks4H3UFT2CxHx15IWAV+lyMrwIPD7EbF1POc3M2uH4u5LWcaEGdMS+tMm29bnniQ3CL0duFLSHwC/TZHC54TcN6kn6UiKAHQMRdfv70j6Viq7LiIukHQOcA7w0fG8h5lZW1SqUDJOqHT79PET4CcRE7vXnXW1IuJ+Sasongv9AjghIn41zvd8GXBTRGwDkPRD4BTgJOD1aZ9LgR/gIGRmU4p7x9X5IrBU0m3AvwH/j+K7/ZlWTlI2qd0djMyCvwioAjdLIiJe3lqdAbgTOF/SYuBXwJspmm4HRMQmgIjYJGn/Meq0hmJwFAcffPA43t7MbJyyxgnNjNtxEbFC0nyKu1qvBj5IkV3nl8C/RcQf5ZynrCX0lolVc3cRcY+kC4F1wHMUTbqBFo5fSzEwiqOXL3eXJzPrnIxnQjOoYwLpjtYPJN0K3Ay8huJxzcrcc5TNJ/TQhGo49nkvAi4CkPTfgI3AZkkHplbQgcBj7XhvM7Nxy0rbkxeEJK0E/ifF3aUvRsQFE61eJ0n6jxQtoKMoxo8OBaLfjohf5p6n7HbcbRFx9ET3aXDM/hHxmKQXAW8DXgUsBVYDF6R/r2rlnGZmbafqpMysKqkK/D1wPMUf4bdKujoi7p6EWnbKWuDfgX8AboiIn43nJGW3414m6adNtgvYexzv+7X0TKgfOCsitkq6ALhC0pnAw8Cp4zivmVnb5HXRzmoJHQNsiIj7ASRdTtE5q5eC0N7Ab1K0hj6RkpluAm4EboyI7+ecpCwIvbRkOxQpfVoSEa9tUPYkcFyr5zIz65zs23GHSqofK7M2Pc8echBFT+MhG4FXTlIlOyIiBoHb0vJ3kg6gGM7zYeCTFLcZS3XlmZC1rmnamfptM6d7qE3QDBpUOWlyUoWl7fdHxGlNdmt0kp7qaCXp5RStoKFlNsUM3H9L0WU7y4wZVWVmNlG1CAZrzWNF2fZkI1A/xmQJ8Oj4a9YVl1AEm2uA/zzeRov/bDYza0FkLBluBZZJWippNrAKuLoN1W2nTwL3RsRlEfGQpJsl3Z+Wt+eeJCsIpXE9pWVmZtNZLcqXnCQ2ETEAfAC4lmI20isi4q721n7S/SkjezHPAX6LIvPN+3NPktsSOr5BWdaERWZm00UERETpkneu+HZE/HpEvDgizm9z1dthdkTUd67414h4MiIeBvbIPUnZOKH3A39E0dOjvqv2XrTw4MnMbDqoAYMlMaZs+zSysP5FRHyg7uV+uScp65jwzxQPnf47RVbrIc9GxJbcNzEzmw4i3XJrJq9fwrRws6T3RMQX6gslvRe4JfckZUGoCjwDnDV6g6RFDkRmNpPk3W6bMVHow8A3U/qe21LZcopnQyfnnqQsCP2I4Ss6ul97AIfmvpGZWa8LiltyzZRtny4i4jHg1ZLeAByRir+VmylhSNlg1aXjrJ+Z2bRTCxgsiTJl26ebFHRaCjz1sgerSloILAPm1r35DeN948mgiOaZBGaKmGG/9WPp0WwRzlzQO3Jux01wotEZJysISfpPwNkUo3pvB46lSFL3hvZVzcxsasm5HecQ1JrcPx3PphiE9FBE/AfgFcDjbauVmdkUFJG3WL7c23HbI2K7JCTNiYh/T2m7zcxmjKDIH9dM2XYbKTcIbZS0D/BNYJ2krfResj0zswmpRflg1Bk0TmhSZAWhiDglrX5C0vUUkxl9p221MjObgnJutzkGtablqRwi4oftqIiZ2VQXBLWSMFO23UbqSp9WSR+WdJekOyVdJmmupE9IekTS7Wl5czfqZmY2lsAdEyZbxye1k3QQ8EHg8Ij4laQrKObSAPhsRHym03UyM8tRq5UPRq152F5LujW6rw+YJ6kPmI87OZhZD6il23HNF2tFx4NQRDwCfAZ4GNgEPB0R302bPyDpp5IuThkadiNpjaT1ktY//sQTHaq19YTa4K4lpJ5ZrLf4dtzk6ngQSsHlJGAp8EJgD0nvAD4PvBg4iiI4/WWj4yNibUSsiIgV++27b4dqbWZWjAHKWSxfN27HvRF4ICIej4h+4OvAqyNic0QMRkQN+AJwTBfqZmY2psEa9A9G02XAA4Va0o0g9DBwrKT5kgQcB9wj6cC6fU4B7uxC3czMxlSLYDBjsXwd7x0XETdLupJiEqQB4MfAWuCLko6i6AX5IPDeTtfNzKyZnNttzqLdmo4HIYCIOA84b1TxGd2oi5lZLs8nNPm6EoTMzHrRYAT9JQOBBjy/V0schMzMMtWiPEGp+yW0xkHIzCxTLYLBkijjjAmtcRAyM8sU4fmEJpuDkJlZpoFa0F/SEvI4odY4CNm0Udn+zK716nNjzz4/uMfiXesxa15b62TTS87tOI8Tao2DkJlZJt+Om3wOQmZmmWp4eu/J1q2pHMzMek4nEphKOjVN+lmTtGLUtnMlbZB0r6QT68qXS7ojbfublBINSXMkfTWV3yzpkLpjVku6Ly2r68qXpn3vS8fOntAHKuEgZGaWabDWPHnpJCUwvRN4G3BDfaGkwykmAD0CWAl8TlI1bf48sAZYlpaVqfxMYGtEvAT4LHBhOtciiqw1r6RIFn1e3fQ5F1JMMLoM2JrO0TYOQmZmmSKrJTTh97gnIu5tsOkk4PKI2BERDwAbgGNS8ucFEXFjFInrvgScXHfMpWn9SuC41Eo6EVgXEVsiYiuwDliZtr0h7Us6duhcbeFnQmZmmWpR3vstbT9U0vq64rURsXaCb38QcFPd642prD+tjy4fOuYXABExIOlpYHF9+ahjFgNPRcRAg3O1hYOQmVmmWgS1kqZOFNvvj4jTxtpH0veAFzTY9PGIuGqswxq9XZPy8RzT7Fxt4SBkZpapSGDa/Du5bDtARLxxHG+/ETi47vUS4NFUvqRBef0xGyX1AXsDW1L560cd8wPgCWAfSX2pNVR/rrbwMyEzs0xDt+OaLbX2NRyuBlalHm9LKTog3BIRm4BnJR2bnum8E7iq7pihnm9vB76fnhtdC5wgaWHqkHACcG3adn3al3TsWC2zSeGW0BRWfeaXu9ZVGxhzv9qcPeoOmoTelD2Uir5Slxmhsv3ZXeu1559ptDsAVR5oa51sd7F92671gc0Pd7EmE5OVMWGCPRMknQL8LbAf8C1Jt0fEiRFxl6QrgLspJgQ9KyIG02HvBy4B5gHXpAXgIuDLkjZQtIBWAUTEFkmfAm5N+30yIrak9Y8Cl0v6NMWkoxdN6AOVcBAyM8s0WMvJoj2xIBQR3wC+Mca284HzG5SvB45sUL4dOHWMc10MXNyg/H6Kbtsd0ZXbcZI+nAZj3SnpMklzJS2StC4NkFpX12fdzGxKGBgMdg7Umi9lKRVshI4HIUkHAR8EVkTEkUCVool4DnBdGiB1XXptZjZlDN2Oa7ZMtCU003SrY0IfMC/11phP0fuiflBV2wdImZm1qhaUBqGJPhOaaToehCLiEeAzwMPAJuDpiPgucEDq5UH6d/9Gx0taI2m9pPWPP/FEp6ptZkYtIwA5i3ZrunE7biFFq2cp8EJgD0nvyD0+ItZGxIqIWLHfvvu2q5pmZrsZqJU8Dxqo0T/QO71Lp4Ju9I57I/BARDwOIOnrwKuBzZIOjIhNKRfSY12om5nZmGq18i7Yg45BLenGM6GHgWMlzU8Dq44D7mHkoKq2D5AyM2vVYBRZspstnlm1NR1vCUXEzZKuBG6jGHD1Y2AtsCdwhaQzKQJVw77tZmbdUsvoeOCOCa3pymDViDiPYi6LejsoWkVmZlPSUO+4ZhyEWtPTGROev/dn3PL63eOWKiMTwe770sWdqlLLVJ34HdE5++y1a706q6d/pB2x46nh9D4D23dmHbN96/YRr5/d9Nyk1WfOgjkjXu+x//xd69XZ1dG7WxcN1GrsLHno01/zQ6FW+BvLzCxTTscED1ZtjYOQmVmmwawEph2qzDThIGRmlqkWNQZLbrfVeigL/VTgIGRmlmkogWkzZdttJAchM7NMQ+OEmu7jZ0ItcRAyM8uUlTHBg1Vb4iBkZpYpZ2ZV99BujYOQmVmmnGdCTmDaGgchM7NMOS0h345rTU8HoSe39fPPtzxavuNNj7S/MtYVi+oyCswelSljPCbzC6SqidenmUPmz9q1Pq869nuNVY9OflnW/2z2XDh31/rcfeaO2G+vF+45ofcZ3DncCtmyYcuEztVI1IIoCUJl222kng5CZmadFFGeEcEZE1rjIGRmlilqNWolKRHcEmqNg5CZWaZaLdwSmmQOQmZmmSKgLCuPs/a0xkHIzCxXQJR06CjbbiM5CJmZZapF+e04B6HWdDwISToM+Gpd0aHAfwH2Ad4DPJ7KPxYR3+5w9czMxhSDQa1kMGrZdhup40EoIu4FjgKQVAUeAb4BvAv4bER8ptN1MjPLERHUSlo6ZdttpG7fjjsO+HlEPKQ2D+wzM5uoiIzBqo5BLel2EFoFXFb3+gOS3gmsBz4SEVtHHyBpDbAG4AVz5/Du3zl0t5PG4Mjfguc3Pz+JVbaZqG/eyP8qcxbMmbRzD+4cHPF6xzM7dq2P/l2ebvqf7x/xest9u/2XH7fqrGr5Ti0qese1N2OCpL8Afg/YCfwceFdEPJW2nQucCQwCH4yIa1P5cuASYB7wbeDsiAhJc4AvAcuBJ4HTIuLBdMxq4M/T2346Ii5N5UuBy4FFwG3AGRGxc0IfqolKu05cRtJs4K3A/05FnwdeTHGrbhPwl42Oi4i1EbEiIlbsM3t2R+pqZgZQq9UYHGy+lA1mzbAOODIiXg78DDgXQNLhFH+4HwGsBD6XHmlA8f25BliWlpWp/Exga0S8BPgscGE61yLgPOCVwDHAeZIWpmMupHg0sgzYms7RNl0LQsCbgNsiYjNARGyOiMGIqAFfoLgwZmZTRtQylgk2XiPiuxExkF7eBCxJ6ycBl0fEjoh4ANgAHCPpQGu0LfcAAAlKSURBVGBBRNwYRde8LwEn1x1zaVq/EjhOxbOPE4F1EbEl3XFaB6xM296Q9iUdO3SutuhmEDqdultx6UIOOQW4s+M1MjNrYih3XNkCHCppfd2yZpxv+W7gmrR+EPCLum0bU9lBaX10+YhjUmB7Gljc5FyLgafqgmD9udqiK8+EJM0HjgfeW1f8PyQdBQTw4KhtZmZd10IW7fsj4rSx9pH0PeAFDTZ9PCKuSvt8HBgAvjJ0WKO3a1I+nmOanastuhKEImIbRcStLzujG3UxM8s19Eyo+T7l39kR8cZm21OngbcAx8Xw6NeNwMF1uy0BHk3lSxqU1x+zUVIfsDewJZW/ftQxPwCeAPaR1JdaQ/Xnaotu3o4zM+stMdwaarZMhKSVwEeBt6Y/2IdcDaySNCf1YFsG3BIRm4BnJR2bnum8E7iq7pjVaf3twPdTULsWOEHSwtQh4QTg2rTt+rQv6dihc7VFt7tom5n1jKxxQhPPov13wBxgXRo/eVNEvC8i7pJ0BXA3xW26syJiqH//+xnuon0Nw8+RLgK+LGkDRQtoVfocWyR9Crg17ffJiBiaBfCjwOWSPg38OJ2jbRyEzMwyRc5UDhPsHpe6U4+17Xzg/Abl64EjG5RvB04d41wXAxc3KL+fDvZOdhAyM8tUq0X5pHYTHyc0ozgImZlliows2p7UrjU9HYTmLJjHIW/8jXEfr+rU7pfR7C+q+rpP9b+8cq/zVP8cU82Op57btT7YP9Bkz6lrcPuOEa93PrNtjD3z1P+uzV2898iNd9zKREWtRtQGm+/jWe1a0tNByMyskyIyglDJdhvJQcjMLFdOS6jmllArHITMzDLVBgeo9TdPKB0DbUs4PS05CJmZ5cq6HeeWUCschMzMMmV1TPAzoZY4CJmZZcrqmODecS1xEDIzy1Qb7Gew5JlPbaC/6XYbyUHIzCyTu2hPPgchM7NcWYNVHYRa0dtBSKI6d3a3a9E+syZ5v6luCnyOyezZpEp7M3Ls/eLhCS8rs3vzv3Klb+QPXfP2mNgJ6wLE4PPPTuxcDUStRgyWBCFn/mhJx/PWSDpM0u11yzOSPiRpkaR1ku5L/y7sdN3MzJqJwQFqAzubL4N+JtSKjgehiLg3Io6KiKOA5cA24BvAOcB1EbEMuC69NjObMoaeCZUtlq/bGTyPA34eEQ8BJwGXpvJLgZO7ViszswaKSe1qTRc8WLUl3b6RvAq4LK0fkKapJSI2Sdq/e9UyM2sgo6XjcUKt6VoQkjQbeCtwbovHrQHWACxZuFcbamZm1tjg4ADyOKFJ1c3bcW8CbouIzen1ZkkHAqR/H2t0UESsjYgVEbFi8Z7zO1RVMzOKLtqDg80XPxNqSTeD0OkM34oDuBpYndZXA1d1vEZmZk3kdUzw7bhWdOV2nKT5wPHAe+uKLwCukHQm8DBwajfqZmY2lpyMCbgl1JKuBKGI2AYsHlX2JEVvOTOzqah/cPMdDG59oPlegzsB/GAokyKi23UYN0mPAw8B+wJPdLk6U4GvQ8HXoeDrMGxfYI+I2G+8J0h3cF6WufuD6Q9rK9HTQWiIpPURsaLb9eg2X4eCr0PB12GYr8XU1e3BqmZmNoM5CJmZWddMlyC0ttsVmCJ8HQq+DgVfh2G+FlPUtHgmZGZmvWm6tITMzKwHOQiZmVnXTIsgJOlPJIWkfevKzpW0QdK9kk7sZv3aTdJfSPp3ST+V9A1J+9RtmzHXAUDSyvRZN0iaMXNSSTpY0vWS7pF0l6SzU/mMnCxSUlXSjyX9S3o9I69DL+j5ICTpYIoUQA/XlR1OMU3EEcBK4HOSqt2pYUesA46MiJcDPyNlJp9p1yF9tr+nSI57OHB6ugYzwQDwkYh4GXAscFb67DN1ssizgXvqXs/U6zDl9XwQAj4L/BlQ38PiJODyiNgREQ8AG4BjulG5ToiI70bEQHp5E7Akrc+o60Dx2TZExP0RsRO4nOIaTHsRsSkibkvrz1J8AR/EDJwsUtIS4HeBL9YVz7jr0Ct6OghJeivwSET8ZNSmg4Bf1L3emMpmgncD16T1mXYdZtrnbUjSIcArgJsZNVkkMBMmi/xrij9M69NZz8Tr0BO6PbNqKUnfA17QYNPHgY8BJzQ6rEFZT/dFb3YdIuKqtM/HKW7LfGXosAb79/R1KDHTPu9uJO0JfA34UEQ8IzW6JNOXpLcAj0XEjyS9vtv1sXJTPghFxBsblUv6DWAp8JP0H20JcJukYyj+Aj64bvclwKNtrmpbjXUdhkhaDbwFOC6GB39Nu+tQYqZ93hEkzaIIQF+JiK+n4s2SDoyITc0mi5xGXgO8VdKbgbnAAkn/xMy7Dj2jZ2/HRcQdEbF/RBwSEYdQfAEdHRG/pJggb5WkOZKWAsuAW7pY3baStBL4KPDWNE3GkBl1HYBbgWWSlqbp41dRXINpT8VfYhcB90TEX9VtmlGTRUbEuRGxJH0nrAK+HxHvYIZdh14y5VtC4xERd0m6Arib4vbUWRExnWea+jtgDrAutQpvioj3zbTrEBEDkj4AXAtUgYsj4q4uV6tTXgOcAdwh6fZU9jE8WeQQX4cpyml7zMysa3r2dpyZmfU+ByEzM+saByEzM+saByEzM+saByEzM+saByHrSZI+IelPul0PM5sYByEzM+saByHrGZI+nuYK+h5wWCp7j6RbJf1E0tckzZe0l6QHUhobJC2Q9KCkWZI+KOnuNPfS5V39QGbmIGS9QdJyijQsrwDeBvxW2vT1iPitiPhNiukLzkxTGfyAIp0/6bivRUQ/xTwyr0hzL72vgx/BzBpwELJe8VrgGxGxLSKeYTgn3JGS/q+kO4A/oJjAD4q5ZN6V1t8F/K+0/lPgK5LeQZHKyMy6yEHIekmjHFOXAB+IiN8A/itF5mQi4t+AQyT9DlCNiDvT/r9LMfvqcuBHkqZl/kSzXuEgZL3iBuAUSfMk7QX8XirfC9iUnv/8wahjvgRcRmoFSaoAB0fE9RSTnu0D7NmJyptZY05gaj0jTdr3TuAhiqk77gaepwgoDwF3AHtFxB+m/V8APAAcGBFPpUB1PbA3xQR4/xQRF3T6c5jZMAchm7YkvR04KSLO6HZdzKwx3w+3aUnS3wJvAt7c7bqY2djcEjIzs65xxwQzM+saByEzM+saByEzM+saByEzM+saByEzM+ua/w/5Bvp4r/cgCwAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "sea_gwd.sel(alt = slice(70,110)).plot(x='days',robust=True)" ] }, { "cell_type": "code", "execution_count": null, "id": "6e2915f0", "metadata": { "papermill": { "duration": 0.012963, "end_time": "2021-06-07T20:43:23.017994", "exception": false, "start_time": "2021-06-07T20:43:23.005031", "status": "completed" }, "tags": [] }, "outputs": [], "source": [] } ], "metadata": { "celltoolbar": "Tags", "kernelspec": { "display_name": "Python 3", "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.6.10" }, "papermill": { "default_parameters": {}, "duration": 19.637191, "end_time": "2021-06-07T20:43:23.441666", "environment_variables": {}, "exception": null, "input_path": "GWD_analysis.ipynb", "output_path": "GWD_analysis_Andenes.ipynb", "parameters": { "station": "Andenes" }, "start_time": "2021-06-07T20:43:03.804475", "version": "2.2.0" } }, "nbformat": 4, "nbformat_minor": 5 }