{ "cells": [ { "cell_type": "markdown", "id": "3241e248", "metadata": { "papermill": { "duration": 0.021183, "end_time": "2021-06-08T09:49:48.680834", "exception": false, "start_time": "2021-06-08T09:49:48.659651", "status": "completed" }, "tags": [] }, "source": [ "# Analysis of the Gravity wave drag" ] }, { "cell_type": "code", "execution_count": 1, "id": "87f41c6e", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:48.703595Z", "iopub.status.busy": "2021-06-08T09:49:48.703131Z", "iopub.status.idle": "2021-06-08T09:49:49.815199Z", "shell.execute_reply": "2021-06-08T09:49:49.814772Z" }, "papermill": { "duration": 1.126452, "end_time": "2021-06-08T09:49:49.815314", "exception": false, "start_time": "2021-06-08T09:49:48.688862", "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": "45590ae2", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.835463Z", "iopub.status.busy": "2021-06-08T09:49:49.835046Z", "iopub.status.idle": "2021-06-08T09:49:49.836843Z", "shell.execute_reply": "2021-06-08T09:49:49.837171Z" }, "papermill": { "duration": 0.013598, "end_time": "2021-06-08T09:49:49.837273", "exception": false, "start_time": "2021-06-08T09:49:49.823675", "status": "completed" }, "tags": [ "parameters" ] }, "outputs": [], "source": [ "station = \"Leipzig\"" ] }, { "cell_type": "code", "execution_count": 3, "id": "ac49c265", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.856288Z", "iopub.status.busy": "2021-06-08T09:49:49.855866Z", "iopub.status.idle": "2021-06-08T09:49:49.857601Z", "shell.execute_reply": "2021-06-08T09:49:49.857927Z" }, "papermill": { "duration": 0.012977, "end_time": "2021-06-08T09:49:49.858028", "exception": false, "start_time": "2021-06-08T09:49:49.845051", "status": "completed" }, "tags": [ "injected-parameters" ] }, "outputs": [], "source": [ "# Parameters\n", "station = \"RioGrande\"" ] }, { "cell_type": "code", "execution_count": 4, "id": "31ea4aed", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.883155Z", "iopub.status.busy": "2021-06-08T09:49:49.882713Z", "iopub.status.idle": "2021-06-08T09:49:49.884719Z", "shell.execute_reply": "2021-06-08T09:49:49.884335Z" }, "papermill": { "duration": 0.019318, "end_time": "2021-06-08T09:49:49.884819", "exception": false, "start_time": "2021-06-08T09:49:49.865501", "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": "aaf201ee", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.906553Z", "iopub.status.busy": "2021-06-08T09:49:49.906141Z", "iopub.status.idle": "2021-06-08T09:49:49.933478Z", "shell.execute_reply": "2021-06-08T09:49:49.934180Z" }, "papermill": { "duration": 0.040691, "end_time": "2021-06-08T09:49:49.934389", "exception": false, "start_time": "2021-06-08T09:49:49.893698", "status": "completed" }, "tags": [] }, "outputs": [], "source": [ "ds = xr.open_dataset(dir)" ] }, { "cell_type": "code", "execution_count": 6, "id": "19c4f8e1", "metadata": { "execution": { "iopub.execute_input": "2021-06-08T09:49:49.970687Z", "iopub.status.busy": "2021-06-08T09:49:49.967214Z", "iopub.status.idle": "2021-06-08T09:49:49.973694Z", "shell.execute_reply": "2021-06-08T09:49:49.974203Z" }, "papermill": { "duration": 0.024004, "end_time": "2021-06-08T09:49:49.974365", "exception": false, "start_time": "2021-06-08T09:49:49.950361", "status": "completed" }, "tags": [] }, "outputs": [ { "data": { "text/html": [ "
[113976 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(['2008-02-01T00:00:00.000000000', '2008-02-01T01:00:00.000000000',\n",
" '2008-02-01T02:00:00.000000000', ..., '2021-01-31T21:00:00.000000000',\n",
" '2021-01-31T22:00:00.000000000', '2021-01-31T23:00:00.000000000'],\n",
" dtype='datetime64[ns]')[2963376 values with dtype=float64]