pdf_cru.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 19,
   "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": 20,
   "id": "9b4e84e2",
   "metadata": {},
   "outputs": [],
   "source": [
    "rain_minimum = 1.0\n",
    "rain_maximum = 120.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",
    "\n",
    "\n",
    " "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "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:  (lon: 720, lat: 360, time: 1452)\n",
       "Coordinates:\n",
       "  * lon      (lon) float32 -179.8 -179.2 -178.8 -178.2 ... 178.8 179.2 179.8\n",
       "  * lat      (lat) float32 -89.75 -89.25 -88.75 -88.25 ... 88.75 89.25 89.75\n",
       "  * time     (time) datetime64[ns] 1901-01-16 1901-02-15 ... 2021-12-16\n",
       "Data variables:\n",
       "    pre      (time, lat, lon) float32 ...\n",
       "    stn      (time, lat, lon) float64 ...\n",
       "Attributes:\n",
       "    Conventions:  CF-1.4\n",
       "    title:        CRU TS4.06 Precipitation\n",
       "    institution:  Data held at British Atmospheric Data Centre, RAL, UK.\n",
       "    source:       Run ID = 2205201912. Data generated from:pre.2204291209.dtb\n",
       "    history:      Fri 20 May 21:52:05 BST 2022 : User f098 : Program makegrid...\n",
       "    references:   Information on the data is available at http://badc.nerc.ac...\n",
       "    comment:      Access to these data is available to any registered CEDA user.\n",
       "    contact:      support@ceda.ac.uk</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-240463b4-a99c-4cc7-b503-94ff1e993909' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-240463b4-a99c-4cc7-b503-94ff1e993909' 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'>lon</span>: 720</li><li><span class='xr-has-index'>lat</span>: 360</li><li><span class='xr-has-index'>time</span>: 1452</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-f68fae46-36e3-4ade-9bc5-9f7f872af3a8' class='xr-section-summary-in' type='checkbox'  checked><label for='section-f68fae46-36e3-4ade-9bc5-9f7f872af3a8' class='xr-section-summary' >Coordinates: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>-179.8 -179.2 ... 179.2 179.8</div><input id='attrs-3dbb9449-0027-49c4-90b0-00d9be9b888d' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-3dbb9449-0027-49c4-90b0-00d9be9b888d' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-299ff837-77d3-4aba-9d46-bb8b3951ea1a' class='xr-var-data-in' type='checkbox'><label for='data-299ff837-77d3-4aba-9d46-bb8b3951ea1a' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>array([-179.75, -179.25, -178.75, ...,  178.75,  179.25,  179.75],\n",
       "      dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>-89.75 -89.25 ... 89.25 89.75</div><input id='attrs-4ac02725-1984-49af-af3f-a5d3eab863e9' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-4ac02725-1984-49af-af3f-a5d3eab863e9' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-a4ba41f4-7adf-46b1-b4a7-410a1a2f1310' class='xr-var-data-in' type='checkbox'><label for='data-a4ba41f4-7adf-46b1-b4a7-410a1a2f1310' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd></dl></div><div class='xr-var-data'><pre>array([-89.75, -89.25, -88.75, ...,  88.75,  89.25,  89.75], dtype=float32)</pre></div></li><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'>1901-01-16 ... 2021-12-16</div><input id='attrs-99acc102-7a72-4c53-88a5-d7da60f3696b' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-99acc102-7a72-4c53-88a5-d7da60f3696b' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-94c66d7b-f176-4e1d-99e1-1c427f70809f' class='xr-var-data-in' type='checkbox'><label for='data-94c66d7b-f176-4e1d-99e1-1c427f70809f' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;1901-01-16T00:00:00.000000000&#x27;, &#x27;1901-02-15T00:00:00.000000000&#x27;,\n",
       "       &#x27;1901-03-16T00:00:00.000000000&#x27;, ..., &#x27;2021-10-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;2021-11-16T00:00:00.000000000&#x27;, &#x27;2021-12-16T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-d02b9834-dad9-49bc-bd29-51f9498e0453' class='xr-section-summary-in' type='checkbox'  checked><label for='section-d02b9834-dad9-49bc-bd29-51f9498e0453' class='xr-section-summary' >Data variables: <span>(2)</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>pre</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-65f20a19-429d-4eb1-b811-363c069d4cd5' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-65f20a19-429d-4eb1-b811-363c069d4cd5' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-d154657d-5961-4734-aee4-fdcf8c63d141' class='xr-var-data-in' type='checkbox'><label for='data-d154657d-5961-4734-aee4-fdcf8c63d141' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>precipitation</dd><dt><span>units :</span></dt><dd>mm/month</dd><dt><span>correlation_decay_distance :</span></dt><dd>450.0</dd></dl></div><div class='xr-var-data'><pre>[376358400 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>stn</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-66c61a3b-386f-4331-b8b3-8411fb71eff2' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-66c61a3b-386f-4331-b8b3-8411fb71eff2' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-6bb31239-6f6c-40ea-8397-c39711399670' class='xr-var-data-in' type='checkbox'><label for='data-6bb31239-6f6c-40ea-8397-c39711399670' 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>description :</span></dt><dd>number of stations contributing to each datum</dd></dl></div><div class='xr-var-data'><pre>[376358400 values with dtype=float64]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-09239e3b-a0d5-4b8a-8d40-b56b271a4a52' class='xr-section-summary-in' type='checkbox'  checked><label for='section-09239e3b-a0d5-4b8a-8d40-b56b271a4a52' class='xr-section-summary' >Attributes: <span>(8)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>Conventions :</span></dt><dd>CF-1.4</dd><dt><span>title :</span></dt><dd>CRU TS4.06 Precipitation</dd><dt><span>institution :</span></dt><dd>Data held at British Atmospheric Data Centre, RAL, UK.</dd><dt><span>source :</span></dt><dd>Run ID = 2205201912. Data generated from:pre.2204291209.dtb</dd><dt><span>history :</span></dt><dd>Fri 20 May 21:52:05 BST 2022 : User f098 : Program makegridsauto.for called by update.for</dd><dt><span>references :</span></dt><dd>Information on the data is available at http://badc.nerc.ac.uk/data/cru/</dd><dt><span>comment :</span></dt><dd>Access to these data is available to any registered CEDA user.</dd><dt><span>contact :</span></dt><dd>support@ceda.ac.uk</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:  (lon: 720, lat: 360, time: 1452)\n",
       "Coordinates:\n",
       "  * lon      (lon) float32 -179.8 -179.2 -178.8 -178.2 ... 178.8 179.2 179.8\n",
       "  * lat      (lat) float32 -89.75 -89.25 -88.75 -88.25 ... 88.75 89.25 89.75\n",
       "  * time     (time) datetime64[ns] 1901-01-16 1901-02-15 ... 2021-12-16\n",
       "Data variables:\n",
       "    pre      (time, lat, lon) float32 ...\n",
       "    stn      (time, lat, lon) float64 ...\n",
       "Attributes:\n",
       "    Conventions:  CF-1.4\n",
       "    title:        CRU TS4.06 Precipitation\n",
       "    institution:  Data held at British Atmospheric Data Centre, RAL, UK.\n",
       "    source:       Run ID = 2205201912. Data generated from:pre.2204291209.dtb\n",
       "    history:      Fri 20 May 21:52:05 BST 2022 : User f098 : Program makegrid...\n",
       "    references:   Information on the data is available at http://badc.nerc.ac...\n",
       "    comment:      Access to these data is available to any registered CEDA user.\n",
       "    contact:      support@ceda.ac.uk"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\n",
    "pr = xr.open_dataset('/home/mkibria/Academic/pr_of_thesis/observation/cru/cru_ts4.06.1901.2021.pre.dat.nc')\n",
    "\n",
    "pr\n",
    "\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "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": 23,
   "id": "f02ba2d8",
   "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:  (lon: 16, lat: 18, time: 431)\n",
       "Coordinates:\n",
       "  * lon      (lon) float32 86.25 86.75 87.25 87.75 ... 92.25 92.75 93.25 93.75\n",
       "  * lat      (lat) float32 18.25 18.75 19.25 19.75 ... 25.25 25.75 26.25 26.75\n",
       "  * time     (time) datetime64[ns] 1980-02-15 1980-03-16 ... 2015-12-16\n",
       "Data variables:\n",
       "    pre      (time, lat, lon) float32 ...\n",
       "    stn      (time, lat, lon) float64 ...\n",
       "Attributes:\n",
       "    Conventions:  CF-1.4\n",
       "    title:        CRU TS4.06 Precipitation\n",
       "    institution:  Data held at British Atmospheric Data Centre, RAL, UK.\n",
       "    source:       Run ID = 2205201912. Data generated from:pre.2204291209.dtb\n",
       "    history:      Fri 20 May 21:52:05 BST 2022 : User f098 : Program makegrid...\n",
       "    references:   Information on the data is available at http://badc.nerc.ac...\n",
       "    comment:      Access to these data is available to any registered CEDA user.\n",
       "    contact:      support@ceda.ac.uk</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-af61938d-dde4-4c54-a1d7-f1dd478d33cc' class='xr-section-summary-in' type='checkbox' disabled ><label for='section-af61938d-dde4-4c54-a1d7-f1dd478d33cc' 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'>lon</span>: 16</li><li><span class='xr-has-index'>lat</span>: 18</li><li><span class='xr-has-index'>time</span>: 431</li></ul></div><div class='xr-section-details'></div></li><li class='xr-section-item'><input id='section-6d0a3631-30cb-4e95-b495-0f0860e488f6' class='xr-section-summary-in' type='checkbox'  checked><label for='section-6d0a3631-30cb-4e95-b495-0f0860e488f6' class='xr-section-summary' >Coordinates: <span>(3)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><ul class='xr-var-list'><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lon</span></div><div class='xr-var-dims'>(lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>86.25 86.75 87.25 ... 93.25 93.75</div><input id='attrs-0a7d98f0-9d9f-414b-aa8c-7815b8ab6f12' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-0a7d98f0-9d9f-414b-aa8c-7815b8ab6f12' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-de9bcc1f-0005-4875-a573-d9f11d766f42' class='xr-var-data-in' type='checkbox'><label for='data-de9bcc1f-0005-4875-a573-d9f11d766f42' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>longitude</dd><dt><span>units :</span></dt><dd>degrees_east</dd></dl></div><div class='xr-var-data'><pre>array([86.25, 86.75, 87.25, 87.75, 88.25, 88.75, 89.25, 89.75, 90.25, 90.75,\n",
       "       91.25, 91.75, 92.25, 92.75, 93.25, 93.75], dtype=float32)</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span class='xr-has-index'>lat</span></div><div class='xr-var-dims'>(lat)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>18.25 18.75 19.25 ... 26.25 26.75</div><input id='attrs-7eed74c6-b1a5-44f1-aa02-b38f9e444da8' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-7eed74c6-b1a5-44f1-aa02-b38f9e444da8' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-4581d21b-4c53-49bd-b8de-1d85972b5fef' class='xr-var-data-in' type='checkbox'><label for='data-4581d21b-4c53-49bd-b8de-1d85972b5fef' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>latitude</dd><dt><span>units :</span></dt><dd>degrees_north</dd></dl></div><div class='xr-var-data'><pre>array([18.25, 18.75, 19.25, 19.75, 20.25, 20.75, 21.25, 21.75, 22.25, 22.75,\n",
       "       23.25, 23.75, 24.25, 24.75, 25.25, 25.75, 26.25, 26.75], dtype=float32)</pre></div></li><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'>1980-02-15 ... 2015-12-16</div><input id='attrs-94caacdd-1777-408d-b9a1-953cceafceff' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-94caacdd-1777-408d-b9a1-953cceafceff' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-594c4825-623d-4558-88c5-0a56c6bac551' class='xr-var-data-in' type='checkbox'><label for='data-594c4825-623d-4558-88c5-0a56c6bac551' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>time</dd></dl></div><div class='xr-var-data'><pre>array([&#x27;1980-02-15T00:00:00.000000000&#x27;, &#x27;1980-03-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;1980-04-16T00:00:00.000000000&#x27;, ..., &#x27;2015-10-16T00:00:00.000000000&#x27;,\n",
       "       &#x27;2015-11-16T00:00:00.000000000&#x27;, &#x27;2015-12-16T00:00:00.000000000&#x27;],\n",
       "      dtype=&#x27;datetime64[ns]&#x27;)</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-b98abc9c-c9b6-45b3-992d-a2a5eb2709ef' class='xr-section-summary-in' type='checkbox'  checked><label for='section-b98abc9c-c9b6-45b3-992d-a2a5eb2709ef' class='xr-section-summary' >Data variables: <span>(2)</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>pre</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float32</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d852fc60-a529-4c66-9166-d426b674660e' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d852fc60-a529-4c66-9166-d426b674660e' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-c5ec672d-8300-4368-95e8-0a50aaa816b6' class='xr-var-data-in' type='checkbox'><label for='data-c5ec672d-8300-4368-95e8-0a50aaa816b6' title='Show/Hide data repr'><svg class='icon xr-icon-database'><use xlink:href='#icon-database'></use></svg></label><div class='xr-var-attrs'><dl class='xr-attrs'><dt><span>long_name :</span></dt><dd>precipitation</dd><dt><span>units :</span></dt><dd>mm/month</dd><dt><span>correlation_decay_distance :</span></dt><dd>450.0</dd></dl></div><div class='xr-var-data'><pre>[124128 values with dtype=float32]</pre></div></li><li class='xr-var-item'><div class='xr-var-name'><span>stn</span></div><div class='xr-var-dims'>(time, lat, lon)</div><div class='xr-var-dtype'>float64</div><div class='xr-var-preview xr-preview'>...</div><input id='attrs-d72d5846-6e0f-490a-92c2-80f47d1df258' class='xr-var-attrs-in' type='checkbox' ><label for='attrs-d72d5846-6e0f-490a-92c2-80f47d1df258' title='Show/Hide attributes'><svg class='icon xr-icon-file-text2'><use xlink:href='#icon-file-text2'></use></svg></label><input id='data-16e623e3-3a8e-4e4e-8fc4-0a8e02491c34' class='xr-var-data-in' type='checkbox'><label for='data-16e623e3-3a8e-4e4e-8fc4-0a8e02491c34' 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>description :</span></dt><dd>number of stations contributing to each datum</dd></dl></div><div class='xr-var-data'><pre>[124128 values with dtype=float64]</pre></div></li></ul></div></li><li class='xr-section-item'><input id='section-f4b8afe7-6c58-41aa-994c-ada15f06bd26' class='xr-section-summary-in' type='checkbox'  checked><label for='section-f4b8afe7-6c58-41aa-994c-ada15f06bd26' class='xr-section-summary' >Attributes: <span>(8)</span></label><div class='xr-section-inline-details'></div><div class='xr-section-details'><dl class='xr-attrs'><dt><span>Conventions :</span></dt><dd>CF-1.4</dd><dt><span>title :</span></dt><dd>CRU TS4.06 Precipitation</dd><dt><span>institution :</span></dt><dd>Data held at British Atmospheric Data Centre, RAL, UK.</dd><dt><span>source :</span></dt><dd>Run ID = 2205201912. Data generated from:pre.2204291209.dtb</dd><dt><span>history :</span></dt><dd>Fri 20 May 21:52:05 BST 2022 : User f098 : Program makegridsauto.for called by update.for</dd><dt><span>references :</span></dt><dd>Information on the data is available at http://badc.nerc.ac.uk/data/cru/</dd><dt><span>comment :</span></dt><dd>Access to these data is available to any registered CEDA user.</dd><dt><span>contact :</span></dt><dd>support@ceda.ac.uk</dd></dl></div></li></ul></div></div>"
      ],
      "text/plain": [
       "<xarray.Dataset>\n",
       "Dimensions:  (lon: 16, lat: 18, time: 431)\n",
       "Coordinates:\n",
       "  * lon      (lon) float32 86.25 86.75 87.25 87.75 ... 92.25 92.75 93.25 93.75\n",
       "  * lat      (lat) float32 18.25 18.75 19.25 19.75 ... 25.25 25.75 26.25 26.75\n",
       "  * time     (time) datetime64[ns] 1980-02-15 1980-03-16 ... 2015-12-16\n",
       "Data variables:\n",
       "    pre      (time, lat, lon) float32 ...\n",
       "    stn      (time, lat, lon) float64 ...\n",
       "Attributes:\n",
       "    Conventions:  CF-1.4\n",
       "    title:        CRU TS4.06 Precipitation\n",
       "    institution:  Data held at British Atmospheric Data Centre, RAL, UK.\n",
       "    source:       Run ID = 2205201912. Data generated from:pre.2204291209.dtb\n",
       "    history:      Fri 20 May 21:52:05 BST 2022 : User f098 : Program makegrid...\n",
       "    references:   Information on the data is available at http://badc.nerc.ac...\n",
       "    comment:      Access to these data is available to any registered CEDA user.\n",
       "    contact:      support@ceda.ac.uk"
      ]
     },
     "execution_count": 23,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "pr= pr.sel(lon=slice(86,94),lat=slice(18,27))\n",
    "pr"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 24,
   "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": 25,
   "id": "2100d68b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(431, 18, 16)\n"
     ]
    }
   ],
   "source": [
    "size = pr.pre.shape\n",
    "print(size)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "77d6e633",
   "metadata": {},
   "outputs": [],
   "source": [
    "for time in range(size[0]):\n",
    "        data = pr.pre[time,Ellipsis] / mpy\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]\n",
    "        #print(rain_bin)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 27,
   "id": "307f9480",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYUAAAEWCAYAAACJ0YulAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAAsTAAALEwEAmpwYAAAn40lEQVR4nO3de7RcdX338fcnCQiHSzQhbRXICRC8oKJiQKytteqyoSVClSXgAbEgKbZofVov8MRlfVrjo3ipWmg1IpAmEfD2WEAFrBXRikjwBoo0EQkErYSAEAwSSL7PH789J/tM5rJnzuwzs898XmvNOjN79uW3z0nmO7/b96eIwMzMDGBGvwtgZmaDw0HBzMzGOSiYmdk4BwUzMxvnoGBmZuMcFMzMbJyDgvWdpDFJ1/Z63w6u/7Ckg3t5zibXuUTSewru+xVJp2XPXy/pW+WWzixxULCWJN0p6ZHsg/NXki6WtHcvrxERayLiFd3sKykkLSx6LUnXSXpD3Tn3jog7ipe4fBFxTESsnOx5st/Pb7K/3z2SPixpZvZe7W+7RdKvJX1b0lmSZuSOv0TStuz42uPEyZbLBpeDghWxJCL2Bo4AjgTeWb+DpFlTXior6jnZ3+9lwGuBM3PvLYmIfYBR4H3AO4BP1R1/XhY4a4/LJ1MY/1sZbA4KVlhE3AN8BXgWjH8L/WtJ64B12bZjJf0g983z8Nrxkg6U9AVJmyRtlnR+tn1C80h23jdLukPSfZI+UPv2mt9X0vXZIT+sfYOV9CRJV2XXeCB7fkC2/3LgD4Hzs/3Pz11vYfZ8tqR/y47fIOmd9deW9MHs3D+XdEyz35ek50n6XvZN/HJgj9x7TcuZvb9LjSbbfoGkD9Vtu1LSW1r/9SAifgp8k+zvV/fegxFxBXAicJqkXfZpR9JMSf9b0s+ye75Z0oHZexP+rUhakG2blTu+4T3b1HJQsMKy/+B/Cnw/t/l44AXAYZKOAC4C/hKYC3wCuELSE7Imi6uADcACYH/gshaX+3NgEal2chxwev0OEfHi7Olzct9gZwAXk775zgceAc7P9l9G+lA8O9v/7AbX/WdgNnAw8EfA64C/yL3/AuB2YD/gPOBTklR/Ekm7A18EVgFzgM8Cr87t0rScbawETs4Fqv1INYBL2x0o6TBSUPx+s30i4rvAxmy/Tv0tcDLp38i+pL/Z1tz7x5P9W+ni3DZFHBSsiC9K+jXwLeAbwHtz7/3fiLg/Ih4hNUt8IiJujIjtWZv4o8DRwFHAU4C3RcRvIuK3EdGq8/T92XnvAj5C+rBpKyI2R8TnI2JrRGwBlpM+3NvKAteJwLkRsSUi7gQ+BJya221DRHwyIraTPqCfDPxug9MdDewGfCQiHouIzwE3Tbac2Yf2g6RAAHAScF1E/KrFYd+T9ABwJXAhKRi18gtSIKt5a1bz+7Wk+1oc9wbgnRFxeyQ/jIjNuffz/1ZsQLltz4o4PiL+o8l7d+eej5KaHt6U27Y7KRhsJ32gPl7wmvnzbsjO0ZakEeCfgMXAk7LN+0iamX2Qt7JfVt4NddfeP/f6f2pPImJrVklo1PH+FOCemJhxcvy8kyznSuAU4KvZz4+22f+IiFjfZp+8/YH7c68/GBET+pEkjZFqggDfjIhjgAOBn7U4790t3rMB4ZqCTVb+Q+9uYHlEPDH3GImIS7P35nfQyXhg7vl80rfXIv4OeBrwgojYF6g1MdWaeFqlBb4PeIwU3PLXvqfgtfN+Cexf17Q0v4NytrIaOE7Sc4BnkJqpekLSkaSg0HIIbDYKrNbxXOtXuRs4pNVhuee/yX6O5Lb9Xqfltd5zULBe+iRwlqQXKNlL0p9J2gf4LumD8n3Z9j0kvajFud6WdcYeCPwN0GzEy69I7f81+5Da538taQ7w9232H5d9Q/8MsFzSPpJGSe3kq1vedWM3AI8Db5Y0S9KrSE1oRcvZVERsJDVFrQI+34vmGEn7SjqW1M+zOiJu6eI0FwL/KOnQ7O9/uKS5jXaMiE2kYHtK1kF9Oq0Dik0RBwXrmYhYS+pXOB94AFgPvD57bzuwBFgI3EXqzGw13v3fgZuBHwBfYtdhkjXvBlZm7d2vIfU/7En61v8d4Oq6/T8KnJCN+PlYg/O9ifQt9g7St+VPkzrPOxIR24BXke7/AdK9fiG3S7tytrMSeDYpMEzGlZK2kL7lLwM+zMSO9U58mBRUrwUeIv3N9myx/5nA24DNwDOBb3d5XesheZEdGzSSAji0w3bwoSLpxaQazIKI2NHv8tj04ZqCWcVI2o3UpHahA4L1moOCWYVIegbwa9JQ2I/0tTA2Lbn5yMzMxrmmYGZm4yo9eW2//faLBQsW9LsYZmaVcvPNN98XEfMavTdQQUHSXsD1wN9HxFXt9l+wYAFr164tv2BmZtOIpA3N3iu1+UjSRZLulXRr3fbFkm6XtF7SObm33kEa52xmZn1Qdp/CJaTcLuOypGMXAMeQsiWeLOkwSS8HfkKacWpmZn1QavNRRFwvaUHd5qOA9bWVriRdRkqNvDewFylQPCLpyx6DbWY2tfrRp7A/E7MlbiQlBTsb0kImwH3NAoKkpcBSgPnz5zfaxczMutSPIamNskCOT5aIiEtadTJHxIqIWBQRi+bNa9h5bmZmXepHUNjIxLTIB1A8LTIAkpZIWvHggw92fPE1a2DBApgxI/1cs6bjU5iZTVv9CAo3AYdKOihbsvAk4IpOThARV0bE0tmzZ3d04TVrYOlS2LABItLPpUsdGMzMasoeknopKa/80yRtlHRGtvLW2cA1wG3AZyLix2WWo2bZMti6deK2rVvTdjMzq2juI0lLgCULFy48c926dYWPmzEj1RB2PR/s8DgnMxsSkm6OiEWN3qtk7qNum4+aDVbyICYzs6SSQaFby5fDyMjEbSMjabuZmVU0KHQ7+mhsDFasgNHR1GQ0Oppej42VVFAzs4qpZJ9CzaJFi8IJ8czMOjPt+hTMzKwclQwKk5m8ZmZmzVUyKHQ7+sjMzFqrZFAwM7NyDHVQcB4kM7OJBmo5zqJyM5q7PkctD1It7UUtDxJ4iKqZDa+hHZK6YEEKBPVGR+HOOydVLDOzgeYhqQ3cdVdn283MhsHQBgXnQTIz29XQBoVu8yDVOqclmDUr/XQntZlNF5UMCr2YvJbPgwQwc+bOtRXqP+DzgeDUU3f2RWzfnn5u2JC2O0CYWdVVMij0avLa2NjOGkP+Az6/Glt+tbZ07WZlany8mVmVDO3oo5pmo5AA5s6FBx7obgEej2Iys0Hl0UcttBpttHlz9yuybdjg2oKZVc/QB4UyRxu5GcnMqmbog0KjUUjtSOnnzJkTX9erdVybmVVFJYNCL1Nn149CamfmTFi1KnUsP/54+rlqVfP9N2zwiCQzq45KBoVep84eG0udwqtXt641jIzAypW75kYaG2sdVDwiycyqopJBoSz1azjPnZseRdZzbtcM5aYkM6uCSmZJLdPYWHdZUmvHLFvWfIir8yqZ2aBzTaGHas1QzZqSnFfJzAadg0IJGjUlSe50NrPB56BQgvoRTZLTYJhZNTgolCTflFSfSWTrVjjtNAcGMxs8Dgola9a5vH27awxmNngqGRR6OXmtbK06lz1M1cwGTSWDQq8nr5Wp3fwFD1M1s0HieQolq81fOO20nWs25HmYqpkNkkrWFKpmbCylx+hm+U8zs6nkoDBFGqXQ2HPPtIyn5y6Y2aBwUJhCtWGqq1bBI4+kRXwiPHfBzAaHg0IfLFuWRh7leSSSmQ0CB4U+aDbiyCORzKzfHBT6oNmIoxkz3IRkZv3loNAHzeYueJazmfWbg0If1EYi1dZ4znPfgpn1k4NCn4yNwY4djd9zim0z65eBCQqSniHp45I+J+mN/S7PVGg1m9nDVM2sH0oNCpIuknSvpFvrti+WdLuk9ZLOAYiI2yLiLOA1wKIyyzUovK6zmQ2asmsKlwCL8xskzQQuAI4BDgNOlnRY9t4rgW8BXyu5XAOhfjGeRjxM1cymUqlBISKuB+6v23wUsD4i7oiIbcBlwHHZ/ldExO8DY83OKWmppLWS1m7atKmsok8Zr+tsZoOkH30K+wN3515vBPaX9BJJH5P0CeDLzQ6OiBURsSgiFs2bN6/ssk4Zr+tsZoOgH6mz1WBbRMR1wHVTW5TBUUuxvWxZCgSN1nXO72dmVoZ+1BQ2AgfmXh8A/KKTE1Rp5bVOtFvX2Z3OZla2fgSFm4BDJR0kaXfgJOCKTk5QpZXXuuHcSGbWL2UPSb0UuAF4mqSNks6IiMeBs4FrgNuAz0TEjzs877SsKdQ061x2p7OZla3s0UcnR8STI2K3iDggIj6Vbf9yRDw1Ig6JiI7XHpvuNYVm8xceftgdzmZWroGZ0Ww71eYvzJ07cfvmzZ7lbGblqmRQmO7NR5ACw95777p961Y45RQPUzWzclQyKEz35qOaVh3Lzo1kZmVoOU9B0hEFzvFYRNzSo/JYzvz56cO/mdowVc9dMLNeaTd57RukIaSNJpzVHAQs6FWBbKfly1NtoH495zwPUzWzXmoXFG6KiJe22kHSf/awPIVIWgIsWbhw4VRfekrVz3JupLaEp2sLZtYLLfsU2gWEovv02rD0KcDOWc6rV3sJTzMrX6GOZkmL2+9lZfISnmY2FVoGBUlHZesfvDe3bVXppbKGvISnmZWtXU3hZFJn88GS3i/ptUCREUmlGoZ5Cs14CU8zK1O7oPCOiPgD4C7gS8Ac4PckfUfS5aWXrolh6lOo5yU8zaxM7UYfXSPpcWAeMJe0+M3pEXG0pANKL53tosiIJA9TNbNutRt99MfAK4GHgUOA9wALJX2RlPLa+sBLeJpZWdqOPoqIR4C7I+KDEfFa4GfAmcDPyy6cteZsqmbWa4WGpNbNRTg/IjZFxOdLKlNbw9zRnOdsqmbWa4r6dR8rZNGiRbF27dp+F6PvFixo3L8wOpqamczM8iTdHBGLGr3Xbp7CigInb7uPlcvLd5pZr7QbfXS8pN+2eF/AH/ewPNaFZtlU589PTUjLlqUAMX9+6odwniQza6ZdUHhbgXN8sxcFse41yqYqpUBx6qlQayGsTW4DBwYza6xlUIiIlVNVEOte/dwFaWcgqO8y8hoMZtZKJVde8+ijXeXnLrQbO+C+BjNrppJBYZjTXLRT5AO/tgaDmVm9tkFB0kxJH5iKwtjkFZnN7DUYzKyZIjOatwPPl9RqSU4bEO0S5tU4cZ6ZNVK0+ej7wL9LOlXSq2qPMgtm3anNch4dTR3OzfIjgddgMLNdtRuSWjMH2Azk010E8IWel8gmbWxs4uiiZjOewcNUzWwip7kYAmvW7DqPoZ5TYpgNj67TXORO8FRJX5N0a/b6cEnv7GUhrTz5JqVmNmyA/fZLI5PcpGQ2vIr2KXwSOBd4DCAifoTXU6iUdmswQMquGuFlPc2GWdGgMBIR363b9nivC1OUJ691z6OTzKyVokHhPkmHkDqXkXQC8MvSStWGJ691r9aUVIRnPpsNn6JB4a+BTwBPl3QP8BbgrLIKZeUaG2vdjFTjZT3Nhk/RldfuiIiXA/OAp0fEH0REk0GOVgXtmpFGRtI+ZjZcio4+mivpY6Q02ddJ+qikue2Os8FVP8lt7tz0qE14W7HC8xbMhlHRyWuXAdcDr85ejwGXAy8vo1A2NeonuZmZFe1TmBMR/xgRP88e7wGeWGK5rE/WrEnzFDxfwWw4FQ0KX5d0kqQZ2eM1wJfKLJhNvdrM5w0bds5XOPXU1KTkAGE2HFqmuZC0hTQMVcBewI7srRnAwxGxb+klbMFpLnqrVY4kSJ3P7mswq76u01xExD4RsW/2c0ZEzMoeM/odEKz32s1L2LoVTjkFZs1y7cFsuira0Yykw4EF+WMiwllSp5H581vXFGq2b08/nWHVbPopOiT1IuAi0uijJdnj2BLLZX1QNAVGntNhmE0vRWsKR0fEYaWWBJB0PPBnwO8AF0TEtWVf03aqfdtftizVAqTU4dyO02GYTR9FRx/dIKmroCDpIkn31tJu57YvlnS7pPWSzgGIiC9GxJnA64ETu7meTU4tm2oErFrldBhmw6ZoUFhJCgy3S/qRpFsk/ajgsZcAi/MbJM0ELgCOAQ4DTq4LOu/M3rc+qgWI1aubNytJqVbhzmez6aFo89FFwKnALewcllpIRFwvaUHd5qOA9RFxB4Cky4DjJN0GvA/4SkR8r9H5JC0FlgLM91fUKVHfrDRzZupszjcvufPZbHooWlO4KyKuyGYzb6g9JnHd/YG7c683ZtveREqdcYKkhllYI2JFRCyKiEXz5s2bRBGsE/lmpccfT81Kzfob3PlsVl1Fawo/lfRp4Erg0drGSQxJVYNtEREfAz7W5TltCrXrXHbns1k1Fa0p7EkKBq+gN0NSNwIH5l4fAPyi6MFeea3/2rXcRbh/wayKWqa56NlFUp/CVRHxrOz1LOC/gZcB9wA3Aa+NiB93cl6nueifWp6krVtb7+fUGGaDp1Wai0LNR5IuJluKMy8iTi9w7KXAS4D9JG0E/j4iPiXpbOAaYCZwUScBQdISYMnChQuLHmI91qzzuV6tf8FBwawaCtUUJL0693IP4M+BX0TEm8sqWBGuKQyWGTOadz6PjqZ+hjlz0uv7709NUMuXO2CYTbVJ1xQi4vN1J7wU+I8elM2mkVa5k2rbN2+euM3DV80GS9GO5nqHAn2bJOCO5sHk3Elm1Vc0Id4WSQ/VHqShqe8ot2jNRcSVEbF09uzZ/SqCNZBf97kTHr5qNjiKNh/tU3ZBbHqorfvcqn+hniemmw2OlkFBUsv/rhHh73jWUNG1GUZGUrOTmQ2GdjWFL7FzOc6aAOaR0lvPLKlcLXlI6uBbvrzxPIa99oI99vDoI7NB1W45zmdHxOHZz2eTZjL/F/Aw8JYpKF+zcrlPYcDl+xek9HP1anj4YbjvvpSWG+DUUz3z2WyQFJ2ncCiwDHgB8CFgZUQ8VnLZ2vI8hWpqNBu6lnF1dNS1B7OydT1PQdKzSMHgmcB5wBkR0WDeqllxy5bt2qxU+27iuQtm/dWypiBpOynF9ZeAXYJBv2Y05/oUzly3bl0/imCTUGRk0uhoStVtZr03mRnNbXMb9UNEXAlcuWjRojP7XRbrXJGRSZ67YNYfLYNCRKycqoLY8Gg2MilvzpzUAX3XXR6lZDaVuk1zYda1+pnPqltyabfdYMuWVJuI2NnP4BFKZuVzULC+yC/vuWrVxKGr++4L27ZN3N85ksymRtHcRy8qsm2qOCHe9FILEDt2pJ/33994P/czmJWvaE3hnwtumxKevDa9NcuF5BxJZuVrN0/hhcDvA/Mk/W3urX3pU4oLm/4adUQ7R5LZ1GhXU9gd2JsUPPbJPR4CTii3aDas6lNkzJ0Le+7plBhmU6HdkNRvAN+QdElEFMh5ada9NWtSZ3J+GCpMrDV4xrNZuYrmPnoq8FZgAblAEhEvLa1kBTj30fTRKB/SyEiqIeSX8KzxjGez7k16jWbgs8DHgQtpkO5iqjl19vTTKB/S1q3NJ7h5JJJZOYoGhccj4l9LLUkHnOZi+un0Q94jkczKUXRI6pWS/krSkyXNqT1KLZkNlWYf8nPnpmakPI9EMitP0aBwGvA24NvAzdnDjfnWM8uXN/7w/+hHd12sZ8WK1p3Ma9akUUozZni0klmnCjUfRcRBZRfEhlvtQ75+9FFte9GRRvUd1h6tZNaZomkuRiS9U9KK7PWhko4tt2g2bOrTXdR/iBepATTrsHbeJLNiijYfXQxsI81uBtgIvKeUEpk1UKsBtMuc2qzD2qOVzIopGhQOiYjzgMcAIuIRQK0PMeudojUA500ym5yiQWGbpD2BAJB0CPBoaaUyq9Psm/6GDRObkpp1WHu0klkxRYPCu4GrgQMlrQG+Bry9rEK149TZw6fVN/18U5LzJplNTqE0FwCS5gJHk5qNvhMR95VZsCKc5mJ4NEqDUa8+9UWz1BnthrSaTXet0lwUHX10BfAK4LqIuGoQAoINl/olPBupb2LySCSzzhVtPvoQ8IfATyR9VtIJkvYosVxmu6gNWW0WGObMmThkdUOTvL4eiWTWXKGgEBHfiIi/Ag4GVgCvAe4ts2BmzTTqTN5tN9iyZeKQVTUZH+eRSGbNFa0pkI0+ejVwFnAksLKsQpm1Ut+ZPDoK++4L27ZN3C9i18DgkUhmrRXtU7gcuA14KXABad7Cm8osmFkr9bOf77+/8X4RneVNMht2RVNnXwy8NiL6vpaCWSPz5zfuQ/BiPGadaVlTkPR2gIi4GnhV3XvvLbFcZh3xpDWz3mjXfHRS7vm5de8t7nFZzLrWqJ/BTUVmnWvXfKQmzxu9Nptya9Y0T7dtZp1rFxSiyfNGr82mlNdOMOu9ds1Hz5H0kKQtwOHZ89rrZ/eyIJIOlvQpSZ/r5Xlt+up2xrJXZjNrrmVQiIiZEbFvROwTEbOy57XXu7U7uaSLJN0r6da67Ysl3S5pvaRzsmvdERFnTO52bJh0s3ZC0XUZzIZV4clrXbqEug5pSTNJcx2OAQ4DTpZ0WMnlsGmom7UTitQuXJOwYVZqUIiI64H6aUVHAeuzmsE24DLguKLnlLRU0lpJazdt2tTD0lrVdDMMtV3twjUJG3Zl1xQa2R+4O/d6I7C/pLmSPg48T1L98NdxEbEiIhZFxKJ58+aVXVYbYN0MQ21Xu3BmVRt2RWc091KjoawREZtJeZXan0BaAixZuHBhTwtm1TM21tlIo+XLG6+xUKtdeI1nG3b9qClsBA7MvT4A+EUnJ4iIKyNi6ezZs3taMJv+2tUuvMazDbt+BIWbgEMlHSRpd9Ks6Sv6UA4bUvXJ9PI1DafLsGFXalCQdClwA/A0SRslnRERjwNnA9eQMq9+JiJ+3OF5vUazdazIqCKny7BhV3iN5kHkNZqtKK/XbLbTpNdoNqs6jyoyK6aSQcHNR9apIvMTPGHNrKJBwaOPrFOtRhV5wprZTpUMCmadajWqyE1LZjtVMii4+cg61WpUkSesme1UyaDg5iPrRqP5CWvWpH6ERvJNTu5zsGHRjzQXZgOh1pewffuu7+UnrHkxHxsmnqdgQ2vBgvQBX2/mTFi5cucHfrP9RkdTjcOsaqbdPAX3KVgvNOsz2LFjYg3AfQ42TCoZFNynYL1QNPmdk+TZMKlkUDDrhSLJ79asgYcf3vXYRvt12xHtTmwbKBFR2cfzn//8MJuM1asjRkcjpPRz9eqJ742MRKQpbTsfc+e2329kZOI+ra7f7bFm3QLWRpPPVXc0mzVRtIN5Mh3R7sS2fnBHs1kXinYwT6Yj2p3YNmgqGRTCHc02BaaiI9qd2DZoKhkUzKZC0VXYJrNam1d6s0HjoGDWRNFV2CazWptXerNB4zQXZi2Mje36Ab1mTcqgumFDmv28fXv6MF++vP2Hee3Yu+6COXPStvvvT81Fq1Y5GFj/OSiYdaA+D1Itb1KRfEj1x27evPM951OyQeEhqWYdaDaEtKbVUNJ2x7Y73qxXPCTVrEfaDRVt9f5khqiaTZVKBgUPSbV+aTdUtNX7kxmiajZVKhkUzPql0RDSmnZDSVsdW+R4s6ngoGDWgfwQUkijj6DYUNL64adz56aHh6LaIHFHs5nZkJl2Hc1mZlYOBwWzHurn2gidXLvMcnp9iIprllO7Cg+vp2CDpJ9rI3Ry7TLL6fUhqgGvp2BWvn6ujdDJtcssp9eHqIZWfQqVDAqSlgBLFi5ceOa6dev6XRwzIDWXNPrvJMGOHYNz7TLL2c/fgRU37Tqaw5PXbAD1c22ETq5dZjm9PkT1VTIomA2ifq6N0Mm1yyyn14eoPgcFsx7p59oInVy7zHJ6fYjqq2SfQo07ms3MOjft+hTMzKwcDgpmZjbOQcHMzMY5KJiZ2TgHBTMzG+egYGZm4xwUzMxs3MAEBUl7SVop6ZOSPNXFKm+yKaQbHV/bJsGsWennfvulR7PrtCpHuzK2KkO7bd1cO7+t3X0VvYde/S4me+1eKf2azdKn9uIBXATcC9xat30xcDuwHjgn23YqsCR7fnmR8zt1tg2qyaaQbnT8brtF7L77xG2NHvnrtCpHuzIWLUOjbSMjEW98Y2fXbnd/jX5/vUoZ3s3fqx9pwnt1TVqkzi47KLwYOCIfFICZwM+Ag4HdgR8ChwHnAs/N9vl0kfM7KNigGh1t/ME2Ojq544s+atdpVY52ZZxsGWbO7PzaRe+rm9/zZH4XZfyNu9Gra7YKCqWnuZC0ALgqIp6VvX4h8O6I+JPs9bnZrhuBByLiKkmXRcRJTc63FFgKMH/+/OdvaJS83azPJptCutnxRdWu06oc0LqMky1Dq7I1u3aRY/O/v16lDG9WnlZ/r36kCe/VNQctzcX+wN251xuzbV8AXi3pX4Ermx0cESsiYlFELJo3b165JTXr0mRTSE821XTt+FblaFfGyZZh5szOr91O/XG9Shnezd+rH2nCp+Ka/QgKarAtIuI3EfEXEfHGiGjZdSJpiaQVDz74YElFNJucyaaQbnT8brvB7ru3PzZ/nVblaFfGomVotG1kBJYu7eza7e6v0e+vVynDu/l79SNN+JRcs1m7Uq8ewAIm9im8ELgm9/pc4Nxuzu0+BRtkq1entl4p/ey0M7DR8bVt+Tb7uXPTo9l1WpWjXRlblaHdtm6und/W7r6K3kOvfheTvXav9OKaDFifwizgv4GXAfcANwGvjYgfd3pup842M+tc3/oUJF0K3AA8TdJGSWdExOPA2cA1wG3AZzoNCG4+MjMrhxfZMTMbMoM2+sjMzAZUJYOCm4/MzMpRyaAQEVdGxNLZs2f3uyhmZtNKpfsUJG0Cup3SvB9wXw+L02/T6X6m072A72eQTad7geL3MxoRDWf/VjooTIaktc06WqpoOt3PdLoX8P0Msul0L9Cb+6lk85GZmZXDQcHMzMYNc1BY0e8C9Nh0up/pdC/g+xlk0+leoAf3M7R9CmZmtqthrimYmVkdBwUzMxs3lEFB0mJJt0taL+mcfpenE5IOlPR1SbdJ+rGkv8m2z5H0VUnrsp9P6ndZOyFppqTvS7oqe13J+5H0REmfk/TT7G/0wqreC4Ck/5X9O7tV0qWS9qjS/Ui6SNK9km7NbWtafknnZp8Lt0v6k/6UurEm9/KB7N/ajyT9P0lPzL3X1b0MXVCQNBO4ADiGtDb0yZIO62+pOvI48HcR8QzgaOCvs/KfA3wtIg4Fvpa9rpK/IWXNranq/XwUuDoing48h3RPlbwXSfsDbwYWZanvZwInUa37uQRYXLetYfmz/0cnAc/MjvmX7PNiUFzCrvfyVeBZEXE4aUmCc2Fy9zJ0QQE4ClgfEXdExDbgMuC4PpepsIj4ZUR8L3u+hfShsz/pHlZmu60Eju9LAbsg6QDgz4ALc5srdz+S9gVeDHwKICK2RcSvqeC95MwC9szWQRkBfkGF7icirgfur9vcrPzHAZdFxKMR8XNgPenzYiA0upeIuDZbjgDgO8AB2fOu72UYg0KzNaIrJ1vA6HnAjcDvRsQvIQUO4Hf6WLROfQR4O5BferyK93MwsAm4OGsKu1DSXlTzXoiIe4APAncBvwQejIhrqej95DQrf9U/G04HvpI97/pehjEoNFwjespLMUmS9gY+D7wlIh7qd3m6JelY4N6IuLnfZemBWcARwL9GxPOA3zDYTSstZW3txwEHAU8B9pJ0Sn9LVarKfjZIWkZqWq6tb9/1vQxjUNgIHJh7fQCpSlwZknYjBYQ1EfGFbPOvJD05e//JwL39Kl+HXgS8UtKdpKa8l0paTTXvZyOwMSJuzF5/jhQkqngvAC8Hfh4RmyLiMeALwO9T3fupaVb+Sn42SDoNOBYYi50Tz7q+l2EMCjcBh0o6SNLupM6YK/pcpsIkidRmfVtEfDj31hXAadnz04B/n+qydSMizo2IAyJiAelv8Z8RcQoVvJ+I+B/gbklPyza9DPgJFbyXzF3A0ZJGsn93LyP1YVX1fmqalf8K4CRJT5B0EHAo8N0+lK8wSYuBdwCvjIitube6v5eIGLoH8KeknvqfAcv6XZ4Oy/4HpGrgj4AfZI8/BeaSRlKsy37O6XdZu7i3lwBXZc8reT/Ac4G12d/ni8CTqnov2f38H+CnwK3AKuAJVbof4FJSf8hjpG/PZ7QqP7As+1y4HTim3+UvcC/rSX0Htc+Cj0/2XpzmwszMxg1j85GZmTXhoGBmZuMcFMzMbJyDgpmZjXNQMDOzcQ4K1heStkv6QZZ987OSRnpwzn+Q9PI2+5wl6XXZ89dLekqB807YL0tfMekkipLeLekeSf8w2XO1uMbVWWK7/LaX1LLRdnG+E7PMm10db4PPQcH65ZGIeG6k7JvbgLPyb3aTnTIi3hUR/9Fmn49HxL9lL19PSt/QzoT9IuINEfGTTsvXxD9FxLt6dK4JJO1JGoN/T6/OGRGXA2/o1fls8Dgo2CD4JrAw+wb7dUmfBm5RWmPhA5JuyvLF/2XtAElvl3SLpB9Kel+27RJJJ2TP75T0fknfzR4Ls+3vlvTWbL9FwJqsxrKnpHdl17pV0goljfa7TtKi7HwnZ+W4VdL7c+V7WNLyrHzfkfS77X4JWdlWSro2K/+rJJ2Xnf/qLL1J7d7eK+kGSWslHSHpGkk/k5QPri8BrsuOWayUd/9bwKty1zxK0reVEvh9uzYbW9I3JT03t99/STq8+J/UqspBwfpKKSXzMcAt2aajSLPMDyPN2HwwIo4EjgTOzNKTHENKd/yCiHgOcF6T0z8UEUcB55MysY6LiM+RZh6PZTWWR4DzI+LIrPayJ3Bsk/1qZX8K8H7gpaSZzEdKOj57ey/gO1n5rgfOLPgrOYSURvw4YDXw9Yh4NvBItr3m7oh4ISmgXgKcQFpfI98UdQxwtaQ9gE8CS4A/BH4vt89PgRdHSuD3LuC92fYLSTUkJD0VeEJE/KjgPViFOShYv+wp6QekD9y7yNYgAL4bKf87wCuA12X73UhKT3AoKVHbxZHleomI+nz5NZfmfr6wQJn+WNKNkm4hfdA/s83+RwLXRUoYV8tQ+eLsvW1Ard39ZmBBgesDfCVS8rlbSIvaXJ1tv6XuHFfktt8YEVsiYhPwW+1cfetFwLeAp5MS262LlMJgde48s4HPKq3m9U/svOfPAsdmtZPTSYHHhsCsfhfAhtYjEfHc/IaUc43f5DcBb4qIa+r2W0yxNMDR5Pkusm/T/0JaZexuSe8G9mhz/kbpiWsei505ZLZT/P/aowARsUNS/hw76s7xaG77o7ntO4BZkg4m1Sa2Zb/XZvf/j6TayJ8rrc9xXXb9rZK+SqqxvIbUhGZDwDUFG2TXAG/MtaU/VWnRmmuB05WNWJI0p8nxJ+Z+3tDg/S3APtnzWgC4T2mtihOa7Jd3I/BHkvbLOsZPBr5R6M7Kdww7axk/BQ6SdEj2+uTcfrOBWkf06+vOcSHwMeCmFrUxm2ZcU7BBdiGpyeR7Sl93NwHHR8TVWSfoWknbgC8D/7vB8U+QdCPpy8/JDd6/BPi4pEdIzUufJDXH3ElKsd5sPyCt2iXpXODrpFrDlyNiUNJILwbeBBARv5W0FPiSpPtITUrPyvY7D1gp6W+B/8yfICJulvQQcPHUFdv6zVlSbVpSWrRnUUTc1++yNJM1UT0cER/s8XmfAPxXREyqySfrSL8OeHpE7Mhtfwnw1og4djLnt8Hk5iOz/nkYWKoeT16LtFj7ZAPC60jNY8vqAsKJpL6XByZXShtUrimYmdk41xTMzGycg4KZmY1zUDAzs3EOCmZmNs5BwczMxv1/LTCwwywZZJAAAAAASUVORK5CYII=\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, 'bo')\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-cru')\n",
    "plt.savefig(\"/home/mkibria/Academic/pr_of_thesis/figures/pdf_cru.png\")\n",
    "\n",
    "\n",
    "\n",
    "plt.show( )"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "a43f36ef",
   "metadata": {},
   "outputs": [],
   "source": [
    " \n",
    " "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "2f5defd5",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "b3c21546",
   "metadata": {},
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "0f0c6eb5",
   "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
}