{ "cells": [ { "cell_type": "markdown", "id": "885c8ea0", "metadata": { "papermill": { "duration": 0.018911, "end_time": "2021-06-08T09:49:48.666285", "exception": false, "start_time": "2021-06-08T09:49:48.647374", "status": "completed" }, "tags": [] }, "source": [ "# Analysis of the Gravity wave drag" ] }, { "cell_type": "code", "execution_count": 1, "id": "5fd6fb79", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:48.688646Z", "iopub.status.busy": "2021-06-08T09:49:48.688202Z", "iopub.status.idle": "2021-06-08T09:49:49.777567Z", "shell.execute_reply": "2021-06-08T09:49:49.778462Z" }, "papermill": { "duration": 1.104506, "end_time": "2021-06-08T09:49:49.778775", "exception": false, "start_time": "2021-06-08T09:49:48.674269", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "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\n", "\n", "from definitions import * # imports all functions from definitions.py" ] }, { "cell_type": "code", "execution_count": 2, "id": "152f2dcb", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.809345Z", "iopub.status.busy": "2021-06-08T09:49:49.808920Z", "iopub.status.idle": "2021-06-08T09:49:49.810668Z", "shell.execute_reply": "2021-06-08T09:49:49.810963Z" }, "papermill": { "duration": 0.014605, "end_time": "2021-06-08T09:49:49.811063", "exception": false, "start_time": "2021-06-08T09:49:49.796458", "status": "completed" }, "tags": [ "parameters" ] }, "outputs": [], "source": [ "station = \"Leipzig\"" ] }, { "cell_type": "code", "execution_count": 3, "id": "2da7832b", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.831826Z", "iopub.status.busy": "2021-06-08T09:49:49.831400Z", "iopub.status.idle": "2021-06-08T09:49:49.833314Z", "shell.execute_reply": "2021-06-08T09:49:49.833615Z" }, "papermill": { "duration": 0.014082, "end_time": "2021-06-08T09:49:49.833710", "exception": false, "start_time": "2021-06-08T09:49:49.819628", "status": "completed" }, "tags": [ "injected-parameters" ] }, "outputs": [], "source": [ "# Parameters\n", "station = \"Davis\"" ] }, { "cell_type": "code", "execution_count": 4, "id": "6e2346c1", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.856826Z", "iopub.status.busy": "2021-06-08T09:49:49.856280Z", "iopub.status.idle": "2021-06-08T09:49:49.857721Z", "shell.execute_reply": "2021-06-08T09:49:49.858053Z" }, "papermill": { "duration": 0.01623, "end_time": "2021-06-08T09:49:49.858160", "exception": false, "start_time": "2021-06-08T09:49:49.841930", "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 == \"Sodankyla_Kiruna\":\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": 5, "id": "b803fc91", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.879234Z", "iopub.status.busy": "2021-06-08T09:49:49.878817Z", "iopub.status.idle": "2021-06-08T09:49:49.903829Z", "shell.execute_reply": "2021-06-08T09:49:49.903278Z" }, "papermill": { "duration": 0.038018, "end_time": "2021-06-08T09:49:49.903968", "exception": false, "start_time": "2021-06-08T09:49:49.865950", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "ds = xr.open_dataset(dir)" ] }, { "cell_type": "code", "execution_count": 6, "id": "c2cbc6bf", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.944241Z", "iopub.status.busy": "2021-06-08T09:49:49.941897Z", "iopub.status.idle": "2021-06-08T09:49:49.946593Z", "shell.execute_reply": "2021-06-08T09:49:49.947149Z" }, "papermill": { "duration": 0.029379, "end_time": "2021-06-08T09:49:49.947319", "exception": false, "start_time": "2021-06-08T09:49:49.917940", "status": "completed" }, "tags": [] }, "outputs": [ { "data": { "text/html": [ "
[140256 values with dtype=int64]
array([ 70, 72, 74, 76, 78, 80, 82, 84, 86, 88, 90, 92, 94, 96,\n",
" 98, 100, 102, 104, 106, 108, 110, 112, 114, 116, 118, 120])array(['2005-01-01T00:00:00.000000000', '2005-01-01T01:00:00.000000000',\n",
" '2005-01-01T02:00:00.000000000', ..., '2020-12-31T21:00:00.000000000',\n",
" '2020-12-31T22:00:00.000000000', '2020-12-31T23:00:00.000000000'],\n",
" dtype='datetime64[ns]')[3646656 values with dtype=float64]