pdf_regcm_eraint.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: 467, bnds: 2, y: 337, x: 429)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 1979-01-16T12:00:00 ... 2017-11-16\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) datetime64[ns] 1979-01-01 1979-02-01 ... 2017-12-01\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:                                         Wed Jul 13 17:38:34 2022...\n",
       "    source:                                          RegCM Model output file\n",
       "    institution:                                     International Centre for...\n",
       "    project_id:                                      CORDEX\n",
       "    ...                                              ...\n",
       "    experiment_id:                                   evaluation\n",
       "    driving_experiment:                              ECMWF-ERAINT, evaluation...\n",
       "    driving_model_id:                                ECMWF-ERAINT\n",
       "    driving_experiment_name:                         evaluation\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-3c8638d8-f8bd-4c2a-b174-65248727167e' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-3c8638d8-f8bd-4c2a-b174-65248727167e' 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>: 467</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-91ab8938-f109-4a17-b33c-b95759a2e324' class='xr-section-summary-in' type='checkbox'  checked><label for='section-91ab8938-f109-4a17-b33c-b95759a2e324' 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'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>1979-01-16T12:00:00 ... 2017-11-16</div><input id='attrs-358e24ed-91f7-4fe2-af28-65134f5ebaab' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-358e24ed-91f7-4fe2-af28-65134f5ebaab' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-14525554-9180-4a10-87b0-7998a720abe0' class='xr-var-data-in' type='checkbox'><label for='data-14525554-9180-4a10-87b0-7998a720abe0' 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([&#x27;1979-01-16T12:00:00.000000000&#x27;, &#x27;1979-02-15T00:00:00.000000000&#x27;,\n",
       "       &#x27;1979-03-16T12:00:00.000000000&#x27;, ..., &#x27;2017-09-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;2017-10-16T12:00:00.000000000&#x27;, &#x27;2017-11-16T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</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-b1fc86a9-e15c-4f3a-93b7-8e7b17725e3a' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-b1fc86a9-e15c-4f3a-93b7-8e7b17725e3a' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-80deb1f0-d8f5-48dc-9de3-919263767248' class='xr-var-data-in' type='checkbox'><label for='data-80deb1f0-d8f5-48dc-9de3-919263767248' 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-41e0e9f2-70a6-41ee-9224-38442f60b737' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-41e0e9f2-70a6-41ee-9224-38442f60b737' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-7a6507a6-a686-4c65-a54d-78a74c1147c2' class='xr-var-data-in' type='checkbox'><label for='data-7a6507a6-a686-4c65-a54d-78a74c1147c2' 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-70390108-cdc8-4928-b765-b0538d96690b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-70390108-cdc8-4928-b765-b0538d96690b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0b015430-0b0b-42d7-9a8f-16733eab16d7' class='xr-var-data-in' type='checkbox'><label for='data-0b015430-0b0b-42d7-9a8f-16733eab16d7' 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-58f98ab6-e260-457f-8623-94814ad35766' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-58f98ab6-e260-457f-8623-94814ad35766' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-068e0bde-9ffe-40f2-94ea-6f1e564aee35' class='xr-var-data-in' type='checkbox'><label for='data-068e0bde-9ffe-40f2-94ea-6f1e564aee35' 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-1baf69d9-ab1b-429b-9310-20b9bb991926' class='xr-section-summary-in' type='checkbox'  checked><label for='section-1baf69d9-ab1b-429b-9310-20b9bb991926' 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'>datetime64[ns]</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-b0ef7fa1-2d12-4230-9516-f2bf24a52fbf' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-b0ef7fa1-2d12-4230-9516-f2bf24a52fbf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-0c82e20b-4e39-4a85-a2c2-50f6cc8841b8' class='xr-var-data-in' type='checkbox'><label for='data-0c82e20b-4e39-4a85-a2c2-50f6cc8841b8' 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([[&#x27;1979-01-01T00:00:00.000000000&#x27;, &#x27;1979-02-01T00:00:00.000000000&#x27;],\n",
       "       [&#x27;1979-02-01T00:00:00.000000000&#x27;, &#x27;1979-03-01T00:00:00.000000000&#x27;],\n",
       "       [&#x27;1979-03-01T00:00:00.000000000&#x27;, &#x27;1979-04-01T00:00:00.000000000&#x27;],\n",
       "       ...,\n",
       "       [&#x27;2017-09-01T00:00:00.000000000&#x27;, &#x27;2017-10-01T00:00:00.000000000&#x27;],\n",
       "       [&#x27;2017-10-01T00:00:00.000000000&#x27;, &#x27;2017-11-01T00:00:00.000000000&#x27;],\n",
       "       [&#x27;2017-11-01T00:00:00.000000000&#x27;, &#x27;2017-12-01T00:00:00.000000000&#x27;]],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</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-6a17f1a1-fbec-41ec-84de-978be25de899' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-6a17f1a1-fbec-41ec-84de-978be25de899' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c4ea2d1d-f403-4dd8-830d-d74bd78dc8ee' class='xr-var-data-in' type='checkbox'><label for='data-c4ea2d1d-f403-4dd8-830d-d74bd78dc8ee' 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-fe0c2347-f3fb-43c1-9b7a-612114a481f2' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-fe0c2347-f3fb-43c1-9b7a-612114a481f2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4266cdc4-cfe5-4a82-ba9e-88d4516e6a3e' class='xr-var-data-in' type='checkbox'><label for='data-4266cdc4-cfe5-4a82-ba9e-88d4516e6a3e' 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>[67515591 values with dtype=float32]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-fdfa2bfe-94d5-477d-bb8d-e84ed88432c0' class='xr-section-summary-in' type='checkbox'  ><label for='section-fdfa2bfe-94d5-477d-bb8d-e84ed88432c0' 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>Wed Jul 13 17:38:34 2022: cdo mergetime pr_WAS-22_ECMWF-ERAINT_evaluation_r1i1p1_ICTP-RegCM4-7_v0_mon_197901-198012.nc pr_WAS-22_ECMWF-ERAINT_evaluation_r1i1p1_ICTP-RegCM4-7_v0_mon_198101-199012.nc pr_WAS-22_ECMWF-ERAINT_evaluation_r1i1p1_ICTP-RegCM4-7_v0_mon_199101-200012.nc pr_WAS-22_ECMWF-ERAINT_evaluation_r1i1p1_ICTP-RegCM4-7_v0_mon_200101-201012.nc pr_WAS-22_ECMWF-ERAINT_evaluation_r1i1p1_ICTP-RegCM4-7_v0_mon_201101-201711.nc pr_was-22-ecmwf_197901-201711.nc\n",
       "2018-02-13 16:42:58 : Created by RegCM RegCM Model program</dd><dt><span>source :</span></dt><dd>RegCM Model output file</dd><dt><span>institution :</span></dt><dd>International Centre for Theoretical Physics</dd><dt><span>project_id :</span></dt><dd>CORDEX</dd><dt><span>ipcc_scenario_code :</span></dt><dd>EVALUATION</dd><dt><span>institute_id :</span></dt><dd>ICTP</dd><dt><span>note :</span></dt><dd>The domain is larger than WAS-22</dd><dt><span>comment :</span></dt><dd>RegCM CORDEX India25 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-12T15:13:03Z</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>none</dd><dt><span>contact :</span></dt><dd>esp@ictp.it</dd><dt><span>product :</span></dt><dd>output</dd><dt><span>tracking_id :</span></dt><dd>f1ed9934-3b9a-11e9-8719-0894ef50c218</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>EIN75</dd><dt><span>model_sst_data_source :</span></dt><dd>EIN75</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> 1979-01-01 00:00:00 UTC</dd><dt><span>model_simulation_start :</span></dt><dd> 1979-01-01 00:00:00 UTC</dd><dt><span>model_simulation_end :</span></dt><dd> 2016-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>evaluation</dd><dt><span>driving_experiment :</span></dt><dd>ECMWF-ERAINT, evaluation, r1i1p1</dd><dt><span>driving_model_id :</span></dt><dd>ECMWF-ERAINT</dd><dt><span>driving_experiment_name :</span></dt><dd>evaluation</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: 467, bnds: 2, y: 337, x: 429)\n",
       "Coordinates:\n",
       "  * time       (time) datetime64[ns] 1979-01-16T12:00:00 ... 2017-11-16\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) datetime64[ns] ...\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:                                         Wed Jul 13 17:38:34 2022...\n",
       "    source:                                          RegCM Model output file\n",
       "    institution:                                     International Centre for...\n",
       "    project_id:                                      CORDEX\n",
       "    ...                                              ...\n",
       "    experiment_id:                                   evaluation\n",
       "    driving_experiment:                              ECMWF-ERAINT, evaluation...\n",
       "    driving_model_id:                                ECMWF-ERAINT\n",
       "    driving_experiment_name:                         evaluation\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/Era_intrim/pr_was-22-ecmwf_197901-201711.nc')\n",
    "\n",
    "pr\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "72e85834",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "-4207954"
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#pr.lon.values\n",
    "\n",
    "\n",
    " "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "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": 7,
   "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": 8,
   "id": "2100d68b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(431, 337, 429)\n"
     ]
    }
   ],
   "source": [
    "size = pr.pr.shape\n",
    "print(size)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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": 10,
   "id": "307f9480",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAwHUlEQVR4nO3de5icdX338fcnu+xyBtngQgJJCGxCUo+IKLVoan1soFCo5VFiRERKBIrFyx7E0sdHW+OjVmtjOTUiECUF8dAWFEXbmkJsVAIeAJccjATCwhKCAjGYsLvf54/7ns1kMufDzs7M53Vde2Xmnvvwu3c2853f6ftTRGBmZgYwpdkFMDOzycNBwczMxjkomJnZOAcFMzMb56BgZmbjHBTMzGycg4KNk7RY0rfrvW8F198uaXY9z1ngOjdK+miZ+35T0nnp43dJWt3Y0lmtJF0r6f80uxytykGhRUh6WNLz6QfnsKQbJB1Yz2tExMqIeHM1+0oKSceVey1JqyT9Sc45D4yITeWXuPEi4tSIWFHredLfz6/T9+8xSf8gqSt9LfPePifpV5L+R9JFkqZkHX+jpF3p8Zmft9VarnYUERdFxN+Vs28lXxA6hYNCazkjIg4ETgBeDfxN7g6Suie8VFaul6fv3+8BbwcuzHrtjIg4CJgJfBz4APD5nOM/mQbOzM+XyrloO/1NtNO9TFYOCi0oIh4Dvgm8BMa/hf6ppA3AhnTb6ZJ+nPXN82WZ4yUdLelrkrZK2ibpynT7Hs0j6Xn/TNImSU9J+vvMt9fsfSXdlR7yk8w3WEkvkvT19Bq/TB8fle6/FDgFuDLd/8qs6x2XPj5E0hfS4zdL+pvca0v6VHruX0g6tdDvS9IrJd2XfhP/ErBv1msFy5m+vleNJt1+laRP52y7XdL7ir97EBEPAXeTvn85rz0TEbcBbwPOk7TXPqVI+rCkr0i6SdKzwLvS3+fnJT2e1lQ+mlVT6ZL06fQ9/oWkS9P3ojt9/bC0ZjqU/o7+Ld2+QNIWSX8l6cn03GdJOk3SeklPS/rrMso7RdLlkn6e/j3eKumw9LVZaVkukPQI8F/p9i9LekLSM5LukvRbWecb//afVcY/zyrj+elrS4DFwF+lf4e3V/q7bkcOCi1I0tHAacCPsjafBbwGmC/pBOB64D1AH/DPwG2SetMPgq8Dm4FZwHTgliKX+yPgRJLayZnAu3N3iIjXpw9fnvUNdgpwA8k33xnA88CV6f5XkHwoXpruf2me6/4TcAgwG3gD8E7g/KzXXwOsA6YCnwQ+L0m5J5HUA/wb8EXgMODLwB9n7VKwnCWsABZlBaqpJDWAm0sdKGk+SVD8UaF9IuKHwJZ0v2qcCXwFOBRYmZZ3BDgOeCXwZiAT7C4ETgVeQfI+n5Vzri8C+wO/BbwY+EzWa0eQBNnpwIeAzwHvAF6Vlv1DKt1P9GfpNd8ATAN+CVyVs88bgHnA76fPvwkMpOW5L73HQo4g+VuaDlwAXCXpRRGxPD0uUwM7o0Q5O0NE+KcFfoCHge3Ar0g+0K8G9ktfC+CNWfteA/xdzvHrSP5jnQxsBbrzXONdwOqs5wEszHp+CfCfRfY9rkj5XwH8Muv5KuBPcvYJkg+tLmAnMD/rtfcAq7KuvTHrtf3TY4/Ic93XA0OAsrb9D/DRSsuZ554Hgf+VPr4UuKPI/QfwLMkH3s+BjwJTst7bN+U55vvAFenjG4HfpO//r4Cnilzrw8BdWc/709/nflnbFgHfTR//F/CerNfelJa3GzgSGANelOc6C0iCaFf6/KD0uNdk7XMvcFaJv+1B4Peynh8JvJBef1Z6ztlFjj803eeQrN/VR3PK2J21/5PAa3P39U/y4/a51nJWRPxHgdcezXo8k6Tp4b1Z23pIvoWNApsjYqTMa2afd3N6jpIk7U/yjXIh8KJ080GSuiJitMThU9Pybs659vSs509kHkTEjrSSkK/jfRrwWKSfAFnnqkc5V5B8K/5O+u+yEvufEBEbS+yTbTrwdNbzT0XEHv1IkhaT1AQB7o6ITDNa7t/DPsDjWZWpKVn7TMvZP/vx0cDTEfHLAmXclvV7ej79dzjr9efJ/75kmwn8q6SxrG2jJMFsrzKltd2lwP8GDicJWpD83TxToIzZf+87yihTx3LzUfvI/tB7FFgaEYdm/ewfETenr81Q+R12R2c9nkHyrbscfw7MJfnWeDDJN3aAzKdSsfS8T5F8U5yZc+3Hyrx2tseB6TlNSzMqKGcxNwFnSno5SdPGv1VRvrwkvZokKBQdAhvJKLBMx3N2v0ru38NOYGrW38PBEZFph38cOCpr/6Nzjj1M0qHV3ksZHgVOzfl73TeSvrOM7Pt5O0nz2JtImoVmpdvLec9yOU10DgeF9vQ54CJJr1HiAEl/IOkg4IckHwIfT7fvK+l1Rc71l2ln7NHAZUChES/DJO3/GQeRfEv8Vdpp+H9L7D8u/eZ5K7BU0kGSZgLvJ/kQrtQakrb0P5PULektwEkVlLOgiNgC3EPS5v7ViHi+xCElSTpY0ukk/Tw3RcT9tZ4zIh4Hvg18Oj3/FEnHSnpDusutwGWSpqcf/h/IOfabwNXp38E+kl6fe40aXUvyXs8EkHS4pDOL7H8QSZDbRtJ0+LEarl3w77BTOSi0oYhYS9J5eCVJG/ZGkvbwzAfuGSRt94+QdGYWG+/+7yTtwj8GvsHewyQzPgysUDLa6a3APwL7kXzr/z7wrZz9lwFnp6NZPpvnfO8Ffg1sIvm2/C8knecViYhdwFtI7v+XJPf6taxdSpWzlBXAS0kCQy1ul/QcybfmK4B/YM+O9Vq9k6RJ7mckv4evkLTdQ/Il4tvAT0k6v+8gCaSZZqFzSWpuD5G0x7+vjuWC5G/hNuDb6e/g+yQDCQr5AkkT4GMk9/P9Gq79eZLBGb/KjKrqdNqzqdVsN0kBDFTYDt5R0m/NNwGzImKs1P6tQMnw3msjYmbJna3tuKZgViVJ+5A0qV3XygFB0n7p3IJuSdNJmtD+tdnlsuZwUDCrgqR5JENDjyRpgmplAj5C0qz0I5Ihoh+q6wWSHFLb8/yUnNxmE8vNR2ZmNs41BTMzG9fSk9emTp0as2bNanYxzMxayr333vtURBye77WWDgqzZs1i7dq1zS6GmVlLkbS50GtuPjIzs3EOCmZmNs5BwczMxk2aPgVJp5AseNFNkjL5t5tcJDOzjtPQmoKk69PVjh7I2b5Q0jpJGyVdDhARd0fERSQLwNS8Jq6ZmVWu0c1HN5LkqR+X5kK/imSlp/kkq1fNz9rl7ZSxelW1hlcOs2bWGlZNWcWaWWsYXjlc+iAzsw7R0KAQEXex5yIhkKQt3hgRm9IMlreQ5EZH0gzgmYh4thHlGV45zLol69i5eScE7Ny8k3VL1jkwmJmlmtHRPJ09V3bawu4VtS4gWS+3IElLJK2VtHbr1q0VXXjTFZsY27Fn3rKxHWNsumJTRecxM2tXzehozrc6UrKIbUTJBU4iYrmkx4Ezenp6XlXJhXc+srOi7WZmnaYZNYUt7Lnc31GUv8RjTXpn9Fa03cys0zQjKNwDDEg6RlIPcA7Jqktli4jbI2LJIYccUtGFZy+dzZT997zlKftPYfZSr8ZnZgaNH5J6M8kauXMlbZF0QUSMAJcCd5Lkbb81Ih6s8LxnSFr+zDPPVFSe/sX9zF0+l66+rvFtU/bz/D0zs4yG9ilExKIC2+8gWQe22vPeDtx+4oknXljV8c/vXkNiZNsI65asA5KgYWbWyVrya3K1NQXwCCQzs2JaMihU26cAHoFkZlZMSwaFWmoKBUcaTcGT2Mys47VkUKilppBvBBIAo3h2s5l1vJYMCrXIjECia+/X3LdgZp2uJYNCLc1HkI4yGsv/2s7N7lsws87VkkGhluajjIJ9C3Lfgpl1rpYMCvUwe+nsglmY3IRkZp2qY4NC/+L+NA3f3tyEZGadqiWDQq19Chm9M92EZGaWrSWDQj36FKB4E9L6y9bXdG4zs1bUkkGhXoo1IY1uG3Vtwcw6TkcHBSjShAQMnjfowGBmHaXjg0LRtRQ8y9nMOkxLBoV6dTRD0oTU3Vc4g/jYjjH3L5hZx2jJoFCvjuaMgWUD+fMhpdy/YGadoiWDQr0Vy4eU4dqCmXUCB4VU/+J+5q2YV/B11xbMrBM4KGQp1b/g2oKZtbtJExQkTZG0VNI/STqvWeUYWDZQ8DXXFsys3TU0KEi6XtKTkh7I2b5Q0jpJGyVdnm4+E5gOvABsaWS5iilVW/DcBTNrZ42uKdwILMzeIKkLuAo4FZgPLJI0H5gLrImI9wMXN7hcRRWrLXjugpm1s4YGhYi4C3g6Z/NJwMaI2BQRu4BbSGoJW4BfpvuMFjqnpCWS1kpau3Xr1kYU23MXzKxjNaNPYTrwaNbzLem2rwG/L+mfgLsKHRwRy4GPAPf19PQ0rJCeu2BmnagZQSFvXtKI2BERF0TEeyPiqmInqPfktXw8d8HMOlEzgsIW4Ois50cBQ5WcoJ5pLorx3AUz6zTNCAr3AAOSjpHUA5wD3FbJCSaippBRqn/BS3eaWTtp9JDUm4E1wFxJWyRdEBEjwKXAncAgcGtEPFjheSekppBRbDSSl+40s3bS6NFHiyLiyIjYJyKOiojPp9vviIg5EXFsRCxtZBnqoWhtwUt3mlkbmTQzmisxkc1HGQPLBrx0p5m1vZYMCs3gpTvNrBO0ZFCY6D6FDC/daWbtriWDQjOaj8BLd5pZ+2vJoNAs5aS/8BBVM2tlLRkUmtV8BKXTX3iIqpm1spYMCs1qPoIy0l94iKqZtbCWDArNNp7+wkNUzazNtGRQaGbzUYaHqJpZO2rJoNDM5qNsxYaousPZzFpRSwaFyaLYEFV3OJtZK3JQqIFzIplZu3FQqFGxnEie5WxmraYlg8Jk6GjOKNbh7FnOZtZqWjIoTJaO5oxiHc6e5WxmraQlg8JkM3vpbM9yNrO24KBQB57lbGbtwkGhTjzL2czawaQJCpIWSLpb0rWSFjS7PNXwLGcza3UNDQqSrpf0pKQHcrYvlLRO0kZJl6ebA9gO7AtsaWS5GqlYp7NrC2Y22TW6pnAjsDB7g6Qu4CrgVGA+sEjSfODuiDgV+ADwkQaXq2GKzXJ2bcHMJrvCK8YAkk4o4xwvRMT9+V6IiLskzcrZfBKwMSI2pde4BTgzIn6Wvv5LoPDX7Umuf3E/Gy7bwMi2kbyvb7piU9LMZGY2CRUNCsB/A/eQv/s04xhgVgXXnA48mvV8C/AaSW8Bfh84FLiy0MGSlgBLAGbMmFHBZSfOwLIBBt8xmPc1D081s8msVFC4JyLeWGwHSf9V4TXzjs+JiK8BXyt1cEQsl/Q4cEZPT8+rKrz2hChaW0iHp7q2YGaTUdE+hVIBodx9cmwBjs56fhQwVOE5Jj3nRDKzVlRWR7OkhaX3Kts9wICkYyT1AOcAt1VygsmW5iIf50Qys1ZUNChIOikdLfSxrG1fLPfkkm4G1gBzJW2RdEFEjACXAncCg8CtEfFgJYWeTAnxiimVE8lDVM1ssilVU1hE0tk8W9InJL0dKGdEEgARsSgijoyIfSLiqIj4fLr9joiYExHHRsTSSgvdCjUFKJ0TyUNUzWyyKRUUPhARvwM8AnwDOAw4QtL3JX2p4aUroFVqCiVzIuEJbWY2uZQKCndK+g5wONAH3AFsjojXAn/e6MIV0io1BcjKiVSAawtmNpmUGn30u8AfkqSfOBb4KHCcpH8j6SBuilapKWQUXbYT1xbMbPIoOfooIp4HHo2IT0XE24GfAxcCv2h04YqUqWVqChkDywYKvubagplNFmUNSc2Zi3BlRGyNiK82qExtqVRtwXMXzGwyqDghXmYEUTO1WvNRRrHagucumNlkUGqewvJSJyhnn3prxeYjKF1b8NwFM2u2UrmPzpL0myKvC/jdOpan7Q0sG2DdknWM7RjL+3qmf8G5kcysGUoFhb8s4xx316MgnSLzYT943iCM5t9n/WXrHRTMrCmKBoWIWDFRBamEpDOAM4477rhmF6Uq44GhQHpt1xbMrFkmzRrNlWjVPoVsnrtgZpNRSwaFdlFq7sLqqas9GsnMJlTJoCCpS9LfT0RhOk2p2sLIthEPUzWzCVXOjOZR4FWSii3JaVUqOncBD1M1s4lVbvPRj4B/l3SupLdkfhpZsGJadfJaPqVqC5A0Ja2/xIHBzBqv3KBwGLANeCNwRvpzeqMKVUo7dDRnG1g2UHTdBYCha4fcjGRmDVdqngIAEXF+owvSyTJDT9dftp7RbQUmL4TnL5hZ45W7RvMcSf8p6YH0+csk/U1ji9ZZ+hf3c8pTpxRtSnIzkpk1WrnNR58DPgi8ABARP6WJ6ym0s4FlA0nykALcjGRmjVRuUNg/In6Ys22k3oWRdICkeyU1rb+i2foX9zPtommFdwin2Tazxik3KDwl6VggACSdDTxe6iBJ10t6MtPslLV9oaR1kjZKujzrpQ8At5ZZprY15+o5xUckjcLguYNuSjKzuis3KPwp8M/A8ZIeA94HXFTGcTcCC7M3SOoCrgJOBeYDiyTNl/Qm4GeAvwJTuhmJgKFrhhwYzKyuyh19tAl4k6QDgCkR8VyZx90laVbO5pOAjek5kXQLcCZwIHAASaB4XtIdEbFXfmlJS4AlADNmzCinGC2pf3E/z3zvGYauHUrrZ/kNXTMEJLULM7NalTv6qE/SZ0nSZK+StExSX5XXnA48mvV8CzA9Iq6IiPcB/wJ8Ll9AAIiI5cBHgPt6enqqLEJrmHP1HOZ9cR50Fd9v6Joh50kys7oot/noFmAr8MfA2enjL1V5zXyNIuPfhSPixoj4erETtNvktWL6F/czb8W84k1JOE+SmdVH2TOaI+LvIuIX6c9HgUOrvOYW4Ois50cBQ5WcoJ3SXJSj5IiklPMkmVmtyg0K35V0jqQp6c9bgW9Uec17gAFJx0jqIZnvcFuV5+oYc66ew7SLSwcGT3Azs1ooonAvpqTnSJp2RNIJnGnnnwJsj4iDi55cuhlYAEwlGVX0fyPi85JOA/6RpLX8+ohYWk3hTzzxxFi7dm01h7as9ZesL9n5DDDt4mnufDazvCTdGxEn5n2tWFCYrLKW47xww4YNzS7OhBteOVw8T1LKgcHM8qlLUJD0MmAWWcNYI+Jr9ShgtTqxppBt9dTVjGwrPrG8u6+bgWUDTqRnZuOKBYVyh6ReD1xPMvqo6amzLVFyghvJqCTPfjazcpU1eQ14bUTMb2hJKpDVfNTsojTV+AS3a0oM3gpPcjOz8pQ7+miNpEkTFDppnkIp5Y5KAqfFMLPSyg0KK0gCwzpJP5V0v6SfNrJgxXTaPIVSxgNDGatoe/azmRVTVkezpI3A+4H72T0slYjY3LiildbpHc25yh2VBIBg2kUenWTWiWruaAYeiYjb0tnMmzM/dSyj1UFm9baympOcZdXM8ii3o/khSf8C3A7szGxs9pBUyy/z7b9kBzTugDazPZUbFPYjCQZvztoWQFOCgkcflTYeGMqY/Tx0zRBD1wx5ToOZteaM5gz3KZRWUT8DoB5x/PXHOzCYtbGaZzRLuoE83zcj4t21F696DgrlW3/J+rKakwC6+ro45alTGlwiM2uWenQ0f50kK+o3gP8EDga216d4NhEqmc/gTKtmnavc5Ti/mv08zX76Hw0pkTVMNf0MvTN7mb10tpuTzDpEuTWFXANA0xZI9uS16mWW+OzqK7HGZ2rn5p3OnWTWQcpNiPecpGczPyRDUz/Q2KIV5jQXtcnMZ5h307zy/gI8p8GsY5TbfHRQowtiEy/TJDR47mDJ5iTwnAazTlA0KEgq2kQUEY/Utzg20crOtJrynAaz9laqpvANdi/HmRHA4cCLSZbTtBZXSQd0xsi2EQbfMcgz33vGNQezNlK0RTkiXhoRL0v/fSnJ4jrfIxmO+r56FkTSPEnXSvqKpIvreW4rLdMB3Tuzt6LjnHXVrL2U29E8IOlG4JvAvcD8iPinMo67XtKTkh7I2b4wTcO9UdLlABExGBEXAW8F8k6qsMbqX9zPyQ+fzIJYUPacBkhqDeuWrHNgMGsDRYOCpJekcxK+SjIv4SURcV1EvFDm+W8EFuacswu4CjgVmA8syizgI+kPgdUkE+SsiSpZowFgbMcYg+cNOjCYtbhSNYWfACcDdwMnAZ+R9NnMT6mTR8RdwNM5m08CNkbEpojYBdwCnJnuf1tE/DawuML7sAaodE4DozD4jkHuPuhuBwezFlWqo7kRuY2mA49mPd8CvEbSAuAtQC9wR6GDJS0BlgDMmNG0+XMdo39xP/2L+xleOcxD73mI+HXpnujR7aPuhDZrUUWDQkSsaMA18zVIRESsAlaVOjgilkt6HDijp6fnVXUumxWQCQ6rp65mZNtIWccMXTPEk7c+6aGrZi2k2jQXtdgCHJ31/CigvEHy1nQDywaYsn/5fzaZoatuUjJrDc0ICvcAA5KOkdQDnAPcVskJnOaiefoX9zN3+dyKh65mmpRWaRVrZq1xgDCbpModkvq6crbl2edmYA0wV9IWSRdExAhwKXAnMAjcGhEPVlJoJ8RrrmqHrmbs3LzTQ1jNJqlyF9m5LyJOKLVtonmRncmhkk7oPUyBeV+Y5/4GswlW9cprkk4Gfptk9vJnsl46GPijiHh5HctZtqw1mi/csGFDM4pgeVS69GdG14FdzLl2joOD2QSpZeW1HuBAklFKB2X9PAucXc9CVsJ9CpNTJiV3pU1Ko9tHeejdD7k5yWwSKLf5aGZEbJ6A8pTFNYXJr+omJXAGVrMGq7r5KOsEc4C/AGaRNbchIt5YpzJWxX0Kk9/wymE2XbGJnZt3VnzstIunefKbWQPUIyj8BLiWJBneeINxRNxbr0JWw0GhdQyvHGbw/EEoN2tWyrUGs/orFhTKWnkNGImIa+pYpppkNR81uyhWpsyHeqVNSpkMrNnnMLPGKXfy2u2SLpF0pKTDMj8NLVkR7mhuTf2L+3nD9jcw76YKkuyRZmB9pzOwmk2EcpuPfpFnc0TE7PoXqXxuPmpt1TQpdR3YxeivR+md0cvspbNdezCrQs3NRxFxTH2LVBs3H7WHzAd6JXMbRrcn++3cvJPBc52J1azeyq0p7A+8H5gREUskDQBzI+LrjS5gMa4ptJdaJr+59mBWvlomr2XcAOwimd0MSabTj9ahbGbjMpPfuvvKHf+QGN0+CuGcSmb1UG5QODYiPkna+hsRz1P2Qo1mlRlYNgD7VHfs2I4x1l+2vr4FMusg5QaFXZL2AwJA0rFA5bORzMrQv7ifeTfMQwdU971jdNso6y9xYDCrRrl9Cm8GrgDmA98GXge8K10tbcI5zUXnWH/JeoauHUq/jlShCxiF3pnubzDLqHlGc3qSPuC1JM1G34+Ip+pXxOq4o7kzjKfKeGQn2l9V5VMCQDDtIqfOMKtHmovbgJuB2yLi13UuX9UcFDpTzbUHXHOwzlaP0UefBk4Bfibpy5LOlrRv3UpoVoE5V89h3hfnVbwkaLadm3d6eVCzPMoKChHx3xFxCTAbWA68FXiykQUzK6bWJUGzeSir2W7l1hRIRx/9MXAR8GpgRaMKZVaJOVfPSQJDDYOkx3aMMfiOQdcarOOVFRQkfQkYBN4IXEUyb+G99S6MpLMkfU7Sv6cjnszKUo8mJXCtwazcjuaFwHciorL8A8mx1wOnA09GxEtyzrmMZNDgdRHx8azXXgR8KiIuKHZudzRbIbUs7jMuZzgrMD4Kyik1rJVV3dEs6a8AIuJbwFtyXvtYmde/EViYc2wXSY3jVJK5D4skzc/a5W/S182qkulzqKnmkH4FynRKD75jMAkyTqlhbaxU89E5WY8/mPPaQsoQEXcBT+dsPgnYGBGbImIXcAtwphKfAL4ZEfflO5+kJZLWSlq7devWcopgHWz20tlM2b/srrOKjO0YY9MVmxpybrNmKfW/RQUe53teienAo1nPt6Tb3gu8CThb0kX5DoyI5cBHgPt6enpqKIJ1gv7F/cxdPjepMQi6+rpQT/3Sdu18xNlerL2UCgpR4HG+55XI978yIuKzEfGqiLgoIq4tWCivvGYVGB++OraAU546heOvP353s1Kt8SHwiCVrK6WCwsslPSvpOeBl6ePM85fWcN0twNFZz48Chso9WNIZkpY/88wzNRTBOlX2HId5X6xsadB8Mn0Oq6eudnCwllc0cX1E1Pa/pbB7gAFJxwCPkfRdvL1B1zIrqH9xP/2L++uSOmNk2wiD5yYd0r0ze+k7rY9td2zzaCVrKY3pgcsi6WZgDTBX0hZJF0TECHApcCfJ/IdbI+LBcs/p5iOrtz3mOdTS95AGlZ2bdzJ0zZBHK1nLKTtL6mTi1Nk2EbKzs/bO6GV0+ygj20ZqO6kYDxzdfd0MLBtw7cEmXF1SZ09GnrxmE2l45TCD5w7WNsQij2kXO523Tax6ZEk163j9i/uZdlFtOZbyGbpmiNVTV7NqijO2WvO1ZFDw6CNrlty+h3r9DxrZNuK+B5sU3HxkVoPhlcMMnj8IL9T3vL0zezn54ZP36tfwCCarh7brU3BHs00mwyuHWX/Zeka3VZwvsqiuvq69z7kPdB/czcjTIw4SVrW2CwoZrinYZNOoAFHIlP2nMHf5XAcGq4g7ms0mSP/ifk556hTm3VRgprRA+9avp7pQUr7hlcOsmbUmb+d1sdfMXFMwa6B8fQKNGNbaO7N3/Bp9p/XxxIonGNsxtsc+3X3dvPitL97rNdc2Ok/bNR+5T8Fa2ZpZa2pb/KcBMh3b1hnarvnIaS6sleVd46HOcx8q5RTgltGSQcGsleWu8dA7szeZFLdPEws1JWnqcn+DtWTzUYb7FKxd5Bu1pANEPB8wVuTAetoHJBG7dn8mTFR/g+djTKy2az4yayfDK4dZt2TdHgFhyv5TOP6fj2feF+Y1bDnRvbzAHgEBJmbJ0cz9O6Ps5NCSQcFpLqydbLpi014jhTIfxvmamg79vUMntHyN7m8odv828dx8ZNZkq6asyj9EVbBgbEHeY1ZPXV17Gu9yCboPS2ZRdx3WhVBdZ1RXc/9WGzcfmU1ivTN6K9oOMPJ0eQGh68B0Al0to5tid8K+0W2jVSXvK9aBXc39W+M4KJg1Wb4hqlP2n8LspbMLHlPwA7OL8WameTfN45TnThlfi5oGLK47tmOMwfMGiwaGfH0Gg+cOskpJgOg7ra/i+7fGcVAwa7J8/QalRvwUCiTzVsxjwdgCTn745D2O71/c37hRTKMUrDEMrxxm8LzBvfoMcpctDQXdfd1l3781zqTpU5A0G7gCOCQizi7nGPcpWCerdBjnRMyk7p25uxyZGsJeAaEIL1E6MYr1KXQ3+MLXA6cDT0bES7K2LwSWkVRor4uIj0fEJuACSV9pZJnM2kX/4v6KPjxnL51d8Yd0pTL9DJB/VFEpI9tGxo93YGiORjcf3QgszN4gqQu4CjgVmA8skjS/weUw63j5mqkaYWzHGIPvGKy6VlLJcFTPwK6/htYUIuIuSbNyNp8EbExrBki6BTgT+Fk555S0BFgCMGPGjPoV1qwD5NYuCjYpdQETsyREXjs372SVVo1ndt12x7a9mslym6eyaynl1jLyNcEB49saMQR3smtGR/N04NGs51uA6ZL6JF0LvFLSBwsdHBHLI+LEiDjx8MMPb3RZzdpaoQ7raUuanIspNbJthKFrhvLOdq510lveUVHnD/LQux8a31btENxW1oygkG/EdETEtoi4KCKOjYj/V/QEntFsVheFRj7NuXoO3Qc3tCGhaplhsIWap8qdgZ23zyNPqo/ca7f7TOtmvOtbgKOznh8FDDWhHGZG4Q7rcifINUWRpq1yJ71Vm76j3dOMN6OmcA8wIOkYST3AOcBtlZzA6ymYNd6Ezyiuw6dRvklvwyuHuXvq3azSKlZpFaunrmZ45XDV99fuM60bGhQk3QysAeZK2iLpgogYAS4F7gQGgVsj4sEKz+vmI7MGy7sY0D6gngI5M/YBeqq82D7UvNBQd1/3XpPehlcOM3j+4B4ZaEe2jfDQux/KO5O66P3RGTOtJ83ktWp48ppZYxUdnbN55/gopa6+3aN0xkfsbBvZPYopd7+s13pn9jK6fbTmBH+ZJUXHy1xqSGwXTFsyLRnZVOpeijxuxVFJXqPZzBom38zlShfnKZgptRKCeV+cV9EEvSn7T+GI847giRVPFC1/qdnZE7UYUb20XVDIcE3BrPkKzXXIfHOv5RyVyEzGq/g8BeZkZJe/nPJVcr/N1naps92nYDZ5FBqNU8konbz9FxXa9cSu6gJLgZFMOzfvHJ+TUM69tMuopJYMCh59ZDZ51GM9hHzzJaZdPI2uvqx835lPqwIpwGNn/Vs9Hnr3Q2WPVGqXUUmTc3aKmbWMfIn2qhmlk2++xJyr5+Tdt6rmpn2AFyo7JHYFm67YVDKZYDuNSmrJoJDV0dzsoph1vMwHeSVpvGtVaVNNJqX34DsGq76W9hPsSLbpANG1b9f4KKqxHWOsv2w9Gy7bMD4iqe+0vr1GNmWnFq9UpanSq+WOZjNrOWXXFNIRSZkPz2pqGF19XcTzsVdNKN+opXJUM1KpHiO8srVdR7OZdbayO6aDPXIVVdqhrR4hlDfx3tDyoarWpqgmf1Ktyf8q0ZJBwaOPzDpbbsf0Hh3SObKbmjLHlaO7r5vjrz++cA6oGlKLV9r8VY8RXuVqyaDg0Udm1r+4n5MfPpkFYws45alTCi4alDsqqH9xf+F9Z/ayIBawIBbwO0/9TrJvoVFFheNQSZWOVKrHCK9ytWRQMDPLVShX0+j20b1WZus7rW+vXEuFRhAV2vfQBYdWna9pdPtowXUZ8q0mV2jdi0aMeHJQMLO2kK9JSdJei+Ssv2Q9T6x4Ys+0GoIjzjtir07b4ZXDefc9+OSDeXbNs4VTcwj2m79fwaCRWYs6NzDkW/gns5pcvnUvPPooh0cfmVkhlS41mi9NRU3LlZaxT+4165EypBxtN/rIHc1mVkrBTthCaS3y7F/pOSrdJ/f8E9mhXEhLBgV3NJtZKZV2EOfbv6ZO5jL2yT3/RHYoF9KSQcHMrJRCnbPTlkwru9O2knNUs0/uNSeyQ7lguSbsSmZmEyhfkr25y+cy5+o5ZXfalnuOrr4uuvu6K94n95qFrjeR6zS4o9nMrMO0XUezmZk1xqTJkirpAOBqYBewKiJWNrlIZtbBCmUlLZatdHjlMOsvW8/otmToUXdfNwPLBoDdWWS7Duti7DdjxK/TVpopwFjSVDSeWTXfmtgTlIG2oc1Hkq4HTgeejIiXZG1fCCwj6Z+/LiI+Lulc4FcRcbukL0XE20qd381HZtYIhbKSFlvPGWDw/MG912zoAnWJ2FXFZ+0+IO15bD3Wg25m89GNwMKcwnQBVwGnAvOBRZLmA0cBj6a71ZBqysysNoWykubLjJrJVrrpik35F/EZpbqAAPDC3sc2KjtqRkODQkTcBTyds/kkYGNEbIqIXcAtwJnAFpLAULRckpZIWitp7datWxtRbDPrcNVMfJvICWaNvFYzOpqns7tGAEkwmA58DfhjSdcAtxc6OCKWAx8B7uvp6WlkOc2sQ1Uz8W0iJ5g18lrNCAr5UkRFRPw6Is6PiItLdTJ7RrOZNVI1E99mL52drAOdqytZrKcq++x9bKMnszUjKGwBjs56fhQwVMkJnPvIzBqpmolv/Yv7mXfDvD0W/Onu62beinkcf/3xe2ZvPSDrgz79FO6d2cu0i6ftce55N+x57ERMZmv45DVJs4CvZ0YfSeoG1gO/BzwG3AO8PSIerOCcZwBnHHfccRdu2LCh/oU2M2tjTRt9JOlmYA0wV9IWSRdExAhwKXAnMAjcWklAADcfmZk1SkMnr0XEogLb7wDuqPa8WTWFak9hZmZ5tGSaC9cUzMwaoyWDgjuazcwaoyWDgmsKZmaN0dKpsyVtBTZXefhU4Kk6FqdVdOJ9+547RyfedzX3PDMiDs/3QksHhVpIWltoSFY768T79j13jk6873rfc0s2H5mZWWM4KJiZ2bhODgrLm12AJunE+/Y9d45OvO+63nPH9imYmdneOrmmYGZmORwUzMxsXEcGBUkLJa2TtFHS5c0uT6NIeljS/ZJ+LGltuu0wSd+RtCH990XNLmetJF0v6UlJD2RtK3ifkj6YvvfrJP1+c0pdmwL3/GFJj6Xv948lnZb1Wjvc89GSvitpUNKDki5Lt7fte13knhv3XkdER/2QrJ30c2A20AP8BJjf7HI16F4fBqbmbPskcHn6+HLgE80uZx3u8/XACcADpe6TZF3wnwC9wDHp30JXs++hTvf8YeAv8uzbLvd8JHBC+vggkhT889v5vS5yzw17rzuxplBojehOcSawIn28AjireUWpj8i/Fnih+zwTuCUidkbEL4CNJH8TLaXAPRfSLvf8eETclz5+jiT1/nTa+L0ucs+F1HzPnRgUCq0R3Y4C+LakeyUtSbf1R8TjkPzBAS9uWukaq9B9tvv7f6mkn6bNS5lmlLa753TxrlcCP6BD3uuce4YGvdedGBTyrhE94aWYGK+LiBOAU4E/lfT6ZhdoEmjn9/8a4FjgFcDjwKfT7W11z5IOBL4KvC8ini22a55tLXnfee65Ye91JwaFmteIbhURMZT++yTwryTVyGFJRwKk/z7ZvBI2VKH7bNv3PyKGI2I0IsaAz7G72aBt7lnSPiQfjisj4mvp5rZ+r/PdcyPf604MCvcAA5KOkdQDnAPc1uQy1Z2kAyQdlHkMvBl4gORez0t3Ow/49+aUsOEK3edtwDmSeiUdAwwAP2xC+eou88GY+iOS9xva5J4lCfg8MBgR/5D1Utu+14XuuaHvdbN715vUo38aSS/+z4Erml2eBt3jbJJRCD8BHszcJ9AH/CewIf33sGaXtQ73ejNJFfoFkm9KFxS7T+CK9L1fB5za7PLX8Z6/CNwP/DT9cDiyze75d0iaQn4K/Dj9Oa2d3+si99yw99ppLszMbFwnNh+ZmVkBDgpmZjbOQcHMzMY5KJiZ2TgHBTMzG+egYE0haTTN7viApC9L2r8O5/xbSW8qsc9Fkt6ZPn6XpGllnHeP/SRdJ2l+HcqbyXT5t7Weq8g1viVpes62BZK+XuX53pZm4KzqeJv8HBSsWZ6PiFdExEuAXcBF2S9K6qr0hBHxoYj4jxL7XBsRX0ifvgsoGRRy94uIP4mIn1VavgI+ExEfqtO59iBpP5Ix+4/V65wR8SXgT+p1Ppt8HBRsMrgbOC79BvtdSf8C3C+pS9LfS7onTfz1nswBkv5KyVoRP5H08XTbjZLOTh8/LOkTkn6Y/hyXbv+wpL9I9zsRWJnWWPaT9KH0Wg9IWq5Evv1WSToxPd+itBwPSPpEVvm2S1qalu/7kvpL/RLSsq2Q9O20/G+R9Mn0/N9K0x1k7u1jktZIWivpBEl3Svq5pOzgugBYlR6zUNJDklYDb8m65kmS/kfSj9J/56bb75b0iqz9vifpZeW/pdaqHBSsqSR1kyTsuz/ddBLJ7Ov5JLN0n4mIVwOvBi5M05OcSpIe+TUR8XKSfPr5PBsRJwFXAv+Y/UJEfAVYCyxOayzPA1dGxKvT2st+wOkF9suUfRrwCeCNJInJXi3prPTlA4Dvp+W7C7iwzF/JscAfkKRAvgn4bkS8FHg+3Z7xaEScTBJQbwTOBl4LZDdFnQp8S9K+JPlxzgBOAY7I2uch4PUR8UrgQ8DH0u3XkdSQkDQH6I2In5Z5D9bCHBSsWfaT9GOSD9xHSPK7APwwkjzwkORreme63w9I0hkMAG8CboiIHQARUWhdgZuz/j25jDL9rqQfSLqf5IP+t0rs/2pgVURsjYgRYCXJ4jeQNIll2t3vBWaVcX2Ab0bECyRBsgv4Vrr9/pxz3Ja1/QcR8VxEbAV+I+nQ9LXXAauB44FfRMSGSFIY3JR1nkOALytZwe0z7L7nLwOnp7WTd5MEHusA3c0ugHWs5yPiFdkbktxf/Dp7E/DeiLgzZ7+FlJcOOAo83kv6bfpq4MSIeFTSh4F9S5w/X5rijBdidw6ZUcr/v7YTICLGJGWfYyznHDuztu/M2j4GdEuaTVKb2JX+Xgvd/9+R1Eb+SEm+/lXp9XdI+g5JjeWtJE1o1gFcU7DJ7E7g4qy29DlKMr5+G3i30hFLkg4rcPzbsv5dk+f150iWOITdAeApJbnrzy6wX7YfAG+QNDXtGF8E/HdZd9Z4p7K7lvEQcIykY9Pni7L2OwTIdES/K+cc1wGfBe4pUhuzNuOagk1m15E0mdyn5OvuVuCsiPhW2gm6VtIu4A7gr/Mc3yvpByRffhblef1G4FpJz5M0L32OpDnmYZIU64X2A5JVviR9EPguSa3hjoiYLKnIFwLvBYiI3yhZee8bkp4iaVJ6SbrfJ4EVkt4P/Ff2CSLiXknPAjdMXLGt2Zwl1dqSpIdJmoKeanZZCkmbqLZHxKfqfN5e4HsRUVOTT9qRvgo4PpLFXDLbF5AsGn96Lee3ycnNR2bNsx1YojpPXotk0fZaA8I7SZrHrsgJCG8j6Xv5ZW2ltMnKNQUzMxvnmoKZmY1zUDAzs3EOCmZmNs5BwczMxjkomJnZuP8PMRMjm2dKUcIAAAAASUVORK5CYII=\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_eraint')\n",
    "plt.savefig(\"/home/mkibria/Academic/pr_of_thesis/figures/pdf_regcm_eraint.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
}