{ "cells": [ { "cell_type": "markdown", "id": "7866ac9a", "metadata": {}, "source": [ "# FVCOM 3D" ] }, { "cell_type": "markdown", "id": "11f260c9", "metadata": {}, "source": [ "[](https://colab.research.google.com/github/asascience-open/xarray-subset-grid/blob/main/docs/examples/fvcom_3d.ipynb)" ] }, { "cell_type": "code", "execution_count": null, "id": "66ea390e", "metadata": {}, "outputs": [], "source": [ "# Install required libraries\n", "# !pip install -q xarray_subset_grid@git+https://github.com/asascience-open/xarray-subset-grid.git\n", "# !pip install -q s3fs cftime xarray cf-xarray fsspec dask h5netcdf" ] }, { "cell_type": "code", "execution_count": 1, "id": "29314777-4fff-4f63-ac75-65620ba193e0", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset> Size: 113GB\n",
"Dimensions: (time: 255, four: 4, nele: 569405, node: 303714,\n",
" three: 3, maxnode: 10, maxelem: 8, siglev: 41,\n",
" siglay: 40)\n",
"Coordinates:\n",
" lat (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" latc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" lon (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" lonc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" siglay (siglay, node) float32 49MB dask.array<chunksize=(40, 303714), meta=np.ndarray>\n",
" siglev (siglev, node) float32 50MB dask.array<chunksize=(41, 303714), meta=np.ndarray>\n",
" * time (time) datetime64[ns] 2kB 2024-06-18T21:00:00 ... 202...\n",
"Dimensions without coordinates: four, nele, node, three, maxnode, maxelem\n",
"Data variables: (12/42)\n",
" Itime (time) datetime64[ns] 2kB dask.array<chunksize=(1,), meta=np.ndarray>\n",
" Times (time) |S26 7kB dask.array<chunksize=(1,), meta=np.ndarray>\n",
" a1u (four, nele) float32 9MB dask.array<chunksize=(4, 569405), meta=np.ndarray>\n",
" a2u (four, nele) float32 9MB dask.array<chunksize=(4, 569405), meta=np.ndarray>\n",
" art1 (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" art2 (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" ... ...\n",
" ww (time, siglay, nele) float32 23GB dask.array<chunksize=(1, 40, 569405), meta=np.ndarray>\n",
" x (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" xc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" y (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" yc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" zeta (time, node) float32 310MB dask.array<chunksize=(1, 303714), meta=np.ndarray>\n",
"Attributes: (12/14)\n",
" Conventions: CF-1.0\n",
" CoordinateProjection: init=nad83:4205\n",
" CoordinateSystem: GeoReferenced\n",
" GroundWater_Forcing: GROUND WATER FORCING IS OFF!\n",
" River_Forcing: THERE ARE 63 RIVERS IN THIS MODEL.\\nRIVER IN...\n",
" Surface_Heat_Forcing: FVCOM variable surface heat forcing file:\\nF...\n",
" ... ...\n",
" Tidal_Forcing: TIDAL ELEVATION FORCING IS OFF!\n",
" history: model started at: 18/06/2024 21:15\n",
" institution: School for Marine Science and Technology\n",
" references: http://fvcom.smast.umassd.edu, http://codfis...\n",
" source: FVCOM_4.3\n",
" title: NGOFS2<xarray.Dataset> Size: 113GB\n",
"Dimensions: (time: 255, four: 4, nele: 569405, node: 303714,\n",
" three: 3, maxnode: 10, maxelem: 8, siglev: 41,\n",
" siglay: 40)\n",
"Coordinates:\n",
" lat (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" latc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" lon (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" lonc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" siglay (siglay, node) float32 49MB dask.array<chunksize=(40, 303714), meta=np.ndarray>\n",
" siglev (siglev, node) float32 50MB dask.array<chunksize=(41, 303714), meta=np.ndarray>\n",
" * time (time) datetime64[ns] 2kB 2024-06-18T21:00:00 ... 202...\n",
"Dimensions without coordinates: four, nele, node, three, maxnode, maxelem\n",
"Data variables: (12/42)\n",
" Itime (time) datetime64[ns] 2kB dask.array<chunksize=(1,), meta=np.ndarray>\n",
" Times (time) |S26 7kB dask.array<chunksize=(1,), meta=np.ndarray>\n",
" a1u (four, nele) float32 9MB dask.array<chunksize=(4, 569405), meta=np.ndarray>\n",
" a2u (four, nele) float32 9MB dask.array<chunksize=(4, 569405), meta=np.ndarray>\n",
" art1 (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" art2 (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" ... ...\n",
" ww (time, siglay, nele) float32 23GB dask.array<chunksize=(1, 40, 569405), meta=np.ndarray>\n",
" x (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" xc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" y (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" yc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" zeta (time, node) float32 310MB dask.array<chunksize=(1, 303714), meta=np.ndarray>\n",
"Attributes: (12/14)\n",
" Conventions: CF-1.0\n",
" CoordinateProjection: init=nad83:4205\n",
" CoordinateSystem: GeoReferenced\n",
" GroundWater_Forcing: GROUND WATER FORCING IS OFF!\n",
" River_Forcing: THERE ARE 63 RIVERS IN THIS MODEL.\\nRIVER IN...\n",
" Surface_Heat_Forcing: FVCOM variable surface heat forcing file:\\nF...\n",
" ... ...\n",
" Tidal_Forcing: TIDAL ELEVATION FORCING IS OFF!\n",
" history: model started at: 18/06/2024 21:15\n",
" institution: School for Marine Science and Technology\n",
" references: http://fvcom.smast.umassd.edu, http://codfis...\n",
" source: FVCOM_4.3\n",
" title: NGOFS2<xarray.Dataset> Size: 12GB\n",
"Dimensions: (nele: 569405, node: 303714, three: 3, time: 255, siglay: 40)\n",
"Coordinates:\n",
" lonc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" latc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" lat (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" lon (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" siglay (siglay, node) float32 49MB dask.array<chunksize=(40, 303714), meta=np.ndarray>\n",
" * time (time) datetime64[ns] 2kB 2024-06-18T21:00:00 ... 2024-07-20T...\n",
"Dimensions without coordinates: nele, node, three\n",
"Data variables:\n",
" fvcom_mesh int32 4B ...\n",
" nv (three, nele) int32 7MB dask.array<chunksize=(3, 569405), meta=np.ndarray>\n",
" nbe (three, nele) int32 7MB dask.array<chunksize=(3, 569405), meta=np.ndarray>\n",
" salinity (time, siglay, node) float32 12GB dask.array<chunksize=(1, 40, 303714), meta=np.ndarray>\n",
"Attributes: (12/14)\n",
" Conventions: CF-1.0\n",
" CoordinateProjection: init=nad83:4205\n",
" CoordinateSystem: GeoReferenced\n",
" GroundWater_Forcing: GROUND WATER FORCING IS OFF!\n",
" River_Forcing: THERE ARE 63 RIVERS IN THIS MODEL.\\nRIVER IN...\n",
" Surface_Heat_Forcing: FVCOM variable surface heat forcing file:\\nF...\n",
" ... ...\n",
" Tidal_Forcing: TIDAL ELEVATION FORCING IS OFF!\n",
" history: model started at: 18/06/2024 21:15\n",
" institution: School for Marine Science and Technology\n",
" references: http://fvcom.smast.umassd.edu, http://codfis...\n",
" source: FVCOM_4.3\n",
" title: NGOFS2<xarray.Dataset> Size: 332MB\n",
"Dimensions: (nele: 569405, node: 303714, three: 3, time: 255)\n",
"Coordinates:\n",
" lonc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" latc (nele) float32 2MB dask.array<chunksize=(569405,), meta=np.ndarray>\n",
" lat (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" lon (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" siglay (node) float32 1MB dask.array<chunksize=(303714,), meta=np.ndarray>\n",
" * time (time) datetime64[ns] 2kB 2024-06-18T21:00:00 ... 2024-07-20T...\n",
"Dimensions without coordinates: nele, node, three\n",
"Data variables:\n",
" fvcom_mesh int32 4B ...\n",
" nv (three, nele) int32 7MB dask.array<chunksize=(3, 569405), meta=np.ndarray>\n",
" nbe (three, nele) int32 7MB dask.array<chunksize=(3, 569405), meta=np.ndarray>\n",
" salinity (time, node) float32 310MB dask.array<chunksize=(1, 303714), meta=np.ndarray>\n",
"Attributes: (12/14)\n",
" Conventions: CF-1.0\n",
" CoordinateProjection: init=nad83:4205\n",
" CoordinateSystem: GeoReferenced\n",
" GroundWater_Forcing: GROUND WATER FORCING IS OFF!\n",
" River_Forcing: THERE ARE 63 RIVERS IN THIS MODEL.\\nRIVER IN...\n",
" Surface_Heat_Forcing: FVCOM variable surface heat forcing file:\\nF...\n",
" ... ...\n",
" Tidal_Forcing: TIDAL ELEVATION FORCING IS OFF!\n",
" history: model started at: 18/06/2024 21:15\n",
" institution: School for Marine Science and Technology\n",
" references: http://fvcom.smast.umassd.edu, http://codfis...\n",
" source: FVCOM_4.3\n",
" title: NGOFS2<xarray.Dataset> Size: 59MB\n",
"Dimensions: (nele: 99246, node: 54024, three: 3, time: 255)\n",
"Coordinates:\n",
" lonc (nele) float32 397kB dask.array<chunksize=(99246,), meta=np.ndarray>\n",
" latc (nele) float32 397kB dask.array<chunksize=(99246,), meta=np.ndarray>\n",
" lat (node) float32 216kB dask.array<chunksize=(54024,), meta=np.ndarray>\n",
" lon (node) float32 216kB dask.array<chunksize=(54024,), meta=np.ndarray>\n",
" siglay (node) float32 216kB dask.array<chunksize=(54024,), meta=np.ndarray>\n",
" * time (time) datetime64[ns] 2kB 2024-06-18T21:00:00 ... 2024-07-20T...\n",
"Dimensions without coordinates: nele, node, three\n",
"Data variables:\n",
" fvcom_mesh int32 4B ...\n",
" nv (three, nele) int32 1MB dask.array<chunksize=(3, 99246), meta=np.ndarray>\n",
" nbe (three, nele) int32 1MB dask.array<chunksize=(3, 99246), meta=np.ndarray>\n",
" salinity (time, node) float32 55MB dask.array<chunksize=(1, 54024), meta=np.ndarray>\n",
"Attributes: (12/14)\n",
" Conventions: CF-1.0\n",
" CoordinateProjection: init=nad83:4205\n",
" CoordinateSystem: GeoReferenced\n",
" GroundWater_Forcing: GROUND WATER FORCING IS OFF!\n",
" River_Forcing: THERE ARE 63 RIVERS IN THIS MODEL.\\nRIVER IN...\n",
" Surface_Heat_Forcing: FVCOM variable surface heat forcing file:\\nF...\n",
" ... ...\n",
" Tidal_Forcing: TIDAL ELEVATION FORCING IS OFF!\n",
" history: model started at: 18/06/2024 21:15\n",
" institution: School for Marine Science and Technology\n",
" references: http://fvcom.smast.umassd.edu, http://codfis...\n",
" source: FVCOM_4.3\n",
" title: NGOFS2