{ "cells": [ { "cell_type": "markdown", "id": "8111d1f3-d047-4a0c-9025-1cee328426bf", "metadata": {}, "source": [ "# Notebook to check if the individual glacier interannual variability coincides between glacier models\n", "- related to [PartB_1_annual_variability.ipynb](PartB_1_annual_variability.ipynb), but here we look at individual glaciers to check if the interannual variability is similar between glacier models" ] }, { "cell_type": "markdown", "id": "16312966-e1b0-4fac-92d8-f1f4ab23ceee", "metadata": {}, "source": [ "- attention: area is in m2 * 1e6 in the raw GloGEMFlow per glacier files" ] }, { "cell_type": "code", "execution_count": 1, "id": "b380d17d-14f3-495e-b243-6b677eca02af", "metadata": {}, "outputs": [], "source": [ "save_path = f'/home/www/lschuster/runs_glacierMIP3_oggm_v16/output/RGI11/'\n", "#save_path = os.path.join(OUTPUT_DIR, 'RGI' + rgi_reg)\n", "import os \n", "import numpy as np\n", "import xarray as xr\n", "import pandas as pd\n", "\n", "import matplotlib.pyplot as plt\n", "files = os.listdir(path=save_path)" ] }, { "cell_type": "code", "execution_count": null, "id": "43856e59-b39c-4521-b701-ff8e63602da0", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 2, "id": "d4829ff6-85f9-4eb7-9592-5921f88b8de2", "metadata": {}, "outputs": [], "source": [ "fpath = '/home/www/lschuster/glacierMIP3_analysis/GloGEMflow/glogemflow_individual_glacier_check/'\n", "\n", "test_fs = ['Compagno_rgi11_dyn_1995-2014_mri-esm2-0_hist.nc',\n", " 'Compagno_rgi11_dyn_2021-2040_gfdl-esm4_ssp370.nc',\n", " 'Compagno_rgi11_dyn_2061-2080_mpi-esm1-2-hr_ssp126.nc']\n", "\n", "fpath_sum = '/home/www/lschuster/glacierMIP3_analysis/GloGEMflow/'\n", "test_fs_sum = ['GloGEMflow_rgi11_sum_1995-2014_mri-esm2-0_hist.nc',\n", " 'GloGEMflow_rgi11_sum_2021-2040_gfdl-esm4_ssp370.nc',\n", " 'GloGEMflow_rgi11_sum_2061-2080_mpi-esm1-2-hr_ssp126.nc']\n", "\n", "# PyGEM\n", "fpath_pygem = '/home/www/drounce/per_glacier-filled-norunaway/11/'\n", "test_fs_pygem = ['Rounce_rgi11_glaciers_1995-2014_mri-esm2-0_hist_Batch-1001-2000.nc',\n", " 'Rounce_rgi11_glaciers_2021-2040_gfdl-esm4_ssp370_Batch-1001-2000.nc',\n", " 'Rounce_rgi11_glaciers_2061-2080_mpi-esm1-2-hr_ssp126_Batch-1001-2000.nc']\n", "\n", "# OGGM\n", "oggm_fpath = '/home/www/lschuster/runs_glacierMIP3_oggm_v16/output/RGI11/'\n", "oggm_test_fs= ['mri-esm2-0_historical_1995-2014.nc',\n", " 'gfdl-esm4_ssp370_2021-2040.nc',\n", " 'mpi-esm1-2-hr_ssp126_2061-2080.nc']" ] }, { "cell_type": "code", "execution_count": 3, "id": "0f40f984-775b-4f49-a88e-1d204081318e", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
<xarray.Dataset>\n",
"Dimensions: (y: 5000, x: 1)\n",
"Dimensions without coordinates: y, x\n",
"Data variables:\n",
" area_m2 (y, x) float64 ...\n",
" volume_m3 (y, x) float64 ...\n",
" simulation_year (y, x) float64 ...\n",
"Attributes:\n",
" contributor: Compagno\n",
" contributor_mail: compagno@vaw.baug.ethz.ch\n",
" creation_date: 07-Mar-2022 19:46:26\n",
" reg_region: rgi11\n",
" period: 1995-2014\n",
" gcm: mri-esm2-0\n",
" ssp: hist\n",
" information: no info<xarray.DataArray 'volume_m3' (y: 5000)>\n",
"[5000 values with dtype=float64]\n",
"Dimensions without coordinates: y\n",
"Attributes:\n",
" units: m3\n",
" long_name: regional glacier volume at timestamp<xarray.DataArray 'volume_m3' (y: 5000)>\n",
"[5000 values with dtype=float64]\n",
"Dimensions without coordinates: y\n",
"Attributes:\n",
" units: m3\n",
" long_name: regional glacier volume at timestamp