{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Making a diagnostic\n", "This notebook aims to help you make a new diagnostic function." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from pathlib import Path\n", "from datatree import DataTree\n", "import xarray as xr\n", "import pandas as pd\n", "import valenspy as vp\n", "import numpy as np" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Demo data\n", "The demo data set is already in CF convention so no input conversion and preprocessing is needed." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "True" ] }, "execution_count": 2, "metadata": {}, "output_type": "execute_result" } ], "source": [ "demo_data = vp.demo_data_CF\n", "demo_ds = xr.open_dataset(demo_data)\n", "vp.cf_checks.is_cf_compliant(demo_ds)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "If more intrecate/complex data is needed to test the new diagnostic we recommend using CMIP6 data from Google cloud." ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "True\n" ] }, { "data": { "text/html": [ "
<xarray.Dataset> Size: 1GB\n", "Dimensions: (lat: 256, bnds: 2, lon: 512, time: 1980)\n", "Coordinates:\n", " height float64 8B ...\n", " * lat (lat) float64 2kB -89.46 -88.77 -88.07 ... 88.07 88.77 89.46\n", " lat_bnds (lat, bnds) float64 4kB dask.array<chunksize=(256, 2), meta=np.ndarray>\n", " * lon (lon) float64 4kB 0.0 0.7031 1.406 2.109 ... 357.9 358.6 359.3\n", " lon_bnds (lon, bnds) float64 8kB dask.array<chunksize=(512, 2), meta=np.ndarray>\n", " * time (time) datetime64[ns] 16kB 1850-01-16T12:00:00 ... 2014-12-16T...\n", " time_bnds (time, bnds) datetime64[ns] 32kB dask.array<chunksize=(1980, 2), meta=np.ndarray>\n", "Dimensions without coordinates: bnds\n", "Data variables:\n", " tas (time, lat, lon) float32 1GB dask.array<chunksize=(126, 256, 512), meta=np.ndarray>\n", "Attributes: (12/49)\n", " Conventions: CF-1.7 CMIP-6.2\n", " activity_id: CMIP\n", " branch_method: standard\n", " branch_time_in_child: 0.0\n", " branch_time_in_parent: 29219.0\n", " cmor_version: 3.4.0\n", " ... ...\n", " table_info: Creation Date:(09 May 2019) MD5:cde93...\n", " title: EC-Earth3-Veg output prepared for CMIP6\n", " tracking_id: hdl:21.14100/009fd018-f3fe-4316-ac94-...\n", " variable_id: tas\n", " variant_label: r1i1p1f1\n", " version_id: v20211207
<xarray.DatasetView> Size: 0B\n", "Dimensions: ()\n", "Data variables:\n", " *empty*