cre_pex_regional_extremes.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "import xarray as xr\n",
    "import numpy as np\n",
    "from pathlib import Path\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "from cycler import cycler\n",
    "import palettable\n",
    "import colorcet as cc"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "def _load(p):\n",
    "    if not isinstance(p, list):\n",
    "        p = [p]\n",
    "    if len(p) == 1:\n",
    "        ds = xr.open_dataset(p[0])\n",
    "    else:\n",
    "        ds = xr.open_mfdataset(p)\n",
    "    dv = ds.data_vars\n",
    "    if 'pr' in dv:\n",
    "        return dv['pr']\n",
    "    elif 'PRECT' in dv:\n",
    "        return dv['PRECT']\n",
    "    else:\n",
    "        raise IOError(\"Could not find pr or PRECT\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "#\n",
    "# We have 3 pieces of information, each in its own file: zonal mean, quantiles, histogram\n",
    "#\n",
    "region_string = \"tropical\" # \"nh_extratropics\"  #  \n",
    "p0 = Path(\"/Volumes/Jedha/\")\n",
    "p1 = p0 / \"CMIP6\"\n",
    "p2 = p0 / \"cloud_locking\"\n",
    "\n",
    "models = {\"CESM2\": ('aqua-control', 'aqua-control-lwoff', 'amip', 'amip-lwoff', 'f-cntl', 'f-lock', 'c-cntl', 'c-lock'),\n",
    "         'CNRM-CM6-1': ('aqua-control', 'aqua-control-lwoff', 'amip', 'amip-lwoff'),\n",
    "         'IPSL-CM6A-LR': ('aqua-control', 'aqua-control-lwoff', 'amip', 'amip-lwoff'),\n",
    "         'HadGEM3-GC31-LL': ('aqua-control', 'aqua-control-lwoff', 'amip', 'amip-lwoff'),\n",
    "         'MRI-ESM2-0': ('amip', 'amip-lwoff'),\n",
    "         }\n",
    "\n",
    "\n",
    "\n",
    "use_member = {}\n",
    "for m in models:\n",
    "    if m == 'CNRM-CM6-1':\n",
    "        use_member[m] = 'r1i1p1f2'\n",
    "    elif m == 'HadGEM3-GC31-LL':\n",
    "        use_member[m] = '*'\n",
    "    else:\n",
    "        use_member[m] = 'r1i1p1f1'\n",
    "\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We want to construct a stacked bar plot that shows the precipitation distribution broken into regions. I don't think we need to worry about area weighting for this."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 104,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "pr_day_MRI-ESM2-0_amip_r1i1p1f1_*\n",
      "Number of files: 1\n",
      "[PosixPath('/Volumes/Jedha/CMIP6/pr_day_MRI-ESM2-0_amip_r1i1p1f1_gn_19790101-20141231.nc')]\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;pr&#x27; (time: 13149, lat: 160, lon: 320)&gt;\n",
       "array([[[2.67101198e-01, 2.35772356e-01, 2.30664775e-01, ...,\n",
       "         2.52189100e-01, 2.32562289e-01, 2.26287246e-01],\n",
       "        [1.31279632e-01, 1.35269061e-01, 1.30600467e-01, ...,\n",
       "         1.40292794e-01, 1.37421206e-01, 1.33562788e-01],\n",
       "        [4.71350811e-02, 4.46746834e-02, 4.45041358e-02, ...,\n",
       "         5.06040566e-02, 4.89875786e-02, 4.85395156e-02],\n",
       "        ...,\n",
       "        [7.56014436e-02, 8.67677927e-02, 1.02531180e-01, ...,\n",
       "         1.13704197e-01, 9.00612846e-02, 7.05437735e-02],\n",
       "        [7.91385863e-03, 7.42563559e-03, 9.16410051e-03, ...,\n",
       "         1.25986254e-02, 1.44611746e-02, 9.40386206e-03],\n",
       "        [1.35917086e-02, 1.43885324e-02, 1.55960834e-02, ...,\n",
       "         1.23632764e-02, 1.23326853e-02, 1.27605749e-02]],\n",
       "\n",
       "       [[1.33665621e-01, 1.42159492e-01, 1.20417401e-01, ...,\n",
       "         1.33317068e-01, 1.34539649e-01, 1.33625612e-01],\n",
       "        [1.13974206e-01, 1.15015261e-01, 1.16569333e-01, ...,\n",
       "         1.18133925e-01, 1.13970421e-01, 1.23785287e-01],\n",
       "        [1.73828006e-01, 1.60179630e-01, 1.93370670e-01, ...,\n",
       "         1.75499454e-01, 1.66666463e-01, 1.70476928e-01],\n",
       "...\n",
       "        [2.32361421e-01, 2.29877830e-01, 2.33837441e-01, ...,\n",
       "         2.85293251e-01, 2.46625215e-01, 2.36508369e-01],\n",
       "        [4.62087721e-01, 4.59472626e-01, 4.53511268e-01, ...,\n",
       "         4.82449055e-01, 4.81340319e-01, 4.69268948e-01],\n",
       "        [9.66097474e-01, 1.16981566e+00, 1.08266509e+00, ...,\n",
       "         9.62051332e-01, 9.97630119e-01, 9.77471828e-01]],\n",
       "\n",
       "       [[2.04111841e-02, 2.02876385e-02, 2.02328227e-02, ...,\n",
       "         2.09000967e-02, 2.08847765e-02, 2.07346268e-02],\n",
       "        [6.02674820e-02, 6.05076551e-02, 6.12820685e-02, ...,\n",
       "         6.53327033e-02, 6.66729286e-02, 6.27168566e-02],\n",
       "        [1.30121157e-01, 1.39833197e-01, 1.40520871e-01, ...,\n",
       "         1.18493840e-01, 1.30778000e-01, 1.31784067e-01],\n",
       "        ...,\n",
       "        [8.84294659e-02, 1.06997639e-01, 6.63979277e-02, ...,\n",
       "         1.01991758e-01, 1.31017476e-01, 1.00867175e-01],\n",
       "        [1.24816984e-01, 1.19286813e-01, 1.16639100e-01, ...,\n",
       "         1.24225348e-01, 1.27329558e-01, 1.35236725e-01],\n",
       "        [5.14209978e-02, 5.40016480e-02, 5.38464859e-02, ...,\n",
       "         5.76767549e-02, 5.57630956e-02, 5.42576984e-02]]], dtype=float32)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1979-01-01T12:00:00 ... 2014-12-31T12:00:00\n",
       "  * lat      (lat) float64 -89.14 -88.03 -86.91 -85.79 ... 86.91 88.03 89.14\n",
       "  * lon      (lon) float64 0.0 1.125 2.25 3.375 4.5 ... 355.5 356.6 357.8 358.9\n",
       "Attributes:\n",
       "    standard_name:  precipitation_flux\n",
       "    long_name:      Precipitation\n",
       "    comment:        includes both liquid and solid phases\n",
       "    units:          mm d-1\n",
       "    original_name:  PRECIPI\n",
       "    cell_methods:   area: time: mean\n",
       "    cell_measures:  area: areacella\n",
       "    history:        2019-05-22T02:27:29Z altered by CMOR: replaced missing va...</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'pr'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 13149</li><li><span class='xr-has-index'>lat</span>: 160</li><li><span class='xr-has-index'>lon</span>: 320</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-a930523c-5a01-48a2-9a1e-bea60dfe28ef' class='xr-array-in' type='checkbox' checked><label for='section-a930523c-5a01-48a2-9a1e-bea60dfe28ef' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>0.2671 0.2358 0.2307 0.2427 0.2456 ... 0.06071 0.05768 0.05576 0.05426</span></div><div class='xr-array-data'><pre>array([[[2.67101198e-01, 2.35772356e-01, 2.30664775e-01, ...,\n",
       "         2.52189100e-01, 2.32562289e-01, 2.26287246e-01],\n",
       "        [1.31279632e-01, 1.35269061e-01, 1.30600467e-01, ...,\n",
       "         1.40292794e-01, 1.37421206e-01, 1.33562788e-01],\n",
       "        [4.71350811e-02, 4.46746834e-02, 4.45041358e-02, ...,\n",
       "         5.06040566e-02, 4.89875786e-02, 4.85395156e-02],\n",
       "        ...,\n",
       "        [7.56014436e-02, 8.67677927e-02, 1.02531180e-01, ...,\n",
       "         1.13704197e-01, 9.00612846e-02, 7.05437735e-02],\n",
       "        [7.91385863e-03, 7.42563559e-03, 9.16410051e-03, ...,\n",
       "         1.25986254e-02, 1.44611746e-02, 9.40386206e-03],\n",
       "        [1.35917086e-02, 1.43885324e-02, 1.55960834e-02, ...,\n",
       "         1.23632764e-02, 1.23326853e-02, 1.27605749e-02]],\n",
       "\n",
       "       [[1.33665621e-01, 1.42159492e-01, 1.20417401e-01, ...,\n",
       "         1.33317068e-01, 1.34539649e-01, 1.33625612e-01],\n",
       "        [1.13974206e-01, 1.15015261e-01, 1.16569333e-01, ...,\n",
       "         1.18133925e-01, 1.13970421e-01, 1.23785287e-01],\n",
       "        [1.73828006e-01, 1.60179630e-01, 1.93370670e-01, ...,\n",
       "         1.75499454e-01, 1.66666463e-01, 1.70476928e-01],\n",
       "...\n",
       "        [2.32361421e-01, 2.29877830e-01, 2.33837441e-01, ...,\n",
       "         2.85293251e-01, 2.46625215e-01, 2.36508369e-01],\n",
       "        [4.62087721e-01, 4.59472626e-01, 4.53511268e-01, ...,\n",
       "         4.82449055e-01, 4.81340319e-01, 4.69268948e-01],\n",
       "        [9.66097474e-01, 1.16981566e+00, 1.08266509e+00, ...,\n",
       "         9.62051332e-01, 9.97630119e-01, 9.77471828e-01]],\n",
       "\n",
       "       [[2.04111841e-02, 2.02876385e-02, 2.02328227e-02, ...,\n",
       "         2.09000967e-02, 2.08847765e-02, 2.07346268e-02],\n",
       "        [6.02674820e-02, 6.05076551e-02, 6.12820685e-02, ...,\n",
       "         6.53327033e-02, 6.66729286e-02, 6.27168566e-02],\n",
       "        [1.30121157e-01, 1.39833197e-01, 1.40520871e-01, ...,\n",
       "         1.18493840e-01, 1.30778000e-01, 1.31784067e-01],\n",
       "        ...,\n",
       "        [8.84294659e-02, 1.06997639e-01, 6.63979277e-02, ...,\n",
       "         1.01991758e-01, 1.31017476e-01, 1.00867175e-01],\n",
       "        [1.24816984e-01, 1.19286813e-01, 1.16639100e-01, ...,\n",
       "         1.24225348e-01, 1.27329558e-01, 1.35236725e-01],\n",
       "        [5.14209978e-02, 5.40016480e-02, 5.38464859e-02, ...,\n",
       "         5.76767549e-02, 5.57630956e-02, 5.42576984e-02]]], dtype=float32)</pre></div></div></li><li class='xr-section-item'><input id='section-5fb44d1b-8611-4f02-92d9-08e05641b0b4' class='xr-section-summary-in' type='checkbox'  checked><label for='section-5fb44d1b-8611-4f02-92d9-08e05641b0b4' class='xr-section-summary' >Coordinates: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>time</span></div><div class='xr-var-dims'>(time)</div><div class='xr-var-dtype'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>1979-01-01T12:00:00 ... 2014-12-...</div><input id='attrs-6184bdcb-a5bc-4c8c-9071-d5c8a0b23a03' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6184bdcb-a5bc-4c8c-9071-d5c8a0b23a03' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-63dcffe9-c44b-4511-9db5-96f411e65869' class='xr-var-data-in' type='checkbox'><label for='data-63dcffe9-c44b-4511-9db5-96f411e65869' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>time_bnds</dd><dt><span>axis :</span></dt><dd>T</dd><dt><span>long_name :</span></dt><dd>time</dd><dt><span>standard_name :</span></dt><dd>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;1979-01-01T12:00:00.000000000&#x27;, &#x27;1979-01-02T12:00:00.000000000&#x27;,\n",
       "       &#x27;1979-01-03T12:00:00.000000000&#x27;, ..., &#x27;2014-12-29T12:00:00.000000000&#x27;,\n",
       "       &#x27;2014-12-30T12:00:00.000000000&#x27;, &#x27;2014-12-31T12:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-89.14 -88.03 ... 88.03 89.14</div><input id='attrs-32f0c9eb-20fd-4c67-8ef8-90e65bb26cb6' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-32f0c9eb-20fd-4c67-8ef8-90e65bb26cb6' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f81698ff-834b-4bd2-a957-a6ae549ef7d0' class='xr-var-data-in' type='checkbox'><label for='data-f81698ff-834b-4bd2-a957-a6ae549ef7d0' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>lat_bnds</dd><dt><span>units :</span></dt><dd>degrees_north</dd><dt><span>axis :</span></dt><dd>Y</dd><dt><span>long_name :</span></dt><dd>Latitude</dd><dt><span>standard_name :</span></dt><dd>latitude</dd></dl></div><div class='xr-var-data'><pre>array([-89.14152, -88.02943, -86.91077, -85.79063, -84.66992, -83.54895,\n",
       "       -82.42782, -81.30659, -80.18531, -79.06398, -77.94262, -76.82124,\n",
       "       -75.69984, -74.57843, -73.45701, -72.33558, -71.21414, -70.09269,\n",
       "       -68.97124, -67.84978, -66.72833, -65.60686, -64.4854 , -63.36393,\n",
       "       -62.24246, -61.12099, -59.99952, -58.87804, -57.75657, -56.63509,\n",
       "       -55.51361, -54.39214, -53.27066, -52.14917, -51.02769, -49.90621,\n",
       "       -48.78473, -47.66325, -46.54176, -45.42028, -44.29879, -43.17731,\n",
       "       -42.05582, -40.93434, -39.81285, -38.69137, -37.56988, -36.44839,\n",
       "       -35.32691, -34.20542, -33.08393, -31.96244, -30.84096, -29.71947,\n",
       "       -28.59798, -27.47649, -26.355  , -25.23351, -24.11203, -22.99054,\n",
       "       -21.86905, -20.74756, -19.62607, -18.50458, -17.38309, -16.2616 ,\n",
       "       -15.14011, -14.01862, -12.89713, -11.77564, -10.65415,  -9.53266,\n",
       "        -8.41117,  -7.28968,  -6.16819,  -5.0467 ,  -3.92521,  -2.80372,\n",
       "        -1.68223,  -0.56074,   0.56074,   1.68223,   2.80372,   3.92521,\n",
       "         5.0467 ,   6.16819,   7.28968,   8.41117,   9.53266,  10.65415,\n",
       "        11.77564,  12.89713,  14.01862,  15.14011,  16.2616 ,  17.38309,\n",
       "        18.50458,  19.62607,  20.74756,  21.86905,  22.99054,  24.11203,\n",
       "        25.23351,  26.355  ,  27.47649,  28.59798,  29.71947,  30.84096,\n",
       "        31.96244,  33.08393,  34.20542,  35.32691,  36.44839,  37.56988,\n",
       "        38.69137,  39.81285,  40.93434,  42.05582,  43.17731,  44.29879,\n",
       "        45.42028,  46.54176,  47.66325,  48.78473,  49.90621,  51.02769,\n",
       "        52.14917,  53.27066,  54.39214,  55.51361,  56.63509,  57.75657,\n",
       "        58.87804,  59.99952,  61.12099,  62.24246,  63.36393,  64.4854 ,\n",
       "        65.60686,  66.72833,  67.84978,  68.97124,  70.09269,  71.21414,\n",
       "        72.33558,  73.45701,  74.57843,  75.69984,  76.82124,  77.94262,\n",
       "        79.06398,  80.18531,  81.30659,  82.42782,  83.54895,  84.66992,\n",
       "        85.79063,  86.91077,  88.02943,  89.14152])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0 1.125 2.25 ... 357.8 358.9</div><input id='attrs-dff3538c-7ad3-4381-ba8b-bd1fdd237612' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-dff3538c-7ad3-4381-ba8b-bd1fdd237612' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-14f7316b-aac6-4fc3-a693-a61227d6a33d' class='xr-var-data-in' type='checkbox'><label for='data-14f7316b-aac6-4fc3-a693-a61227d6a33d' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>lon_bnds</dd><dt><span>units :</span></dt><dd>degrees_east</dd><dt><span>axis :</span></dt><dd>X</dd><dt><span>long_name :</span></dt><dd>Longitude</dd><dt><span>standard_name :</span></dt><dd>longitude</dd></dl></div><div class='xr-var-data'><pre>array([  0.   ,   1.125,   2.25 , ..., 356.625, 357.75 , 358.875])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-76b768f3-60a5-4b4a-ad6e-3d6e03d4ca96' class='xr-section-summary-in' type='checkbox'  checked><label for='section-76b768f3-60a5-4b4a-ad6e-3d6e03d4ca96' class='xr-section-summary' >Attributes: <span>(8)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>standard_name :</span></dt><dd>precipitation_flux</dd><dt><span>long_name :</span></dt><dd>Precipitation</dd><dt><span>comment :</span></dt><dd>includes both liquid and solid phases</dd><dt><span>units :</span></dt><dd>mm d-1</dd><dt><span>original_name :</span></dt><dd>PRECIPI</dd><dt><span>cell_methods :</span></dt><dd>area: time: mean</dd><dt><span>cell_measures :</span></dt><dd>area: areacella</dd><dt><span>history :</span></dt><dd>2019-05-22T02:27:29Z altered by CMOR: replaced missing value flag (-9.99e+33) with standard missing value (1e+20).</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'pr' (time: 13149, lat: 160, lon: 320)>\n",
       "array([[[2.67101198e-01, 2.35772356e-01, 2.30664775e-01, ...,\n",
       "         2.52189100e-01, 2.32562289e-01, 2.26287246e-01],\n",
       "        [1.31279632e-01, 1.35269061e-01, 1.30600467e-01, ...,\n",
       "         1.40292794e-01, 1.37421206e-01, 1.33562788e-01],\n",
       "        [4.71350811e-02, 4.46746834e-02, 4.45041358e-02, ...,\n",
       "         5.06040566e-02, 4.89875786e-02, 4.85395156e-02],\n",
       "        ...,\n",
       "        [7.56014436e-02, 8.67677927e-02, 1.02531180e-01, ...,\n",
       "         1.13704197e-01, 9.00612846e-02, 7.05437735e-02],\n",
       "        [7.91385863e-03, 7.42563559e-03, 9.16410051e-03, ...,\n",
       "         1.25986254e-02, 1.44611746e-02, 9.40386206e-03],\n",
       "        [1.35917086e-02, 1.43885324e-02, 1.55960834e-02, ...,\n",
       "         1.23632764e-02, 1.23326853e-02, 1.27605749e-02]],\n",
       "\n",
       "       [[1.33665621e-01, 1.42159492e-01, 1.20417401e-01, ...,\n",
       "         1.33317068e-01, 1.34539649e-01, 1.33625612e-01],\n",
       "        [1.13974206e-01, 1.15015261e-01, 1.16569333e-01, ...,\n",
       "         1.18133925e-01, 1.13970421e-01, 1.23785287e-01],\n",
       "        [1.73828006e-01, 1.60179630e-01, 1.93370670e-01, ...,\n",
       "         1.75499454e-01, 1.66666463e-01, 1.70476928e-01],\n",
       "...\n",
       "        [2.32361421e-01, 2.29877830e-01, 2.33837441e-01, ...,\n",
       "         2.85293251e-01, 2.46625215e-01, 2.36508369e-01],\n",
       "        [4.62087721e-01, 4.59472626e-01, 4.53511268e-01, ...,\n",
       "         4.82449055e-01, 4.81340319e-01, 4.69268948e-01],\n",
       "        [9.66097474e-01, 1.16981566e+00, 1.08266509e+00, ...,\n",
       "         9.62051332e-01, 9.97630119e-01, 9.77471828e-01]],\n",
       "\n",
       "       [[2.04111841e-02, 2.02876385e-02, 2.02328227e-02, ...,\n",
       "         2.09000967e-02, 2.08847765e-02, 2.07346268e-02],\n",
       "        [6.02674820e-02, 6.05076551e-02, 6.12820685e-02, ...,\n",
       "         6.53327033e-02, 6.66729286e-02, 6.27168566e-02],\n",
       "        [1.30121157e-01, 1.39833197e-01, 1.40520871e-01, ...,\n",
       "         1.18493840e-01, 1.30778000e-01, 1.31784067e-01],\n",
       "        ...,\n",
       "        [8.84294659e-02, 1.06997639e-01, 6.63979277e-02, ...,\n",
       "         1.01991758e-01, 1.31017476e-01, 1.00867175e-01],\n",
       "        [1.24816984e-01, 1.19286813e-01, 1.16639100e-01, ...,\n",
       "         1.24225348e-01, 1.27329558e-01, 1.35236725e-01],\n",
       "        [5.14209978e-02, 5.40016480e-02, 5.38464859e-02, ...,\n",
       "         5.76767549e-02, 5.57630956e-02, 5.42576984e-02]]], dtype=float32)\n",
       "Coordinates:\n",
       "  * time     (time) datetime64[ns] 1979-01-01T12:00:00 ... 2014-12-31T12:00:00\n",
       "  * lat      (lat) float64 -89.14 -88.03 -86.91 -85.79 ... 86.91 88.03 89.14\n",
       "  * lon      (lon) float64 0.0 1.125 2.25 3.375 4.5 ... 355.5 356.6 357.8 358.9\n",
       "Attributes:\n",
       "    standard_name:  precipitation_flux\n",
       "    long_name:      Precipitation\n",
       "    comment:        includes both liquid and solid phases\n",
       "    units:          mm d-1\n",
       "    original_name:  PRECIPI\n",
       "    cell_methods:   area: time: mean\n",
       "    cell_measures:  area: areacella\n",
       "    history:        2019-05-22T02:27:29Z altered by CMOR: replaced missing va..."
      ]
     },
     "execution_count": 104,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "model_name = 'MRI-ESM2-0'\n",
    "experiment_name = 'amip'\n",
    "member = use_member[model_name]\n",
    "print(f'pr_day_{model_name}_{experiment_name}_{member}_*')\n",
    "fils = sorted(list(Path(\"/Volumes/Jedha/CMIP6\").glob(f'pr_day_{model_name}_{experiment_name}_{member}_*')))\n",
    "print(f\"Number of files: {len(fils)}\")\n",
    "print(fils)\n",
    "if len(fils) > 1:\n",
    "    ds = xr.open_mfdataset(fils)\n",
    "else:\n",
    "    ds = xr.open_dataset(fils[0])\n",
    "pr = ds['pr']\n",
    "# convert kg m-2 s-1 to mm/day \n",
    "# * 1/1000kg/m3 * 1000 mm/m * 86400 s/d \n",
    "pr *= 86400.0\n",
    "pr.attrs['units'] = 'mm d-1'\n",
    "pr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 105,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div><svg style=\"position: absolute; width: 0; height: 0; overflow: hidden\">\n",
       "<defs>\n",
       "<symbol id=\"icon-database\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M16 0c-8.837 0-16 2.239-16 5v4c0 2.761 7.163 5 16 5s16-2.239 16-5v-4c0-2.761-7.163-5-16-5z\"></path>\n",
       "<path d=\"M16 17c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "<path d=\"M16 26c-8.837 0-16-2.239-16-5v6c0 2.761 7.163 5 16 5s16-2.239 16-5v-6c0 2.761-7.163 5-16 5z\"></path>\n",
       "</symbol>\n",
       "<symbol id=\"icon-file-text2\" viewBox=\"0 0 32 32\">\n",
       "<path d=\"M28.681 7.159c-0.694-0.947-1.662-2.053-2.724-3.116s-2.169-2.030-3.116-2.724c-1.612-1.182-2.393-1.319-2.841-1.319h-15.5c-1.378 0-2.5 1.121-2.5 2.5v27c0 1.378 1.122 2.5 2.5 2.5h23c1.378 0 2.5-1.122 2.5-2.5v-19.5c0-0.448-0.137-1.23-1.319-2.841zM24.543 5.457c0.959 0.959 1.712 1.825 2.268 2.543h-4.811v-4.811c0.718 0.556 1.584 1.309 2.543 2.268zM28 29.5c0 0.271-0.229 0.5-0.5 0.5h-23c-0.271 0-0.5-0.229-0.5-0.5v-27c0-0.271 0.229-0.5 0.5-0.5 0 0 15.499-0 15.5 0v7c0 0.552 0.448 1 1 1h7v19.5z\"></path>\n",
       "<path d=\"M23 26h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 22h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "<path d=\"M23 18h-14c-0.552 0-1-0.448-1-1s0.448-1 1-1h14c0.552 0 1 0.448 1 1s-0.448 1-1 1z\"></path>\n",
       "</symbol>\n",
       "</defs>\n",
       "</svg>\n",
       "<style>/* CSS stylesheet for displaying xarray objects in jupyterlab.\n",
       " *\n",
       " */\n",
       "\n",
       ":root {\n",
       "  --xr-font-color0: var(--jp-content-font-color0, rgba(0, 0, 0, 1));\n",
       "  --xr-font-color2: var(--jp-content-font-color2, rgba(0, 0, 0, 0.54));\n",
       "  --xr-font-color3: var(--jp-content-font-color3, rgba(0, 0, 0, 0.38));\n",
       "  --xr-border-color: var(--jp-border-color2, #e0e0e0);\n",
       "  --xr-disabled-color: var(--jp-layout-color3, #bdbdbd);\n",
       "  --xr-background-color: var(--jp-layout-color0, white);\n",
       "  --xr-background-color-row-even: var(--jp-layout-color1, white);\n",
       "  --xr-background-color-row-odd: var(--jp-layout-color2, #eeeeee);\n",
       "}\n",
       "\n",
       "html[theme=dark],\n",
       "body.vscode-dark {\n",
       "  --xr-font-color0: rgba(255, 255, 255, 1);\n",
       "  --xr-font-color2: rgba(255, 255, 255, 0.54);\n",
       "  --xr-font-color3: rgba(255, 255, 255, 0.38);\n",
       "  --xr-border-color: #1F1F1F;\n",
       "  --xr-disabled-color: #515151;\n",
       "  --xr-background-color: #111111;\n",
       "  --xr-background-color-row-even: #111111;\n",
       "  --xr-background-color-row-odd: #313131;\n",
       "}\n",
       "\n",
       ".xr-wrap {\n",
       "  display: block;\n",
       "  min-width: 300px;\n",
       "  max-width: 700px;\n",
       "}\n",
       "\n",
       ".xr-text-repr-fallback {\n",
       "  /* fallback to plain text repr when CSS is not injected (untrusted notebook) */\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-header {\n",
       "  padding-top: 6px;\n",
       "  padding-bottom: 6px;\n",
       "  margin-bottom: 4px;\n",
       "  border-bottom: solid 1px var(--xr-border-color);\n",
       "}\n",
       "\n",
       ".xr-header > div,\n",
       ".xr-header > ul {\n",
       "  display: inline;\n",
       "  margin-top: 0;\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-obj-type,\n",
       ".xr-array-name {\n",
       "  margin-left: 2px;\n",
       "  margin-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-obj-type {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-sections {\n",
       "  padding-left: 0 !important;\n",
       "  display: grid;\n",
       "  grid-template-columns: 150px auto auto 1fr 20px 20px;\n",
       "}\n",
       "\n",
       ".xr-section-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-section-item input {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-item input + label {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label {\n",
       "  cursor: pointer;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-item input:enabled + label:hover {\n",
       "  color: var(--xr-font-color0);\n",
       "}\n",
       "\n",
       ".xr-section-summary {\n",
       "  grid-column: 1;\n",
       "  color: var(--xr-font-color2);\n",
       "  font-weight: 500;\n",
       "}\n",
       "\n",
       ".xr-section-summary > span {\n",
       "  display: inline-block;\n",
       "  padding-left: 0.5em;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label {\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in + label:before {\n",
       "  display: inline-block;\n",
       "  content: '►';\n",
       "  font-size: 11px;\n",
       "  width: 15px;\n",
       "  text-align: center;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:disabled + label:before {\n",
       "  color: var(--xr-disabled-color);\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label:before {\n",
       "  content: '▼';\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked + label > span {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-section-summary,\n",
       ".xr-section-inline-details {\n",
       "  padding-top: 4px;\n",
       "  padding-bottom: 4px;\n",
       "}\n",
       "\n",
       ".xr-section-inline-details {\n",
       "  grid-column: 2 / -1;\n",
       "}\n",
       "\n",
       ".xr-section-details {\n",
       "  display: none;\n",
       "  grid-column: 1 / -1;\n",
       "  margin-bottom: 5px;\n",
       "}\n",
       "\n",
       ".xr-section-summary-in:checked ~ .xr-section-details {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-array-wrap {\n",
       "  grid-column: 1 / -1;\n",
       "  display: grid;\n",
       "  grid-template-columns: 20px auto;\n",
       "}\n",
       "\n",
       ".xr-array-wrap > label {\n",
       "  grid-column: 1;\n",
       "  vertical-align: top;\n",
       "}\n",
       "\n",
       ".xr-preview {\n",
       "  color: var(--xr-font-color3);\n",
       "}\n",
       "\n",
       ".xr-array-preview,\n",
       ".xr-array-data {\n",
       "  padding: 0 5px !important;\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-array-data,\n",
       ".xr-array-in:checked ~ .xr-array-preview {\n",
       "  display: none;\n",
       "}\n",
       "\n",
       ".xr-array-in:checked ~ .xr-array-data,\n",
       ".xr-array-preview {\n",
       "  display: inline-block;\n",
       "}\n",
       "\n",
       ".xr-dim-list {\n",
       "  display: inline-block !important;\n",
       "  list-style: none;\n",
       "  padding: 0 !important;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list li {\n",
       "  display: inline-block;\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "}\n",
       "\n",
       ".xr-dim-list:before {\n",
       "  content: '(';\n",
       "}\n",
       "\n",
       ".xr-dim-list:after {\n",
       "  content: ')';\n",
       "}\n",
       "\n",
       ".xr-dim-list li:not(:last-child):after {\n",
       "  content: ',';\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-has-index {\n",
       "  font-weight: bold;\n",
       "}\n",
       "\n",
       ".xr-var-list,\n",
       ".xr-var-item {\n",
       "  display: contents;\n",
       "}\n",
       "\n",
       ".xr-var-item > div,\n",
       ".xr-var-item label,\n",
       ".xr-var-item > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-even);\n",
       "  margin-bottom: 0;\n",
       "}\n",
       "\n",
       ".xr-var-item > .xr-var-name:hover span {\n",
       "  padding-right: 5px;\n",
       "}\n",
       "\n",
       ".xr-var-list > li:nth-child(odd) > div,\n",
       ".xr-var-list > li:nth-child(odd) > label,\n",
       ".xr-var-list > li:nth-child(odd) > .xr-var-name span {\n",
       "  background-color: var(--xr-background-color-row-odd);\n",
       "}\n",
       "\n",
       ".xr-var-name {\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-var-dims {\n",
       "  grid-column: 2;\n",
       "}\n",
       "\n",
       ".xr-var-dtype {\n",
       "  grid-column: 3;\n",
       "  text-align: right;\n",
       "  color: var(--xr-font-color2);\n",
       "}\n",
       "\n",
       ".xr-var-preview {\n",
       "  grid-column: 4;\n",
       "}\n",
       "\n",
       ".xr-var-name,\n",
       ".xr-var-dims,\n",
       ".xr-var-dtype,\n",
       ".xr-preview,\n",
       ".xr-attrs dt {\n",
       "  white-space: nowrap;\n",
       "  overflow: hidden;\n",
       "  text-overflow: ellipsis;\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-var-name:hover,\n",
       ".xr-var-dims:hover,\n",
       ".xr-var-dtype:hover,\n",
       ".xr-attrs dt:hover {\n",
       "  overflow: visible;\n",
       "  width: auto;\n",
       "  z-index: 1;\n",
       "}\n",
       "\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  display: none;\n",
       "  background-color: var(--xr-background-color) !important;\n",
       "  padding-bottom: 5px !important;\n",
       "}\n",
       "\n",
       ".xr-var-attrs-in:checked ~ .xr-var-attrs,\n",
       ".xr-var-data-in:checked ~ .xr-var-data {\n",
       "  display: block;\n",
       "}\n",
       "\n",
       ".xr-var-data > table {\n",
       "  float: right;\n",
       "}\n",
       "\n",
       ".xr-var-name span,\n",
       ".xr-var-data,\n",
       ".xr-attrs {\n",
       "  padding-left: 25px !important;\n",
       "}\n",
       "\n",
       ".xr-attrs,\n",
       ".xr-var-attrs,\n",
       ".xr-var-data {\n",
       "  grid-column: 1 / -1;\n",
       "}\n",
       "\n",
       "dl.xr-attrs {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  display: grid;\n",
       "  grid-template-columns: 125px auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt,\n",
       ".xr-attrs dd {\n",
       "  padding: 0;\n",
       "  margin: 0;\n",
       "  float: left;\n",
       "  padding-right: 10px;\n",
       "  width: auto;\n",
       "}\n",
       "\n",
       ".xr-attrs dt {\n",
       "  font-weight: normal;\n",
       "  grid-column: 1;\n",
       "}\n",
       "\n",
       ".xr-attrs dt:hover span {\n",
       "  display: inline-block;\n",
       "  background: var(--xr-background-color);\n",
       "  padding-right: 10px;\n",
       "}\n",
       "\n",
       ".xr-attrs dd {\n",
       "  grid-column: 2;\n",
       "  white-space: pre-wrap;\n",
       "  word-break: break-all;\n",
       "}\n",
       "\n",
       ".xr-icon-database,\n",
       ".xr-icon-file-text2 {\n",
       "  display: inline-block;\n",
       "  vertical-align: middle;\n",
       "  width: 1em;\n",
       "  height: 1.5em !important;\n",
       "  stroke-width: 0;\n",
       "  stroke: currentColor;\n",
       "  fill: currentColor;\n",
       "}\n",
       "</style><pre class='xr-text-repr-fallback'>&lt;xarray.DataArray &#x27;sftlf&#x27; (lat: 160, lon: 320)&gt;\n",
       "array([[100., 100., 100., ..., 100., 100., 100.],\n",
       "       [100., 100., 100., ..., 100., 100., 100.],\n",
       "       [100., 100., 100., ..., 100., 100., 100.],\n",
       "       ...,\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.],\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.],\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.]], dtype=float32)\n",
       "Coordinates:\n",
       "  * lat      (lat) float64 -89.14 -88.03 -86.91 -85.79 ... 86.91 88.03 89.14\n",
       "  * lon      (lon) float64 0.0 1.125 2.25 3.375 4.5 ... 355.5 356.6 357.8 358.9\n",
       "    type     |S4 b&#x27;land&#x27;\n",
       "Attributes:\n",
       "    standard_name:   land_area_fraction\n",
       "    long_name:       Percentage of the grid  cell occupied by land (including...\n",
       "    comment:         Please express &#x27;X_area_fraction&#x27; as the percentage of ho...\n",
       "    units:           %\n",
       "    original_name:   lfrac\n",
       "    original_units:  1\n",
       "    history:         2019-05-24T05:34:57Z altered by CMOR: Converted units fr...\n",
       "    cell_methods:    area: mean\n",
       "    cell_measures:   area: areacella</pre><div class='xr-wrap' hidden><div class='xr-header'><div class='xr-obj-type'>xarray.DataArray</div><div class='xr-array-name'>'sftlf'</div><ul class='xr-dim-list'><li><span class='xr-has-index'>lat</span>: 160</li><li><span class='xr-has-index'>lon</span>: 320</li></ul></div><ul class='xr-sections'><li class='xr-section-item'><div class='xr-array-wrap'><input id='section-7d7d2cef-2568-4f67-b21b-4e546e6ac06d' class='xr-array-in' type='checkbox' checked><label for='section-7d7d2cef-2568-4f67-b21b-4e546e6ac06d' title='Show/hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-array-preview xr-preview'><span>100.0 100.0 100.0 100.0 100.0 100.0 100.0 ... 0.0 0.0 0.0 0.0 0.0 0.0</span></div><div class='xr-array-data'><pre>array([[100., 100., 100., ..., 100., 100., 100.],\n",
       "       [100., 100., 100., ..., 100., 100., 100.],\n",
       "       [100., 100., 100., ..., 100., 100., 100.],\n",
       "       ...,\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.],\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.],\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.]], dtype=float32)</pre></div></div></li><li class='xr-section-item'><input id='section-21464fd2-c1d4-4a86-86a3-1df4838e635c' class='xr-section-summary-in' type='checkbox'  checked><label for='section-21464fd2-c1d4-4a86-86a3-1df4838e635c' class='xr-section-summary' >Coordinates: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-89.14 -88.03 ... 88.03 89.14</div><input id='attrs-7993aa72-6b26-42ef-aec4-da5f5accd8d4' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7993aa72-6b26-42ef-aec4-da5f5accd8d4' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-ce80b180-3989-4e59-892f-4bc60f36772f' class='xr-var-data-in' type='checkbox'><label for='data-ce80b180-3989-4e59-892f-4bc60f36772f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>lat_bnds</dd><dt><span>units :</span></dt><dd>degrees_north</dd><dt><span>axis :</span></dt><dd>Y</dd><dt><span>long_name :</span></dt><dd>Latitude</dd><dt><span>standard_name :</span></dt><dd>latitude</dd></dl></div><div class='xr-var-data'><pre>array([-89.14152, -88.02943, -86.91077, -85.79063, -84.66992, -83.54895,\n",
       "       -82.42782, -81.30659, -80.18531, -79.06398, -77.94262, -76.82124,\n",
       "       -75.69984, -74.57843, -73.45701, -72.33558, -71.21414, -70.09269,\n",
       "       -68.97124, -67.84978, -66.72833, -65.60686, -64.4854 , -63.36393,\n",
       "       -62.24246, -61.12099, -59.99952, -58.87804, -57.75657, -56.63509,\n",
       "       -55.51361, -54.39214, -53.27066, -52.14917, -51.02769, -49.90621,\n",
       "       -48.78473, -47.66325, -46.54176, -45.42028, -44.29879, -43.17731,\n",
       "       -42.05582, -40.93434, -39.81285, -38.69137, -37.56988, -36.44839,\n",
       "       -35.32691, -34.20542, -33.08393, -31.96244, -30.84096, -29.71947,\n",
       "       -28.59798, -27.47649, -26.355  , -25.23351, -24.11203, -22.99054,\n",
       "       -21.86905, -20.74756, -19.62607, -18.50458, -17.38309, -16.2616 ,\n",
       "       -15.14011, -14.01862, -12.89713, -11.77564, -10.65415,  -9.53266,\n",
       "        -8.41117,  -7.28968,  -6.16819,  -5.0467 ,  -3.92521,  -2.80372,\n",
       "        -1.68223,  -0.56074,   0.56074,   1.68223,   2.80372,   3.92521,\n",
       "         5.0467 ,   6.16819,   7.28968,   8.41117,   9.53266,  10.65415,\n",
       "        11.77564,  12.89713,  14.01862,  15.14011,  16.2616 ,  17.38309,\n",
       "        18.50458,  19.62607,  20.74756,  21.86905,  22.99054,  24.11203,\n",
       "        25.23351,  26.355  ,  27.47649,  28.59798,  29.71947,  30.84096,\n",
       "        31.96244,  33.08393,  34.20542,  35.32691,  36.44839,  37.56988,\n",
       "        38.69137,  39.81285,  40.93434,  42.05582,  43.17731,  44.29879,\n",
       "        45.42028,  46.54176,  47.66325,  48.78473,  49.90621,  51.02769,\n",
       "        52.14917,  53.27066,  54.39214,  55.51361,  56.63509,  57.75657,\n",
       "        58.87804,  59.99952,  61.12099,  62.24246,  63.36393,  64.4854 ,\n",
       "        65.60686,  66.72833,  67.84978,  68.97124,  70.09269,  71.21414,\n",
       "        72.33558,  73.45701,  74.57843,  75.69984,  76.82124,  77.94262,\n",
       "        79.06398,  80.18531,  81.30659,  82.42782,  83.54895,  84.66992,\n",
       "        85.79063,  86.91077,  88.02943,  89.14152])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>0.0 1.125 2.25 ... 357.8 358.9</div><input id='attrs-c7b98a9c-7431-4ec7-bb10-00ea44e84571' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-c7b98a9c-7431-4ec7-bb10-00ea44e84571' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f5a8286f-028c-4e01-b6e1-6cfdd21fa412' class='xr-var-data-in' type='checkbox'><label for='data-f5a8286f-028c-4e01-b6e1-6cfdd21fa412' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>bounds :</span></dt><dd>lon_bnds</dd><dt><span>units :</span></dt><dd>degrees_east</dd><dt><span>axis :</span></dt><dd>X</dd><dt><span>long_name :</span></dt><dd>Longitude</dd><dt><span>standard_name :</span></dt><dd>longitude</dd></dl></div><div class='xr-var-data'><pre>array([  0.   ,   1.125,   2.25 , ..., 356.625, 357.75 , 358.875])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>type</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>|S4</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-56ea4bd8-6873-4ca6-98d4-296b530522f9' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-56ea4bd8-6873-4ca6-98d4-296b530522f9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3225d6f9-57e4-4ecd-af09-534471694e05' class='xr-var-data-in' type='checkbox'><label for='data-3225d6f9-57e4-4ecd-af09-534471694e05' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>Land area type</dd><dt><span>standard_name :</span></dt><dd>area_type</dd></dl></div><div class='xr-var-data'><pre>array(b&#x27;land&#x27;, dtype=&#x27;|S4&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-695edd54-0e27-421b-96dc-aa64f0459ea6' class='xr-section-summary-in' type='checkbox'  checked><label for='section-695edd54-0e27-421b-96dc-aa64f0459ea6' class='xr-section-summary' >Attributes: <span>(9)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>standard_name :</span></dt><dd>land_area_fraction</dd><dt><span>long_name :</span></dt><dd>Percentage of the grid  cell occupied by land (including lakes)</dd><dt><span>comment :</span></dt><dd>Please express &#x27;X_area_fraction&#x27; as the percentage of horizontal area occupied by X.</dd><dt><span>units :</span></dt><dd>%</dd><dt><span>original_name :</span></dt><dd>lfrac</dd><dt><span>original_units :</span></dt><dd>1</dd><dt><span>history :</span></dt><dd>2019-05-24T05:34:57Z altered by CMOR: Converted units from &#x27;1&#x27; to &#x27;%&#x27;. 2019-05-24T05:34:57Z altered by CMOR: Treated scalar dimension: &#x27;type&#x27;. 2019-05-24T05:34:57Z altered by CMOR: replaced missing value flag (-9.99e+33) with standard missing value (1e+20).</dd><dt><span>cell_methods :</span></dt><dd>area: mean</dd><dt><span>cell_measures :</span></dt><dd>area: areacella</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.DataArray 'sftlf' (lat: 160, lon: 320)>\n",
       "array([[100., 100., 100., ..., 100., 100., 100.],\n",
       "       [100., 100., 100., ..., 100., 100., 100.],\n",
       "       [100., 100., 100., ..., 100., 100., 100.],\n",
       "       ...,\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.],\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.],\n",
       "       [  0.,   0.,   0., ...,   0.,   0.,   0.]], dtype=float32)\n",
       "Coordinates:\n",
       "  * lat      (lat) float64 -89.14 -88.03 -86.91 -85.79 ... 86.91 88.03 89.14\n",
       "  * lon      (lon) float64 0.0 1.125 2.25 3.375 4.5 ... 355.5 356.6 357.8 358.9\n",
       "    type     |S4 ...\n",
       "Attributes:\n",
       "    standard_name:   land_area_fraction\n",
       "    long_name:       Percentage of the grid  cell occupied by land (including...\n",
       "    comment:         Please express 'X_area_fraction' as the percentage of ho...\n",
       "    units:           %\n",
       "    original_name:   lfrac\n",
       "    original_units:  1\n",
       "    history:         2019-05-24T05:34:57Z altered by CMOR: Converted units fr...\n",
       "    cell_methods:    area: mean\n",
       "    cell_measures:   area: areacella"
      ]
     },
     "execution_count": 105,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# get the correct land mask\n",
    "fils_land = list(Path(\"/Volumes/Jedha/CMIP6\").glob(f'sftlf_fx_{model_name}_{experiment_name}_{member}_*.nc'))\n",
    "if len(fils_land) > 1:\n",
    "    print(fils_land)\n",
    "ds_land = xr.open_dataset(fils_land[0])\n",
    "land = ds_land['sftlf']\n",
    "land"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 106,
   "metadata": {},
   "outputs": [],
   "source": [
    "pr_tropics = pr.sel(lat=slice(-30,30))\n",
    "pr_tropics_land = pr_tropics.where(land > 0)\n",
    "pr_tropics_ocean = pr_tropics.where(land <= 0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 107,
   "metadata": {},
   "outputs": [],
   "source": [
    "pr_nh = pr.sel(lat=slice(30,None))\n",
    "pr_sh = pr.sel(lat=slice(None,-30))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 108,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "tropical land\n",
      "tropical ocean\n",
      "NH\n",
      "SH\n"
     ]
    }
   ],
   "source": [
    "# flatten the array and make histogram using specified bins\n",
    "\n",
    "prbins = np.linspace(0, 150, 51)\n",
    "\n",
    "print(\"tropical land\")\n",
    "h_trop_land = np.histogram(pr_tropics_land.values.flatten(), density=False, bins=prbins)\n",
    "print(\"tropical ocean\")\n",
    "h_trop_ocean = np.histogram(pr_tropics_ocean.values.flatten(), density=False, bins=prbins)\n",
    "print(\"NH\")\n",
    "h_nh = np.histogram(pr_nh.values.flatten(), density=False, bins=prbins)\n",
    "print(\"SH\")\n",
    "h_sh = np.histogram(pr_sh.values.flatten(), density=False, bins=prbins)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 109,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAk10lEQVR4nO3df5xd853H8dfbCBGRlCTbLRETojQhIoZqa4n6FUsEqyVi0VXpKqstsapaiaput4m1LF1NSYNmKW2qorF+tRpdrEyIH/EzIWVCiagQpGF99o9z5rq5uXdm7uSeuffOvJ+Px32Y+z0/7ufemPuZ74/zOYoIzMzMADaqdgBmZlY7nBTMzCzHScHMzHKcFMzMLMdJwczMcpwUzMwsx0nBzMxynBSsW5G0TNJaSQML2hdJCkmNkmal+6yW9IakuyTtnLfvyZL+0MZrTJX0fnp86+PNvO3j09d7S9Lrku6R1Jh3bEg6s+CcX0/bp6bP907jekPSCkk3S/pEO+/9AElPS3pX0u8kbVfGR2cGOClY9/QCMKH1iaRdgc0K9vlhRPQFtgGWA9eU+Ro/j4i+eY+Ppa81DLgOOBvoDwwFfgR8mHfss8BJBec7MW1vtSUwA2gEtgPeBn5aKpg0Cc4BvgNsBTQDPy/zPZk5KVi3dD3Jl2yrk0i+qNcTEe8BNwGjKvTao4AXIuKeSLwdEb+MiBfz9lkA9JE0AiD972Zpe2tct0fEzRHxVkS8C1wBfK6N1z0aWJweswaYCuyW3wMy6wgnBeuOHgT6SfqUpAbgWOBnxXaUtDlJr2JJhV77YWBnSZdK2l9S3xL75Seukkkrz77A4ja2jwAebX0SEe8AS9N2sw5zUrDuqvVL9yDgaZIhonyT03mAt4F9gL8v8/xflPRm3uN3ABHxPDCGZFjqJuD1dA6jMDn8DJggqRdwHCWSFoCkkcAFwDltxNMXWFXQtgrYooz3ZOakYN3W9cDxwMkU/yt8ejoP0Ai8B+xU7CSSJuZNJt+et+mmiPhY3mP/1g0R8WBEfDEiBgF/Q/JX/vn5502Hk5YA3weei4iXSrz+MOB24GsRcV/aNiR/kjvddTXQr+DwfiRJz6zDnBSsW4qIP5JMOP8tyQRsqf1eBL4GXCapcDKaiJidN5l8aCfiWJC+/i5FNrdOSBcdOkpXD90NXBQR1+fHnD/JnTYvBnbLO3ZzYAfaHnIyW4+TgnVnpwCfT8fXS4qIu4CXgUkb+oKS9pF0qqS/Sp/vDBxBMs9R6OfAwSTDTIXn2Qb4LXBlRFzVgZf+FbCLpL+T1JtkuOmxiHi6k2/FeignBeu2ImJpRDR3cPdpwD9L2rSD+x9bcJ3C6jQRvEmSBB5Ph3b+m+QL+4dF4nsvIu5OV0AV+jKwPTClyFDReiJiBfB3wMXAn4FPk8xVmJVFvsmOmZm1ck/BzMxynBTMzCzHScHMzHKcFMzMLGfjagfQStJwknotK4F7IuIX7R0zcODAaGxszDgyM7PuZeHCha+nF1euJ9OkIGkmcDjwWkTsktc+FrgMaACujogfAIcC/xER90m6FWg3KTQ2NtLc3NEVh2ZmBiDpj6W2ZT18NAsYWxBMA3AlSRIYTlL/ZThJWYLjJE0DBmQcl5mZFZFpUoiI+cAbBc17AUsi4vmIWAvcCIyPiNci4nTgm8DrWcZlZmbFVWNOYRsgv/hXC/Dp9M5U3wI2J7m6tChJk0jLEQwZMiS7KM3MeqBqJAUVaYuIWEYHas9ExAySO1LR1NTky7HNepD333+flpYW1qxZU+1Q6kLv3r0ZPHgwvXr16vAx1UgKLcC2ec8HkxQj6zBJ44Bxw4YNq2RcZlbjWlpa2GKLLWhsbEQq9veltYoIVq5cSUtLC0OHDu3wcdW4TmEBsKOkoZI2ISnadWs5J4iIuRExqX///pkEaGa1ac2aNQwYMMAJoQMkMWDAgLJ7VZkmBUk3AA8AO0lqkXRKRHwAnAHcATxFcrOSsmq+SxonacaqVYU3mjKz7s4JoeM681llOnwUERNKtM8D5m3AeecCc5uamk7t7DnMzGx9NXNFczk8p2BmAI9/dbuKnm/XH5W8pitHEmeddRaXXHIJANOnT2f16tVMnTqVqVOn0rdvXyZPnpzbv/Ui24EDB1Y01qzUZe0jzymYWbVsuummzJkzh9df756XU9VlUjAzq5aNN96YSZMmcemll1Y7lEzUZVLwRLOZVdPpp5/O7NmzKfYddOmllzJq1Kjc4+WXy1pxX3V1mRQ8fGRm1dSvXz9OPPFELr/88vW2feMb32DRokW5x9Zbb12FCDuvLieaK+GQi35TtP2O7xzWxZGYWT36+te/zujRo/nSl75U7VAqqi57Ch4+MrNq22qrrfjiF7/INddcU+1QKqouewq+TsHMoGNLSLN09tlnc8UVV1Q1hkqry6RgZlYtq1evzv388Y9/nHfffTf3fOrUqevtv2zZsi6IqnJ6bFJYstmZJbZ4TsHMei7PKZiZWU5dJgUvSTUzy0aPHT468Zzni2+YXLzZzKwnqMuegpmZZcNJwczMcupy+Mils80M4MIK33BnSnTstu8XX3wx//Vf/0VDQwMbbbQRP/7xjzn33HOZPn06TU1NQLIU9fDDD+eJJ56oaIxZq8uk4IvXzKxaHnjgAW677TYefvhhNt10U15//XXWrl1b7bAqpi6TgplZtbzyyisMHDiQTTfdFKBubp7TUU4KZmZlOPjgg/nud7/LJz/5SQ488ECOPfZY9ttvPwAmTpzIZpttBsDatWvZaKP6m7atv4jNzKqob9++LFy4kBkzZjBo0CCOPfZYZs2aBcDs2bNzJbPnzev0beiryj0FM7MyNTQ0MGbMGMaMGcOuu+7KtddeW+2QKsY9BTOzMjzzzDM899xzueeLFi1iu+22q2JElVWXPQUvSTUz6PgS0kpavXo1//RP/8Sbb77JxhtvzLBhw5gxYwbHHHNMl8eShbpMCl6SambVsscee3D//fev137vvfeu87yxsbHurlGAOk0KlXD0aUOqHYKZWc3xnIKZmeU4KZiZWY6TgpmZ5TgpmJlZjpOCmZnl9NjVR2ZW/3aYvkNFz7d08tI2t69cuZIDDjgAgD/96U80NDQwaNAgAB566CE22WQTZs2aRXNzM1dccUVFY2tsbKS5uTnzAnw1kxQkDQGuAF4Hno2IH1Q5JDOzdQwYMIBFixYBMHXqVPr27cvkyR/dw/eDDz6oUmSVk2lSkDQTOBx4LSJ2yWsfC1wGNABXpwngk8BvIuLHkq7LMi4zs0o5+eST2WqrrXjkkUcYPXo0u+66a9H9TjvtNBYsWMB7773HMcccw4UXXggkPYCTTjqJuXPn8v7773PzzTez8847s3LlSiZMmMCKFSvYa6+9iC66ejvrOYVZwNj8BkkNwJXAocBwYIKk4cAjwHGSfgv8LuO4zMwq5tlnn+Xuu+/mkksuKbnPxRdfTHNzM4899hi///3veeyxx3LbBg4cyMMPP8xpp53G9OnTAbjwwgvZZ599eOSRRzjiiCN48cUXM38fkHFSiIj5wBsFzXsBSyLi+YhYC9wIjAe+BEyJiM8Dh5U6p6RJkpolNa9YsSKr0M3MOuwLX/gCDQ0Nbe5z0003MXr0aHbffXcWL17Mk08+mdt29NFHA0kJjWXLlgEwf/58TjjhBAAOO+wwttxyy2yCL1CN1UfbAC/lPW9J2/4bOFPSVcCyUgdHxIyIaIqIptYJHjOzatp8883b3P7CCy8wffp07rnnHh577DEOO+ww1qxZk9veehe3hoaGdeYlVOF7UHdENZJCsXcZEfFERBwTEf8YEZOL7PPRCaRxkmasWrUqoxDNzCrnrbfeYvPNN6d///68+uqr3H777e0es++++zJ79mwAbr/9dv785z9nHSZQndVHLcC2ec8HAy+XcwJXSTUzaH8JabXMmjWLW265Jff8wQcfZPfdd2fEiBFsv/32fO5zn2v3HFOmTGHChAmMHj2a/fbbjyFDuqaIp7Ke0ZbUCNzWuvpI0sbAs8ABwHJgAXB8RCwu45yt91M4Nf9mF+V4/KvFb4qx64/+2KnzmVn2nnrqKT71qU9VO4y6Uuwzk7QwIpqK7Z/p8JGkG4AHgJ0ktUg6JSI+AM4A7gCeAm4qJyFA0lOIiEn9+/evfNBmZj1YpsNHETGhRPs8oD7vam1m1o3VZe0jTzSbmWWjLpOCh4/MzLJRM7WPypE30Vzxcx9y0W+Ktt/xnZLX05mZdRvuKZiZWU5d9hTMzKB0z76z2hsR6EjpbIBbbrmFCy64gLVr19KrVy8uuugijjzyyIrGmhUnBTOzDupI6ezFixczefJk7rrrLoYOHcoLL7zAQQcdxPbbb8/IkSOrFHnH1eXwkVcfmVmtOPnkkznrrLPYf//9Offcc5k+fTrf+ta3GDp0KABDhw7lvPPOY9q0aQAsWbKEAw88kN12243Ro0ezdGlyVfa0adPYc889GTlyJFOmTMmd/8gjj2SPPfZgxIgRzJgxI9fet29fzj//fHbbbTf23ntvXn311Yq8n7pMClnOKSzZ7MyiDzOzUvJLZy9evJg99thjne1NTU0sXpxcoztx4kROP/10Hn30Ue6//34+8YlPcOedd/Lcc8/x0EMPsWjRIhYuXMj8+fMBmDlzJgsXLqS5uZnLL7+clStXAvDOO++w99578+ijj7Lvvvvyk5/8pCLvpS6TgplZLckvnR0R61U3bW17++23Wb58OUcddRQAvXv3pk+fPtx5553ceeed7L777owePZqnn36a1hI+l19+ea438NJLL+XaN9lkEw4//HBg3ZLbG6ou5xSyXJJ64jnPF9/QZt1WM+vJ8ktnjxgxgubm5nXmDx5++GGGDx9e8u5pEcF5553HV77ylXXa7733Xu6++24eeOAB+vTpw5gxY3Ilt3v16pVLPoUltzdEXfYUvCTVzGrV5MmT+Zd/+ZfcX+7Lli3j+9//PmeffTb9+vVj8ODBuQqqf/nLX3j33Xc55JBDmDlzJqtXrwZg+fLlvPbaa6xatYott9ySPn368PTTT/Pggw9mHn9d9hTMzKA2LyodNWoU//qv/8q4ceN4//336dWrFz/84Q8ZNWoUANdffz1f+cpXuOCCC+jVqxc333wzBx98ME899RSf+cxngGQS+Wc/+xljx47lqquuYuTIkey0007svffemcefeensLDU1NUVzc3Onji1VOnvOfxa/D+qUOv6czLoLl84uX02VzjYzs/ripGBmZjl1OaeQ5eqjUtq6nL4WxzXNuqtiSz6tuM5MD9RlT8Grj8x6pt69e7Ny5cpOfdn1NBHBypUr6d27d1nH1WVPwcx6psGDB9PS0sKKFSuqHUpd6N27N4MHDy7rGCcFM6sbvXr1ytUUsmzU5fCRmZllw0nBzMxyPHzUQW1XSvXqIzPrHuqyp+D7KZiZZaMuewoRMReY29TUdGqlz330aUOKtl9X6RcyM6tBdZkUqqFkSW1wWW0z6zacFCqg1NXOvtLZzOpNXc4pmJlZNpwUzMwsx0nBzMxynBTMzCzHScHMzHJqZvWRpL8BJpLENDwiPlvlkMzMepxMk4KkmcDhwGsRsUte+1jgMqABuDoifhAR9wH3SToSWJBlXJVWugSGl6SaWX3JuqcwC7iCvAuCJTUAVwIHAS3AAkm3RsST6S7HA1/OOK4u4esXzKzeZDqnEBHzgTcKmvcClkTE8xGxFrgRGA8gaQiwKiLeKnVOSZMkNUtq9o02zMwqqxoTzdsAL+U9b0nbAE4BftrWwRExIyKaIqJp0KBBGYVoZtYzVWOiudgdtwMgIqZ06ATSOGDcsGHDKhlXp5Wqi3T/d7s4EDOzDVSNnkILsG3e88HAy+WcICLmRsSk/v37VzQwM7Oerho9hQXAjpKGAsuB40gmlzusGj2FUiW1Aeb854tdFoeZWZYy7SlIugF4ANhJUoukUyLiA+AM4A7gKeCmiFhcznndUzAzy0amPYWImFCifR4wr7PnrbU5hXJ5qaqZ1aqauaK5HFneea2SSl3UNuy9y7s4EjOzjnHtIzMzy2kzKUg6o6sCKYekcZJmrFq1qtqhmJl1K+31FP6hS6Iokyeazcyy4eEjMzPLaW+ieaSkYnWIBERE9MsgpnbV++ojM7NapYgovVF6JCJ278J4ytLU1BTNzc2dOvbxr25XsTjKvXjt/u/eVtb+XqpqZpUkaWFENBXb5uEjMzPLaW/46OYuiaLOlSqBUaoH4esXzKxWtZcUekm6oMS2iIiLKh1QR3hOwcwsG+0NH60G3il4BMl9D87NNrTSvCTVzCwbbfYUIuKS1p8lbQF8jeTahRuBS0odZ2Zm9and2keStgLOAiYC1wKjI+LPWQdmZmZdr82kIGkacDQwA9g1IlZ3SVS2DldVNbOu0t6cwtnA1sC3gZclvZU+3i5xUVuXcO0jM7NstDenUJPXMdRL6exSSt3T+bppXqpqZtXV3vDRVgVNAbwZbV0GbWZmdau9ieaFJIlAeW1bSFoEfDkilmUUV7dQ7kVtZmbV1t7w0dBi7ZKOBq4CxmYRlJmZVUenbscZEXMkfbvSwVh5vCrJzCqtUxPJkvp29lgzM6td7U00n1WkeUvgCOCKTCLqANc+MjPLRnvDR1sUPA/gT8AJEfF4NiG1r96XpJbLVVXNrKu0N9F8YevP6ZBRRMQ7mUfVQ5W+fmH7Lo7EzHqqducFJJ0m6UXgj8CLkv4o6avZh2ZmZl2tvTmFbwOfBcZExPNp2/bAZZK2iojvdUGM3U7W1y94VZKZdVZ7PYW/B45uTQgA6c9fBE7MMjAzM+t67Q4fRcSaIm3vAR9mEpGZmVVNe0mhRdIBhY1p2yvZhGRmZtXS3pLUM4FfS/oDH9VB2hP4HDA+49isHV6qamaV1l5S+AtwMvBJYARJYbz5wDXAesNKG0LSRsBFQD+gOSKureT561mllqp6AtrM2tPe8NG/A29FxMyIODsizoqIa4B3021tkjRT0muSnihoHyvpGUlLJH0zbR4PbAO8D7SU+0bMzGzDtZcUGiPiscLGiGgGGjtw/lkUVFKV1ABcCRwKDAcmSBoO7AQ8EBFnAad14NxmZlZh7Q0f9W5j22btnTwi5ktqLGjeC1iSd93DjSS9hJeAtek+/1fqnJImAZMAhgwpvt6/Xvn+C2ZWbe31FBZIWq++kKRTSCaeO2MbkgTQqiVtmwMcIuk/SOYtioqIGRHRFBFNgwYN6mQIZmZWTHs9ha8Dv5I0kY+SQBOwCXBUJ19TRdoiIt4FTunQCVwltU1elWRmndVeQbxXgc9K2h/YJW3+TUT8dgNeswXYNu/5YODlck7Q06qkZs2rksysVYfuvBYRvwN+V6HXXADsKGkosBw4Dji+nBO4p2Bmlo1O3Y6zoyTdAIwBBkpqAaZExDWSzgDuABqAmRGxuJzzuqeQcKltM6u0TJNCREwo0T4PmJfla5uZWfkyTQpZ6WnDR16qamZdpS6TgoePOqfcVUmegDbredotnW1mZj1HXSYFSeMkzVi1alW1QzEz61Y8fNQNeVWSmXVWXfYUzMwsG3XZU+hpq49qjSegzbqvukwKHj5KVGqpqmslmVkrDx+ZmVmOk4KZmeXUZVLwklQzs2x4TqEHyXqpaqkJaPAktFm9qMukYF3DE9BmPY+TQjfkAnpm1ll1OadgZmbZqMuk4IlmM7Ns1OXwkSea64+vgjarD3WZFKyyXEDPzFo5KVjZvCrJrPuqyzkFMzPLhnsKPUgtLlX1XINZbXFPwczMcuqyp+D7KXQNT0Cb9Tx1mRS8JLU2lZqABk9Cm9ULDx+ZmVlOXfYUrPvzBLRZdTgpWE2uSjKz6vDwkZmZ5binYGXrzKokXwVtVh/cUzAzsxz3FKyueALaLFs1kxQkjQEuAhYDN0bEvdWMxzwBbdYTZTp8JGmmpNckPVHQPlbSM5KWSPpm2hzAaqA30JJlXGZmVlzWPYVZwBXAda0NkhqAK4GDSL78F0i6FbgvIn4v6ePAvwETM47NKqzUBDS4NIZZvcg0KUTEfEmNBc17AUsi4nkASTcC4yPiyXT7n4FNS51T0iRgEsCQIcWHN6x+eFWSWW2pxpzCNsBLec9bgE9LOho4BPgYSe+iqIiYAcwAaGpqiuzCtHriCWizyqhGUlCRtoiIOcCcDp3AVVLNzDJRjaTQAmyb93ww8HI5J3CV1OryqiSz7qsaSWEBsKOkocBy4Djg+HJO4J5C/fG9GczqQ6ZJQdINwBhgoKQWYEpEXCPpDOAOoAGYGRGLyzmvewrdX6UmoD3XYFaerFcfTSjRPg+Y19nzuqdgZpaNuqx9FBFzI2JS//79qx2KmVm3UjNlLqz+lZqABk9Cm9WLuuwpSBonacaqVauqHYqZWbeiiPq9/qupqSmam5s7dezjX92uwtFYW8rtKZS7KqlSV0B7Atp6AkkLI6Kp2La67CmYmVk26jIpePjIzCwbdZkUvPrIzCwbXn1kXSLr0hi+2M2sMpwUrCa5LIZZddTl8JHnFMzMslGXScFzCmZm2fDwkVkHeK7BegonBauqciegy51r8O0+zcpTl8NHZmaWjbpMCp5oNjPLRl0OH/kmO1YrPNdg3U1d9hTMzCwbddlTsO7PE9Bm1eGegpmZ5binYJaBUnMN4PkGq23uKZiZWU5d9hQkjQPGDRs2rNqhWDfjuQbr6Xw7TusW6uV2n23xsJJ1Fd+O08zMOsRJwczMcupyTsGsq3muwXoKJwXrFrK+2K0ruGSG1QIPH5mZWY6TgpmZ5Xj4yKzGeVjJulJNJQVJmwPzgSkRcVu14zFrjyegrbvJNClImgkcDrwWEbvktY8FLgMagKsj4gfppnOBm7KMyXqW7jABbdaVsu4pzAKuAK5rbZDUAFwJHAS0AAsk3QpsDTwJ9M44JrNuwcNKloVMk0JEzJfUWNC8F7AkIp4HkHQjMB7oC2wODAfekzQvIj7MMj6zrJQaVgIPLVltq8acwjbAS3nPW4BPR8QZAJJOBl4vlRAkTQImAQwZUnxowMzMOqcaSUFF2nJV+SJiVlsHR8QMYAYkBfEqGplZCfU01+BhJdsQ1UgKLcC2ec8HAy+XcwKXzrYNVe4EtFlPUY2ksADYUdJQYDlwHHB8OSeIiLnA3KamplMziM8sU17GarUs6yWpNwBjgIGSWkiuP7hG0hnAHSRLUmdGxOIyz+ueglmZPKxkHZH16qMJJdrnAfM24LzuKVhNqKe5BrOOqKkrmjvKPQXLiucarKery6TgnoJ1R9Waa/CwkuVzlVQzM8upy56Ch4/MsuceRM9Ul0nBw0fW1SpVWA88CW21TRH1e1GwpBXAHzt5+EDg9QqGk4Vaj7HW4wPHWAm1Hh/Ufoy1Ft92ETGo2Ia6TgobQlJzRDRVO4621HqMtR4fOMZKqPX4oPZjrPX48nmi2czMcpwUzMwspycnhRnVDqADaj3GWo8PHGMl1Hp8UPsx1np8OT12TsHMzNbXk3sKZmZWwEnBzMxyemRSkDRW0jOSlkj6Zg3Es62k30l6StJiSV9L27eSdJek59L/blnlOBskPSLpthqN72OSfiHp6fSz/EwNxviN9N/4CUk3SOpd7RglzZT0mqQn8tpKxiTpvPR35xlJh1Qpvmnpv/Njkn4l6WPViq9UjHnbJksKSQOrGWNH9bikIKkBuBI4FBgOTJA0vLpR8QFwdkR8CtgbOD2N6ZvAPRGxI3BP+ryavgY8lfe81uK7DPjviNgZ2I0k1pqJUdI2wJlAU0TsQnI/keNqIMZZwNiCtqIxpf9fHgeMSI/5Ufo71dXx3QXsEhEjgWeB86oYX6kYkbQtcBDwYl5btWLskB6XFIC9gCUR8XxErAVuBMZXM6CIeCUiHk5/fpvky2ybNK5r092uBY6sSoCApMHAYcDVec21FF8/YF/gGoCIWBsRb1JDMaY2BjaTtDHQh+RWtFWNMSLmA28UNJeKaTxwY0T8JSJeAJaQ/E51aXwRcWdEfJA+fZDktr5Via9UjKlLgX8m7z701Yqxo3piUtgGeCnveUvaVhMkNQK7A/8LfDwiXoEkcQB/VcXQ/p3kf+4P89pqKb7tgRXAT9MhrqslbV5LMUbEcmA6yV+NrwCrIuLOWooxT6mYavH35x+A29OfayY+SUcAyyPi0YJNNRNjMT0xKahIW02sy5XUF/gl8PWIeKva8bSSdDjwWkQsrHYsbdgYGA38Z0TsDrxD9Yez1pGOy48HhgJbA5tLOqG6UZWtpn5/JJ1PMvw6u7WpyG5dHp+kPsD5wAXFNhdpq4nvIOiZSaEF2Dbv+WCSLnxVSepFkhBmR8SctPlVSZ9It38CeK1K4X0OOELSMpLhts9L+lkNxQfJv2tLRPxv+vwXJEmilmI8EHghIlZExPvAHOCzNRZjq1Ix1czvj6STgMOBifHRBVe1Et8OJMn/0fT3ZjDwsKS/pnZiLKonJoUFwI6ShkrahGTC59ZqBiRJJGPhT0XEv+VtuhU4Kf35JODXXR0bQEScFxGDI6KR5PP6bUScUCvxAUTEn4CXJO2UNh0APEkNxUgybLS3pD7pv/kBJPNHtRRjq1Ix3QocJ2lTSUOBHYGHujo4SWOBc4EjIuLdvE01EV9EPB4RfxURjenvTQswOv3/tCZiLCkietwD+FuSFQtLgfNrIJ59SLqPjwGL0sffAgNIVn48l/53qxqIdQxwW/pzTcUHjAKa08/xFmDLGozxQuBp4AngemDTascI3EAyx/E+yZfXKW3FRDIsshR4Bji0SvEtIRmXb/19uapa8ZWKsWD7MmBgNWPs6MNlLszMLKcnDh+ZmVkJTgpmZpbjpGBmZjlOCmZmluOkYGZmOU4K1mGS/k/SorTC583pVZuF7XNbK1ZKapT0Xrqt9XFiuq2vpB9LWppWDZ0v6dPpttVFjn9S0lWSNkrbn5B0SN55V6cVJxdJui49/qi0OuXO6fP/Tbe/KGlF3rGNkpa1VrGUNFjSr5VUCF0q6bL0mhYkjUnPOS7vc7lN0pgin9csSS+kr/GopAMKtn9D0hpJ/SUNyIvnT5KW5z3fJO8zbn2sd7W2pHvTz+CIDf7HLkJSL0nrXdUuaaqkyZ0857T0/XbqeKu8jasdgNWV9yJiFICk2cA/Av9W0H4tcDpwcXrM0tZtBa4GXgB2jIgPJW0PfKrIfksjYpSSAnK/JSnM1lo88A7gjvR17wUmR0Rz3rETgD+QXHA3NSJak87JJJVKz2jdMbmWLHch4RySchnjlVSvnJG+n3PS3VtI1pnPbevDSp0TEb+QtH96nh0L4lsAHBURs0ius0DSVGB1REzPi++9Ep9joYkFn0El7QPcX8kTRsQ5kt6p5Dltw7inYJ11HzCsSPsDtFPcS9IOwKeBb0fEhwCRVK39TaljIqmIeX+J1yz2Gn1JynOcQpIUOurzwJqI+Gn6uv8HfAP4h9aeEfAosErSQWWcd53PJf0M+gLfJkkOFZf2HC5Ne2FPSdpT0py0B/S9dJ9GJfcluDrtfc2WdKCk/0n3y6/eOZa08Jyk89Neyd3ATnmveaqkBWnP6JdKrt7eIu0x9Ur36Zf2zHpl8b5twzgpWNnSv9oPBR4vaG8gKd2QXzZkh4Jhj78hqSO/KP3C7ehr9knP/Xh7+6aOJLm3wrPAG5JGd/C4EcA6QySRFCd8kXUT0vdIvtA7aizJVdatJpBcBXsfsJOk9iqjblbwOR7bwdddGxH7AleRlKo4HdgFOFnSgHSfYST3ohgJ7AwcT9IrmAx8K+9c+wP3StqDJNHuDhwN7Jm3z5yI2DMiWu9ncUok5eDvJSm9TnrsLyOp/2Q1xsNHVo7NJC1Kf76P9N4Fee2NJF+od+Uds97wUZlj3juk5w7g1xFxu5Ly4u2ZQFLuG5IifhNIh53aIYpXrFynPSLuk0Sa5NoyTdIPSUpP753XfhzJsNGHkuYAXyC5+VMpHR0+KtSaoB8HFkdaDlvS8yRF2d4kKdL3eNq+mOTmOiHpcZJ/UyRtDbwREe+m7/lXkdYckpT/R8AuaS/kYyQ9oTvS9qtJSq/fAnwJOLUT78W6gJOClaPUF9N76bh/f+A2kr9GL2/jPIuB3SRt1Dp81IZScxIlpX8Bf57kCypI7nAWkv452q/rshj4u4Lz9SP5Al1KUhOo1cUkcwsfUNo5JHMUZ5LcrGYPSSNJ5hbuSucyNgGep+2k0Fl/Sf/7Yd7Prc83LtincL/8fQ7loy94KF3qeRZwZEQ8ms7djAGIiP9Jh6r2AxoiYr3bVlpt8PCRVUxErCL58pvc1nhxRCwlKVx3YTqxi6QdJVXqDnjHANdFxHaRVKnclmRSe58OHHsP0EcfrZJqAC4BZsW61TiJ5AY5W5Lc+rOkNPFdBmyk5H68E0gmvhvTx9bANpK2K+9tdqncfAIwHzhK0maStgDG5e23BfBK+u8/seAc15EMmf0062Ct85wUrKIi4hGSidjWyd3COYUz0/YvA38NLEmHKX5C5WrKTwB+VdD2S5Kx8jalPYmjgC9Ieo6kmu4a1h1bz3cxH90Ksr3zfo9kCOW4IvH9irYnxAvnFH7Q3mtWSpoYd4yIpwEiuXXsz0mqk/6SZCix1XdI7hp4F0k12HyzSZLoDRmHbBvAVVLNuokSy3Ircd59gBMi4h838DzHAOMj4u8L2qdSsATXqsc9BbPu4w1gVpkT+e2KiD9UICH8B/AD4KKC9mnACSS3T7Ua4J6CmZnluKdgZmY5TgpmZpbjpGBmZjlOCmZmluOkYGZmOf8PBswsN5q56G0AAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "bctr = prbins[0:-1] + 0.5*(prbins[1:]-prbins[0:-1])\n",
    "width = 3\n",
    "fig, ax = plt.subplots()\n",
    "ax.bar(bctr, h_nh[0], color='chocolate', width=width, label='NH')\n",
    "ax.bar(bctr, h_sh[0], color='maroon', bottom=h_nh[0], width=width, label='SH')\n",
    "ax.bar(bctr, h_trop_land[0], color='forestgreen', bottom=h_nh[0]+h_sh[0], width=width, label='TrLand')\n",
    "ax.bar(bctr, h_trop_ocean[0], color='steelblue', bottom=h_nh[0]+h_sh[0]+h_trop_land[0], width=width, label='TrOcean')\n",
    "ax.semilogy()\n",
    "# ax.set_ylim([0, 1e5])\n",
    "ax.legend()\n",
    "ax.set_ylabel('COUNT')\n",
    "ax.set_xlabel('PRECIPITATION RATE [mm/day]')\n",
    "ax.set_title(model_name)\n",
    "fig.savefig(f\"/Users/brianpm/Dropbox/Manuscripts/M_cre_pex/REVISION1/pr_histogram_{model_name}_{experiment_name}.png\", dpi=144)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 110,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYIAAAEWCAYAAABrDZDcAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/Il7ecAAAACXBIWXMAAAsTAAALEwEAmpwYAAAjTklEQVR4nO3de5xd473H8c/XSCQRkQTVykWCuN9KXFrqrqJUcNpKqFI0pW5VFAelh/alxVFFRUqaUgdFEJcgbkWVJlpBxCUJYlBycU0I4Xf+WGtiZ2fv2ZPJXrP3zPq+X695zaxnPXvt397J7N88z7PWbykiMDOz/Fqu1gGYmVltORGYmeWcE4GZWc45EZiZ5ZwTgZlZzjkRmJnlnBOBmVnOORFYuyfpFUmfSFq1qP0pSSFpgKQxaZ8PJc2VNEHS+gV9D5X0aDPPcbakT9PHN329W7B/aPp870uaLel+SQMKHhuSjis65k/T9rPT7W3TuOZKmiXpRklfqfDad5X0vKT5kh6UtOZSvHVmgBOBdRwvA8ObNiRtAnQt6vPbiOgO9AFeB65ayue4ISK6F3z1TJ9rHeBq4ERgZWAg8Afg84LHvggcUnS8H6TtTXoBo4ABwJrAB8CfygWTJr6xwJlAb2AScMNSviYzJwLrMK4h+WBtcgjJh/MSIuIj4K/A5lV67s2BlyPi/kh8EBE3R8TMgj4TgW6SNgJIv3dN25viGh8RN0bE+xExH7gU2K6Z590fmJI+5mPgbGCzwpGOWUs4EVhH8TjQQ9IGkhqAA4C/lOooaUWS0cO0Kj33v4D1JV0kaWdJ3cv0K0xWZRNVgR2AKc3s3wiY3LQREfOA6Wm7WYs5EVhH0vRBuzvwPMn0T6GT0nn9D4DtgYOX8vjfk/RuwdeDABExA9iJZMrpr8DsdE2iOCH8BRguqRMwjDKJCkDSpsAvgJObiac78F5R23vASkvxmsycCKxDuQY4EDiU0n9tX5DO6w8APgLWK3UQSQcVLAiPL9j114joWfC1c9OOiHg8Ir4XEasB3yD5a/70wuOmU0XTgF8DL0XEa2Wefx1gPHB8RDyStvUvXKhOu34I9Ch6eA+SRGfWYk4E1mFExKski8bfIllELddvJnA8cLGk4gVlIuLaggXhPVsRx8T0+TcusbtpUbnktFB61s99wDkRcU1hzIUL1WnzFGCzgseuCKxN89NJZktwIrCO5nBgl3S+vKyImAC8AYxY1ieUtL2kH0n6Urq9PrAPybpFsRuAb5JMIRUfpw/wAHBZRIxswVPfAmws6b8kdSGZSno6Ip5v5UuxnHIisA4lIqZHxKQWdj8f+LmkFVrY/4Ci6wg+TD/83yX54H8mnba5m+RD+rcl4vsoIu5Lz1wqdgSwFnBWiWmgJUTELOC/gF8B7wDbkKw9mC0V+cY0Zmb55hGBmVnOORGYmeWcE4GZWc5llggkjZb0tqRny+yXpN9LmibpaUlbZBWLmZmVt3yGxx5DUiul3GX0ewKD0q9tgMvT781addVVY8CAAdWJ0MwsJ5588snZ6QWPS8gsEUTEw01leMsYClwdyWlLj0vqKekrEfFmc8cdMGAAkya19OxAMzMDkPRquX21XCPoAxReYt+YtpmZWRvKcmqoEpVoK3lRg6QRpFeA9u/fv9VPuMc5d7b6sWZmtXbPmXtlctxajggagX4F231JLvlfQkSMiojBETF4tdVKTnGZmVkr1XJEMA44RtL1JIvE71VaHzCzfFqx83J8d/NefLlHJ1RyMiEfpk6dWrFPly5d6Nu3L506dWrxcTNLBJKuI6nRvqqkRuAsoBNAWlDrLpIqkdOA+cAPs4rFzNq3727ei40GrsEKK66ElN9EsO4aPZvdHxHMmTOHxsZGBg4c2OLjZnnW0PAK+wM4OqvnN7OO48s9OuU+CbSEJFZZZRVmzZq1VI/zlcVmVveEnARaqDXvkxOBmVnO1XKx2MysVRacu1nlTkthhTMmN7t/vT69+OGIozn1rHMBuGrkJcyfN49jTzyVSy48j24rrsjhRx67qP8u22zKTeMfpHfvVaoaZ1Y8IjAzq6DzCitw7/jbmTt3Tq1DyYQTgZlZBcs3LM8BBx3Cn0f9odahZMJTQ2ZmLXDQoUewz27bc8RPjlti35g/Xs64m7+4DfXbb/2nLUNbZk4EZmYt0H2lHgz9zjCuvmoUXbp2WWzfoT86aok1gvbEU0NmZi10yBFHcfP11/DR/Pm1DqWqnAjMzFqoZ69eDPn2vtx03TW1DqWqPDVkZu1OpdM9s3TYj4/h2j9dWbPnz4ITgZlZBf9+qXHRz6uu9iUmT/+iUPKxJ566RP8Hnni6TeKqFk8NmZnlnBOBmVnOORGYmeWcE4GZWc45EZiZ5ZwTgZlZzvn0UTNrd67r06uqxxv++jsV+1x+8QXccevNLNewHMtpOf7nNxdxwa/P5udnnsMmm30VgMbXZnLkIQdwxwP/qGp8WXMiMDOr4N+T/slD993LLXc/ROcVVmDu3Dl8+skntQ6rapwIzMwqmPX2W/Tq3ZvOK6wA0G5uONNSTgRmZhVst+POXHbRb9lj+8F87Rs78q199mfrr20HwEnHjKBLl6Qa6aeffspyy7W/eys7EZiZVbDiit0Ze/dDTHriHzzx2COccNRhnHjaWQBccOmoJdYI2hsnAjOzFmhoaGCbr2/PNl/fnnXX35Bbb7yu1iFVjU8fNTOrYMa0l3hlxvRF21OnPMMaffvVMKLq8ojAzNqdlpzuWU3z58/j3DN+zvvvv0/D8g2sOWAt/ue3v+P4EYe0aRxZcSIwM6tg40035/px9y7Rfs1Ndyy23bdf/3Z3DQF4asjMLPecCMzMcs6JwMws55wIzMxyzonAzCznnAjMzHLOp4+aWbuz5/9tWdXjjT/wybL73pk7l0MPGArA7Flvs1xDw6KiczfeeT+dO3dm7A3/x7NP/5tf/Or8qsa1yzabctP4BzMvcudEYGbWjF69e3PbhEcAuOTC8+i24oocfuSxi/YvXLiwVqFVTaaJQNIQ4GKgAbgyIs4r2r8y8BegfxrLBRHxpyxjMjNbVqf+9Ces3LMXzz37NBttshnrrr9hyX5nnfoznpn8bxZ8/DF77LUPx510GpD8pb/vd4fz4IS7WbjwU353xRjWXmdd3pk7lxOPPoK5c2az6eZbEBFt8noyWyOQ1ABcBuwJbAgMl1T8bh0NPBcRmwE7ARdK6pxVTGZm1fLKjGmMueFWTj3r3LJ9TjjlTMaOf5Bx9z3KxMf/zvPPPbtoX6/eq3DLPX9j2MGHMXrkpQBcdtFv2GLrbbn13ofZ5Zt78sbrjZm/Dsh2sXhrYFpEzIiIT4DrgaFFfQJYSZKA7sBcoP2Ps8yswxuy9740NDQ022f87bew3x47su8eO/DSC88z/aUXFu375p57A0n5itdfmwnAxMcfY+j+3wNgp932YOWePbMJvkiWU0N9gNcKthuBbYr6XAqMA94AVgIOiIjPM4zJzKwqunbr1uz+12a+yugrLuWmOx9g5Z49OfWnP2HBxwsW7e+U3u1suYYGPvus4O9ftf2NbbIcEZR6NcUTXnsATwFrAJsDl0rqscSBpBGSJkmaNGvWrGrHaWZWdfM+eJ+uXbuxUo8ezJ71Ng8/eF/Fx2y17de5feyNAPztgQm89+67GUeZyHJE0AgUFuzuS/KXf6EfAudFsiIyTdLLwPrAPws7RcQoYBTA4MGD22b1xMzqVnOne9bKLX+9jvvuvmvR9l9vv5cNN96UvXb+Gv36r8kWWxVPiCzp6BNO4cSjj2C/PXZkq223Y40+fbMMeRFltSotaXngRWBX4HVgInBgREwp6HM58FZEnC1pdeBfwGYRMbvccQcPHhyTJk1qVUx7nHNnqx5nZrV12m5fYY0116p1GDW37ho9W9Rv6tSpbLDBBou1SXoyIgaX6p/ZiCAiFko6BriH5PTR0RExRdKR6f6RwDnAGEnPkEwlndJcEjAzs+rL9DqCiLgLuKuobWTBz28A38wyhpaa1vW4ku3rfPT7No7EzKxt5erK4nIf9mZmeeaic2ZmOedEYGaWc04EZmY5l6s1AjPrGI696u9VPd4lh29Xdl9LylAD3Hf3nVx8wa/59JNP6dRpeY4/+XR2G7JXVePMihOBmVkzWlKGetoLU/nNOWcy+rpb6Nd/TV6b+SqHDd+Pvv3XZP0NN65V6C3mRGBmtpSKy1DPmT2LHx/7M/r1XxOAfv3XZMQxJ3DV5Zdw/iVX8OrLMzjr1J8xd85sGhoauPiKMfQfMJArL/8942+/lU8+WcDuQ/ZeVKb6J4cdxH/eeJ0FCxbwg8N/zAHfPxSA7t27c/zxx3PHHXfQtWtXbrvtNlZfffVlfj1eIzAza4XCMtTTXnyejTfZbLH9m2y6OdNefB6Ak44dwUGHHsG4+x7l+tvuYbXVV+fRvz3Aqy/P4KY77+e2ex9hytNPMfHxZMrr1xdeyti7H+Lmux7gmtFX8M7cuQDMmzePbbfdlsmTJ7PDDjvwxz/+sSqvxSMCM7NWKCxDHRFLVA1NmsSHH37AW2++ye5p2ekVunQB4O9/e5C//+0B9v3mDgDMnz+PV16ewVbbbsc1o69gwvg7AHjzjdd59eXp9Ordm86dO7P33slxttxySyZMmFCV1+JEYGbWCoVlqNdZb32effqpxdYDpjw7mbUHrZdkhBIighHHnMCwg3+4WPsTjz3KY488xA2330vXrt04+Dt7s2BBUr66U6dOKE04DQ0NVbtNpqeGKpjW9biSX2ZmTQ7/8bGMuuQiGtMbzDS+NpMrLvlfDvvxMXRfqQdf/soa3Hd3UvTykwUL+Oij+Wy/0y7cfMO1zJv3IQBvvfkGc2bP4oMP3mfllXvStWs3pk97kaf+1boim0vDIwIza3eaO92zFjbYeBNOOv1sjjp0GJ9+upBOnZbn5NN/yQYbbwLAb38/kl+ccgIXX/BrOi3fiYuvGMP2O+7C9JdeZNg+Sbm1bt26c/4lV7DDTrty/TWj+fZu2zFwrUFsvkXJgqFVlVkZ6qwsSxnqtS9Yu2pxuBidWdtxGepEVmWoPTVkZpZzTgRmZjnnRGBmdS8I2ts0dq205n1yIjCzuvef9z9lwbz3nQwqiAjmzJlDl/RahZbyWUOpH5w8o2T71ed7gcqs1m586h0AvtxjNkIVendcn73XtWKfLl260Lfv0t303onAzOrevE8+Z8w/59Q6jJq758xsqpl6asjMLOdyNSIoN/1jZpZnuUoE1VSuzIQvNDOz9sZTQ2ZmOedEYGaWc04EZmY550RgZpZzTgRmZjnnRGBmlnNOBGZmOefrCKrM1xeYWXvjEYGZWc45EZiZ5ZwTgZlZzjkRmJnlXKaLxZKGABcDDcCVEXFeiT47Ab8DOgGzI2LHLGOqFS8im1m9yiwRSGoALgN2BxqBiZLGRcRzBX16An8AhkTETElfyioeMzMrLcsRwdbAtIiYASDpemAo8FxBnwOBsRExEyAi3s4wnlbxLSzNrKPLco2gD/BawXZj2lZoXaCXpIckPSnpBxnGY2ZmJWQ5Iih1h+ko8fxbArsCXYF/SHo8Il5c7EDSCGAEQP/+/TMItXa8dmBmtZbliKAR6Few3Rd4o0SfuyNiXkTMBh4GNis+UESMiojBETF4tdVWyyxgM7M8yjIRTAQGSRooqTMwDBhX1Oc24BuSlpfUDdgGmJphTGZmViSzqaGIWCjpGOAektNHR0fEFElHpvtHRsRUSXcDTwOfk5xi+mxWMbUnnjIys7aS6XUEEXEXcFdR28ii7fOB87OMw8zMymt2aij9i97MzDqwSiOCw4BL2yIQaxlPGZlZtbnWkJlZzlUaEWwq6f0S7QIiInpkEJOZmbWhSongmYj4aptEYsvEU0Zm1lqeGjIzy7lKieDGNonCzMxqptLUUCdJvyizLyLinGoHZNXlKSMzq6RSIviwRFs34AhgFcCJwMysnWs2EUTEhU0/S1oJOJ7k2oLrgQvLPc7MzNqPiiUmJPUGfgYcBPwZ2CIi3sk6sHrX3m9Y4ykjM2vSbCKQdD6wPzAK2CQiSk0VmZlZO1ZpRHAisAA4AzhdWnSvGV9Q1kGVGymARwtmHVWlNQJfZ2Bm1sFVmhrqXdQUwLsRUXzLSTMza6cqTQ09SfLhX3j/4ZUkPQUcERGvZBRXJg558MHMn2OPtbtk/hz1putng2odgpktg0pTQwNLtUvaHxgJDMkiKDMzazutukNZRIyVdEa1g8nan3feOfPnaC+nj1aTF5HN2sY9Z+6VyXFblQgkdccF6yzlaxLM2rdKi8U/K9HcC9gH37nMzKxDqDQiWKloO4D/AN+PiGeyCck6Co8UzNqHSovFv2z6OZ0OioiYl3lUZmbWZlpSa+go4DRgxXT7Q+A3EfGHjGOzDsojBbP6UmmN4Azg68BOETEjbVsLuFhS74g4tw1itJxwgjCrjUojgoOBzSLi46aGiJgh6XvAZMCJoEh7r0paj5qrf1SKE4fZ0qk4NVSYBAraPpL0eTYhmS0bjyzMlk6lRNAoadeIuL+wUdKuwJvZhWVWfR5ZmJVWKREcB9wm6VG+qDu0FbAdMDTj2MxqyiW5LS8qJYIFwKHAusBGJMXnHgauApaYMjIzs/anUiL4HfDfETG6sFHS4HTft7MJy6y+eZrJOpJK9YIGRMTTxY0RMQkYkElEZmbWpiolguaK63etZiBmZlYblaaGJkr6UUT8sbBR0uEki8dm1gI+pdXqWaVE8FPgFkkH8cUH/2CgM7BfhnF1OL7QzMzqVaWic28BX5e0M7Bx2nxnRDyQeWRmOeCRgtWDFt1cJiIejIhL0q8WJwFJQyS9IGmapFOb6beVpM8kfaelxzYzs+po1R3KWkJSA3AZsDvQSLLeMC4inivR7zfAPVnFYtbeeKRgbSmzRABsDUwrqFp6PcnVyM8V9TsWuJnkimUza4YThGUhy/sO9wFeK9huTNsWkdSHZNF5ZHMHkjRC0iRJk2bNmlX1QM3M8izLEYFKtEXR9u+AUyLiM6lU9/RBEaOAUQCDBw8uPka75rOJrBo8UrBlkWUiaAT6FWz3Bd4o6jMYuD5NAqsC35K0MCJuzTAus9xwgrCWyDIRTAQGSRoIvA4MAw4s7BARA5t+ljQGuMNJwMysbWWWCCJioaRjSM4GagBGR8QUSUem+5tdFzAzs7aR5YiAiLgLuKuorWQCiIhDs4zFzL7gKSMrlGkiMLP2xQkin5wI6pTPJrL2wHdx6xiyvI7AzMzaAY8IzKyipb0jm7UvHhGYmeVcrkYE+x/Vv+y+sZfPbMNIWs9rB9ZeeOG5/chVIjCz2nOCqD9OBGZWF5wgaseJoIPwlJGZtZYXi83Mcs6JwMws5zw11MF5ysg6Kq8pVI8TgZlZCXlKNE4EqXLXGLSX6wvMOqqlvao5Tx/g1eJEkFOeMjKzJk4EZpZrrqPkRGBFyo0UwKMFs47KicDMLGP1vm7hRGAt5nUFa888BVSeE4GZWRW050TjRFCBTyutzCMFy5P2/IFfjhOBZcYJwqx9cCKwNucEYVZfnAisbjR36motODFZXjgRtJLXDjo+X1NheeEy1GZmOecRgVkreJ3DOhInArMqcoKwaih/iupemTyfE0GVee3AzNobJwIzsxqpl4vTnAjaiEcK+eYpI6tnPmvIzCznPCKoMY8U8s0jBasHHhGYmeVcpiMCSUOAi4EG4MqIOK9o/0HAKenmh8BRETE5y5jaC48U8s0jBWtLmSUCSQ3AZcDuQCMwUdK4iHiuoNvLwI4R8Y6kPYFRwDZZxWTW3jlBWBayHBFsDUyLiBkAkq4HhgKLEkFEPFbQ/3Ggb4bxdAgeKZhZtWW5RtAHeK1guzFtK+dwYHypHZJGSJokadKsWbOqGKKZmWU5IlCJtijZUdqZJBFsX2p/RIwimTZi8ODBJY+Rdx4p5JunjGxZZJkIGoF+Bdt9gTeKO0naFLgS2DMi5mQYTy45QZhZJVkmgonAIEkDgdeBYcCBhR0k9QfGAgdHxIsZxmKWSx4pWEtklggiYqGkY4B7SE4fHR0RUyQdme4fCfwCWAX4gySAhRExOKuY7AseKeSbE4QVyvQ6goi4C7irqG1kwc9HAEdkGYMtnXIJApwkzDoql5iwFvMoouPzSCGfnAhsmTlBmLVvTgSWmeammUpx4jCrDScCqxtLmzhqKW9Jq9yUUXM8ndR+uPqomVnOeURg1gqe9qrMC8/thxOBWRvwgrrVMycCsxrKY4LwSKH+OBGYtSN5TByWPScCszrUns6gsvbPicDM6oKnjGrHicCsA8jjlNHSXtvghFKeE4FZB9YREkRrLmazpeNEYJZDHSFBZC1PU1VOBGa2SB4ThEccTgRmtgza0/0r/IFfnhOBmVXk01mXTb1PMzkRmFkm8jjN1F45EZhZm2rvCaIjTjE5EZhZXahWgmjviaYWnAjMrK7l8YO97KjjpGyez4nAzNqljrCAXS/TTL5DmZlZznlEYGa54LvKledEYGZWQp7WJpwIzMyqoDVrFvWSVJwIzMyWQlssUrf1QrgXi83Mcs4jAjOzGqmXU2A9IjAzyzknAjOznHMiMDPLOScCM7OccyIwM8u5TBOBpCGSXpA0TdKpJfZL0u/T/U9L2iLLeMzMbEmZJQJJDcBlwJ7AhsBwSRsWddsTGJR+jQAuzyoeMzMrLcsRwdbAtIiYERGfANcDQ4v6DAWujsTjQE9JX8kwJjMzK5LlBWV9gNcKthuBbVrQpw/wZmEnSSNIRgwAH0p6oZUxrQrMbuVj24pjXHb1Hh/Uf4z1Hh/Uf4zVj+9yLcuj1yy3I8tEUCriaEUfImIUMGqZA5ImRcTgZT1Olhzjsqv3+KD+Y6z3+KD+Y6z3+AplOTXUCPQr2O4LvNGKPmZmlqEsE8FEYJCkgZI6A8OAcUV9xgE/SM8e2hZ4LyLeLD6QmZllJ7OpoYhYKOkY4B6gARgdEVMkHZnuHwncBXwLmAbMB36YVTypZZ5eagOOcdnVe3xQ/zHWe3xQ/zHWe3yLKGKJKXkzM8sRX1lsZpZzTgRmZjmXm0RQqdxFDeLpJ+lBSVMlTZF0fNreW9IESS+l33vVQawNkv4t6Y56jFFST0k3SXo+fT+/Vk8xSjoh/Td+VtJ1krrUOj5JoyW9LenZgrayMUk6Lf3deUHSHjWK7/z03/hpSbdI6lmr+MrFWLDvJEkhadVaxthSuUgELSx30dYWAidGxAbAtsDRaUynAvdHxCDg/nS71o4HphZs11uMFwN3R8T6wGYksdZFjJL6AMcBgyNiY5ITJ4bVQXxjgCFFbSVjSv9fDgM2Sh/zh/R3qq3jmwBsHBGbAi8Cp9UwvnIxIqkfsDsws6CtVjG2SC4SAS0rd9GmIuLNiPhX+vMHJB9efdK4/px2+zOwb00CTEnqC+wFXFnQXDcxSuoB7ABcBRARn0TEu9RRjCRn53WVtDzQjeRamZrGFxEPA3OLmsvFNBS4PiIWRMTLJGf5bd3W8UXEvRGxMN18nOS6o5rEVy7G1EXAz1n84tiaxNhSeUkE5UpZ1AVJA4CvAk8AqzddS5F+/1INQwP4Hcl/6s8L2uopxrWAWcCf0umrKyWtWC8xRsTrwAUkfx2+SXKtzL31El+RcjHV4+/PYcD49Oe6iU/SPsDrETG5aFfdxFhKXhJBi0pZ1IKk7sDNwE8j4v1ax1NI0t7A2xHxZK1jacbywBbA5RHxVWAetZ+qWiSdZx8KDATWAFaU9P3aRrXU6ur3R9LpJFOr1zY1lejW5vFJ6gacDvyi1O4SbXXxGQT5SQR1WcpCUieSJHBtRIxNm99qqsCafn+7VvEB2wH7SHqFZDptF0l/ob5ibAQaI+KJdPsmksRQLzHuBrwcEbMi4lNgLPD1OoqvULmY6ub3R9IhwN7AQfHFRVD1Et/aJAl/cvo70xf4l6QvUz8xlpSXRNCSchdtSpJI5rWnRsT/FuwaBxyS/nwIcFtbx9YkIk6LiL4RMYDkPXsgIr5PfcX4H+A1SeulTbsCz1E/Mc4EtpXULf0335VkPahe4itULqZxwDBJK0gaSHL/kH+2dXCShgCnAPtExPyCXXURX0Q8ExFfiogB6e9MI7BF+n+0LmIsKyJy8UVSyuJFYDpweh3Esz3J0PBp4Kn061vAKiRnbLyUfu9d61jTeHcC7kh/rqsYgc2BSel7eSvQq55iBH4JPA88C1wDrFDr+IDrSNYsPiX5wDq8uZhIpjymAy8Ae9Yovmkk8+xNvy8jaxVfuRiL9r8CrFrLGFv65RITZmY5l5epITMzK8OJwMws55wIzMxyzonAzCznnAjMzHLOicCaJekzSU+llTNvTK+eLG6/vakSpKQBkj5K9zV9/SDd113SFZKmp9U4H5a0TbrvwxKPf07SSEnLpe3PStqj4LgfppUcn5J0dfr4/dKqj+un20+k+2dKmlXw2AGSXmmqDimpr6TblFTenC7p4vSaEyTtlB7z2wXvyx2Sdirxfo2R9HL6HJMl7Vq0/wRJH0taWdIqBfH8R9LrBdudC97jpq8lrpiW9FD6HuyzzP/YJUjqJGmJK8slnS3ppFYe8/z09bbq8VZ9md2q0jqMjyJicwBJ1wJHAv9b1P5n4GjgV+ljpjftK3Il8DIwKCI+l7QWsEGJftMjYnMlRdoeICl+1lSg7x6S258i6SHgpIiYVPDY4cCjJBfAnR0RTYnmUJIKoMc0dUyu71p0cd9YkjIVQ5VUhRyVvp6T0+6NJOeB397cm5U6OSJukrRzepxBRfFNBPaLiDEk10Ag6Wzgw4i4oCC+j8q8j8UOKnoPqml74LFqHjAiTpY0r5rHtGXjEYEtjUeAdUq0/4MKBbQkrQ1sA5wREZ8DRFIN9s5yj4mk0uRjZZ6z1HN0JymLcThJImipXYCPI+JP6fN+BpwAHNY0AgImA+9J2n0pjrvY+5K+B92BM0gSQtWlI4SL0tHWVElbSRqbjnTOTfsMUFLX/8p0lHWtpN0k/T3tV1gVcwhpcTdJp6ejj/uA9Qqe80eSJqYjoJuVXEW9Ujoy6pT26ZGOwDpl8bpt2TgRWIukf53vCTxT1N5AUjahsGTH2kVTGt8gqcP+VPoh29Ln7JYe+5lKfVP7ktyX4EVgrqQtWvi4jYDFpj8iKQA4k8WT0LkkH+ItNYTkSucmw0muRn0EWE9SpYqjXYvexwNa+LyfRMQOwEiSMhFHAxsDh0paJe2zDsl9HDYF1gcOJPnr/yTgvwuOtTPwkKQtSZLrV4H9ga0K+oyNiK0iouleEIdHUlr9IZIS5qSPvTmSektWZzw1ZJV0lfRU+vMjpHX/C9oHkHyITih4zBJTQ0s5h712euwAbouI8UpKdVcynKRsNiRF8oaTTilVIEpXglysPSIekUSa2JpzvqTfkpRx3ragfRjJlNDnksYC3yW5YVI5LZ0aKtaUlJ8BpkRaWlrSDJLCZ++SFMJ7Jm2fQnJDmpD0DMm/KZLWAOZGxPz0Nd8SaY0fSYWJf+N0tNGTZMRzT9p+JUkJ81uBHwI/asVrsTbgRGCVlPsw+iidx18ZuIPkr87fN3OcKcBmkpZrmhpqRrk1hrLSv3R3IflQCpI7gYWkn0flOipTgP8qOl4Pkg/N6SQ1eJr8imStYCHlnUyy5nAcyQ1etpS0KclawYR0baIzMIPmE0FrLUi/f17wc9P28kV9ivsV9tmTLz7UoXzZ5DHAvhExOV2L2QkgIv6eTkPtCDRExBK3dLT64KkhWyYR8R7JB95Jzc3/RsR0ksJwv0wXZ5E0SFK17hT3HeDqiFgzkuqP/UgWprdvwWPvB7rpi7ObGoALgTGxeJVLIrmpTC+SW2KWlSa7i4HllNyfdjjJ4vWA9GsNoI+kNZfuZbapResDwMPAfpK6SloJ+HZBv5WAN9N//4OKjnE1yXTYn7IO1lrPicCWWUT8m2QxtWmBtniN4Li0/Qjgy8C0dArij1SvJvtw4JaitptJ5r6blY4Y9gO+K+klkiq1H7P4XHmhX/HFbRIrHfdckumRYSXiu4XmF7WL1wjOq/Sc1ZImw0ER8TxAJLdVvYGk6ufNJNOETc4kubveBJIqq4WuJUmc12Ucsi0DVx81a8fKnEJbjeNuD3w/Io5cxuN8BxgaEQcXtZ9N0emyVjseEZi1b3OBMUu5GF9RRDxahSRwCXAecE5R+/nA90luK2p1wCMCM7Oc84jAzCznnAjMzHLOicDMLOecCMzMcs6JwMws5/4fuAG/wBhwwLMAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# alternatively, show each bin's fractional contributions\n",
    "\n",
    "bctr = prbins[0:-1] + 0.5*(prbins[1:]-prbins[0:-1])\n",
    "\n",
    "tot = h_nh[0]+h_sh[0]+h_trop_land[0]+h_trop_ocean[0]\n",
    "\n",
    "width = 3\n",
    "fig, ax = plt.subplots()\n",
    "ax.bar(bctr, h_nh[0]/tot, color='chocolate', width=width, label='NH')\n",
    "ax.bar(bctr, h_sh[0]/tot, color='maroon', bottom=h_nh[0]/tot, width=width, label='SH')\n",
    "ax.bar(bctr, h_trop_land[0]/tot, color='forestgreen', bottom=(h_nh[0]+h_sh[0])/tot, width=width, label='TrLand')\n",
    "ax.bar(bctr, h_trop_ocean[0]/tot, color='steelblue', bottom=(h_nh[0]+h_sh[0]+h_trop_land[0])/tot, width=width, label='TrOcean')\n",
    "# ax.semilogy()\n",
    "ax.set_ylim([0,1.01])\n",
    "ax.axhline(.5, color='w')\n",
    "ax.legend()\n",
    "ax.set_ylabel('COUNT')\n",
    "ax.set_xlabel('PRECIPITATION RATE [mm/day]')\n",
    "ax.set_title(model_name)\n",
    "fig.savefig(f\"/Users/brianpm/Dropbox/Manuscripts/M_cre_pex/REVISION1/pr_histogram_fractional_{model_name}_{experiment_name}.png\", dpi=144)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "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.1"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}