From 61bf3d157f21e82dbe7575eda71500859533729e Mon Sep 17 00:00:00 2001
From: Nicolas PIERRE <nicolas.pierre@artelys.com>
Date: Wed, 14 Dec 2022 16:21:00 +0100
Subject: [PATCH 1/2] dynawaltz example

Signed-off-by: Nicolas PIERRE <nicolas.pierre@artelys.com>
---
 dynawaltz.ipynb       | 287 ++++++++++++++++++++++++++++++++++++++++++
 dynawaltz/IEEE14.iidm | 224 +++++++++++++++++++++++++++++++++
 dynawaltz/config.yaml |  18 +++
 dynawaltz/models.par  |  66 ++++++++++
 dynawaltz/network.par |  25 ++++
 dynawaltz/solver.par  |  12 ++
 6 files changed, 632 insertions(+)
 create mode 100644 dynawaltz.ipynb
 create mode 100644 dynawaltz/IEEE14.iidm
 create mode 100644 dynawaltz/config.yaml
 create mode 100644 dynawaltz/models.par
 create mode 100644 dynawaltz/network.par
 create mode 100644 dynawaltz/solver.par

diff --git a/dynawaltz.ipynb b/dynawaltz.ipynb
new file mode 100644
index 0000000..e3a72e4
--- /dev/null
+++ b/dynawaltz.ipynb
@@ -0,0 +1,287 @@
+{
+ "cells": [
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "# Pypowsybl dynawaltz simulation"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 1,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import logging\n",
+    "import pypowsybl as pp\n",
+    "import pypowsybl.dynamic as dyn\n",
+    "import pandas as pd"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 2,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "network = pp.network.load(\"./dynawaltz/IEEE14.iidm\")"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "You can add dynamic mappings with dataframes like this :"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 3,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "df_alpha_beta = pd.DataFrame.from_dict({\"static_id\": [network.get_loads().loc[l].name for l in network.get_loads().index],\n",
+    "                             \"parameter_set_id\": [\"LAB\" for l in network.get_loads().index]})\n",
+    "df_GSTWPR = pd.DataFrame.from_dict({\"static_id\": [network.get_generators().loc[l].name for l in network.get_generators().index],\n",
+    "                             \"parameter_set_id\": [\"GSTWPR\" for l in network.get_generators().index]})\n",
+    "df_omega_ref = pd.DataFrame.from_dict({\"generator_id\": [network.get_generators().loc[l].name for l in network.get_generators().index]})"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Don't forget to index the dataframe on the network element id column"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 4,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "model_mapping = dyn.ModelMapping()\n",
+    "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.ALPHA_BETA_LOAD, df_alpha_beta.set_index(\"static_id\"))\n",
+    "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.GENERATOR_SYNCHRONOUS_THREE_WINDINGS_PROPORTIONAL_REGULATIONS, df_GSTWPR.set_index(\"static_id\"))\n",
+    "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.OMEGA_REF, df_omega_ref.set_index(\"generator_id\"))"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Adding events"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 5,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "events = dyn.EventMapping()\n",
+    "events.add_event(\"EQD\", dyn.EventType.BRANCH_DISCONNECTION, \"_BUS____1-BUS____5-1_AC\")"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Adding curves, you can batch curves creation for a given id"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 6,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "timeseries = dyn.CurveMapping()\n",
+    "timeseries.add_curves(\"_LOAD___2_EC\", [\"load_PPu\", \"load_QPu\"])"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Running the simulation (will load the default config file in ~/.itools folder)"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "sim = dyn.Simulation()\n",
+    "res = sim.run(network, model_mapping, events, timeseries, 0, 30)"
+   ]
+  },
+  {
+   "attachments": {},
+   "cell_type": "markdown",
+   "metadata": {},
+   "source": [
+    "Display data of the run:"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 9,
+   "metadata": {},
+   "outputs": [
+    {
+     "data": {
+      "text/html": [
+       "<div>\n",
+       "<style scoped>\n",
+       "    .dataframe tbody tr th:only-of-type {\n",
+       "        vertical-align: middle;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe tbody tr th {\n",
+       "        vertical-align: top;\n",
+       "    }\n",
+       "\n",
+       "    .dataframe thead th {\n",
+       "        text-align: right;\n",
+       "    }\n",
+       "</style>\n",
+       "<table border=\"1\" class=\"dataframe\">\n",
+       "  <thead>\n",
+       "    <tr style=\"text-align: right;\">\n",
+       "      <th></th>\n",
+       "      <th>_LOAD___2_EC_load_QPu</th>\n",
+       "      <th>_LOAD___2_EC_load_PPu</th>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>timestamp</th>\n",
+       "      <th></th>\n",
+       "      <th></th>\n",
+       "    </tr>\n",
+       "  </thead>\n",
+       "  <tbody>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>0.126992</td>\n",
+       "      <td>0.216992</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>0.126992</td>\n",
+       "      <td>0.216992</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>0.126992</td>\n",
+       "      <td>0.216992</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>0.126992</td>\n",
+       "      <td>0.216992</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>0</th>\n",
+       "      <td>0.126992</td>\n",
+       "      <td>0.216992</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>...</th>\n",
+       "      <td>...</td>\n",
+       "      <td>...</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>9900</th>\n",
+       "      <td>0.125207</td>\n",
+       "      <td>0.215157</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>12710</th>\n",
+       "      <td>0.125208</td>\n",
+       "      <td>0.215157</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>18331</th>\n",
+       "      <td>0.125208</td>\n",
+       "      <td>0.215158</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>28331</th>\n",
+       "      <td>0.125208</td>\n",
+       "      <td>0.215158</td>\n",
+       "    </tr>\n",
+       "    <tr>\n",
+       "      <th>30000</th>\n",
+       "      <td>0.125208</td>\n",
+       "      <td>0.215158</td>\n",
+       "    </tr>\n",
+       "  </tbody>\n",
+       "</table>\n",
+       "<p>252 rows × 2 columns</p>\n",
+       "</div>"
+      ],
+      "text/plain": [
+       "           _LOAD___2_EC_load_QPu  _LOAD___2_EC_load_PPu\n",
+       "timestamp                                              \n",
+       "0                       0.126992               0.216992\n",
+       "0                       0.126992               0.216992\n",
+       "0                       0.126992               0.216992\n",
+       "0                       0.126992               0.216992\n",
+       "0                       0.126992               0.216992\n",
+       "...                          ...                    ...\n",
+       "9900                    0.125207               0.215157\n",
+       "12710                   0.125208               0.215157\n",
+       "18331                   0.125208               0.215158\n",
+       "28331                   0.125208               0.215158\n",
+       "30000                   0.125208               0.215158\n",
+       "\n",
+       "[252 rows x 2 columns]"
+      ]
+     },
+     "execution_count": 9,
+     "metadata": {},
+     "output_type": "execute_result"
+    }
+   ],
+   "source": [
+    "res.status()\n",
+    "res.curves()"
+   ]
+  }
+ ],
+ "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.8.10 (default, Jun 22 2022, 20:18:18) \n[GCC 9.4.0]"
+  },
+  "orig_nbformat": 4,
+  "vscode": {
+   "interpreter": {
+    "hash": "e7370f93d1d0cde622a1f8e1c04877d8463912d04d973331ad4851f04de6915a"
+   }
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 2
+}
diff --git a/dynawaltz/IEEE14.iidm b/dynawaltz/IEEE14.iidm
new file mode 100644
index 0000000..df38208
--- /dev/null
+++ b/dynawaltz/IEEE14.iidm
@@ -0,0 +1,224 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<iidm:network xmlns:iidm="http://www.itesla_project.eu/schema/iidm/1_0" id="ieee14bus" caseDate="2017-06-09T10:14:24.146+02:00" forecastDistance="0" sourceFormat="CIM1">
+    <iidm:substation id="_BUS___10_SS" name="BUS   10_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS___10_VL" name="BUS   10_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS___10_TN" v="14.5036" angle="-15.0972"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD__10_EC" name="LOAD  10" loadType="UNDEFINED" p0="9.0" q0="5.8" bus="_BUS___10_TN" connectableBus="_BUS___10_TN" p="9.0" q="5.8"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS___11_SS" name="BUS   11_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS___11_VL" name="BUS   11_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS___11_TN" v="14.5853" angle="-14.7906"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD__11_EC" name="LOAD  11" loadType="UNDEFINED" p0="3.5" q0="1.8" bus="_BUS___11_TN" connectableBus="_BUS___11_TN" p="3.5" q="1.8"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS___12_SS" name="BUS   12_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS___12_VL" name="BUS   12_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS___12_TN" v="14.5616" angle="-15.0755"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD__12_EC" name="LOAD  12" loadType="UNDEFINED" p0="6.1" q0="1.6" bus="_BUS___12_TN" connectableBus="_BUS___12_TN" p="6.1" q="1.6"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS___13_SS" name="BUS   13_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS___13_VL" name="BUS   13_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS___13_TN" v="14.4952" angle="-15.15652"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD__13_EC" name="LOAD  13" loadType="UNDEFINED" p0="13.5" q0="5.8" bus="_BUS___13_TN" connectableBus="_BUS___13_TN" p="13.5" q="5.8"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS___14_SS" name="BUS   14_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS___14_VL" name="BUS   14_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS___14_TN" v="14.306159" angle="-16.0336"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD__14_EC" name="LOAD  14" loadType="UNDEFINED" p0="14.9" q0="5.0" bus="_BUS___14_TN" connectableBus="_BUS___14_TN" p="14.9" q="5.0"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS____1_SS" name="BUS    1_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS____1_VL" name="BUS    1_VL" nominalV="69.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____1_TN" v="73.14" angle="0.0"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="_GEN____1_SM" name="GEN    1" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="232.3463" targetV="73.14" targetQ="-16.759" bus="_BUS____1_TN" connectableBus="_BUS____1_TN" p="-232.39" q="16.55">
+                <iidm:minMaxReactiveLimits minQ="-999.0" maxQ="999.0"/>
+            </iidm:generator>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS____2_SS" name="BUS    2_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS____2_VL" name="BUS    2_VL" nominalV="69.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____2_TN" v="72.11" angle="-4.98"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="_GEN____2_SM" name="GEN    2" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="40.0" targetV="72.105" targetQ="42.4" bus="_BUS____2_TN" connectableBus="_BUS____2_TN" p="-40.0" q="-43.56">
+                <iidm:minMaxReactiveLimits minQ="-40.0" maxQ="50.0"/>
+            </iidm:generator>
+            <iidm:load id="_LOAD___2_EC" name="LOAD   2" loadType="UNDEFINED" p0="21.7" q0="12.7" bus="_BUS____2_TN" connectableBus="_BUS____2_TN" p="21.7" q="12.7"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS____3_SS" name="BUS    3_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS____3_VL" name="BUS    3_VL" nominalV="69.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____3_TN" v="69.69" angle="-12.73"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="_GEN____3_SM" name="GEN    3" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="0.0" targetV="69.69" targetQ="23.4" bus="_BUS____3_TN" connectableBus="_BUS____3_TN" p="-0.0" q="-25.07">
+                <iidm:minMaxReactiveLimits minQ="0.0" maxQ="40.0"/>
+            </iidm:generator>
+            <iidm:load id="_LOAD___3_EC" name="LOAD   3" loadType="UNDEFINED" p0="94.2" q0="19.0" bus="_BUS____3_TN" connectableBus="_BUS____3_TN" p="94.2" q="19.0"/>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:substation id="_BUS____4_SS" name="BUS    4_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS____9_VL" name="BUS    9_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____9_TN" v="14.5719" angle="-14.9385"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD___9_EC" name="LOAD   9" loadType="UNDEFINED" p0="29.5" q0="16.6" bus="_BUS____9_TN" connectableBus="_BUS____9_TN" p="29.5" q="16.6"/>
+            <iidm:shunt id="_BANK___9_SC" name="BANK   9" bPerSection="0.099769" maximumSectionCount="1" currentSectionCount="1" bus="_BUS____9_TN" connectableBus="_BUS____9_TN" q="-21.256718"/>
+        </iidm:voltageLevel>
+        <iidm:voltageLevel id="_BUS____7_VL" name="BUS    7_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____7_TN" v="14.649" angle="-13.3596"/>
+            </iidm:busBreakerTopology>
+        </iidm:voltageLevel>
+        <iidm:voltageLevel id="_BUS____4_VL" name="BUS    4_VL" nominalV="69.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____4_TN" v="70.2193" angle="-10.3129"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD___4_EC" name="LOAD   4" loadType="UNDEFINED" p0="47.8" q0="-3.9" bus="_BUS____4_TN" connectableBus="_BUS____4_TN" p="47.8" q="-3.9"/>
+        </iidm:voltageLevel>
+        <iidm:twoWindingsTransformer id="_BUS____4-BUS____9-1_PT" name="BUS    4-BUS    9-1" r="0.0" x="1.0591881" g="0.0" b="0.0" ratedU1="69.0" ratedU2="13.8" bus1="_BUS____4_TN" connectableBus1="_BUS____4_TN" voltageLevelId1="_BUS____4_VL" bus2="_BUS____9_TN" connectableBus2="_BUS____9_TN" voltageLevelId2="_BUS____9_VL" p1="16.299362" q1="1.270369" p2="-16.299362" q2="0.050373">
+            <iidm:ratioTapChanger lowTapPosition="1" tapPosition="4" regulating="false" loadTapChangingCapabilities="false">
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.1111112"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0834236"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0570825"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0319917"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0157440"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9803922"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9615385"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9433963"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9259259"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9090909"/>
+            </iidm:ratioTapChanger>
+            <iidm:currentLimits1 permanentLimit="836.74"/>
+            <iidm:currentLimits2 permanentLimit="4183.7"/>
+        </iidm:twoWindingsTransformer>
+        <iidm:twoWindingsTransformer id="_BUS____4-BUS____7-1_PT" name="BUS    4-BUS    7-1" r="0.0" x="0.39824802" g="0.0" b="0.0" ratedU1="69.0" ratedU2="13.8" bus1="_BUS____4_TN" connectableBus1="_BUS____4_TN" voltageLevelId1="_BUS____4_VL" bus2="_BUS____7_TN" connectableBus2="_BUS____7_TN" voltageLevelId2="_BUS____7_VL" p1="28.129929" q1="-10.561864" p2="-28.129929" q2="12.3099">
+            <iidm:ratioTapChanger lowTapPosition="1" tapPosition="5" regulating="false" loadTapChangingCapabilities="false">
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.1111112"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0875476"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0649627"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0432966"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0224948"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.98039216"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9615385"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9433963"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9259259"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9090909"/>
+            </iidm:ratioTapChanger>
+            <iidm:currentLimits1 permanentLimit="836.74"/>
+            <iidm:currentLimits2 permanentLimit="4183.7"/>
+        </iidm:twoWindingsTransformer>
+    </iidm:substation>
+    <iidm:substation id="_BUS____5_SS" name="BUS    5_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS____6_VL" name="BUS    6_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____6_TN" v="14.77" angle="-14.22"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="_GEN____6_SM" name="GEN    6" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="0.0" targetV="14.766" targetQ="12.2" bus="_BUS____6_TN" connectableBus="_BUS____6_TN" p="-0.0" q="-12.73">
+                <iidm:minMaxReactiveLimits minQ="-6.0" maxQ="24.0"/>
+            </iidm:generator>
+            <iidm:load id="_LOAD___6_EC" name="LOAD   6" loadType="UNDEFINED" p0="11.2" q0="7.5" bus="_BUS____6_TN" connectableBus="_BUS____6_TN" p="11.2" q="7.5"/>
+        </iidm:voltageLevel>
+        <iidm:voltageLevel id="_BUS____5_VL" name="BUS    5_VL" nominalV="69.0" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____5_TN" v="70.3464" angle="-8.77381"/>
+            </iidm:busBreakerTopology>
+            <iidm:load id="_LOAD___5_EC" name="LOAD   5" loadType="UNDEFINED" p0="7.6" q0="1.6" bus="_BUS____5_TN" connectableBus="_BUS____5_TN" p="7.6" q="1.6"/>
+        </iidm:voltageLevel>
+        <iidm:twoWindingsTransformer id="_BUS____5-BUS____6-1_PT" name="BUS    5-BUS    6-1" r="0.0" x="0.47994804" g="0.0" b="0.0" ratedU1="69.0" ratedU2="13.8" bus1="_BUS____5_TN" connectableBus1="_BUS____5_TN" voltageLevelId1="_BUS____5_VL" bus2="_BUS____6_TN" connectableBus2="_BUS____6_TN" voltageLevelId2="_BUS____6_VL" p1="43.804256" q1="9.096129" p2="-43.804256" q2="-4.821185">
+            <iidm:ratioTapChanger lowTapPosition="1" tapPosition="3" regulating="false" loadTapChangingCapabilities="false">
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.1111112"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0917031"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0729614"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0474860"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0231924"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="1.0"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9803922"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9615385"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9433963"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9259259"/>
+                <iidm:step r="0.0" x="0.0" g="0.0" b="0.0" rho="0.9090909"/>
+            </iidm:ratioTapChanger>
+            <iidm:currentLimits1 permanentLimit="836.74"/>
+            <iidm:currentLimits2 permanentLimit="4183.7"/>
+        </iidm:twoWindingsTransformer>
+    </iidm:substation>
+    <iidm:substation id="_BUS____8_SS" name="BUS    8_SS" country="AF" geographicalTags="_SGR_01">
+        <iidm:voltageLevel id="_BUS____8_VL" name="BUS    8_VL" nominalV="13.8" topologyKind="BUS_BREAKER">
+            <iidm:busBreakerTopology>
+                <iidm:bus id="_BUS____8_TN" v="15.04" angle="-13.36"/>
+            </iidm:busBreakerTopology>
+            <iidm:generator id="_GEN____8_SM" name="GEN    8" energySource="OTHER" minP="-9999.0" maxP="9999.0" voltageRegulatorOn="true" targetP="0.0" targetV="15.042" targetQ="17.4" bus="_BUS____8_TN" connectableBus="_BUS____8_TN" p="-0.0" q="-17.62">
+                <iidm:minMaxReactiveLimits minQ="-6.0" maxQ="24.0"/>
+            </iidm:generator>
+        </iidm:voltageLevel>
+    </iidm:substation>
+    <iidm:line id="_BUS___10-BUS___11-1_AC" name="BUS   10-BUS   11-1" r="0.156256" x="0.365778" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS___10_TN" connectableBus1="_BUS___10_TN" voltageLevelId1="_BUS___10_VL" bus2="_BUS___11_TN" connectableBus2="_BUS___11_TN" voltageLevelId2="_BUS___11_VL" p1="-3.628976" q1="-1.291152" p2="3.639966" q2="1.316878">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS___12-BUS___13-1_AC" name="BUS   12-BUS   13-1" r="0.42072" x="0.380651" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS___12_TN" connectableBus1="_BUS___12_TN" voltageLevelId1="_BUS___12_VL" bus2="_BUS___13_TN" connectableBus2="_BUS___13_TN" voltageLevelId2="_BUS___13_VL" p1="1.584024" q1="0.715153" p2="-1.578033" q2="-0.709732">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS___13-BUS___14-1_AC" name="BUS   13-BUS   14-1" r="0.325519" x="0.662769" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS___13_TN" connectableBus1="_BUS___13_TN" voltageLevelId1="_BUS___13_VL" bus2="_BUS___14_TN" connectableBus2="_BUS___14_TN" voltageLevelId2="_BUS___14_VL" p1="5.526893" q1="1.540022" p2="-5.47592" q2="-1.43624">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____1-BUS____2-1_AC" name="BUS    1-BUS    2-1" r="0.922682" x="2.81708" g1="0.0" b1="5.54505E-4" g2="0.0" b2="5.54505E-4" bus1="_BUS____1_TN" connectableBus1="_BUS____1_TN" voltageLevelId1="_BUS____1_VL" bus2="_BUS____2_TN" connectableBus2="_BUS____2_TN" voltageLevelId2="_BUS____2_VL" p1="156.78983" q1="-20.382833" p2="-152.49738" q2="27.639011">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____1-BUS____5-1_AC" name="BUS    1-BUS    5-1" r="2.57237" x="10.6189" g1="0.0" b1="5.167E-4" g2="0.0" b2="5.167E-4" bus1="_BUS____1_TN" connectableBus1="_BUS____1_TN" voltageLevelId1="_BUS____1_VL" bus2="_BUS____5_TN" connectableBus2="_BUS____5_TN" voltageLevelId2="_BUS____5_VL" p1="75.579735" q1="3.118322" p2="-72.81625" q2="2.96058">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____2-BUS____3-1_AC" name="BUS    2-BUS    3-1" r="2.23719" x="9.42535" g1="0.0" b1="4.599875E-4" g2="0.0" b2="4.599875E-4" bus1="_BUS____2_TN" connectableBus1="_BUS____2_TN" voltageLevelId1="_BUS____2_VL" bus2="_BUS____3_TN" connectableBus2="_BUS____3_TN" voltageLevelId2="_BUS____3_VL" p1="73.19019" q1="3.564935" p2="-70.86989" q2="1.58502">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____2-BUS____4-1_AC" name="BUS    2-BUS    4-1" r="2.76662" x="8.3946" g1="0.0" b1="3.57068E-4" g2="0.0" b2="3.57068E-4" bus1="_BUS____2_TN" connectableBus1="_BUS____2_TN" voltageLevelId1="_BUS____2_VL" bus2="_BUS____4_TN" connectableBus2="_BUS____4_TN" voltageLevelId2="_BUS____4_VL" p1="56.126595" q1="-2.020396" p2="-54.450264" q2="3.486913">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____2-BUS____5-1_AC" name="BUS    2-BUS    5-1" r="2.71139" x="8.27843" g1="0.0" b1="3.63369E-4" g2="0.0" b2="3.63369E-4" bus1="_BUS____2_TN" connectableBus1="_BUS____2_TN" voltageLevelId1="_BUS____2_VL" bus2="_BUS____5_TN" connectableBus2="_BUS____5_TN" voltageLevelId2="_BUS____5_VL" p1="41.48059" q1="0.250869" p2="-40.580875" q2="-1.196797">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____3-BUS____4-1_AC" name="BUS    3-BUS    4-1" r="3.19035" x="8.14274" g1="0.0" b1="1.344255E-4" g2="0.0" b2="1.344255E-4" bus1="_BUS____3_TN" connectableBus1="_BUS____3_TN" voltageLevelId1="_BUS____3_VL" bus2="_BUS____4_TN" connectableBus2="_BUS____4_TN" voltageLevelId2="_BUS____4_VL" p1="-23.330109" q1="4.002362" p2="23.701889" q2="-4.37021">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____4-BUS____5-1_AC" name="BUS    4-BUS    5-1" r="0.635593" x="2.00486" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____4_TN" connectableBus1="_BUS____4_TN" voltageLevelId1="_BUS____4_VL" bus2="_BUS____5_TN" connectableBus2="_BUS____5_TN" voltageLevelId2="_BUS____5_VL" p1="-61.48091" q1="14.074792" p2="61.99287" q2="-12.459913">
+        <iidm:currentLimits1 permanentLimit="836.74"/>
+    </iidm:line>
+    <iidm:line id="_BUS____6-BUS___11-1_AC" name="BUS    6-BUS   11-1" r="0.18088" x="0.378785" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____6_TN" connectableBus1="_BUS____6_TN" voltageLevelId1="_BUS____6_VL" bus2="_BUS___11_TN" connectableBus2="_BUS___11_TN" voltageLevelId2="_BUS___11_VL" p1="7.191498" q1="3.224791" p2="-7.139966" q2="-3.116878">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____6-BUS___12-1_AC" name="BUS    6-BUS   12-1" r="0.23407" x="0.487165" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____6_TN" connectableBus1="_BUS____6_TN" voltageLevelId1="_BUS____6_VL" bus2="_BUS___12_TN" connectableBus2="_BUS___12_TN" voltageLevelId2="_BUS___12_VL" p1="7.755102" q1="2.463086" p2="-7.684024" q2="-2.315153">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____6-BUS___13-1_AC" name="BUS    6-BUS   13-1" r="0.125976" x="0.248086" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____6_TN" connectableBus1="_BUS____6_TN" voltageLevelId1="_BUS____6_VL" bus2="_BUS___13_TN" connectableBus2="_BUS___13_TN" voltageLevelId2="_BUS___13_VL" p1="17.657656" q1="7.041473" p2="-17.44886" q2="-6.630291">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____7-BUS____8-1_AC" name="BUS    7-BUS    8-1" r="0.0" x="0.33546" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____7_TN" connectableBus1="_BUS____7_TN" voltageLevelId1="_BUS____7_VL" bus2="_BUS____8_TN" connectableBus2="_BUS____8_TN" voltageLevelId2="_BUS____8_VL" p1="-0.0" q1="-16.861053" p2="0.0" q2="17.305046">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____7-BUS____9-1_AC" name="BUS    7-BUS    9-1" r="0.0" x="0.209503" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____7_TN" connectableBus1="_BUS____7_TN" voltageLevelId1="_BUS____7_VL" bus2="_BUS____9_TN" connectableBus2="_BUS____9_TN" voltageLevelId2="_BUS____9_VL" p1="28.129929" q1="4.551154" p2="-28.129929" q2="-3.759173">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____9-BUS___10-1_AC" name="BUS    9-BUS   10-1" r="0.060579" x="0.160922" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____9_TN" connectableBus1="_BUS____9_TN" voltageLevelId1="_BUS____9_VL" bus2="_BUS___10_TN" connectableBus2="_BUS___10_TN" voltageLevelId2="_BUS___10_VL" p1="5.385146" q1="4.546363" p2="-5.371024" q2="-4.508848">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+    <iidm:line id="_BUS____9-BUS___14-1_AC" name="BUS    9-BUS   14-1" r="0.242068" x="0.514912" g1="0.0" b1="0.0" g2="0.0" b2="0.0" bus1="_BUS____9_TN" connectableBus1="_BUS____9_TN" voltageLevelId1="_BUS____9_VL" bus2="_BUS___14_TN" connectableBus2="_BUS___14_TN" voltageLevelId2="_BUS___14_VL" p1="9.544144" q1="3.819154" p2="-9.42408" q2="-3.56376">
+        <iidm:currentLimits1 permanentLimit="4183.7"/>
+    </iidm:line>
+</iidm:network>
diff --git a/dynawaltz/config.yaml b/dynawaltz/config.yaml
new file mode 100644
index 0000000..7b5176b
--- /dev/null
+++ b/dynawaltz/config.yaml
@@ -0,0 +1,18 @@
+# copy to ~/.itools/config.yaml
+dynamic-simulation-default-parameters:
+  startTime: 0
+  stopTime: 30
+
+dynawaltz:
+  homeDir: PATH/TO/DYNAWO # root directory of dynawo installation
+  debug: true
+
+dynawaltz-default-parameters:
+  # All the pathes can be relative to the root directory
+  parametersFile: ./dynawaltz/models.par
+  network.parametersFile: ./dynawaltz/network.par
+  network.parametersId: "1"
+  solver.type: IDA
+  solver.parametersFile: ./dynawaltz/solver.par
+  solver.parametersId: "1"
+
diff --git a/dynawaltz/models.par b/dynawaltz/models.par
new file mode 100644
index 0000000..6054bdd
--- /dev/null
+++ b/dynawaltz/models.par
@@ -0,0 +1,66 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<parametersSet xmlns="http://www.rte-france.com/dynawo">
+     <set id="LAB">
+        <par type="DOUBLE" name="load_alpha" value="1.5"/>
+        <par type="DOUBLE" name="load_beta" value="2.5"/>
+        <reference type="DOUBLE" name="load_P0Pu" origData="IIDM" origName="p_pu"/>
+        <reference type="DOUBLE" name="load_Q0Pu" origData="IIDM" origName="q_pu"/>
+        <reference type="DOUBLE" name="load_U0Pu" origData="IIDM" origName="v_pu"/>
+        <reference type="DOUBLE" name="load_UPhase0" origData="IIDM" origName="angle_pu"/>
+    </set>
+    <set id="GSTWPR">
+        <par type="INT" name="generator_ExcitationPu" value="1" />
+        <par type="DOUBLE" name="generator_md" value="0.16" />
+        <par type="DOUBLE" name="generator_mq" value="0.16" />
+        <par type="DOUBLE" name="generator_nd" value="5.7" />
+        <par type="DOUBLE" name="generator_nq" value="5.7" />
+        <par type="DOUBLE" name="generator_MdPuEfd" value="0" />
+        <par type="DOUBLE" name="generator_DPu" value="0" />
+        <par type="DOUBLE" name="generator_H" value="4.9750000000000001" />
+        <par type="DOUBLE" name="generator_RaPu" value="0.004" />
+        <par type="DOUBLE" name="generator_XlPu" value="0.102" />
+        <par type="DOUBLE" name="generator_XdPu" value="0.75" />
+        <par type="DOUBLE" name="generator_XpdPu" value="0.225" />
+        <par type="DOUBLE" name="generator_XppdPu" value="0.154" />
+        <par type="DOUBLE" name="generator_Tpd0" value="3" />
+        <par type="DOUBLE" name="generator_Tppd0" value="0.04" />
+        <par type="DOUBLE" name="generator_XqPu" value="0.45" />
+        <par type="DOUBLE" name="generator_XppqPu" value="0.2" />
+        <par type="DOUBLE" name="generator_Tppq0" value="0.04" />
+        <par type="DOUBLE" name="generator_UNom" value="15" />
+        <par type="DOUBLE" name="generator_PNomTurb" value="74.4" />
+        <par type="DOUBLE" name="generator_PNomAlt" value="74.4" />
+        <par type="DOUBLE" name="generator_SNom" value="80" />
+        <par type="DOUBLE" name="generator_SnTfo" value="80" />
+        <par type="DOUBLE" name="generator_UNomHV" value="15" />
+        <par type="DOUBLE" name="generator_UNomLV" value="15" />
+        <par type="DOUBLE" name="generator_UBaseHV" value="15" />
+        <par type="DOUBLE" name="generator_UBaseLV" value="15" />
+        <par type="DOUBLE" name="generator_RTfPu" value="0.0" />
+        <par type="DOUBLE" name="generator_XTfPu" value="0.0" />
+        <par type="DOUBLE" name="PmPu_ValueIn" value="0" />
+        <par type="DOUBLE" name="EfdPu_ValueIn" value="0" />
+        <par type="DOUBLE" name="voltageRegulator_LagEfdMax" value="0" />
+        <par type="DOUBLE" name="voltageRegulator_LagEfdMin" value="0" />
+        <par type="DOUBLE" name="voltageRegulator_EfdMinPu" value="-5" />
+        <par type="DOUBLE" name="voltageRegulator_EfdMaxPu" value="5" />
+        <par type="DOUBLE" name="voltageRegulator_UsRefMinPu" value="0.8" />
+        <par type="DOUBLE" name="voltageRegulator_UsRefMaxPu" value="1.2" />
+        <par type="DOUBLE" name="voltageRegulator_Gain" value="20" />
+        <par type="DOUBLE" name="governor_KGover" value="5" />
+        <par type="DOUBLE" name="governor_PMin" value="0" />
+        <par type="DOUBLE" name="governor_PMax" value="74.4" />
+        <par type="DOUBLE" name="governor_PNom" value="74.4" />
+        <par type="DOUBLE" name="URef_ValueIn" value="0" />
+        <par type="DOUBLE" name="Pm_ValueIn" value="0" />
+        <reference name="generator_P0Pu" origData="IIDM" origName="p_pu" type="DOUBLE" />
+        <reference name="generator_Q0Pu" origData="IIDM" origName="q_pu" type="DOUBLE" />
+        <reference name="generator_U0Pu" origData="IIDM" origName="v_pu" type="DOUBLE" />
+        <reference name="generator_UPhase0" origData="IIDM" origName="angle_pu" type="DOUBLE" />
+    </set>
+    <set id="EQD">
+        <par type="DOUBLE" name="event_tEvent" value="1"/>
+        <par type="BOOL" name="event_disconnectOrigin" value="false"/>
+        <par type="BOOL" name="event_disconnectExtremity" value="true"/>
+    </set>
+</parametersSet>
diff --git a/dynawaltz/network.par b/dynawaltz/network.par
new file mode 100644
index 0000000..91c6cda
--- /dev/null
+++ b/dynawaltz/network.par
@@ -0,0 +1,25 @@
+<?xml version="1.0" encoding="UTF-8"?>
+<parametersSet xmlns="http://www.rte-france.com/dynawo">
+    <set id="1">
+        <par type="DOUBLE" name="capacitor_no_reclosing_delay" value="300" />
+        <par type="DOUBLE" name="dangling_line_currentLimit_maxTimeOperation" value="90" />
+        <par type="DOUBLE" name="line_currentLimit_maxTimeOperation" value="90" />
+        <par type="DOUBLE" name="load_Tp" value="90" />
+        <par type="DOUBLE" name="load_Tq" value="90" />
+        <par type="DOUBLE" name="load_alpha" value="1" />
+        <par type="DOUBLE" name="load_alphaLong" value="0" />
+        <par type="DOUBLE" name="load_beta" value="2" />
+        <par type="DOUBLE" name="load_betaLong" value="0" />
+        <par type="BOOL" name="load_isControllable" value="false" />
+        <par type="BOOL" name="load_isRestorative" value="false" />
+        <par type="DOUBLE" name="load_zPMax" value="100" />
+        <par type="DOUBLE" name="load_zQMax" value="100" />
+        <par type="DOUBLE" name="reactance_no_reclosing_delay" value="0" />
+        <par type="DOUBLE" name="transformer_currentLimit_maxTimeOperation" value="90" />
+        <par type="DOUBLE" name="transformer_t1st_HT" value="60" />
+        <par type="DOUBLE" name="transformer_t1st_THT" value="30" />
+        <par type="DOUBLE" name="transformer_tNext_HT" value="10" />
+        <par type="DOUBLE" name="transformer_tNext_THT" value="10" />
+        <par type="DOUBLE" name="transformer_tolV" value="0.014999999700000001" />
+    </set>
+</parametersSet>
diff --git a/dynawaltz/solver.par b/dynawaltz/solver.par
new file mode 100644
index 0000000..339c440
--- /dev/null
+++ b/dynawaltz/solver.par
@@ -0,0 +1,12 @@
+<?xml version='1.0' encoding='UTF-8'?>
+<parametersSet xmlns="http://www.rte-france.com/dynawo">
+    <!-- IDA order 2 solver -->
+    <set id="1">
+        <par type="INT" name="order" value="2" />
+        <par type="DOUBLE" name="initStep" value="0.000001" />
+        <par type="DOUBLE" name="minStep" value="0.000001" />
+        <par type="DOUBLE" name="maxStep" value="10" />
+        <par type="DOUBLE" name="absAccuracy" value="1e-4" />
+        <par type="DOUBLE" name="relAccuracy" value="1e-4" />
+    </set>
+</parametersSet>

From 1446716465ab1e6258d88de81a08f4f905a964d2 Mon Sep 17 00:00:00 2001
From: Nicolas PIERRE <nicolas.pierre@artelys.com>
Date: Thu, 15 Dec 2022 11:10:38 +0100
Subject: [PATCH 2/2] remove omega ref

---
 dynawaltz.ipynb | 140 +++++-------------------------------------------
 1 file changed, 12 insertions(+), 128 deletions(-)

diff --git a/dynawaltz.ipynb b/dynawaltz.ipynb
index e3a72e4..a690d23 100644
--- a/dynawaltz.ipynb
+++ b/dynawaltz.ipynb
@@ -10,7 +10,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 1,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -22,7 +22,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 2,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -39,15 +39,14 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 3,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
     "df_alpha_beta = pd.DataFrame.from_dict({\"static_id\": [network.get_loads().loc[l].name for l in network.get_loads().index],\n",
     "                             \"parameter_set_id\": [\"LAB\" for l in network.get_loads().index]})\n",
     "df_GSTWPR = pd.DataFrame.from_dict({\"static_id\": [network.get_generators().loc[l].name for l in network.get_generators().index],\n",
-    "                             \"parameter_set_id\": [\"GSTWPR\" for l in network.get_generators().index]})\n",
-    "df_omega_ref = pd.DataFrame.from_dict({\"generator_id\": [network.get_generators().loc[l].name for l in network.get_generators().index]})"
+    "                             \"parameter_set_id\": [\"GSTWPR\" for l in network.get_generators().index]})"
    ]
   },
   {
@@ -60,14 +59,13 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
     "model_mapping = dyn.ModelMapping()\n",
     "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.ALPHA_BETA_LOAD, df_alpha_beta.set_index(\"static_id\"))\n",
-    "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.GENERATOR_SYNCHRONOUS_THREE_WINDINGS_PROPORTIONAL_REGULATIONS, df_GSTWPR.set_index(\"static_id\"))\n",
-    "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.OMEGA_REF, df_omega_ref.set_index(\"generator_id\"))"
+    "model_mapping.add_all_dynamic_mappings(dyn.DynamicMappingType.GENERATOR_SYNCHRONOUS_THREE_WINDINGS_PROPORTIONAL_REGULATIONS, df_GSTWPR.set_index(\"static_id\"))"
    ]
   },
   {
@@ -80,7 +78,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -98,7 +96,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 6,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -116,7 +114,7 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 8,
+   "execution_count": null,
    "metadata": {},
    "outputs": [],
    "source": [
@@ -134,123 +132,9 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 9,
+   "execution_count": null,
    "metadata": {},
-   "outputs": [
-    {
-     "data": {
-      "text/html": [
-       "<div>\n",
-       "<style scoped>\n",
-       "    .dataframe tbody tr th:only-of-type {\n",
-       "        vertical-align: middle;\n",
-       "    }\n",
-       "\n",
-       "    .dataframe tbody tr th {\n",
-       "        vertical-align: top;\n",
-       "    }\n",
-       "\n",
-       "    .dataframe thead th {\n",
-       "        text-align: right;\n",
-       "    }\n",
-       "</style>\n",
-       "<table border=\"1\" class=\"dataframe\">\n",
-       "  <thead>\n",
-       "    <tr style=\"text-align: right;\">\n",
-       "      <th></th>\n",
-       "      <th>_LOAD___2_EC_load_QPu</th>\n",
-       "      <th>_LOAD___2_EC_load_PPu</th>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>timestamp</th>\n",
-       "      <th></th>\n",
-       "      <th></th>\n",
-       "    </tr>\n",
-       "  </thead>\n",
-       "  <tbody>\n",
-       "    <tr>\n",
-       "      <th>0</th>\n",
-       "      <td>0.126992</td>\n",
-       "      <td>0.216992</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>0</th>\n",
-       "      <td>0.126992</td>\n",
-       "      <td>0.216992</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>0</th>\n",
-       "      <td>0.126992</td>\n",
-       "      <td>0.216992</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>0</th>\n",
-       "      <td>0.126992</td>\n",
-       "      <td>0.216992</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>0</th>\n",
-       "      <td>0.126992</td>\n",
-       "      <td>0.216992</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>...</th>\n",
-       "      <td>...</td>\n",
-       "      <td>...</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>9900</th>\n",
-       "      <td>0.125207</td>\n",
-       "      <td>0.215157</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>12710</th>\n",
-       "      <td>0.125208</td>\n",
-       "      <td>0.215157</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>18331</th>\n",
-       "      <td>0.125208</td>\n",
-       "      <td>0.215158</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>28331</th>\n",
-       "      <td>0.125208</td>\n",
-       "      <td>0.215158</td>\n",
-       "    </tr>\n",
-       "    <tr>\n",
-       "      <th>30000</th>\n",
-       "      <td>0.125208</td>\n",
-       "      <td>0.215158</td>\n",
-       "    </tr>\n",
-       "  </tbody>\n",
-       "</table>\n",
-       "<p>252 rows × 2 columns</p>\n",
-       "</div>"
-      ],
-      "text/plain": [
-       "           _LOAD___2_EC_load_QPu  _LOAD___2_EC_load_PPu\n",
-       "timestamp                                              \n",
-       "0                       0.126992               0.216992\n",
-       "0                       0.126992               0.216992\n",
-       "0                       0.126992               0.216992\n",
-       "0                       0.126992               0.216992\n",
-       "0                       0.126992               0.216992\n",
-       "...                          ...                    ...\n",
-       "9900                    0.125207               0.215157\n",
-       "12710                   0.125208               0.215157\n",
-       "18331                   0.125208               0.215158\n",
-       "28331                   0.125208               0.215158\n",
-       "30000                   0.125208               0.215158\n",
-       "\n",
-       "[252 rows x 2 columns]"
-      ]
-     },
-     "execution_count": 9,
-     "metadata": {},
-     "output_type": "execute_result"
-    }
-   ],
+   "outputs": [],
    "source": [
     "res.status()\n",
     "res.curves()"
@@ -273,7 +157,7 @@
    "name": "python",
    "nbconvert_exporter": "python",
    "pygments_lexer": "ipython3",
-   "version": "3.8.10 (default, Jun 22 2022, 20:18:18) \n[GCC 9.4.0]"
+   "version": "3.8.10"
   },
   "orig_nbformat": 4,
   "vscode": {