{ "cells": [ { "cell_type": "markdown", "id": "ae3913b6-68cc-4d5f-9dc0-4bb7e172846c", "metadata": {}, "source": [ "# Global Water budget" ] }, { "cell_type": "code", "execution_count": 1, "id": "626c70e5-b71e-49a2-989b-e59866a513ea", "metadata": { "tags": [] }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "/work/scratch/env/blazqueza/.conda/envs/lenapy/lib/python3.11/site-packages/esmpy/interface/loadESMF.py:92: VersionWarning: ESMF installation version 8.6.1, ESMPy version 8.6.0\n", " warnings.warn(\"ESMF installation version {}, ESMPy version {}\".format(\n" ] } ], "source": [ "import lenapy\n", "import xarray as xr\n", "import numpy as np\n", "import matplotlib.pyplot as plt" ] }, { "cell_type": "markdown", "id": "6d98be82-9040-4aba-beed-02fb680bd3c2", "metadata": {}, "source": "### Open the ocean mass and its contributors time series" }, { "cell_type": "code", "execution_count": 2, "id": "01adf90a-3452-4f6e-907f-404376a84344", "metadata": { "tags": [] }, "outputs": [], "source": [ "ds=xr.open_dataset('../../data/global_water_budget.nc',engine='lenapyNetcdf')\n", "data=ds.MassChange" ] }, { "cell_type": "markdown", "id": "53237e20-eeb1-4006-bd95-3de1f3e20c59", "metadata": { "tags": [] }, "source": "## Barystatic Sea Level analysis" }, { "cell_type": "code", "execution_count": 3, "id": "2bad0e23-4bca-4229-906b-666c49b4f637", "metadata": { "tags": [] }, "outputs": [ { "data": { "text/html": [ "
<xarray.DataArray 'MassChange' (product_leakage: 1, product_earthquake: 1,\n",
" product_filter: 2, product_center: 5,\n",
" product_GIA: 2, product_C20: 3,\n",
" product_geocenter: 2, coeffs: 6)> Size: 6kB\n",
"array([[[[[[[[-9.58647476e+00, -1.84482885e-01, -3.74030088e-01,\n",
" -2.92337168e-01, -9.37846338e+00, 4.43607432e-03],\n",
" [-1.02043911e+01, -2.36859232e-01, -4.45448962e-01,\n",
" -2.28616066e-02, -1.31640674e+01, 5.86477881e-03]],\n",
"\n",
" [[-9.70688635e+00, 4.22223490e-01, -2.57326598e-01,\n",
" -4.18650785e-01, -9.44069816e+00, 4.43874449e-03],\n",
" [-1.03248027e+01, 3.69847144e-01, -3.28745472e-01,\n",
" -1.49175224e-01, -1.32263021e+01, 5.86744898e-03]],\n",
"\n",
" [[-9.56854614e+00, -5.21655048e-01, -4.38472098e-01,\n",
" -2.45271156e-01, -9.39839179e+00, 4.43790857e-03],\n",
" [-1.01864625e+01, -5.74031395e-01, -5.09890972e-01,\n",
" 2.42044053e-02, -1.31839958e+01, 5.86661306e-03]]],\n",
"\n",
"\n",
" [[[-9.58647476e+00, -1.84482885e-01, -3.74030088e-01,\n",
" -2.92337168e-01, -1.06199305e+01, 4.89885482e-03],\n",
" [-1.02043911e+01, -2.36859232e-01, -4.45448962e-01,\n",
" -2.28616066e-02, -1.44055344e+01, 6.32755931e-03]],\n",
"...\n",
" [[-8.99560348e+00, -4.35502432e-01, -3.23369275e-01,\n",
" -2.40402003e-01, -9.02240595e+00, 4.30502257e-03],\n",
" [-9.61351983e+00, -4.87878778e-01, -3.94788149e-01,\n",
" 2.90735581e-02, -1.28080099e+01, 5.73372707e-03]]],\n",
"\n",
"\n",
" [[[-9.01353210e+00, -9.83302688e-02, -2.58927265e-01,\n",
" -2.87468015e-01, -1.02439446e+01, 4.76596882e-03],\n",
" [-9.63144845e+00, -1.50706615e-01, -3.30346139e-01,\n",
" -1.79924538e-02, -1.40295486e+01, 6.19467331e-03]],\n",
"\n",
" [[-9.13394369e+00, 5.08376106e-01, -1.42223776e-01,\n",
" -4.13781633e-01, -1.03061794e+01, 4.76863899e-03],\n",
" [-9.75186004e+00, 4.55999760e-01, -2.13642650e-01,\n",
" -1.44306071e-01, -1.40917834e+01, 6.19734348e-03]],\n",
"\n",
" [[-8.99560348e+00, -4.35502432e-01, -3.23369275e-01,\n",
" -2.40402003e-01, -1.02638730e+01, 4.76780307e-03],\n",
" [-9.61351983e+00, -4.87878778e-01, -3.94788149e-01,\n",
" 2.90735581e-02, -1.40494770e+01, 6.19650757e-03]]]]]]]])\n",
"Coordinates:\n",
" * product_center (product_center) <U6 120B 'CSR' 'JPL' ... 'CNES' 'GFZ'\n",
" * product_filter (product_filter) <U4 32B 'DDK6' 'DDK3'\n",
" * product_earthquake (product_earthquake) <U7 28B 'Earthq0'\n",
" * product_leakage (product_leakage) <U13 52B 'IceSheetsOnly'\n",
" * product_GIA (product_GIA) <U10 80B 'ICE6GD' 'Caron_2018'\n",
" * product_C20 (product_C20) <U7 84B 'Lemoine' 'Cheng' 'Loomis'\n",
" * product_geocenter (product_geocenter) <U7 56B 'Lemoine' 'Sun'\n",
" maskGlobalBudget <U25 100B 'Ocean'\n",
" * coeffs (coeffs) <U13 312B 'cosAnnual' 'sinAnnual' ... 'order_1'