{ "cells": [ { "cell_type": "markdown", "id": "e323ca80-b90e-4c8f-9d95-055be2df4d72", "metadata": {}, "source": [ "# Check 2300 GCMs" ] }, { "cell_type": "markdown", "id": "884c89b7-8b81-4188-bd9e-490e872b31a9", "metadata": {}, "source": [ "- todo: check if the two missing GCMs are now inside!!!" ] }, { "cell_type": "code", "execution_count": 1, "id": "771128b9-649c-45c7-b1d1-735f64e9f9f8", "metadata": {}, "outputs": [], "source": [ "import xarray as xr\n", "import os\n", "import numpy as np\n", "import pandas as pd\n", "import glob\n" ] }, { "cell_type": "code", "execution_count": 5, "id": "f0e2d2b7-b1ce-4dfa-b479-62cae8092e02", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "47.112724743783474\n" ] } ], "source": [ "dw5e5_ref = xr.open_dataset('/home/www/oggm/climate/gswp3-w5e5/unflattened/monthly/gswp3-w5e5_obsclim_pr_global_monthly_1901_2019.nc').pr#.groupby('time.year').mean().sel(year=2000).mean()\n", "print(dw5e5_ref.groupby('time.year').mean().max().values *60*60*24)" ] }, { "cell_type": "code", "execution_count": 6, "id": "7982ccb9-0a9b-4f7a-aceb-08f2cf0316d9", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'pr' (time: 1428, lat: 360, lon: 720)>\n",
       "[370137600 values with dtype=float32]\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1901-01-01 1901-02-01 ... 2019-12-01\n",
       "  * lon      (lon) float32 -179.8 -179.2 -178.8 -178.2 ... 178.8 179.2 179.8\n",
       "  * lat      (lat) float32 89.75 89.25 88.75 88.25 ... -88.75 -89.25 -89.75\n",
       "Attributes:\n",
       "    standard_name:  precipitation_flux\n",
       "    long_name:      Precipitation\n",
       "    units:          kg m-2 s-1
" ], "text/plain": [ "\n", "[370137600 values with dtype=float32]\n", "Coordinates:\n", " * time (time) datetime64[ns] 1901-01-01 1901-02-01 ... 2019-12-01\n", " * lon (lon) float32 -179.8 -179.2 -178.8 -178.2 ... 178.8 179.2 179.8\n", " * lat (lat) float32 89.75 89.25 88.75 88.25 ... -88.75 -89.25 -89.75\n", "Attributes:\n", " standard_name: precipitation_flux\n", " long_name: Precipitation\n", " units: kg m-2 s-1" ] }, "execution_count": 6, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dw5e5_ref" ] }, { "cell_type": "code", "execution_count": 6, "id": "539a5007-4739-467c-851d-d00af9878bef", "metadata": {}, "outputs": [], "source": [ "dcc = xr.open_dataset('/home/www/rchizzola/cmip5_done/pr/pr_mon_CanESM2_rcp26_r1i1p1.nc')" ] }, { "cell_type": "code", "execution_count": 7, "id": "f7b9cd47-78e8-4f13-86b9-9bb5a598534e", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'lat' (lat: 64)>\n",
       "array([-87.863801, -85.096529, -82.312915, -79.525609, -76.736902, -73.947518,\n",
       "       -71.157755, -68.367759, -65.57761 , -62.787354, -59.997023, -57.206634,\n",
       "       -54.416202, -51.625736, -48.835243, -46.044729, -43.254197, -40.463651,\n",
       "       -37.673092, -34.882523, -32.091946, -29.301362, -26.510772, -23.720176,\n",
       "       -20.929577, -18.138973, -15.348367, -12.557759,  -9.767148,  -6.976536,\n",
       "        -4.185923,  -1.395309,   1.395309,   4.185923,   6.976536,   9.767148,\n",
       "        12.557759,  15.348367,  18.138973,  20.929577,  23.720176,  26.510772,\n",
       "        29.301362,  32.091946,  34.882523,  37.673092,  40.463651,  43.254197,\n",
       "        46.044729,  48.835243,  51.625736,  54.416202,  57.206634,  59.997023,\n",
       "        62.787354,  65.57761 ,  68.367759,  71.157755,  73.947518,  76.736902,\n",
       "        79.525609,  82.312915,  85.096529,  87.863801])\n",
       "Coordinates:\n",
       "  * lat      (lat) float64 -87.86 -85.1 -82.31 -79.53 ... 79.53 82.31 85.1 87.86\n",
       "Attributes:\n",
       "    bounds:         lat_bnds\n",
       "    units:          degrees_north\n",
       "    axis:           Y\n",
       "    long_name:      latitude\n",
       "    standard_name:  latitude
" ], "text/plain": [ "\n", "array([-87.863801, -85.096529, -82.312915, -79.525609, -76.736902, -73.947518,\n", " -71.157755, -68.367759, -65.57761 , -62.787354, -59.997023, -57.206634,\n", " -54.416202, -51.625736, -48.835243, -46.044729, -43.254197, -40.463651,\n", " -37.673092, -34.882523, -32.091946, -29.301362, -26.510772, -23.720176,\n", " -20.929577, -18.138973, -15.348367, -12.557759, -9.767148, -6.976536,\n", " -4.185923, -1.395309, 1.395309, 4.185923, 6.976536, 9.767148,\n", " 12.557759, 15.348367, 18.138973, 20.929577, 23.720176, 26.510772,\n", " 29.301362, 32.091946, 34.882523, 37.673092, 40.463651, 43.254197,\n", " 46.044729, 48.835243, 51.625736, 54.416202, 57.206634, 59.997023,\n", " 62.787354, 65.57761 , 68.367759, 71.157755, 73.947518, 76.736902,\n", " 79.525609, 82.312915, 85.096529, 87.863801])\n", "Coordinates:\n", " * lat (lat) float64 -87.86 -85.1 -82.31 -79.53 ... 79.53 82.31 85.1 87.86\n", "Attributes:\n", " bounds: lat_bnds\n", " units: degrees_north\n", " axis: Y\n", " long_name: latitude\n", " standard_name: latitude" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dcc.lat" ] }, { "cell_type": "code", "execution_count": 9, "id": "bf2bff0a-2223-4de0-b13c-eb244c394066", "metadata": {}, "outputs": [], "source": [ "dce = xr.open_dataset('/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc')" ] }, { "cell_type": "code", "execution_count": 11, "id": "bb13ffaf-54cc-4390-abf2-2251c2bb6fa3", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'lat' (lat: 192)>\n",
       "array([-90.      , -89.057594, -88.115181, -87.172775, -86.23037 , -85.287956,\n",
       "       -84.345551, -83.403145, -82.460732, -81.518326, -80.575912, -79.633507,\n",
       "       -78.691101, -77.748688, -76.806282, -75.863876, -74.921463, -73.979057,\n",
       "       -73.036652, -72.094238, -71.151833, -70.209427, -69.267014, -68.324608,\n",
       "       -67.382202, -66.439789, -65.497383, -64.554977, -63.612564, -62.670158,\n",
       "       -61.727749, -60.785339, -59.842934, -58.900524, -57.958115, -57.015705,\n",
       "       -56.073299, -55.13089 , -54.18848 , -53.246075, -52.303665, -51.361256,\n",
       "       -50.41885 , -49.47644 , -48.534031, -47.591621, -46.649216, -45.706806,\n",
       "       -44.764397, -43.821991, -42.879581, -41.937172, -40.994766, -40.052357,\n",
       "       -39.109947, -38.167538, -37.225132, -36.282722, -35.340313, -34.397907,\n",
       "       -33.455498, -32.513088, -31.570681, -30.628273, -29.685863, -28.743456,\n",
       "       -27.801046, -26.858639, -25.916231, -24.973822, -24.031414, -23.089005,\n",
       "       -22.146597, -21.204189, -20.26178 , -19.319372, -18.376963, -17.434555,\n",
       "       -16.492147, -15.549738, -14.607329, -13.664922, -12.722513, -11.780105,\n",
       "       -10.837696,  -9.895288,  -8.95288 ,  -8.010471,  -7.068063,  -6.125654,\n",
       "        -5.183246,  -4.240838,  -3.298429,  -2.356021,  -1.413613,  -0.471204,\n",
       "         0.471204,   1.413613,   2.356021,   3.298429,   4.240838,   5.183246,\n",
       "         6.125654,   7.068063,   8.010471,   8.95288 ,   9.895288,  10.837696,\n",
       "        11.780105,  12.722513,  13.664922,  14.607329,  15.549738,  16.492147,\n",
       "        17.434555,  18.376963,  19.319372,  20.26178 ,  21.204189,  22.146597,\n",
       "        23.089005,  24.031414,  24.973822,  25.916231,  26.858639,  27.801046,\n",
       "        28.743456,  29.685863,  30.628273,  31.570681,  32.513088,  33.455498,\n",
       "        34.397907,  35.340313,  36.282722,  37.225132,  38.167538,  39.109947,\n",
       "        40.052357,  40.994766,  41.937172,  42.879581,  43.821991,  44.764397,\n",
       "        45.706806,  46.649216,  47.591621,  48.534031,  49.47644 ,  50.41885 ,\n",
       "        51.361256,  52.303665,  53.246075,  54.18848 ,  55.13089 ,  56.073299,\n",
       "        57.015705,  57.958115,  58.900524,  59.842934,  60.785339,  61.727749,\n",
       "        62.670158,  63.612564,  64.554977,  65.497383,  66.439789,  67.382202,\n",
       "        68.324608,  69.267014,  70.209427,  71.151833,  72.094238,  73.036652,\n",
       "        73.979057,  74.921463,  75.863876,  76.806282,  77.748688,  78.691101,\n",
       "        79.633507,  80.575912,  81.518326,  82.460732,  83.403145,  84.345551,\n",
       "        85.287956,  86.23037 ,  87.172775,  88.115181,  89.057594,  90.      ])\n",
       "Coordinates:\n",
       "  * lat      (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 87.17 88.12 89.06 90.0\n",
       "Attributes:\n",
       "    bounds:         lat_bnds\n",
       "    units:          degrees_north\n",
       "    axis:           Y\n",
       "    long_name:      latitude\n",
       "    standard_name:  latitude
" ], "text/plain": [ "\n", "array([-90. , -89.057594, -88.115181, -87.172775, -86.23037 , -85.287956,\n", " -84.345551, -83.403145, -82.460732, -81.518326, -80.575912, -79.633507,\n", " -78.691101, -77.748688, -76.806282, -75.863876, -74.921463, -73.979057,\n", " -73.036652, -72.094238, -71.151833, -70.209427, -69.267014, -68.324608,\n", " -67.382202, -66.439789, -65.497383, -64.554977, -63.612564, -62.670158,\n", " -61.727749, -60.785339, -59.842934, -58.900524, -57.958115, -57.015705,\n", " -56.073299, -55.13089 , -54.18848 , -53.246075, -52.303665, -51.361256,\n", " -50.41885 , -49.47644 , -48.534031, -47.591621, -46.649216, -45.706806,\n", " -44.764397, -43.821991, -42.879581, -41.937172, -40.994766, -40.052357,\n", " -39.109947, -38.167538, -37.225132, -36.282722, -35.340313, -34.397907,\n", " -33.455498, -32.513088, -31.570681, -30.628273, -29.685863, -28.743456,\n", " -27.801046, -26.858639, -25.916231, -24.973822, -24.031414, -23.089005,\n", " -22.146597, -21.204189, -20.26178 , -19.319372, -18.376963, -17.434555,\n", " -16.492147, -15.549738, -14.607329, -13.664922, -12.722513, -11.780105,\n", " -10.837696, -9.895288, -8.95288 , -8.010471, -7.068063, -6.125654,\n", " -5.183246, -4.240838, -3.298429, -2.356021, -1.413613, -0.471204,\n", " 0.471204, 1.413613, 2.356021, 3.298429, 4.240838, 5.183246,\n", " 6.125654, 7.068063, 8.010471, 8.95288 , 9.895288, 10.837696,\n", " 11.780105, 12.722513, 13.664922, 14.607329, 15.549738, 16.492147,\n", " 17.434555, 18.376963, 19.319372, 20.26178 , 21.204189, 22.146597,\n", " 23.089005, 24.031414, 24.973822, 25.916231, 26.858639, 27.801046,\n", " 28.743456, 29.685863, 30.628273, 31.570681, 32.513088, 33.455498,\n", " 34.397907, 35.340313, 36.282722, 37.225132, 38.167538, 39.109947,\n", " 40.052357, 40.994766, 41.937172, 42.879581, 43.821991, 44.764397,\n", " 45.706806, 46.649216, 47.591621, 48.534031, 49.47644 , 50.41885 ,\n", " 51.361256, 52.303665, 53.246075, 54.18848 , 55.13089 , 56.073299,\n", " 57.015705, 57.958115, 58.900524, 59.842934, 60.785339, 61.727749,\n", " 62.670158, 63.612564, 64.554977, 65.497383, 66.439789, 67.382202,\n", " 68.324608, 69.267014, 70.209427, 71.151833, 72.094238, 73.036652,\n", " 73.979057, 74.921463, 75.863876, 76.806282, 77.748688, 78.691101,\n", " 79.633507, 80.575912, 81.518326, 82.460732, 83.403145, 84.345551,\n", " 85.287956, 86.23037 , 87.172775, 88.115181, 89.057594, 90. ])\n", "Coordinates:\n", " * lat (lat) float64 -90.0 -89.06 -88.12 -87.17 ... 87.17 88.12 89.06 90.0\n", "Attributes:\n", " bounds: lat_bnds\n", " units: degrees_north\n", " axis: Y\n", " long_name: latitude\n", " standard_name: latitude" ] }, "execution_count": 11, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dce.lat" ] }, { "cell_type": "code", "execution_count": 113, "id": "c2558a03-dfbd-47da-b6da-caab081835f0", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:                       (time: 301, rgi_id: 567, month_2d: 12)\n",
       "Coordinates:\n",
       "  * time                          (time) float64 2e+03 2.001e+03 ... 2.3e+03\n",
       "  * rgi_id                        (rgi_id) object 'RGI60-06.00001' ... 'RGI60...\n",
       "    hydro_year                    (time) int64 2000 2001 2002 ... 2298 2299 2300\n",
       "    hydro_month                   (time) int64 4 4 4 4 4 4 4 4 ... 4 4 4 4 4 4 4\n",
       "    calendar_year                 (time) int64 2000 2001 2002 ... 2298 2299 2300\n",
       "    calendar_month                (time) int64 1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1\n",
       "  * month_2d                      (month_2d) int64 1 2 3 4 5 6 7 8 9 10 11 12\n",
       "    calendar_month_2d             (month_2d) int64 1 2 3 4 5 6 7 8 9 10 11 12\n",
       "Data variables: (12/29)\n",
       "    volume                        (time, rgi_id) float32 3.386e+08 ... 0.0\n",
       "    volume_bsl                    (time, rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0\n",
       "    volume_bwl                    (time, rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0\n",
       "    area                          (time, rgi_id) float32 4.894e+06 ... 0.0\n",
       "    length                        (time, rgi_id) float32 2.952e+03 ... 0.0\n",
       "    calving                       (time, rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0\n",
       "    ...                            ...\n",
       "    snowfall_on_glacier_monthly   (time, month_2d, rgi_id) float32 1.358e+09 ...\n",
       "    snow_bucket_monthly           (time, month_2d, rgi_id) float32 0.0 ... 0.0\n",
       "    residual_mb_monthly           (time, month_2d, rgi_id) float32 0.0 ... nan\n",
       "    water_level                   (rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n",
       "    glen_a                        (rgi_id) float32 6.436e-24 ... 6.436e-24\n",
       "    fs                            (rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n",
       "Attributes:\n",
       "    description:    OGGM model output\n",
       "    oggm_version:   1.6.1.dev26+gf8a1745\n",
       "    calendar:       365-day no leap\n",
       "    creation_date:  2023-07-26 11:12:41
" ], "text/plain": [ "\n", "Dimensions: (time: 301, rgi_id: 567, month_2d: 12)\n", "Coordinates:\n", " * time (time) float64 2e+03 2.001e+03 ... 2.3e+03\n", " * rgi_id (rgi_id) object 'RGI60-06.00001' ... 'RGI60...\n", " hydro_year (time) int64 ...\n", " hydro_month (time) int64 ...\n", " calendar_year (time) int64 ...\n", " calendar_month (time) int64 ...\n", " * month_2d (month_2d) int64 1 2 3 4 5 6 7 8 9 10 11 12\n", " calendar_month_2d (month_2d) int64 ...\n", "Data variables: (12/29)\n", " volume (time, rgi_id) float32 3.386e+08 ... 0.0\n", " volume_bsl (time, rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0\n", " volume_bwl (time, rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0\n", " area (time, rgi_id) float32 4.894e+06 ... 0.0\n", " length (time, rgi_id) float32 2.952e+03 ... 0.0\n", " calving (time, rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0\n", " ... ...\n", " snowfall_on_glacier_monthly (time, month_2d, rgi_id) float32 1.358e+09 ...\n", " snow_bucket_monthly (time, month_2d, rgi_id) float32 0.0 ... 0.0\n", " residual_mb_monthly (time, month_2d, rgi_id) float32 0.0 ... nan\n", " water_level (rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n", " glen_a (rgi_id) float32 6.436e-24 ... 6.436e-24\n", " fs (rgi_id) float32 0.0 0.0 0.0 ... 0.0 0.0 0.0\n", "Attributes:\n", " description: OGGM model output\n", " oggm_version: 1.6.1.dev26+gf8a1745\n", " calendar: 365-day no leap\n", " creation_date: 2023-07-26 11:12:41" ] }, "execution_count": 113, "metadata": {}, "output_type": "execute_result" } ], "source": [ "xr.open_dataset('/home/users/lschuster/www_lschuster/runs_oggm_v16/runs_2023.3/output/RGI06/run_hydro_gcm_from_2000_endyr2300_CMIP6_MRI-ESM2-0_ssp585_bc_2000_2019_rgi06_0_1000.nc').dropna(dim='rgi_id', how='all')" ] }, { "cell_type": "code", "execution_count": 107, "id": "50590d5c-a4a2-49f4-93b9-6fb588121f61", "metadata": { "tags": [] }, "outputs": [], "source": [ "gcms_cmip6 = pd.read_csv('/home/www/oggm/cmip6/all_gcm_list.csv', index_col='index')\n", "pd_t = pd.read_html('/home/www/oggm/cmip6/all_gcm_table.html')[0]\n", "folder_p = '/home/www/oggm/cmip6/GCM/'\n", "pd_t = pd_t.loc[pd_t.y1>=2299]\n", "pd_t['path'] = folder_p + pd_t['folder'] + '/' + pd_t['fname'] \n", "# now we select those from the total list \n", "gcms_cmip6.index = gcms_cmip6.path\n", "gcms_cmip6= gcms_cmip6.loc[pd_t.path] \n", "gcms_cmip6 = gcms_cmip6.reset_index(drop=True)\n", "gcms_cmip6 = gcms_cmip6.drop(columns=['Unnamed: 0'])\n", "gcms_cmip6.to_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300.csv')\n", "gcms_cmip6 = pd.read_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300.csv', index_col= 'Unnamed: 0')" ] }, { "cell_type": "code", "execution_count": 109, "id": "434dab65-f5b7-4e86-930a-9d311c41c003", "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", " \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", " \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", "
pathfolderfnamevalidgcmsspreavar
0/home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2...ACCESS-CM2ACCESS-CM2_ssp126_r1i1p1f1_pr.ncTrueACCESS-CM2ssp126r1i1p1f1pr
1/home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2...ACCESS-CM2ACCESS-CM2_ssp126_r1i1p1f1_tas.ncTrueACCESS-CM2ssp126r1i1p1f1tas
2/home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2...ACCESS-CM2ACCESS-CM2_ssp585_r1i1p1f1_pr.ncTrueACCESS-CM2ssp585r1i1p1f1pr
3/home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2...ACCESS-CM2ACCESS-CM2_ssp585_r1i1p1f1_tas.ncTrueACCESS-CM2ssp585r1i1p1f1tas
4/home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-...ACCESS-ESM1-5ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.ncTrueACCESS-ESM1-5ssp126r1i1p1f1pr
5/home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-...ACCESS-ESM1-5ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.ncTrueACCESS-ESM1-5ssp126r1i1p1f1tas
6/home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-...ACCESS-ESM1-5ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.ncTrueACCESS-ESM1-5ssp585r1i1p1f1pr
7/home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-...ACCESS-ESM1-5ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.ncTrueACCESS-ESM1-5ssp585r1i1p1f1tas
8/home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC...CESM2-WACCMCESM2-WACCM_ssp126_r1i1p1f1_pr.ncTrueCESM2-WACCMssp126r1i1p1f1pr
9/home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC...CESM2-WACCMCESM2-WACCM_ssp126_r1i1p1f1_tas.ncTrueCESM2-WACCMssp126r1i1p1f1tas
10/home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC...CESM2-WACCMCESM2-WACCM_ssp534-over_r1i1p1f1_pr.ncTrueCESM2-WACCMssp534-overr1i1p1f1pr
11/home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC...CESM2-WACCMCESM2-WACCM_ssp534-over_r1i1p1f1_tas.ncTrueCESM2-WACCMssp534-overr1i1p1f1tas
12/home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC...CESM2-WACCMCESM2-WACCM_ssp585_r1i1p1f1_pr.ncTrueCESM2-WACCMssp585r1i1p1f1pr
13/home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC...CESM2-WACCMCESM2-WACCM_ssp585_r1i1p1f1_tas.ncTrueCESM2-WACCMssp585r1i1p1f1tas
14/home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp12...CanESM5CanESM5_ssp126_r1i1p1f1_pr.ncTrueCanESM5ssp126r1i1p1f1pr
15/home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp12...CanESM5CanESM5_ssp126_r1i1p1f1_tas.ncTrueCanESM5ssp126r1i1p1f1tas
16/home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp53...CanESM5CanESM5_ssp534-over_r1i1p1f1_pr.ncTrueCanESM5ssp534-overr1i1p1f1pr
17/home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp53...CanESM5CanESM5_ssp534-over_r1i1p1f1_tas.ncTrueCanESM5ssp534-overr1i1p1f1tas
18/home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp58...CanESM5CanESM5_ssp585_r1i1p1f1_pr.ncTrueCanESM5ssp585r1i1p1f1pr
19/home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp58...CanESM5CanESM5_ssp585_r1i1p1f1_tas.ncTrueCanESM5ssp585r1i1p1f1tas
20/home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6...IPSL-CM6A-LRIPSL-CM6A-LR_ssp126_r1i1p1f1_pr.ncTrueIPSL-CM6A-LRssp126r1i1p1f1pr
21/home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6...IPSL-CM6A-LRIPSL-CM6A-LR_ssp126_r1i1p1f1_tas.ncTrueIPSL-CM6A-LRssp126r1i1p1f1tas
22/home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6...IPSL-CM6A-LRIPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.ncTrueIPSL-CM6A-LRssp534-overr1i1p1f1pr
23/home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6...IPSL-CM6A-LRIPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.ncTrueIPSL-CM6A-LRssp534-overr1i1p1f1tas
24/home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6...IPSL-CM6A-LRIPSL-CM6A-LR_ssp585_r1i1p1f1_pr.ncTrueIPSL-CM6A-LRssp585r1i1p1f1pr
25/home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6...IPSL-CM6A-LRIPSL-CM6A-LR_ssp585_r1i1p1f1_tas.ncTrueIPSL-CM6A-LRssp585r1i1p1f1tas
26/home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0...MRI-ESM2-0MRI-ESM2-0_ssp126_r1i1p1f1_pr.ncTrueMRI-ESM2-0ssp126r1i1p1f1pr
27/home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0...MRI-ESM2-0MRI-ESM2-0_ssp126_r1i1p1f1_tas.ncTrueMRI-ESM2-0ssp126r1i1p1f1tas
28/home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0...MRI-ESM2-0MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.ncTrueMRI-ESM2-0ssp534-overr1i1p1f1pr
29/home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0...MRI-ESM2-0MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.ncTrueMRI-ESM2-0ssp534-overr1i1p1f1tas
30/home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0...MRI-ESM2-0MRI-ESM2-0_ssp585_r1i1p1f1_pr.ncTrueMRI-ESM2-0ssp585r1i1p1f1pr
31/home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0...MRI-ESM2-0MRI-ESM2-0_ssp585_r1i1p1f1_tas.ncTrueMRI-ESM2-0ssp585r1i1p1f1tas
\n", "
" ], "text/plain": [ " path folder \\\n", "0 /home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2... ACCESS-CM2 \n", "1 /home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2... ACCESS-CM2 \n", "2 /home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2... ACCESS-CM2 \n", "3 /home/www/oggm/cmip6/GCM/ACCESS-CM2/ACCESS-CM2... ACCESS-CM2 \n", "4 /home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-... ACCESS-ESM1-5 \n", "5 /home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-... ACCESS-ESM1-5 \n", "6 /home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-... ACCESS-ESM1-5 \n", "7 /home/www/oggm/cmip6/GCM/ACCESS-ESM1-5/ACCESS-... ACCESS-ESM1-5 \n", "8 /home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC... CESM2-WACCM \n", "9 /home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC... CESM2-WACCM \n", "10 /home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC... CESM2-WACCM \n", "11 /home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC... CESM2-WACCM \n", "12 /home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC... CESM2-WACCM \n", "13 /home/www/oggm/cmip6/GCM/CESM2-WACCM/CESM2-WAC... CESM2-WACCM \n", "14 /home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp12... CanESM5 \n", "15 /home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp12... CanESM5 \n", "16 /home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp53... CanESM5 \n", "17 /home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp53... CanESM5 \n", "18 /home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp58... CanESM5 \n", "19 /home/www/oggm/cmip6/GCM/CanESM5/CanESM5_ssp58... CanESM5 \n", "20 /home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6... IPSL-CM6A-LR \n", "21 /home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6... IPSL-CM6A-LR \n", "22 /home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6... IPSL-CM6A-LR \n", "23 /home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6... IPSL-CM6A-LR \n", "24 /home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6... IPSL-CM6A-LR \n", "25 /home/www/oggm/cmip6/GCM/IPSL-CM6A-LR/IPSL-CM6... IPSL-CM6A-LR \n", "26 /home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0... MRI-ESM2-0 \n", "27 /home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0... MRI-ESM2-0 \n", "28 /home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0... MRI-ESM2-0 \n", "29 /home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0... MRI-ESM2-0 \n", "30 /home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0... MRI-ESM2-0 \n", "31 /home/www/oggm/cmip6/GCM/MRI-ESM2-0/MRI-ESM2-0... MRI-ESM2-0 \n", "\n", " fname valid gcm \\\n", "0 ACCESS-CM2_ssp126_r1i1p1f1_pr.nc True ACCESS-CM2 \n", "1 ACCESS-CM2_ssp126_r1i1p1f1_tas.nc True ACCESS-CM2 \n", "2 ACCESS-CM2_ssp585_r1i1p1f1_pr.nc True ACCESS-CM2 \n", "3 ACCESS-CM2_ssp585_r1i1p1f1_tas.nc True ACCESS-CM2 \n", "4 ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc True ACCESS-ESM1-5 \n", "5 ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc True ACCESS-ESM1-5 \n", "6 ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc True ACCESS-ESM1-5 \n", "7 ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc True ACCESS-ESM1-5 \n", "8 CESM2-WACCM_ssp126_r1i1p1f1_pr.nc True CESM2-WACCM \n", "9 CESM2-WACCM_ssp126_r1i1p1f1_tas.nc True CESM2-WACCM \n", "10 CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc True CESM2-WACCM \n", "11 CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc True CESM2-WACCM \n", "12 CESM2-WACCM_ssp585_r1i1p1f1_pr.nc True CESM2-WACCM \n", "13 CESM2-WACCM_ssp585_r1i1p1f1_tas.nc True CESM2-WACCM \n", "14 CanESM5_ssp126_r1i1p1f1_pr.nc True CanESM5 \n", "15 CanESM5_ssp126_r1i1p1f1_tas.nc True CanESM5 \n", "16 CanESM5_ssp534-over_r1i1p1f1_pr.nc True CanESM5 \n", "17 CanESM5_ssp534-over_r1i1p1f1_tas.nc True CanESM5 \n", "18 CanESM5_ssp585_r1i1p1f1_pr.nc True CanESM5 \n", "19 CanESM5_ssp585_r1i1p1f1_tas.nc True CanESM5 \n", "20 IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc True IPSL-CM6A-LR \n", "21 IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc True IPSL-CM6A-LR \n", "22 IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc True IPSL-CM6A-LR \n", "23 IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc True IPSL-CM6A-LR \n", "24 IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc True IPSL-CM6A-LR \n", "25 IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc True IPSL-CM6A-LR \n", "26 MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc True MRI-ESM2-0 \n", "27 MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc True MRI-ESM2-0 \n", "28 MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc True MRI-ESM2-0 \n", "29 MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc True MRI-ESM2-0 \n", "30 MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc True MRI-ESM2-0 \n", "31 MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc True MRI-ESM2-0 \n", "\n", " ssp rea var \n", "0 ssp126 r1i1p1f1 pr \n", "1 ssp126 r1i1p1f1 tas \n", "2 ssp585 r1i1p1f1 pr \n", "3 ssp585 r1i1p1f1 tas \n", "4 ssp126 r1i1p1f1 pr \n", "5 ssp126 r1i1p1f1 tas \n", "6 ssp585 r1i1p1f1 pr \n", "7 ssp585 r1i1p1f1 tas \n", "8 ssp126 r1i1p1f1 pr \n", "9 ssp126 r1i1p1f1 tas \n", "10 ssp534-over r1i1p1f1 pr \n", "11 ssp534-over r1i1p1f1 tas \n", "12 ssp585 r1i1p1f1 pr \n", "13 ssp585 r1i1p1f1 tas \n", "14 ssp126 r1i1p1f1 pr \n", "15 ssp126 r1i1p1f1 tas \n", "16 ssp534-over r1i1p1f1 pr \n", "17 ssp534-over r1i1p1f1 tas \n", "18 ssp585 r1i1p1f1 pr \n", "19 ssp585 r1i1p1f1 tas \n", "20 ssp126 r1i1p1f1 pr \n", "21 ssp126 r1i1p1f1 tas \n", "22 ssp534-over r1i1p1f1 pr \n", "23 ssp534-over r1i1p1f1 tas \n", "24 ssp585 r1i1p1f1 pr \n", "25 ssp585 r1i1p1f1 tas \n", "26 ssp126 r1i1p1f1 pr \n", "27 ssp126 r1i1p1f1 tas \n", "28 ssp534-over r1i1p1f1 pr \n", "29 ssp534-over r1i1p1f1 tas \n", "30 ssp585 r1i1p1f1 pr \n", "31 ssp585 r1i1p1f1 tas " ] }, "execution_count": 109, "metadata": {}, "output_type": "execute_result" } ], "source": [ "gcms_cmip6" ] }, { "cell_type": "code", "execution_count": 108, "id": "a418f414-aaa5-45e5-b264-58a13041888f", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 76, "id": "c309b075-eaaa-48d3-9991-aeecbf3c4f49", "metadata": {}, "outputs": [], "source": [ "assert len(gcms_cmip6) == 32" ] }, { "cell_type": "code", "execution_count": 80, "id": "816c2995-b60d-4fb4-b871-6baefe21ae94", "metadata": {}, "outputs": [], "source": [ "gcms_cmip6_2300_bef_merge = pd.read_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_vJul2023.csv', index_col='Unnamed: 0')\n" ] }, { "cell_type": "code", "execution_count": 85, "id": "6ad04b8d-ea46-4107-b994-1a448570bd8e", "metadata": {}, "outputs": [], "source": [ "for fname in gcms_cmip6_2300_bef_merge.fname.values:\n", " if fname in gcms_cmip6.fname.values:\n", " pass\n", " else:\n", " print(f'{fname} is missing')" ] }, { "cell_type": "code", "execution_count": null, "id": "396cbfee-63fc-4d6a-a48e-0bae0fbaaebf", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "e00874fa-cf16-46bf-9a55-5d747229c60d", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 18, "id": "39311c3e-1493-4811-914a-053b47849f5b", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "2100-12-15 12:00:00\n", "2299-12-15 12:00:00\n", "2100-12-15 12:00:00\n", "2299-12-15 12:00:00\n" ] } ], "source": [ "for v in ['pr', 'tas']:\n", " do = xr.open_dataset(f'/home/www/oggm/cmip6/GCM/CESM2-WACCM//CESM2-WACCM_ssp126_r1i1p1f1_{v}.nc')[v]\n", " dn = xr.open_dataset(f'/home/www/oggm/cmip6/redundant/redundant_nico/CESM2-WACCM_ssp126_r1i1p1f1_{v}.nc')\n", " print(do.isel(time=-1).time.values)\n", " print(dn.isel(time=-1).time.values)\n", " dn = dn.sel(time=slice('1850','2100'))[v]\n", " np.testing.assert_allclose(do, dn)" ] }, { "cell_type": "code", "execution_count": 40, "id": "194a69ea-4692-4411-8d83-896a15380397", "metadata": {}, "outputs": [], "source": [ "end_yr=2100" ] }, { "cell_type": "code", "execution_count": 49, "id": "80c00839-ec82-47c5-a10a-729461e73352", "metadata": {}, "outputs": [], "source": [ "df1 = df1.loc[df1.rcp!='historicalGHG']\n", "df1 = df1.loc[df1.rcp!='historicalNat']" ] }, { "cell_type": "code", "execution_count": 50, "id": "5d3e3672-b5b7-409d-ab41-f0e0431eeee9", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "_CMIP5_tas_mon_CanESM2_rcp26_r1i1p1_g025.nc\n", "_CMIP5_tas_mon_CanESM2_rcp85_r1i1p1_g025.nc\n", "_CMIP5_tas_mon_CanESM2_rcp45_r1i1p1_g025.nc\n" ] } ], "source": [ " gcms_cmip5 = pd.read_csv('/home/www/oggm/cmip5-ng/all_gcm_list.csv', index_col=0)\n", " if end_yr==2300:\n", " raise InvalidParamsError('No CMIP5 GCMs until 2300!')\n", " for gcm in ['CanESM2']: #gcms_cmip5.gcm.unique(): ##TEST:\n", " df1 = gcms_cmip5.loc[gcms_cmip5.gcm == gcm]\n", " df1 = df1.loc[df1.rcp!='historicalGHG']\n", " df1 = df1.loc[df1.rcp!='historicalNat']\n", " for rcp in df1.rcp.unique():\n", " df2 = df1.loc[df1.rcp == rcp]\n", " assert len(df2) == 2\n", " ft = df2.loc[df2['var'] == 'tas'].iloc[0]\n", " rid = ft.fname.replace('_r1i1p1f1_tas.nc', '')\n", " rid='_CMIP5_' + rid\n", " print(rid)" ] }, { "cell_type": "code", "execution_count": 42, "id": "52d9ef2a-73b0-484f-bcc4-c90062cab1d5", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'_CMIP5_tas_mon_GFDL-ESM2G_rcp45_r1i1p1_g025.nc'" ] }, "execution_count": 42, "metadata": {}, "output_type": "execute_result" } ], "source": [ "rid" ] }, { "cell_type": "code", "execution_count": 23, "id": "2d96c439-0a96-42e4-9724-2d2ed1c07c7d", "metadata": {}, "outputs": [], "source": [ " gcms_cmip5 = pd.read_csv('/home/www/oggm/cmip5-ng/all_gcm_list.csv', index_col=0)\n" ] }, { "cell_type": "code", "execution_count": null, "id": "2cd156cf-b51e-4be7-8a59-306f913685ad", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 37, "id": "00020ce4-0a39-47fd-bf3a-a2f7925e7959", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array(['rcp45', 'historicalGHG', 'rcp60', 'historicalNat', 'rcp26',\n", " 'rcp85'], dtype=object)" ] }, "execution_count": 37, "metadata": {}, "output_type": "execute_result" } ], "source": [ "gcms_cmip5.rcp.unique()" ] }, { "cell_type": "code", "execution_count": 38, "id": "557e0f7b-ee0d-4f8a-9d90-49ad56ff665a", "metadata": {}, "outputs": [ { "ename": "ValueError", "evalue": "The truth value of a Series is ambiguous. Use a.empty, a.bool(), a.item(), a.any() or a.all().", "output_type": "error", "traceback": [ "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m", "\u001b[0;31mValueError\u001b[0m Traceback (most recent call last)", "Input \u001b[0;32mIn [38]\u001b[0m, in \u001b[0;36m\u001b[0;34m()\u001b[0m\n\u001b[0;32m----> 1\u001b[0m gcms_cmip5\u001b[38;5;241m.\u001b[39mloc[\u001b[43mgcms_cmip5\u001b[49m\u001b[38;5;241;43m.\u001b[39;49m\u001b[43mrcp\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;129;43;01min\u001b[39;49;00m\u001b[43m \u001b[49m\u001b[43m[\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[38;5;124;43mrcp26\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[43m,\u001b[49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[38;5;124;43mrcp45\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[38;5;124;43mrcp60\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[43m,\u001b[49m\u001b[43m \u001b[49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[38;5;124;43mrcp85\u001b[39;49m\u001b[38;5;124;43m'\u001b[39;49m\u001b[43m]\u001b[49m]\n", "File \u001b[0;32m~/.local/lib/python3.10/site-packages/pandas/core/generic.py:1527\u001b[0m, in \u001b[0;36mNDFrame.__nonzero__\u001b[0;34m(self)\u001b[0m\n\u001b[1;32m 1525\u001b[0m \u001b[38;5;129m@final\u001b[39m\n\u001b[1;32m 1526\u001b[0m \u001b[38;5;28;01mdef\u001b[39;00m \u001b[38;5;21m__nonzero__\u001b[39m(\u001b[38;5;28mself\u001b[39m):\n\u001b[0;32m-> 1527\u001b[0m \u001b[38;5;28;01mraise\u001b[39;00m \u001b[38;5;167;01mValueError\u001b[39;00m(\n\u001b[1;32m 1528\u001b[0m \u001b[38;5;124mf\u001b[39m\u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mThe truth value of a \u001b[39m\u001b[38;5;132;01m{\u001b[39;00m\u001b[38;5;28mtype\u001b[39m(\u001b[38;5;28mself\u001b[39m)\u001b[38;5;241m.\u001b[39m\u001b[38;5;18m__name__\u001b[39m\u001b[38;5;132;01m}\u001b[39;00m\u001b[38;5;124m is ambiguous. \u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 1529\u001b[0m \u001b[38;5;124m\"\u001b[39m\u001b[38;5;124mUse a.empty, a.bool(), a.item(), a.any() or a.all().\u001b[39m\u001b[38;5;124m\"\u001b[39m\n\u001b[1;32m 1530\u001b[0m )\n", "\u001b[0;31mValueError\u001b[0m: The truth value of a Series is ambiguous. Use a.empty, a.bool(), a.item(), a.any() or a.all()." ] } ], "source": [ "gcms_cmip5.loc[gcms_cmip5.rcp in ['rcp26','rcp45','rcp60', 'rcp85']]" ] }, { "cell_type": "code", "execution_count": 35, "id": "07c035a3-a546-445f-80f8-b581aefdabe1", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:  (time: 1632, lat: 72, lon: 144)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1870-01-15T12:00:00 ... 2005-12-15T12:00:00\n",
       "  * lat      (lat) float64 -88.75 -86.25 -83.75 -81.25 ... 83.75 86.25 88.75\n",
       "  * lon      (lon) float64 1.25 3.75 6.25 8.75 11.25 ... 351.2 353.8 356.2 358.8\n",
       "Data variables:\n",
       "    pr       (time, lat, lon) float32 ...\n",
       "Attributes:\n",
       "    source_model:        CCSM4\n",
       "    source_ensemble:     r1i1p1\n",
       "    source_experiment:   historicalNat\n",
       "    source_institution:  NCAR (National Center for Atmospheric Research) Boul...\n",
       "    source_contact:      cesm_data@ucar.edu\n",
       "    source_files:        pr_Amon_CCSM4_historicalNat_r1i1p1_185001-200512.nc\n",
       "    source_md5:          1d48dde930ef3772a4128a3b03ff8e78\n",
       "    freq:                monthly\n",
       "    interpolation_grid:  2.5 x 2.5 degrees bilinear interpolation\n",
       "    contact:             cmip5-archive@env.ethz.ch\n",
       "    modifications:       nothing
" ], "text/plain": [ "\n", "Dimensions: (time: 1632, lat: 72, lon: 144)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1870-01-15T12:00:00 ... 2005-12-15T12:00:00\n", " * lat (lat) float64 -88.75 -86.25 -83.75 -81.25 ... 83.75 86.25 88.75\n", " * lon (lon) float64 1.25 3.75 6.25 8.75 11.25 ... 351.2 353.8 356.2 358.8\n", "Data variables:\n", " pr (time, lat, lon) float32 ...\n", "Attributes:\n", " source_model: CCSM4\n", " source_ensemble: r1i1p1\n", " source_experiment: historicalNat\n", " source_institution: NCAR (National Center for Atmospheric Research) Boul...\n", " source_contact: cesm_data@ucar.edu\n", " source_files: pr_Amon_CCSM4_historicalNat_r1i1p1_185001-200512.nc\n", " source_md5: 1d48dde930ef3772a4128a3b03ff8e78\n", " freq: monthly\n", " interpolation_grid: 2.5 x 2.5 degrees bilinear interpolation\n", " contact: cmip5-archive@env.ethz.ch\n", " modifications: nothing" ] }, "execution_count": 35, "metadata": {}, "output_type": "execute_result" } ], "source": [ "xr.open_dataset(gcms_cmip5['path'][3])" ] }, { "cell_type": "code", "execution_count": 27, "id": "9758548b-4192-4da3-bf48-3922d63b418f", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "'/home/www/oggm/cmip5-ng/pr/pr_mon_NorESM1-M_rcp45_r1i1p1_g025.nc'" ] }, "execution_count": 27, "metadata": {}, "output_type": "execute_result" } ], "source": [] }, { "cell_type": "markdown", "id": "6193c4e7-783e-4e2f-91d3-9223a477ac3f", "metadata": {}, "source": [ "## CMIP5 until 2300" ] }, { "cell_type": "code", "execution_count": 16, "id": "bec9b27b-bf16-4ba5-a945-10024018c377", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "['tas', 'mon', 'CESM1-CAM5', 'rcp60', 'r1i1p1.nc']" ] }, "execution_count": 16, "metadata": {}, "output_type": "execute_result" } ], "source": [ "fname.split('_')" ] }, { "cell_type": "code", "execution_count": 64, "id": "d2bfd5f7-da32-4726-bcad-4fedd6c85513", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "['/home/www/rchizzola/cmip5_done/ccsm4/tas_mon_CCSM4_rcp60_r1i1p1.nc',\n", " '/home/www/rchizzola/cmip5_done/ccsm4/tas_mon_CCSM4_rcp85_r1i1p1.nc']" ] }, "execution_count": 64, "metadata": {}, "output_type": "execute_result" } ], "source": [ "files_l_tas2" ] }, { "cell_type": "code", "execution_count": 63, "id": "40578208-4f04-4790-a789-0aa8f706523a", "metadata": {}, "outputs": [], "source": [ "files_l_tas = glob.glob(\"/home/www/rchizzola/cmip5_done/tas/*\")\n", "files_l_pr = glob.glob(\"/home/www/rchizzola/cmip5_done/pr/*\")\n", "files_l_tas2 = glob.glob(\"/home/www/rchizzola/cmip5_done/ccsm4/tas*\")\n", "files_l_pr2 = glob.glob(\"/home/www/rchizzola/cmip5_done/ccsm4/pr*\")\n", "files_l_tas = files_l_tas2 + files_l_tas\n", "files_l_pr = files_l_pr2 + files_l_pr\n", "\n", "files_l_rebekka = files_l_tas + files_l_pr \n", "\n", "load = True\n", "if load:\n", " gcms_cmip5_2300 = pd.DataFrame()\n", "\n", " for fpath in files_l_rebekka:\n", " if '.nc' in fpath:\n", " fname = fpath.split('/')[-1]\n", " folder_p = fpath[:-len(fname)]\n", " var, _, gcm, rcp, realisation = fname.split('_')\n", " gcms_cmip5_2300.loc[fname, 'rcp'] = rcp\n", "\n", " gcms_cmip5_2300.loc[fname, 'gcm'] = gcm\n", " gcms_cmip5_2300.loc[fname, 'rea'] = realisation\n", " gcms_cmip5_2300.loc[fname, 'var'] = var\n", " gcms_cmip5_2300.loc[fname, 'path'] = fpath\n", " gcms_cmip5_2300.loc[fname, 'fname'] = fname\n", "\n", " gcms_cmip5_2300.to_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_cmip5_vAug2023.csv')\n", " gcms_cmip6_2300.to_html('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_cmip5_vAug2023.html')\n", "\n", "else:\n", " gcms_cmip5_2300 = pd.read_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_cmip5_vAug2023.csv', index_col='Unnamed: 0')\n", " #gcms_cmip6_2300 = pd.read_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_update.csv', index_col='Unnamed: 0')" ] }, { "cell_type": "code", "execution_count": 65, "id": "76dbc9e4-f46c-4b6b-94e6-f25d1a418ea9", "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", " \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", "
rcpgcmreavarpathfname
pr_mon_CCSM4_rcp60_r1i1p1.ncrcp60CCSM4r1i1p1.ncpr/home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CC...pr_mon_CCSM4_rcp60_r1i1p1.nc
pr_mon_CCSM4_rcp85_r1i1p1.ncrcp85CCSM4r1i1p1.ncpr/home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CC...pr_mon_CCSM4_rcp85_r1i1p1.nc
pr_mon_CESM1-CAM5_rcp26_r1i1p1.ncrcp26CESM1-CAM5r1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1...pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc
pr_mon_CESM1-CAM5_rcp45_r1i1p1.ncrcp45CESM1-CAM5r1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1...pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc
pr_mon_CESM1-CAM5_rcp60_r1i1p1.ncrcp60CESM1-CAM5r1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1...pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc
pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.ncrcp45CSIRO-Mk3-6-0r1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO...pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc
pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.ncrcp85CSIRO-Mk3-6-0r1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO...pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc
pr_mon_CanESM2_rcp26_r1i1p1.ncrcp26CanESM2r1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_CanES...pr_mon_CanESM2_rcp26_r1i1p1.nc
pr_mon_MPI-ESM-LR_rcp26_r1i1p1.ncrcp26MPI-ESM-LRr1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-E...pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc
pr_mon_MPI-ESM-LR_rcp45_r1i1p1.ncrcp45MPI-ESM-LRr1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-E...pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc
pr_mon_MPI-ESM-LR_rcp85_r1i1p1.ncrcp85MPI-ESM-LRr1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-E...pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc
pr_mon_NorESM1-M_rcp45_r1i1p1.ncrcp45NorESM1-Mr1i1p1.ncpr/home/www/rchizzola/cmip5_done/pr/pr_mon_NorES...pr_mon_NorESM1-M_rcp45_r1i1p1.nc
tas_mon_CCSM4_rcp60_r1i1p1.ncrcp60CCSM4r1i1p1.nctas/home/www/rchizzola/cmip5_done/ccsm4/tas_mon_C...tas_mon_CCSM4_rcp60_r1i1p1.nc
tas_mon_CCSM4_rcp85_r1i1p1.ncrcp85CCSM4r1i1p1.nctas/home/www/rchizzola/cmip5_done/ccsm4/tas_mon_C...tas_mon_CCSM4_rcp85_r1i1p1.nc
tas_mon_CESM1-CAM5_rcp26_r1i1p1.ncrcp26CESM1-CAM5r1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_CES...tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc
tas_mon_CESM1-CAM5_rcp45_r1i1p1.ncrcp45CESM1-CAM5r1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_CES...tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc
tas_mon_CESM1-CAM5_rcp60_r1i1p1.ncrcp60CESM1-CAM5r1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_CES...tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc
tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.ncrcp45CSIRO-Mk3-6-0r1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_CSI...tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc
tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.ncrcp85CSIRO-Mk3-6-0r1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_CSI...tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc
tas_mon_CanESM2_rcp26_r1i1p1.ncrcp26CanESM2r1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_Can...tas_mon_CanESM2_rcp26_r1i1p1.nc
tas_mon_MPI-ESM-LR_rcp26_r1i1p1.ncrcp26MPI-ESM-LRr1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_MPI...tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc
tas_mon_MPI-ESM-LR_rcp45_r1i1p1.ncrcp45MPI-ESM-LRr1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_MPI...tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc
tas_mon_MPI-ESM-LR_rcp85_r1i1p1.ncrcp85MPI-ESM-LRr1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_MPI...tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc
tas_mon_NorESM1-M_rcp45_r1i1p1.ncrcp45NorESM1-Mr1i1p1.nctas/home/www/rchizzola/cmip5_done/tas/tas_mon_Nor...tas_mon_NorESM1-M_rcp45_r1i1p1.nc
\n", "
" ], "text/plain": [ " rcp gcm rea var \\\n", "pr_mon_CCSM4_rcp60_r1i1p1.nc rcp60 CCSM4 r1i1p1.nc pr \n", "pr_mon_CCSM4_rcp85_r1i1p1.nc rcp85 CCSM4 r1i1p1.nc pr \n", "pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc rcp26 CESM1-CAM5 r1i1p1.nc pr \n", "pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc rcp45 CESM1-CAM5 r1i1p1.nc pr \n", "pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc rcp60 CESM1-CAM5 r1i1p1.nc pr \n", "pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc rcp45 CSIRO-Mk3-6-0 r1i1p1.nc pr \n", "pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc rcp85 CSIRO-Mk3-6-0 r1i1p1.nc pr \n", "pr_mon_CanESM2_rcp26_r1i1p1.nc rcp26 CanESM2 r1i1p1.nc pr \n", "pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc rcp26 MPI-ESM-LR r1i1p1.nc pr \n", "pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc rcp45 MPI-ESM-LR r1i1p1.nc pr \n", "pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc rcp85 MPI-ESM-LR r1i1p1.nc pr \n", "pr_mon_NorESM1-M_rcp45_r1i1p1.nc rcp45 NorESM1-M r1i1p1.nc pr \n", "tas_mon_CCSM4_rcp60_r1i1p1.nc rcp60 CCSM4 r1i1p1.nc tas \n", "tas_mon_CCSM4_rcp85_r1i1p1.nc rcp85 CCSM4 r1i1p1.nc tas \n", "tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc rcp26 CESM1-CAM5 r1i1p1.nc tas \n", "tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc rcp45 CESM1-CAM5 r1i1p1.nc tas \n", "tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc rcp60 CESM1-CAM5 r1i1p1.nc tas \n", "tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc rcp45 CSIRO-Mk3-6-0 r1i1p1.nc tas \n", "tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc rcp85 CSIRO-Mk3-6-0 r1i1p1.nc tas \n", "tas_mon_CanESM2_rcp26_r1i1p1.nc rcp26 CanESM2 r1i1p1.nc tas \n", "tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc rcp26 MPI-ESM-LR r1i1p1.nc tas \n", "tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc rcp45 MPI-ESM-LR r1i1p1.nc tas \n", "tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc rcp85 MPI-ESM-LR r1i1p1.nc tas \n", "tas_mon_NorESM1-M_rcp45_r1i1p1.nc rcp45 NorESM1-M r1i1p1.nc tas \n", "\n", " path \\\n", "pr_mon_CCSM4_rcp60_r1i1p1.nc /home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CC... \n", "pr_mon_CCSM4_rcp85_r1i1p1.nc /home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CC... \n", "pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1... \n", "pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1... \n", "pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1... \n", "pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO... \n", "pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO... \n", "pr_mon_CanESM2_rcp26_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_CanES... \n", "pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-E... \n", "pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-E... \n", "pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-E... \n", "pr_mon_NorESM1-M_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/pr/pr_mon_NorES... \n", "tas_mon_CCSM4_rcp60_r1i1p1.nc /home/www/rchizzola/cmip5_done/ccsm4/tas_mon_C... \n", "tas_mon_CCSM4_rcp85_r1i1p1.nc /home/www/rchizzola/cmip5_done/ccsm4/tas_mon_C... \n", "tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_CES... \n", "tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_CES... \n", "tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_CES... \n", "tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_CSI... \n", "tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_CSI... \n", "tas_mon_CanESM2_rcp26_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_Can... \n", "tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_MPI... \n", "tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_MPI... \n", "tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_MPI... \n", "tas_mon_NorESM1-M_rcp45_r1i1p1.nc /home/www/rchizzola/cmip5_done/tas/tas_mon_Nor... \n", "\n", " fname \n", "pr_mon_CCSM4_rcp60_r1i1p1.nc pr_mon_CCSM4_rcp60_r1i1p1.nc \n", "pr_mon_CCSM4_rcp85_r1i1p1.nc pr_mon_CCSM4_rcp85_r1i1p1.nc \n", "pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc \n", "pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc \n", "pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc \n", "pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc \n", "pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc \n", "pr_mon_CanESM2_rcp26_r1i1p1.nc pr_mon_CanESM2_rcp26_r1i1p1.nc \n", "pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc \n", "pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc \n", "pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc \n", "pr_mon_NorESM1-M_rcp45_r1i1p1.nc pr_mon_NorESM1-M_rcp45_r1i1p1.nc \n", "tas_mon_CCSM4_rcp60_r1i1p1.nc tas_mon_CCSM4_rcp60_r1i1p1.nc \n", "tas_mon_CCSM4_rcp85_r1i1p1.nc tas_mon_CCSM4_rcp85_r1i1p1.nc \n", "tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc \n", "tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc \n", "tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc \n", "tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc \n", "tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc \n", "tas_mon_CanESM2_rcp26_r1i1p1.nc tas_mon_CanESM2_rcp26_r1i1p1.nc \n", "tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc \n", "tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc \n", "tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc \n", "tas_mon_NorESM1-M_rcp45_r1i1p1.nc tas_mon_NorESM1-M_rcp45_r1i1p1.nc " ] }, "execution_count": 65, "metadata": {}, "output_type": "execute_result" } ], "source": [ "\n", "gcms_cmip5_2300.sort_values(by='fname')" ] }, { "cell_type": "code", "execution_count": 66, "id": "79c5e289-21cb-41c1-a8a1-ee82a0e4e2e8", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "12.0" ] }, "execution_count": 66, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# amount of GCMS * SSPs\n", "len(gcms_cmip5_2300)/2" ] }, { "cell_type": "code", "execution_count": 67, "id": "afbb808f-cf48-4bbb-ab16-800d616dec58", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array(['CCSM4', 'CESM1-CAM5', 'MPI-ESM-LR', 'NorESM1-M', 'CanESM2',\n", " 'CSIRO-Mk3-6-0'], dtype=object)" ] }, "execution_count": 67, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# amount of GCMs\n", "gcms_cmip5_2300.gcm.unique()" ] }, { "cell_type": "code", "execution_count": 68, "id": "8280771d-9489-422c-b97c-5af3a1ba4ee4", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array(['rcp60', 'rcp85', 'rcp26', 'rcp45'], dtype=object)" ] }, "execution_count": 68, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# amount of SSPs\n", "gcms_cmip5_2300.rcp.unique()\n" ] }, { "cell_type": "code", "execution_count": 69, "id": "e63dbb3f-673a-4c14-a55d-cf1721c6a8cb", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "87.77165226638317" ] }, "execution_count": 69, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#xr.open_dataset('/home/www/oggm/cmip6/GCM/CAMS-CSM1-0/CAMS-CSM1-0_ssp119_r2i1p1f1_tas.nc').time\n", "dold_pr = xr.open_dataset('/home/www/oggm/cmip5-ng/pr/pr_mon_CCSM4_rcp26_r1i1p1_g025.nc')\n", "#dold_tas = xr.open_dataset('/home/www/oggm/cmip6/GCM/GFDL-ESM4/GFDL-ESM4_ssp119_r1i1p1f1_tas.nc')\n", "pr_max_other_cmip5 = dold_pr.pr.max().values*60*60*24 # kg m-2 day-1\n", "pr_max_other_cmip5" ] }, { "cell_type": "code", "execution_count": 70, "id": "3a247a3d-37ab-4c1f-befb-349756c1127c", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:  (time: 2772, lat: 72, lon: 144)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1870-01-15T12:00:00 ... 2100-12-15T12:00:00\n",
       "  * lat      (lat) float64 -88.75 -86.25 -83.75 -81.25 ... 83.75 86.25 88.75\n",
       "  * lon      (lon) float64 1.25 3.75 6.25 8.75 11.25 ... 351.2 353.8 356.2 358.8\n",
       "Data variables:\n",
       "    pr       (time, lat, lon) float32 1.514e-06 1.579e-06 ... 7.333e-06\n",
       "Attributes:\n",
       "    source_model:        CCSM4\n",
       "    source_ensemble:     r1i1p1\n",
       "    source_experiment:   historical,rcp26\n",
       "    source_institution:  NCAR (National Center for Atmospheric Research) Boul...\n",
       "    source_contact:      cesm_data@ucar.edu\n",
       "    source_files:        pr_Amon_CCSM4_historical_r1i1p1_185001-200512.nc,pr_...\n",
       "    source_md5:          cfc87ae596265b9eba109c0d8b29ace8,436d6b5447194e7fbb6...\n",
       "    freq:                monthly\n",
       "    interpolation_grid:  2.5 x 2.5 degrees bilinear interpolation\n",
       "    contact:             cmip5-archive@env.ethz.ch\n",
       "    modifications:       nothing
" ], "text/plain": [ "\n", "Dimensions: (time: 2772, lat: 72, lon: 144)\n", "Coordinates:\n", " * time (time) datetime64[ns] 1870-01-15T12:00:00 ... 2100-12-15T12:00:00\n", " * lat (lat) float64 -88.75 -86.25 -83.75 -81.25 ... 83.75 86.25 88.75\n", " * lon (lon) float64 1.25 3.75 6.25 8.75 11.25 ... 351.2 353.8 356.2 358.8\n", "Data variables:\n", " pr (time, lat, lon) float32 1.514e-06 1.579e-06 ... 7.333e-06\n", "Attributes:\n", " source_model: CCSM4\n", " source_ensemble: r1i1p1\n", " source_experiment: historical,rcp26\n", " source_institution: NCAR (National Center for Atmospheric Research) Boul...\n", " source_contact: cesm_data@ucar.edu\n", " source_files: pr_Amon_CCSM4_historical_r1i1p1_185001-200512.nc,pr_...\n", " source_md5: cfc87ae596265b9eba109c0d8b29ace8,436d6b5447194e7fbb6...\n", " freq: monthly\n", " interpolation_grid: 2.5 x 2.5 degrees bilinear interpolation\n", " contact: cmip5-archive@env.ethz.ch\n", " modifications: nothing" ] }, "execution_count": 70, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dold_pr" ] }, { "cell_type": "code", "execution_count": 71, "id": "c0bc8583-8feb-427b-bfea-650771aa902a", "metadata": {}, "outputs": [], "source": [ "weight = np.cos(np.deg2rad(dold_pr.lat.astype(np.float64))).clip(0)\n", "weight = dold_pr.pr.isel(time=0) * 0. + weight\n", "dold_pr['weight'] = (('lat', 'lon'), weight.values / weight.sum().values)\n", "global_avg_pr_daily_old = ((dold_pr.weight * dold_pr.pr).sum(dim='lon').sum(dim='lat')*60*60*24).max().values" ] }, { "cell_type": "code", "execution_count": 72, "id": "3154205a-7407-47c7-a35b-b023cffd5dc0", "metadata": {}, "outputs": [], "source": [ "dold_pr_lon_mean_daily = dold_pr.pr.mean(dim='lon')" ] }, { "cell_type": "code", "execution_count": 73, "id": "f9304950-e0b7-4a4a-b783-8ec1f8ee59b1", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "5412 192 288 /home/www/rchizzola/cmip5_done/ccsm4/tas_mon_CCSM4_rcp60_r1i1p1.nc\n", "5412 192 288 /home/www/rchizzola/cmip5_done/ccsm4/tas_mon_CCSM4_rcp85_r1i1p1.nc\n", "5412 192 288 /home/www/rchizzola/cmip5_done/tas/tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc\n", "5412 96 192 /home/www/rchizzola/cmip5_done/tas/tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc\n", "5412 96 144 /home/www/rchizzola/cmip5_done/tas/tas_mon_NorESM1-M_rcp45_r1i1p1.nc\n", "5412 192 288 /home/www/rchizzola/cmip5_done/tas/tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc\n", "5412 64 128 /home/www/rchizzola/cmip5_done/tas/tas_mon_CanESM2_rcp26_r1i1p1.nc\n", "5412 96 192 /home/www/rchizzola/cmip5_done/tas/tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc\n", "5412 96 192 /home/www/rchizzola/cmip5_done/tas/tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc\n", "5412 192 288 /home/www/rchizzola/cmip5_done/tas/tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc\n", "5412 96 192 /home/www/rchizzola/cmip5_done/tas/tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc\n", "5412 96 192 /home/www/rchizzola/cmip5_done/tas/tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc\n" ] } ], "source": [ "for p in gcms_cmip5_2300.loc[gcms_cmip5_2300['var']=='tas'].path:\n", " dtest = xr.open_dataset(p, use_cftime=True) \n", " time, lat, lon = np.shape(dtest.tas)\n", " # either from 1850 until 2300 or until 2999\n", " try:\n", " assert time==5412\n", " except:\n", " assert time==5400\n", " print( time, lat,lon, p)" ] }, { "cell_type": "code", "execution_count": 74, "id": "5f13b89b-9d8d-4f06-8a29-7866f32c4a75", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "5412 192 288 /home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CCSM4_rcp60_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 192 288 /home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CCSM4_rcp85_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 96 192 /home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 64 128 /home/www/rchizzola/cmip5_done/pr/pr_mon_CanESM2_rcp26_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 96 192 /home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 96 192 /home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 96 192 /home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 192 288 /home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 96 144 /home/www/rchizzola/cmip5_done/pr/pr_mon_NorESM1-M_rcp45_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 96 192 /home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 192 288 /home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc 1850-01-16 12:00:00\n", "5412 192 288 /home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc 1850-01-16 12:00:00\n" ] } ], "source": [ "for p in gcms_cmip5_2300.loc[gcms_cmip5_2300['var']=='pr'].path:\n", " dtest = xr.open_dataset(p, use_cftime=True) \n", " time, lat, lon = np.shape(dtest.pr)\n", " # either from 1850 until 2300 or until 2999\n", " try:\n", " assert time==5412\n", " except:\n", " assert time==5400\n", " print(time, lat,lon, p, dtest.time[0].values )" ] }, { "cell_type": "code", "execution_count": 75, "id": "47d4d118-d649-4cf2-b67b-a880c46549af", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/home/www/rchizzola/cmip5_done/ccsm4/tas_mon_CCSM4_rcp60_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -72.0°C, tas_max = 48.1°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/ccsm4/tas_mon_CCSM4_rcp85_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -71.3°C, tas_max = 53.6°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_CESM1-CAM5_rcp60_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -80.8°C, tas_max = 50.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_MPI-ESM-LR_rcp26_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, tas_min = -76.8°C, tas_max = 47.8°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_NorESM1-M_rcp45_r1i1p1.nc: lat_resol=1.9°, lon_resol=2.5°, tas_min = -71.0°C, tas_max = 44.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_CESM1-CAM5_rcp45_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -80.8°C, tas_max = 49.0°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_CanESM2_rcp26_r1i1p1.nc: lat_resol=2.8°, lon_resol=2.8°, tas_min = -79.5°C, tas_max = 46.8°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_MPI-ESM-LR_rcp45_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, tas_min = -76.8°C, tas_max = 49.6°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_MPI-ESM-LR_rcp85_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, tas_min = -76.8°C, tas_max = 59.0°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_CESM1-CAM5_rcp26_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -80.8°C, tas_max = 47.0°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, tas_min = -82.6°C, tas_max = 63.2°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/tas/tas_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, tas_min = -82.6°C, tas_max = 49.6°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CCSM4_rcp60_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = 0.0 kg m-2 day-1, pr_max = 141.3 kg m-2 day-1, pr_global_avg_max = 3.4 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/ccsm4/pr_mon_CCSM4_rcp85_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = 0.0 kg m-2 day-1, pr_max = 157.5 kg m-2 day-1, pr_global_avg_max = 3.6 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO-Mk3-6-0_rcp45_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 64.2 kg m-2 day-1, pr_global_avg_max = 3.3 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_CanESM2_rcp26_r1i1p1.nc: lat_resol=2.8°, lon_resol=2.8°, pr_min = 0.0 kg m-2 day-1, pr_max = 91.0 kg m-2 day-1, pr_global_avg_max = 3.0 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_CSIRO-Mk3-6-0_rcp85_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 100.5 kg m-2 day-1, pr_global_avg_max = 3.9 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-ESM-LR_rcp45_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 70.7 kg m-2 day-1, pr_global_avg_max = 3.3 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-ESM-LR_rcp85_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 79.5 kg m-2 day-1, pr_global_avg_max = 3.6 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp26_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = -0.0 kg m-2 day-1, pr_max = 148.0 kg m-2 day-1, pr_global_avg_max = 3.3 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_NorESM1-M_rcp45_r1i1p1.nc: lat_resol=1.9°, lon_resol=2.5°, pr_min = 0.0 kg m-2 day-1, pr_max = 88.5 kg m-2 day-1, pr_global_avg_max = 3.1 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_MPI-ESM-LR_rcp26_r1i1p1.nc: lat_resol=1.8°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 69.4 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp45_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = -0.0 kg m-2 day-1, pr_max = 166.7 kg m-2 day-1, pr_global_avg_max = 3.4 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/cmip5_done/pr/pr_mon_CESM1-CAM5_rcp60_r1i1p1.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = -0.0 kg m-2 day-1, pr_max = 153.7 kg m-2 day-1, pr_global_avg_max = 3.5 kg m-2 day-1\n", "1850-01-16 12:00:00\n" ] } ], "source": [ "pd_max_prcp_l = []\n", "for f in files_l_rebekka:\n", " if '.nc' in f:\n", " dtest = xr.open_dataset(f,use_cftime=True) \n", " # Check resolution:\n", " lat_resol = np.unique((dtest.lat.values[1:] - dtest.lat.values[:-1]).round(1))[0]\n", " lon_resol = np.unique((dtest.lon.values[1:] - dtest.lon.values[:-1]).round(1))[0]\n", " assert lon_resol<3\n", " assert lat_resol<3\n", " \n", " # Check units\n", " if 'pr' in f:\n", " pr_max = (dtest.pr.max()*60*60*24).values # kg m-2 day-1\n", " pr_min = (dtest.pr.min()*60*60*24).values\n", " try:\n", " assert pr_max=dold_pr.pr.max()*4).pr\n", " dmax = dmax.dropna(dim='lon', how='all').dropna(dim='lat', how='all').dropna(dim='time', how='all')\n", " _pd_max_prcp = (dmax*60*60*24).to_dataframe('daily_prcp_max').dropna()\n", " _pd_max_prcp['file'] = f\n", " pd_max_prcp_l.append(_pd_max_prcp)\n", " # global avg. precipitation \n", " weight = np.cos(np.deg2rad(dtest.lat.astype(np.float64))).clip(0)\n", " weight = dtest.pr.isel(time=0) * 0. + weight\n", " dtest['weight'] = (('lat', 'lon'), weight.values / weight.sum().values)\n", " global_avg_pr_daily_test = (dtest.weight * dtest.pr).sum(dim='lon').sum(dim='lat')*60*60*24\n", " max_global_avg_pr_daily_test = global_avg_pr_daily_test.max().values\n", "\n", " try:\n", " assert max_global_avg_pr_daily_test < global_avg_pr_daily_old *4\n", " except:\n", " print(f'{f}: very large global average precipitation values')\n", "\n", " np.testing.assert_allclose(pr_min,0, atol=1e-2)\n", " print(f + f': lat_resol={lat_resol}°, lon_resol={lon_resol}°, pr_min = {pr_min:.1f} kg m-2 day-1, pr_max = {pr_max:.1f} kg m-2 day-1, pr_global_avg_max = {max_global_avg_pr_daily_test:.1f} kg m-2 day-1')\n", "\n", " elif 'tas' in f:\n", " tas_min = (dtest.tas.min()-273.15).values\n", " tas_max = (dtest.tas.max()-273.15).values\n", " assert tas_min>-100\n", " assert tas_min<-40\n", " assert tas_max<100\n", " assert tas_max>40\n", " print(f + f': lat_resol={lat_resol}°, lon_resol={lon_resol}°, tas_min = {tas_min:.1f}°C, tas_max = {tas_max:.1f}°C')\n", " print(dtest.time.isel(time=0).values)" ] }, { "cell_type": "code", "execution_count": null, "id": "66f25c50-c50e-4c80-a7df-d416052f711d", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "3ad940ad-c724-414d-9cbb-52c89718487a", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": null, "id": "32bb10fa-d935-4de9-870c-fcb2578799d2", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "markdown", "id": "b4ad3f06-be86-411a-a168-af1fb1b0221d", "metadata": {}, "source": [ "## CMIP6 until 2300" ] }, { "cell_type": "code", "execution_count": 25, "id": "6b62f083-d5c0-4720-ab1b-6b072d3e1b55", "metadata": {}, "outputs": [], "source": [ "files_l_tas = glob.glob(\"/home/www/rchizzola/data_final_redo/tas/*\")\n", "files_l_pr = glob.glob(\"/home/www/rchizzola/data_final_redo/pr/*\")\n", "\n" ] }, { "cell_type": "code", "execution_count": 26, "id": "9375fb53-004b-40e8-88af-0299bda0d4a0", "metadata": {}, "outputs": [], "source": [ "files_l_nico_tas = glob.glob(\"/home/www/rchizzola/data_final_redo/nico_data/tas/*\")\n", "files_l_nico_pr = glob.glob(\"/home/www/rchizzola/data_final_redo/nico_data/pr/*\")" ] }, { "cell_type": "code", "execution_count": 27, "id": "fe60e323-a23e-4803-b32c-7394bd4ae557", "metadata": {}, "outputs": [], "source": [ "files_l_rebekka = files_l_tas + files_l_pr \n", "files_l_nico = files_l_nico_tas+ files_l_nico_pr\n", "files_l = files_l_rebekka + files_l_nico" ] }, { "cell_type": "code", "execution_count": 28, "id": "e8df6e06-3ea3-4641-87bc-52e1df1a92ee", "metadata": {}, "outputs": [], "source": [ "load = True\n", "if load:\n", " gcms_cmip6_2300 = pd.DataFrame()\n", "\n", " for fpath in files_l_rebekka:\n", " if '.nc' in fpath:\n", " fname = fpath.split('/')[-1]\n", " folder_p = fpath[:-len(fname)]\n", " gcm, ssp, realisation, var = fname.split('_')\n", " gcms_cmip6_2300.loc[fname, 'ssp'] = ssp\n", "\n", " gcms_cmip6_2300.loc[fname, 'gcm'] = gcm\n", " gcms_cmip6_2300.loc[fname, 'rea'] = realisation\n", " gcms_cmip6_2300.loc[fname, 'var'] = var.split('.')[0]\n", " gcms_cmip6_2300.loc[fname, 'path'] = fpath\n", " gcms_cmip6_2300.loc[fname, 'fname'] = fname\n", " for fpath in files_l_nico:\n", " if '.nc' in fpath:\n", " fname = fpath.split('/')[-1]\n", " folder_p = fpath[:-len(fname)]\n", " # put that only in the list if not available from Rebekka's files\n", " if fname not in gcms_cmip6_2300.index:\n", " gcm, ssp, realisation, var = fname.split('_')\n", " gcms_cmip6_2300.loc[fname, 'ssp'] = ssp\n", " gcms_cmip6_2300.loc[fname, 'gcm'] = gcm\n", " gcms_cmip6_2300.loc[fname, 'rea'] = realisation\n", " gcms_cmip6_2300.loc[fname, 'var'] = var.split('.')[0]\n", " gcms_cmip6_2300.loc[fname, 'path'] = fpath\n", " gcms_cmip6_2300.loc[fname, 'fname'] = fname\n", "\n", " gcms_cmip6_2300.to_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_vJul2023.csv')\n", "\n", "else:\n", " gcms_cmip6_2300 = pd.read_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_vJul2023.csv', index_col='Unnamed: 0')\n", " #gcms_cmip6_2300 = pd.read_csv('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_update.csv', index_col='Unnamed: 0')" ] }, { "cell_type": "code", "execution_count": 29, "id": "7ea80ddb-c5c0-40bf-a27c-063b3aac99e8", "metadata": {}, "outputs": [], "source": [ "gcms_cmip6_2300.to_html('/home/www/lschuster/runs_oggm_v16/all_gcm_list_2300_vJul2023.html')" ] }, { "cell_type": "code", "execution_count": 165, "id": "bcbecf45-fb87-485e-aaec-eca5e5369d16", "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", " \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", " \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", "
sspgcmreavarpathfname
ACCESS-CM2_ssp126_r1i1p1f1_pr.ncssp126ACCESS-CM2r1i1p1f1pr/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-CM2_ssp126_r1i1p1f1_pr.nc
ACCESS-CM2_ssp126_r1i1p1f1_tas.ncssp126ACCESS-CM2r1i1p1f1tas/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-CM2_ssp126_r1i1p1f1_tas.nc
ACCESS-CM2_ssp585_r1i1p1f1_pr.ncssp585ACCESS-CM2r1i1p1f1pr/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-CM2_ssp585_r1i1p1f1_pr.nc
ACCESS-CM2_ssp585_r1i1p1f1_tas.ncssp585ACCESS-CM2r1i1p1f1tas/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-CM2_ssp585_r1i1p1f1_tas.nc
ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.ncssp126ACCESS-ESM1-5r1i1p1f1pr/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc
ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.ncssp126ACCESS-ESM1-5r1i1p1f1tas/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc
ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.ncssp585ACCESS-ESM1-5r1i1p1f1pr/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc
ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.ncssp585ACCESS-ESM1-5r1i1p1f1tas/home/www/rchizzola/data_final_redo/nico_data/...ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc
CESM2-WACCM_ssp126_r1i1p1f1_pr.ncssp126CESM2-WACCMr1i1p1f1pr/home/www/rchizzola/data_final_redo/nico_data/...CESM2-WACCM_ssp126_r1i1p1f1_pr.nc
CESM2-WACCM_ssp126_r1i1p1f1_tas.ncssp126CESM2-WACCMr1i1p1f1tas/home/www/rchizzola/data_final_redo/nico_data/...CESM2-WACCM_ssp126_r1i1p1f1_tas.nc
CESM2-WACCM_ssp534-over_r1i1p1f1_pr.ncssp534-overCESM2-WACCMr1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/CESM2-W...CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc
CESM2-WACCM_ssp534-over_r1i1p1f1_tas.ncssp534-overCESM2-WACCMr1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/CESM2-...CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc
CESM2-WACCM_ssp585_r1i1p1f1_pr.ncssp585CESM2-WACCMr1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/CESM2-W...CESM2-WACCM_ssp585_r1i1p1f1_pr.nc
CESM2-WACCM_ssp585_r1i1p1f1_tas.ncssp585CESM2-WACCMr1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/CESM2-...CESM2-WACCM_ssp585_r1i1p1f1_tas.nc
CanESM5_ssp126_r1i1p1f1_pr.ncssp126CanESM5r1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/CanESM5...CanESM5_ssp126_r1i1p1f1_pr.nc
CanESM5_ssp126_r1i1p1f1_tas.ncssp126CanESM5r1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/CanESM...CanESM5_ssp126_r1i1p1f1_tas.nc
CanESM5_ssp534-over_r1i1p1f1_pr.ncssp534-overCanESM5r1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/CanESM5...CanESM5_ssp534-over_r1i1p1f1_pr.nc
CanESM5_ssp534-over_r1i1p1f1_tas.ncssp534-overCanESM5r1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/CanESM...CanESM5_ssp534-over_r1i1p1f1_tas.nc
CanESM5_ssp585_r1i1p1f1_pr.ncssp585CanESM5r1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/CanESM5...CanESM5_ssp585_r1i1p1f1_pr.nc
CanESM5_ssp585_r1i1p1f1_tas.ncssp585CanESM5r1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/CanESM...CanESM5_ssp585_r1i1p1f1_tas.nc
IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.ncssp126IPSL-CM6A-LRr1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/IPSL-CM...IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc
IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.ncssp126IPSL-CM6A-LRr1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/IPSL-C...IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc
IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.ncssp534-overIPSL-CM6A-LRr1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/IPSL-CM...IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc
IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.ncssp534-overIPSL-CM6A-LRr1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/IPSL-C...IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc
IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.ncssp585IPSL-CM6A-LRr1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/IPSL-CM...IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc
IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.ncssp585IPSL-CM6A-LRr1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/IPSL-C...IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc
MRI-ESM2-0_ssp126_r1i1p1f1_pr.ncssp126MRI-ESM2-0r1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/MRI-ESM...MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc
MRI-ESM2-0_ssp126_r1i1p1f1_tas.ncssp126MRI-ESM2-0r1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/MRI-ES...MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc
MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.ncssp534-overMRI-ESM2-0r1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/MRI-ESM...MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc
MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.ncssp534-overMRI-ESM2-0r1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/MRI-ES...MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc
MRI-ESM2-0_ssp585_r1i1p1f1_pr.ncssp585MRI-ESM2-0r1i1p1f1pr/home/www/rchizzola/data_final_redo/pr/MRI-ESM...MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc
MRI-ESM2-0_ssp585_r1i1p1f1_tas.ncssp585MRI-ESM2-0r1i1p1f1tas/home/www/rchizzola/data_final_redo/tas/MRI-ES...MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc
\n", "
" ], "text/plain": [ " ssp gcm \\\n", "ACCESS-CM2_ssp126_r1i1p1f1_pr.nc ssp126 ACCESS-CM2 \n", "ACCESS-CM2_ssp126_r1i1p1f1_tas.nc ssp126 ACCESS-CM2 \n", "ACCESS-CM2_ssp585_r1i1p1f1_pr.nc ssp585 ACCESS-CM2 \n", "ACCESS-CM2_ssp585_r1i1p1f1_tas.nc ssp585 ACCESS-CM2 \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc ssp126 ACCESS-ESM1-5 \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc ssp126 ACCESS-ESM1-5 \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc ssp585 ACCESS-ESM1-5 \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc ssp585 ACCESS-ESM1-5 \n", "CESM2-WACCM_ssp126_r1i1p1f1_pr.nc ssp126 CESM2-WACCM \n", "CESM2-WACCM_ssp126_r1i1p1f1_tas.nc ssp126 CESM2-WACCM \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc ssp534-over CESM2-WACCM \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc ssp534-over CESM2-WACCM \n", "CESM2-WACCM_ssp585_r1i1p1f1_pr.nc ssp585 CESM2-WACCM \n", "CESM2-WACCM_ssp585_r1i1p1f1_tas.nc ssp585 CESM2-WACCM \n", "CanESM5_ssp126_r1i1p1f1_pr.nc ssp126 CanESM5 \n", "CanESM5_ssp126_r1i1p1f1_tas.nc ssp126 CanESM5 \n", "CanESM5_ssp534-over_r1i1p1f1_pr.nc ssp534-over CanESM5 \n", "CanESM5_ssp534-over_r1i1p1f1_tas.nc ssp534-over CanESM5 \n", "CanESM5_ssp585_r1i1p1f1_pr.nc ssp585 CanESM5 \n", "CanESM5_ssp585_r1i1p1f1_tas.nc ssp585 CanESM5 \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc ssp126 IPSL-CM6A-LR \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc ssp126 IPSL-CM6A-LR \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc ssp534-over IPSL-CM6A-LR \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc ssp534-over IPSL-CM6A-LR \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc ssp585 IPSL-CM6A-LR \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc ssp585 IPSL-CM6A-LR \n", "MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc ssp126 MRI-ESM2-0 \n", "MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc ssp126 MRI-ESM2-0 \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc ssp534-over MRI-ESM2-0 \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc ssp534-over MRI-ESM2-0 \n", "MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc ssp585 MRI-ESM2-0 \n", "MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc ssp585 MRI-ESM2-0 \n", "\n", " rea var \\\n", "ACCESS-CM2_ssp126_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "ACCESS-CM2_ssp126_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "ACCESS-CM2_ssp585_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "ACCESS-CM2_ssp585_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "CESM2-WACCM_ssp126_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "CESM2-WACCM_ssp126_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "CESM2-WACCM_ssp585_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "CESM2-WACCM_ssp585_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "CanESM5_ssp126_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "CanESM5_ssp126_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "CanESM5_ssp534-over_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "CanESM5_ssp534-over_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "CanESM5_ssp585_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "CanESM5_ssp585_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc r1i1p1f1 pr \n", "MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc r1i1p1f1 tas \n", "\n", " path \\\n", "ACCESS-CM2_ssp126_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-CM2_ssp126_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-CM2_ssp585_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-CM2_ssp585_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "CESM2-WACCM_ssp126_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "CESM2-WACCM_ssp126_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/nico_data/... \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/CESM2-W... \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/CESM2-... \n", "CESM2-WACCM_ssp585_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/CESM2-W... \n", "CESM2-WACCM_ssp585_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/CESM2-... \n", "CanESM5_ssp126_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/CanESM5... \n", "CanESM5_ssp126_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/CanESM... \n", "CanESM5_ssp534-over_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/CanESM5... \n", "CanESM5_ssp534-over_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/CanESM... \n", "CanESM5_ssp585_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/CanESM5... \n", "CanESM5_ssp585_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/CanESM... \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/IPSL-CM... \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/IPSL-C... \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/IPSL-CM... \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/IPSL-C... \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/IPSL-CM... \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/IPSL-C... \n", "MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/MRI-ESM... \n", "MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/MRI-ES... \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/MRI-ESM... \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/MRI-ES... \n", "MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc /home/www/rchizzola/data_final_redo/pr/MRI-ESM... \n", "MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc /home/www/rchizzola/data_final_redo/tas/MRI-ES... \n", "\n", " fname \n", "ACCESS-CM2_ssp126_r1i1p1f1_pr.nc ACCESS-CM2_ssp126_r1i1p1f1_pr.nc \n", "ACCESS-CM2_ssp126_r1i1p1f1_tas.nc ACCESS-CM2_ssp126_r1i1p1f1_tas.nc \n", "ACCESS-CM2_ssp585_r1i1p1f1_pr.nc ACCESS-CM2_ssp585_r1i1p1f1_pr.nc \n", "ACCESS-CM2_ssp585_r1i1p1f1_tas.nc ACCESS-CM2_ssp585_r1i1p1f1_tas.nc \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc \n", "ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc \n", "ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc \n", "CESM2-WACCM_ssp126_r1i1p1f1_pr.nc CESM2-WACCM_ssp126_r1i1p1f1_pr.nc \n", "CESM2-WACCM_ssp126_r1i1p1f1_tas.nc CESM2-WACCM_ssp126_r1i1p1f1_tas.nc \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc \n", "CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc \n", "CESM2-WACCM_ssp585_r1i1p1f1_pr.nc CESM2-WACCM_ssp585_r1i1p1f1_pr.nc \n", "CESM2-WACCM_ssp585_r1i1p1f1_tas.nc CESM2-WACCM_ssp585_r1i1p1f1_tas.nc \n", "CanESM5_ssp126_r1i1p1f1_pr.nc CanESM5_ssp126_r1i1p1f1_pr.nc \n", "CanESM5_ssp126_r1i1p1f1_tas.nc CanESM5_ssp126_r1i1p1f1_tas.nc \n", "CanESM5_ssp534-over_r1i1p1f1_pr.nc CanESM5_ssp534-over_r1i1p1f1_pr.nc \n", "CanESM5_ssp534-over_r1i1p1f1_tas.nc CanESM5_ssp534-over_r1i1p1f1_tas.nc \n", "CanESM5_ssp585_r1i1p1f1_pr.nc CanESM5_ssp585_r1i1p1f1_pr.nc \n", "CanESM5_ssp585_r1i1p1f1_tas.nc CanESM5_ssp585_r1i1p1f1_tas.nc \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc \n", "IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc \n", "IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc \n", "IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc \n", "MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc \n", "MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc \n", "MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc \n", "MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc \n", "MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc " ] }, "execution_count": 165, "metadata": {}, "output_type": "execute_result" } ], "source": [ "gcms_cmip6_2300.sort_values(by='fname')" ] }, { "cell_type": "code", "execution_count": 166, "id": "25f661c9-bb12-4d99-a2f6-92d7fa4d2184", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "16.0" ] }, "execution_count": 166, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# amount of GCMS * SSPs\n", "len(gcms_cmip6_2300)/2" ] }, { "cell_type": "code", "execution_count": 167, "id": "ada0c0ba-9df0-4f92-895e-f6d0817025e1", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array(['CanESM5', 'IPSL-CM6A-LR', 'CESM2-WACCM', 'MRI-ESM2-0',\n", " 'ACCESS-CM2', 'ACCESS-ESM1-5'], dtype=object)" ] }, "execution_count": 167, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# amount of GCMs\n", "gcms_cmip6_2300.gcm.unique()" ] }, { "cell_type": "code", "execution_count": 168, "id": "8da8fbfb-7203-40de-ade4-0782c890983e", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array(['ssp585', 'ssp126', 'ssp534-over'], dtype=object)" ] }, "execution_count": 168, "metadata": {}, "output_type": "execute_result" } ], "source": [ "# amount of SSPs\n", "gcms_cmip6_2300.ssp.unique()" ] }, { "cell_type": "code", "execution_count": 280, "id": "59bc6299-27df-4630-852c-2802d15dfd29", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ssp126 915\n", "ssp534-over 912\n", "ssp585 918\n" ] } ], "source": [ "for ssp in ['ssp126', 'ssp534-over', 'ssp585']:\n", " dt = xr.open_dataset(f'/home/www/lschuster/runs_oggm_v16/runs_2023.3/output/RGI11/run_hydro_w5e5_gcm_merged_endyr2300_CMIP6_MRI-ESM2-0_{ssp}_bc_2000_2019_rgi11_3000_4000.nc')\n", " dt = dt.dropna(dim='rgi_id', how='all')\n", " print(ssp, len(dt.rgi_id.values))" ] }, { "cell_type": "code", "execution_count": 57, "id": "fea396fd-5ad2-4302-8d1e-2a9476b2c779", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'time' (time: 301)>\n",
       "array([2000., 2001., 2002., ..., 2298., 2299., 2300.])\n",
       "Coordinates:\n",
       "  * time            (time) float64 2e+03 2.001e+03 ... 2.299e+03 2.3e+03\n",
       "    hydro_year      (time) int64 2000 2001 2002 2003 ... 2297 2298 2299 2300\n",
       "    hydro_month     (time) int64 4 4 4 4 4 4 4 4 4 4 4 ... 4 4 4 4 4 4 4 4 4 4 4\n",
       "    calendar_year   (time) int64 2000 2001 2002 2003 ... 2297 2298 2299 2300\n",
       "    calendar_month  (time) int64 1 1 1 1 1 1 1 1 1 1 1 ... 1 1 1 1 1 1 1 1 1 1 1\n",
       "Attributes:\n",
       "    description:  Floating year
" ], "text/plain": [ "\n", "array([2000., 2001., 2002., ..., 2298., 2299., 2300.])\n", "Coordinates:\n", " * time (time) float64 2e+03 2.001e+03 ... 2.299e+03 2.3e+03\n", " hydro_year (time) int64 ...\n", " hydro_month (time) int64 ...\n", " calendar_year (time) int64 ...\n", " calendar_month (time) int64 ...\n", "Attributes:\n", " description: Floating year" ] }, "execution_count": 57, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dt.time" ] }, { "cell_type": "markdown", "id": "d304bbfc-4987-49a3-a38d-887aeb399075", "metadata": {}, "source": [ "### Check if format is similar to the old one:" ] }, { "cell_type": "code", "execution_count": 221, "id": "b894055d-6d3b-4d1a-9564-32f32e00ec2f", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "84.27310455590487" ] }, "execution_count": 221, "metadata": {}, "output_type": "execute_result" } ], "source": [ "#xr.open_dataset('/home/www/oggm/cmip6/GCM/CAMS-CSM1-0/CAMS-CSM1-0_ssp119_r2i1p1f1_tas.nc').time\n", "dold_pr = xr.open_dataset('/home/www/oggm/cmip6/GCM/GFDL-ESM4/GFDL-ESM4_ssp585_r1i1p1f1_pr.nc')\n", "#dold_tas = xr.open_dataset('/home/www/oggm/cmip6/GCM/GFDL-ESM4/GFDL-ESM4_ssp119_r1i1p1f1_tas.nc')\n", "pr_max_other_cmip6 = dold_pr.pr.max().values*60*60*24 # kg m-2 day-1\n", "pr_max_other_cmip6" ] }, { "cell_type": "code", "execution_count": 254, "id": "917b562a-802b-422b-a710-22fe46904fdc", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.Dataset>\n",
       "Dimensions:    (time: 3012, bnds: 2, lon: 288, lat: 180)\n",
       "Coordinates:\n",
       "  * time       (time) object 1850-01-16 12:00:00 ... 2100-12-16 12:00:00\n",
       "  * lon        (lon) float64 0.625 1.875 3.125 4.375 ... 355.6 356.9 358.1 359.4\n",
       "  * lat        (lat) float64 -89.5 -88.5 -87.5 -86.5 ... 86.5 87.5 88.5 89.5\n",
       "Dimensions without coordinates: bnds\n",
       "Data variables:\n",
       "    time_bnds  (time, bnds) object 1850-01-01 00:00:00 ... 2101-01-01 00:00:00\n",
       "    lon_bnds   (lon, bnds) float64 0.0 1.25 1.25 2.5 ... 357.5 358.8 358.8 360.0\n",
       "    lat_bnds   (lat, bnds) float64 -90.0 -89.0 -89.0 -88.0 ... 89.0 89.0 90.0\n",
       "    pr         (time, lat, lon) float32 1.209e-06 1.209e-06 ... 1.26e-05\n",
       "Attributes: (12/48)\n",
       "    CDI:                    Climate Data Interface version 1.9.5 (http://mpim...\n",
       "    history:                Wed Apr 08 14:40:22 2020: cdo mergetime /scratch_...\n",
       "    source:                 GFDL-ESM4 (2018):\\natmos: GFDL-AM4.1 (Cubed-spher...\n",
       "    institution:            National Oceanic and Atmospheric Administration, ...\n",
       "    Conventions:            CF-1.7 CMIP-6.0 UGRID-1.0\n",
       "    external_variables:     areacella\n",
       "    ...                     ...\n",
       "    tracking_id:            hdl:21.14100/29468e1c-b66b-40d6-92aa-9851fba964ee\n",
       "    variable_id:            pr\n",
       "    variant_info:           N/A\n",
       "    references:             see further_info_url attribute\n",
       "    variant_label:          r1i1p1f1\n",
       "    CDO:                    Climate Data Operators version 1.9.5 (http://mpim...
" ], "text/plain": [ "\n", "Dimensions: (time: 3012, bnds: 2, lon: 288, lat: 180)\n", "Coordinates:\n", " * time (time) object 1850-01-16 12:00:00 ... 2100-12-16 12:00:00\n", " * lon (lon) float64 0.625 1.875 3.125 4.375 ... 355.6 356.9 358.1 359.4\n", " * lat (lat) float64 -89.5 -88.5 -87.5 -86.5 ... 86.5 87.5 88.5 89.5\n", "Dimensions without coordinates: bnds\n", "Data variables:\n", " time_bnds (time, bnds) object 1850-01-01 00:00:00 ... 2101-01-01 00:00:00\n", " lon_bnds (lon, bnds) float64 0.0 1.25 1.25 2.5 ... 357.5 358.8 358.8 360.0\n", " lat_bnds (lat, bnds) float64 -90.0 -89.0 -89.0 -88.0 ... 89.0 89.0 90.0\n", " pr (time, lat, lon) float32 1.209e-06 1.209e-06 ... 1.26e-05\n", "Attributes: (12/48)\n", " CDI: Climate Data Interface version 1.9.5 (http://mpim...\n", " history: Wed Apr 08 14:40:22 2020: cdo mergetime /scratch_...\n", " source: GFDL-ESM4 (2018):\\natmos: GFDL-AM4.1 (Cubed-spher...\n", " institution: National Oceanic and Atmospheric Administration, ...\n", " Conventions: CF-1.7 CMIP-6.0 UGRID-1.0\n", " external_variables: areacella\n", " ... ...\n", " tracking_id: hdl:21.14100/29468e1c-b66b-40d6-92aa-9851fba964ee\n", " variable_id: pr\n", " variant_info: N/A\n", " references: see further_info_url attribute\n", " variant_label: r1i1p1f1\n", " CDO: Climate Data Operators version 1.9.5 (http://mpim..." ] }, "execution_count": 254, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dold_pr" ] }, { "cell_type": "code", "execution_count": 274, "id": "df340362-ab36-421c-a397-735b8ef48e81", "metadata": {}, "outputs": [], "source": [ "weight = np.cos(np.deg2rad(dold_pr.lat.astype(np.float64))).clip(0)\n", "weight = dold_pr.pr.isel(time=0) * 0. + weight\n", "dold_pr['weight'] = (('lat', 'lon'), weight.values / weight.sum().values)\n", "global_avg_pr_daily_old = ((dold_pr.weight * dold_pr.pr).sum(dim='lon').sum(dim='lat')*60*60*24).max().values" ] }, { "cell_type": "code", "execution_count": 275, "id": "a83a25dc-026f-472b-a629-496199463bc8", "metadata": {}, "outputs": [ { "data": { "text/plain": [ "array(3.19732436)" ] }, "execution_count": 275, "metadata": {}, "output_type": "execute_result" } ], "source": [] }, { "cell_type": "code", "execution_count": 260, "id": "edc55928-237d-4844-8382-cff18630fb8c", "metadata": {}, "outputs": [ { "data": { "text/html": [ "
\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "\n", "
<xarray.DataArray 'weight' ()>\n",
       "array(1.)
" ], "text/plain": [ "\n", "array(1.)" ] }, "execution_count": 260, "metadata": {}, "output_type": "execute_result" } ], "source": [ "dold_pr['weight'].sum()" ] }, { "cell_type": "code", "execution_count": 263, "id": "3119533f-e063-4cb8-bee9-2112c5b6d25d", "metadata": {}, "outputs": [], "source": [] }, { "cell_type": "code", "execution_count": 250, "id": "9beff0d8-125f-4057-99a9-f06a1daf4a5d", "metadata": {}, "outputs": [], "source": [ "dold_pr_lon_mean_daily = dold_pr.pr.mean(dim='lon')" ] }, { "cell_type": "code", "execution_count": 224, "id": "854adeed-ad64-49fe-9445-b0707a3e4afc", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "500 km 5412 64 128 /home/www/rchizzola/data_final_redo/tas/CanESM5_ssp585_r1i1p1f1_tas.nc\n", "250 km 5412 143 144 /home/www/rchizzola/data_final_redo/tas/IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc\n", "500 km 5412 64 128 /home/www/rchizzola/data_final_redo/tas/CanESM5_ssp126_r1i1p1f1_tas.nc\n", "500 km 5412 64 128 /home/www/rchizzola/data_final_redo/tas/CanESM5_ssp534-over_r1i1p1f1_tas.nc\n", "100 km 5400 192 288 /home/www/rchizzola/data_final_redo/tas/CESM2-WACCM_ssp585_r1i1p1f1_tas.nc\n", "250 km 5412 143 144 /home/www/rchizzola/data_final_redo/tas/IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc\n", "250 km 5412 143 144 /home/www/rchizzola/data_final_redo/tas/IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc\n", "100 km 5412 160 320 /home/www/rchizzola/data_final_redo/tas/MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc\n", "100 km 5400 192 288 /home/www/rchizzola/data_final_redo/tas/CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc\n", "100 km 5412 160 320 /home/www/rchizzola/data_final_redo/tas/MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc\n", "100 km 5412 160 320 /home/www/rchizzola/data_final_redo/tas/MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc\n", "100 km 5400 192 288 /home/www/rchizzola/data_final_redo/nico_data/tas/CESM2-WACCM_ssp126_r1i1p1f1_tas.nc\n", "250 km 5412 144 192 /home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-CM2_ssp126_r1i1p1f1_tas.nc\n", "250 km 5412 145 192 /home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc\n", "250 km 5412 144 192 /home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-CM2_ssp585_r1i1p1f1_tas.nc\n", "250 km 5412 145 192 /home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc\n" ] } ], "source": [ "for p in gcms_cmip6_2300.loc[gcms_cmip6_2300['var']=='tas'].path:\n", " dtest = xr.open_dataset(p, use_cftime=True) \n", " time, lat, lon = np.shape(dtest.tas)\n", " # either from 1850 until 2300 or until 2999\n", " try:\n", " assert time==5412\n", " except:\n", " assert time==5400\n", " print(dtest.attrs['nominal_resolution'], time, lat,lon, p)" ] }, { "cell_type": "code", "execution_count": 225, "id": "e4858d2a-3139-4c44-becf-f117e5604702", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "100 km 5412 160 320 /home/www/rchizzola/data_final_redo/pr/MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "250 km 5412 143 144 /home/www/rchizzola/data_final_redo/pr/IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "100 km 5412 160 320 /home/www/rchizzola/data_final_redo/pr/MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "100 km 5412 160 320 /home/www/rchizzola/data_final_redo/pr/MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "500 km 5412 64 128 /home/www/rchizzola/data_final_redo/pr/CanESM5_ssp534-over_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "100 km 5400 192 288 /home/www/rchizzola/data_final_redo/pr/CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc 1850-01-15 12:00:00\n", "250 km 5412 143 144 /home/www/rchizzola/data_final_redo/pr/IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "500 km 5412 64 128 /home/www/rchizzola/data_final_redo/pr/CanESM5_ssp585_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "250 km 5412 143 144 /home/www/rchizzola/data_final_redo/pr/IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "100 km 5400 192 288 /home/www/rchizzola/data_final_redo/pr/CESM2-WACCM_ssp585_r1i1p1f1_pr.nc 1850-01-15 12:00:00\n", "500 km 5412 64 128 /home/www/rchizzola/data_final_redo/pr/CanESM5_ssp126_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "250 km 5412 144 192 /home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp585_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "250 km 5412 145 192 /home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "250 km 5412 144 192 /home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp126_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "250 km 5412 145 192 /home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc 1850-01-16 12:00:00\n", "100 km 5400 192 288 /home/www/rchizzola/data_final_redo/nico_data/pr/CESM2-WACCM_ssp126_r1i1p1f1_pr.nc 1850-01-15 12:00:00\n" ] } ], "source": [ "for p in gcms_cmip6_2300.loc[gcms_cmip6_2300['var']=='pr'].path:\n", " dtest = xr.open_dataset(p, use_cftime=True) \n", " time, lat, lon = np.shape(dtest.pr)\n", " # either from 1850 until 2300 or until 2999\n", " try:\n", " assert time==5412\n", " except:\n", " assert time==5400\n", " print(dtest.attrs['nominal_resolution'], time, lat,lon, p, dtest.time[0].values )" ] }, { "cell_type": "markdown", "id": "d706970f-40f8-4a64-a206-74d223f6308f", "metadata": {}, "source": [ "## Check units:" ] }, { "cell_type": "code", "execution_count": 173, "id": "0785e235-c349-4d08-8c7c-905b809a072d", "metadata": {}, "outputs": [], "source": [ "test_convert = xr.open_dataset(files_l[-4],use_cftime=True)\n", "test_convert.time[0].values #'cftime.DatetimeGregorian'\n", "#test_convert = xr.open_dataset(files_l[2])\n", "\n", "#test_convert['time'] = xr.open_dataset(files_l[0]).time.values\n", "#test_convert = test_convert.drop('time_bounds')\n", "#test_convert.to_netcdf('test_convert2.nc')" ] }, { "cell_type": "code", "execution_count": 276, "id": "0037dc54-74a9-437a-ad38-83033b81590c", "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "/home/www/rchizzola/data_final_redo/tas/CanESM5_ssp585_r1i1p1f1_tas.nc: lat_resol=2.8°, lon_resol=2.8°, tas_min = -93.1°C, tas_max = 66.6°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, tas_min = -84.6°C, tas_max = 62.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/CanESM5_ssp126_r1i1p1f1_tas.nc: lat_resol=2.8°, lon_resol=2.8°, tas_min = -93.1°C, tas_max = 45.3°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/CanESM5_ssp534-over_r1i1p1f1_tas.nc: lat_resol=2.8°, lon_resol=2.8°, tas_min = -93.1°C, tas_max = 46.3°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/CESM2-WACCM_ssp585_r1i1p1f1_tas.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -69.2°C, tas_max = 61.7°C\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, tas_min = -84.6°C, tas_max = 42.5°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/IPSL-CM6A-LR_ssp534-over_r1i1p1f1_tas.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, tas_min = -84.6°C, tas_max = 43.2°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc: lat_resol=1.1°, lon_resol=1.1°, tas_min = -68.2°C, tas_max = 47.7°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/CESM2-WACCM_ssp534-over_r1i1p1f1_tas.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -69.2°C, tas_max = 46.4°C\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc: lat_resol=1.1°, lon_resol=1.1°, tas_min = -68.2°C, tas_max = 57.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/tas/MRI-ESM2-0_ssp534-over_r1i1p1f1_tas.nc: lat_resol=1.1°, lon_resol=1.1°, tas_min = -68.2°C, tas_max = 47.5°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc: lat_resol=1.1°, lon_resol=1.1°, pr_min = 0.0 kg m-2 day-1, pr_max = 117.8 kg m-2 day-1, pr_global_avg_max = 3.7 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/IPSL-CM6A-LR_ssp534-over_r1i1p1f1_pr.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, pr_min = 0.0 kg m-2 day-1, pr_max = 157.9 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc: lat_resol=1.1°, lon_resol=1.1°, pr_min = 0.0 kg m-2 day-1, pr_max = 70.4 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/MRI-ESM2-0_ssp534-over_r1i1p1f1_pr.nc: lat_resol=1.1°, lon_resol=1.1°, pr_min = 0.0 kg m-2 day-1, pr_max = 73.3 kg m-2 day-1, pr_global_avg_max = 3.3 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/CanESM5_ssp534-over_r1i1p1f1_pr.nc: lat_resol=2.8°, lon_resol=2.8°, pr_min = 0.0 kg m-2 day-1, pr_max = 77.6 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/CESM2-WACCM_ssp534-over_r1i1p1f1_pr.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = 0.0 kg m-2 day-1, pr_max = 110.4 kg m-2 day-1, pr_global_avg_max = 3.3 kg m-2 day-1\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, pr_min = 0.0 kg m-2 day-1, pr_max = 282.8 kg m-2 day-1, pr_global_avg_max = 3.8 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/CanESM5_ssp585_r1i1p1f1_pr.nc: lat_resol=2.8°, lon_resol=2.8°, pr_min = 0.0 kg m-2 day-1, pr_max = 111.2 kg m-2 day-1, pr_global_avg_max = 4.0 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, pr_min = 0.0 kg m-2 day-1, pr_max = 161.9 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/CESM2-WACCM_ssp585_r1i1p1f1_pr.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = 0.0 kg m-2 day-1, pr_max = 289.8 kg m-2 day-1, pr_global_avg_max = 3.8 kg m-2 day-1\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/pr/CanESM5_ssp126_r1i1p1f1_pr.nc: lat_resol=2.8°, lon_resol=2.8°, pr_min = 0.0 kg m-2 day-1, pr_max = 82.7 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/CESM2-WACCM_ssp126_r1i1p1f1_tas.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -69.2°C, tas_max = 45.7°C\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/CanESM5_ssp585_r1i1p1f1_tas.nc: lat_resol=2.8°, lon_resol=2.8°, tas_min = -93.1°C, tas_max = 66.6°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/IPSL-CM6A-LR_ssp126_r1i1p1f1_tas.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, tas_min = -84.6°C, tas_max = 42.5°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/MRI-ESM2-0_ssp126_r1i1p1f1_tas.nc: lat_resol=1.1°, lon_resol=1.1°, tas_min = -68.2°C, tas_max = 47.7°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/IPSL-CM6A-LR_ssp585_r1i1p1f1_tas.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, tas_min = -84.6°C, tas_max = 62.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/CanESM5_ssp126_r1i1p1f1_tas.nc: lat_resol=2.8°, lon_resol=2.8°, tas_min = -93.1°C, tas_max = 45.3°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/CESM2-WACCM_ssp585_r1i1p1f1_tas.nc: lat_resol=0.9°, lon_resol=1.2°, tas_min = -69.2°C, tas_max = 61.7°C\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/MRI-ESM2-0_ssp585_r1i1p1f1_tas.nc: lat_resol=1.1°, lon_resol=1.1°, tas_min = -68.2°C, tas_max = 57.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-CM2_ssp126_r1i1p1f1_tas.nc: lat_resol=1.2°, lon_resol=1.9°, tas_min = -81.9°C, tas_max = 45.4°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-ESM1-5_ssp585_r1i1p1f1_tas.nc: lat_resol=1.2°, lon_resol=1.9°, tas_min = -70.4°C, tas_max = 55.2°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-CM2_ssp585_r1i1p1f1_tas.nc: lat_resol=1.2°, lon_resol=1.9°, tas_min = -81.9°C, tas_max = 59.5°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/tas/ACCESS-ESM1-5_ssp126_r1i1p1f1_tas.nc: lat_resol=1.2°, lon_resol=1.9°, tas_min = -70.4°C, tas_max = 43.7°C\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp585_r1i1p1f1_pr.nc: very large precipitation values\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp585_r1i1p1f1_pr.nc: lat_resol=1.2°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 1103.5 kg m-2 day-1, pr_global_avg_max = 4.1 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc: very large precipitation values\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-ESM1-5_ssp585_r1i1p1f1_pr.nc: lat_resol=1.2°, lon_resol=1.9°, pr_min = -0.0 kg m-2 day-1, pr_max = 406.1 kg m-2 day-1, pr_global_avg_max = 3.8 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp126_r1i1p1f1_pr.nc: very large precipitation values\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp126_r1i1p1f1_pr.nc: lat_resol=1.2°, lon_resol=1.9°, pr_min = 0.0 kg m-2 day-1, pr_max = 501.4 kg m-2 day-1, pr_global_avg_max = 3.5 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-ESM1-5_ssp126_r1i1p1f1_pr.nc: lat_resol=1.2°, lon_resol=1.9°, pr_min = -0.0 kg m-2 day-1, pr_max = 100.1 kg m-2 day-1, pr_global_avg_max = 3.5 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/CESM2-WACCM_ssp585_r1i1p1f1_pr.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = 0.0 kg m-2 day-1, pr_max = 289.8 kg m-2 day-1, pr_global_avg_max = 3.8 kg m-2 day-1\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/MRI-ESM2-0_ssp126_r1i1p1f1_pr.nc: lat_resol=1.1°, lon_resol=1.1°, pr_min = 0.0 kg m-2 day-1, pr_max = 70.4 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/IPSL-CM6A-LR_ssp126_r1i1p1f1_pr.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, pr_min = 0.0 kg m-2 day-1, pr_max = 161.9 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/CanESM5_ssp585_r1i1p1f1_pr.nc: lat_resol=2.8°, lon_resol=2.8°, pr_min = 0.0 kg m-2 day-1, pr_max = 111.2 kg m-2 day-1, pr_global_avg_max = 4.0 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/MRI-ESM2-0_ssp585_r1i1p1f1_pr.nc: lat_resol=1.1°, lon_resol=1.1°, pr_min = 0.0 kg m-2 day-1, pr_max = 117.8 kg m-2 day-1, pr_global_avg_max = 3.7 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/CESM2-WACCM_ssp126_r1i1p1f1_pr.nc: lat_resol=0.9°, lon_resol=1.2°, pr_min = 0.0 kg m-2 day-1, pr_max = 126.6 kg m-2 day-1, pr_global_avg_max = 3.3 kg m-2 day-1\n", "1850-01-15 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/CanESM5_ssp126_r1i1p1f1_pr.nc: lat_resol=2.8°, lon_resol=2.8°, pr_min = 0.0 kg m-2 day-1, pr_max = 82.7 kg m-2 day-1, pr_global_avg_max = 3.2 kg m-2 day-1\n", "1850-01-16 12:00:00\n", "/home/www/rchizzola/data_final_redo/nico_data/pr/IPSL-CM6A-LR_ssp585_r1i1p1f1_pr.nc: lat_resol=1.2999999523162842°, lon_resol=2.5°, pr_min = 0.0 kg m-2 day-1, pr_max = 282.8 kg m-2 day-1, pr_global_avg_max = 3.8 kg m-2 day-1\n", "1850-01-16 12:00:00\n" ] } ], "source": [ "pd_max_prcp_l = []\n", "for f in files_l:\n", " if '.nc' in f:\n", " dtest = xr.open_dataset(f,use_cftime=True) \n", " # Check resolution:\n", " lat_resol = np.unique((dtest.lat.values[1:] - dtest.lat.values[:-1]).round(1))[0]\n", " lon_resol = np.unique((dtest.lon.values[1:] - dtest.lon.values[:-1]).round(1))[0]\n", " assert lon_resol<3\n", " assert lat_resol<3\n", " \n", " # Check units\n", " if 'pr' in f:\n", " pr_max = (dtest.pr.max()*60*60*24).values # kg m-2 day-1\n", " pr_min = (dtest.pr.min()*60*60*24).values\n", " try:\n", " assert pr_max=dold_pr.pr.max()*4).pr\n", " dmax = dmax.dropna(dim='lon', how='all').dropna(dim='lat', how='all').dropna(dim='time', how='all')\n", " _pd_max_prcp = (dmax*60*60*24).to_dataframe('daily_prcp_max').dropna()\n", " _pd_max_prcp['file'] = f\n", " pd_max_prcp_l.append(_pd_max_prcp)\n", " # global avg. precipitation \n", " weight = np.cos(np.deg2rad(dtest.lat.astype(np.float64))).clip(0)\n", " weight = dtest.pr.isel(time=0) * 0. + weight\n", " dtest['weight'] = (('lat', 'lon'), weight.values / weight.sum().values)\n", " global_avg_pr_daily_test = (dtest.weight * dtest.pr).sum(dim='lon').sum(dim='lat')*60*60*24\n", " max_global_avg_pr_daily_test = global_avg_pr_daily_test.max().values\n", "\n", " try:\n", " assert max_global_avg_pr_daily_test < global_avg_pr_daily_old *4\n", " except:\n", " print(f'{f}: very large global average precipitation values')\n", "\n", " np.testing.assert_allclose(pr_min,0, atol=1e-10)\n", " print(f + f': lat_resol={lat_resol}°, lon_resol={lon_resol}°, pr_min = {pr_min:.1f} kg m-2 day-1, pr_max = {pr_max:.1f} kg m-2 day-1, pr_global_avg_max = {max_global_avg_pr_daily_test:.1f} kg m-2 day-1')\n", "\n", " elif 'tas' in f:\n", " tas_min = (dtest.tas.min()-273.15).values\n", " tas_max = (dtest.tas.max()-273.15).values\n", " assert tas_min>-100\n", " assert tas_min<-40\n", " assert tas_max<100\n", " assert tas_max>40\n", " print(f + f': lat_resol={lat_resol}°, lon_resol={lon_resol}°, tas_min = {tas_min:.1f}°C, tas_max = {tas_max:.1f}°C')\n", " print(dtest.time.isel(time=0).values)" ] }, { "cell_type": "code", "execution_count": 239, "id": "a613f8cf-1a01-4236-83e5-bcf8447d9026", "metadata": {}, "outputs": [], "source": [ "pd_max_prcp = pd.concat(pd_max_prcp_l)" ] }, { "cell_type": "markdown", "id": "7adac452-7e42-4e5d-a8f8-5d0204dead58", "metadata": {}, "source": [ "The **ACCESS** GCMs have sometimes very large precipitation values. Not sure if this is an error. " ] }, { "cell_type": "code", "execution_count": 246, "id": "16fad317-c0c3-4ad4-b6d0-1fdc227dc685", "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", "
daily_prcp_maxfile
timelatlon
2086-09-16 00:00:0014.375265.3125493.053680/home/www/rchizzola/data_final_redo/nico_data/...
2091-11-16 00:00:0010.62575.9375467.192627/home/www/rchizzola/data_final_redo/nico_data/...
2095-10-16 12:00:008.125276.5625421.974792/home/www/rchizzola/data_final_redo/nico_data/...
2099-05-16 12:00:000.625336.5625431.707214/home/www/rchizzola/data_final_redo/nico_data/...
2099-10-16 12:00:006.875278.4375337.138184/home/www/rchizzola/data_final_redo/nico_data/...
...............
2298-09-16 00:00:0020.62592.8125455.371948/home/www/rchizzola/data_final_redo/nico_data/...
2299-01-16 12:00:00-4.375149.0625424.989990/home/www/rchizzola/data_final_redo/nico_data/...
2299-02-15 00:00:00-5.625147.1875678.665710/home/www/rchizzola/data_final_redo/nico_data/...
-4.375145.3125573.085571/home/www/rchizzola/data_final_redo/nico_data/...
2299-09-16 00:00:008.125107.8125354.976074/home/www/rchizzola/data_final_redo/nico_data/...
\n", "

271 rows × 2 columns

\n", "
" ], "text/plain": [ " daily_prcp_max \\\n", "time lat lon \n", "2086-09-16 00:00:00 14.375 265.3125 493.053680 \n", "2091-11-16 00:00:00 10.625 75.9375 467.192627 \n", "2095-10-16 12:00:00 8.125 276.5625 421.974792 \n", "2099-05-16 12:00:00 0.625 336.5625 431.707214 \n", "2099-10-16 12:00:00 6.875 278.4375 337.138184 \n", "... ... \n", "2298-09-16 00:00:00 20.625 92.8125 455.371948 \n", "2299-01-16 12:00:00 -4.375 149.0625 424.989990 \n", "2299-02-15 00:00:00 -5.625 147.1875 678.665710 \n", " -4.375 145.3125 573.085571 \n", "2299-09-16 00:00:00 8.125 107.8125 354.976074 \n", "\n", " file \n", "time lat lon \n", "2086-09-16 00:00:00 14.375 265.3125 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2091-11-16 00:00:00 10.625 75.9375 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2095-10-16 12:00:00 8.125 276.5625 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2099-05-16 12:00:00 0.625 336.5625 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2099-10-16 12:00:00 6.875 278.4375 /home/www/rchizzola/data_final_redo/nico_data/... \n", "... ... \n", "2298-09-16 00:00:00 20.625 92.8125 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2299-01-16 12:00:00 -4.375 149.0625 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2299-02-15 00:00:00 -5.625 147.1875 /home/www/rchizzola/data_final_redo/nico_data/... \n", " -4.375 145.3125 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2299-09-16 00:00:00 8.125 107.8125 /home/www/rchizzola/data_final_redo/nico_data/... \n", "\n", "[271 rows x 2 columns]" ] }, "execution_count": 246, "metadata": {}, "output_type": "execute_result" } ], "source": [ "pd_max_prcp.loc[pd_max_prcp.file =='/home/www/rchizzola/data_final_redo/nico_data/pr/ACCESS-CM2_ssp585_r1i1p1f1_pr.nc']" ] }, { "cell_type": "code", "execution_count": 244, "id": "3e47c04a-57b0-4ec0-a34e-3ac26b3eb9a9", "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", "
daily_prcp_maxfile
timelatlon
2144-12-16 12:00:00-1.875119.06251103.514771/home/www/rchizzola/data_final_redo/nico_data/...
2252-04-16 00:00:00-4.375145.31251017.436768/home/www/rchizzola/data_final_redo/nico_data/...
2270-09-16 00:00:0015.625265.31251059.056641/home/www/rchizzola/data_final_redo/nico_data/...
\n", "
" ], "text/plain": [ " daily_prcp_max \\\n", "time lat lon \n", "2144-12-16 12:00:00 -1.875 119.0625 1103.514771 \n", "2252-04-16 00:00:00 -4.375 145.3125 1017.436768 \n", "2270-09-16 00:00:00 15.625 265.3125 1059.056641 \n", "\n", " file \n", "time lat lon \n", "2144-12-16 12:00:00 -1.875 119.0625 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2252-04-16 00:00:00 -4.375 145.3125 /home/www/rchizzola/data_final_redo/nico_data/... \n", "2270-09-16 00:00:00 15.625 265.3125 /home/www/rchizzola/data_final_redo/nico_data/... " ] }, "execution_count": 244, "metadata": {}, "output_type": "execute_result" } ], "source": [ "pd_max_prcp[pd_max_prcp.daily_prcp_max>1000]" ] }, { "cell_type": "code", "execution_count": null, "id": "830931ea-d851-43ee-96c4-882b73b484c2", "metadata": {}, "outputs": [], "source": [] } ], "metadata": { "kernelspec": { "display_name": "Python 3 (ipykernel)", "language": "python", "name": "python3" }, "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.10.12" } }, "nbformat": 4, "nbformat_minor": 5 }