pdf_regcm_noresm.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:36:42 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_name:                         historical\n",
       "    model_id:                                        ICTP-RegCM4-7\n",
       "    driving_experiment:                              NCC-NorESM1-M, historica...\n",
       "    driving_model_id:                                NCC-NorESM1-M\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-32de1a66-0376-4a44-954d-390e992a4dc4' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-32de1a66-0376-4a44-954d-390e992a4dc4' 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-6843facd-3bc7-4c01-8869-492a11103b6a' class='xr-section-summary-in' type='checkbox'  checked><label for='section-6843facd-3bc7-4c01-8869-492a11103b6a' 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-f9dabfec-19a9-4855-b2da-b85fa7048f27' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-f9dabfec-19a9-4855-b2da-b85fa7048f27' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c3e91ffc-f462-422b-8386-ef46c77f723f' class='xr-var-data-in' type='checkbox'><label for='data-c3e91ffc-f462-422b-8386-ef46c77f723f' 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-73227030-438b-485d-8e80-3d50d8801910' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-73227030-438b-485d-8e80-3d50d8801910' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-00c8c743-284d-452f-88c1-49ee9abf2106' class='xr-var-data-in' type='checkbox'><label for='data-00c8c743-284d-452f-88c1-49ee9abf2106' 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-b80b81d6-467e-47bf-83c7-5602b4654113' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-b80b81d6-467e-47bf-83c7-5602b4654113' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-1ea517d9-7a92-48c8-b1fb-99a06b6ec996' class='xr-var-data-in' type='checkbox'><label for='data-1ea517d9-7a92-48c8-b1fb-99a06b6ec996' 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-865359b8-407e-49e2-8d3a-8e50b1c48baf' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-865359b8-407e-49e2-8d3a-8e50b1c48baf' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-3989da18-6922-41dc-8937-38f7a684f3f2' class='xr-var-data-in' type='checkbox'><label for='data-3989da18-6922-41dc-8937-38f7a684f3f2' 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-8acd5445-f76e-4d42-ad81-ae8c069ccf9d' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-8acd5445-f76e-4d42-ad81-ae8c069ccf9d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d8740705-3b7a-49ae-b0ce-618d49c0f80c' class='xr-var-data-in' type='checkbox'><label for='data-d8740705-3b7a-49ae-b0ce-618d49c0f80c' 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-55b20c55-8024-4112-aaeb-f47e546022ba' class='xr-section-summary-in' type='checkbox'  checked><label for='section-55b20c55-8024-4112-aaeb-f47e546022ba' 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-0a322e44-b23d-41ab-be27-4d00fec9a5c0' class='xr-var-attrs-in' type='checkbox' disabled><label for='attrs-0a322e44-b23d-41ab-be27-4d00fec9a5c0' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-038a1295-89e5-4045-af1a-54078dc3443e' class='xr-var-data-in' type='checkbox'><label for='data-038a1295-89e5-4045-af1a-54078dc3443e' 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-16fe40fd-acc3-475d-ae42-c5a221b07c71' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-16fe40fd-acc3-475d-ae42-c5a221b07c71' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6160d961-8e2b-4e16-ad05-d8894130bc12' class='xr-var-data-in' type='checkbox'><label for='data-6160d961-8e2b-4e16-ad05-d8894130bc12' 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-cdf54b32-5c63-4827-8cbb-13bc0bb8cfc1' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-cdf54b32-5c63-4827-8cbb-13bc0bb8cfc1' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-f0c015e5-5a82-4eca-9dfb-9cabd8555f8d' class='xr-var-data-in' type='checkbox'><label for='data-f0c015e5-5a82-4eca-9dfb-9cabd8555f8d' 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-eeee7704-487b-4946-9c7b-acdc7fa66684' class='xr-section-summary-in' type='checkbox'  ><label for='section-eeee7704-487b-4946-9c7b-acdc7fa66684' 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:36:42 2022: cdo mergetime pr_WAS-22_NCC-NorESM1-M_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_197001-198012.nc pr_WAS-22_NCC-NorESM1-M_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_198101-199012.nc pr_WAS-22_NCC-NorESM1-M_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_199101-200012.nc pr_WAS-22_NCC-NorESM1-M_historical_r1i1p1_ICTP-RegCM4-7_v0_mon_200101-200512.nc noresm_pr_197001-200512.nc\n",
       "2018-03-21 03:01:05 : 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-19T16:13:52Z</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>e18baf90-8c3e-11e9-a8fd-f04da273d2fb</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>NO_85</dd><dt><span>model_sst_data_source :</span></dt><dd>NO_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> 2005-12-31 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_name :</span></dt><dd>historical</dd><dt><span>model_id :</span></dt><dd>ICTP-RegCM4-7</dd><dt><span>driving_experiment :</span></dt><dd>NCC-NorESM1-M, historical, r1i1p1</dd><dt><span>driving_model_id :</span></dt><dd>NCC-NorESM1-M</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:36:42 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_name:                         historical\n",
       "    model_id:                                        ICTP-RegCM4-7\n",
       "    driving_experiment:                              NCC-NorESM1-M, historica...\n",
       "    driving_model_id:                                NCC-NorESM1-M\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/NorESM1-M/noresm_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/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAq+klEQVR4nO3debxcdX3/8dc7NyTsoQQ6bcISQhaSuqC9gGiRiNSG1ghFHkK8ICASweLS/lrF4sNiBX/Vai2WrZElKBFE9KegCBZqBDQiiUUBLwlp2K7BEAIBAjEhN5/fH+fMZTLMnTkzd+bO9n4+HvPIzFk/M+dmPvNdjyICMzMzgDHNDsDMzFqHk4KZmQ1xUjAzsyFOCmZmNsRJwczMhjgpmJnZECeFLiWpT9KP6r1tFeffKGlqPY85zHkWSbog47Y/lHRq+vw0SXc3Njqz1uOk0IIkPSppU/rFuVbS1ZJ2rec5ImJxRLyjlm0lhaRpWc8laYmkDxQdc9eIWJ094saLiGMi4pqRHif9fF5Mr99vJf2bpJ50Xf7aviBpg6SfSTpL0piC/RdJ2pLun3+cONK4zLJwUmhd8yJiV+CNwCHAp4o3kDR21KOyrF6fXr+3A+8FzixYNy8idgP2B/4F+ARwZdH+X0gTZ/7xzSwn7aa/iXyitfpyUmhxEfFb4IfAa2DoV+jfSHoYeDhd9k5J9xX88nxdfn9J+0r6jqR1ktZLujhdvl31SHrcj0haLelpSf+a//VauK2kO9NdfpX/BSvpDyR9Pz3Hs+nzfdLtLwSOAC5Ot7+44HzT0ucTJH0t3f8xSZ8qPrekL6bHfkTSMcN9XpLeIOmX6S/xbwI7FqwbNs50/atKNOnySyR9qWjZzZI+Vv7qQUQ8BNxFev2K1j0XETcBJwKnSnrVNpVIOl/SjZKulfQ8cFr6eV4p6cm0pHJBQUmlR9KX0mv8iKRz0msxNl2/Z1oyXZN+Rt9Nl8+RNCDp45KeSo99nKS/lLRS0jOS/jFjvDek1/sFSQ9K6i1YPyu9DhvSde8qWLdI0mWSbpH0IvA2SZMkfTu9po9I+kjB9odKWibpeSUl7n9Ll09J3/Ppkp5I3+dZkg6R9Ov03BdXey06RkT40WIP4FHg6PT5vsCDwGfT1wH8F7AnsBNJSeIp4DCgBzg13X98+vpXwJeBXUi+IP8sPc5pwN0F5wzgx+lx9wNWAh8os+20gtcTgXcDOwO7Ad8Cvluwfkn+WKWOAXwN+F6675T03GcUnPtlkl/aPcDZwBpAJT63ccBjwN8COwAnpPteUG2che8ZODQ955j09V7AS0BumOtX+N5mA78reD9D17Zon8eBs9Pni/IxZ/hbOT99j8eR/MjbCfgu8J/pNf9D4BfAB9PtzwJ+A+wD/AFwexrv2HT9D4Bvput2AI5Ml88BtgKfTpefCawDvpF+ln8C/B6YmiHe3wN/mV7P/wv8PF23A7AK+Mf0Wh4FvADMLPhcngPekr7XnYHlaUzjgKnAauAv0u2XAqekz3cF3pQ+n5K+58tJ/k+8I43pu+nnNZnk/9SRzf4uaMr3T7MD8KPERUm+ODYCG0i+5C4FdkrXBXBUwbaXkSaMgmUrgCOBw9P/uGNLnOM0Xv1FP7fg9YeAO8psO61M/AcDzxa8XsIwSSH9YtgMzC5Y90FgScG5VxWs2znd949KnPetFCUM4GcM8wVbLs4S77kf+PP0+TnALWXefwDPA88C/wtcwCsJ5VFKJ4WfA+elzxelX1Ib0sfTZc51PnBnwetc+nnuVLBsPvDj9Pl/kyaI9PXRabxjgT8GtgF/UOI8c4BNQE/6erd0v8MKtlkOHFfhb/t84PaC17OBTenzI0gS6JiC9dcB5xd8Ll8rWHcY8HjR8T8JXJ0+vxP4DLBX0TZT0tgnFyxbD5xY8PrbwMdG+n+5HR9dU//Yho6LiNuHWfdEwfP9SaoePlywbBwwCRgEHouIrRnPWXjcx9JjVCRpZ5LSyFySX5gAu0nqiYjBCrvvxSu/8AvPPbng9e/yTyLiJUmQ/PIrNgn4baT/qwuOVY84rwFOJimlnQxcVGH7N0bEqgrbFJoMPFPw+osRsV07kqQ+khIAwF0Rka9GK/572AF4Mv2cIPlVnd9mUtH2hc/3BZ6JiGeHiXF9wee0Kf13bcH6TZS+LsV+V/D8JWDHtPpqEvBERGwrWF/8t1D8XidJ2lCwrIekug7gDOCfgYckPQJ8JiK+X7Btcey1vJeO46TQngq/9J4ALoyIC4s3knQ4sJ+ksRkTQ76qCpIqpDUZ4/k/wEySX42/k3Qw8D9A/lup3FS8T5NUf+xPUq2RP/dvM5670JPAZEkqSAz7kfxazxJnOdcCD0h6PTCLpKqhLiQdQvLFV7YLbEQsBhaXWlXw/AmSksJew1zzJ0mqjvL2Ldp3T0l7RMSGDKHX2xpgX0ljChJDviozr/i9PhIR00sdLCIeBuan7VPHAzdKmtiAuDuKG5rb31eBsyQdpsQukv5K0m4kdclPAv+SLt9R0lvKHOsf0sbYfYGPktQtl7KWpP42bzeSX1YbJO0J/FOF7YekvzxvAC6UtJuk/YG/I/kSrtZSknrvj0gaK+l4kvaArHEOKyIGgHuBrwPfjohNFXapSNLukt4JXA9cGxH3j/SYEfEk8CPgS+nxx0g6UNKR6SY3AB+VNFnSHiQ9nwr3/SFwafp3sIOkt440pircA7wIfDw99xxgHsnnU8ovgOclfULSTmkj+mvSJIukkyXtnSaYDek+lUqEXc9Joc1FxDKSRr+LSeqwV5HUh+e/cOeR1N0/DgyQ9HQZzvdI6oXvI2lwLO4mmXc+cE3aS+M9wL+TNHA+TVI3fmvR9hcBJ6S9PL5S4ngfJvkyWE3ya/kbwFVl4iwpIraQ/CI8jeSzOBH4TsEmleKs5BrgtSSJYSRulvQCyS/d84B/A04f4TELvY+kSu43JJ/DjSTtBZD8iPgR8GuSUtItJIk0/2V5CknJ7SGSxtaP1TGustLr9y7gGJJrdCnwvkh6cJXaPv/3fTDwSLrPFcCEdJO5wIOSNpL8DZ4UEb9v5HvoBNq++tW6laQApldZD95V0l/N1wJTiuq925aS7r2XR8T+zY7FWoNLCmYZSNqBpErtinZOCGk1y1+m1WuTSarQ/l+z47LW4aRgVoGkWSR10n9MUgXVzkTSTfNZkuqjfpJ+/vU7QTKH1MYSj4qD26z5XH1kZmZDXFIwM7MhTgpmZjakrQev7bXXXjFlypRmh2Fm1laWL1/+dETsXWpdWyeFKVOmsGzZsmaHYWbWViQ9Ntw6Vx+ZmdmQlikpSDoC6COJaXZEvLnJIZmZdZ2GlhQkXaXkhhwPFC2fK2mFpFWSzgWIiLsi4izg+yTTCZiZ2ShrdPXRIpL5R4YouQPUJSTzm8wmmcVwdsEm7yWZQ93MzEZZQ5NCRNzJ9nPEQzJr5aqIWJ1OgHU9cCyApP2A5yLi+UbGZWZmpTWjoXky298oY4BXbqJxBnB1uZ0lLUjvu7ps3bp1DQrRzKw7NaOhudQNTZJ7GEZUnN8+IhZKehKYN27cuD+td3BmZt2sGSWFAba/29M+ZL/DFwARcXNELJgwYULljc3MLLNmJIV7gemSDpA0DjgJuKmaA0iaJ2nhc88915AAzcy6VaO7pF5HcovEmZIGJJ2R3jf2HOA2kml7b4iIB8sdp5hLCmZmjdHQNoWImD/M8ltIbgNYE0nzgHnTpk2r9RBmZlZCW05zMZKSwtrFa1k6ZSlLxixh6ZSlrF28tgERmpm1p7ZMCrW2KaxdvJYVC1aw+bHNELD5sc2sWLDCicHMLNWWSaHWksLq81az7aXtb6+77aVtrD5vdT3DMzNrW22ZFGq1+fHNVS03M+s2bZkUaq0+Gr/f+KqWm5l1m7ZMCrVWH029cCpjdt7+LY/ZeQxTL5xaz/Bq5kZwM2u2lrmfwmjI9eWApG1h8+ObGb/feKZeOHVoeTPlG8HzbR75RnCgJeIzs+6giGh2DFUrGKdw5sMPP9zscOpi6ZSlSa+oIuP3H8/hjx7ehIjMrFNJWh4RvaXWdVX1UStrZCO4q6XMLKu2TAqdqFGN4B6bYWbVcFJoEY1qBPfYDDOrRlsmhU6cJTXXl2PmwpmM3388KGlLmLlw5ogbmetRLeXqJ7Pu0ZYNzXm9vb2xbNmyZofR0kbagF3cKwqSEkw9EpaZNUfHNTRbdiOtlqql+sklC7P21VXjFLrRSMdmVFv95PEWZu3NSaEL5PpyNX8hj99vfOnqp2F6RZUrWRTHsHbx2pYcSGjWzVx9ZGVVW/2UtWSRtausq6LMRldbJoVO7H3UqqrtFZV1vEWWtgqPsTAbfe59ZHWVtbfSkjFLoNSfnmDOtjlA5Z5Trn4yq417H9moyVqyyFKiKFcV5VKEWWO4odnqLkvD9tQLp5YsURS2VZRr5K7UoO1ShFltXFKwpshSoijXyO1ShFljuE3BWtpwv/jLtTcAJdf1TOxh7K5jXXqwrleuTaFlqo8kjQE+C+wOLIuIa5ockrWA4aqiylU/9Z/SX/JYg+sHGVw/CHhQndlwGlp9JOkqSU9JeqBo+VxJKyStknRuuvhYYDLwMjDQyLis/ZWrfso63Xi+DcJjIcxe0dDqI0lvBTYCX4uI16TLeoCVwJ+TfPnfC8wH3gU8GxH/KenGiDih0vFdfWSllOoWW86Yncd4wj/rKk3rkhoRdwLPFC0+FFgVEasjYgtwPUkpYQB4Nt1mcLhjSlogaZmkZevWrWtE2NbmSpUixk4cpqa0h5K9mPpP7XfJwbpSM9oUJgNPFLweAA4DLgL+Q9IRwJ3D7RwRC4GFkJQUGhintbHitojhBtUNW5pIf5a47cG6TTO6pKrEsoiIlyLijIj4cERcUvYAnubCqjRcG0S+t1I5217axsqPrnS7g3WFZpQUBoB9C17vA6yp5gARcTNwc29v75n1DMw623A9mbK0P7jnknWLZpQU7gWmSzpA0jjgJOCmag7gkoLVS3EJgp5s+7ndwTpVo7ukXgcsBWZKGpB0RkRsBc4BbgP6gRsi4sFGxmFWTq4vx+GPHs6cbXOYdc2sV42iHtYgHjFtHccjms2KFI+iHtw4yNb1Wyvul/W+12bNVq5LqpOCWQVVjXsQnkLDWl7HTZ3tNgUbTVW1O7g6ydqcSwpmVaqm5DB+f5carPW4pGBWR68qOZThUoO1G5cUzEZouGm8t9MDbHN7g7WGjispmLWSUjcDehV3X7U20ZZJwdVH1kq2q07KYNtL2+g/ud+D3qwltWVSiIibI2LBhAkTmh2KGfDKALhZ12Yf/OZSg7WitkwKZq2q2mkz8jf6MWsVTgpmdVbttBmbH9vsqiRrGW2ZFNymYO0ia3uDq5KsVbhLqtkoyTrozQPerNHcJdWsBbjUYO3AScFsFOXbGyolhvz9GpwYbLQ5KZg1QdYBby4x2GhzUjBrgqxVSe6yaqOtLZOCex9ZJ8g64M1dVm00tWVS8Ihm6yT5UkO5gW5ufLbRMrbcSklvzHCMlyPi/jrFY9aV8t1Py3VZzTc+F25vVm9lkwLwE+Beys8afwAwpV4BmXWr/Bf96vNWDz8Vd9r4XLi9WT2VHbwm6b8j4qiyB8iwTaN48Jp1qiz3aPAgN6tVzYPXsnzZNyshmHWyLF1W3c5gjZCpoVnS3EYHImmOpLskXS5pTqPPZ9bKsjQ+g7usWv2VTQqSDpXUA3yuYNnXsx5c0lWSnpL0QNHyuZJWSFol6dx0cQAbgR2BgczvwKxD5fpymWdZdWnB6qVSSWE+SWPzVEmfl/ReIEuPpLxFwHaljDTJXAIcA8wG5kuaDdwVEccAnwA+U8U5zDpW1kFurkayeqmUFD4REX8GPA78ANgT+CNJP5f0zUoHj4g7gWeKFh8KrIqI1RGxBbgeODYi8v3wngWy3dfQrAtkGeTmuZKsXip1Sb1N0lZgb2AicAvw/oh4k6R9ajznZOCJgtcDwGGSjgf+AtgDuHi4nSUtABYA7LfffjWGYNZ+8r2M+k/uL72Bu6taHVTqffQ24F0kdf0HAhcA0yR9FzipxnOWGvMQEfGdiPhgRJwYEUvKxLSQpHrpl+PGjasxBLP2lOvLla1KconBRqpi76OI2AQ8ERFfjIj3Av8LnAk8UuM5B4B9C17vA6yp5gCe5sK6WcXuqp5d1UYgU5fUorEIF0fEuoj4do3nvBeYLukASeNIShw3VXMAT4hn3SxLd1WXGKxWVU+IFxFXZt1W0nXAUmCmpAFJZ0TEVuAc4DagH7ghIh6sMgaXFKyrZequ6hKD1aDShHgLI2JBrdtExPxhlt9C0mhdE0nzgHnTpk2r9RBmbW+o4fnUfhgsvY0n0bNqVep9dJyk35dZL+BtdYwnk4i4Gbi5t7f3zNE+t1kryTK7qnslWTUqJYV/yHCMu+oRSDVcUjB7hUsMVk9lZ0ltdZ4l1ewVaxevLV9iAMbsPIaZC2c6MXS5crOkViopmFmbcInB6qEtb8fpLqlmpblXko1UxaQgqUfSv45GMFm5S6rZ8DyOwUYiy4jmQeBPJZW7JaeZtRCXGKxWWauP/gf4nqRTJB2ffzQysHJcfWRWWdYSw8qPrhy9oKzlZU0KewLrgaOAeenjnY0KqhJXH5llk6XEMLh+0KUFG5Kp91FEnN7oQMysMbL0SnKPJMvLeo/mGZLuyN9WU9LrJH2qsaGZWb3kSwzDGkzu03D3Xne71NDlslYffRX4JPAyQET8mtrvpzBiblMwq16uL8fYieUrB7au3+rG5y6XNSnsHBG/KFq2td7BZOU2BbPaTL9oevkeSbjxudtlTQpPSzoQCABJJwBPNiwqM2uILD2SIGl8dlVSd8o6zcXfAAuBgyT9luSua30Ni8rMGibTzKq8UpVUuI91vqx3XlsdEUcDewMHRcSfRcRjjQ3NzBolX2LomVi+yOCqpO6TtffRRElfIZkme4mkiyRNbGxoZtZIub4cRzx9RMXGZ1cldZesbQrXA+uAdwMnpM+/2aigKnHvI7P6ydL47F5J3SPziOaI+GxEPJI+LgD2aGBcZbn3kVn9uCrJCmVNCj+WdJKkMenjPcAPGhmYmY2eaqqSVn7IiaGTlU0Kkl6Q9DzwQeAbwJb0cT3wt40Pz8xGU5aqpDWXrXEbQwcre/UjYreI2D39d0xEjE0fYyJi99EK0sxGR9aqpK3rt9J/Sr9LDR0o8+04Jb0OmFK4T0R8pwExmVkT5fpy5Ppy3L3X3WxdX2bigkhKDQAzLp0xStFZo2XtknoVcBVJ76OGTZ0taRdJyyU1bVpuM0tMv2g6ZLi11prL1rjE0EGyNjS/KSJ6I+LUiDg9fby/0k6SrpL0VH521YLlcyWtkLRK0rkFqz4B3FBF/GbWILm+HJPOmpQ5MbidoTNkTQpLJc2u4fiLgLmFCyT1AJcAxwCzgfmSZks6GvgN4L8qsxYx49IZzPr6rIptDOB2hk6RtU3hGpLE8DtgM8lvh4iI15XbKSLulDSlaPGhwKqIWA0g6XrgWGBXYBeSRLFJ0i0R8aqJWSQtABYA7LfffhnDN7Na5dsYVn5o5VAbwrDcztD2siaFq4BTgPuB4WfQymYy8ETB6wHgsIg4B0DSacDTpRICQEQsJJmcj97e3hhhLGaWUf5LvmJiSLdZ+/W1zLh8hifTazNZk8LjEXFTnc5ZqoZy6Ms9IhZVPIA0D5g3bdq0OoVkZlkMJYbL1xT8ry1tcOMg/Sf389xPn3OpoY1kbVN4SNI3JM2XdHz+UeM5B4B9C17vA1T+6VHA01yYNU817Qzg3kntJmtS2ImkLeEdjLxL6r3AdEkHSBpHclvPqkohnhDPrLny02JMOntSpu2dGNpHpuqjiDi9loNLug6YA+wlaQD4p4i4UtI5wG0k93+6KiIerOX4ZtZc1bYzFO5jrUkRldtqJV1NiRrELGMVGqm3tzeWLVvWzBDMDLL1TEr17NrjBugmk7Q8InpLrctaffR9kllRfwDcAewObKxPeGbW7mZcOoNZ185Cu1Qe6ZZvgL5rt7s82K0FZSopvGonaQxwe0QcVf+QMp0/3/vozIcffrgZIZjZMKopNYBLDs1Qj5JCselA00aOufeRWeuacemMzA3Q8ErJwQ3RrSFTQ7OkF9i+TeF3JPMUNYXHKZi1tmoaoPPcEN0aMpUUCu6rkH/MiIhvNzq4MvG4pGDW4qotMYC7rraCsiUFSWWriCLi8fqGY2adZMalM5jwlgk89MGHiBeztV+6xNBcZRuaJd1PUm1U2KUggL2BP4yIbEMa68wNzWbtZ+3itVUlBzdAN07NDc0R8dqIeF3672tJRjL/lKQ76sfqHmlGrj4yaz+5vhxHbjyy6q6rrk4aXVnvvDZd0iLgh8ByYHZE/EcjAzOzzpRPDtVMkeExDaOnbFKQ9Jp0qopvA7cDr4mIKyLi5VGJzsw6VjUN0S41jJ5KXVJ/RXLvgx+Q3BznUOmVYl9EfKRxoQ3PXVLNOkO1XVfdCN14lRqaTy23c0RcU/eIquC5j8w6g0dBj65yDc1lSwrN/tI3s+5QbYnBN/BpnFqnuTAzq6tqJtXL82C3+nNSMLOWUdhtteSNe0twYqivrF1S35JlmZlZPeT6csz6+qzkNlwZrLlsDUvGLHFyqIOsJYVSYxKaNk7Bt+M063y5vhyzrqmiOimS5HDf0fc1NK5OV6n30eHAm0lGL3+5YNXuwF9HxOsbGl0F7n1k1h2q7Z20x9v34ODbD25cQG1uJPdTGAfsStJLabeCx/PACfUM0sxsOPlG6KztDBvu2OASQ40qdUn9CfATSYsi4rFRisnM7FXyYxL6T+2Hwcrbb7hjAys/tNJdVquU6SY7wHhJC4Ephfs063acZtad8okh62yray5bw0srX3JVUhUy3aNZ0q+Ay0kmwxvK0RGxvHGhVeY2BbPutXbxWvrf3w9bKm/rNobt1eMezVsj4rKI+EVELM8/6hijmVlVcn055myewx5v36Pithvu2OCZVjPKmhRulvQhSX8sac/8o56BSJol6XJJN0o6u57HNrPOdfDtB6MdM96f4dR+J4YKsiaFU4F/AH5GUoW0HKhYbyPpKklPSXqgaPlcSSskrZJ0LkBE9EfEWcB7gJLFGjOzUg664qBsA90GSaqcbFiZkkJEHFDiMTXDrouAuYULJPUAlwDHALOB+ZJmp+veBdwN3FHFezCzLpcf6Ma4DBtvgXv+5J6Gx9Susk5zsbOkT6U9kPJ3Yntnpf0i4k7gmaLFhwKrImJ1RGwBrgeOTbe/KSLeDPSViWWBpGWSlq1bty5L+GbWBfJtDFmqkjb9ZpPbGIaRtfroapI2/jenrweAC2o852SSG/fkDQCTJc2R9BVJ/wncMtzOEbEQ+Azwy3HjsvwsMLNukrUqaXDjIP2nu42hWNZxCgdGxImS5gNExCYV3oKtOqX2i4hYAizJcoCIuBm4ube398waYzCzDlXVWIaXof/k/u3263ZZSwpbJO0EBICkA4HNNZ5zANi34PU+QPZJTfCEeGZWXn4K7izdVQHf/7lA1qRwPnArsK+kxSQNwR+v8Zz3AtMlHSBpHHAScFM1B4iImyNiwYQJE2oMwcy6wcG3H8xOs3fKtK3vy5DI2vvoR8DxwGnAdUBvWt1TlqTrgKXATEkDks6IiK3AOcBtQD9wQ0Q8WE3QLimYWVaHPXgYPbtmuzHDmsvWdH0bQ9ZpLm4iSQY3RcSLDY8qI09zYWZZrF28lv7T++HlyttqF3HkxiMbH1QT1WOaiy8BRwC/kfQtSSdI2rFuEVbJJQUzq0auL8esq7ONY4gXo6urkTKVFIY2TgaeHQWcCcyNiN0bFVgWLimYWbXu2u0uBjdWnnt70tmTOnba7XqUFEh7H70bOAs4BLimPuGZmY2eGZfPgB0qb9etDc9ZRzR/k6RR+CiSKSoOjIgPNzKwCvG4+sjMapKvSspy7+duvOdz1obmucB/RUSG+x2NHlcfmdlIdGtVUs3VR5I+DhARt5J0SS1c97m6RWhm1gQzLs/2Rd9NXVUrVR+dVPD8k0Xr5tIkrj4ys3rI9eUYOzHbbD8Pvf+hrkgMlZKChnle6vWo8YhmM6uX6RdNz9TwHFuChz74UOMDarJKSSGGeV7qtZlZ26mm4TlejI4vLVRKCq+X9LykF4DXpc/zr187CvGZmTVcfgK9SWdPqrht/8n9LJ2ytGOTQ9mkEBE9EbF7ROwWEWPT5/nXGQpcjeE2BTNrhKw9jDY/tpkVC1Z0ZGLIPHitlbhNwcwapWditsnztr20jdXnrW5wNKOvLZOCmVmjzLgo24hnSEoMncZJwcysQL7hOWuJodPu9eykYGZWJNeX44inj8jU8X5w42BHjWFwUjAzG8aksyr3RoJkDEOntC+0ZVJw7yMzGw0zLp2Bxmcbp9sp7QttmRTc+8jMRstBVx6Uedu797q77auR2jIpmJmNllxfjlnXzsrUvrB1/da2b1/INhOUmVkXy/XlgGQ0cyWxJeg/pX+7/dqJSwpmZhnk+nL07JqtmyqRJJB2rE5yUjAzy2jG5TPQuOwTRG9dv7XtpsNoqaQg6ThJX5X0PUnvaHY8ZmaFcn05Droqe8MztN90GA1PCpKukvSUpAeKls+VtELSKknnAkTEdyPiTOA04MRGx2ZmVq1cX47x+4+vap/Nj7dPd9XRKCksougubZJ6gEuAY4DZwHxJsws2+VS63sys5Uy9cCpjdq7i6zNgSc8SVn5oZaNCqpuGJ4WIuBN4pmjxocCqiFgdEVuA64Fjlfg88MOI+GWjYzMzq0WuL8fMhTMhY7szANuSez3fd/R9jQqrLprVpjAZeKLg9UC67MPA0cAJks4qtaOkBZKWSVq2bt26xkdqZlZCri/HrGtmVVdiADbcsaGlG56bNU6hVPN9RMRXgK+U2zEiFkp6Epg3bty4P21IdGZmGVQzfqFQK49jaFZJYQDYt+D1PsCarDt7mgszaxW1NDznxzEsGdN67QzNSgr3AtMlHSBpHHAScFPWnT0hnpm1kqkXTq1tx0jaGVopMYxGl9TrgKXATEkDks6IiK3AOcBtQD9wQ0Q8mPWYLimYWSvJ9eUYO7H22vg1l61pmdHPDW9TiIj5wyy/BbillmNKmgfMmzZt2khCMzOrm+kXTWfFghVse2lbTfvnJ9OD5rY1tNSI5qxcUjCzVlNTN9UisSXoP7W/qSWGtkwKblMws1ZUazfV7QzS1PmS2jIpuKRgZq0qX2KoukdSgWbOl9SWScElBTNrZbm+HIc/ejizrq291NCs23u2ZVJwScHM2sGISg0jaJsYibZMCmZm7SJfaqg6MQw2557PbZkUXH1kZu2m6plVac49n9syKbj6yMzazVCX1SrFlhjVRue2TApmZu2opnmSGN2b9DgpmJmNolqqkRjDqFUhOSmYmY2imnokjeKAtrZMCm5oNrN2lu+RNCfmMOnsSZn22fbSNvpPafwUGG2ZFNzQbGadYsalM9COpe47VkJA//samxjaMimYmXWKtYvXMmaXKr6Kt8HDH324YfE063acZmZdb+3itTVNt711/dYGReSSgplZ06w+b3XN919olLZMCm5oNrNOMJLxB0unLG1I20JbJgU3NJtZJxi/X+3Ta29+bHNDuqm2ZVIwM+sEJQey7ZB9/0bcd8FJwcysSbYbyCYYv/94Zl09CzL2UIX6T4Hh3kdmZk2U68uR68ttt+y5nz7HmsvWZNp/JFVQpbikYGbWYmZcOoM93r5Hpm2nXji1rud2UjAza0GbVm2quI120atKGSPVMklB0lRJV0q6sdmxmJk1W6a2gsH6z57a0KQg6SpJT0l6oGj5XEkrJK2SdC5ARKyOiDMaGY+ZWbvI0lYQvw8eOq2+d2ZrdElhETC3cIGkHuAS4BhgNjBf0uwGx2Fm1laythXE1vrema2hSSEi7gSeKVp8KLAqLRlsAa4Hjs16TEkLJC2TtGzdunV1jNbMrHXk+nKMnZitg2g9u6U2o01hMvBEwesBYLKkiZIuB94g6ZPD7RwRCyOiNyJ6995770bHambWNNMvmp7pLm317JbajHEKpYZlRESsB87KdABpHjBv2rRpdQ3MzKyV5HsWrfzoSgbXDw67XT27pTYjKQwA+xa83gfINkrDzKzL5BND/8n9Fbeph2ZUH90LTJd0gKRxwEnATdUcwBPimVk3KdeQ3DOxp67nanSX1OuApcBMSQOSzoiIrcA5wG1AP3BDRDxY5XE9dbaZdY1yDckzLppR13MpIup6wNHU29sby5Yta3YYZmYNtXTKUjY/9urE0DOxhyOePqLq40laHhG9pda1zIjmarikYGbdpNQU22N2HlP3UgK0aVJwm4KZdZNSU2zPXDiz7vMegafONjNrC6Wm2G6EtiwpuPrIzKwx2jIpuPrIzKwx2jIpmJlZY7RlUnD1kZlZY7RlUnD1kZlZY7T14DVJ64DH0pcTgFJFh+GW7wU83aDQRmq4mJt93Gr3z7p9lu3KbVPtOl/7xu/va1+90bz2+0dE6WmmI6IjHsDCKpcva3bM1b6XZh+32v2zbp9lu3LbVLvO197XvhUfrXLt27L6aBg3V7m8lTUq5pEet9r9s26fZbty29S6rhX52le3ja99nY/b1tVHIyFpWQwz94d1Nl/77uVrX1knlRSqtbDZAVjT+Np3L1/7Crq2pGBmZq/WzSUFMzMr4qRgZmZDnBTMzGyIk0JK0i6SrpH0VUl9zY7HRo+kqZKulHRjs2Ox0SXpuPT//PckvaPZ8bSCjk4Kkq6S9JSkB4qWz5W0QtIqSeemi48HboyIM4F3jXqwVlfVXPuIWB0RZzQnUqu3Kq/9d9P/86cBJzYh3JbT0UkBWATMLVwgqQe4BDgGmA3MlzQb2Ad4It1scBRjtMZYRPZrb51lEdVf+0+l67teRyeFiLgTeKZo8aHAqvTX4RbgeuBYYIAkMUCHfy7doMprbx2kmmuvxOeBH0bEL0c71lbUjV9+k3mlRABJMpgMfAd4t6TLaL/h8ZZNyWsvaaKky4E3SPpkc0KzBhvu//2HgaOBEySd1YzAWk033qNZJZZFRLwInD7awdioGu7arwf8hdDZhrv2XwG+MtrBtLJuLCkMAPsWvN4HWNOkWGx0+dp3L1/7jLoxKdwLTJd0gKRxwEnATU2OyUaHr3338rXPqKOTgqTrgKXATEkDks6IiK3AOcBtQD9wQ0Q82Mw4rf587buXr/3IeEI8MzMb0tElBTMzq46TgpmZDXFSMDOzIU4KZmY2xEnBzMyGOCmYmdkQJwVrCkmDku6T9ICkb0nauQ7H/GdJR1fY5ixJ70ufnyZpUobjbredpCvqMbuqpPMl/VbSP4/0WGXOcaukyUXL5kj6fo3HOzGderqm/a31OSlYs2yKiIMj4jXAFormHkqnOq5KRHw6Im6vsM3lEfG19OVpQMWkULxdRHwgIn5TbXzD+HJEfLpOx9qOpJ2APSPit/U6ZkR8E/hAvY5nrcdJwVrBXcC09BfsjyV9A7hfUo+kf5V0r6RfS/pgfgdJH5d0v6RfSfqXdNkiSSekzx+V9HlJv0gf09Ll50v6+3S7XmBxWmLZSdKn03M9IGlhOq1yqe2WSOpNjzc/jeOBdArmfHwbJV2YxvdzSblKH0Ia2zWSfpTGf7ykL6THv1XSDgXv7XOSlkpaJumNkm6T9L9FM33OAZak+8yV9JCku0luKJU/56GSfibpf9J/Z6bL75J0cMF2P5X0uuyX1NqVk4I1laSxJDc+uT9ddChwXkTMBs4AnouIQ4BDgDPTuWuOAY4DDouI1wNfGObwz0fEocDFwL8XroiIG4FlQF9aYtkEXBwRh6Sll52Adw6zXT72ScDngaOAg4FDJB2Xrt4F+Hka353AmRk/kgOBvyK5z8O1wI8j4rXApnR53hMRcThJQl0EnAC8CSisijoGuFXSjsBXgXnAEcAfFWzzEPDWiHgD8Gngc+nyK0hKSEiaAYyPiF9nfA/WxpwUrFl2knQfyRfu48CV6fJfRMQj6fN3AO9Lt7sHmAhMJ5n//uqIeAkgIopvqJJ3XcG/h2eI6W2S7pF0P8kX/Z9U2P4QYElErEvn1lkMvDVdtwXI17svB6ZkOD8kN3t5mSRJ9gC3psvvLzrGTQXL74mIFyJiHfB7SXuk694C3A0cBDwSEQ9HMq/NtQXHmQB8S8mtK7/MK+/5W8A709LJ+0kSj3WBbryfgrWGTRFxcOECSQAvFi4CPhwRtxVtNxfIMmlXDPP8VdJf05cCvRHxhKTzgR0rHL/UHP15L8crE4sNkv3/2maAiNgmqfAY24qOsblg+eaC5duAsZKmkpQmtqSf63Dv/7MkpZG/ljSFtLopIl6S9F8kJZb3kFShWRdwScFa2W3A2QV16TMk7QL8CHi/0h5LkvYcZv8TC/5dWmL9C8Bu6fN8Anha0q4k1TGltit0D3CkpL3ShvH5wE8yvbPGO4ZXShkPAQdIOjB9Pb9guwlAviH6tKJjXEFyA5p7y5TGrMO4pGCt7AqSKpNfKvm5uw44LiJuTRtBl0naAtwC/GOJ/cdLuofkx8/8EusXAZdL2kRSvfRVkuqYR0nm3x9uOwAi4kklt+/8MUmp4ZaI+F6tb7bO5pLcapKI+L2kBcAPJD1NUqX0mnS7LwDXSPo74L8LDxARyyU9D1w9emFbs3nqbOtIkh4lqQp6utmxDCetotoYEV+s83HHAz+NiBFV+aQN6UuAgyJiW8HyOcDfR8Q7R3J8a02uPjJrno3AAtV58FpEbK5DQngfSfXYeUUJ4USStpdnRxaltSqXFMzMbIhLCmZmNsRJwczMhjgpmJnZECcFMzMb4qRgZmZDnBTMzGzI/wcFa/PabZJi0AAAAABJRU5ErkJggg==\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_noresm')\n",
    "plt.savefig(\"/home/mkibria/Academic/pr_of_thesis/figures/pdf_regcm_noresm.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
}