diff --git a/examples/stofs_2d.ipynb b/examples/stofs_2d.ipynb new file mode 100644 index 0000000..e420b8b --- /dev/null +++ b/examples/stofs_2d.ipynb @@ -0,0 +1,2273 @@ +{ + "cells": [ + { + "cell_type": "markdown", + "metadata": {}, + "source": [ + "# Example: Subset STOFS dataset given a polygon" + ] + }, + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + "
<xarray.Dataset> Size: 810MB\n", + "Dimensions: (time: 3, node: 12785004, nele: 24875336, nvertex: 3,\n", + " mesh: 1, nbou: 262, nvel: 12421)\n", + "Coordinates:\n", + " * time (time) datetime64[ns] 24B 2023-08-19 ... 2023-08-26T12:...\n", + " x (node) float64 102MB dask.array<chunksize=(511401,), meta=np.ndarray>\n", + " y (node) float64 102MB dask.array<chunksize=(511401,), meta=np.ndarray>\n", + "Dimensions without coordinates: node, nele, nvertex, mesh, nbou, nvel\n", + "Data variables:\n", + " element (nele, nvertex) int32 299MB dask.array<chunksize=(2763927, 1), meta=np.ndarray>\n", + " adcirc_mesh (mesh) int32 4B dask.array<chunksize=(1,), meta=np.ndarray>\n", + " nvell (nbou) int32 1kB dask.array<chunksize=(262,), meta=np.ndarray>\n", + " max_nvell int32 4B ...\n", + " ibtype (nbou) int32 1kB dask.array<chunksize=(262,), meta=np.ndarray>\n", + " nbvv (nvel) int32 50kB dask.array<chunksize=(12421,), meta=np.ndarray>\n", + " depth (node) float64 102MB dask.array<chunksize=(511401,), meta=np.ndarray>\n", + " wind_max (node) float64 102MB dask.array<chunksize=(511401,), meta=np.ndarray>\n", + " time_of_wind_max (node) float64 102MB dask.array<chunksize=(511401,), meta=np.ndarray>\n", + "Attributes: (12/49)\n", + " _FillValue: -99999.0\n", + " model: ADCIRC\n", + " version: 56.dev\n", + " grid_type: Triangular\n", + " description: 2023081900 :-6 hr nowcast and +180 hr forecast ! 32 C...\n", + " agrid: OceanMesh2D\n", + " ... ...\n", + " sfea0: 45.0\n", + " cf: 0.001\n", + " eslm: -0.2\n", + " cori: 0.0\n", + " ntif: 8\n", + " nbfr: 0
<xarray.Dataset> Size: 23MB\n", + "Dimensions: (time: 3, node: 356230, nele: 686355, nvertex: 3,\n", + " mesh: 1, nbou: 262, nvel: 12421)\n", + "Coordinates:\n", + " * time (time) datetime64[ns] 24B 2023-08-19 ... 2023-08-26T12:...\n", + " x (node) float64 3MB dask.array<chunksize=(79,), meta=np.ndarray>\n", + " y (node) float64 3MB dask.array<chunksize=(79,), meta=np.ndarray>\n", + "Dimensions without coordinates: node, nele, nvertex, mesh, nbou, nvel\n", + "Data variables:\n", + " element (nele, nvertex) int32 8MB dask.array<chunksize=(44090, 1), meta=np.ndarray>\n", + " adcirc_mesh (mesh) int32 4B dask.array<chunksize=(1,), meta=np.ndarray>\n", + " nvell (nbou) int32 1kB dask.array<chunksize=(262,), meta=np.ndarray>\n", + " max_nvell int32 4B 1772\n", + " ibtype (nbou) int32 1kB dask.array<chunksize=(262,), meta=np.ndarray>\n", + " nbvv (nvel) int32 50kB dask.array<chunksize=(12421,), meta=np.ndarray>\n", + " depth (node) float64 3MB dask.array<chunksize=(79,), meta=np.ndarray>\n", + " wind_max (node) float64 3MB dask.array<chunksize=(79,), meta=np.ndarray>\n", + " time_of_wind_max (node) float64 3MB dask.array<chunksize=(79,), meta=np.ndarray>\n", + "Attributes: (12/49)\n", + " _FillValue: -99999.0\n", + " model: ADCIRC\n", + " version: 56.dev\n", + " grid_type: Triangular\n", + " description: 2023081900 :-6 hr nowcast and +180 hr forecast ! 32 C...\n", + " agrid: OceanMesh2D\n", + " ... ...\n", + " sfea0: 45.0\n", + " cf: 0.001\n", + " eslm: -0.2\n", + " cori: 0.0\n", + " ntif: 8\n", + " nbfr: 0