Plot_Event_hours_timeseries.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e88a707a",
   "metadata": {},
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import scipy.io as scio\n",
    "import matplotlib.pyplot as plt\n",
    "from scipy.signal import savgol_filter\n",
    "import matplotlib as mpl\n",
    "import cartopy.feature as cfeature\n",
    "import cartopy.crs as ccrs\n",
    "import netCDF4 as nc\n",
    "import numpy as np\n",
    "from Load_data import Data_from_nc\n",
    "import xarray as xr\n",
    "import math\n",
    "import gc\n",
    "\n",
    "def get_data_from_mat(file, variable):\n",
    "    var_dict = scio.loadmat(file)\n",
    "    var = var_dict[variable]\n",
    "    return var\n",
    "\n",
    "def get_data_from_mat_for_calcu(file, variable):\n",
    "    var_dict = scio.loadmat(file)\n",
    "    var = var_dict[variable]\n",
    "    var = var[:, 29:]\n",
    "    return var.T\n",
    "    \n",
    "def plot_line(ax, index, data_mean, data_min, data_max, title, title1, label, color, ylabel, xlabel, marker):\n",
    "    ax.text(-0.05, 1.05, index, color='dimgrey', fontsize=12, transform=ax.transAxes, weight='bold')\n",
    "    ax.plot(range(1901, 2015), data_mean, label = label, color = color, linewidth = 2, marker = marker, markersize=1.2, alpha=0.8)\n",
    "    ax.fill_between(range(1901, 2015), data_min, data_max, color = color, alpha=0.2)\n",
    "    plt.title(title, loc='right')\n",
    "    plt.title(title1, loc='left')\n",
    "    plt.ylabel(ylabel, fontsize = 14)\n",
    "    plt.xlabel(xlabel, fontsize = 14)\n",
    "    plt.yticks(fontsize = 12)\n",
    "    plt.xticks(fontsize=12)\n",
    "    plt.xlim(1901,2015)\n",
    "    plt.legend(loc = 'upper left')\n",
    "    \n",
    "    \n",
    "    \n",
    "ar6_region = get_data_from_mat_for_calcu('/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/ar6_region.mat','ar6_region')\n",
    "irr_diff = get_data_from_mat_for_calcu('/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/irr_diff_out.mat', 'irr_diff_out')\n",
    "\n",
    "AREA = get_data_from_mat('/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/input_data/AREA.mat', 'AREA') # Area for CLM_data\n",
    "AREA = AREA[:, 29:]\n",
    "AREA = AREA.T\n",
    "\n",
    "AREA[irr_diff < 40] = 0\n",
    "area_more_40 = np.sum(AREA)\n",
    "\n",
    "\n",
    "\n",
    "AREA = get_data_from_mat('/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/input_data/AREA.mat', 'AREA') # Area for CLM_data\n",
    "AREA = AREA[:, 29:]\n",
    "AREA = AREA.T\n",
    "\n",
    "AREA[irr_diff < 30] = 0\n",
    "AREA[irr_diff >= 40] = 0\n",
    "area_30_40 = np.sum(AREA)\n",
    "\n",
    "\n",
    "\n",
    "AREA = get_data_from_mat('/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/input_data/AREA.mat', 'AREA') # Area for CLM_data\n",
    "AREA = AREA[:, 29:]\n",
    "AREA = AREA.T\n",
    "\n",
    "AREA[irr_diff < 20] = 0\n",
    "AREA[irr_diff >= 30] = 0\n",
    "area_20_30 = np.sum(AREA)\n",
    "\n",
    "\n",
    "\n",
    "AREA = get_data_from_mat('/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/input_data/AREA.mat', 'AREA') # Area for CLM_data\n",
    "AREA = AREA[:, 29:]\n",
    "AREA = AREA.T\n",
    "\n",
    "AREA[irr_diff < 0.1] = 0\n",
    "AREA[irr_diff >= 20] = 0\n",
    "area_00_20 = np.sum(AREA)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "4cd6926c",
   "metadata": {},
   "outputs": [],
   "source": [
    "def all_total(metric, var_show, ylim1, ylim2, ylim3):\n",
    "    \n",
    "    fig = plt.figure(figsize=(12, 9),dpi=300)\n",
    "    fig.subplots_adjust(hspace=0.4, wspace=0.2, left = 0.05, right = 0.95, top = 0.95, bottom = 0.05)\n",
    "    \n",
    "    var = metric\n",
    "    percentile = 0.99\n",
    "    irr_low = 0.01\n",
    "    irr_high = 20\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    \n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(341, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'a'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_00_20,15,2), savgol_filter((data_NOI_25)/area_00_20,15,2), savgol_filter((data_NOI_75)/area_00_20,15,2), r'$\\Delta$irr: ' +  '0%-20%', var_show + ' 99.0%', '1901irr', 'orange', 'hours/yr', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_00_20,15,2), savgol_filter((data_IRR_25)/area_00_20,15,2), savgol_filter((data_IRR_75)/area_00_20,15,2), r'$\\Delta$irr: ' +  '0%-20%', var_show + ' 99.0%', 'tranirr', 'dodgerblue', 'hours/yr', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim1])\n",
    "    \n",
    "    percentile = 0.995\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(345, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'e'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_00_20,15,2), savgol_filter((data_NOI_25)/area_00_20,15,2), savgol_filter((data_NOI_75)/area_00_20,15,2), r'$\\Delta$irr: ' +  '0%-20%', var_show + ' 99.5%', '1901irr', 'orange', 'hours/yr', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_00_20,15,2), savgol_filter((data_IRR_25)/area_00_20,15,2), savgol_filter((data_IRR_75)/area_00_20,15,2), r'$\\Delta$irr: ' +  '0%-20%', var_show + ' 99.5%', 'tranirr', 'dodgerblue', 'hours/yr', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim2])\n",
    "    percentile = 0.999\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(349, frameon=True)\n",
    "    \n",
    "    ax1.spines[\"top\"].set_visible(True)\n",
    "    ax1.spines[\"right\"].set_visible(True)\n",
    "    index = 'i'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_00_20,15,2), savgol_filter((data_NOI_25)/area_00_20,15,2), savgol_filter((data_NOI_75)/area_00_20,15,2), r'$\\Delta$irr: ' +  '0%-20%', var_show + ' 99.9%', '1901irr', 'orange', 'hours/yr', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_00_20,15,2), savgol_filter((data_IRR_25)/area_00_20,15,2), savgol_filter((data_IRR_75)/area_00_20,15,2), r'$\\Delta$irr: ' +  '0%-20%', var_show + ' 99.9%', 'tranirr', 'dodgerblue', 'hours/yr', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim3])\n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    percentile = 0.99\n",
    "    irr_low = 20\n",
    "    irr_high = 30\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(342, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'b'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_20_30,15,2), savgol_filter((data_NOI_25)/area_20_30,15,2), savgol_filter((data_NOI_75)/area_20_30,15,2), r'$\\Delta$irr: ' +  '20%-30%', var_show + ' 99.0%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_20_30,15,2), savgol_filter((data_IRR_25)/area_20_30,15,2), savgol_filter((data_IRR_75)/area_20_30,15,2), r'$\\Delta$irr: ' +  '20%-30%', var_show + ' 99.0%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim1])\n",
    "    \n",
    "    percentile = 0.995\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(346, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'f'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_20_30,15,2), savgol_filter((data_NOI_25)/area_20_30,15,2), savgol_filter((data_NOI_75)/area_20_30,15,2), r'$\\Delta$irr: ' +  '20%-30%', var_show + ' 99.5%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_20_30,15,2), savgol_filter((data_IRR_25)/area_20_30,15,2), savgol_filter((data_IRR_75)/area_20_30,15,2), r'$\\Delta$irr: ' +  '20%-30%', var_show + ' 99.5%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim2])\n",
    "    percentile = 0.999\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(3,4,10, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'j'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_20_30,15,2), savgol_filter((data_NOI_25)/area_20_30,15,2), savgol_filter((data_NOI_75)/area_20_30,15,2), r'$\\Delta$irr: ' +  '20%-30%', var_show + ' 99.9%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_20_30,15,2), savgol_filter((data_IRR_25)/area_20_30,15,2), savgol_filter((data_IRR_75)/area_20_30,15,2), r'$\\Delta$irr: ' +  '20%-30%', var_show + ' 99.9%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim3])\n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    percentile = 0.99\n",
    "    irr_low = 30\n",
    "    irr_high = 40\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(343, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'c'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_30_40,15,2), savgol_filter((data_NOI_25)/area_30_40,15,2), savgol_filter((data_NOI_75)/area_30_40,15,2), r'$\\Delta$irr: ' +  '30%-40%', var_show + ' 99.0%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_30_40,15,2), savgol_filter((data_IRR_25)/area_30_40,15,2), savgol_filter((data_IRR_75)/area_30_40,15,2), r'$\\Delta$irr: ' +  '30%-40%', var_show + ' 99.0%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim1])\n",
    "    \n",
    "    percentile = 0.995\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(347, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'g'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_30_40,15,2), savgol_filter((data_NOI_25)/area_30_40,15,2), savgol_filter((data_NOI_75)/area_30_40,15,2), r'$\\Delta$irr: ' +  '30%-40%', var_show + ' 99.5%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_30_40,15,2), savgol_filter((data_IRR_25)/area_30_40,15,2), savgol_filter((data_IRR_75)/area_30_40,15,2), r'$\\Delta$irr: ' +  '30%-40%', var_show + ' 99.5%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim2])\n",
    "    percentile = 0.999\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(3,4,11, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'k'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_30_40,15,2), savgol_filter((data_NOI_25)/area_30_40,15,2), savgol_filter((data_NOI_75)/area_30_40,15,2), r'$\\Delta$irr: ' +  '30%-40%', var_show + ' 99.9%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_30_40,15,2), savgol_filter((data_IRR_25)/area_30_40,15,2), savgol_filter((data_IRR_75)/area_30_40,15,2), r'$\\Delta$irr: ' +  '30%-40%', var_show + ' 99.9%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim3])\n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    percentile = 0.99\n",
    "    irr_low = 40\n",
    "    irr_high = 100\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(344, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'd'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_more_40,15,2), savgol_filter((data_NOI_25)/area_more_40,15,2), savgol_filter((data_NOI_75)/area_more_40,15,2), r'$\\Delta$irr: ' +  '>40%', var_show + ' 99.0%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_more_40,15,2), savgol_filter((data_IRR_25)/area_more_40,15,2), savgol_filter((data_IRR_75)/area_more_40,15,2), r'$\\Delta$irr: ' +  '>40%', var_show + ' 99.0%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim1])\n",
    "    \n",
    "    percentile = 0.995\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(348, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'h'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_more_40,15,2), savgol_filter((data_NOI_25)/area_more_40,15,2), savgol_filter((data_NOI_75)/area_more_40,15,2), r'$\\Delta$irr: ' +  '>40%', var_show + ' 99.5%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_more_40,15,2), savgol_filter((data_IRR_25)/area_more_40,15,2), savgol_filter((data_IRR_75)/area_more_40,15,2), r'$\\Delta$irr: ' +  '>40%', var_show + ' 99.5%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim2])\n",
    "    percentile = 0.999\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(3,4,12, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'l'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_more_40,15,2), savgol_filter((data_NOI_25)/area_more_40,15,2), savgol_filter((data_NOI_75)/area_more_40,15,2), r'$\\Delta$irr: ' +  '>40%', var_show + ' 99.9%', '1901irr', 'orange', '', '', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_more_40,15,2), savgol_filter((data_IRR_25)/area_more_40,15,2), savgol_filter((data_IRR_75)/area_more_40,15,2), r'$\\Delta$irr: ' +  '>40%', var_show + ' 99.9%', 'tranirr', 'dodgerblue', '', '', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim3])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "75e13063",
   "metadata": {},
   "outputs": [],
   "source": [
    "all_total('TSA', '$\\mathregular{T_{2m}}$', 350, 200, 80)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b9f59f74",
   "metadata": {},
   "outputs": [],
   "source": [
    "def all_total_for_paper(ylim1, ylim2, ylim3):\n",
    "    \n",
    "    fig = plt.figure(figsize=(12, 3),dpi=300)\n",
    "    fig.subplots_adjust(hspace=0.4, wspace=0.2, left = 0.05, right = 0.95, top = 0.95, bottom = 0.05)\n",
    "    \n",
    "    metric = 'TSA'\n",
    "    var_show = '$\\mathregular{T_{2m}}$'\n",
    "    var = metric\n",
    "    percentile = 0.999\n",
    "    irr_low = 40\n",
    "    irr_high = 100\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    \n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(131, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'a'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_more_40,15,2), savgol_filter((data_NOI_25)/area_more_40,11,2), savgol_filter((data_NOI_75)/area_more_40,11,2), r'$\\Delta$irr: ' +  '>40%', '$\\mathregular{T_{2m}}$ 99.9%', '1901irr', 'orange', 'hours/yr', 'year', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_more_40,11,2), savgol_filter((data_IRR_25)/area_more_40,11,2), savgol_filter((data_IRR_75)/area_more_40,11,2), r'$\\Delta$irr: ' +  '>40%', '$\\mathregular{T_{2m}}$ 99.9%', 'tranirr', 'dodgerblue', 'hours/yr', 'year', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim1])\n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    percentile = 0.999\n",
    "    metric = 'HUMIDEX'\n",
    "    var = metric\n",
    "    var_show = '$\\mathregular{HU}$'\n",
    "    irr_low = 40\n",
    "    irr_high = 100\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(132, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'b'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_more_40,15,2), savgol_filter((data_NOI_25)/area_more_40,11,2), savgol_filter((data_NOI_75)/area_more_40,11,2), r'$\\Delta$irr: ' +  '>40%', 'HU 99.9%', '1901irr', 'orange', '', 'year', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_more_40,11,2), savgol_filter((data_IRR_25)/area_more_40,11,2), savgol_filter((data_IRR_75)/area_more_40,11,2), r'$\\Delta$irr: ' +  '>40%', 'HU 99.9%', 'tranirr', 'dodgerblue', '', 'year', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim2])\n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    \n",
    "    metric = 'T_w'\n",
    "    var = metric\n",
    "    var_show = '$\\mathregular{T_{w}}$'\n",
    "    percentile = 0.999\n",
    "    irr_low = 40\n",
    "    irr_high = 100\n",
    "    data_IRR_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CESM2_fam = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/after_process/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_CNRM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/CNRM/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_IPSL = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/IPSL-CM6/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_IRR_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "    data_NOI_E3SM = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/E3SM_cdo/data_NOI_\" + \\\n",
    "                         var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "#     data_IRR_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_IRR_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_IRR')\n",
    "#     data_NOI_MIROC = get_data_from_mat(\"/dodrio/scratch/projects/2022_200/project_output/cesm/yi_yao_IRRMIP/MIROC/data_NOI_\" + \\\n",
    "#                          var + '_' + str(percentile) + '_' + str(irr_low) + '_' + str(irr_high) + '.mat', 'data_NOI')\n",
    "    \n",
    "    data_IRR = np.vstack((data_IRR_CESM2_fam,\n",
    "                               data_IRR_CNRM,\n",
    "                               data_IRR_IPSL,\n",
    "                               data_IRR_E3SM))\n",
    "    \n",
    "    data_NOI = np.vstack((data_NOI_CESM2_fam,\n",
    "                               data_NOI_CNRM,\n",
    "                               data_NOI_IPSL,\n",
    "                               data_NOI_E3SM))\n",
    "    data_IRR_75 = np.percentile(data_IRR, 75, axis = 0)\n",
    "    data_IRR_50 = np.percentile(data_IRR, 50, axis = 0)\n",
    "    data_IRR_25 = np.percentile(data_IRR, 25, axis = 0)\n",
    "\n",
    "    data_NOI_75 = np.percentile(data_NOI, 75, axis = 0)\n",
    "    data_NOI_50 = np.percentile(data_NOI, 50, axis = 0)\n",
    "    data_NOI_25 = np.percentile(data_NOI, 25, axis = 0)\n",
    "\n",
    "    \n",
    "    ax1 = plt.subplot(133, frameon=True)\n",
    "    \n",
    "\n",
    "    index = 'c'\n",
    "    plot_line(ax1, index, savgol_filter((data_NOI_50)/area_more_40,15,2), savgol_filter((data_NOI_25)/area_more_40,11,2), savgol_filter((data_NOI_75)/area_more_40,11,2), r'$\\Delta$irr: ' +  '>40%', '$\\mathregular{T_{w}}$ 99.9%', '1901irr', 'orange', '', 'year', '')\n",
    "    plot_line(ax1, index, savgol_filter((data_IRR_50)/area_more_40,11,2), savgol_filter((data_IRR_25)/area_more_40,11,2), savgol_filter((data_IRR_75)/area_more_40,11,2), r'$\\Delta$irr: ' +  '>40%', '$\\mathregular{T_{w}}$ 99.9%', 'tranirr', 'dodgerblue', '', 'year', '')\n",
    "    plt.grid(linestyle = '--')\n",
    "    plt.ylim([0, ylim3])\n",
    "    \n",
    "    "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9b9f2a11",
   "metadata": {},
   "outputs": [],
   "source": [
    "all_total_for_paper(30, 120, 240)"
   ]
  }
 ],
 "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.9.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}