{ "cells": [ { "cell_type": "markdown", "id": "645bc418-a8d8-4ea2-9fa0-890e29bcca48", "metadata": {}, "source": [ "## Maximum vertial wind map\n", "\n", "Use data provided by the Unified Model's .pe0 file to geographically locate the site of the largest vertial winds within a model domain. Animate this data as the model progresses. Useful for diagnosing model instabilities." ] }, { "cell_type": "code", "execution_count": 1, "id": "691460e8-3f7f-4ead-abd7-6aaefe11144c", "metadata": { "tags": [] }, "outputs": [], "source": [ "import xarray as xr\n", "import datetime\n", "import cartopy.crs as ccrs\n", "from matplotlib import pyplot as plt, animation" ] }, { "cell_type": "code", "execution_count": 2, "id": "6cfd36d7-e5b2-402c-90af-ffc63086ad14", "metadata": { "tags": [] }, "outputs": [], "source": [ "ts1=30\n", "starttime=datetime.datetime(year=2022,month=2,day=21,hour=18)\n", "td1=datetime.timedelta(seconds=ts1)\n", "\n", "ts2=75\n", "td2=datetime.timedelta(seconds=ts2)\n", "\n", "tdg=datetime.timedelta(seconds=15)" ] }, { "cell_type": "code", "execution_count": 3, "id": "0847e72c-f217-4667-b901-c865c8f58522", "metadata": { "tags": [] }, "outputs": [], "source": [ "ds=xr.open_dataset(f\"/g/data/hh5/tmp/WACI-Hackathon-2023/AUS2200/data/surf/orog_1hr_{(starttime + datetime.timedelta(hours=7)).strftime('%Y%m%d')}_20220307.nc\")" ] }, { "cell_type": "code", "execution_count": 4, "id": "f66e93c6-7fbc-42d4-b990-b35a02ccdb44", "metadata": { "tags": [] }, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n",
"Dimensions: (latitude: 2120, longitude: 2600, time: 336)\n",
"Coordinates:\n",
" * latitude (latitude) float32 -48.79 -48.77 -48.75 ... -6.871 -6.852 -6.832\n",
" * longitude (longitude) float32 107.5 107.5 107.6 107.6 ... 158.9 159.0 159.0\n",
" * time (time) datetime64[ns] 2022-02-22T00:00:30 ... 2022-03-08\n",
"Data variables:\n",
" orog (time, latitude, longitude) float32 ...\n",
"Attributes: (12/25)\n",
" Conventions: CF-1.6, ACDD-1.3\n",
" title: AUS2200 UM amip simulation over Australian region f...\n",
" summary: A limited area model study of the entire Australian...\n",
" source: UM v12.2 ECMWF ERA5 forecast fields used as forcings\n",
" license: http://creativecommons.org/licenses/by/4.0/\n",
" product_version: 0.1\n",
" ... ...\n",
" geospatial_lon_min: 107.5101\n",
" geospatial_lon_max: 0.0198\n",
" contributor_name: Sam Green, Paola Petrelli\n",
" contributor_email: sam.green@unsw.edu.au, paola.petrelli@utas.edu.au\n",
" contributor_url: https://orcid.org/0000-0003-1129-4676, https://orci...\n",
" contributor_role: post-processing, post-processing<xarray.DataArray 'longitude' ()>\n",
"array(158.9802, dtype=float32)\n",
"Coordinates:\n",
" longitude float32 159.0\n",
"Attributes:\n",
" standard_name: longitude\n",
" units: degrees