{ "cells": [ { "cell_type": "markdown", "id": "f9fae6f4", "metadata": {}, "source": [ "# OMB Analysis" ] }, { "cell_type": "code", "execution_count": 1, "id": "31822224", "metadata": {}, "outputs": [], "source": [ "%config InlineBackend.figure_format = 'svg'\n", "\n", "import xarray as xr\n", "import matplotlib.pyplot as plt\n", "import numpy as np\n", "import pandas as pd\n", "\n", "from openMWR.models import Model\n", "from openMWR.omb import Omb_Analysis, get_input_noise\n", "\n", "DATA_DIR = '../../../../data'" ] }, { "cell_type": "code", "execution_count": 2, "id": "d96664eb", "metadata": {}, "outputs": [], "source": [ "site = 'munich_G5'\n", "\n", "data_source = 'radiosonde'\n", "\n", "model = Model.load('NN_opt', site=site, data_dir=DATA_DIR)\n", "\n", "input_noise = get_input_noise(model)" ] }, { "cell_type": "code", "execution_count": 3, "id": "3a7c7053", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.DataArray 'TB' (time: 770, frq: 14)> Size: 86kB\n",
"array([[-1.40297411, -1.26841401, -1.00572544, ..., 0.93969751,\n",
" 1.066543 , 1.09720991],\n",
" [-2.58776186, -2.37139069, -2.01889836, ..., 0.90047798,\n",
" 1.04439825, 1.07300073],\n",
" [ 0.39260703, 0.67978445, 1.03742362, ..., -0.23038844,\n",
" -0.22769913, -0.23773495],\n",
" ...,\n",
" [-1.37367536, -0.92333656, -0.9081708 , ..., 0.33742529,\n",
" 0.38558368, 0.32665472],\n",
" [-1.31882046, -0.74580405, -0.80185115, ..., 0.90980798,\n",
" 1.08574595, 1.10787999],\n",
" [-1.84070434, -1.36453684, -1.31819138, ..., 0.88988553,\n",
" 1.08184554, 1.1346612 ]], shape=(770, 14))\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 6kB 2024-11-01T10:45:00 ... 2026-01-19T22:...\n",
" synop (time) datetime64[ns] 6kB 2024-11-01T12:00:00 ... 2026-01-20\n",
" * frq (frq) float64 112B 22.24 23.04 23.84 25.44 ... 56.66 57.3 58.0<xarray.DataArray 'TB' (time: 344, frq: 14)> Size: 39kB\n",
"array([[-1.40297411, -1.26841401, -1.00572544, ..., 0.93969751,\n",
" 1.066543 , 1.09720991],\n",
" [-2.58776186, -2.37139069, -2.01889836, ..., 0.90047798,\n",
" 1.04439825, 1.07300073],\n",
" [-0.54014704, -0.19074624, -0.45763188, ..., 0.27013049,\n",
" 0.34309084, 0.25480603],\n",
" ...,\n",
" [-1.89182856, -1.39838717, -1.25472159, ..., 1.03117385,\n",
" 1.25703493, 1.28493482],\n",
" [-1.31882046, -0.74580405, -0.80185115, ..., 0.90980798,\n",
" 1.08574595, 1.10787999],\n",
" [-1.84070434, -1.36453684, -1.31819138, ..., 0.88988553,\n",
" 1.08184554, 1.1346612 ]], shape=(344, 14))\n",
"Coordinates:\n",
" * time (time) datetime64[ns] 3kB 2024-11-01T10:45:00 ... 2026-01-19T22:...\n",
" synop (time) datetime64[ns] 3kB 2024-11-01T12:00:00 ... 2026-01-20\n",
" * frq (frq) float64 112B 22.24 23.04 23.84 25.44 ... 56.66 57.3 58.0<xarray.DataArray (frq: 14)> Size: 112B\n",
"array([0.33823873, 0.2018027 , 0.18984346, 0.09704575, 0.13247507,\n",
" 0.27405771, 0.53062621, 0.26204186, 0.24931903, 0.50929831,\n",
" 0.41307335, 0.11798895, 0.08452118, 0.13577296])\n",
"Coordinates:\n",
" * frq (frq) float64 112B 22.24 23.04 23.84 25.44 ... 56.66 57.3 58.0