{ "cells": [ { "cell_type": "code", "execution_count": 1, "id": "22be82f8-f882-4407-a6b3-d468d14c42ce", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import numpy as np\n", "import pandas as pd\n", "from oggm import utils\n", "import oggm" ] }, { "cell_type": "code", "execution_count": 2, "id": "3f9dc5ab-6625-4c45-90c7-e73f4ac1e03e", "metadata": {}, "outputs": [], "source": [ "df_itmix = pd.read_hdf(oggm.utils.get_demo_file('rgi62_itmix_df.h5'))\n", "df_rgi6g = pd.read_hdf(utils.file_downloader('https://cluster.klima.uni-bremen.de/~oggm/rgi/rgi62_stats.h5'))\n", "# level 2 glaciers should be excluded (as was the case for GlacierMIP2)\n", "df_rgi6g = df_rgi6g.loc[df_rgi6g.Connect != 2]\n", "rgi_regs = list(df_rgi6g['O1Region'].unique())\n", "rgi_regs.append('All')\n", "# directly copy/paste from Farinotti et al. (2019) – Table 1 \n", "rgidf_vol_1e3_km3_farinotti={\n", " \"01\": 18.98,\n", " \"02\": 1.06,\n", " \"03\": 28.33,\n", " \"04\": 8.61,\n", " \"05\": 15.69,\n", " \"06\": 3.77,\n", " \"07\": 7.47,\n", " \"08\": 0.30,\n", " \"09\": 14.64,\n", " \"10\": 0.14,\n", " \"11\": 0.13,\n", " \"12\": 0.06,\n", " \"13\": 3.27,\n", " \"14\": 2.87,\n", " \"15\": 0.88,\n", " \"16\": 0.10,\n", " \"17\": 5.34,\n", " \"18\": 0.07,\n", " \"19\": 46.47,\n", " \"All\": 158.17\n", "}" ] }, { "cell_type": "code", "execution_count": 3, "id": "24521596-ee5e-4295-9fec-5de05b6065a7", "metadata": {}, "outputs": [], "source": [ "rgidf_mass_gt = pd.DataFrame()\n", "for rgi_reg in rgi_regs:\n", " if rgi_reg != 'All':\n", " rgidf = df_rgi6g.loc[df_rgi6g.O1Region == rgi_reg]\n", " else:\n", " rgidf = df_rgi6g\n", " rgidf_mass_gt.loc[rgi_reg,'oggm_version'] = (df_itmix.loc[rgidf.index]['vol_itmix_m3'].sum()*900 *1e-12).round()\n", " # 1e3 km³ --> into GT \n", " rgidf_mass_gt_farinotti_reg = rgidf_vol_1e3_km3_farinotti[rgi_reg] * 900\n", " rgidf_mass_gt.loc[rgi_reg,'original_farinotti_2019_version'] = rgidf_mass_gt_farinotti_reg\n", " _perc = (100 * rgidf_mass_gt.loc[rgi_reg,'oggm_version']/rgidf_mass_gt_farinotti_reg).round(2)\n", " rgidf_mass_gt.loc[rgi_reg,'perc_oggm_original_farinotti_2019_version'] = _perc" ] }, { "cell_type": "code", "execution_count": 4, "id": "f99e2429-9d8f-4410-877f-2e1e63719dbf", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", " \n", "
oggm_versionoriginal_farinotti_2019_versionperc_oggm_original_farinotti_2019_version
0117081.017082.099.99
02942.0954.098.74
0325498.025497.0100.00
047750.07749.0100.01
0514123.014121.0100.01
063393.03393.0100.00
076723.06723.0100.00
08269.0270.099.63
0913176.013176.0100.00
10122.0126.096.83
11115.0117.098.29
1257.054.0105.56
132944.02943.0100.03
142579.02583.099.85
15790.0792.099.75
1689.090.098.89
174806.04806.0100.00
1866.063.0104.76
1941820.041823.099.99
All142341.0142353.099.99
\n", "
" ], "text/plain": [ " oggm_version original_farinotti_2019_version \\\n", "01 17081.0 17082.0 \n", "02 942.0 954.0 \n", "03 25498.0 25497.0 \n", "04 7750.0 7749.0 \n", "05 14123.0 14121.0 \n", "06 3393.0 3393.0 \n", "07 6723.0 6723.0 \n", "08 269.0 270.0 \n", "09 13176.0 13176.0 \n", "10 122.0 126.0 \n", "11 115.0 117.0 \n", "12 57.0 54.0 \n", "13 2944.0 2943.0 \n", "14 2579.0 2583.0 \n", "15 790.0 792.0 \n", "16 89.0 90.0 \n", "17 4806.0 4806.0 \n", "18 66.0 63.0 \n", "19 41820.0 41823.0 \n", "All 142341.0 142353.0 \n", "\n", " perc_oggm_original_farinotti_2019_version \n", "01 99.99 \n", "02 98.74 \n", "03 100.00 \n", "04 100.01 \n", "05 100.01 \n", "06 100.00 \n", "07 100.00 \n", "08 99.63 \n", "09 100.00 \n", "10 96.83 \n", "11 98.29 \n", "12 105.56 \n", "13 100.03 \n", "14 99.85 \n", "15 99.75 \n", "16 98.89 \n", "17 100.00 \n", "18 104.76 \n", "19 99.99 \n", "All 99.99 " ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "rgidf_mass_gt" ] }, { "cell_type": "code", "execution_count": null, "id": "24532809-8f7e-48f8-8821-0c85ac2f62d3", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "de976f91-3dca-4947-8aa3-fe876f32e420", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python [conda env:oggm_env_2025]", "language": "python", "name": "conda-env-oggm_env_2025-py" }, "language_info": { "codemirror_mode": { "name": "ipython", "version": 3 }, "file_extension": ".py", "mimetype": "text/x-python", "name": "python", "nbconvert_exporter": "python", "pygments_lexer": "ipython3", "version": "3.11.14" } }, "nbformat": 4, "nbformat_minor": 5 }