tutorials/building_a_controller.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Building a Controller"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The following documents the development of a new controller.\n",
    "In this case we are going to implement an arbitrary controllable storage unit. This\n",
    "may be a battery, an electrically powered car or some sort of reservoir storage."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Modelling a Battery\n",
    "\n",
    "In order to simulate a storage system we use the storage element of pandapower. The power of a storage can be positive or negative. To overcome this issue, a storage element can be created. \n",
    "\n",
    "For storage elements the signing is based on the consumer viewpoint (positive active power means power consumption and therefore charging of the battery).\n",
    "\n",
    "As pandapower is not a time dependend simulation tool and there is no time domain parameter in default power flow calculations, the state of charge (SoC) is not updated during any power flow calculation. \n",
    "In Order to update the SoC we build our own storage controller and keep track of the SoC.\n",
    "\n",
    "State of charge (SoC [\\%]) is the level of charge of an electric battery relative to its capacity.\n"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Controller init\n",
    "First we start by creating a new file *control/storage_control.py*, containing our new class."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Import and inherit from the parent class `Controller` and override methods you would like to use. Next we write the actual code for the methods."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "import pandapower as pp\n",
    "import control\n",
    "import pandas as pd\n",
    "import timeseries as ts\n",
    "\n",
    "# importing a grid from the library\n",
    "from pandapower.networks import mv_oberrhein\n",
    "\n",
    "\n",
    "class Storage(control.basic_controller.Controller):\n",
    "    \"\"\"\n",
    "        Example class of a Storage-Controller. Models an abstract energy storage.\n",
    "    \"\"\"\n",
    "    def __init__(self, net, gid, data_source=None, p_profile=None, in_service=True,\n",
    "                 recycle=False, order=0, level=0, **kwargs):\n",
    "        super().__init__(net, in_service=in_service, recycle=recycle, order=order, level=level,\n",
    "                    initial_powerflow = True, **kwargs)\n",
    "        \n",
    "        # read generator attributes from net\n",
    "        self.gid = gid  # index of the controlled storage\n",
    "        self.bus = net.storage.at[gid, \"bus\"]\n",
    "        self.p_mw = net.storage.at[gid, \"p_mw\"]\n",
    "        self.q_mvar = net.storage.at[gid, \"q_mvar\"]\n",
    "        self.sn_mva = net.storage.at[gid, \"sn_mva\"]\n",
    "        self.name = net.storage.at[gid, \"name\"]\n",
    "        self.gen_type = net.storage.at[gid, \"type\"]\n",
    "        self.in_service = net.storage.at[gid, \"in_service\"]\n",
    "        self.applied = False\n",
    "\n",
    "        # specific attributes\n",
    "        self.max_e_mwh = net.storage.at[gid, \"max_e_mwh\"]\n",
    "        self.soc_percent = net.storage.at[gid, \"soc_percent\"] = 0\n",
    "        \n",
    "        # profile attributes\n",
    "        self.data_source = data_source\n",
    "        self.p_profile = p_profile\n",
    "        self.last_time_step = None\n",
    "        \n",
    "    # We choose to represent the storage-unit as a storage element in pandapower. \n",
    "    # We start with a function calculating the amout of stored energy:    \n",
    "    def get_stored_ernergy(self):\n",
    "        # calculating the stored energy\n",
    "        return self.max_e_mwh * self.soc_percent / 100        \n",
    "    \n",
    "    # convergence check\n",
    "    # Also remember that 'is_converged()' returns the boolean value of convergence:\n",
    "    def is_converged(self):\n",
    "        # check if controller already was applied\n",
    "        return self.applied\n",
    "    \n",
    "    # Also a first step we want our controller to be able to write its P and Q and state of charge values back to the\n",
    "    # data structure net.\n",
    "    def write_to_net(self):\n",
    "        # write p, q and soc_percent to bus within the net\n",
    "        self.net.storage.at[self.gid, \"p_mw\"] = self.p_mw\n",
    "        self.net.storage.at[self.gid, \"q_mvar\"] = self.q_mvar\n",
    "        self.net.storage.at[self.gid, \"soc_percent\"]= self.soc_percent\n",
    "        \n",
    "    # In case the controller is not yet converged, the control step is executed. In the example it simply\n",
    "    # adopts a new value according to the previously calculated target and writes back to the net.\n",
    "    def control_step(self):\n",
    "        # Call write_to_net and set the applied variable True\n",
    "        self.write_to_net()\n",
    "        self.applied = True\n",
    "        \n",
    "    # In a time-series simulation the battery should read new power values from a profile and keep track\n",
    "    # of its state of charge as depicted below.\n",
    "    def time_step(self, time):\n",
    "        # keep track of the soc (assuming time is given in 15min values)\n",
    "        if self.last_time_step is not None:\n",
    "            # The amount of Energy produce or consumed in the last timestep is added relative to the \n",
    "            # maximum of the possible stored energy\n",
    "            self.soc_percent += (self.p_mw * (time-self.last_time_step) * 15 / 60) / self.max_e_mwh * 100\n",
    "        self.last_time_step = time\n",
    "\n",
    "        # read new values from a profile\n",
    "        if self.data_source:\n",
    "            if self.p_profile is not None:\n",
    "                self.p_mw = self.data_source.get_time_step_value(time_step=time,\n",
    "                                                                profile_name=self.p_profile)\n",
    "                \n",
    "        self.applied = False # reset applied variable"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We are now ready to create objects of our newly implemented class and simulate with it!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "hp.pandapower.control.util.auxiliary - INFO: Creating controller 0 of type <class '__main__.Storage'> \n"
     ]
    }
   ],
   "source": [
    "# importing a grid from the library\n",
    "from pandapower.networks import mv_oberrhein\n",
    "\n",
    "# loading the network with the usecase 'generation'\n",
    "net = mv_oberrhein()\n",
    "pp.runpp(net)\n",
    "\n",
    "# creating a simple time series\n",
    "framedata = pd.DataFrame([0.1, .05, 0.1, .005, -0.2, 0], columns=['P'])\n",
    "datasource = ts.DFData(framedata)\n",
    "\n",
    "# creating storage unit in the grid, which will be controlled by our controller\n",
    "store_el = pp.create_storage(net, 30, p_mw = .1, q_mvar = 0, max_e_mwh = 0.1, )\n",
    "\n",
    "# creating an Object of our new build storage controller, controlling the storage unit\n",
    "ctrl = Storage(net=net, gid=store_el, data_source=datasource,p_profile='P')"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we run a small time-series-simulation and track the results using the outputwriter:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Progress: |██████████████████████████████████████████████████| 100.0% Complete\n",
      "\n"
     ]
    }
   ],
   "source": [
    "# defining an OutputWriter to track certain variables\n",
    "ow = ts.OutputWriter(net)\n",
    "\n",
    "ow.log_variable(\"res_storage\", \"p_mw\")\n",
    "ow.log_variable(\"storage\", \"soc_percent\")\n",
    "\n",
    "# starting time series simulation\n",
    "ts.run_timeseries(net, time_steps=range(0, 6), output_writer=ow)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "To visualize the results we plot directly with the dataframe:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.legend.Legend at 0x207346c0dd8>"
      ]
     },
     "execution_count": 6,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX4AAAEGCAYAAABiq/5QAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO3dd3xUdfb/8ddJI/QSQk1C7z0kseDiuih2sCFBcXXVr6Ku4s91FV3rrutXXXVta2HdglIFVLAhimL5Wkih9yKQQIBQEgKk5/z+mItGDGGAzNwp5/l4zGPmTmbufU+Ukzv3fu75iKpijDEmfES4HcAYY4x/WeE3xpgwY4XfGGPCjBV+Y4wJM1b4jTEmzES5HcAbLVu21I4dO7odwxhjgkpWVtZuVY0/8vmgKPwdO3YkMzPT7RjGGBNURGRLTc/boR5jjAkzVviNMSbMWOE3xpgwExTH+GtSXl5Obm4uJSUlbkcJSLGxsSQkJBAdHe12FGNMgAnawp+bm0vjxo3p2LEjIuJ2nICiquzZs4fc3Fw6derkdhxjTIAJ2kM9JSUlxMXFWdGvgYgQFxdn34aMMTUK2sIPWNGvhf1ujDFHE7SHeowxoaWqSikqrWB/cTn7S8opLC5nf3EF+0vK2V9cTnmlcu3pHWgQY2XrZNlv8CRERkbSr18/VJXIyEheeuklTj/99KO+vqCggKlTp3LrrbcCsHnzZr755huuuuoqf0U2xmdUleLySvYXV3iKtlOwPQW8nP0lFT8tl5T/4nVFpRUca3qQmKgIbjjDzludLCv8J6F+/fosWbIEgI8//pj77ruPL7744qivLygo4OWXX/5Z4Z86depxFf7KykoiIyNPLrgxR1FWUfXLol1y9L3w/Uf8vLyy9srdICaSpvWjaRIbTZP6UbRrFkvP2MY0qR/tucVGOffRntfVj/I8bhDNNf9axPRFW7l+iA3oOFlW+OvI/v37ad68OQAHDhxg5MiR7Nu3j/Lych577DFGjhzJhAkT2LhxIwMHDuScc87hq6++YvXq1QwcOJBrr72WO+64gwkTJrBw4UJKS0u57bbbuPnmm1m4cCGPPvoobdu2ZcmSJVx++eW0bNmS8ePHA/CnP/2J1q1bc8cdd7j5KzABoLJKKSr5qTj/tLf986J9tL3wkvKqWtcfExnhFOkomtaPplmDGJLiGtIkNsop1DUUbef5xrFRREee+GnFMamJTHh7OdlbCxjcofkJr8eESOF/9L2VrNq+v07X2btdEx6+uE+trykuLmbgwIGUlJSQl5fHZ599BnjG0L/zzjs0adKE3bt3c+qppzJixAieeOIJVqxY8eO3hIULF/L000/z/vvvAzBx4kSaNm1KRkYGpaWlDBkyhOHDhwOwaNEiVqxYQadOndi8eTOXXXYZ48ePp6qqiunTp7No0aI6/fwmcJVWVPLsJ+vYlH/wx6Jd5BTwotKKWt8bIfyiOLdq3Kha0fbscVffK//pcTSx0e5927x4QDv+8v4qpi/aaoX/JIVE4XdL9UM93377Lb/97W9ZsWIFqsr999/Pl19+SUREBNu2bWPnzp3HXN/8+fNZtmwZs2bNAqCwsJD169cTExNDWlraj2PyO3bsSFxcHIsXL2bnzp0MGjSIuLg4331QE1AembuKaYu20rON5xBJQvMGNe5h/2wv3FluVC8qaA+TNKwXxcUD2jFnyXYeurg3jWPt4sQTFRKF/1h75v5w2mmnsXv3bvLz8/nwww/Jz88nKyuL6OhoOnbs6NWYelXlxRdf5Nxzz/3Z8wsXLqRhw4Y/e+7GG2/kv//9Lzt27OD666+v089iAtdbGTlMW7SVcWd2YcL5Pd2O43fpaUlMz8hh7tLtXH1KB7fjBK2gHscfSNasWUNlZSVxcXEUFhbSqlUroqOj+fzzz9myxdMZtXHjxhQVFf34niOXzz33XF555RXKy8sBWLduHQcPHqxxe5deeinz5s0jIyPjF38oTGhallvAA3NWMKRrHHcP7+52HFcMSGhKzzaNmZGR43aUoBYSe/xuOXyMHzx765MmTSIyMpKrr76aiy++mJSUFAYOHEjPnp49s7i4OIYMGULfvn05//zzefzxx4mKimLAgAFcd911jB8/ns2bN5OcnIyqEh8fz7vvvlvjtmNiYjjrrLNo1qyZjfIJA3sPlnHL5GxaNozhhfRBRJ3ESdJgJiKkpybyyHurWLm9kD7tmrodKSiJHmvgbABISUnRIydiWb16Nb169XIpkfuqqqpITk5m5syZdOvWrcbXhPvvKFRUVinX/nsRi37Yy8xxpzEgsZnbkVxVcKiMtMcXkJ6ayJ9H9nU7TkATkSxVTTnyeZ/uNohIMxGZJSJrRGS1iJwmIi1E5BMRWe/c2+n547Rq1Sq6du3KsGHDjlr0Teh4Zv5avt6wmz+P7BP2RR+gWYMYLujbhncWb6O4rNLtOEHJ198XnwfmqWpPYACwGpgALFDVbsACZ9kch969e7Np0yaeeeYZt6MYH/t45Q5eXriR9NRE0tOS3I4TMEanJlFUUsFHK/LcjhKUfFb4RaQJMBT4F4CqlqlqATASmOS8bBJwyYluIxgOU7nFfjfBb2P+Af7w1lL6JzTlkRHuj1wLJKd2bkHHuAZMX2QneU+EL/f4OwP5wH9EZLGIvC4iDYHWqpoH4Ny3qunNInKTiGSKSGZ+fv4vfh4bG8uePXuswNXgcD/+2NhYt6OYE3SwtIJxb2YRExXBK2MHu3rhVCASEUanJrFo81425h9wO07Q8eWoniggGbhdVb8Xkec5jsM6qjoRmAiek7tH/jwhIYHc3Fxq+qNgfpqBywQfVeWeWcvYmH+AN284hfbN6rsdKSBdPrg9z8xfy4yMHO6/wAYxHA9fFv5cIFdVv3eWZ+Ep/DtFpK2q5olIW2DXiaw8OjraZpcyIen1r37gg+V53HteT4Z0bel2nIDVqnEsw3q1YnZWLncP70FMVHgOcT0RPvtNqeoOIEdEejhPDQNWAXOBa53nrgXm+CqDMcHm2417eGLeGs7r04ZxZ3Z2O07AS09LYs/BMj5dfeyWKOYnvr6A63ZgiojEAJuA3+H5Y/OWiNwAbAVG+TiDMUEhr7CY30/NpmNcA/42qn/Q9tTxp6Hd4mnXNJbpGTlc0K+t23GChk8Lv6ouAX5x8QCevX9jjKO0opJbJmdTUl7Ja9ecag3IvBQZIYxKSeSFz9aTs/cQiS0auB0pKNhBMWMCwJ/fW8WSnAL+NmoAXVs1djtOULkyNRGAmVm5LicJHlb4jXHZzMwcpny/lZvP7GyHK05A+2b1GdotnpmZOVRW2fBub1jhN8ZFK7YV8qd3V3B6lzj+OLzHsd9gajQmLZG8whK+XGfDu71hhd8Yl+w7WMbNb2bRsmEML44J346bdeE3PVvTslEM0xZtdTtKULD/04xxQWWVcsf0xeQXlfLK2MHENarndqSgFhMVweWDE1iwZhe7io496VG4s8JvjAv+/sk6vlq/m0et42adGZ2SSGWVMstO8h6TFX5j/Gz+yh289PkGRqckMsY6btaZzvGNOKVTC2Zk5FBlJ3lrZYXfGD/aVK3j5qMjreNmXUtPS2TLnkN898Met6MENCv8xvjJwdIKxk3OIipSePnqZOu46QPn921Lk9goa9d8DFb4jfEDVeXe2cvYsOsAL45JJqG5XWHqC7HRkVw6qD3zVuxg38Eyt+MELCv8xvjBv77+gfeX5XH3uT04o5t13PSl9LQkyiqreGfxNrejBCwr/Mb42Heb9vC/H63h3D6tueXMLm7HCXm92jZhQEJTZmTk2ERNR2GF3xgf2lFYwu+nZtMhrgFPjxpgHTf9JD0tibU7i1icU+B2lIBkhd8YHymrqOKWKVkcKqvktbGDreOmH108oB0NYiKZYSd5a2SF3xgf+cv7q1i8tYC/XTGAbq2t46Y/NaoXxcX92/Hesu0cKK1wO07AscJvjA/Mzsrlze+2cNPQzlzY3zpuumF0WiKHyip5b+l2t6MEHCv8xtSxFdsKuf+d5ZzWOY57zrWOm24ZlNiMHq0bM90at/2CFX5j6lDBoTJumZJFi4YxvHiVddx0k4gwOjWRpbmFrNq+3+04AcX+rzSmjlRWKeOnL2FnYSkvX51MS+u46brLktsTExXBjAzb66/OCr8xdeT5T9fxxbp8Hh7Rm0FJzd2OY4BmDWI4r08b3lm8jZLySrfjBAwr/MbUgU9X7eSFzzYwanACV1nHzYCSnpbI/pIKPlqR53aUgGGF35iT9MPug/y/GUvo274Jf7mkr12kFWBO7RRHh7gG1ritGp8WfhHZLCLLRWSJiGQ6z7UQkU9EZL1zb9+JTdA6VFbBuDeziIwUXrl6sHXcDEAREZ6TvN//sJdN+QfcjhMQ/LHHf5aqDlTVFGd5ArBAVbsBC5xlY4KOp+PmctbtKuKF9EEktrCOm4HqiuQEIiOEGZm21w/uHOoZCUxyHk8CLnEhgzEn7d//t5n3lm7n7uE9GNo93u04phatmsQyrGcrZmflUlZR5XYc1/m68CswX0SyROQm57nWqpoH4Ny3qumNInKTiGSKSGZ+fr6PYxpzfL7ftIfHP1zNOb2t42awSE9LZPeBMj5bs9PtKK7zdeEfoqrJwPnAbSIy1Ns3qupEVU1R1ZT4eNubMoFj5/4Sbpu6mA4tGvDMlQOIiLCTucHgzO6taNs0lml2kte3hV9Vtzv3u4B3gDRgp4i0BXDud/kygzF1qayiilsmZ3GorIJXrxlME+u4GTQiI4RRKYl8uT6fbQXFbsdxlc8Kv4g0FJHGhx8Dw4EVwFzgWudl1wJzfJXBmLr22AeryN5awFNX9Ke7ddwMOlemJADwVkZ47/X7co+/NfC1iCwFFgEfqOo84AngHBFZD5zjLBsT8N7OzuWNb7dw4xmduKh/O7fjmBOQ0LwBv+oWz8zMHCqrwnd2rihfrVhVNwEDanh+DzDMV9s1xhdWbi/kvreXc0qnFkw4v6fbccxJSE9N5NYp2Xy5Pp+zetQ4tiTkeb3HLyKxInKDiNwuInG+DGVMICk4VMa4yVk0bxDDS1clW8fNIHd2r9bENYwJ69m5juf/4OfxfEMoAd71TRxjAktVlXLnjCXsKCzh5bHJxDe2jpvBLiYqgssHJ/Dp6p3kF5W6HccVRy38IjJVRKoPUG4BTAGmAdZmwYSF5xasZ+HafB66uA/J1nEzZIxOTaSiSpmdnet2FFfUtsf/APAXEXlaRJoCT+MZkTMfeMQP2Yxx1YLVO3lhwXouT05g7CnWcTOUdIlvRFrHFszIyEE1/E7yHrXwq+omVb0Kz2GdGXjG4J+jqqer6ix/BTTGDZt3H+TOGUvo064Jf73UOm6GovS0RH7YfZDvf9jrdhS/q+1QT3MRuQ3oDVwJFAIfi8hF/gpnjBsOlVUwbnIWESK8OtY6boaq8/u2pXFsVFjOyVvboZ53gVIgFnhTVd8ALgYGi8hcf4Qzxt9UlfveXs7anUW8MMY6boay+jGRXDqoPR+u2EHhoXK34/hVbYU/DpiK54RuewBVLVbVR4Gb/ZDNGL/77zebmbNkO3ed3Z0zreNmyBudmkhZRRXvLA6vk7y1Ff6HgE/w9Nj5Wc/8w901jQklGZv38tcPVnN2r1bcdlZXt+MYP+jTrin9E5oyPcxO8tZ2cvdtVR2iqkNV9VN/hjLG33btL+HWKdkkNK/PM1cOtI6bYWR0aiJrdhSxNLfQ7Sh+Y5cgmrBXVlHFrVOyOVBSwWvXpNC0vnXcDCcjBrSjfnRkWJ3ktcJvwt7jH64mc8s+nrqiPz3aWMfNcNM4NpqL+rdl7tLtHCitcDuOX1jhN2HtncW5/PebzdxwRicuHmAdN8NVeloSh8oqeX/pdrej+MUxu3OKSDzwP0DH6q9X1et9F8sY31u1fT/3vb2cNOu4GfaSk5rRrVUjpmfkkJ4W+ldpe7PHPwdoCnwKfFDtZkzQKjxUzrjJWTStH80/rkom2jpuhjURIT0tiSU5BazZsd/tOD7nzf/tDVT1XlV9S1VnH775PJkxPuLpuLmYvMJiXr56sHXcNABcOqg9MZERTA+Dds3eFP73ReQCnycxxk9e+Gw9n6/N56GLejO4g3XcNB4tGsZwbt82vLN4GyXllW7H8SlvCv94PMW/WET2i0iRiIT+dyETkj5fs4vnF6znsuT2jD21g9txTIBJT02ksLicj1fucDuKTx2z8KtqY1WNUNX6qtrEWW7ij3DG1KUtew4yfvpierVpwuOX9rOOm+YXTuscR1KLBkwL8TH9tXXn7OncJ9d0819EY05ecVkl4yZnIyK8do113DQ1i4gQRqcm8t2mvfyw+6DbcXymtj3+u5z7Z2q4Pe3jXMbUGVXl/neWs2bHfp5PH2gdN02trhicQGSEMCMjdE/yHnUcv6re5Nyf5b84xtS9N77dwjuLt3HXOd35dY9WbscxAa51k1jO6tGKWVm5/GF495Ac6uvzTyQikSKyWETed5Y7icj3IrJeRGaISIyvM5jwlbl5L395fxVn92rF763jpvHSmLREdh8oZcHqXW5H8Ql//CkbD6yutvwk8HdV7QbsA27wQwYThnYVWcdNc2LO7B5PmyaxzMgIzZO8Pi38IpIAXAi87iwL8Bvg8Jy9k4BLfJnBhKfyyip+P2UxRSUVvHrNYOu4aY5LVGQEo1IS+GJdPtsLit2OU+eOWfjFY6yIPOQsJ4lImpfrfw64B6hyluOAAlU93AIvF2d2rxq2e5OIZIpIZn5+vpebM8bj8Q9Xs2jzXp64vB8929joY3P8rkxJRIGZmaE3O5c3e/wvA6cBY5zlIuAfx3qTMyn7LlXNqv50DS+tcdobVZ2oqimqmhIfb1PgGe/NWbKN//zfZn43pCMjB9a4X2HMMSW2aMAZXVvyVmYOlVWhNTuXN4X/FFW9DSgBUNV9gDcnZIcAI0RkMzAdzyGe54BmInJ4NFECEB59UI1frNmxnwmzl5PWsQX3X9DL7TgmyKWnJrGtoJivN+x2O0qd8qbwl4tIJM6eudOmuar2t4Cq3qeqCaraEUgHPlPVq4HPgSucl12Lp/unMSetsLicm9/MonFsFC9dPSgkh+EZ/zq7dytaNIwJudm5vPmX8QKeCddbichfga+Bx09im/cCd4nIBjzH/P91EusyBvB03LxrxhK27SvmlbHJtGoc63YkEwLqRUVyeXJ7Plm1k90HSt2OU2e86dUzBc8J2v8F8oBLVHXm8WxEVReq6kXO402qmqaqXVV1lKqGzm/TuObFzzawYM0uHryoN4M7tHA7jgkho1MTqahSZmeFzkleb0b1tAB2AdOAqcBOEbGxcSZgfL52F88tWMelg9rz29Os46apW11bNSa1Y3NmZOSgGhoneb051JMN5APrgPXO4x9EJFtEBvsynDG1qaxS/v31D9w2JZue1nHT+NDo1CQ27T7Ioh/2uh2lTnhT+OcBF6hqS1WNA84H3gJuxTPU0xi/W7V9P5e9/H/8+f1VpHVqwX+uS6V+jHXcNL5xYb+2NK4XFTKN27wp/Cmq+vHhBVWdDwxV1e8Am7PO+FVJeSVPzlvDiJe+JndfMc+nD+Q/16XSpqmdzDW+Uz8mkpGD2vHB8jwKD5W7HeekeVP494rIvSLSwbndA+xzhngec1inMXXlmw27Oe+5L3ll4UYuGdSeT+86k5ED29vhHeMX6alJlFZUMWfpNrejnDRvCv9VeC60ete5JTrPRQJX+i6aMR4Fh8r448ylXPX69ygw5cZTeHrUAJo3tMauxn/6tm9K3/ZNmLYo+E/yHrUfP3haKgP3qurtR3nJhrqPZIyHqvLesjz+/N5K9h0qZ9yZXbjz7G42e5ZxTXpqEg+8u4Ll2wrpn9DM7TgnrNY9flWtBGzkjvG7bQXF3DApkzumLaZds/rM/f0QJpzf04q+cdWIge2oHx3JtEXBfZK31j1+x2IRmQvMBH6chFJV3/ZZKhO2KquUSd9s5un5awF48KLeXHd6RyKtl74JAE1io7mwf1vmLtnGAxf2omE9b0po4PEmdQtgD54ma4cpYIXf1Kk1O/Zz7+zlLM0p4Nc94nnskr4kNLf5cU1gSU9NZFZWLh8sy+PK1ES345yQYxZ+Vf2dP4KY8FVSXsmLn63ntS820bR+NM+nD2TEgHY2WscEpMEdmtO1VSOmZ2wN3cIvIrF4pkfsA/w4WFpVr/dhLhMmvt24h/vfWc4Puw9yxeAE/nRBLxutYwKaiJCemshjH6xm3c4iurdu7Hak4+bNcM43gTbAucAXeIZ2FvkylAl9hYfKuXfWMsb88zuqVG2IpgkqlyUnEB0pTAvSds3eFP6uqvogcFBVJ+GZQ7efb2OZUKWqvL9sO8Oe/YJZ2bmMO7ML88YPZUjXlm5HM8ZrLRrGMLxPG95ZvI2S8kq34xw3ryZice4LRKQv0BTo6LNEJmRtLyjmxkmZ/H7qYto1i/1xiKb12DHBaExqEgWHyvl45Q63oxw3b0b1TBSR5sCDwFygEfCQT1OZkFJZpbz57Wb+9vFaqhQeuLAX153ekSibIcsEsdO7xJHYoj4zMnKCbm5nb0b1vO48/ALo7Ns4JtQcngN3SU4BZ3b3DNFMbGFDNE3wi4gQRqck8vT8dWzZc5AOcQ3djuQ1b0b11AMux3N458fXq+qffRfLBLvqQzSb2BBNE6KuGJzIs5+sY0ZGDvec19PtOF7z5lDPHKAQyAJsmkRzTNWHaF6enMADF9oQTROa2jSN5Tc9WzEzK5f/d053ooPk8KU3hT9BVc/zeRIT9AoPlfP4h6uZkZlDUosGvHlDGr/qFu92LGN8anRqEp+uzuTzNbsY3qeN23G84k3h/0ZE+qnqcp+nMUFJVflgeR6PzF3FvkNl3Dy0M3ee3d1G65iwcFaPeFo3qcf0jJzgL/wishxPT54o4HcisgnPoR4BVFX7+yeiCWTbC4p58N0VLFizi77tm/Df36XSt31Tt2MZ4zdRkRGMGpzIyws3kFdYTNum9d2OdEy17fFfdDIrdlo9fIlnesYoYJaqPiwinYDpeJq/ZQPXqGrZyWzL+N+RQzT/dEEvfjfEhmia8HRlSiIvfb6BmZm53DGsm9txjqm2f6WtgN6quqX6DegPeHOZZSnwG1UdAAwEzhORU4Engb+rajdgH54+QCaIrN1RxBWvfsMj760iuUNz5v+/ofzP0M5W9E3YSoprwBldWzIjI4eqqsCfnau2f6l/A1bX8Pxq52e1Uo8DzmK0c1M87Z1nOc9PAi7xOq1xVUl5JU9/vJYLX/iKLXsO8ffRA3jj+jQbl28MMDo1kW0FxXy9YbfbUY6ptkM9caq6+cgnVXWDiMR5s3Jn6sYsoCvwD2AjUKCqFc5LcoEaL3kTkZuAmwCSkpK82Zzxoe827eH+t5ezafdBLhvUngcu6k0LG6JpzI+G92lN8wbRzMjIYWj3wB7NVlvhr+0MhVeXqDlTNw4UkWbAO0Cvml52lPdOBCYCpKSkBP53pxBVWFzOEx+tZtqiHBJb1OeN69MC/n9qY9xQLyqSy5ITeOPbzew5UEpco3puRzqq2g71fCoif5UjLrUUkUeBz45nI6paACwETgWaicjhPzgJwPbjWZfxD1Xlw+V5nP3sF7yVmcvNQzvz8Z1DregbU4v01ETKK5W3s7e5HaVWtRX+P+DpzbNBRGY7tw1AD+CuY61YROKdPX1EpD5wNp7zA58DVzgvuxbPlcEmgOQVFvM/b2Rx65RsWjepx5zbhnDfBb1oEBOc84sa4y/dWjdmcIfmTMvYimrgHqg46r9kVT0IjBGRznhm3wJYqaqbvFx3W2CSc5w/AnhLVd8XkVXAdBF5DFgM/OvE45u6VFWlTP5+C0/NW0tFVZUN0TTmBKSnJvLHWcvI3LKP1I4t3I5TI2+6c24CvC321d+3DBh0lPWlHe/6jG+t21nEhNnLyN5awK+6teTxS/vZaB1jTsCF/dvy5/dWMW3R1oAt/LYrF+ZKyit5dr5niOYPuw/aEE1jTlKDmChGDGzHh8vzKCwuP/YbXHDUwu9cYWtC2KIf9nLBC1/xwmcbuLh/Oxb84ddcOijBWicbc5LSU5MoKa9i7pLAPMlb2x7/LAARWeCnLMZPCovLue/t5Vz52reUV1bxxvVpPDt6oI3LN6aO9EtoSp92TZiekeN2lBrVdow/QkQeBrqLyC9G8ajqs76LZXxBVZm3YgcPz13J7gOl3DS0M3ee3c1G6xjjA+mpiTw4ZyXLcwvplxBYjQtr2+NPB0rw/HFoXMPNBJG8wmJuejOLW6ZkE9+4HnNuO4P7bYimMT4zYmB7YqMjmJ6x1e0ov1DbcM61wJMiskxVP/JjJlOHjhyied/5PbnhjE42RNMYH2taP5oL+rVlzpLt/OnCwNrJ8uZf/zci8qyIZDq3Z0QksL63mBqt2+npovnQnJUMTGzG/DvP5OYzu1jRN8ZPxqQlcaC0gg+W5bkd5We8qQD/BoqAK53bfuA/vgxlTk5pxc+HaD4zagBv3pBGUpwN0TTGn1I6NKdLfMOAO8nrzXePLqp6ebXlR0Vkia8CmZOTtWUf98xaysb8g1wysB0PXtQ7oJtFGRPKRIT01CT++uFq1u8solvrwDg96s0ef7GInHF4QUSGAMW+i2RO1NodRYx9/XtKyquYdH0az6UPsqJvjMsuS25PdKQE1F6/N3v844A3qh3X34enuZoJIPtLyhk3OYtGsVG8fevptG4S63YkYwwQ16gew3u34e3sXO45rwf1oiLdjnTsPX5VXepMn9gf6K+qg5w+PCZAVFUpd81YSs7eQ/zjqmQr+sYEmNGpiew7VM78lTvdjgIcR68eVd2vqvt9GcacmJcXbuDT1Tu5/4JepHUKzKZQxoSzM7q2pH2z+swIkMM9Nq4vyH2xLp9nPlnHyIHt+N2Qjm7HMcbUICJCGJ2ayNcbdrN1zyG341jhD2Y5ew8xfvpierRuzP9e1s+aqxkTwEalJBAhMCPT/St5j1n4RaSBiDwoIv90lruJyEW+j2ZqU1JeybjJWVRWKa+OHRxQVwUaY36pbdP6/LpHK2Zm5lJRWeVqFm/2+P8DlAKnOcu5wGM+S2SOSVV54N0VrNy+nzkcj/AAABHGSURBVOdGD6Rjy4ZuRzLGeCE9NZFdRaV8vjbf1RzeFP4uqvoUUA6gqsWAHVNw0ZTvtzIrK5c7hnVjWK/WbscxxnjprJ6tiG9cjxkuN27zpvCXOZOlK4CIdMHzDcC4IHvrPh59byW/7hHPncO6uR3HGHMcoiMjGDU4gc/W7GJHYYlrObwp/I8A84BEEZkCLADu9WUoU7P8olJunZxNm6axPDd6IBER9sXLmGAzOjWRKoVZWe4N7fTmAq75wGXAdcA0IEVVP/dxLnOEisoqbp+Wzb5DZbw6djDNGthsWcYEow5xDTm9SxwzMnOoqlJXMngzqmeBqu5R1Q9U9X1V3W3TMfrfk/PW8N2mvfzvZf3o0866YhsTzNLTksjZW8w3G/e4sv3aJluPFZEWQEsRaS4iLZxbR6DdsVYsIoki8rmIrBaRlSIy3nm+hYh8IiLrnfvmdfVhQtX7y7bzz69+4LendeCy5AS34xhjTtLw3q1p1iCaaS6d5K1tj/9mIAvo6dwfvs0B/uHFuiuAP6hqL+BU4DYR6Q1MABaoajc85wsmnHj80LduZxH3zFrG4A7NeeDC3m7HMcbUgdjoSC4blMD8lTvYe7DM79s/auFX1edVtRNwt6p2VtVOzm2Aqr50rBWrap6qZjuPi4DVQHtgJDDJedkk4JKT/hQhan9JOePezKJBTBQvX51MTJRdaG1MqEhPS6S8Unk7O9fv2z7m5Z6q+qKI9AV6A7HVnn/D2404h4cGAd8DrVU1z1lHnoi0Osp7bgJuAkhKSvJ2UyGjqkq5+62lbNl7iKk3nmIdN40JMd1bNyY5qRnTM3K44YxOfm254s3J3YeBF53bWcBTwAhvNyAijYDZwJ3H091TVSeqaoqqpsTHx3v7tpDxyhcbmb/K03HzlM5xbscxxvhAemoSG3YdIGvLPr9u15tjB1cAw4Adqvo7YADg1bROIhKNp+hPUdW3nad3ikhb5+dtgV3HnTrEfbU+n2fmr+XiAe243jpuGhOyLuzflkb1ovw+O5dXUy+qahVQISJN8BTqzsd6k3i+t/wLWK2qz1b70Vx+msHrWjwni40jd98h7pi2mG6tGvPk5dZx05hQ1rBeFBcPaMf7y7azv6Tcb9v1pvBnikgz4J94RvVkA4u8eN8Q4BrgNyKyxLldADwBnCMi64FznGWDp+PmLZOzqahUXr3GOm4aEw7GpCVSUl7F3CXb/bZNb07u3uo8fFVE5gFNvJl6UVW/5ujN3IZ5HzE8qCoPzVnB8m2FvP7bFDpZx01jwkK/9k3p1bYJ0zO2MvbUDn7ZpldX7h5+rKqbVXWZXblb96YtyuGtzFxu/01Xzu5tHTeNCRciwpi0RFZs28+KbYV+2abPrtw13luSU8Ajc1cytHs8d57d3e04xhg/GzmgPfWiIpjupyt5fXnlrvHC7gOl3DI5i1ZN6vFC+kAireOmMWGnaYNoLuzXljmLt1NcVunz7fnsyl1zbBWVVdw+dTF7D1rHTWPC3ejURIpKK/hgeZ7Pt1XboZ5UEWmjqi86y78VkTki8oJzCMicpKc+Xsu3m/bw10v70be9ddw0JpyldWpB55YN/TI7V22Hel4DygBEZCieYZdvAIXARJ8nC3EfLMtj4pebGHtqElcMto6bxoQ7EWF0aiIZm/exYVeRT7dVW+GPVNW9zuPRwERVna2qDwJdfZoqxK3fWcQfZy1lUFIzHrqoj9txjDEB4vLBCURFCDN8fCVvrYVfRA6P8x8GfFbtZ3Zl0QkqKinn5jezaBATyStXD7aOm8aYH7VsVI9zerdmdvY2Sit8d5K3tqozDfhCROYAxcBXACLSFc/hHnOcVJW7Z3o6br50VTJtmlrHTWPMz6WnJbH3YBmfrvJdG7PaRvX8FfgD8F/gDFXVau+53WeJQtgrX2zk45U7ue/8npxqHTeNMTU4o2tL2jer79Mx/bUeZ1DV71T1HVU9WO25dYcnWDHe+3r9bp7+eC0X9W/LDWd0cjuOMSZARUYIV6Yk8tX63eTsPeSTbdgBZj/I3XeI26dl07VVI568vL913DTG1GpUSgIRAm9l+uYkrxV+H/tZx82xg2lYz86LG2Nq165Zfc7sHs/MzFwqKqvqfP1W+H3s4TkrWb6tkGeuHEDn+EZuxzHGBInRqZ6TvKvz6n5Mv+1++tC0RVuZkZnDbWd1YXifNm7HMcYEkWG9WvHd/cNo0bDuW7nYHr+PLMkp4OE5K/lVt5bcdU4Pt+MYY4JMdGSET4o+WOH3iT0HSrn1x46bg6zjpjEmoNihnjpWUVnF7dMWs+dgGbNvOZ3mPvqLbYwxJ8oKfx372/y1fLNxD3+7or913DTGBCQ71FOHPlqex2tfbOLqU5IYlZLodhxjjKmRFf46smFXEXfPXMrAxGY8dHFvt+MYY8xRWeGvAwdKK7j5zSzqx0Tyythk6kVFuh3JGGOOymeFX0T+LSK7RGRFtedaiMgnIrLeuW/uq+37i6ryx5lL2bznEC+OSaZt0/puRzLGmFr5co//v8B5Rzw3AVigqt2ABc5yUHvty018tGIHE87ryWldrOOmMSbw+azwq+qXwN4jnh4JTHIeTwIu8dX2/eH/NuzmqXlruLB/W278lXXcNMYEB38f42+tqnkAzn2ro71QRG4SkUwRyczPz/dbQG9tKyjm9mmL6RLfiKes46YxJogE7MldVZ2oqimqmhIfH+92nJ8pKa/k1slZlFVU8eo11nHTGBNc/F34d4pIWwDn3ndzi/nQo++tZGmup+NmF+u4aYwJMv4u/HOBa53H1wJz/Lz9kzYjYyvTFuVw66+7cK513DTGBCFfDuecBnwL9BCRXBG5AXgCOEdE1gPnOMtBY1luAQ86HTf/MNw6bhpjgpPPDk6r6pij/GiYr7bpS3sPlnHL5GziG9Xjeeu4aYwJYnZW0guVVcod0xaTf6CU2eNO91mPbGOM8YeAHdUTSJ6ev5avN+zmsZF96ZdgHTeNMcHNCv8xzFuxg1cWbmRMWhJXplrHTWNM8LPCX4sNuw5w98ylDEhsxiMjrOOmMSY0WOE/igOlFYybnEW9qAheudo6bhpjQoed3K2BqnLPrKVsyj/A5BtOoV0z67hpjAkdtsdfg39+tYkPl+/g3vN6cnrXlm7HMcaYOmWF/wjfbNzNEx+t4YJ+bbhpaGe34xhjTJ2zwl/N9oJibp+6mM7xjXjqigHWcdMYE5Ks8DtKKyq5ZUo2pRVVvDp2MI2s46YxJkRZdXM8+t4qluYU8OrYZLq2so6bxpjQZXv8wFsZOUz9fivjzuzCeX3buh3HGGN8KuwL/7LcAh6Ys4IhXeO4e3h3t+MYY4zPhXXhP9xxs2XDGF5IH0RUZFj/OowxYSJsj/H/2HGzqJSZ404jrlE9tyMZY4xfhG3hf8bpuPnEZf0YkNjM7TjGGOM3YXls4+OVO3h54UbGpCWSnpbkdhxjjPGrsCv8G/MP8Ie3ljIgoSmPjOjjdhxjjPG7sCr8B0srGPdmFjFREbwydrB13DTGhKWwKfyejpvL2Jh/gJfGDLKOm8aYsBU2hf/1r37gg+V53GMdN40xYS4sCv+3G/fwxLw1nN+3DTdbx01jTJhzpfCLyHkislZENojIBF9uK6+wmN9PzaZjXAP+Nso6bhpjjN8Lv4hEAv8Azgd6A2NExCcT2pZWVHLL5GxKyit57ZoU67hpjDG4s8efBmxQ1U2qWgZMB0b6YkN/fm8VS3IKeHrUAOu4aYwxDjcKf3sgp9pyrvPcz4jITSKSKSKZ+fn5x70RVaVTy4bcdlYXzu9nHTeNMeYwN4591HSQXX/xhOpEYCJASkrKL35+zI2IcOOv7ESuMcYcyY09/lwgsdpyArDdhRzGGBOW3Cj8GUA3EekkIjFAOjDXhRzGGBOW/H6oR1UrROT3wMdAJPBvVV3p7xzGGBOuXBnfqKofAh+6sW1jjAl3YXHlrjHGmJ9Y4TfGmDBjhd8YY8KMFX5jjAkzonrc10b5nYjkA1tO8O0tgd11GCcY2GcOD/aZQ9/Jft4Oqhp/5JNBUfhPhohkqmqK2zn8yT5zeLDPHPp89XntUI8xxoQZK/zGGBNmwqHwT3Q7gAvsM4cH+8yhzyefN+SP8RtjjPm5cNjjN8YYU40VfmOMCTMhXfj9Oal7IBCRf4vILhFZ4XYWfxCRRBH5XERWi8hKERnvdiZfE5FYEVkkIkudz/yo25n8RUQiRWSxiLzvdhZ/EJHNIrJcRJaISGadrjtUj/E7k7qvA87BM/lLBjBGVVe5GsyHRGQocAB4Q1X7up3H10SkLdBWVbNFpDGQBVwS4v+NBWioqgdEJBr4Ghivqt+5HM3nROQuIAVooqoXuZ3H10RkM5CiqnV+wVoo7/H7bVL3QKGqXwJ73c7hL6qap6rZzuMiYDU1zN8cStTjgLMY7dxCc++tGhFJAC4EXnc7SygI5cLv1aTuJjSISEdgEPC9u0l8zznksQTYBXyiqiH/mYHngHuAKreD+JEC80UkS0RuqssVh3Lh92pSdxP8RKQRMBu4U1X3u53H11S1UlUH4pmvOk1EQvqwnohcBOxS1Sy3s/jZEFVNBs4HbnMO5daJUC78Nql7GHCOc88Gpqjq227n8SdVLQAWAue5HMXXhgAjnGPe04HfiMhkdyP5nqpud+53Ae/gOXxdJ0K58Nuk7iHOOdH5L2C1qj7rdh5/EJF4EWnmPK4PnA2scTeVb6nqfaqaoKod8fw7/kxVx7ocy6dEpKEzYAERaQgMB+pstF7IFn5VrQAOT+q+Gngr1Cd1F5FpwLdADxHJFZEb3M7kY0OAa/DsAS5xbhe4HcrH2gKfi8gyPDs3n6hqWAxvDDOtga9FZCmwCPhAVefV1cpDdjinMcaYmoXsHr8xxpiaWeE3xpgwY4XfGGPCjBV+Y4wJM1b4jTEmzFjhN34nInHVhl/uEJFt1Za/8cH2UkTkheN8z19FJEdEDhzx/HUikl8t743Hud5xIvLb43h9hIi8ICIrnE6NGSLSyfnZ/cezbWMOs+GcxlUi8ghwQFWfdjtLdSJyKrAFWK+qjao9fx2ejom/91OOMcDlwJWqWuU0KzuoqvtE5ED1bMZ4y/b4TUA5vIctIr8WkS9E5C0RWSciT4jI1U4v+uUi0sV5XbyIzHb2hDNEZEgN6/z14R7uIvKIM2/BQhHZJCJ31JRDVb9T1bzjyO1t3kdE5G7n8UIRedJ5zToR+VUNq24L5KlqlZMr1yn6TwD1nW8dU5z1jXXWtUREXnNakyMiB0TkGRHJFpEFIhLvPH+HiKwSkWUiMt3bz2qCnxV+E8gGAOOBfniu0O2uqml4WvPe7rzmeeDvqpqKZ8/Ym7a9PYFz8fQ+edjp93M8LneK5SwRqd4Pypu8R4pyXnMn8HANP38LuNgp5s+IyCAAVZ0AFKvqQFW9WkR6AaPxNPYaCFQCVzvraAhkOw2/vqi2nQnAIFXtD4w7zt+BCWJW+E0gy3B67pcCG4H5zvPLgY7O47OBl5w2xXOBJod7nNTiA1UtdSa42IXn8nhvvQd0dIrlp8Ck48x7pMON5bJqeo2q5gI9gPvwtCReICLDaljPMGAwkOH8LoYBnZ2fVQEznMeTgTOcx8uAKSIyFqg4Sj4TgqLcDmBMLUqrPa6qtlzFT//vRgCnqWrxCa63kuP4d6Cqe6ot/hN48ijrPVreo2U5ag7nD8lHwEcishO4BFhwxMsEmKSq9x3rM/BTe/ILgaHACOBBEenj9LgyIc72+E2wm4+nGR8AIjLQlxsTz3SPh43A0wDQl9tLFpF2zuMIoD+ek84A5dUOUy0ArhCRVs5rW4hIB+dnEcAVzuOr8DT/igASVfVzPBOcNAPsRHGYsD1+E+zuAP7hdKuMAr6kDo5Xi8hTeIpkAxHJBV5X1UeAO0RkBJ5DI3uB6052W8fQCviniNRzlhcBLzmPJwLLRCTbOc7/AJ4ZmyKAcuA2PH8kDgJ9RCQLKMRzLiASmCwiTfF8W/i709/fhAEbzmlMiLNhn+ZIdqjHGGPCjO3xG2NMmLE9fmOMCTNW+I0xJsxY4TfGmDBjhd8YY8KMFX5jjAkz/x+AdHONxVYzLQAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "# plotting the state of charge\n",
    "ow.output['storage.soc_percent'].columns = ['Battery']\n",
    "ax = ow.output['storage.soc_percent'].plot()\n",
    "ax.set_xlabel('Time in 15min Steps')\n",
    "ax.set_ylabel('State of Charge in %')\n",
    "ax.legend()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "The state of charge grows until a certain point and drops in the end following the give values from the timeseries. The time windows between each timestep and the reactive power values define the amount of stored or consumed energy.\n",
    "\n",
    "In the shown case a SoC of more than 100\\% would be possible, because the SoC just gets summed up.\n",
    "To make the controller more realistic you could implement tresholds for the SoC and different initial SoC-values. Try to experiment!"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}