pdf_regcm_miroc5.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "b2254f18",
   "metadata": {},
   "outputs": [],
   "source": [
    "#writing some conda package or module\n",
    "#!pip install netCDF4\n",
    "%matplotlib inline\n",
    "import warnings\n",
    "from netCDF4 import Dataset\n",
    "import numpy as np\n",
    "import xarray as xr\n",
    "import pandas as pd\n",
    "import matplotlib.pyplot as plt\n",
    "#!pip install cartopy\n",
    "import cartopy.feature as cfeature\n",
    "import cartopy.crs as ccrs\n",
    "import cartopy.mpl.ticker as cticker\n",
    "from cartopy.util import add_cyclic_point"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "9b4e84e2",
   "metadata": {},
   "outputs": [],
   "source": [
    "rain_minimum = 1.0\n",
    "rain_maximum = 250.0\n",
    "rain_step = 1.0\n",
    "spd = 86400.0\n",
    "sph = 3600.0\n",
    "mpy = 30.4167\n",
    "\n",
    "nrainbin = int((rain_maximum-rain_minimum)/rain_step) + 1\n",
    "\n",
    "rain_pdf = np.linspace(rain_minimum,rain_maximum, num=nrainbin, endpoint=True)\n",
    "rain_bin = np.zeros((nrainbin-1))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "3e1b1bb0",
   "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 !important;\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.Dataset&gt;\n",
       "Dimensions:    (time: 432, bnds: 2, y: 337, x: 429)\n",
       "Coordinates:\n",
       "  * time       (time) object 1970-01-16 12:00:00 ... 2005-12-16 12:00:00\n",
       "    lon        (y, x) float64 ...\n",
       "    lat        (y, x) float64 ...\n",
       "  * x          (x) float64 -5.35e+06 -5.325e+06 -5.3e+06 ... 5.325e+06 5.35e+06\n",
       "  * y          (y) float64 -4.2e+06 -4.175e+06 -4.15e+06 ... 4.175e+06 4.2e+06\n",
       "Dimensions without coordinates: bnds\n",
       "Data variables:\n",
       "    time_bnds  (time, bnds) object 1970-01-01 00:00:00 ... 2006-01-01 00:00:00\n",
       "    crs        |S1 b&#x27;&#x27;\n",
       "    pr         (time, y, x) float32 ...\n",
       "Attributes: (12/165)\n",
       "    CDI:                                             Climate Data Interface v...\n",
       "    Conventions:                                     CF-1.7\n",
       "    history:                                         Mon Jun 27 15:34:28 2022...\n",
       "    source:                                          RegCM Model output file\n",
       "    institution:                                     Oak Ridge National Labor...\n",
       "    project_id:                                      CORDEX\n",
       "    ...                                              ...\n",
       "    experiment_id:                                   historical\n",
       "    driving_experiment:                              MIROC-MIROC5, historical...\n",
       "    driving_model_id:                                MIROC-MIROC5\n",
       "    driving_experiment_name:                         historical\n",
       "    model_id:                                        ICTP-RegCM4-7\n",
       "    CDO:                                             Climate Data Operators v...</pre><div class='xr-wrap' style='display:none'><div class='xr-header'><div class='xr-obj-type'>xarray.Dataset</div></div><ul class='xr-sections'><li class='xr-section-item'><input id='section-8ca5d43e-9edf-4aad-bb90-af752b8af609' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-8ca5d43e-9edf-4aad-bb90-af752b8af609' class='xr-section-summary'  title='Expand/collapse section'>Dimensions:</label><div class='xr-section-inline-details'><ul class='xr-dim-list'><li><span class='xr-has-index'>time</span>: 432</li><li><span>bnds</span>: 2</li><li><span class='xr-has-index'>y</span>: 337</li><li><span class='xr-has-index'>x</span>: 429</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-ea355394-4d78-4fa6-abf5-fe29af47000a' class='xr-section-summary-in' type='checkbox'  checked><label for='section-ea355394-4d78-4fa6-abf5-fe29af47000a' class='xr-section-summary' >Coordinates: <span>(5)</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'>object</div><div class='xr-var-preview xr-preview'>1970-01-16 12:00:00 ... 2005-12-...</div><input id='attrs-68db560f-d9a5-4f7a-83f4-56cb5a4288ef' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-68db560f-d9a5-4f7a-83f4-56cb5a4288ef' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-67cc8b69-a8e5-4f19-97c3-1d21d63a790f' class='xr-var-data-in' type='checkbox'><label for='data-67cc8b69-a8e5-4f19-97c3-1d21d63a790f' 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>standard_name :</span></dt><dd>time</dd><dt><span>long_name :</span></dt><dd>time</dd><dt><span>bounds :</span></dt><dd>time_bnds</dd><dt><span>axis :</span></dt><dd>T</dd></dl></div><div class='xr-var-data'><pre>array([cftime.DatetimeNoLeap(1970, 1, 16, 12, 0, 0, 0, has_year_zero=True),\n",
       "       cftime.DatetimeNoLeap(1970, 2, 15, 0, 0, 0, 0, has_year_zero=True),\n",
       "       cftime.DatetimeNoLeap(1970, 3, 16, 12, 0, 0, 0, has_year_zero=True),\n",
       "       ...,\n",
       "       cftime.DatetimeNoLeap(2005, 10, 16, 12, 0, 0, 0, has_year_zero=True),\n",
       "       cftime.DatetimeNoLeap(2005, 11, 16, 0, 0, 0, 0, has_year_zero=True),\n",
       "       cftime.DatetimeNoLeap(2005, 12, 16, 12, 0, 0, 0, has_year_zero=True)],\n",
       "      dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lon</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4a459c7e-f5ff-4de6-9fb1-e6ff3c523e9b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4a459c7e-f5ff-4de6-9fb1-e6ff3c523e9b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c0adbb1b-e8c4-41a3-9463-1df474e0ceb9' class='xr-var-data-in' type='checkbox'><label for='data-c0adbb1b-e8c4-41a3-9463-1df474e0ceb9' 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>standard_name :</span></dt><dd>longitude</dd><dt><span>long_name :</span></dt><dd>Longitude on Cross Points</dd><dt><span>units :</span></dt><dd>degrees_east</dd><dt><span>_CoordinateAxisType :</span></dt><dd>Lon</dd></dl></div><div class='xr-var-data'><pre>[144573 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>lat</span></div><div class='xr-var-dims'>(y, x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-90576e18-cf64-497a-a6c3-96ed3d883cb7' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-90576e18-cf64-497a-a6c3-96ed3d883cb7' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-62d9dd18-70c6-4a64-b6e2-d81866060e49' class='xr-var-data-in' type='checkbox'><label for='data-62d9dd18-70c6-4a64-b6e2-d81866060e49' 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>standard_name :</span></dt><dd>latitude</dd><dt><span>long_name :</span></dt><dd>Latitude on Cross Points</dd><dt><span>units :</span></dt><dd>degrees_north</dd><dt><span>_CoordinateAxisType :</span></dt><dd>Lat</dd></dl></div><div class='xr-var-data'><pre>[144573 values with dtype=float64]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>x</span></div><div class='xr-var-dims'>(x)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-5.35e+06 -5.325e+06 ... 5.35e+06</div><input id='attrs-263efd44-605f-4b4d-af97-ae131a5972cf' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-263efd44-605f-4b4d-af97-ae131a5972cf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c040cf31-87a0-47fc-bbc7-ef9718914af2' class='xr-var-data-in' type='checkbox'><label for='data-c040cf31-87a0-47fc-bbc7-ef9718914af2' 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>standard_name :</span></dt><dd>projection_x_coordinate</dd><dt><span>long_name :</span></dt><dd>x-coordinate in Cartesian system</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>axis :</span></dt><dd>X</dd></dl></div><div class='xr-var-data'><pre>array([-5350000., -5325000., -5300000., ...,  5300000.,  5325000.,  5350000.])</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>y</span></div><div class='xr-var-dims'>(y)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>-4.2e+06 -4.175e+06 ... 4.2e+06</div><input id='attrs-d8afe54a-2177-44d7-a09f-50a1d036066a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d8afe54a-2177-44d7-a09f-50a1d036066a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-b2ce507f-1ce0-4fbc-9c43-367dc9ac4ac0' class='xr-var-data-in' type='checkbox'><label for='data-b2ce507f-1ce0-4fbc-9c43-367dc9ac4ac0' 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>standard_name :</span></dt><dd>projection_y_coordinate</dd><dt><span>long_name :</span></dt><dd>y-coordinate in Cartesian system</dd><dt><span>units :</span></dt><dd>m</dd><dt><span>axis :</span></dt><dd>Y</dd></dl></div><div class='xr-var-data'><pre>array([-4200000., -4175000., -4150000., ...,  4150000.,  4175000.,  4200000.])</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-e13d5911-1d34-48bb-9a80-0209e23bd234' class='xr-section-summary-in' type='checkbox'  checked><label for='section-e13d5911-1d34-48bb-9a80-0209e23bd234' class='xr-section-summary' >Data variables: <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>time_bnds</span></div><div class='xr-var-dims'>(time, bnds)</div><div class='xr-var-dtype'>object</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-9828e24a-0fb6-4e8a-ae0e-d93aba53021a' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-9828e24a-0fb6-4e8a-ae0e-d93aba53021a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1cada5de-95c1-4ee1-a99c-f2e3743936b3' class='xr-var-data-in' type='checkbox'><label for='data-1cada5de-95c1-4ee1-a99c-f2e3743936b3' 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'></dl></div><div class='xr-var-data'><pre>array([[cftime.DatetimeNoLeap(1970, 1, 1, 0, 0, 0, 0, has_year_zero=True),\n",
       "        cftime.DatetimeNoLeap(1970, 2, 1, 0, 0, 0, 0, has_year_zero=True)],\n",
       "       [cftime.DatetimeNoLeap(1970, 2, 1, 0, 0, 0, 0, has_year_zero=True),\n",
       "        cftime.DatetimeNoLeap(1970, 3, 1, 0, 0, 0, 0, has_year_zero=True)],\n",
       "       [cftime.DatetimeNoLeap(1970, 3, 1, 0, 0, 0, 0, has_year_zero=True),\n",
       "        cftime.DatetimeNoLeap(1970, 4, 1, 0, 0, 0, 0, has_year_zero=True)],\n",
       "       ...,\n",
       "       [cftime.DatetimeNoLeap(2005, 10, 1, 0, 0, 0, 0, has_year_zero=True),\n",
       "        cftime.DatetimeNoLeap(2005, 11, 1, 0, 0, 0, 0, has_year_zero=True)],\n",
       "       [cftime.DatetimeNoLeap(2005, 11, 1, 0, 0, 0, 0, has_year_zero=True),\n",
       "        cftime.DatetimeNoLeap(2005, 12, 1, 0, 0, 0, 0, has_year_zero=True)],\n",
       "       [cftime.DatetimeNoLeap(2005, 12, 1, 0, 0, 0, 0, has_year_zero=True),\n",
       "        cftime.DatetimeNoLeap(2006, 1, 1, 0, 0, 0, 0, has_year_zero=True)]],\n",
       "      dtype=object)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>crs</span></div><div class='xr-var-dims'>()</div><div class='xr-var-dtype'>|S1</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-56473792-8600-48e7-af77-17ad08b27cf8' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-56473792-8600-48e7-af77-17ad08b27cf8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-643374e0-e791-4e5d-88f1-9b6caeeb2de3' class='xr-var-data-in' type='checkbox'><label for='data-643374e0-e791-4e5d-88f1-9b6caeeb2de3' 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>proj4_params :</span></dt><dd>+proj=omerc +lat_0=13.00 +alpha=89.999999 +lonc=70.00 +x_0=-12500. +y_0=-12500. +ellps=sphere +a=6371229. +b=6371229. +units=m +no_defs</dd><dt><span>semi_major_axis :</span></dt><dd>6371229.0</dd><dt><span>inverse_flattening :</span></dt><dd>0.0</dd><dt><span>azimuth_of_central_line :</span></dt><dd>89.999999</dd><dt><span>latitude_of_projection_origin :</span></dt><dd>13.0</dd><dt><span>longitude_of_projection_origin :</span></dt><dd>70.0</dd><dt><span>scale_factor_at_projection_origin :</span></dt><dd>1.0</dd><dt><span>false_easting :</span></dt><dd>-12500.0</dd><dt><span>false_northing :</span></dt><dd>-12500.0</dd><dt><span>grid_mapping_name :</span></dt><dd>oblique_mercator</dd></dl></div><div class='xr-var-data'><pre>array(b&#x27;&#x27;, dtype=&#x27;|S1&#x27;)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>pr</span></div><div class='xr-var-dims'>(time, y, x)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-4d232a78-87cf-4dfe-ad77-14609e5c2370' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4d232a78-87cf-4dfe-ad77-14609e5c2370' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0a62a612-f3af-4b29-94f9-e51376e66f73' class='xr-var-data-in' type='checkbox'><label for='data-0a62a612-f3af-4b29-94f9-e51376e66f73' 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>standard_name :</span></dt><dd>precipitation_flux</dd><dt><span>long_name :</span></dt><dd>Precipitation</dd><dt><span>units :</span></dt><dd>kg m-2 s-1</dd><dt><span>grid_mapping :</span></dt><dd>crs</dd><dt><span>cell_methods :</span></dt><dd>time: mean</dd></dl></div><div class='xr-var-data'><pre>[62455536 values with dtype=float32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-4895f1a7-32ce-4d71-99bd-dbaea2a4da60' class='xr-section-summary-in' type='checkbox'  ><label for='section-4895f1a7-32ce-4d71-99bd-dbaea2a4da60' class='xr-section-summary' >Attributes: <span>(165)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>CDI :</span></dt><dd>Climate Data Interface version 1.9.9rc1 (https://mpimet.mpg.de/cdi)</dd><dt><span>Conventions :</span></dt><dd>CF-1.7</dd><dt><span>history :</span></dt><dd>Mon Jun 27 15:34:28 2022: cdo mergetime pr_WAS-22_MIROC-MIROC5_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_197001-198012.nc pr_WAS-22_MIROC-MIROC5_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_198101-199012.nc pr_WAS-22_MIROC-MIROC5_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_199101-200012.nc pr_WAS-22_MIROC-MIROC5_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_200101-200512.nc miroc5_pr_197001-200512.nc\n",
       "2018-05-11 16:23:33 : Created by RegCM RegCM Model program</dd><dt><span>source :</span></dt><dd>RegCM Model output file</dd><dt><span>institution :</span></dt><dd>Oak Ridge National Laboratory</dd><dt><span>project_id :</span></dt><dd>CORDEX</dd><dt><span>ipcc_scenario_code :</span></dt><dd>HISTORICAL</dd><dt><span>institute_id :</span></dt><dd>ORNL</dd><dt><span>note :</span></dt><dd>The domain is larger than WAS-22</dd><dt><span>comment :</span></dt><dd>RegCM CORDEX SAsia25 run</dd><dt><span>experiment :</span></dt><dd>WAS-22</dd><dt><span>driving_model_ensemble_member :</span></dt><dd>r1i1p1</dd><dt><span>creation_date :</span></dt><dd>2019-09-12T16:23:22Z</dd><dt><span>CORDEX_domain :</span></dt><dd>WAS-22</dd><dt><span>rcm_version_id :</span></dt><dd>v0</dd><dt><span>ICTP_version_note :</span></dt><dd>Model run at ORNL on Titan</dd><dt><span>contact :</span></dt><dd>mashfaq@ornl.gov</dd><dt><span>product :</span></dt><dd>output</dd><dt><span>tracking_id :</span></dt><dd>635194b8-8309-11e9-827b-f04da273f023</dd><dt><span>title :</span></dt><dd>ICTP Regional Climatic model V4</dd><dt><span>references :</span></dt><dd>http://gforge.ictp.it/gf/project/regcm</dd><dt><span>model_revision :</span></dt><dd>tag-4.7.0</dd><dt><span>projection :</span></dt><dd>ROTMER</dd><dt><span>grid_size_in_meters :</span></dt><dd>25000.0</dd><dt><span>latitude_of_projection_origin :</span></dt><dd>13.0</dd><dt><span>longitude_of_projection_origin :</span></dt><dd>70.0</dd><dt><span>grid_north_pole_latitude :</span></dt><dd>13.0</dd><dt><span>grid_north_pole_longitude :</span></dt><dd>70.0</dd><dt><span>grid_factor :</span></dt><dd>0.0</dd><dt><span>landsurface_model :</span></dt><dd>clm4.5</dd><dt><span>model_icbc_data_source :</span></dt><dd>MI_85</dd><dt><span>model_sst_data_source :</span></dt><dd>MI_85</dd><dt><span>boundary_nspgx :</span></dt><dd>30</dd><dt><span>boundary_nspgd :</span></dt><dd>30</dd><dt><span>boundary_high_nudge :</span></dt><dd>8.0</dd><dt><span>boundary_medium_nudge :</span></dt><dd>6.0</dd><dt><span>boundary_low_nudge :</span></dt><dd>4.0</dd><dt><span>boundary_nm :</span></dt><dd>-1.0</dd><dt><span>boundary_dm :</span></dt><dd>-1.0</dd><dt><span>model_is_restarted :</span></dt><dd>No</dd><dt><span>model_simulation_initial_start :</span></dt><dd> 1970-01-01 00:00:00 UTC</dd><dt><span>model_simulation_start :</span></dt><dd> 1970-01-01 00:00:00 UTC</dd><dt><span>model_simulation_end :</span></dt><dd> 1975-01-01 00:00:00 UTC</dd><dt><span>atmosphere_time_step_in_seconds :</span></dt><dd>40.0</dd><dt><span>surface_interaction_time_step_in_seconds :</span></dt><dd>480.0</dd><dt><span>convection_time_step_in_seconds :</span></dt><dd>240.0</dd><dt><span>radiation_scheme_time_step_in_seconds :</span></dt><dd>1800.0</dd><dt><span>absorption_emission_time_step_in_seconds :</span></dt><dd>64800.0</dd><dt><span>dynamical_core :</span></dt><dd>1</dd><dt><span>asselin_filter_nu_1 :</span></dt><dd>0.0625</dd><dt><span>asselin_filter_nu_2 :</span></dt><dd>0.0625</dd><dt><span>diffusion_hgt_factor :</span></dt><dd>1</dd><dt><span>upstream_mode :</span></dt><dd>Yes</dd><dt><span>off_centering_max :</span></dt><dd>0.25</dd><dt><span>stability_enhance :</span></dt><dd>Yes</dd><dt><span>temperature_extreme_gradient :</span></dt><dd>5.0</dd><dt><span>cloud_extreme_gradient_fraction :</span></dt><dd>0.5</dd><dt><span>vapor_extreme_gradient_fraction :</span></dt><dd>0.2</dd><dt><span>tracer_extreme_gradient_fraction :</span></dt><dd>0.5</dd><dt><span>lateral_boundary_condition_scheme :</span></dt><dd>5</dd><dt><span>semi_lagrangian_advection_scheme :</span></dt><dd>0</dd><dt><span>boundary_layer_scheme :</span></dt><dd>2</dd><dt><span>cumulus_convection_scheme_lnd :</span></dt><dd>4</dd><dt><span>cumulus_convection_scheme_ocn :</span></dt><dd>5</dd><dt><span>moisture_scheme :</span></dt><dd>1</dd><dt><span>ocean_flux_scheme :</span></dt><dd>2</dd><dt><span>zeng_ocean_roughness_formula :</span></dt><dd>1</dd><dt><span>zeng_ocean_roughness_method :</span></dt><dd>1</dd><dt><span>pressure_gradient_scheme :</span></dt><dd>0</dd><dt><span>surface_emissivity_factor_computed :</span></dt><dd>0</dd><dt><span>lake_model_activated :</span></dt><dd>0</dd><dt><span>chemical_aerosol_scheme_activated :</span></dt><dd>0</dd><dt><span>diurnal_cycle_sst_scheme :</span></dt><dd>0</dd><dt><span>ocean_precipitaion_cooling_effect :</span></dt><dd>0</dd><dt><span>ocean_albedo_from_whitecapping :</span></dt><dd>0</dd><dt><span>simple_sea_ice_scheme :</span></dt><dd>0</dd><dt><span>seasonal_desert_albedo :</span></dt><dd>0</dd><dt><span>convective_lwp_as_large_scale :</span></dt><dd>0</dd><dt><span>large_scale_cloud_fraction_scheme :</span></dt><dd>0</dd><dt><span>ocean_stratiform_clouds_scheme :</span></dt><dd>1</dd><dt><span>rrtm_radiation_scheme_activated :</span></dt><dd>0</dd><dt><span>climatic_ozone_input_dataset :</span></dt><dd>0</dd><dt><span>climatic_aerosol_input_dataset :</span></dt><dd>0</dd><dt><span>static_solar_constant_used :</span></dt><dd>0</dd><dt><span>cumulus_cloud_model :</span></dt><dd>1</dd><dt><span>sun_year_offset :</span></dt><dd>0.0</dd><dt><span>fix_solar_parameters :</span></dt><dd>0</dd><dt><span>cloud_bottom_level_with_no_clouds :</span></dt><dd>1</dd><dt><span>cloud_rh_with_fcc_one :</span></dt><dd>1.01</dd><dt><span>cloud_rh_threshold_for_land :</span></dt><dd>0.8</dd><dt><span>cloud_rh_threshold_for_ocean :</span></dt><dd>0.9</dd><dt><span>cloud_rh1_limit_temperature :</span></dt><dd>238.0</dd><dt><span>cloud_cloud_fraction_maximum :</span></dt><dd>0.75</dd><dt><span>cloud_condensation_threshold :</span></dt><dd>1.0</dd><dt><span>cloud_cloud_fraction_max_for_convection :</span></dt><dd>1.0</dd><dt><span>cloud_cloud_liqwat_max_for_convection :</span></dt><dd>0.0003</dd><dt><span>cloud_shallow_convective_cf_param :</span></dt><dd>0.07</dd><dt><span>cloud_deep_convective_cf_param :</span></dt><dd>0.14</dd><dt><span>cloud_k2_convective_cf_param :</span></dt><dd>500.0</dd><dt><span>cloud_surface_radiation_hack :</span></dt><dd>No</dd><dt><span>cloud_arctic_cloud_correction :</span></dt><dd>Yes</dd><dt><span>subex_auto_conversion_rate_for_land :</span></dt><dd>0.00025</dd><dt><span>subex_auto_conversion_rate_for_ocean :</span></dt><dd>0.00025</dd><dt><span>subex_gultepe_factor_when_rain_for_land :</span></dt><dd>0.4</dd><dt><span>subex_gultepe_factor_when_rain_for_ocean :</span></dt><dd>0.4</dd><dt><span>subex_land_raindrop_evaporation_rate :</span></dt><dd>0.0002</dd><dt><span>subex_ocean_raindrop_evaporation_rate :</span></dt><dd>0.0001</dd><dt><span>subex_land_raindrop_accretion_rate :</span></dt><dd>3.0</dd><dt><span>subex_ocean_raindrop_accretion_rate :</span></dt><dd>3.0</dd><dt><span>mit_lowest_convection_sigma :</span></dt><dd>0.95</dd><dt><span>mit_autoconversion_threshold_mixing_over_ocean :</span></dt><dd>0.00011</dd><dt><span>mit_autoconversion_threshold_mixing_over_land :</span></dt><dd>0.00011</dd><dt><span>mit_autoconversion_threshold_temperature :</span></dt><dd>-55.0</dd><dt><span>mit_mixing_coefficient_in_entrainment :</span></dt><dd>1.5</dd><dt><span>mit_fractional_area_unsaturated_downdraft :</span></dt><dd>0.05</dd><dt><span>mit_fractional_precip_outside_cloud :</span></dt><dd>0.12</dd><dt><span>mit_pressure_velocity_of_rain :</span></dt><dd>50.0</dd><dt><span>mit_pressure_velocity_of_snow :</span></dt><dd>5.5</dd><dt><span>mit_rain_evaporation_coefficient :</span></dt><dd>1.0</dd><dt><span>mit_snow_evaporation_coefficient :</span></dt><dd>0.8</dd><dt><span>mit_momentum_transport_coefficient :</span></dt><dd>0.7</dd><dt><span>mit_downdraft_velocity_coefficient :</span></dt><dd>10.0</dd><dt><span>mit_max_parcel_neg_temp_perturbation :</span></dt><dd>0.9</dd><dt><span>mit_approach_rate_quasi_eq_coeff_a :</span></dt><dd>0.2</dd><dt><span>mit_approach_rate_quasi_eq_coeff_d :</span></dt><dd>0.1</dd><dt><span>mit_maximum_land_precipitation_efficiency :</span></dt><dd>0.999</dd><dt><span>mit_maximum_ocean_precipitation_efficiency :</span></dt><dd>0.999</dd><dt><span>tiedtke_actual_scheme :</span></dt><dd>4</dd><dt><span>tiedtke_entrainment_rate_downdraft :</span></dt><dd>0.0003</dd><dt><span>tiedtke_entrainment_rate_deep_land :</span></dt><dd>0.00075</dd><dt><span>tiedtke_entrainment_rate_deep_ocean :</span></dt><dd>0.00075</dd><dt><span>tiedtke_penetrative :</span></dt><dd>Yes</dd><dt><span>tiedtke_midlevel :</span></dt><dd>Yes</dd><dt><span>tiedtke_shallow :</span></dt><dd>Yes</dd><dt><span>tiedtke_cumulus_downdraft :</span></dt><dd>Yes</dd><dt><span>tiedtke_prognostic_cloud :</span></dt><dd>Yes</dd><dt><span>tiedtke_cumulus_friction :</span></dt><dd>Yes</dd><dt><span>tiedtke_ke_dissipation :</span></dt><dd>Yes</dd><dt><span>tiedtke_tracer_transport :</span></dt><dd>Yes</dd><dt><span>tiedtke_tracer_smooth_massflux :</span></dt><dd>No</dd><dt><span>tiedtke_shallow_wstar_closure :</span></dt><dd>No</dd><dt><span>tiedtke_detrainment_rate_deep_land :</span></dt><dd>7.5e-05</dd><dt><span>tiedtke_detrainment_rate_deep_ocean :</span></dt><dd>7.5e-05</dd><dt><span>tiedtke_shallow_entrainment :</span></dt><dd>2.0</dd><dt><span>tiedtke_cloud_cover_evap_over_land :</span></dt><dd>0.05</dd><dt><span>tiedtke_cloud_cover_evap_over_ocean :</span></dt><dd>0.05</dd><dt><span>tiedtke_coeff_evap_over_land :</span></dt><dd>5.55e-05</dd><dt><span>tiedtke_coeff_evap_over_ocean :</span></dt><dd>5.55e-05</dd><dt><span>tiedtke_critical_rh_over_land :</span></dt><dd>0.9</dd><dt><span>tiedtke_critical_rh_over_ocean :</span></dt><dd>0.9</dd><dt><span>tiedtke_cloud_water_conv_over_land :</span></dt><dd>0.03</dd><dt><span>tiedtke_cloud_water_conv_over_ocean :</span></dt><dd>0.0014</dd><dt><span>tiedtke_cape_adjustment_timescale :</span></dt><dd>2400.0</dd><dt><span>uwpbl_advection_scheme :</span></dt><dd>0</dd><dt><span>uwpbl_cloud_evap_entr_incr_efficiency :</span></dt><dd>15.0</dd><dt><span>uwpbl_eddy_LS_stable_PBL_scaling :</span></dt><dd>1.5</dd><dt><span>uwpbl_czero :</span></dt><dd>5.869</dd><dt><span>uwpbl_nuk :</span></dt><dd>5.0</dd><dt><span>frequency :</span></dt><dd>mon</dd><dt><span>experiment_id :</span></dt><dd>historical</dd><dt><span>driving_experiment :</span></dt><dd>MIROC-MIROC5, historical, r1i1p1</dd><dt><span>driving_model_id :</span></dt><dd>MIROC-MIROC5</dd><dt><span>driving_experiment_name :</span></dt><dd>historical</dd><dt><span>model_id :</span></dt><dd>ICTP-RegCM4-7</dd><dt><span>CDO :</span></dt><dd>Climate Data Operators version 1.9.9rc1 (https://mpimet.mpg.de/cdo)</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:    (time: 432, bnds: 2, y: 337, x: 429)\n",
       "Coordinates:\n",
       "  * time       (time) object 1970-01-16 12:00:00 ... 2005-12-16 12:00:00\n",
       "    lon        (y, x) float64 ...\n",
       "    lat        (y, x) float64 ...\n",
       "  * x          (x) float64 -5.35e+06 -5.325e+06 -5.3e+06 ... 5.325e+06 5.35e+06\n",
       "  * y          (y) float64 -4.2e+06 -4.175e+06 -4.15e+06 ... 4.175e+06 4.2e+06\n",
       "Dimensions without coordinates: bnds\n",
       "Data variables:\n",
       "    time_bnds  (time, bnds) object ...\n",
       "    crs        |S1 ...\n",
       "    pr         (time, y, x) float32 ...\n",
       "Attributes: (12/165)\n",
       "    CDI:                                             Climate Data Interface v...\n",
       "    Conventions:                                     CF-1.7\n",
       "    history:                                         Mon Jun 27 15:34:28 2022...\n",
       "    source:                                          RegCM Model output file\n",
       "    institution:                                     Oak Ridge National Labor...\n",
       "    project_id:                                      CORDEX\n",
       "    ...                                              ...\n",
       "    experiment_id:                                   historical\n",
       "    driving_experiment:                              MIROC-MIROC5, historical...\n",
       "    driving_model_id:                                MIROC-MIROC5\n",
       "    driving_experiment_name:                         historical\n",
       "    model_id:                                        ICTP-RegCM4-7\n",
       "    CDO:                                             Climate Data Operators v..."
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "pr = xr.open_dataset('/home/mkibria/Academic/pr_of_thesis/regcm driven by gcm/MIROC5/miroc5_pr_197001-200512.nc')\n",
    "\n",
    "pr\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "c6e12424",
   "metadata": {},
   "outputs": [],
   "source": [
    "#pr = pr.sel(time=slice('1980-01-16T12:00:00.000000000','2015-12-16T00:00:00.000000000'))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "a7bc5951",
   "metadata": {},
   "outputs": [],
   "source": [
    "\n",
    "#interpolating to new lat and lon\n",
    "\n",
    "#new_lon = np.arange(85,94,0.5)\n",
    "#new_lat = np.arange(17,28,0.5)\n",
    "#pr = pr.interp(latitude=new_lat,longitude=new_lon)\n",
    "#pr\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "2100d68b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(432, 337, 429)\n"
     ]
    }
   ],
   "source": [
    "size = pr.pr.shape\n",
    "print(size)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "77d6e633",
   "metadata": {},
   "outputs": [],
   "source": [
    "for time in range(size[0]):\n",
    "        data = pr.pr[time,Ellipsis] * spd\n",
    "        nzd  = data.values[tuple(np.where(data >= rain_minimum))]\n",
    "        hist = np.histogram(nzd,\n",
    "                            bins = nrainbin-1,\n",
    "                            range = (rain_minimum, rain_maximum))\n",
    "        rain_bin = rain_bin + hist[0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "307f9480",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEaCAYAAAD+E0veAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAq6UlEQVR4nO3de5zcdX3v8dc7GxIJl1ACHUu4hJBsSKqIuoDoUVNLbWiJUKVKXBQEiWixeHpaL8WHpdV41GotCgEjl0RBENCjoAjWS+RiVBalAi65NBKIwRDCHWJCNp/zx+83y2QyO/Ob3bnP+/l47CMzv/ldPru/yXzme1dEYGZmBjCu2QGYmVnrcFIwM7NhTgpmZjbMScHMzIY5KZiZ2TAnBTMzG+ak0OUk9Uv6fq33reL6z0iaXstzjnCdpZI+kXHf70k6LX18uqTb6xudVVJ4T6y+nBRamKQHJG1JPzg3SrpC0p61vEZEXBURbxzNvpJC0oys15K0XNK7i865Z0SszR5x/UXE8RGxbKznSf8+z6b373eS/kNST/pa/t4+LekJST+VdLakcQXHL5W0LT0+//O2scbVjmp1T4ql78k/FPx9V9b6Gu3GSaH1zY+IPYFXAEcBHy3eQdL4hkdlWb0svX9/DrwdOKvgtfkRsRdwCPAp4EPAZUXHfyZNnPmfr2e5aDe9J2rwu55T8PedVZOg2piTQpuIiN8B3wNeAsPfQv9O0mpgdbrtBEl3F3zzPCJ/vKSDJH1T0iZJmyVdmG7fqXokPe/fS1or6VFJ/57/9lq4r6Rb00P+O/8NVtIfSfpOeo3H08cHpvsvAl4LXJjuf2HB9WakjydL+kp6/DpJHy2+tqTPpuf+raTjR/p7SXq5pF+m38S/Dryo4LUR40xf36VEk26/SNLnirbdKOkD5e8eRMT9wG2k96/otScj4gbgbcBpknbZpxJJ50u6XtKVkp4CTk//npdJejgtqXyioKTSI+lz6T3+raRz0nsxPn1937RkuiH9G30r3T5X0npJH5T0SHrukyT9laRVkh6T9M8Z470ujfdpSfdI6pX0kfS8D0kqLJUO35P0vXCHpM9Legw4v9x7Jz3mLEmD6bV+I+kV1f6Nu4WTQpuQdBDwV8CvCjafBBwDzEnf5JcD7wGmAF8CbpA0Mf0g+A6wDpgGTAWuKXO5vwH6SEonJwJnFO8QEa9LH76s4BvsOOAKkm++BwNbgAvT/c8j+VDMfys7p8R1vwhMBqYDrwfeCbyr4PVjgJXAfsBngMskqfgkkiYA3wK+CuwLXAe8pWCXEeOsYBmwoCBR7UdSAri60oGS5pAkxV+NtE9E/AJYn+43GicC1wP7AFel8W4HZgAvB94I5JPdWcDxwJEk9/mkonN9FZgE/Cnwx8DnC157MUmSnQp8DPgycCrwyjT2jylbO9H89Dp/RPJ3uYXk3kwF/o3kPTySY4C1aWyLKPPekfS3wPnptr2BNwGbC871f9PkeIekuRni7mwR4Z8W/QEeAJ4BniD5QF8M7J6+FsAbCva9GPh40fErSf6DHAtsAsaXuMbpwO0FzwOYV/D8fcAPy+w7o0z8RwKPFzxfDry7aJ8g+dDqAbYCcwpeew+wvODaawpem5Qe++IS130dsAFQwbafAp+oNs4Sv/Mg8Bfp43OAm8r8/gE8BTwO/A/wCWBcwb09rsQxPwPOSx8vBf6Q3v8ngEfLXOt84NaC57n077l7wbYFwI/Txz8C3lPw2nFpvOOBPwF2AH9U4jpzSZJoT/p8r/S4Ywr2uQs4qcJ7+3zgvwqezyd5rxefd58R7smDBcdWeu/cApw7QhzHpNeaCJwGPA0cVuv/y+300zX1jm3spIj4wQivPVTw+BCSqof3F2ybABwADAHrImJ7xmsWnnddeo6KJE0i+UY5j+TbH8BeknoiYqjC4ful8a4ruvbUgue/zz+IiOfSQkKphvcDgN9F+r++4Fy1iHMZybfi/0r/vaDC/q+IiDUV9ik0FXis4PlnI2KndiRJ/bzwLfq2iMhXoxW/H3YDHi4oTI0r2OeAov0LHx8EPBYRj48Q4+aCv9OW9N+NBa9vofR9KVZ8zKMlzrsnSUIsVhhvpffOQSRJeRcR8fOCp8skLSApkX8xQ/wdydVH7a3wQ+8hYFFE7FPwMykirk5fO1jZG+QOKnh8MMm37iz+DzCL5Fvj3iTf2AHyn0rlpuR9FHie5MOs8Nq/y3jtQg8DU4uqlg6uIs5yrgROlPQyYDZJNVVNSDqK5IOsbBfYSHqB5RtGC9tVit8PW4H9Ct4Pe0fEn6avPwwcWLD/QUXH7itpn9H+Lg1Q+LtWeu88BBxWxXmzvA86lpNC5/gycLakY5TYQ9JfS9oL+AXJh8Cn0u0vkvSaMuf6p7Qx9iDgXGCkHi8bSepw8/Yi+Yb3hKR9gX+psP+w9BvitcAiSXtJOgT4B5IP4WqtIKlL/3tJ4yW9GTi6ijhHFBHrgTtJ6sK/ERFbKhxSkaS9JZ1A0s5zZUTcM9ZzRsTDwPeBz6XnHyfpMEmvT3e5FjhX0tT0w/9DRcd+D1icvg92k/S64mu0igzvnUuBf5T0yvT/xgxJh0jaR9Jfpv8fxqclsNeRVDd1LSeFDhERAySNhxeS1GGvIal7zf+nmU9Sd/8gSWNmuf7u3yapF74b+C67dpPMO5+kyP2EpLcC/wnsTvLN7WfAzUX7XwCcnPZm+UKJ870feJakAfF24GskjedViYhtwJtJfv/HSX7XbxbsUinOSpYBLyVJDGNxo6SnSb7Jngf8Bzs3rI/VO0mqVX5D8ne4nqS9AJIvEd8Hfk3SyHsTSSLNV9+8g+Tb9/3AI8AHahhXPYz43omI60gao79G0mbwLZIOCLuRtPNsInkvvJ+kurarxypo52pX63aSAphZZT14V0m/NV8JTIuIHc2OpxaUdO+9JCIOqbizdTSXFMyqIGk3kiq1S9s5IUjaPR1bMF7SVJIqtP/X7Lis+ZwUzDKSNJukJ8yfkFRBtTMB/0pSrfQrkq62H6vpBZL5ip4p8VNxcJs1j6uPzMxsmEsKZmY2zEnBzMyGtfWI5v322y+mTZvW7DDMzNrKXXfd9WhE7F/qtbZOCtOmTWNgYKDZYZiZtRVJ60Z6zdVHZmY2rGVKCpJeC/STxDQnIl7d5JDMzLpOXUsKki5PF8y4t2j7PEkrJa2R9GGAiLgtIs4mmfe/5svumZlZZfWuPlpKMj3xsHTBl4tIFviYQ7JoyZyCXd5OhkVLzMys9uqaFCLiVnaeGx6S2SrXRMTadOKya0hWjELSwcCTEfHUSOeUtFDSgKSBTZs21St0M7Ou1IyG5qnsvEDGel5YDONMkmUSRxQRSyKiLyL69t+/ZI8qMzMbpWY0NJdawCJZuzAi07z2kuYD82fMmFHLuMzMul4zSgrr2XmVpwPJvrKXmZnVUTOSwp3ATEmHSpoAnALcUM0JIuLGiFg4efLkugRoZtat6t0l9WqSpRFnSVov6cx08fhzSJa8GwSujYj7qjzvfElLnnzyydoHbWbWxdp66uy+vr7wNBdmZtWRdFdE9JV6rS2nuXBJwcysPtoyKYylTWHjVRtZMW0Fy8ctZ8W0FWy8amMdIjQza08tM/dRI2y8aiMrF65kx3PJ0rpb121l5cKVAOT6c80MzcysJbRlSWG01Udrz1s7nBDydjy3g7Xnra1leGZmbastk8Joq4+2Pri1qu1mZt2mLZPCaE08eGJV2xvN7R1m1mxtmRRGW300fdF0xk3a+VceN2kc0xdNr2V4o5Jv79i6bivEC+0dTgxm1khtmRRGW32U688xa8ksJh4yEQQTD5nIrCWzWqKR2e0dZtYKuqr3ESSJoRWSQLF6tndsvGoja89by9YHtzLx4IlMXzS9Jf8GZtZ8bVlS6ET1au9wtZSZVaMtk0InjmiuV3uHq6XMrBptmRQ6cZbUerV31KJayr2izLpH17UptLJ6tHdMPHhiUnVUYnsWHgVu1l3asqRg2Y21Wmo01U8uWZi1L5cUOlz+2/xoex9VW/3kkoVZe3NS6AJjqZaqtvqpXMmiOAZ3lTVrPW1ZfdSJvY9aVbXVT1lLFu4qa9aa2jIpdGLvo1ZVba+orOMtsrZVuH3CrLFcfWQVVVP9NH3R9J3aFKB0ySJLiaJS+4Srn8xqz0nBaiprw3aWtopKpQk3aJvVniKi2TGMWl9fXwwMDDQ7DBuF4lIAJCWKwqqp5eOWQ6m3p8oklUMmcuwDx7oUYVaGpLsioq/Uay4pWFNkKVGUK02Uq34aqdrpyTueZPNNm50ozMpomZKCpHHAx4G9gYGIWFbpGJcUOlu50sTa89aOWFIASr6G2KnkUVwyMesW5UoKde19JOlySY9Iurdo+zxJKyWtkfThdPOJwFTgeWB9PeOy9lCu51O5rrIjzutU9P0n3z7hHk5mL6hrSUHS64BngK9ExEvSbT3AKuAvSD787wQWAG8CHo+IL0m6PiJOrnR+lxS620jtBiumrShdUhjBuEnjdm7Q3g3G7z2e7Y9tdzWTdaSmtSlExK2SphVtPhpYExFr0+CuISklPARsS/cZqmdc1hlG6ipbqltscdXRsB526eHE87B983bAvZqs+zRj8NpUkgSQtz7d9k3gLyV9Ebh1pIMlLZQ0IGlg06ZN9Y3U2lKpaqcDzj6gZHVTlq8fO57bwapzV7mKybpCM3ofqcS2iIjngDMrHRwRSyQ9DMyfMGHCK2senXWEUqWIya+ZvEt100gN1sWGNg8xtDnJIC49WCdrRlJYDxxU8PxAYEM1J4iIG4Eb+/r6zqplYNbZRqpu2qWqKYMdz+1g8LRBBt8x6HYH6yjNqD66E5gp6VBJE4BTgBuqOYEnxLNaKa5q6pnSgyaUKsyWMMTwZH6D7xrk9v1ud/WStb169z66GpgL7AdsBP4lIi6T9FfAfwI9wOURsWg053fvI6uH4l5NQ88MDTc8Z+UxENbKyvU+apnBa9WQNB+YP2PGjLNWr17d7HCsw5UaRJfVxENctWStp+OSQp5LCtYohaUHxlFdp2mPe7AW07QRzfXiNgVrtFx/jmMfOJa5O+Yye9nsXbq3lpUf9+DFhKwNtGVS8CI71kxjapwm7bl06qAbpK0leZZUs1Eo7t46XL1UxfQaHu9gragtSwquPrJWk69emn1ldVVLLjVYq2nLpODqI2tVo61acluDtQr3PjKrs2qrltyN1erNvY/MmqjaqqWt67YyeGoyQtolB2u0tkwKrj6ydrRT1VIG2zdvd5WSNVxbJgWzdlVtqcEN0dZoTgpmTVBtqcEN0dYobZkU3KZgnWBUpYbTBp0YrK7aMim4TcE6Sb7U0DOlp/LOQ7jEYHXVlknBrNPk+nO89tHXMvvK2RWrlNzOYPXkpGDWQqqpUnI7g9WD5z4ya0H5gWuDpw2WnaY7385QeIzZWLikYNaicv25bNN0u53Baqgtk4J7H1m3yNp11T2TrFbaMim495F1k8ztDC4xWA2UbVOQ9IoM53g+Iu6pUTxmNoIs7QxuY7CxqtTQ/BPgTqDc3L+HAtNqFZCZjSz/Qb9y4Up2PLej9E5DMHjqIKvPXc3MC2Y6OVhVKiWFOyPiDeV2kPSjGsZjZhVk7ZmUn1Cv8BizSsq2KVRKCFn3MbPaytozyQ3QVq1MDc2S5tU7EElzJd0m6RJJc+t9PbN2l++ZRKXZMdwAbVUomxQkHS2pB/hkwbavZj25pMslPSLp3qLt8yStlLRG0ofTzQE8A7wIWJ/5NzDrYi4xWK1VKiksIGlsni7p05LeDmTpkZS3FNiplJEmmYuA44E5wAJJc4DbIuJ44EPAv1ZxDbOulnlCvbQB2iu6WTmVksKHIuJ/AQ8C3wX2BV4s6WeSvl7p5BFxK/BY0eajgTURsTYitgHXACdGRL4rxePAiCN1JC2UNCBpYNOmTZVCMOsKhRPqVapO8opuVk6l3ke3SNoO7A9MAW4CzoiIV0k6cJTXnAo8VPB8PXCMpDcDfwnsA1w40sERsQRYAtDX1xejjMGsI2XqsorHM9jIyiaFiPgzSbsDvwYOA/4WmCHpW8DtwGdHcc1SYx4iIr4JfDPTCaT5wPwZM2aM4vJmnS1rl9V8A3ThMWYVex9FxBbgoYj4bES8Hfgf4Czgt6O85nrgoILnBwIbqjmBp7kwK6+aBuhV565qUFTWDjJ1SS0ai3BhRGyKiG+M8pp3AjMlHSppAnAKcEM1J/CEeGaVZW2AHto85PYFG1b1hHgRcVnWfSVdDawAZklaL+nMiNgOnAPcAgwC10bEfVXG4JKCWQZZG6DdXdXyKk2ItyQiFo52n4hYMML2m0garUfFbQpm1RluZzh1sPQObl+wlCJG7sAj6RGSLqMj7gLMi4iZtQ4si76+vhgYGGjGpc3a0u373c72zdtH3qEHZi+b7cTQ4STdFRF9pV6r1CX1nzKc/7bqQxoblxTMRmfmBTMrz7D6jkGevONJehf3NjY4awllSwqtziUFs+ptvGpj5e6qgtlfdYmhU5UrKbTlymvufWQ2epm6q4Ybn7tVWyYF9z4yG5tMM6x6dtWuVDEpSOqR9O+NCMbMGidfYii3rqIHt3WfLCOah4BXSiq3JGdDufrIrDZy/TkOOPuAsolhaPOQZ1btIpkamiV9DpgJXAc8m9+ezlfUNG5oNquNLI3P4yaNY9aSWW587gBj6ZKaty+wGSic7iLIOIGdmbW2ioPbeKEqyUmhs2VKChHxrnoHYmbNlevPsfrc1WUHtw1tHmLV+1Z5DEMHy7pGc6+kH+aX1ZR0hKSP1jc0M2u0mRfMrDiz6oZLNrh9oYNl7ZL6ZeAjwPMAEfFrktlNm8INzWb1kWlmVY9h6GhZk8KkiPhF0bYyE6jUl8cpmNVPfmbV8VPK1C6n02Gsep+7q3aarEnhUUmHkTQuI+lk4OG6RWVmTTfzgpllu6oSsOHiDU4MHSZrUvg74EvA4ZJ+B3wAOLteQZlZ82UZwwBuY+g0WVdeWxsRxwH7A4dHxP+KiHX1Dc3Mmq13cS+zv1p+gR4Cj3ruIFl7H02R9AWSabKXS7pA0pT6hlY2Hjc0mzVIlukwPOq5c2StProG2AS8BTg5ffz1egVViRuazRpruCqpjO2bt7vxuQNkHtEcER8veP4JSSfVIR4za1H5AWsbLt4w8k7xwuse4NaespYUfizpFEnj0p+3At+tZ2Bm1np6F/eW76qacq+k9lU2KUh6WtJTwHuArwHb0p9rgP9d//DMrNVkGfUMTgztqmzKj4i9GhWImbWH/IR4q85dxdDmcmt6uiqpHWVtU0DSEcC0wmOaPXW2mTVHrj9Hrj/HqvetKt/GgBNDu8mUFCRdDhwB3AfsSDfXfOpsSXsAtwL/EhHfqeW5zaz2MjU+48TQTrKWFF4VEXOqPXmaTE4AHomIlxRsnwdcQDIk5tKI+FT60oeAa6u9jpk1jxNDZ8na+2iFpKqTArAUmFe4QVIPcBFwPDAHWCBpjqTjgN8AHv1i1mZ6F/dywHvLj2MANz63g6wlhWUkieH3wFaSsY0REUeUOygibpU0rWjz0cCaiFgLIOka4ERgT2APkkSxRdJNEbGj6FgkLQQWAhx88MEZwzezenOJoTNkTQqXA+8A7uGFNoXRmgo8VPB8PXBMRJwDIOl04NFSCQEgIpYASyBZo3mMsZhZDTkxtL+sSeHBiLihRtcsNYPK8Id7RCyteAJpPjB/xowZNQrJzGrFiaG9ZU0K90v6GnAjSfURMOouqeuBgwqeHwiUf/eYWVtxYmhfWRuadydJBm8E5qc/J4zymncCMyUdKmkCybKeVZVCPCGeWetz43N7ylRSiIh3jebkkq4G5gL7SVpPMv7gMknnALeQdEm9PCLuq/K8rj4yawMuMbQfRVRuq5V0BQX1/nkRcUY9gsqqr68vBgYGmhmCmWWQZeQzwAHvPcCJoQEk3RURfaVey1p99B2SWVG/C/wQ2Bt4pjbhVc+L7Ji1l2qqkm7b6zYv1tNEmUoKuxwkjQN+EBFvqH1I2bmkYNZespYYwKWGeqpFSaHYTMAjx8ysKllLDOAG6GbJukbz05Keyv+QdE39UH1DKxuPq4/M2lS1icFVSY2VKSlExF4RsXfBT29EfKPewZWJx11SzdpYNYlh8B2DTgwNVGnltYPL/TQqSDPrPJkTQ8DgqYNugG6QSuMUvkvSFbVwaooA9gf+mGScQcN5nIJZZ+hd3Mvk10zm/vfcTzxbvtPL0DNDDJ46yJN3POkG6DoqW1KIiJdGxBHpvy8lGcl8B0l31A80IL6R4nL1kVmHyPXneP0zr6dnz2zfMd0AXV9ZG5pnSloKfA+4C5gTEV+sZ2Bm1l16L8n+7d+JoX4qtSm8JJ2q4hvAD4CXRMSlEfF8Q6Izs66R689lbnwGJ4Z6KTt4TdIQydoH3wWGil+PiL+vX2gjK2hTOGv16tXNCMHM6mTjVRsztTHkzb5yNrn+XJ2j6izlBq9VSgqnlTtxRCwbY2xj4hHNZp0r6+hn7SFe/8zrGxBR5yiXFMr2Pmr2h76Zda+sPZPi2eDu4+7myB8c2bjgOthop7kwM6u7rD2TnvjhEx7HUCNOCmbW8nov6a04Kio/jsGNz2OTtUvqa7JsMzOrh1x/jtnLZmfad8PFG/jJ7j9xqWGUspYUSo1JaNo4BU+IZ9Z9qumyGn8Iz5k0SpV6Hx0LvJpk9PLnC17aG/ibiHhZXaOrwL2PzLrP3cfdzRM/fCLTvu6ZVNpY1lOYAOxJ0ktpr4Kfp4CTaxmkmVkWR/7gyOwlhmfDbQxVqtQl9SfATyQtjYh1DYrJzKys/IR4WcYxbLh4A5NfM9kD3DLK2qYwUdISSd+X9KP8T10jMzMro3dxL7OvnJ3UZ1Rw/3vur39AHSJrUrgO+BXwUeCfCn7MzJom159j7ta57PPn+5Tdz9VI2WVNCtsj4uKI+EVE3JX/qWUgkmZLukTS9ZLeW8tzm1lny9LOsOHiDdx93N2NCaiNZU0KN0p6n6Q/kbRv/qfSQZIul/SIpHuLts+TtFLSGkkfBoiIwYg4G3grULJV3MxsJL2LezONfHZiKC9rUjiNpLropyTrKdwFZOkLuhSYV7hBUg9wEXA8MAdYIGlO+tqbgNuBH2aMy8xsWJY1GZ744RMev1BGpqQQEYeW+Jme4bhbgceKNh8NrImItRGxDbgGODHd/4aIeDXQX92vYWaWtDFkWcFt8FQPbBtJ1mkuJkn6qKQl6fOZkk4Y5TWnkqzRkLcemCpprqQvSPoScFOZWBZKGpA0sGnTplGGYGadKss8SYDnSRpB2XEKBa4gqTJ6dfp8PUmPpO+M4poqsS0iYjmwvNLBEbFE0sPA/AkTJrxyFNc3sw6WH48weMYgbCu/74ZLPIahWNY2hcMi4jPA8wARsYXSH+5ZrAcOKnh+IFB5BEqBiLgxIhZOnjx5lCGYWSfLd1XViyp8TAWsOtelhUJZk8I2SbsDASDpMGDrKK95JzBT0qGSJgCnADdUcwJPiGdmWRx+6eEV9xnaPORqpAJZk8L5wM3AQZKuIukd9MFKB0m6GlgBzJK0XtKZEbEdOAe4BRgEro2I+6oJ2iUFM8si15+rOLANkjEMTgyJsrOk7rSjNAV4FUm10c8i4tF6BlYhlvnA/BkzZpy1evXqZoVhZm0i63rPs6+c3RXtC2OZJTV/ghuANwLLI+I7zUwI4JKCmVWnd3EvEw+ZWHE/z5GUvfroc8Brgd9Iuk7SyZJeVMe4ynKbgplVa/qiikOriGej68cvZB289pOIeB8wHVhCMhXFI/UMrEI8LimYWVWyrtzW7b2RspYUSHsfvQU4GzgKWFavoMzM6qF3cS/jp5QfnjW0eairSwtZ2xS+TtJT6A0k8xYdFhHvr2dgFeJx9ZGZjcrMC2bCbuX3GTxjsDHBtKCsJYUrSBLB2RHxo4jYUc+gKnH1kZmNVq4/x+wrZpffaRtd20W1bFKQ9EGAiLgZeHPRa5+sY1xmZnWT689V7I204eIN3L7f7V1XlVSppHBKweOPFL02jyZx9ZGZjVWW3kjbN2/n/jPu76rEUCkpaITHpZ43jKuPzGyssvZGim3B2vPWNiCi1lApKcQIj0s9NzNrK72LezNNg7H1wdFO9dZ+KiWFl0l6StLTwBHp4/zzlzYgPjOzujryB0dW3qlCb6VOUjYpRERPROwdEXtFxPj0cf550/5MblMws1rqmVJhVZ4u6o2UefBaK3GbgpnVUu8Fldd23nDxhq5ocG7LpGBmVktZG527YcI8JwUzM7I1OnfDhHlOCmZmqSyNzp3ePdVJwcysQKWRzlvXdXb31LZMCu59ZGb1Mn3R9IpdUDu5Cqktk4J7H5lZvQxPmFdmzobBUwc7totqWyYFM7N6yrJO84aLN3RkYnBSMDMrYeLBldd07sSxC04KZmYlZJlFFWDlwpUdlRicFMzMSsj15+jZs8L0F8CO53Yw+M7BjkkMLZUUJJ0k6cuSvi3pjc2Ox8y6W+8llae/AGAH3H96Z6y7UPekIOlySY9Iurdo+zxJKyWtkfRhgIj4VkScBZwOvK3esZmZlZPrz6E9si0dE9ujI6bBaERJYSlFq7RJ6gEuAo4H5gALJM0p2OWj6etmZk11+JcOz7xvPBss13KWa3nbLuVZ96QQEbcCjxVtPhpYExFrI2IbcA1wohKfBr4XEb8sdT5JCyUNSBrYtGlTfYM3s66XdbK8Yu26lGez2hSmAg8VPF+fbns/cBxwsqSzSx0YEUsioi8i+vbff//6R2pmXa93cS+zr5xd9XHtuJTn+CZdt1QlXUTEF4AvVDxYmg/MnzFjRs0DMzMrJdefY/Adg1UvRLx13VaWazkA46eMZ+YFMzMNjmuWZpUU1gMHFTw/ENiQ9WBPc2FmzZBlQFs57VCl1KykcCcwU9KhkiYApwA3ZD3YE+KZWTNkmSyvktgWrDq3dafHaESX1KuBFcAsSeslnRkR24FzgFuAQeDaiLiv3rGYmY1FfrK8ims6VzC0eYjlWs6KaStartSgiCoryFpIX19fDAwMNDsMM+tCK6atqMnaCuMmjWPWklkNbWeQdFdE9JWMp2FRmJl1kOmLpjNu0tg/Qnc8t4PBUwdbptTQlknBbQpm1my5/hyzlsxKVmrLNui5rK3rtjL4jsGmVyu5+sjMrAZWvW8VGy7O3IkykwPeewC9izPOv1SFjqs+cknBzFpN7+LeUY18LqcZ6zW0ZVLwOAUza0W9i3uT6qQaanT31bZMCmZmrapWDdB5Q5uHGlpaaMuk4OojM2tVOzVA10gj509qy6Tg6iMza2W5/hzHPnAsc2NuTdoZtq7b2rDSQlsmBTOzdpGfYVUTxtZvtVFrQTspmJnVWa4/x+GXHz48pqFnSg9MqO4cO57bweBp9V8LullTZ4+Jp842s3aT68/tMpVFfkrtzIaSEkP+fPXQliUFtymYWScYTWN0vsSwfFx9Rj63ZVIwM+sE0xdNH92BQ0AkDdC1bmtwUjAza5Jcf47xU8ZWi7/juR017bLqpGBm1kQzL5g55sFuWx8c+xTeeW3Z0Gxm1inyDcaDpw6O+hxjXSa0UFuWFDyi2cw6Sa4/N+pBbuMmjRt920Sp89XsTA3k3kdm1ml6F/eyz5/vU/VxLz7txTXtntqWScHMrBNtWbOl6mM2ftVdUs3MOtJoGoyHnhmqaQxOCmZmLaKWDcaj5aRgZtYiar0Ww2g4KZiZtYhRrcUwtslXd9EySUHSdEmXSbq+2bGYmTVL4VoMmT7wo7bXr2tSkHS5pEck3Vu0fZ6klZLWSPowQESsjYgz6xmPmVk7ydLGUOs1oetdUlgKzCvcIKkHuAg4HpgDLJA0p85xmJm1nemLpsNuGfapobomhYi4FXisaPPRwJq0ZLANuAY4Mes5JS2UNCBpYNOmTTWM1systeT6c8y+YvaIn9Q9U3pqvq5CM9oUpgIPFTxfD0yVNEXSJcDLJX1kpIMjYklE9EVE3/7771/vWM3MmirXn2P2V2bv0itp3KRx9F7QW/PrNWNCvFJNJxERm4GzM53AK6+ZWRfJ9ed48o4n2bBkQ7KWQk/tp7fIa0ZJYT1wUMHzA4ENTYjDzKwtbLxqI79f9vskIQAMwe+X/b4u6zU3IyncCcyUdKikCcApwA3VnMAT4plZN1l73lp2PLdjp221Xlwnr95dUq8GVgCzJK2XdGZEbAfOAW4BBoFrI+K+Ks/rqbPNrGuMNCdSLRfXyatrm0JELBhh+03ATWM4743AjX19fWeN9hxmZu1i4sET2bpu1wRQj7mSWmZEczVcUjCzblJqTqRaL64zfN6an7EB3KZgZt1kpzmRlIxinrVkVl16H3mNZjOzNpDrz9UlCRRry5KCq4/MzOqjLZOCq4/MzOqjLZOCmZnVR1smBVcfmZnVR1smBVcfmZnVhyJqvGxPA0naBKxLn04GShUdRtq+H/BonUIbq5FibvZ5qz0+6/5Z9iu3T7Wv+d7X/3jf++o18t4fEhGlp5mOiI74AZZUuX2g2TFX+7s0+7zVHp91/yz7ldun2td8733vW/GnVe59W1YfjeDGKre3snrFPNbzVnt81v2z7Fdun9G+1op876vbx/e+xudt6+qjsZA0EBF9zY7DGs/3vnv53lfWSSWFai1pdgDWNL733cv3voKuLSmYmdmuurmkYGZmRZwUzMxsmJOCmZkNc1JISdpD0jJJX5bU3+x4rHEkTZd0maTrmx2LNZakk9L/89+W9MZmx9MKOjopSLpc0iOS7i3aPk/SSklrJH043fxm4PqIOAt4U8ODtZqq5t5HxNqIOLM5kVqtVXnvv5X+nz8deFsTwm05HZ0UgKXAvMINknqAi4DjgTnAAklzgAOBh9LdhhoYo9XHUrLfe+ssS6n+3n80fb3rdXRSiIhbgceKNh8NrEm/HW4DrgFOBNaTJAbo8L9LN6jy3lsHqebeK/Fp4HsR8ctGx9qKuvHDbyovlAggSQZTgW8Cb5F0Me03PN6yKXnvJU2RdAnwckkfaU5oVmcj/b9/P3AccLKks5sRWKvpxjWaVWJbRMSzwLsaHYw11Ej3fjPgD4TONtK9/wLwhUYH08q6saSwHjio4PmBwIYmxWKN5XvfvXzvM+rGpHAnMFPSoZImAKcANzQ5JmsM3/vu5XufUUcnBUlXAyuAWZLWSzozIrYD5wC3AIPAtRFxXzPjtNrzve9evvdj4wnxzMxsWEeXFMzMrDpOCmZmNsxJwczMhjkpmJnZMCcFMzMb5qRgZmbDnBSsKSQNSbpb0r2SrpM0qQbn/DdJx1XY52xJ70wfny7pgAzn3Wk/SZfWYnZVSedL+p2kfxvrucpc42ZJU4u2zZX0nVGe723p1NOjOt5an5OCNcuWiDgyIl4CbKNo7qF0quOqRMTHIuIHFfa5JCK+kj49HaiYFIr3i4h3R8Rvqo1vBJ+PiI/V6Fw7kbQ7sG9E/K5W54yIrwPvrtX5rPU4KVgruA2YkX6D/bGkrwH3SOqR9O+S7pT0a0nvyR8g6YOS7pH035I+lW5bKunk9PEDkj4t6Rfpz4x0+/mS/jHdrw+4Ki2x7C7pY+m17pW0JJ1WudR+yyX1pedbkMZxbzoFcz6+ZyQtSuP7maRcpT9CGtsySd9P43+zpM+k579Z0m4Fv9snJa2QNCDpFZJukfQ/RTN9zgWWp8fMk3S/pNtJFpTKX/NoST+V9Kv031np9tskHVmw3x2Sjsh+S61dOSlYU0kaT7LwyT3ppqOB8yJiDnAm8GREHAUcBZyVzl1zPHAScExEvAz4zAinfyoijgYuBP6z8IWIuB4YAPrTEssW4MKIOCotvewOnDDCfvnYDwA+DbwBOBI4StJJ6ct7AD9L47sVOCvjn+Qw4K9J1nm4EvhxRLwU2JJuz3soIo4lSahLgZOBVwGFVVHHAzdLehHwZWA+8FrgxQX73A+8LiJeDnwM+GS6/VKSEhKSeoGJEfHrjL+DtTEnBWuW3SXdTfKB+yBwWbr9FxHx2/TxG4F3pvv9HJgCzCSZ//6KiHgOICKKF1TJu7rg32MzxPRnkn4u6R6SD/o/rbD/UcDyiNiUzq1zFfC69LVtQL7e/S5gWobrQ7LYy/MkSbIHuDndfk/ROW4o2P7ziHg6IjYBf5C0T/raa4DbgcOB30bE6kjmtbmy4DyTgeuULF35eV74na8DTkhLJ2eQJB7rAt24noK1hi0RcWThBkkAzxZuAt4fEbcU7TcPyDJpV4zweBfpt+nFQF9EPCTpfOBFFc5fao7+vOfjhYnFhsj+f20rQETskFR4jh1F59hasH1rwfYdwHhJ00lKE9vSv+tIv//HSUojfyNpGml1U0Q8J+m/SEosbyWpQrMu4JKCtbJbgPcW1KX3StoD+D5whtIeS5L2HeH4txX8u6LE608De6WP8wngUUl7klTHlNqv0M+B10vaL20YXwD8JNNvVn/H80Ip437gUEmHpc8XFOw3Gcg3RJ9edI5LSRagubNMacw6jEsK1souJaky+aWSr7ubgJMi4ua0EXRA0jbgJuCfSxw/UdLPSb78LCjx+lLgEklbSKqXvkxSHfMAyfz7I+0HQEQ8rGT5zh+TlBpuiohvj/aXrbF5JEtNEhF/kLQQ+K6kR0mqlF6S7vcZYJmkfwB+VHiCiLhL0lPAFY0L25rNU2dbR5L0AElV0KPNjmUkaRXVMxHx2RqfdyJwR0SMqconbUhfDhweETsKts8F/jEiThjL+a01ufrIrHmeARaqxoPXImJrDRLCO0mqx84rSghvI2l7eXxsUVqrcknBzMyGuaRgZmbDnBTMzGyYk4KZmQ1zUjAzs2FOCmZmNsxJwczMhv1/9zjp+bltsIkAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig = plt.figure()\n",
    "ax = fig.add_subplot(1, 1, 1)\n",
    "plt.plot(rain_bin, 'mo')\n",
    "ax.set_yscale('log')\n",
    "ax.set_xscale('log')\n",
    "plt.xlabel('Precipitation [mm/day]')\n",
    "plt.ylabel('Event Number [#]')\n",
    "plt.title('Precipitation daily PDF-regcm_miroc5')\n",
    "plt.savefig(\"/home/mkibria/Academic/pr_of_thesis/figures/pdf_regcm_miroc5.png\")\n",
    "\n",
    "\n",
    "\n",
    "plt.show( )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a43f36ef",
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.7"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}