{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Eigenvalue Analysis" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Run Eigenvalue Analysis" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "execution": { "iopub.execute_input": "2021-09-26T22:41:51.305680Z", "iopub.status.busy": "2021-09-26T22:41:51.305313Z", "iopub.status.idle": "2021-09-26T22:41:51.838406Z", "shell.execute_reply": "2021-09-26T22:41:51.838711Z" } }, "outputs": [], "source": [ "import numpy as np\n", "\n", "import andes\n", "from andes.utils.paths import list_cases, get_case" ] }, { "cell_type": "code", "execution_count": 2, "metadata": { "execution": { "iopub.execute_input": "2021-09-26T22:41:51.841517Z", "iopub.status.busy": "2021-09-26T22:41:51.840878Z", "iopub.status.idle": "2021-09-26T22:41:51.844084Z", "shell.execute_reply": "2021-09-26T22:41:51.844562Z" } }, "outputs": [], "source": [ "case_path = get_case('kundur/kundur_full.xlsx')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Pass the routine name `EIG` or `eig` to `andes.run`." ] }, { "cell_type": "code", "execution_count": 3, "metadata": { "execution": { "iopub.execute_input": "2021-09-26T22:41:51.846712Z", "iopub.status.busy": "2021-09-26T22:41:51.846110Z", "iopub.status.idle": "2021-09-26T22:41:52.746189Z", "shell.execute_reply": "2021-09-26T22:41:52.746774Z" } }, "outputs": [ { "name": "stderr", "output_type": "stream", "text": [ "Working directory: \"/home/hacui/repos/andes/examples\"\n", "> Loaded config from file \"/home/hacui/.andes/andes.rc\"\n", "> Loaded generated Python code in \"/home/hacui/.andes/pycode\".\n", "Parsing input file \"/home/hacui/repos/andes/andes/cases/kundur/kundur_full.xlsx\"...\n", "Input file parsed in 0.1985 seconds.\n", "System internal structure set up in 0.0309 seconds.\n", "-> System connectivity check results:\n", " No islanded bus detected.\n", " System is interconnected.\n", " Each island has a slack bus correctly defined and enabled.\n", "\n", "-> Power flow calculation\n", " Numba: Off\n", " Sparse solver: KLU\n", " Solution method: NR method\n", "Power flow initialized in 0.0042 seconds.\n", "0: |F(x)| = 14.9282832\n", "1: |F(x)| = 3.608627841\n", "2: |F(x)| = 0.1701107882\n", "3: |F(x)| = 0.002038626956\n", "4: |F(x)| = 3.745104027e-07\n", "Converged in 5 iterations in 0.0070 seconds.\n", "Report saved to \"kundur_full_out.txt\" in 0.0010 seconds.\n", "Initialization for dynamics completed in 0.0301 seconds.\n", "Initialization was successful.\n", "\n", "-> Eigenvalue Analysis:\n", " Positive 0\n", " Zeros 1\n", " Negative 52\n", "Eigenvalue analysis finished in 0.0051 seconds.\n", "Report saved to \"kundur_full_eig.txt\".\n" ] }, { "name": "stdout", "output_type": "stream", "text": [ "-> Single process finished in 0.5363 seconds.\n" ] } ], "source": [ "ss = andes.run(case_path, routine='eig')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Plotting in S-Domain" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "execution": { "iopub.execute_input": "2021-09-26T22:41:52.768563Z", "iopub.status.busy": "2021-09-26T22:41:52.768058Z", "iopub.status.idle": "2021-09-26T22:41:53.540092Z", "shell.execute_reply": "2021-09-26T22:41:53.540496Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAjkAAAG6CAYAAADqJ1eGAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAABI4UlEQVR4nO3df1BTd74//mdQxJ8kgNhef6HBatfaXQ20XWuLrYT+lNW1IHXvX5epZPfO7H7mzt2aMvfeD8veO5fCdu6d6ezdFnSYuTufrxaJtF1qf0hwi1vrbpVou7ZrbROJv9oKhBOsKEQ53z/YZIkESDi/cpLnY4YZyTl5v18GTJ6e8/5hEEVRBBEREVGCSdG6ACIiIiIlMOQQERFRQmLIISIiooTEkENEREQJiSGHiIiIEhJDDhERESUkhhwiIiJKSAw5RERElJAYcoiIiCghMeQQERFRQpqudQFSNTQ0QBAEmEwmuN1uVFZWwmQyaV0WERERaUzXIaeurg4VFRWhUCMIAnbu3Inm5mZtCyMiIiLN6fp2VVtbW9hVG5PJBEEQNKuHiIiI4oeur+SYTCYUFRWhubkZJpMJHo8HZrN53PMHBwcxODgY+n54eBg+nw9ZWVkwGAxqlExEREQSiaKIq1evYuHChUhJGf96jUEURVHFumQlCALy8vLg8Xiwa9cu5ObmoqKiYtzzf/GLX6C6ulrFComIiEgpFy5cwOLFi8c9ruuQA4wMPG5ra4PD4YDVag1d1Ynk9is5fr8fS5cuxYULF5Cenq5SxURERNG7desWAGDatGmq993c3IzS0lLV+51Mf38/lixZAkEQYDQaxz1P17er7HZ76HaVx+NBaWkp8vLy4Ha7I56flpaGtLS0MY+np6cz5BARUVwSRRFvv/02nn76adX73rBhQ1x/Pk421ES3A489Hg8EQYDVagUAmM1mdHZ2wmQyweFwaFwdERGRPAwGA/x+v9Zl6JKuQ06k21I2m039YoiIiBS0ZMkSDA8Pq97vyZMnVe9TTroNOVarFS6Xa8yU8c7OTpSUlGhTFBERkQIsFkvYmFKKjq7H5DQ3N6OmpgZZWVmhNXJqa2u1LouIiEhWXV1dCAQCWLt2rar9FhcXq9qf3HQdckwmE0MNERElvOzsbHz88ceq93vs2DEUFRWp3q9cdHu7ioiIKFlkZWVh9erVqvfb09Ojep9yYsghIiKKc9OmTcPly5dV7zcjI0P1PuXEkENERKQDwXE5anrkkUdU7U9uDDlEREQ6kJ2dDZ/Pp2qfr7/+uqr9yU3XA4+JiIiSxcaNG7UuQXd4JYeIiEgHent78fvf/17VPu+9915V+5MbQw4REZEOmEwmdHd3q9rnzJkzVe1Pbgw5REREOjB9+nSYzWZV+zx+/Liq/cmNIYeIiEiiyWY9yTUratGiRbh586YsbSUDhhwiIiIJAoEAiouL4fV6Ix73er0oLi6WJeicOXNG1QX6nnrqKdX6UgJDDhERkQSpqamor69HeXn5mKDj9XpRXl6O+vp6pKamSu4rOzsbV65ckdxOtFwul2p9KYEhh4iISKKcnBw0NjaGBZ1gwGlsbEROTo4s/axYsULVcTlff/21an0pgSGHiIhIBqODzpEjR2QPOMDIbCc1r67MmzdPtb6UwMUAiYiIZJKTk4Oqqips3LgRHR0dsgYcADAYDLh48aKsbU7kiSeeUK0vJfBKDhERkQwCgQC8Xi+qq6vR0dGB6urqsDE6cs2wMhqNGB4elqWtyTQ3N6vSj1IYcoiIiCQKBAKwWq3YsWMHGhsbUVBQEDZGR84ZVk888YRqIUfveLuKiIhIosuXLyMQCMBgMIQeC47R2bFjBwwGA/bu3SvLDKtPP/0UKSkpWLNmjeS2JrN69WrF+1ASr+QQERFJEAgEYLPZsG/fPuzdu3fMVHKDwYDp06dj4cKFsvSn5jRyo9GoSj9KYcghIiKSIDU1Fa2trcjJyYk4w2rv3r1wOp2yXMUBgPnz52PRokWytDWZY8eOqdKPUni7ioiISKLRAUbpGVZyhaVkwCs5REREMppohpVcTp06hVu3bsne7u0ee+wxxftQEkMOERGRTEavcnz7DCs5ZWZmore3V9Y2I/nss88U70NJDDlEREQyiLSNQ6TtHuTw8MMPIysrS7b2xqPmwoNKYMghIiKSKDjDKtI2DsGgY7PZZFsQcGBgAJ2dnbK0NZHZs2cr3oeSDKIoiloXoZX+/n4YjUb4/X6kp6drXQ4REelYIBCYcFDwZMdjMTQ0hN/97ncoKSmRpT29ifbzm1dyiIiIZDBZgJFzVtSMGTNgMplka288+/btU7wPJTHkEBER6dD999/P7R0mwZBDRESkQ3/84x8Vn2G1cuVKRdtXGkMOERGRDmVnZ6O7u1vRPhYsWKBo+0pjyCEiItKh3NxcxUPIBx98oGj7SmPIISIi0qF58+bpfh0bpTHkEBER6ZDBYFB8ReJNmzYp2r7SGHKIiIh0KiUlRdEZVm63W7G21cCQQ0REpFPbt29HSopyH+VKbC6qpulaFyAHu92O3NxcACObliXrCpBERJRc/vKXv2D69OlYtWqVIu3PmDFDkXbVouuQIwgCCgsL0d7eDpPJBJfLhby8PCTxThVERJRE0tPTce7cOcVCzjPPPKNIu2rR9e0qu92OsrKy0NLWFosFbW1t2hZFRESkkuzsbMyZM0ex9pubmxVrWw26vpLT0NAAt9sNj8cDj8cDq9UKq9U67vmDg4MYHBwMfd/f369GmURERIqYOXMmFi5cqFj7N2/eVKxtNej2So7H4wEAuFwuCIIAs9kMm80Gp9M57nNqampgNBpDX0uWLFGrXCIiIkW8//77is2wMpvNirSrFt2HHJPJBIvFArPZjNraWpSWlo77nMrKSvj9/tDXhQsX1CqXiIhIERkZGejr61Ok7aVLlyrSrlp0G3KC8vPzQ382mUwQBGHcqzlpaWlIT08P+yIiItKzBx54QLFxOe+//74i7apFtyFnvEtoJpMpdJWHiIgo0RkMBt0v2qcUXYccs9k8JtAIghB2dYeIiCiRpaWl4ezZs4q0XVBQoEi7atFtyAGA2tpaNDU1hb53OBywWq2wWCwaVkVERKSeWbNmKdb2pUuXFGtbDbqeQl5SUgKfz4e6ujoAQG9vL9fJISKipLNlyxaIogiDwSBru263G/fff7+sbapJ1yEHACoqKrQugYiISFPvvfcevv/97yMjI0PWdpXcF0sN+q6eiIiIkJ2djStXrsjebllZmextqokhh4iISOeWLVumyDTylpYW2dtUE0MOERGRzmVlZSEQCMje7uitkPSIIYeIiEjnDAYDjh07Jnu7XPGYiIiI4oIoirK2t3LlSlnbUxtDDhERUQLYvHmz7G1OtOm1HjDkEBERJYBvvvkG58+f17qMuMKQQ0RElABSU1NlX6F4w4YNsranNoYcIiKiBJCdnY2bN2/K2mZPT4+s7amNIYeIiCgBzJkzB3l5ebK2+fnnn8vantoYcoiIiBLEm2++KfsMKz1jyCEiIkoQ6enp6O/vl609butAREREcWHdunWYNm2abO299dZbsrWlBYYcIiKiBGE0GmUdLHzt2jXZ2tICQw4REVEC+eSTT2Rra9GiRbK1pQWGHCIiogQxZ84cfPvtt7K1t2bNGtna0gJDDhERUYIwGAwoKSmRrb333ntPtra0wJBDRESUQD788ENZZ1jpGUMOERFRApk5cya6u7tlaeuBBx6QpR2tMOQQERElkMWLF8u2IKCc43u0wJBDRESUQBYtWoSMjAxZ2vr0009laUcrDDlEREQJRu8DhuXCkENERJRADAYDAMhyy0rOmVpaYMghIiJKMAUFBbKEnEOHDslQjXYYcoiIiBLM8PAwrly5IrkdvU9FZ8ghIiJKMENDQ/B6vZLbueOOO2SoRjsMOURERAkmOzsbfX19ktvJz8+XoRrtMOQQERElmPT0dDz++OOS2zl48KAM1WiHIYeIiCjB3Lx5E83NzeMOPg4EAipXpA2GHCIiogQSCARQXFyMQCCAgYGBMce9Xm/o+GTy8vKUKFE1DDlEREQJJDU1FfX19Xj99dfR1dUVdszr9aK8vBz19fVITU2dtC29X/FhyCEiIkowOTk5+M///E+88MILoVlWwYDT2NiInJycqNr55JNPlCxTcQw5RERECWjBggX4+7//e5SXl+PIkSMxB5xEMF3rAoiIiEh+RqMRw8PDqKqqwsaNG9HR0RFzwNm6dasyxakkoa7kFBUVaV0CERFRXDAYDFizZg2qq6vR0dGB6urqmBcIPHLkiELVqSNhQo7D4YDT6dS6DCIiorjg9XpRV1eH3/zmNygoKEBjYyPKy8tjCjo+n0/BCpWXECFHEATd/yCIiIjkEhxk/Oyzz2LmzJkARgYjxxp0srKylCxTcQkRcvbv34/t27dPet7g4CD6+/vDvoiIiBJJIBCAzWZDY2MjvvOd7+DatWuhY8GgY7PZopoe/tBDDylZquIMohx7sWvI6XTCbDYjMzMTGRkZE24t/4tf/ALV1dVjHvf7/UhPT1eyTCIiItUEAgGkpqZCFEUEAgHMmDEj4vHJ7Nu3Dzt27FCqzCnr7++H0Wic9PNb91dyBEGA2WyO6tzKykr4/f7Q14ULFxSujoiISH2jA8yBAwcmPJ7IdD2FvKGhARUVFVGfn5aWhrS0NAUrIiIiih8Gg2HCOxyTWbt2rXzFaEC3V3JcLpfut4AnIiJS2tq1azE8PDyl56ak6DYmANDxlRyfzweXyxWaNu52uwEAdXV1MJvNKCkp0bI8IiKiuDB//nz4/X5kZGTE/FyXy4VVq1YpUJU6dBtyrFYrrFZr6HuXy4WGhgbs2rVLw6qIiIjiS19fH65evZqUdz/0fR3qrxwOB2pqagAAdrudiwISERH9VXZ2Nrq7u6f03M2bN8tcjbp0P4VciminoBEREemVKIq4ceMGZs2aFfNz29vbUVhYqEBV0iTNFHIiIiIan8FgQEdHx5See+XKFZmrURdDDhERUYIbGBjAjRs3Yn6e0WhUoBr1MOQQEREluOXLl4dt7xCteLxVFQuGHCIiogS3Zs2aKS2G29LSokA16mHIISIiSnD9/f344IMPtC5DdQw5RERECS4jIwN9fX0xP2/NmjUKVKMehhwiIqIEl5KSggceeCDm582ePVuBatTDkENERJQEbty4gcHBwZie89FHHylUjToYcoiIiJJAT0/PlFc+1iuGHCIioiSQnZ2Nnp6emJ7z5JNPKlSNOhhyiIiIksCqVavwve99L6bnfPzxxwpVow6GHCIioiRgMBjwxhtvxPScy5cvK1OMShhyiIiIkoDBYMD169djes7cuXMVqkYdDDlERERJwmw2Y3h4OOrzn3rqKQWrUR5DDhERUZJYvXp1TFdz9u/fr2A1ymPIISIiShIXL17EF198oXUZqmHIISIiShLZ2dkxrZVz9913K1iN8hhyiIiIkkRmZmZM08gzMzMVrEZ5DDlERERJYtq0afjyyy+jPv/DDz9UsBrlMeQQERElkUuXLmFoaEjrMlTBkENERJRE7rzzTvj9/qjOLSoqUrgaZTHkEBERJZGHH34YWVlZUZ175swZhatRFkMOERFRErly5QoOHz4c1bkXLlxQuBplMeQQERElkczMTPT29kZ17qxZsxSuRlkMOURERElk+vTpWLVqVVTnbt26VdliFMaQQ0RElGQyMjIQCAQmPe+1115ToRrlMOQQERElGbfbjZ6enknPE0VRhWqUw5BDRESUZKLd3uGuu+5SoRrlTNe6ACIiIlLXypUro7pKc+edd6pQjXJ4JYeIiCjJzJgxI6otG/7whz+oUI1yGHKIiIiSjMFgwNdff611GYpjyCEiIkpC8+fPx/Dw8ITnPProoypVowyGHCIioiT0yCOP4NatWxOe09XVpU4xCol64HFLSws++ugjGAyGmDsRRREGgwE1NTUxP5eIiIjkd+bMGRgMBtx7773jnnPu3Dl8//vfV7EqeUUdctra2vDKK69MuaPt27dP+bkTqaurAzAy5x8A6uvrFemHiIgokWRnZ+Mvf/nLhOekpqaqVI0yog45FotFUkdKbNdut9tRW1sb+t5ms6GoqAhtbW2y90VERJRI5s+fj2XLlk14TklJiTrFKCTqMTk7d+6U1JHU599OEAS4XC4IghB6zGazwel0wuPxyNoXERFRoklNTcXAwMCE5zgcDpWqUYbkgcctLS1y1DElJ06cCAs0ZrMZAMKCDxEREUX26aef4ubNm+Mej2Z/q3gW1e2ql156adxjbW1t2LZtm2wFRctkMqGvry/sMafTCeBvYed2g4ODGBwcDH3f39+vXIFERERxLisrC729vbjjjjsiHl++fLnKFckrqpAjiiJ6e3uRm5s75pjRaJS9qKmqqalBfX09TCbTuMerq6vVLYqIiChObdy4EdOnjx8FJhuzE+8MYhSbV/j9fjidTjzzzDNjjrW3t6OwsFCR4mJht9uRlZWFXbt2jXtOpCs5S5Ysgd/vR3p6uhplEhERxY2enh6cPXsWDz74YMTj+/btw44dO1SuanL9/f0wGo2Tfn5HdSXHaDRGDDgA4iLgOBwO5ObmoqKiYsLz0tLSkJaWplJVRERE8S09PR2XLl3SugzF6H7F4+A4nGDAEQSBs6uIiEgXJhvYq/TA3xkzZow7HgcAHn74YUX7V9qUQ85Eg5HV4nK54HK5YLFY4PF44PF40NDQgMzMTK1LIyIimlAgEEBxcTG8Xm/E416vF8XFxYoHnTVr1oy7vYPeN/Gccsjp6emRs46YCYKAwsJC2O125Obmhr7sdvu4A4+JiIjiRWpqKurr61FeXj4m6Hi9XpSXl6O+vl7xVYdPnDgBn88X8dgXX3yhaN9Km3LImcoeVnIKTiEXRXHMFxERkR7k5OSgsbExLOgEA05jYyNycnIUr2HBggW4cuVKxGNaf9ZLpfsxOURERHo2OugcOXJE1YADACtWrMDChQsjHnv22WdVqUEpDDlEREQay8nJQVVVFTZu3IiqqirVAg4AzJkzJ7TJ9e3eeOMN1epQAkMOERGRxrxeL6qrq9HR0YHq6upxByMrwWAw4OzZsxGPXb9+XbU6lDDlkMMZTERERNKNHoNTUFAwZoyOGmbMmIHh4eExjy9ZskS1GpQQ1YrHiSraFROJiIiUMN4gY7UHHwejwO0DjXt6ejB//nzF+49VtJ/fvF1FRESkgUAgAJvNFjHIBAcj22w2VXYC//Of/4wzZ86MebytrU3xvpUkKeS0t7ejrKws9H1LSwsOHz4suSgiIqJEl5qaitbW1nGv1OTk5KC1tVXxdXKAkSEo400j1zPJV3JefPFFAMB9990Hu92Ozs5OlJWVob+/X3JxREREiWyyAKNGwAGA7OxsZGRkjHl8vI079SKqDTrH4/f7sXz5cpw7dw4ulwudnZ1Yu3YtAGDPnj147rnn5KiRiIiIFJSWlhZxtwCfz6fqdHa5SbqS43a70dXVFdpKIRhwAERMhERERBSfPvzwwzEzrCKN09ETSSGnoqICL774IjweD/bv3w8A6OrqAqD/paCJiIiSSXC7pEQi+xTyFStWoLS0FPfddx+2bdsmZ9Oy4xRyIiKiEf39/Zg5cyZmzJgReuzWrVuYNm2ahlVFFu3nt6QxOZF8+eWX8Pv9MBqNcjdNRERECrl16xbOnDmD7373u6HH3n77bRQXF2tYlTRR366KZbZUpIDD2VZERETxa+bMmfjyyy/DHvv22281qkYeUYccu90uqSOpzyciIiLlzJo1K+xWFYBxdyfXi6hvV4miiMrKyil1Ioqi7ldNJCIiSnSPPfYYRFEMTR763ve+p3FF0sQ08Njv90vqLN7G6XDgMRER0d+88847eOCBB0KbcO/btw87duzQuKqxFBl4HG8hhYiISE8CgcCEqxhPdlxp2dnZ6O7uDoUcveMGnURERCoIBAIoLi6G1+uNeNzr9aK4uFiVDTnHs3z58rArI/fff79mtciBIYeIiEgFqampqK+vR3l5+Zig4/V6UV5ejvr6ek2v5GRkZITNhh4YGNCsFjkw5BAREakkJycHjY2NYUEnGHAaGxs13ycqJSUFnZ2doe9Pnz6tYTXSMeQQERGpaHTQOXLkSNwEnCCDwQCZN0PQjKRtHbq6urBs2TIZy1EXZ1cREZFWjhw5go0bN6KjowMFBQValxNy/fp1zJw5EwaDAYODg0hLS9O6pDGi/fyWdCXHZrNJeToREVFS8nq9qK6uRkdHB6qrq8cdjKwFr9cb2my7vb1d22IkkhRy3G43XnrpJbS0tHDbBiIioiiMHoNTUFAwZozOZCabfSV1dtbMmTNx+fJlANLXx9OapJDT1taGn//859i2bRuOHz+O3bt3o6WlRa7aiIiIEkqkQcaRBiOPNnqGU6Rp6KOPyzENfcGCBaGdxxcsWDDlduKBpJCzfPny0J8LCwuxc+dOiKKIrKwslJWV4fDhw5ILJCIiSgSBQAA2my3iIONg0LHZbGEBZWBgAEuXLsXRo0cBjJ2GfvToUSxduhQDAwOyTUOfPXs27r77bgD6XydH0sDjw4cPY9OmTTh16hSamprQ0NCAzMxMVFRUoKSkBJmZmdi/fz+ysrKwbds2OeuWBQceExGRmqay4vHRo0dRXFyM1tZWbNiwAcDIFZtt27bB7Xbj4MGDWLx4sayztPbu3YsdO3bgtddeS55tHW5XUlKCrKws+Hw+VFRUwOl0Yt26dWHn7Ny5EwCwZ88ePPfcc1K6IyIi0rXJrrBEOr5hwwa0traGBZ2LFy/C7XYjNzcXt27dkn0autFoTIixtpJCjtlsRm1tLQoLCyc8b8+ePVK6ISIiSmqjg051dTWqqqpw8OBB3Lp1KzQNXc51dvLy8pCamgqLxSJbm1qQNCansrJy0oDj9/tx4sQJ5OfnS+mKiIgoqW3YsAHV1dX42c9+hurqaixevFixaeizZ8/GV199heHhYdna1IKkkHP8+HFUVlZOeI7RaMSrr76KtWvXSumKiIgoqR09ehRVVVV4+eWX8W//9m/Ytm3blKehTyYlJQWnT5/GqVOnZGlPK5K3dSgrK4v4eHAhISIiIpJm9ODjH/zgB8jNzYXb7cbFixcBTD4NPVZz5szR/eacgMSQU1RUBJ/PF/GYw+GQ0jQREREhPODcf//9sNlsaGlpwcGDB1FcXByaXj7eNPSpMBgMKCsrw5YtW+T4K2hG0hTy7du3w+Px4OTJkzCbzTCZTAAAURRx7tw59Pb2ylXnuOrq6kL9CoKAXbt2Rf1cPUwhHz2dMNLUwsmOk7ymMv0zEWuYiuCb7ni1TXY8Hqn5sxj9oRWpzcmOy0mL38Gp9tnT04P58+eP+7zJjmstuE7Om2++GZo+PvrvevToUWzZsgXnz5/H7NmzxxwP8vv9MBqN4/Yz+vj1wHX89uPf4p3P3sHf9f0dLi64iB/e/UP8/b1/j7Tp8bGPlSp7V3k8HlRWVuLQoUN49dVX8eKLL4a+8vLypDQdlbq6OgBARUUFKioqYLFYEmo/rdErW0Za5XL0ypZyrHJJE4v0MxhNjZ9BPNQwFYFAAE8//TSsVmvE2r1eL6xWKzZv3hx3tY9HzZ9FIBDA5s2bUVhYGPE1DL5+mzZtwtNPP51wv4NT7bOnpwcLFizAq6++GvF5r776KhYsWICenh7ZapXb7Nmzcf78+VDAAcJD7IYNG8ICzu3HgZEAk52djQMHDkTs48CBA8jOzobf74f/hh8/bPohZqfORtOPmlCwoAD7S/ZjWBzGM/ufwbWhazL/DRUmSuB0Oqd0TC4mk0ns6+sLeyyWv5Lf7xcBiH6/X+bK5NPV1SVu2rRJ7OrqiurPpKzxXms1fwbxUMNUdHV1ievXrxcffPDBsBrHe1wP1PxZdHV1iQ8++KBosVjE9evXh9oe73ElafE7ONU+X3nlFdFgMIivvPJKVI8nKofDIU6fPl10OBwTPr7zdzvFYxeOhY6/++67oT873U7x/7zzf1SpdzLRfn5LCjmRnDt3Tjxw4IB44MABuZsO43a7IwYaAGJbW1tUbegh5Iji2KDz4IMPht7M4v2DLRHd/ppr8TOIhxqm4vZAo+eAE6Tmz+L2QPPBBx+oHnBG16L27+BU+7w90CRbwAm6PdDc/r1vwCdufW1r2HMGBgbCvn/q/3tKvDp4VZ2CJ6BayPH7/eKBAwfE3bt3h31t375datMTamtrixhyTCaT2NzcHPE5N27cEP1+f+jrwoULugg5ovi3f8wdHR2hD4WOjg5dfLAlotE/D61+BvFQw1QEg43FYhHXrVun64ATpObPIhh07rnnHnHu3LniPffco3rAGV2L2r+DU+0zGGyeeeaZpAw4QcFg87Of/WzMlZ2Org6x9oPasPP37t0b9v3/Pfx/xc7LnarUOpFoQ46kgccnT55EaWkpzGYzfD4fzGYzBEFAX18fmpubsWzZMsm308bjdDpRVFSE28vPzc2F3W5HRUXFmOf84he/QHV19ZjH9+zZg9mzZ2Pbtm1ob2+H3+/HggULcP/99+Ott94CAFgsFgwPD4fWDNiyZQs++OAD9Pb2IjMzEwUFBXjjjTcAAN/97neRmpqKzs5OAMDTTz+NEydO4JtvvkF6ejoee+yx0Oyze+65B3PnzsWf/vQnAMDjjz+O06dP49KlS5gzZw42b96MpqYmACOj3e12O55//nkAwC9/+Uv813/9F6ZNm4a0tDRs27YNTU1NGB4eRm5uLhYtWoQjR44AAB555BGcP38eHo8H06dPR2lpKQ4cOIChoSHk5OQgNzc3tKHqQw89hCtXruDs2bMAgB07duDNN9/EwMAAFi9ejNWrV+PQoUMAgPXr18Pv9+Ozzz4DAJSWluLdd9/F1atXceedd8JiseDtt98GANx33324ceMG/vznPwMAfvjDH+L9999HX18f5s+fj/Xr16O1tRUAQtuDnDx5EgBQXFyMY8eOoaenBxkZGXjkkUfw+uuvAwDuvfdezJw5E8ePHwcAPPXUU3C5XPj6668xb948PPHEE2hubgYArF69GkajEceOHQMAPPbYY/jss89w8eJFzJ49G1u2bMG+ffsAACtXrsSCBQvwwQcfAAA2bdoEt9sNr9eLGTNmIDs7G7/+9a9RUFCA++67D0uXLsX7778PACgoKMClS5fgdruRkpKCsrIytLS0YHBwEEuXLsXKlSvhdDoBjNxT7+npweeffw5gZFmGt956C9euXcOiRYuwZs0avPfeewCABx54AN9++y0+/fRTAMCdd96Jl156CU8++SRWr16N/Px8HDx4EMDIiqWBQACffPIJAGDr1q04cuQIfD4fsrKy8NBDD+HNN98EAKxduxYpKSlwuVwAgM2bN+Ojjz7ClStXYDQaUVhYiJaWFgDAmjVrMHv2bHz00UcAgCeffBIff/wxLl++jLlz5+Kpp57C/v37AQB33303MjMz8eGHHwIYmZF55swZnDp1Cu3t7XjjjTfw61//GpmZmbjrrrtw55134g9/+AMA4NFHH0VXVxfOnTuH1NRUlJSUwOFwIBAIYPny5Vi2bBl+//vfAwAefvhhfP311/jiiy9gMBjw7LPP4o033sD169exZMkS3H333WhrawMAPPjgg/D5fDhz5gyAkQkUb7/9Nr799lssXLgQ3/ve9/DOO+8AGNmccGBgAKdPnwaACd8j0tLS8K//+q/4x3/8R2RlZSn6HvHWW2/h3XffxbVr19Da2oqXX34ZWVlZWLVqFebPnx+abWO1WnH27FmcP39esfeI3t5e/PSnP8Wrr76K9PR0Vd4jent78T//8z/4j//4DwwNDUX9HvHKK6/gq6++giiKqKioUPw94plnnkFzczNu3rwJs9msyXtESUkJDh06hP7+ftxxxx3Iz89HbW0tzp07h4ULF+If/uEfQu8Ri+9fjPbft2PVnFWh94jf/OY3WLRoUeg94v+99/+wJH0JdjyzQ9H3iAsXLmDWrFnYunUrXnvtNYiiGPYeMTAwgOeee27yiUNSktSPf/zj0J9dLpcoCELo+/b2dilNT2qiKzn19fURn8MrOSSXeLiKEg81TAWv5Ejvi1dyeCVnqia6kjN0c0gs/N9CcXh4OPTY6dOnQ3++eeumWPi/heKt4Vuq1hyJKrerRr84giCIu3fvDn2vdMjhmByOydFKPIyHiYcapoJjcqT3xTE5HJMzVZONyRFFUfzvY/8tvnT0pdD3n3/+eejPv/j9L8TdnX/7nNeSKiGnrq5OPHfuXCjclJaWiqdOnRJFURRfeOEFKU1HxWQyiW63O+yxWC5O6SHk3B5wOLtKW1rMKonHGqZivECj56Cj5s/i9oAz+oOes6s4u2oy0c6uGh4eFiudleIOxw7xsOewuOd/94iHvjwkluwvEX/5/i/DrvJoSZWQIwiCaLfbQ4OMOzs7xYyMDDElJUWVkFNbWxt2a6q5uVmsqKiI+vnxHnKGhobExx9/XOzq6gr7c1BXV5f4+OOPi0NDQ2F/JmVE+hmMpsbPIB5qmIqhoSGxqKhILCgoiFh7V1eXWFBQID722GNxV/t41PxZDA0NiY899pj48MMPR3wNg6/fQw89JBYVFSXc7+BU++zu7p4wyASDTnd3t2y1KuHatWuSjguCIKampo4JOEEOh0NMTU0NG3LiFbziL5t+KVb+d6VY84ca8aL/YuyFK0iVgceR+P1+eDye0MBRpdXV1cFsNgMY2TC0trY26udyxWOKVTysNhwPNUwFVzyW3lcQVzyO/jhXPB4Ry4rHQXv37sXmzZvj8vMx2s9v2UNOUEtLC7Zt26ZE07LRQ8ghIqLkdvveVcXFxaivr8fFixdDjwcDkNfrhc1mQ2trq6SgKYoimpubMX/+fGzatEmuv4psVAs5p06dgsfjGbNRZ319fWhKb7xiyCEiIj0YHXQWL16Mbdu2we124+DBg2EBp7y8HI2NjcjJyZHcpyiKeO2117Bjxw7Jbckt2s/v6VI6eeGFF+BwOMI25wRGNsr0eDxSmiYiIqK/2rBhA1pbW1FcXIzq6mq43W7k5uZi8eLFAOQPOJ9//jlu3Lih+wsAkkJOVlYWvvzyy4jHfvWrX0lpmoiIiEbZsGEDqqur8bOf/Qwvv/wyfvCDH6C8vBxVVVWorq6WLeAAQHd3N5YuXYrVq1fL0p5WJO1CHhzwG0lwVV4iIiKS7ujRo6iqqsLLL7+MqqoqXLx4EVVVVdi4cSOqqqpkCzgAYDQakZ2dHVp5W68khZzRy3zf7qWXXpLSNBEREf3V6DE5P/3pT9Ha2oqnn34a//RP/4SOjg5UV1fD6/XK1p/JZMKsWbNka08rkkLOc889h4qKCkybNg133XVX6GvFihWw2+1y1UhERJQQRk/1j/b46IATHGS8ePFi5Obmwu12Y9q0aWhsbER5ebksQUcUxdB+XPfcc4/k9rQkaUwOMDKLKjMzM+wxURTx4osvSm2aiIgoYQQCgdD070i3liJN/x4YGMCWLVvGTBMvLy9HS0sLLl68GFonJxh0pI7NEQQhtGbO3Llzp9xOXJCy4qDT6Rz3mMvlktK0KuJ9xWMiIkosU9meYvSKxpFWfx59XI4VpwcHB0OrH+/du3fK7Sgp2s9vSberCgsLxz2m1orHREREepGTkzPm1tJk07+DKxkDI6tZt7a2hp03+nhOTo7khQDPnj2LW7duTfn58URSyJkIBx4TERGNNTroHDlyJOZbTJMFGKlbanz55ZeYOXMmAODxxx+X1JbWYh6TU1ZWhrKyMmzbtg0rVqyAwWAYc44oijh37hx+/vOfy1IkERFRIsnJyQlN/+7o6JB1+rdUKSkpoatDp0+fRkFBgcYVTV3MIWf58uWh1Y3NZnPEDTFFDjwmIiIal9frRXV1dWj6t5wL+UkhiiKefPLJ0PeXLl3SsBrpJO1dderUKaxduzbisZMnT8b9uBzuXUVERGq7fQyO3FsySCEIAo4ePYqnn34aAPC73/0OP/jBDzStKZJoP78ljcmpqakZ99i6deuwZ88e/OQnP0FLS4uUboiIiBJCpEATaTCyVrq7u5GdnR36fvPmzRpWI52kkFNRUQEA6OrqGnNs9+7daG5uRkVFBdxuN4MOEREltUAgAJvNFvGKTTDo2Gy2SRcMVNLcuXPDtmxqamrSrBY5SAo5LpcLK1asgNlsxrRp07Bnz57QsYaGBtTW1mLdunV4/vnn4Xa7JRdLRESkV5Gmf48mx/Rvqa5duzZmgV89kxRy3G43Ojs7MTw8jFu3bqGvry90Vaevry8sDU60mScREVEyUHr6t1QfffRR2KzpVatWaViNdJJCTn5+fmjpZ2Dk9pXL5QIwEnJGDwaKNNWciIiI4kNwHtLoz+v58+drVY4sJIWcEydOoKWlBV1dXTh16hQqKipgtVoBjIScq1evhs71eDzSKiUiIiJFbdmyJez7o0ePalSJPCSFnNraWrz22mswm83YuXMnrFYrmpqa8JOf/ATNzc147rnncPjwYfzkJz+BxWKRq2YiIiKS2fnz53H+/Hmty5CVpHVyJuP3+9HQ0ACLxTLhPlda4To5REREI44dO4Y77rgjbAzt7VPK40W0n98xr3g8ma6urtC4nG3btuH555+XuwsiIiKSmcFgwIIFC8IeO3v2bFyGnGhJDjn9/f1wOp3w+Xxhj7e1tWHbtm1SmyciIiIVrF69GnPnzg177Pz589iwYYNGFUknKeScPHkSpaWlMJvN8Pl8MJvNEAQBfX19aG5ulqtGIiIiUpAoinjrrbfwox/9KOzxtLQ0jSqSh6SQ09DQgC+//BLASOAxm82hKeWHDx/GsmXLJBdIREREyrp69SrmzZs35nG935GRNLsqOF0cGFnsj1dviIiIojfZFg5qbfEwffp05Ofnj3k8qbd18Hg86Orqwp49e2A0GnHo0CF8/PHHAEbG5BAREVFkgUAAxcXF427K6fV6UVxcrErQ+eabbzB79uwxjw8PDyvet5Ikb9D56quvhgLNCy+8gEcffRTTpk2TpTgiIqJElZqaivr6+oi7jwd3K6+vr1dlq4fTp09H3JkgNzdX8b6VJPs6OX6/Hx6PB+vWrZOzWUVwnRwiItJaMNAEdye//Xs17Nu3D88+++yYoHPp0iUsWrRIlRpiEe3nt6QrOZEYjUasW7cO/f39cjdNRESUcHJyctDY2Ijy8nIcOXJE9YADAGVlZRGv5Bw5ckS1GpQge8gJqqmpUappIiKihJKTk4Oqqips3LgRVVVVqgacq1evor29XbX+1CQp5Jw8eRL5+fnIysoK+8rMzERdXZ1cNRIRESU0r9eL6upqdHR0oLq6etzByEro7u4eswhg0COPPKJaHUqQtE5OTU0NKisrx2y+KYoiXnjhBUmFERERJYPbx+AEb12pdcsqJSUFS5cujXjs/Pnz+Lu/+zvFa1CKpJBTVlaGZ555JuKxyspKKU0TERElvEiDjNUOOrNmzRp3fyqPx4MHHnhA0f6VJOl21e37VY3W19cnpWkiIqKEFggEYLPZIgaZYNCx2WyKr5PT3t4ecdAxMLJIoJ5JnkIevGJz+1z65uZmvPfee1KanlRw3I/b7QYA1NfXx/R8TiEnIiItBQKBCdfBmey4HPbu3Ttmz6p4F+3nt6SI9uMf/xhOpxNmszkUNABAEAR0dnZKaXpSdrsdtbW1oe9tNhuKioq40jIREenGZAFG6YAjiiIeffTRcY8fOHBg3GEpeiAp5OTl5eHVV1+NeOxXv/qVlKYnJAgCXC4XBEGAyWQCMBJy8vLy4PF4YDabFeubiIgoUXzzzTe4cePGuMeHhoZUrEZ+ksbkZGZmjnvs+eefl9L0pE6cOAGPxxP6PhhsBEEY9zmDg4Po7+8P+yIiIkpW58+fx82bN8c9ruZ6PUqQFHJyc3Nx+PDhiMdeeuklKU1PyGQyoa+vL2zqutPpBIAJr+LU1NTAaDSGvpYsWaJYjURERPHu6tWr486sApJ876r8/HwIgoBz586FhQtRFHHu3DncunVLliKjkZeXB5vNhoqKinHPGRwcxODgYOj7/v5+LFmyhAOPiYgoKQ0PDyMlZfzrHfv27cOOHTtUrCg6qgw8BkZmNN1+20oURbz44otRt+FwONDU1DTpeZEWHgRGBiGXlZVNGHAAIC0tDWlpaVHXRURElMiampriMsTIRdKVnPb2dhQWFkY8dvLkSVV2Inc4HPD5fJMGnEg4hZyIiJLVwMAADh06hK1bt457zoULF+JyaIciu5B3dXWFfT9ewAHUuY8XHIcTDDiCIIQNRiYiIqLIbt26hTVr1kx4zpUrV1SqRhkxhRyHwxH1uQ0NDTEXEwuXywWXywWLxQKPxwOPx4OGhoYJZ3wRERHRiIGBgUn3pTp79qxK1SgjpttVK1asQFFRUVTnOp1OfPHFF1MubCKCIGD58uURp4vHcveNt6uIiChZvf3221i/fj0yMjLGPSepBh77fD4cP3580qslPp9vwn2tpApOISciIqKp8fv9oQV1xxOPAScWMYWcysrKqBf5U3LFYyIiIpLmySefHHdjzqA333wTW7ZsUaki+cU0JsdqtSpyLhEREann+vXr+OyzzyY9b2BgQIVqlBNTyIllSrga08eJiIgodt3d3VEt2Lt48WIVqlGOpG0diIiISH9u3LiBhQsXTnre6tWrVahGOQw5RERESWbJkiUT7vUYdOjQIRWqUQ5DDhERUZJ54403YlpyRa8YcoiIiJKMKIoTbswZtH79ehWqUQ5DDhERURIRRRF5eXlRnev3+xWuRlkMOUREREmkr68P8+bNi+rcaKaZxzOGHCIioiTS1dWl+ys00WLIISIiSiLd3d3Izs6O6tzS0lKFq1EWQw4REVES2bRpE7KysqI6991331W4GmUx5BARESWR9957b9I9q4KuXr2qcDXKYsghIiJKEoODgxgcHIz6/DvvvFPBapTHkENERJQkrl+/jhUrVkR9vsViUbAa5THkEBERJYmUlBR85zvfifr8t99+W8FqlMeQQ0RElCT++Mc/Js30cYAhh4iIKGn4fD5kZmZGff59992nYDXKY8ghIiJKEg8++CCmTZsW9fk3btxQsBrlMeQQERElgaGhIfT19cX0nD//+c8KVaMOhhwiIqIk0NPTE3PI0TuGHCIioiTg8/mi3s4h6Ic//KFC1aiDIYeIiCgJ3HPPPVi9enVMz3n//feVKUYlDDlERERJ4MCBAxBFMabn6P32FkMOERFREhgaGkJKSmwf+/Pnz1eoGnUw5BARESU4URRx1113xfy89evXK1CNehhyiIiIEty1a9eQm5sb8/NaW1sVqEY907UugIiIiJTl8XiQkpIS02rHiYBXcoiIiBLclStXYp4+DgDr1q1ToBr18EoOERFRgrNYLDAajVqXoTpeySEiIkpwp0+fjmnPqqCTJ08qUI16GHKIiIgSWCAQwNdff611GZpgyCEiIkpgV69excKFC6f03OLiYpmrURdDDhERUQKbN28eNmzYMKXnHjt2TOZq1JUwIaeoqEjrEoiIiOLOH/7wB3zzzTdTem5PT4/M1agrIUKOw+GA0+nUugwiIqK409PTg6ysrCk9NyMjQ+Zq1KX7kCMIAnw+n9ZlEBERxaXVq1cjNTV1Ss995JFH5C1GZboPOfv378f27dujOndwcBD9/f1hX0RERInq5s2bmDVr1pSf//rrr8tYjfp0HXKcTiesVmvU59fU1MBoNIa+lixZomB1RERE2ggEAgCA3t5enDt3btzjiU7XIUcQBJjN5qjPr6yshN/vD31duHBBweqIiIjUFwgEUFxcDK/Xi+7ubixYsCDsuNfrRXFxcVRB595771WqTFXExbYODocDTU1Nk55XWVkJi8UCAGhoaEBFRUVM/aSlpSEtLW1KNRIREelBamoq6uvrUV5ejvr6eixdujR0zOv1ory8HI2NjVGN05k5c6aSpSrOIIqiqHURsXK5XAAQCjyCICAjIwOx/lX6+/thNBrh9/uRnp4ue51ERERa8Xq9+Pd//3f8y7/8C5YvXx4WcHJycqJqY9++fdixY4fClcYu2s/vuLiSEyufzweXyxWaNu52uwEAdXV1MJvNKCkp0bI8IiIizeXk5GD9+vV47rnnUFVVherq6pgCTiLQ5ZWc27lcLuTl5fFKDhER0V+JoojDhw8jNTUVGzduREdHBwoKCmJqw+/3x+Xu5dF+fut64DEwMp6npqYGAGC327koIBEREUYGIC9ZsgTV1dXo6OhAdXU1vF5vTG0Eh4folS5vV41WUlLC21NERES3+fDDD9HQ0BC6RdXY2BjzmBy9716u+ys5REREFM7r9eK3v/0t/vmf/zkUaEYHnWiv6MybN0/JMhXHkENERJRAAoEAbDYbfvzjH+O73/1u2LFg0LHZbFGtk/PEE08oVaYqGHKIiIgSSGpqKlpbWzFjxoyIa+Hk5OSgtbU1qnVympublShRNQw5RERECSYlJQVnzpwZ9/hUN+zUG4YcIiKiBNPX14fMzEzJ7axevVqGarTDkENERJRg5s+fH9MG1uOJxzVyYsGQQ0RElGCOHj2K7u5uye0cO3ZMhmq0w5BDRESUYL766iuYTCaty9AcQw4REVGCWbhwIdLS0iS389hjj8lQjXYYcoiIiBLI8PAwcnNzZWnrs88+k6UdrTDkEBERJRCfz4dTp07J0tbFixdlaUcrDDlEREQJpLu7G9nZ2bK0NXv2bFna0YruN+gkIiKiv1m0aBFSUuS5hrFlyxZZ2tEKr+QQERElkL/85S+yXYHZt2+fLO1ohSGHiIgogbjdbtmu5OgdXwUiIqIEIYoiZs2aJVt7K1eulK0tLTDkEBERJQhRFLF161bZ2luwYIFsbWmBIYeIiChBnD17FidPnpStvQ8++EC2trTAkENERJQguru7MX/+fK3LiBsMOURERAli/vz5st5i2rRpk2xtaYEhh4iIKEGkpaVh5syZsrXndrtla0sLDDlEREQJQBRF/OlPf5K1Ta/XK2t7amPIISIiSgCCIMBkMsna5owZM2RtT20GURRFrYvQSn9/P4xGI/x+P9LT07Uuh4iIaMpu3ryJGzduYO7cuVqXorhoP795JYeIiCgBnD59Gjdu3JC1zebmZlnbUxtDDhERUQLo6uqSfdfwmzdvytqe2hhyiIiIEsDMmTNlDzlms1nW9tTGkENERKRzoiiioKBA9naXLl0qe5tqYsghIiLSOb/fj/fff1/2dpVoU00MOURERDrX3d2N7OxsrcuIO9O1LoCIiIikMZlMyMrKkr1dJW6BqYlXcoiIiHSup6cHRqNR9nYvXboke5tqYsghIiLSuZMnTyIlRf6PdO5dRURERJoRRREGgwEGg0H2tpUITmrS/bYOdrsdubm5AIDMzEyUlJRE/Vxu60BERHoniiICgYDu95mKRbSf37odeCwIAgoLC9He3g6TyQSXy4W8vDzoPLMRERHFpKurC9euXcOaNWtkb7ulpQXbtm2TvV216PY6lN1uR1lZWWjHVYvFgra2Nm2LIiIiUtlXX32FOXPmKNL24OCgIu2qRbdXchoaGuB2u+HxeODxeGC1WmG1Wid8zuDgYNgPrL+/X+kyiYiIFDVjxgwsWLBAkba54rEGPB4PAMDlckEQBJjNZthsNjidzgmfV1NTA6PRGPpasmSJGuUSEREpZtmyZYpdyVm5cqUi7apF1yHHZDLBYrHAbDajtrYWpaWlEz6vsrISfr8/9HXhwgU1yiUiIlKEKIo4dOiQYu1PdvEg3sXF7SqHw4GmpqZJz6usrITFYgl9n5+fH/qzyWSCIAhwOp3j3rZKS0tDWlqa9IKJiIjiwNWrVzFv3jyty4hbcRFySkpKYpr6Pd7W7yaTKXSVh4iIKNHNmDEDDzzwgGLtb9iwQbG21aDL21Vmsxlms3lMoBEEIezqDhERUSK7cOECUlNTFWu/p6dHsbbVoMuQAwC1tbVht7gcDgesVmvY7SwiIqJEdubMGUyfrtxNmc8//1yxttUQF7erpqKkpAQ+nw91dXUAgN7eXq6TQ0RESeX69euYO3eu1mXELd1v6yAFt3UgIiK9Cn58K7FnVdDw8HBc7l8V7ed3/FVOREREk/r222/xzjvvKNrHW2+9pWj7SmPIISIi0qHu7m5kZGQo2se1a9cUbV9pDDlEREQ6lJqaimXLlinax6JFixRtX2kMOURERDpkMBhwxx13KNqHEjubq4khh4iISIeOHDmi6KBjAHjvvfcUbV9pDDlEREQ6I4oiRFFUPOToHUMOERGRDj3++OOK96HklhFqYMghIiKSQSAQkHQ8Fl999RX6+vpka2883377reJ9KIkhh4iISKJAIIDi4mJ4vd6Ix71eL4qLi2ULOhcvXpSlncl8+umnqvSjFIYcIiIiiVJTU1FfX4/y8vIxQcfr9aK8vBz19fWybaZ548YNLFiwQJa2Ehm3deC2DkREJJNgoGlsbEROTs6Y7+UyNDSEGTNmyNbeeAKBgKK7nE8Vt3UgIiJSWU5ODhobG1FeXo4jR44oEnAAwOFwyNreeA4dOqRKP0rR7S7kRERE8SgnJwdVVVXYuHEjOjo6ZA84165dw+zZs2Vtczz9/f2q9KMUXskhIiKSaPSAYq/Xi+rqanR0dKC6uhper1fWmVWiKGLt2rWytTcRpVdUVhpDDhERkQSjZ1aNHoNTUFCAxsZG/OhHP4LVapUt6Pj9fmRlZcnS1mTy8/NV6UcpDDlEREQSBGdW7dixAz/60Y/GjMERRRE3b97E5cuXZenv448/xs2bN2VpazIHDx5UpR+lcEwOERGRRAsXLkRqauqY21bl5eXYt28fAMBms6G1tVXybCW/3w+TySSpjWTBKeScQk5ERDIIBAK4fPkyysvLUVVVherq6rCrOnJNx7527RrmzJkjuZ1onD17FitXrlSlr1hE+/nNKzlEREQySE1NnXBmlRwB5/r16zhx4gQ2btwoua1oyDlgWgsck0NERCSTSDOr5HTlyhVMmzZN1jYn8sknn6jWlxIYcoiIiGQQaWZVpG0epLh58yYWLVokW3uJjiGHiIhIokjbN4xe/ViuoJOdnY1ly5bJ0lY0tm7dqlpfSmDIISIikiAQCMBms0XcviEYdGw2myzjWw4ePAg15wsdOXJEtb6UwIHHREREEqSmpk44NTwnJ0eWqePAyJo7KSnqXZ/w+Xyq9aUEXskhIiKSaLIAI1fA+f73vy+5nViotbKyUhhyiIiIdKC3t1eWsBSLhx56SNX+5MaQQ0REpAPnz5/HtWvXVO3zzTffVLU/uTHkEBER6UB3dzeys7O1LkNXOPCYiIhIB4qKimAwGFTtc+3atar2JzdeySEiItKBN998U/WQo+ZMLiXou3oiIqIkcOPGDQwPD6ver8vlUr1POTHkEBERxbnBwUGsWrVK6zJ0hyGHiIgozg0PD2PFihWq97t582bV+5QTQw4REVGcO378OK5evap6vx999JHqfcpJ17OrGhoaIAgCTCYT3G43KisrYTKZtC6LiIhIVj6fD5mZmar3e+XKFdX7lJNuQ05dXR0qKipCoUYQBOzcuRPNzc3aFkZERCSzjRs3Ytq0aar3azQaVe9TTrq9XdXW1hZ21cZkMkEQBM3qISIiUsLQ0BAuXbqkSd+FhYWa9CsX3YYck8mEoqKiULDxeDwwm83aFkVERCSz7u5ufPvtt5r03dLSokm/ctHt7ardu3cjLy8PGRkZ2LVrF3Jzc1FfXz/hcwYHBzE4OBj63u/3AwD6+/sVrZWIiGiqLl++jDlz5mjyWTUwMBCXn5HBmkRRnPhEUcfq6+vFkpISEYBotVrFvr6+Cc+vqqoSAfCLX/ziF7/4xa8E+Lpw4cKEn/sGUZwsBinP4XCgqalp0vMqKythsVgAAHa7HUVFRbBarfB4PCgtLYUgCHC73eM+//YrOYIgICcnB+fPn9f94Co19Pf3Y8mSJbhw4QLS09O1Lifu8fWKHV+z2PD1ig1fr9jE8+sliiKuXr2KhQsXTrj1RFzcriopKUFJSUnU53s8HgiCAKvVCgAwm83o7OxEXl4eHA7HuG2lpaUhLS1tzONGozHufoDxLD09na9XDPh6xY6vWWz4esWGr1ds4vX1iubihC4HHns8nojr4dhsNvWLISIioriky5BjtVrhcrnGTBnv7OyM6YoQERERJa64uF01Fc3NzaipqUFWVlZojZza2tqY2khLS0NVVVXEW1g0Fl+v2PD1ih1fs9jw9YoNX6/YJMLrFRcDj4mIiIjkpsvbVURERESTYcghIiKihMSQQ0RERAmJIYeIiIgSkm5nVynFbrcjNzcXAJCZmckp6eNwOp2or69HUVERzGYz2tracN999/H1ikFRURHa2tq0LiNuCYKA/fv3AwDcbjc8Hg92794dcY0sGlFXVwcAoZXfJ9vPj/72e9bc3Mx/j6PU1dWF/q0JgoBdu3ZpW9AUMeT8lSAIKCwsRHt7O0wmE1wuF/Ly8ibf/CtJCYIAp9MJh8MBs9kMu93OgBMDh8MBp9OpdRlxzW63w263w2w2AxhZ7LO0tJQfROOw2+1hy2jYbDYG6Um4XC6cOHECgiDA5/NpXU7cCIbliooKACP/qbXZbLoMzbxd9Vd2ux1lZWWh5GqxWPjmMIlz585BFEW43e7QPwaaHN9Qo+PxeOBwOELf5+bm4sSJExpWFL8EQRizQKrNZoPT6YTH49GusDhnsVhQUVERCtI0oqamJuw93Wq1oqGhQcOKpo4h568aGhpQUlICj8cT+h92cG8sIjnt378f27dv17qMuNfW1hZ2ifz48eP8NzmBEydOhAWa4Af37SvDE00kuDdkpNvCerz6zNtVQOiNweVywWw2w2w2hy6N8011fPv370dmZiZ8Ph/cbnfMK04nI6fTyd+pKXA4HBAEAc3NzVqXEpdMJhP6+vrCHgt+IPEqBcVivCt/wZ0F9IYhB3/7oZpMJlgsFgBAbW0tli9fPuaNg0YEX6fgG2hDQwNKS0v5ITQJQRBgNpt1+WahheCgUEEQUFpaykHHMaipqUF9fT1fM5JF8D+0epOQIcfhcKCpqWnS8yorK0Mf1gCQn58f+nMwtSbD/7yn8nrd/r/D7du3w2azjXuZM9FM5TVraGhI2rFLU/03aTKZQq9ZQ0MDMjIycO7cuYT/HZvq6xUUHGOYTL9vUl8zmpgeAw4AQCTR7XaLAMS+vr6wx00mk1hfX69NUXGuubl5zGMAxM7OTg2qiX+dnZ1hr01fX5/If37j6+vrE3ft2hX2bzL47zTS7x79TXNzM9+3YtTc3CxaLBaty4gLwX9ntwMgtrW1aVCRNBx4DITG4dx+L1IQhLCrOzQieOtg9OsVvP3C+/+R+Xw+OJ1O1NXVoa6uDna7HcDIVM3RM4hohMfjQV1dXdj/HoO/Y4l+FUeK4Dic4BUcQRA4u4piYjabYTKZIv7e6PGuRkLerpqK2tpaNDU1hS5jOhwOWK1WXtaMwGQyYdeuXWGBJjg7jR9AkVmt1rA3CJfLhYaGBt0usKU0i8Uy5ncs+O9Tj2+0anC5XHC5XKFZosDI+1gy3bKaKt3eilFIZWUlnE5n6HdHz79HBlHkandBDQ0Nof8t9vb2crbQBARBCFs3ga9X9IJjBxwOB3bt2oWioiJ+cEdw++9YcAYfg/RYgiBg+fLlEQe08y1+fMG1mJqamuByubBr1y6u3P5XdXV1of9kHD9+XLfv7ww5RERElJA4JoeIiIgSEkMOERERJSSGHCIiIkpIDDlERESUkBhyiIiIKCEx5BAREVFCYsghIiKihMQVj4mIJBi9aCFXsCaKL7ySQ0RJw2azwWazRVwZeKqcTid6e3vDHmtoaIDNZuO+ZEQaY8ghIll5PB7Y7XYYDAbk5uaGbUpqs9lk3TDS5XKhqKgIubm5k54rCAJMJhPq6+tl3RqipKRkTP8VFRWor6/H8ePHZeuHiGLH21VEJCuz2Yza2lq4XC6YzeawWzgejwd5eXlob2+XZfNbi8USCk9ERLfjlRwiUo3ZbEZ+fj527twpW5uZmZmytUVEiYVXcohIVSaTSdZbVkqqq6sbM94GALKysjjImEgHGHKISDWCIMDpdGL37t1jjtXV1cFsNsPj8cBsNqOkpCT0nIaGBpjNZrS1tcFms8lyq2t0TcH2g2N1rFYrAM6WItI7hhwiUozH4wnNMPJ4PHC73RHH45SWlqKsrCwUbIqKimA2m2GxWFBTUwObzRYKPrm5uejs7JRt8LDdbofdbg8FLLvdHgo50XA6nWhra4MgCGHhjIi0xzE5RKSY4Id+SUkJrFYrTpw4AZ/PF3ZOMAiNDgelpaWor68PHXc6nWFtjv5eqhMnTsBut4euIDU3N8f0fKvViubmZrS1tTHgEMUZXskhIlVYLBZUVlaitLQUfX19ocedTidMJlNYcHG73aFxO8HQIQgCPB4PfD7fmKAkRW1tLUpLS5GbmwuLxYLOzk7Z2iYibfFKDhGpxmw2h8JKUPA2j9VqDX3V1taira0NwMhaOKWlpdi/fz/MZjPMZrOsNVmtVvT19aGzsxOCIHABP6IEwpBDRKpzuVyhP1ssloizrQRBgCAIKCwsRGVlJSoqKmAymUKrFcsxQys3NzfUjsVigc1mk3WhQCLSFkMOESki0m2l4Aym4ErADQ0NsFqtyM/PH3MFZf/+/fB4PBAEIWygcrDN0UFpKgRBQElJSdiVoePHj8c06JiI4hvH5BCRrDweD+rr65Gfnw9gZGp48CqMyWRCc3Mzamtr4XA4QgGjra0NdrsdPp8vtLhfRUUFgJFp3Ha7HUVFRQBGxujY7XaUlZXB5XKhpqYmNCuqtrY26jpNJhPuu+++ULjy+XwRp7YTkX4ZRFEUtS6CiEhpgiCgpqYmpiAkVazBi4jkxdtVRERElJAYcoiIiCghMeQQUdJwOByw2WyhGVpKaWhoQGlpqaJ9ENHkOCaHiIiIEhKv5BAREVFCYsghIiKihMSQQ0RERAmJIYeIiIgSEkMOERERJSSGHCIiIkpIDDlERESUkBhyiIiIKCH9/+439JI1sRUtAAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" }, { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 4, "metadata": {}, "output_type": "execute_result" } ], "source": [ "ss.EIG.plot()" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Report\n", "\n", "Report is saved to the file and can be loaded into the notebook." ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "execution": { "iopub.execute_input": "2021-09-26T22:41:53.546062Z", "iopub.status.busy": "2021-09-26T22:41:53.545595Z", "iopub.status.idle": "2021-09-26T22:41:53.548872Z", "shell.execute_reply": "2021-09-26T22:41:53.548508Z" }, "scrolled": true }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "ANDES 1.7.5.post2.dev0+ga53a31fa\n", "Copyright (C) 2015-2022 Hantao Cui\n", "\n", "ANDES comes with ABSOLUTELY NO WARRANTY\n", "Case file: /home/hacui/repos/andes/andes/cases/kundur/kundur_full.xlsx\n", "Report time: 07/05/2022 08:57:50 PM\n", "\n", "Power flow converged in 5 iterations.\n", "Flat-start: No\n", "\n", "\n", "EIGENVALUE ANALYSIS REPORT \n", "\n", "Positives 0\n", "Zeros 1\n", "Negatives 52\n", "\n", "STATISTICS\n", " Most Associated Real Imag. Damped Freq. Frequency Damping [%]\n", "\n", "#1 LL_x EXDC2 1 -1 0 0 0 0\n", "#2 LL_x EXDC2 2 -1 0 0 0 0\n", "#3 LL_x EXDC2 3 -1 0 0 0 0\n", "#4 LL_x EXDC2 4 -1 0 0 0 0\n", "#5 LS_y EXDC2 3 -49.191 0.64233 0.10223 7.8296 99.991\n", "#6 LS_y EXDC2 3 -49.191 -0.64233 0.10223 7.8296 99.991\n", "#7 LS_y EXDC2 1 -49.535 0 0 0 0\n", "#8 LS_y EXDC2 4 -49.515 0 0 0 0\n", "#9 LA_y EXDC2 2 -49.202 0.39349 0.062627 7.831 99.997\n", "#10 LA_y EXDC2 2 -49.202 -0.39349 0.062627 7.831 99.997\n", "#11 LA_y EXDC2 1 -48.885 0 0 0 0\n", "#12 LA_y EXDC2 4 -48.905 0 0 0 0\n", "#13 e2d GENROU 2 -36.896 0 0 0 0\n", "#14 e2d GENROU 4 -36.75 0 0 0 0\n", "#15 e2d GENROU 1 -34.732 0 0 0 0\n", "#16 e2q GENROU 3 -33.89 0 0 0 0\n", "#17 e2q GENROU 2 -33.526 0 0 0 0\n", "#18 e2d GENROU 3 -31.242 0 0 0 0\n", "#19 e2q GENROU 4 -27.613 0 0 0 0\n", "#20 e2q GENROU 1 -24.954 0 0 0 0\n", "#21 delta GENROU 3 -0.6564 7.086 1.1278 1.1326 9.224\n", "#22 delta GENROU 3 -0.6564 -7.086 1.1278 1.1326 9.224\n", "#23 delta GENROU 2 -0.65252 6.8343 1.0877 1.0927 9.5046\n", "#24 delta GENROU 2 -0.65252 -6.8343 1.0877 1.0927 9.5046\n", "#25 delta GENROU 1 -0.19177 4.2247 0.67239 0.67308 4.5346\n", "#26 delta GENROU 1 -0.19177 -4.2247 0.67239 0.67308 4.5346\n", "#27 e1d GENROU 2 -5.4189 0 0 0 0\n", "#28 e1d GENROU 3 -5.3729 0 0 0 0\n", "#29 e1d GENROU 4 -4.1119 0 0 0 0\n", "#30 e1d GENROU 2 -3.6722 0 0 0 0\n", "#31 vp EXDC2 2 -0.48245 1.628 0.25911 0.27025 28.413\n", "#32 vp EXDC2 2 -0.48245 -1.628 0.25911 0.27025 28.413\n", "#33 vp EXDC2 4 -3.0736 0 0 0 0\n", "#34 vp EXDC2 1 -3.0635 0 0 0 0\n", "#35 vp EXDC2 1 -2.3787 0 0 0 0\n", "#36 LAG_y TGOV1 1 -1.5702 0 0 0 0\n", "#37 LAG_y TGOV1 4 -1.9922 0 0 0 0\n", "#38 LAG_y TGOV1 2 -2.0109 0 0 0 0\n", "#39 LAG_y TGOV1 3 -2.0194 0 0 0 0\n", "#40 W_x EXDC2 2 -1.1816 0 0 0 0\n", "#41 W_x EXDC2 4 -0.53573 0.73459 0.11691 0.1447 58.924\n", "#42 W_x EXDC2 4 -0.53573 -0.73459 0.11691 0.1447 58.924\n", "#43* delta GENROU 3 -1.1783e-14 0 0 0 0\n", "#44 LL_x TGOV1 2 -0.30739 0.44664 0.071085 0.086293 56.693\n", "#45 LL_x TGOV1 2 -0.30739 -0.44664 0.071085 0.086293 56.693\n", "#46 W_x EXDC2 1 -0.36003 0.38874 0.061869 0.084328 67.95\n", "#47 W_x EXDC2 1 -0.36003 -0.38874 0.061869 0.084328 67.95\n", "#48 W_x EXDC2 3 -0.35993 0.37956 0.060409 0.083251 68.81\n", "#49 W_x EXDC2 3 -0.35993 -0.37956 0.060409 0.083251 68.81\n", "#50 LL_x TGOV1 1 -0.14157 0 0 0 0\n", "#51 LL_x TGOV1 2 -0.14203 0 0 0 0\n", "#52 LL_x TGOV1 3 -0.14202 0 0 0 0\n", "\n", "EIGENVALUE DATA\n", " #1 #2 #3 #4 #5 #6 #7\n", "\n", "delta GENROU 1 0 0 0 0 0 0 4e-05\n", "delta GENROU 2 0 0 0 0 0 0 6e-05\n", "delta GENROU 3 0 0 0 0 0 0 2e-05\n", "delta GENROU 4 0 0 0 0 0 0 1e-05\n", "omega GENROU 1 0 0 0 0 0 0 2e-05\n", "omega GENROU 2 0 0 0 0 0 0 3e-05\n", "omega GENROU 3 0 0 0 0 0 0 1e-05\n", "omega GENROU 4 0 0 0 0 0 0 0\n", "e1q GENROU 1 0 0 0 0 0.00084 0.00084 0.00306\n", "e1q GENROU 2 0 0 0 0 0.0019 0.0019 0.00224\n", "e1q GENROU 3 0 0 0 0 0.0026 0.0026 0.00069\n", "e1q GENROU 4 0 0 0 0 0.00118 0.00118 0.00063\n", "e1d GENROU 1 0 0 0 0 0 0 0\n", "e1d GENROU 2 0 0 0 0 0 0 0\n", "e1d GENROU 3 0 0 0 0 0 0 0\n", "e1d GENROU 4 0 0 0 0 0 0 0\n", "e2d GENROU 1 0 0 0 0 0.00839 0.00839 0.02012\n", "e2d GENROU 2 0 0 0 0 0.01824 0.01824 0.01364\n", "e2d GENROU 3 0 0 0 0 0.02784 0.02784 0.00467\n", "e2d GENROU 4 0 0 0 0 0.01361 0.01361 0.00474\n", "e2q GENROU 1 0 0 0 0 0.00096 0.00096 0.00036\n", "e2q GENROU 2 0 0 0 0 0.00026 0.00026 0.00114\n", "e2q GENROU 3 0 0 0 0 3e-05 3e-05 0.00031\n", "e2q GENROU 4 0 0 0 0 0.00118 0.00118 9e-05\n", "LAG_y TGOV1 1 0 0 0 0 0 0 0\n", "LAG_y TGOV1 2 0 0 0 0 0 0 0\n", "LAG_y TGOV1 3 0 0 0 0 0 0 0\n", "LAG_y TGOV1 4 0 0 0 0 0 0 0\n", "LL_x TGOV1 1 0 0 0 0 0 0 0\n", "LL_x TGOV1 2 0 0 0 0 0 0 0\n", "LL_x TGOV1 3 0 0 0 0 0 0 0\n", "LL_x TGOV1 4 0 0 0 0 0 0 0\n", "vp EXDC2 1 0 0 0 0 0.00157 0.00157 0.00231\n", "vp EXDC2 2 0 0 0 0 0.00352 0.00352 0.00168\n", "vp EXDC2 3 0 0 0 0 0.00346 0.00346 0.00037\n", "vp EXDC2 4 0 0 0 0 0.00202 0.00202 0.00043\n", "LS_y EXDC2 1 0 0 0 0 0.06526 0.06526 0.53032\n", "LS_y EXDC2 2 0 0 0 0 0.14386 0.14386 0.37893\n", "LS_y EXDC2 3 0 0 0 0 0.30715 0.30715 0.18162\n", "LS_y EXDC2 4 0 0 0 0 0.08575 0.08575 0.10168\n", "LL_x EXDC2 1 0.66046 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0.66046 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 1 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0.65832 0 0 0\n", "LA_y EXDC2 1 0 0 0 0 0.06704 0.06704 0.22072\n", "LA_y EXDC2 2 0 0 0 0 0.20271 0.20271 0.21636\n", "LA_y EXDC2 3 0 0 0 0 0.30051 0.30051 0.072\n", "LA_y EXDC2 4 0 0 0 0 0.1453 0.1453 0.06982\n", "W_x EXDC2 1 0 0 0 0 5e-05 5e-05 0.00015\n", "W_x EXDC2 2 0 0 0 0 0.00012 0.00012 0.00013\n", "W_x EXDC2 3 0 0 0 0 0.00017 0.00017 4e-05\n", "W_x EXDC2 4 0 0 0 0 8e-05 8e-05 4e-05\n", "\n", "PARTICIPATION FACTORS [1/8]\n", " #8 #9 #10 #11 #12 #13 #14\n", "\n", "delta GENROU 1 2e-05 1e-05 1e-05 5e-05 2e-05 0.00044 0.00047\n", "delta GENROU 2 1e-05 0 0 6e-05 1e-05 0.00078 0.00049\n", "delta GENROU 3 5e-05 0 0 2e-05 6e-05 0.00062 0.00061\n", "delta GENROU 4 5e-05 1e-05 1e-05 1e-05 6e-05 0.00032 0.00067\n", "omega GENROU 1 1e-05 1e-05 1e-05 3e-05 1e-05 0.00025 0.00027\n", "omega GENROU 2 1e-05 0 0 4e-05 1e-05 0.00045 0.00028\n", "omega GENROU 3 2e-05 0 0 1e-05 2e-05 0.00025 0.00024\n", "omega GENROU 4 2e-05 0 0 0 2e-05 0.00012 0.00025\n", "e1q GENROU 1 0.00117 0.00225 0.00225 0.0031 0.00115 0.0003 0.00033\n", "e1q GENROU 2 0.00043 0.00285 0.00285 0.00224 0.00041 0.00042 0.00021\n", "e1q GENROU 3 0.00184 0.00167 0.00167 0.00067 0.00185 0.00028 0.00027\n", "e1q GENROU 4 0.00329 0.00134 0.00134 0.00061 0.00333 0.00019 0.00041\n", "e1d GENROU 1 0 0 0 0 0 2e-05 1e-05\n", "e1d GENROU 2 0 0 0 0 0 7e-05 2e-05\n", "e1d GENROU 3 0 0 0 0 0 5e-05 5e-05\n", "e1d GENROU 4 0 0 0 0 0 2e-05 1e-05\n", "e2d GENROU 1 0.00784 0.01923 0.01923 0.02323 0.00867 0.23812 0.24595\n", "e2d GENROU 2 0.00261 0.02337 0.02337 0.01559 0.00284 0.33371 0.16432\n", "e2d GENROU 3 0.01258 0.01507 0.01507 0.00518 0.01432 0.2465 0.23261\n", "e2d GENROU 4 0.02533 0.01309 0.01309 0.00526 0.02886 0.17228 0.34779\n", "e2q GENROU 1 0.00015 0.0011 0.0011 0.00039 0.00016 0.00111 0.00037\n", "e2q GENROU 2 0.00038 0.00045 0.00045 0.00121 0.00039 0.0033 0.00106\n", "e2q GENROU 3 0.00116 7e-05 7e-05 0.00032 0.00123 0.00241 0.00234\n", "e2q GENROU 4 0.00041 0.00079 0.00079 9e-05 0.00044 0.00069 0.00043\n", "LAG_y TGOV1 1 0 0 0 0 0 0 0\n", "LAG_y TGOV1 2 0 0 0 0 0 0 0\n", "LAG_y TGOV1 3 0 0 0 0 0 0 0\n", "LAG_y TGOV1 4 0 0 0 0 0 0 0\n", "LL_x TGOV1 1 0 0 0 0 0 0 0\n", "LL_x TGOV1 2 0 0 0 0 0 0 0\n", "LL_x TGOV1 3 0 0 0 0 0 0 0\n", "LL_x TGOV1 4 0 0 0 0 0 0 0\n", "vp EXDC2 1 0.00094 0.00415 0.00415 0.01364 0.00476 0.00032 0.00034\n", "vp EXDC2 2 0.00034 0.00522 0.00522 0.00979 0.00169 0.00042 0.0002\n", "vp EXDC2 3 0.00105 0.0022 0.0022 0.0021 0.00546 0.0002 0.00019\n", "vp EXDC2 4 0.00242 0.00226 0.00226 0.00246 0.01263 0.00018 0.00038\n", "LS_y EXDC2 1 0.19482 0.203 0.203 0.22143 0.08347 0.00069 0.00074\n", "LS_y EXDC2 2 0.06924 0.25052 0.25052 0.15607 0.02909 0.0009 0.00042\n", "LS_y EXDC2 3 0.46474 0.22878 0.22878 0.07262 0.2043 0.00092 0.00088\n", "LS_y EXDC2 4 0.51255 0.11281 0.11281 0.04075 0.22639 0.00041 0.00085\n", "LL_x EXDC2 1 0 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 0 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0 0 0 0\n", "LA_y EXDC2 1 0.08621 0.2062 0.2062 0.53692 0.19058 0.00079 0.00084\n", "LA_y EXDC2 2 0.04203 0.34908 0.34908 0.51915 0.09112 0.00141 0.00066\n", "LA_y EXDC2 3 0.19588 0.22134 0.22134 0.16772 0.44428 0.001 0.00095\n", "LA_y EXDC2 4 0.37416 0.18904 0.18904 0.16302 0.85269 0.00077 0.0016\n", "W_x EXDC2 1 6e-05 0.00014 0.00014 0.00038 0.00014 0 0\n", "W_x EXDC2 2 2e-05 0.00021 0.00021 0.00031 5e-05 0 0\n", "W_x EXDC2 3 0.00011 0.00013 0.00013 0.0001 0.00026 0 0\n", "W_x EXDC2 4 0.00021 0.00011 0.00011 0.0001 0.0005 0 0\n", "\n", "PARTICIPATION FACTORS [2/8]\n", " #15 #16 #17 #18 #19 #20 #21\n", "\n", "delta GENROU 1 0.0013 0.00058 0.00086 4e-05 1e-05 0.00014 0.04024\n", "delta GENROU 2 0.00034 0.00024 0.00268 6e-05 7e-05 2e-05 0.08541\n", "delta GENROU 3 0.00011 0.00276 0.00076 5e-05 5e-05 0.0002 0.23816\n", "delta GENROU 4 0.00079 0.00199 1e-05 0.00011 1e-05 1e-05 0.16086\n", "omega GENROU 1 0.00075 0.00034 0.0005 3e-05 0 8e-05 0.02377\n", "omega GENROU 2 0.0002 0.00014 0.00155 3e-05 4e-05 1e-05 0.05045\n", "omega GENROU 3 4e-05 0.00109 0.0003 2e-05 2e-05 8e-05 0.0965\n", "omega GENROU 4 0.00029 0.00074 1e-05 4e-05 0 0 0.06107\n", "e1q GENROU 1 0.00054 6e-05 0.00042 0.00094 4e-05 0.0003 0.00067\n", "e1q GENROU 2 0.00065 0.00022 2e-05 0.00115 0.00016 0.00097 0.00141\n", "e1q GENROU 3 0.00038 0.00019 0.00019 0.00144 0.00011 0.00123 0.00194\n", "e1q GENROU 4 0.00034 1e-05 0.0003 0.00112 2e-05 0.00045 0.0034\n", "e1d GENROU 1 0.00197 0.00266 0.00808 1e-05 0.02906 0.04674 0.00263\n", "e1d GENROU 2 0.0033 0.00186 0.01412 0.00441 0.01327 0.03129 0.00637\n", "e1d GENROU 3 0.00138 0.01461 0.00314 0.00479 0.0177 0.02276 0.02746\n", "e1d GENROU 4 0.00094 0.01041 0.00046 5e-05 0.03673 0.03898 0.00549\n", "e2d GENROU 1 0.25874 0.03867 0.05142 0.14321 0.00031 0.0033 0.00066\n", "e2d GENROU 2 0.16397 0.03747 0.0468 0.20085 0.00352 0.02285 0.00139\n", "e2d GENROU 3 0.11321 0.0137 0.07851 0.27348 0.00187 0.03733 0.00213\n", "e2d GENROU 4 0.18053 0.00845 0.08023 0.18763 0.00026 0.00906 0.00371\n", "e2q GENROU 1 0.05939 0.06822 0.19355 0.00017 0.23112 0.22378 0.00134\n", "e2q GENROU 2 0.09922 0.04763 0.3366 0.06848 0.10503 0.1491 0.00325\n", "e2q GENROU 3 0.04116 0.36911 0.07409 0.07356 0.13858 0.10724 0.01383\n", "e2q GENROU 4 0.028 0.26219 0.0109 0.00082 0.28669 0.18311 0.00276\n", "LAG_y TGOV1 1 0 0 0 0 0 0 0.00067\n", "LAG_y TGOV1 2 0 0 0 0 0 0 0.00142\n", "LAG_y TGOV1 3 0 0 0 0 0 0 0.00355\n", "LAG_y TGOV1 4 0 0 0 0 0 0 0.00239\n", "LL_x TGOV1 1 0 0 0 0 0 0 3e-05\n", "LL_x TGOV1 2 0 0 0 0 0 0 7e-05\n", "LL_x TGOV1 3 0 0 0 0 0 0 0.00014\n", "LL_x TGOV1 4 0 0 0 0 0 0 9e-05\n", "vp EXDC2 1 0.00033 3e-05 0.00062 0.00077 9e-05 0.00049 6e-05\n", "vp EXDC2 2 0.00084 0.0003 0.00017 0.00073 0.00025 0.00122 0.00018\n", "vp EXDC2 3 0.00034 0.00023 6e-05 0.00069 0.00014 0.00101 0.00069\n", "vp EXDC2 4 0.00022 4e-05 0.00029 0.0009 3e-05 0.00061 0.00083\n", "LS_y EXDC2 1 0.0006 5e-05 0.00101 0.00103 9e-05 0.00039 1e-05\n", "LS_y EXDC2 2 0.00149 0.00049 0.00026 0.00095 0.00024 0.00095 3e-05\n", "LS_y EXDC2 3 0.00129 0.00084 0.0002 0.00197 0.0003 0.00172 0.00021\n", "LS_y EXDC2 4 0.00041 6e-05 0.00048 0.00123 3e-05 0.00049 0.00012\n", "LL_x EXDC2 1 0 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 0 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0 0 0 0\n", "LA_y EXDC2 1 0.00067 6e-05 0.00113 0.00114 0.0001 0.00042 1e-05\n", "LA_y EXDC2 2 0.00229 0.00074 0.0004 0.00144 0.00037 0.00143 3e-05\n", "LA_y EXDC2 3 0.00137 0.00089 0.00021 0.00207 0.00031 0.00179 0.00018\n", "LA_y EXDC2 4 0.00076 0.00011 0.00089 0.00224 6e-05 0.00089 0.00018\n", "W_x EXDC2 1 0 0 0 0 0 0 0\n", "W_x EXDC2 2 0 0 0 0 0 1e-05 1e-05\n", "W_x EXDC2 3 0 0 0 0 0 1e-05 3e-05\n", "W_x EXDC2 4 0 0 0 1e-05 0 0 3e-05\n", "\n", "PARTICIPATION FACTORS [3/8]\n", " #22 #23 #24 #25 #26 #27 #28\n", "\n", "delta GENROU 1 0.04024 0.17926 0.17926 0.17245 0.17245 0.01877 0.00547\n", "delta GENROU 2 0.08541 0.21636 0.21636 0.08005 0.08005 0.02617 0.00332\n", "delta GENROU 3 0.23816 0.04918 0.04918 0.09817 0.09817 0.00681 0.02482\n", "delta GENROU 4 0.16086 0.08069 0.08069 0.14333 0.14333 0.00315 0.02115\n", "omega GENROU 1 0.02377 0.10604 0.10604 0.10491 0.10491 0.01037 0.00302\n", "omega GENROU 2 0.05045 0.12794 0.12794 0.04872 0.04872 0.01445 0.00183\n", "omega GENROU 3 0.0965 0.01995 0.01995 0.04102 0.04102 0.00257 0.00936\n", "omega GENROU 4 0.06107 0.03068 0.03068 0.05611 0.05611 0.00111 0.00747\n", "e1q GENROU 1 0.00067 0.00511 0.00511 0.00596 0.00596 0.00153 0.00043\n", "e1q GENROU 2 0.00141 0.00145 0.00145 0.00491 0.00491 0.00172 0.00027\n", "e1q GENROU 3 0.00194 0.00018 0.00018 0.00041 0.00041 0.00039 0.0016\n", "e1q GENROU 4 0.0034 0.00215 0.00215 0.00091 0.00091 0.00022 0.00155\n", "e1d GENROU 1 0.00263 0.00462 0.00462 0.00044 0.00044 0.36152 0.09703\n", "e1d GENROU 2 0.00637 0.02796 0.02796 0.00049 0.00049 0.49924 0.07963\n", "e1d GENROU 3 0.02746 0.00894 0.00894 0.00335 0.00335 0.12193 0.48105\n", "e1d GENROU 4 0.00549 0.00157 0.00157 0.00281 0.00281 0.05818 0.38742\n", "e2d GENROU 1 0.00066 0.0043 0.0043 0.00204 0.00204 0.00033 8e-05\n", "e2d GENROU 2 0.00139 0.00146 0.00146 0.00158 0.00158 0.00074 0.00013\n", "e2d GENROU 3 0.00213 0.00012 0.00012 0.00027 0.00027 0.00018 0.00073\n", "e2d GENROU 4 0.00371 0.00209 0.00209 0.00056 0.00056 6e-05 0.0004\n", "e2q GENROU 1 0.00134 0.00229 0.00229 0.00017 0.00017 0.01819 0.00503\n", "e2q GENROU 2 0.00325 0.01381 0.01381 0.00019 0.00019 0.02501 0.00411\n", "e2q GENROU 3 0.01383 0.00437 0.00437 0.00129 0.00129 0.00604 0.02453\n", "e2q GENROU 4 0.00276 0.00076 0.00076 0.00109 0.00109 0.00287 0.01969\n", "LAG_y TGOV1 1 0.00067 0.0032 0.0032 0.00728 0.00728 0.00125 0.00038\n", "LAG_y TGOV1 2 0.00142 0.00386 0.00386 0.00338 0.00338 0.00175 0.00023\n", "LAG_y TGOV1 3 0.00355 0.00079 0.00079 0.00371 0.00371 0.00041 0.00152\n", "LAG_y TGOV1 4 0.00239 0.00129 0.00129 0.0054 0.0054 0.00019 0.00129\n", "LL_x TGOV1 1 3e-05 0.00016 0.00016 0.00064 0.00064 6e-05 2e-05\n", "LL_x TGOV1 2 7e-05 0.0002 0.0002 0.00031 0.00031 8e-05 1e-05\n", "LL_x TGOV1 3 0.00014 3e-05 3e-05 0.00026 0.00026 1e-05 5e-05\n", "LL_x TGOV1 4 9e-05 5e-05 5e-05 0.00038 0.00038 1e-05 5e-05\n", "vp EXDC2 1 6e-05 0.00162 0.00162 0.00482 0.00482 0.00739 0.002\n", "vp EXDC2 2 0.00018 0.00133 0.00133 0.00482 0.00482 0.01079 0.00183\n", "vp EXDC2 3 0.00069 0.00024 0.00024 0.00056 0.00056 0.00174 0.00719\n", "vp EXDC2 4 0.00083 0.00076 0.00076 0.00082 0.00082 0.00101 0.00706\n", "LS_y EXDC2 1 1e-05 0.00022 0.00022 0.00047 0.00047 0.00037 0.0001\n", "LS_y EXDC2 2 3e-05 0.00018 0.00018 0.00046 0.00046 0.00052 9e-05\n", "LS_y EXDC2 3 0.00021 7e-05 7e-05 0.00012 0.00012 0.00018 0.00074\n", "LS_y EXDC2 4 0.00012 0.00011 0.00011 8e-05 8e-05 5e-05 0.00035\n", "LL_x EXDC2 1 0 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 0 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0 0 0 0\n", "LA_y EXDC2 1 1e-05 0.00021 0.00021 0.00039 0.00039 0.00064 0.00017\n", "LA_y EXDC2 2 3e-05 0.00023 0.00023 0.00052 0.00052 0.00126 0.00021\n", "LA_y EXDC2 3 0.00018 6e-05 6e-05 9e-05 9e-05 0.00031 0.00125\n", "LA_y EXDC2 4 0.00018 0.00016 0.00016 0.00011 0.00011 0.00015 0.00102\n", "W_x EXDC2 1 0 5e-05 5e-05 0.00038 0.00038 0.00056 0.00015\n", "W_x EXDC2 2 1e-05 5e-05 5e-05 0.00044 0.00044 0.00093 0.00016\n", "W_x EXDC2 3 3e-05 1e-05 1e-05 7e-05 7e-05 0.00022 0.00094\n", "W_x EXDC2 4 3e-05 3e-05 3e-05 9e-05 9e-05 0.00011 0.00076\n", "\n", "PARTICIPATION FACTORS [4/8]\n", " #29 #30 #31 #32 #33 #34 #35\n", "\n", "delta GENROU 1 0.00184 0.00076 0.0105 0.0105 0.00211 0.00572 0.00175\n", "delta GENROU 2 0.00619 0.00021 8e-05 8e-05 0.0019 0.00509 0.00601\n", "delta GENROU 3 0.00524 0.00086 0.00348 0.00348 0.00522 0.00212 0.00608\n", "delta GENROU 4 0.00254 0 0.00062 0.00062 0.00562 0.00191 0.00119\n", "omega GENROU 1 0.00097 0.00039 0.00811 0.00811 0.00097 0.00262 0.00049\n", "omega GENROU 2 0.00326 0.00011 0.00024 0.00024 0.00087 0.00234 0.00176\n", "omega GENROU 3 0.00188 0.0003 0.00176 0.00176 0.00163 0.00066 0.00119\n", "omega GENROU 4 0.00086 0 0.00023 0.00023 0.00164 0.00056 0.0002\n", "e1q GENROU 1 0.00024 0.00202 0.07076 0.07076 0.00382 0.01402 0.05544\n", "e1q GENROU 2 0.0016 0.0052 0.0874 0.0874 0.00481 0.01018 0.04937\n", "e1q GENROU 3 0.00126 0.00534 0.07194 0.07194 0.01035 0.00282 0.04397\n", "e1q GENROU 4 0.00048 0.00212 0.0524 0.0524 0.01026 0.00535 0.04396\n", "e1d GENROU 1 0.22346 0.14687 0.02655 0.02655 0.0065 0.02563 0.01032\n", "e1d GENROU 2 0.16265 0.15839 0.02676 0.02676 0.00999 0.01879 0.00058\n", "e1d GENROU 3 0.1893 0.12278 0.02515 0.02515 0.02218 0.00485 0.00021\n", "e1d GENROU 4 0.28887 0.12266 0.02446 0.02446 0.02447 0.0144 0.00922\n", "e2d GENROU 1 1e-05 0.00011 0.00296 0.00296 0.00013 0.0005 0.00144\n", "e2d GENROU 2 9e-05 0.00027 0.00367 0.00367 0.00016 0.00032 0.00127\n", "e2d GENROU 3 8e-05 0.0003 0.00331 0.00331 0.00037 0.0001 0.00124\n", "e2d GENROU 4 4e-05 0.00013 0.0025 0.0025 0.00042 0.00022 0.00131\n", "e2q GENROU 1 0.02215 0.0174 0.00812 0.00812 0.00095 0.00378 0.00189\n", "e2q GENROU 2 0.01605 0.01868 0.00815 0.00815 0.00146 0.00276 0.00011\n", "e2q GENROU 3 0.01847 0.01432 0.00757 0.00757 0.00321 0.0007 4e-05\n", "e2q GENROU 4 0.0281 0.01426 0.00734 0.00734 0.00353 0.00208 0.00165\n", "LAG_y TGOV1 1 0.0003 0.00019 0.0023 0.0023 0.00118 0.00327 0.00534\n", "LAG_y TGOV1 2 0.00103 5e-05 7e-05 7e-05 0.00107 0.00292 0.0193\n", "LAG_y TGOV1 3 0.00077 0.0002 0.00065 0.00065 0.0026 0.00108 0.01706\n", "LAG_y TGOV1 4 0.00037 0 9e-05 9e-05 0.00279 0.00096 0.00309\n", "LL_x TGOV1 1 1e-05 1e-05 0.00065 0.00065 5e-05 0.00015 0.00014\n", "LL_x TGOV1 2 5e-05 0 2e-05 2e-05 5e-05 0.00014 0.00054\n", "LL_x TGOV1 3 3e-05 1e-05 0.00015 0.00015 0.0001 4e-05 0.00037\n", "LL_x TGOV1 4 1e-05 0 2e-05 2e-05 0.0001 4e-05 7e-05\n", "vp EXDC2 1 0.0018 0.02792 0.0755 0.0755 0.08166 0.28549 0.19372\n", "vp EXDC2 2 0.01151 0.07037 0.09369 0.09369 0.10072 0.20041 0.17109\n", "vp EXDC2 3 0.00659 0.05154 0.05539 0.05539 0.15513 0.03978 0.10881\n", "vp EXDC2 4 0.00328 0.02678 0.05123 0.05123 0.20162 0.10006 0.1402\n", "LS_y EXDC2 1 4e-05 0.00026 0.00452 0.00452 0.00039 0.00144 0.00445\n", "LS_y EXDC2 2 0.00022 0.00065 0.00551 0.00551 0.00047 0.00099 0.00386\n", "LS_y EXDC2 3 0.00027 0.00103 0.00707 0.00707 0.00158 0.00043 0.00532\n", "LS_y EXDC2 4 7e-05 0.00026 0.00313 0.00313 0.00099 0.00052 0.00329\n", "LL_x EXDC2 1 0 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 0 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0 0 0 0\n", "LA_y EXDC2 1 0.0001 0.00137 0.0025 0.0025 0.00299 0.0104 0.00439\n", "LA_y EXDC2 2 0.0009 0.00465 0.00418 0.00418 0.00497 0.00984 0.00523\n", "LA_y EXDC2 3 0.00078 0.00513 0.00372 0.00372 0.01154 0.00294 0.00501\n", "LA_y EXDC2 4 0.00032 0.00221 0.00285 0.00285 0.01245 0.00614 0.00536\n", "W_x EXDC2 1 0.00026 0.00528 0.03995 0.03995 0.02435 0.08588 0.11819\n", "W_x EXDC2 2 0.00188 0.01518 0.05653 0.05653 0.03424 0.06873 0.119\n", "W_x EXDC2 3 0.00159 0.01641 0.04933 0.04933 0.07785 0.02014 0.11172\n", "W_x EXDC2 4 0.00065 0.00702 0.03754 0.03754 0.08326 0.04168 0.11845\n", "\n", "PARTICIPATION FACTORS [5/8]\n", " #36 #37 #38 #39 #40 #41 #42\n", "\n", "delta GENROU 1 0.00239 0.01045 0.00051 0.00201 0.00333 0.00065 0.00065\n", "delta GENROU 2 0.00236 0.00093 0.00922 0.00128 0.00055 0.00047 0.00047\n", "delta GENROU 3 0.00232 0.00157 0.00622 0.00373 0.00374 0.00048 0.00048\n", "delta GENROU 4 0.00245 0.0139 0.00281 0.00283 8e-05 0.00065 0.00065\n", "omega GENROU 1 0.05646 0.00087 2e-05 7e-05 0.00453 0.00107 0.00107\n", "omega GENROU 2 0.05461 8e-05 0.00045 4e-05 0.00192 0.00065 0.00065\n", "omega GENROU 3 0.03522 8e-05 0.0002 8e-05 0.00528 0.0005 0.0005\n", "omega GENROU 4 0.03404 0.0007 8e-05 6e-05 0.00043 0.00073 0.00073\n", "e1q GENROU 1 0.00021 0.01108 0.00476 0.00142 0.03388 0.08211 0.08211\n", "e1q GENROU 2 0.00026 0.00532 0.00231 0.00021 0.04408 0.06594 0.06594\n", "e1q GENROU 3 0.0001 0.00154 0.00122 0.0004 0.03867 0.06805 0.06805\n", "e1q GENROU 4 7e-05 0.00458 0.00375 0.0018 0.02724 0.07576 0.07576\n", "e1d GENROU 1 0.00019 0.00414 0.00056 0.00018 0.02258 0.00026 0.00026\n", "e1d GENROU 2 0.0002 0.00398 0.0024 0.00016 0.02473 0.00126 0.00126\n", "e1d GENROU 3 0.0001 0.00324 0.00173 0.00029 0.02321 0.00142 0.00142\n", "e1d GENROU 4 0.0001 0.0042 0.00103 0.00044 0.02112 0.00043 0.00043\n", "e2d GENROU 1 0 0.00015 0.00011 2e-05 0.00025 0.00119 0.00119\n", "e2d GENROU 2 0 6e-05 2e-05 2e-05 0.00035 0.00092 0.00092\n", "e2d GENROU 3 0 2e-05 1e-05 6e-05 0.00034 0.00105 0.00105\n", "e2d GENROU 4 0 2e-05 7e-05 2e-05 0.00024 0.00127 0.00127\n", "e2q GENROU 1 4e-05 0.00084 0.00011 4e-05 0.00566 8e-05 8e-05\n", "e2q GENROU 2 4e-05 0.00081 0.00049 3e-05 0.00617 0.00037 0.00037\n", "e2q GENROU 3 2e-05 0.00065 0.00034 6e-05 0.00573 0.00041 0.00041\n", "e2q GENROU 4 2e-05 0.00084 0.00021 9e-05 0.0052 0.00012 0.00012\n", "LAG_y TGOV1 1 0.28771 0.44443 0.03358 0.18229 0.00614 0.0005 0.0005\n", "LAG_y TGOV1 2 0.27828 0.03862 0.60969 0.11639 0.0026 0.0003 0.0003\n", "LAG_y TGOV1 3 0.23415 0.05545 0.34788 0.28645 0.00933 0.0003 0.0003\n", "LAG_y TGOV1 4 0.24068 0.49682 0.15655 0.21672 0.0008 0.00047 0.00047\n", "LL_x TGOV1 1 0.02949 0.00262 0.00012 0.00046 0.00245 0.00046 0.00046\n", "LL_x TGOV1 2 0.02951 0.00024 0.00224 0.0003 0.00107 0.00029 0.00029\n", "LL_x TGOV1 3 0.01928 0.00026 0.00099 0.00058 0.00299 0.00022 0.00022\n", "LL_x TGOV1 4 0.0199 0.00236 0.00045 0.00044 0.00026 0.00035 0.00035\n", "vp EXDC2 1 0.00028 0.01851 0.01021 0.00225 0.02054 0.04771 0.04771\n", "vp EXDC2 2 0.00034 0.00805 0.00341 0.00124 0.02771 0.03765 0.03765\n", "vp EXDC2 3 0.00019 0.0005 0.00126 0.00202 0.01745 0.02788 0.02788\n", "vp EXDC2 4 0.00016 0.00385 0.00608 0.00268 0.01504 0.0403 0.0403\n", "LS_y EXDC2 1 2e-05 0.00065 0.00035 8e-05 0.00189 0.00251 0.00251\n", "LS_y EXDC2 2 2e-05 0.00028 0.00012 4e-05 0.00251 0.00194 0.00194\n", "LS_y EXDC2 3 2e-05 4e-05 9e-05 0.00015 0.00343 0.00312 0.00312\n", "LS_y EXDC2 4 1e-05 0.00014 0.00021 9e-05 0.00142 0.00216 0.00216\n", "LL_x EXDC2 1 0 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 0 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0 0 0 0\n", "LA_y EXDC2 1 0 0.00028 0.00016 4e-05 1e-05 0.00088 0.00088\n", "LA_y EXDC2 2 0 0.00016 7e-05 3e-05 2e-05 0.00094 0.00094\n", "LA_y EXDC2 3 0 2e-05 4e-05 6e-05 2e-05 0.00105 0.00105\n", "LA_y EXDC2 4 0 0.0001 0.00016 7e-05 1e-05 0.00125 0.00125\n", "W_x EXDC2 1 0.0007 0.01967 0.01051 0.00229 0.2114 0.11396 0.11396\n", "W_x EXDC2 2 0.00099 0.00975 0.004 0.00144 0.32512 0.10252 0.10252\n", "W_x EXDC2 3 0.0008 0.00089 0.00219 0.00345 0.30227 0.11206 0.11206\n", "W_x EXDC2 4 0.00057 0.00566 0.00867 0.00377 0.21431 0.13329 0.13329\n", "\n", "PARTICIPATION FACTORS [6/8]\n", " #43* #44 #45 #46 #47 #48 #49\n", "\n", "delta GENROU 1 0.23884 0.00446 0.00446 5e-05 5e-05 4e-05 4e-05\n", "delta GENROU 2 0.24952 0.00469 0.00469 2e-05 2e-05 3e-05 3e-05\n", "delta GENROU 3 0.26107 0.00465 0.00465 2e-05 2e-05 7e-05 7e-05\n", "delta GENROU 4 0.25057 0.00449 0.00449 7e-05 7e-05 5e-05 5e-05\n", "omega GENROU 1 0 0.10227 0.10227 0.00063 0.00063 0.00032 0.00032\n", "omega GENROU 2 0 0.10619 0.10619 0.00022 0.00022 0.00028 0.00028\n", "omega GENROU 3 0 0.07177 0.07177 0.00013 0.00013 0.00029 0.00029\n", "omega GENROU 4 0 0.06484 0.06484 0.00031 0.00031 0.0002 0.0002\n", "e1q GENROU 1 0 2e-05 2e-05 0.07049 0.07049 0.04381 0.04381\n", "e1q GENROU 2 0 2e-05 2e-05 0.05619 0.05619 0.05594 0.05594\n", "e1q GENROU 3 0 3e-05 3e-05 0.0392 0.0392 0.06379 0.06379\n", "e1q GENROU 4 0 4e-05 4e-05 0.04984 0.04984 0.04973 0.04973\n", "e1d GENROU 1 0 0 0 0.00451 0.00451 0.00282 0.00282\n", "e1d GENROU 2 0 0 0 0.00354 0.00354 0.00354 0.00354\n", "e1d GENROU 3 0 0 0 0.00267 0.00267 0.00435 0.00435\n", "e1d GENROU 4 0 1e-05 1e-05 0.00388 0.00388 0.0039 0.0039\n", "e2d GENROU 1 0 0 0 0.00069 0.00069 0.00042 0.00042\n", "e2d GENROU 2 0 0 0 0.00051 0.00051 0.0005 0.0005\n", "e2d GENROU 3 0 0 0 0.00039 0.00039 0.00063 0.00063\n", "e2d GENROU 4 0 0 0 0.00056 0.00056 0.00055 0.00055\n", "e2q GENROU 1 0 0 0 0.00136 0.00136 0.00085 0.00085\n", "e2q GENROU 2 0 0 0 0.00106 0.00106 0.00106 0.00106\n", "e2q GENROU 3 0 0 0 0.00079 0.00079 0.00129 0.00129\n", "e2q GENROU 4 0 0 0 0.00115 0.00115 0.00116 0.00116\n", "LAG_y TGOV1 1 0 0.04715 0.04715 0.00028 0.00028 0.00015 0.00015\n", "LAG_y TGOV1 2 0 0.04895 0.04895 0.0001 0.0001 0.00013 0.00013\n", "LAG_y TGOV1 3 0 0.04317 0.04317 7e-05 7e-05 0.00017 0.00017\n", "LAG_y TGOV1 4 0 0.04147 0.04147 0.0002 0.0002 0.00012 0.00012\n", "LL_x TGOV1 1 0 0.12629 0.12629 0.00092 0.00092 0.00049 0.00049\n", "LL_x TGOV1 2 0 0.13567 0.13567 0.00033 0.00033 0.00045 0.00045\n", "LL_x TGOV1 3 0 0.0929 0.0929 0.00019 0.00019 0.00047 0.00047\n", "LL_x TGOV1 4 0 0.08958 0.08958 0.00052 0.00052 0.00034 0.00034\n", "vp EXDC2 1 0 8e-05 8e-05 0.02856 0.02856 0.01757 0.01757\n", "vp EXDC2 2 0 7e-05 7e-05 0.02194 0.02194 0.02157 0.02157\n", "vp EXDC2 3 0 6e-05 6e-05 0.01098 0.01098 0.01765 0.01765\n", "vp EXDC2 4 0 8e-05 8e-05 0.01848 0.01848 0.01825 0.01825\n", "LS_y EXDC2 1 0 0 0 0.00105 0.00105 0.00064 0.00064\n", "LS_y EXDC2 2 0 0 0 0.00079 0.00079 0.00077 0.00077\n", "LS_y EXDC2 3 0 0 0 0.00086 0.00086 0.00136 0.00136\n", "LS_y EXDC2 4 0 0 0 0.00069 0.00069 0.00068 0.00068\n", "LL_x EXDC2 1 0 0 0 0 0 0 0\n", "LL_x EXDC2 2 0 0 0 0 0 0 0\n", "LL_x EXDC2 3 0 0 0 0 0 0 0\n", "LL_x EXDC2 4 0 0 0 0 0 0 0\n", "LA_y EXDC2 1 0 0 0 0.00049 0.00049 0.0003 0.0003\n", "LA_y EXDC2 2 0 0 0 0.00051 0.00051 0.0005 0.0005\n", "LA_y EXDC2 3 0 0 0 0.00038 0.00038 0.00061 0.00061\n", "LA_y EXDC2 4 0 0 0 0.00054 0.00054 0.00053 0.00053\n", "W_x EXDC2 1 0 0.00025 0.00025 0.11968 0.11968 0.07515 0.07515\n", "W_x EXDC2 2 0 0.00028 0.00028 0.10481 0.10481 0.10518 0.10518\n", "W_x EXDC2 3 0 0.00032 0.00032 0.07742 0.07742 0.12701 0.12701\n", "W_x EXDC2 4 0 0.00035 0.00035 0.10724 0.10724 0.1081 0.1081\n", "\n", "PARTICIPATION FACTORS [7/8]\n", " #50 #51 #52\n", "\n", "delta GENROU 1 0.00304 0.00244 4e-05\n", "delta GENROU 2 0.00157 0.00344 8e-05\n", "delta GENROU 3 0.00179 6e-05 0.00325\n", "delta GENROU 4 0.00266 0 0.00268\n", "omega GENROU 1 0 0 0\n", "omega GENROU 2 0 0 0\n", "omega GENROU 3 0 0 0\n", "omega GENROU 4 0 0 0\n", "e1q GENROU 1 0 3e-05 0\n", "e1q GENROU 2 2e-05 7e-05 0\n", "e1q GENROU 3 0 0 6e-05\n", "e1q GENROU 4 0 0 3e-05\n", "e1d GENROU 1 1e-05 2e-05 0\n", "e1d GENROU 2 0 5e-05 0\n", "e1d GENROU 3 1e-05 0 4e-05\n", "e1d GENROU 4 1e-05 0 3e-05\n", "e2d GENROU 1 0 0 0\n", "e2d GENROU 2 0 0 0\n", "e2d GENROU 3 0 0 0\n", "e2d GENROU 4 0 0 0\n", "e2q GENROU 1 0 1e-05 0\n", "e2q GENROU 2 0 1e-05 0\n", "e2q GENROU 3 0 0 1e-05\n", "e2q GENROU 4 0 0 1e-05\n", "LAG_y TGOV1 1 0.00021 0.00017 0\n", "LAG_y TGOV1 2 0.00011 0.00023 1e-05\n", "LAG_y TGOV1 3 0.0001 0 0.00019\n", "LAG_y TGOV1 4 0.00015 0 0.00015\n", "LL_x TGOV1 1 0.30765 0.38484 0.0055\n", "LL_x TGOV1 2 0.16372 0.56187 0.01273\n", "LL_x TGOV1 3 0.12342 0.00631 0.34609\n", "LL_x TGOV1 4 0.1836 0.00014 0.28568\n", "vp EXDC2 1 3e-05 8e-05 0\n", "vp EXDC2 2 2e-05 6e-05 0\n", "vp EXDC2 3 1e-05 0 4e-05\n", "vp EXDC2 4 2e-05 0 6e-05\n", "LS_y EXDC2 1 0 0 0\n", "LS_y EXDC2 2 0 0 0\n", "LS_y EXDC2 3 0 0 0\n", "LS_y EXDC2 4 0 0 0\n", "LL_x EXDC2 1 0 0 0\n", "LL_x EXDC2 2 0 0 0\n", "LL_x EXDC2 3 0 0 0\n", "LL_x EXDC2 4 0 0 0\n", "LA_y EXDC2 1 0 0 0\n", "LA_y EXDC2 2 0 0 0\n", "LA_y EXDC2 3 0 0 0\n", "LA_y EXDC2 4 0 0 0\n", "W_x EXDC2 1 9e-05 0.00026 0\n", "W_x EXDC2 2 9e-05 0.00023 1e-05\n", "W_x EXDC2 3 5e-05 0 0.00021\n", "W_x EXDC2 4 7e-05 0 0.00029\n", "\n", "\n" ] } ], "source": [ "with open('kundur_full_eig.txt', 'r') as f:\n", " print(f.read())" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Parameter Sweep and Root Loci Plot" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Parameter sweep allows automatically applying a set of parameters to compute the\n", "eigenvalues. It is useful to plot the root loci of the system.\n", "\n", "In this example, we will study the impact of ``EXDC2.KA`` of the device\n", "`1` on the system stability. Let its `KA` vary evenly between 20 to 200 in 10\n", "steps.\n" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [], "source": [ "ret = ss.EIG.sweep(ss.EXDC2.KA, 1, np.linspace(20, 200, 10))" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Plot the 30th to the 50-th eigenvalues (0-based index) on the s-plane:" ] }, { "cell_type": "code", "execution_count": 7, "metadata": {}, "outputs": [ { "data": { "text/plain": [ "(
,\n", " )" ] }, "execution_count": 7, "metadata": {}, "output_type": "execute_result" }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkQAAAG6CAYAAADpkFmuAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjUuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/YYfK9AAAACXBIWXMAAA9hAAAPYQGoP6dpAABMT0lEQVR4nO3dfXQU9b0/8PcmhIQAyeaBitpI3OADohRWEGl7DwkkaBG00oSn21u9ERcrKtAjCbm3EuO1NyTYi9hiScBYT/3xkERaxSfIKsHq0RZYUCvW4j6QKFKEZLKQmBCS7+8PumM22Tzs7uxMNvN+nZNTdmfmO5+djNl3Z77z/RqEEAJEREREOhahdQFEREREWmMgIiIiIt1jICIiIiLdYyAiIiIi3WMgIiIiIt1jICIiIiLdYyAiIiIi3WMgIiIiIt1jICIiIiLdYyAiIiIi3RtSgSgrK0vrEoiIiCgMDZlAVF1dDavVqnUZREREFIaGRCCSJAkNDQ1al0FERERhyjAUZrsvLy/HwoULkZCQgL4+TltbG9ra2uTXnZ2daGhoQFJSEgwGgxqlEhERUZCEEDh37hyuuOIKREQoc21nmCKtaMhqtSIzM3NA6xYXF6OoqCjEFREREZEa6uvr8d3vfleRtsL+ClF1dTWys7MhSZLfV4iamppw1VVXob6+HnFxcWqUS0REYaq2thbp6emq7rO4uBgFBQWq7jMcuN1upKSkQJIkxMfHK9JmWF8hKi8vh8ViGfD60dHRiI6O7vF+XFwcAxEREfXpyJEj+Oqrr7B8+XLV9hkbG8vvpz4o2d0lbDtV22w2TJ06VesyiIhIR06fPq11CRQiYXuFqKGhATabTX7U3m63AwBKS0thMpmQnZ2tZXlERDTEjBo1Cs3Nzarus7OzU9X96VnYBqLMzEyvztQ2mw3l5eXIy8vTsCoiIhqqfv7zn6OqqkrrMihEwvaWWVfV1dUoLi4GAOTn53OARiIiCgmXy6Xq/iZNmqTq/vQsbK8QdZWdnc1bZEREFFKxsbGq7/Ojjz7C3Xffrfp+9WhIXCEiIiJSQ1RUFM6cOaN1GRQCDEREREQDFBkZyX5EQxQDERER0QC1tbWp+uh9RkaGavvSOwYiIiKiAUpISFB1f/v371d1f3o2JDpVExERqeG+++7TugQKEV4hIiIiGqDY2Fhs2LBB6zIoBBiIiIiI/FRfX6/KfiIjI1XZDzEQERER+e2TTz5RZT8dHR2q7IcYiIiIiPxiNpshSZLWZZDC2KmaiIjID3/7299w4cIFVfZlMplU2Q/xChEREZFfrrvuOtX25XA4VNuX3jEQERER+eH222/HqFGjtC6DFMZARERE5IfY2FicP39e6zJIYQxEREREAXjzzTdDvo8lS5aEfB90CQMRERFRANR49H7Hjh0h3wddwqfMiIiI/FRYWIiWlhatyyAF8QoRERGRnzZv3swpPIYYBiIiIiI/ZWVlqbKfqKgoVfZDDERERER+u/baa1XZT3t7uyr7IQYiIiKigMTExODMmTNal0EKYSAiIiIKQGtrK1566aWQ7mPs2LEhbZ++xUBEREQUoNOnT4e0/VOnToW0ffoWH7snIiIKwPTp0xERwesKQwV/k0RERAEwGo14//33tS6DFMJAREREFIBJkyaFfB+FhYUh3wddwkBEREQUgNjY2JDvo6ioKOT7oEsYiIiIiAKUnZ3NR++HCAYiIiKiAFVXV6OiokLrMkgBDEREREQBMhgM+Oabb0LWfkxMTMjaJm8MRERERAEaP358SENLa2tryNombwxEREREAfrxj3+MCxcuaF0GKYCBiIiIKECxsbHo7OwMWfvx8fEha5u8MRAREREFKVRPmjU1NYWkXeqJgYiIiCgI06dPV2VMIgqtsJ7LTJIkVFZWAgDsdjscDge2bt0Ko9GobWFERKQbBw8ehM1mw3/9139pXQoFIayvEOXn5yMzMxMWiwUlJSVITExETk6O1mUREZGOREZGor29PSRtc+oO9YR1IHI4HKiurpZfp6Wl4dChQxpWREREenP99deHrG1O3aGesL5lVlNT4/X64MGDyMzM7HX9trY2tLW1ya/dbnfIaiMiIn1YsGABxwsaAsL6ClFX1dXVkCQJW7du7XWd4uJixMfHyz8pKSkqVkhEREPV8ePHtS6BghT2gUiSJJSXl8PhcCAnJ6fPDtUFBQVoamqSf+rr69UrlIiIhrR33nlH8TZHjRqleJvkW1jfMgMAo9EIi8UCACgvL0dCQgKcTqfPYBQdHY3o6GiVKyQioqFuxIgRIWn3/PnzIWmXegrbK0SSJCE/Px+SJMnvZWZmQpIkWK1W7QojIiLdiY+PxwcffKB1GRSEsA1EDocDpaWlaGhokN/zhCOOQ0RERGpqbm4Oyaz3I0eOVLxN8i1sA5HZbEZeXh5MJpP83q5du2A2m/t80oyIiEhpWVlZIWm3ubk5JO1ST2Hdh6igoAClpaXya0mS8NZbb2lYERER6dFNN93EiVjDXFgHIqPRiLy8PK3LICIiwvPPP48VK1YgOTlZ61IoAGF7y4yIiGiwefnllxVtj1N3qIeBiIiISAERERGKd6zm1B3qYSAiIiJSwLx58xAbG6t1GRQgBiIiIiIFfPLJJ5wBIYwxEBERESlgwYIFireZlJSkeJvkGwMRERGRAmJjYxW/ZXb27FlF26PeMRAREREppKWlBe+9957WZVAAGIiIiIgU9P777yvW1vDhwxVri/rGQERERKSglpYWxdq6cOGCYm1R38J6pGoiIqLBpLCwUNFAROrhFSIiIiKFfPLJJ9iwYYPWZVAAGIiIiIgUMnHiREXbW7FihaLtUe8YiIiIiAapzZs3a12CbjAQERERKSgxMRFvvfWW1mWQnxiIiIiIFNTQ0KDoo/ekDgYiIiIiBUVGRqKjo0ORtq644gpF2qH+8bF7IiIiBd1+++34/PPPFWnr5MmTirRD/eMVIiIiIgVNnToVn332mdZlkJ8YiIiIiAapqKgorUvQDQYiIiKiEPjHP/4RdBvt7e0KVEIDwUBERESksBUrViAmJkbrMsgPDEREREQK27x5M55//nmtyyA/MBAREREpbNSoUYq0s2DBAkXaof4xEBERESls1qxZiIgI/it29+7dClRDA8FAREREpLApU6ZoXQL5iYGIiIgoBDo7O9HS0qJ1GTRADEREREQhUlVVFdT211xzjUKVUH84dQcREVEITJ8+HRMnTgyqjePHjytUDfWHV4iIiIhC4JNPPkFFRYXWZdAAMRARERGFQGJiYtBtKPGkGg0MjzQREVEILFq0KOg2hBAKVEIDwUBEREQUArGxsZgwYUJQbTAQqYeBiIiIKEQ+/fRTHDp0SOsyaAAYiIiIiELowIEDAW87Y8YMBSuhvoT9Y/elpaUAALvdDgAoKyvTshwiIiKZ0WhEUlJSwNu///77mDNnjoIVUW/COhDl5+ejpKREfr18+XJkZWWhpqZGw6qIiIguue6663jLLEyE7S0zSZJgs9kgSZL83vLly2G1WuFwOLQrjIiI6F9Gjx6Njo4OrcugAQjbQAQAhw4d8go/JpMJALxCUldtbW1wu91eP0RERKHygx/8IKjtp02bplAl1J+wDURGoxGNjY0wm83ye1arFcC3wai74uJixMfHyz8pKSmq1EpERPq1ZMmSgLc9ePCggpVQX8I2EPlSXFyMsrIyGI1Gn8sLCgrQ1NQk/9TX16tbIBER6c6OHTuwfft2rcugfoR1p+qu8vPzsWjRIlgsll7XiY6ORnR0tIpVERERAU6nM6DtDAaDwpVQb4ZEIKqurkZaWlqfYYiIiEgLI0aMwMWLFwPaloFIPWF/y8zTb8gThiRJ4lNmREQ0aOTm5gY80WtnZ6fC1VBvwjoQ2Ww22Gw2mM1mOBwOOBwOlJeXKzLDMBERkVL++c9/al0C9SNsb5lJkoTZs2dDkiTk5+d7LcvLy9OoKiKiwaG5uRkjR47sdflXX32F9vZ2XHXVVV7vd3R04Msvv4TRaERcXJzX+5GRkSGrdyhLTk4OeNvrr79ewUqoL2F7hcjz2L0QoscPEZHevPHGGzh69ChOnDiB5uZm3HnnnXC5XHjkkUdw4sQJAJemNnrqqadw9OhRZGdnIycnB/feey8AwOVy4aWXXkJWVhYWL17sNZ6by+XCXXfdxQEGgxATExPQdn//+98VroR6E7aBiIiIvjVhwgQ8+OCDWLx4Mb7++ms899xzyM3NxcKFC5Gbm4t3330X5eXlKCoqwn333YeSkhIMGzYM//jHP/Duu+8iNzcXl112GVpbWyGEkPuuuFwu5Obm4re//S2ef/55uFwur/26XC5s27ZNg08cXtrb27F161aty6A+MBAREQ0Bqamp2L59OwwGgzwQYEVFBdatW4cHHngAd9xxB9ra2nDNNddACIH8/Hxs374dpaWlmD9/PiwWC9atW4ft27djx44dyM3NxYEDB5Cbm4uKigqkpqYiMzMTubm5cijyhKXMzEwNP3l46OjowMmTJ7Uug/rAQERENET4CkVPPPEEli1bhssuuwx1dXXIz8/HiBEjIISAy+XCunXr8Lvf/Q5LlixBYWEhUlNTkZqaisLCQqSnp8vvedqvqKjwGZY8mpubewx6W19fj+bmZrUOw6AUGxsb0Ha33XabwpVQr4SONTU1CQCiqalJ61KIiBTjdDrFjBkzxOTJk8WUKVPEDTfcIACITZs2CQDiwIEDora2VgAQO3bsEBkZGaK2tlZkZGQIp9MpnE5nj/e68mxbW1vbY991dXVe23jaqqurU+GTD17Nzc0Bbff4448rXMnQEIrvbwYiBiIiGoI8oQWAmDx5sti5c6cwGo1i586d4vvf/76YMWOG2LFjhzAajeLPf/6zEOLbIHXrrbf2CDTdX/cWlga6jh4FEm4YiHwLxfc3b5kREQ0xLpcLBQUFuOGGGzBq1Ci0t7dj06ZN2LNnDzZt2oRvvvkGFy5cwDPPPIM9e/Zg3bp1cLlcSElJQUxMDAwGAyIiLn09eG6TPfTQQ7Db7fJtspkzZ8q3z7p3tO7tlhsBLS0tWpdAvWAgIiIaQlwuF5YuXYrW1lbExcWhoqIC9fX1aG1tBQAIIRATE4Nf/epXMBgMSElJkYNNfX09ampqsHPnTq9JslNTU/Hyyy9j//79Xn2GPGHJM2NA1xqKiopQW1uLoqKiHoFJzz777DO/1ufYT+phICIiGiI8YaijowMxMTHYsWMH/vznP+O1115DTEwM7rvvPmzYsEF+umzbtm349a9/LQebI0eOIDIyEldddZXXoIzApS/mZcuW9bjak5qaimXLlsmv6+vrfV5F6t7RWo/S0tLw6aef+rUNx35SDwMREdEQ4AlDQgjs2rULNTU1SE1NxTPPPIMf/vCH2L59O5KSkrBmzRoYDAa88sormDBhAp555hkAl4LN3XffHXQdiYmJeOGFF7yuIr3wwgt9Tqm0cuVKr0f5V65cGXQdg5HT6cTx48e1LoN6EbZTdxAR0bc+/fRTPPvss0hISMC4ceN6LPc8kl9bW4vRo0f3Oa1HMEaOHNmj7ZSUlD63Wb16NXJzc1FYWIiioiJUVFSEpDatjR071u+xiLpPrUKhw0BERDQE/OhHP+p3ndTUVHmqjsGkayfs2traIdsJ++6770ZVVZVf29TV1YWoGuqOt8yIiEhTeumEnZycjNOnT2tdBvWCgYiIiDS1ceNGr07YGzdu1LqkkPr444+1LoF8MAih3+nh3W434uPj0dTU1OOJCiIiIqUVFRVh7NixWL58+YDWr6+v77cPlh6F4vubV4iIiIhUsmbNmgH19/IYqh3MByMGIiIiIpX89re/xfPPP691GeQDAxEREZFKbrnlFq1LoF4wEBEREakkPT3dr/WHDx8emkKoBwYiIiIiFUVGRg54ktcLFy6EuBryYCAiIiJSUUdHB/bs2aN1GdQNAxEREZHK7Hb7gNYbM2ZMiCshD07dQUREpKKZM2ciPj5+QOt+/fXXIa6GPHiFiIiISEWNjY145ZVXtC6DumEgIiIiUtG//du/aV0C+cBAREREpKLk5OQBr1tYWBjCSqgrBiIiIiKVTZ8+HWfOnOl3vaKiIhWqIYCBiIiISHV/+ctf8OKLL2pdBnXBQERERKQBt9utdQnUBQMRERGRylJTUzF69Oh+1xsxYoQK1RDAQERERKS6O+64A+fOnet3vW+++UaFaghgICIiIlJdcnIyhBBal0FdMBARERFppL9JXgdyW42UwUBERESkgTvvvLPfdQZyW42UMeC5zHbv3o2//vWvMBgMfu9ECAGDwYDi4mK/tyUiIhqKXnnlFURHR2Pt2rVal0LwIxDV1NTgd7/7XcA7WrhwYcDb9kaSJFRWVqKqqgo1NTWKt09ERBQqERERaGtr07oM+pcBByKz2RzUjrKysoLavjubzYZDhw5BkiQ0NDQo2jYREVGoffe730VdXV2f63DqDvUMuA/R/fffH9SOgt2+O7PZDIvFApPJpGi7REREavjP//xPjB07ts91OHWHeoLuVL17924l6iAiItKdU6dOaV0C/cuAbpk99dRTvS6rqanBggULFCsolNra2rzu13LYdCIi0prdbkdaWprWZejegK4QCSFw5swZxMfH+/wJF8XFxV51p6SkaF0SERHp2MiRI9Ha2trrco5DpJ4BXSGyWCywWq34yU9+0mNZOPXhKSgowC9+8Qv5tdvtZigiIiLNRERE4PXXX8fEiRN9Luc4ROoZUCCKj4/3GYYAYPbs2YoWFErR0dGIjo7WugwiIiIAwIULF/jo/SDBkaqJiIg0MmnSpD6Xc7Z79QQciPrqaK0mjkFERKSNyspKWK3WPtexWq2orKxUqaLwM3fuXKxZs6bX5ZztXj0BB6IzZ84oWYffHA4HSktLUVZWBpvNhvz8fFRXV2taExGRniQmJmLOnDm4//774XK5eizfvHkz5syZg08++UT94sLIhg0b+p3klUIv4EAUyJxmSjKZTMjLy8Phw4chhEBJSQmys7M1rYmISE8yMzPxm9/8Btu2bUNGRgbsdru8bPPmzXjooYcwbty4sBmaRUt79+7VugTdG/DUHURERN2tWLECjY2NeOyxx3DjjTfikUcewciRI1FYWIjU1FQkJCTgwQcfxD333AOLxaJ1uYNSREQEJEnyuYxTd6iHnaqJiCgov/zlL1FUVITW1lb8+te/RmFhIcaOHYuEhAQAQGdnJ2688UY0NzdrXOng9P3vfx/t7e0+l3HqDvUwEBERUdDWrVuHOXPmoKOjAwBw+vRpnDt3DtHR0diwYQN++ctf8iGYXnz22Wf46quvtC5D9wIORImJiUrWQUREYWzz5s3Yt2+fPHtBZ2cnPv/8cyxevBjr1q3DuHHj5LDkcrmwcuVKLcsdVObPn691CYQgAlFfjwkSEZF+eDpQp6amwmQy4bLLLgMADBs2DKtWrYLFYkFhYSFyc3Nx4MAB5ObmYvXq1RpXPXikpKQgLi7O57IxY8aoXI1+8ZYZEREFrGsY8vQZuvrqq5Gbm4uLFy9i2LBhWLZsGb744gsUFhYiPT1d7nBN33K73V5P6Xl8/fXXGlSjT0EForfeeguLFi2SX+/evRtvv/120EUREdHgZ7Va8fDDD2PJkiW47LLLYDQasXDhQvz617+G0+nE6tWrcfHiRYwePRo/+9nPUFBQgNraWhQVFfkct0jv3njjDa1L0LWgH7tfv349AGDatGmQJAkWiwVlZWXYunVrr5cAiYgo/DU0NGDfvn245ZZbIEkSrrzySkRGRqK5uRkvvPACUlJScPPNN6OtrQ3vvvsuHn30Udxwww2oqKjAxo0bsWnTJq0/wqDS1NTU473hw4drUIk+GYQQItCNd+/ejQULFsDpdGL8+PE4fPgwJk+eDADYtm0bli1bplSdIeF2uxEfH4+mpiaGNyIi0tSZM2eQnJzs9V5RURHHIvIhFN/fQd0ys9vtcLlcyM/Ph9FolMMQAPleMhEREfXt9ddfx+bNm7UuQ9eCCkQWiwXr16+Hw+GQJ+/z3BfWemoPIiKicJGenq51CboXVB+i+Ph4bNmyxeu9zMxM5OTkYNq0aUEVRkREpBexsbE+38/NzVW5Ev1S/LH7zz//HGvXruVkfkRERH4YMWIE3nvvPa/3KioqNKpGfwYciNxu94Ab9YxUGuj2REREwaqsrITVau1zHavVKnf50No333yDd955R+sydGvAgSg/Pz+oHQW7PRERUXd9hZ7ExETMmTMHP/nJT3yGntWrV2POnDmDZioqg8GACxcuaF2Gbg24D5EQAgUFBQHtRAiBmpqagLYlIiLqjSf0rFy5Ehs3bvRalpmZie9973vYvXs3vvjiCyxcuFBeNn/+fLz66qtYtWoVMjMz1S7bp7lz5+LMmTNe71111VUaVaM/Aw5EW7Zs8Tlo1EAFGqaIiIh609DQgLvvvhtPP/00vvzyS68rQQsXLsTRo0cxZswY/PWvf0VGRgb279+PjIwM1NbWYt68eT1ClJZuuOEGbNiwAbfffrv8Xl1dnYYV6YtfT5n56htERESklcTERPzxj3+EyWRCVVUVFi5ciMrKSixcuBBVVVUwmUxwOp2YNGkSamtrERcXh3PnziE9PR179uzRunwvvT1pRurg5K5ERBS2MjMzsXLlSjgcDkRGRqKqqgqXX345qqqqEBkZCYfDgZUrV+LDDz/E6NGjce7cOYwePRr79+/XuvRetbS0yP+OjIzUsBJ9YSAiIqKwtnHjRsybNw8dHR0AgFOnTgEAOjo65NtiGRkZchg6d+4cMjIytCy5V3feeSe++OIL+bXnM1HoMRAREVHY27NnT48poxISErBnzx7Mnz8ftbW1SE9Ph9vtRnp6OmprazF//nyNqu3dK6+8gl27dmldhi4xEBERUdhbuHAhGhsbvd5rbGzE+PHj8eqrr2LevHnybbL9+/dj3rx5ePXVV7F69Wotyu3V8OHD0dnZqXUZuhRUIPLMW0ZERKSV1atXo6qqSp5Dc8SIEfIyu92OyZMn9+hAvWfPHqxatQqbNm3qd/BGNc2YMQNRUVHy69tuu03DavQlqEC0fPlypeogIiLym9VqxaZNm3DVVVdBCIHp06ejpaUF06dPBwAkJyfjww8/9Bl6Nm7ciH379qGhoUHtsnuVnp6Oixcvyq/37t2rYTX6ElQgstvteOqpp7B7925OzUFERKpraGjAj3/8Y9TV1WHevHn44IMPAAAffPAB5s2bhzNnzuDHP/5xr6EnMzPTa8DGwUAIoXUJuhRUIKqpqcGjjz6KBQsW4ODBg9i6dSt2796tVG1ERER9SkxMxJ/+9CesWrWq19tif/rTnwbN9BwD9frrr2tdgu4EFYiuvvpq+d+zZ8/G/fffDyEEkpKSsGjRIrz99ttBF0hERNSbhoYG7Nu3r9cRpwfjbbH+3HbbbUhPTwcAXH/99doWoyN+jVTd3dtvv41Zs2bh6NGj2LVrF8rLy5GYmIi1a9ciOzsbiYmJ2Lp1K5KSkrBgwQKlaiYiIgKAAd3uGixzlQ3U/v37sXfvXhQWFuLvf/+71uXoRlCBKDs7G0lJSWhoaIDFYoHVasWUKVO81rn//vsBANu2bcOyZcuC2R0REdGQN2rUqLC6ojVUBBWITCYTSkpKMHv27D7X27ZtWzC7ISIi0o25c+fixRdfBABERHC4QLUEdaQLCgr6DUNNTU04dOgQpk6dGsyuiIiIdCEtLQ2FhYUA+MSZmoIKRAcPHkRBQUGf68THx2PLli2YPHlyMLsiIiLSHQYi9QR9LW7RokU+3+co1kRERBQuggpEWVlZvXb8qq6uDqZpIiIiXRJCoMZRg/858D+IvDIS/zz/T61L0oWgOlWXlZXB4XDgyJEjMJlMMBqNAC79Mp1OJx599FElauxTaWmpvF9JkpCXlxfyfYajN954A5dffjkSEhKQnJyMkSNH9ljH5XKhtrYWs2fPRkpKigZVDh6rV6/G2LFjkZ+fL79XWVmJxMRE+RHekpISnDp1ymv8E6vVioaGhkE38q0/un9OX8Lxc3r+G4iOjsaECRN8ruNyubB3717k5OSE3UB+pD5ffye68/V3oi87/7YTP3/15zh34RwA4Jfil0jZmILvjf0e3v7Z2xgdPVqR2qmnoK4QORwOFBQUYN++fdiyZQvWr18v/9x8881K1dir0tJSAIDFYoHFYoHZbOb8ar2YMGECHnzwQSxcuBDjx4/HunXr8N5776GjowPApS+CpUuX4qmnnsLkyZN1/2Tg2LFjsXbtWq9hJBITEzFnzhysXr0aU6ZMwdq1azF27Fh5+erVqzFnzpyw/yLt+jl9CdfPOWHCBDzwwAO46aab8MQTT/RY7nK5MGvWLKxYsQL//Cf/Hzn1z9ffia58/Z3oy4sfvYif/fFnONd2DhGIwPCI4QAuXWQ4cvIIrvvtdfim/RvF6qduRBCsVmtAy5RiNBpFY2Oj13v+fKSmpiYBQDQ1NSlc2eDkdDrFjBkzxMSJE8XIkSPFqFGjxK233iry8vKE2WwWEydOFLGxsSImJkaYzWZRVlYmhBCirq5OnD9/XuPq1Td58mQBQCQnJ8vvrVixQgAQAMTkyZOFEJeOz9133y0AiFWrVmlVrqJWrVolAIh58+Z5vT9v3ryw/pxOp1NcffXVAoBYtGiR1/vf+c53BABRVFSkYYUUbnz9nRBCiOTkZK+/E/252HFRjP7f0SLy8UgR8z8xIv5/48UVv75CFD5eKEb+aqQY/sRwEVkUKR7d+2goPkbYCcX3d1CByBen0yleeukl8dJLLyndtBe73e4z/AAQNTU1A2pDb4FICO9Q5Ak/I0aM8Prf2NhYMWXKFOF0OoXT6RQZGRmirq7Or/088sgjwul0yvt85JFHQvBpQs/zRy0qKko4nU4xbNgwORDNnTtX3HvvvSIpKUkAEDk5OVqXqyhP+ElPTxdCCJGenu4zJIWbruFn7ty5wul0yr/DriGJaKA8fyfi4uKEEELExcX5DEl92fm3nSKyKFIMf2K4MBYbxX1/uk80tTaJP+z+gzBtMonYJ2NF5OORwlhsDNXHCCuh+P4Oqg8RALjdbrk/QVc1NTUhna7D4XD4fN9oNEKSJJ/L2tra0NbWJr92u92hKG1QS01Nxfbt27F06VJceeWVsNvt6OzshMFgAHBp/IvRo0cjOjoaJ06cQFFRESoqKvzuU7R69Wrk5uaisLBQbiMcff3110hMTERjY6PX3H2jRo3C66+/jmHDhuHixYvIyclBZWWlhpUqb8+ePcjIyEBtbS3i4uJw7tw5pKen95hAM9ykpqbiL3/5C6ZOnYrXX38d11xzDS5evIi5c+di586dWpdHYejrr79GfHw83G63/Lc0Li4OX3/99YDb2O/cDwCINERidPRobJm/BcMihsH+kR0/z/o5njzwJNoj2tFysQUXOy9iWETQX9/UjUGIwAc5OHLkCHJycmAymdDQ0ACTyQRJktDY2IiqqiqkpqYqWKo3q9WKrKysHmM0pKWlIT8/HxaLpcc2jz/+OIqKinq8v3btWsTExOC2227D3r175fevv/56eR6ZiIgIiEtX1AAA06ZNw8GDB+V1Z8yYgffffx8AYDAYYDAY0NnZCQC45pprcPz4cXndBQsWYPfu3fLrqKgotLe3AwCuuOIKnDx5Ul62YsUKbN68WX49fPhwXLhwAQCQlJSEs2fPyss84cNj5MiRaG5uBnDpC/z8+fPysvT0dOzff+k/QCEEzp49i+TkZADAsGHD0NHRASEEDAZDj3bHjh2LU6dO9agdAJYsWYIdO3bIrz/77DNcd911AIDIyEi5zxIAZGRkyDUAwKRJk/DRRx/JxxuAfAy7LvO1bde2TSaTV2DuXlPXmrt+Fl/HMCYmBq2trRBC4MyZM0hKSpLPhaKiInnwNIPB4HWMux57X+12/d11/Z0CPX/nXc+J7se7+7nU9Vzrfs52PUcB73O4+znb9dwXQuDll1/GXXfdJf+x73q8r7rqKtTV1cnt5ubmegXgrp9vzJgxXl8S3Y/LiBEj8M03l/pIjB49GufOnet13a7Lu27na92u++1ajxACmzdvxooVK+TP1vXzdD9n9fQ3ovu68fHxaGpqAvDtfxe9revP34iu/72G898I4FIwSkpKkn+/3dft629Ewy0NSPhrAgBAQCDBmIAm6dLxvhhxERGdETDg0jm64sEVePbZZ+Vttf4bAfT8byPUfyNaW1uxfv16NDU1IS4uDooI5vLSAw88IP/bZrMJSZLk12+99VYwTferpqbG5y0zo9Eo933prrW1VTQ1Nck/9fX1urtlJsS3t82uueYaERERIQAIg8EgDAaDGD9+vJgyZYq49dZbRW1trcjIyJBvffm7j4yMjKDaGCwSEhLk22Sen1GjRgkA8i20oXa7zMNzm2z06NFet8/CXdfbZJ7f4dy5c7Uui8KY5zaZ58dz+2ygbCdtYtgTw0RUUZSI/9948b3ffU/8vw//n3js8cfE5RsuFyOeHCEiH48U4zaOC80HCDOhuGUW1FNmXR/LNZlMqKqqCqY5v5hMJp/vS5LU67Lo6GjExcV5/eiN52kyt9uNL7/8EsOHD8eIESMQHR2NmJgYfPHFF/jss8/Q1taGcePGoaKiArm5uaivr/drPxs3bkRFRQVmzpyJioqKAT9yOtiMGTMGjY2NiIqKgtPpxLBhly5Tnz9/HnPnzsVPf/pTJCUloaqqKqweQR+I+fPno7a2Funp6XC73UhPT0dtbS3mz5+vdWlBcblcmD59Os6ePYu5c+fi+PHjSEpKwuuvv47FixdrXR6FoTFjxsDtdiMuLg5CCMTFxcHtdmPMmDEDbmPK5VNw+ajL0SE60HKxBcfPHsfDbz6MTnSiqa0J7R3tgAH4VcavQvhJdC6YNFVaWiqcTqfYunWrEEKInJwccfToUSGEEGvXrg06rfXHaDQKu93u9Z4/H0lvnar5lJl/+JQZnzIj6o9ST5kJIYSjwXHpSbOiSBHxeIQY/sRwUfh4oYh4PEJEFkWKH+/4sdLlh61B95SZJEkiPz9fLFy4UAghxOHDh0VCQoKIiIhQJRCVlJR43R6rqqoSFotlwNvrKRB5wtAtt9wixo4dKx577DHx7rvviosXL3otnzhxokhMTJRDrl6tX7++xx+zmpoaYTAYxKpVq+Q/guvXr5eXr1q1ShgMhgE/5ThYdf2cvoTr53Q6nWL69OkiMjLSZ+jxhKXIyEhx7NgxDSqkcOPr70RXvv5O9Oerc1+J+dvnixFPjhBRT0SJwscLxdgNY8Wzf31WqbKHhFB8fwfVqdqXpqYmOByOXgeqUlppaal8i+zgwYMoKSkZ8LZut1vuKDjUb59xpGr/cKRqjlQdbgNPkvpCMVK1x8XOizjbchZv7n4T9/zsnmBLHXJC8f2teCDy2L17d0gfu1eCngIRERGFn65Ps9K3QvH9HfRABkePHoXD4egxDlFZWdmgD0REREREQJCBaO3ataiurvaa2BW49KRXbwMnEhEREQ02QQWipKQkfP755z6XbdiwIZimiYiIdOmdd97B/v37UVhYyDstKgpqHKLexvsBgDVr1gTTNBERkS5dffXVcsf/rqNMU2gFFYjS0tLw9ttv+1z21FNPBdM0ERGRLv3+97+H0+nUugzdCeqW2bJlyyBJEpxOp9fVIiEEnE4nHn300aALJCIi6s0bb7yBCRMm9Dl3psvlwqeffoof/ehH6hUWhM7OTkRFRWldhu4E/ZRZWVlZj/E6hBBYv359sE0TERH16fLLL8fixYuxc+dOn6HI5XJh8eLF2LJli/rFBSErKwsAGIxUFFQgKikpwezZs30uKygoCKZpIiKifkVHR+PQoUOYNWsW3n77ba9Q5HK5MGvWLNTV1SE6Olq7Iv00Y8YM3HTTTQDgNXM9hVZQfYh6C0MAVBupmoiI9Ovjjz/G0qVL4XQ6MWHCBLz77rsAgHfffRfXXnstnE4nli5dio8//tjn9larFZWVlWqW3KczZ87g/fff17oMXQoqEPWFnaqJiCjUEhMT8eKLLyI5ORmtra3IyMjApk2bkJGRIV9d+cMf/uBzKpbVq1djzpw5g2qalpdfftnr9RVXXKFRJfrjdyBatGiR/Bjg+PHjcc011/T4GT9+fJ9zuxARESkhMzMTK1euxJkzZwAAFy9exKpVq3Dx4kWv9TZt2uT1ev78+Xj66aexcuXKPuftU9vo0aMxYsQI+fXJkyc1rEZf/O5DdPXVV8ujUptMJp+TqbJTNRERqWXjxo1wu92oqKjosSw3NxenT5/Gq6++ioyMDOzfvx8ZGRmora3FvHnz/J50NdQcDgeuvfZarcvQJb8DUdegU1paismTJ/tcj52qiYhILY899hh+//vfo7OzU34vIiICjz32GFJTU+UQFBcXh3PnziE9PR179uzRsGLf2tra0NLSonUZuhRUH6Li4uJel02ZMgXbtm3Dz3/+c460SUREIeNyuTB9+nSvMARcGs9n+vTpcLlc2L9/P0aPHo1z585h9OjR2L9/v0bV9q/rdB0rVqzQsBJ9CSoQWSwWAJdOxu62bt2KqqoqWCwW2O12hiIiIlKc59H606dPy+9dffXV8r9Pnz6NWbNm4dZbb5XD0Llz55CRkaFFuf3KyMhAbGys/Hrz5s0aVqMvQQUim82G8ePHw2QyITIyEtu2bZOXlZeXo6SkBFOmTMGaNWtgt9uDLpaIiKirvXv3oq6uTn4aKzc3Fw6HA7m5uQAuPaXldDrxl7/8Benp6XC73UhPT0dtbS3mz5+vZek9vPfee4P6ytVQF1QgstvtOHz4MDo7O9HR0YHGxkb5alFjY6PXdB59TQRLREQUiJycHPz7v/87Tp48idzcXDz33HMAgOeeew65ubnyU1pZWVly2Ni/fz/mzZuHV199FatXr9as9u44/pC2ggpEU6dORXx8vPzaYrHAZrMBuBSI4uLi5GUGgyGYXREREfVgs9nwhz/8AatWrZLDkMeSJUvkf+fl5Xkt27NnD1atWoVNmzbBarWqUmt/2tvbERkZ6fXe8OHDNapGf4KauuPQoUNITEyE2WyGJEkoLi7G1q1bAVwKRJ77tcClRwmJiIiU1NDQgH379vkcS6ihoQE1NTXyv7vbuHEj7rjjDp/LtJCSkoKZM2d6vXfhwgWNqtGfoOcyu//++5GdnY2bb74ZFosFu3btgs1mQ1VVFZYtW4bly5ejqqoKOTk5StVMREQEAFi4cGFAyzwG06CMdrsdEydOREpKital6FJQgSg+Pr7POWAyMzNRXl6O7OxszJo1K5hdERERDXnd5wFNSkrSqBL9CSoQ+eJyueR+RAsWLMCaNWuU3gUREdGQExUV1eO9s2fPalCJPgUdiNxuN6xWa497sDU1NV6DSxEREZFvu3fvRkdHh9Zl6FpQgejIkSPIycmByWRCQ0MDTCYTJElCY2MjqqqqlKqRiIhoSDt27FiPkbZJXUEFovLycnz++ecALoUjk8kkP4b/9ttvIzU1NegCiYhocKqsrERiYmKfHZM9dxAG0sFZzzo6OjBsWM+v5MLCQg2q0aegxiHq+h+ByWTiVSEiIh1JTEzEnDlzcO+99+K9996Tb/k0Nzejvr4emzdvxpw5c3D48GHk5ubi2LFjAC71NV25cqWWpQ86a9as8TlIZFFRkQbV6FNQgcjhcMDlcmHbtm2Ij4/Hvn378OGHHwKAPPYDERENTZmZmfjNb36DF154AZmZmfjhD3+I9evX48MPP8SsWbPw0EMPISoqCq+99hpqa2uxbNkyHDhwALm5uYNqhOjBYMOGDVqXoHtB3TKzWCwoLi6G0+nEsmXLsHbtWmRkZKCpqanHqKBERDT0eGZjf+ihh2Cz2eB2u7FlyxacOHECAHD55Zejrq4Or7/+Ojo6OuR5xNiloqeuk7qS+oIeh2j9+vXya7PZDKfTCYfD0WMsBSIiGpq6hqLPPvtMvnV25ZVX4sSJE9ixYwe++93vIjc3F7W1tSgqKkJFRQVD0b+cOXOm12UjR45UsRJ9C+qWmS/x8fGYMmUK3G630k0TEdEgtWLFCixYsMDr0fEvv/wSTz/9NMrLy+UQNHPmTFRUVGDjxo0aVju4/PWvf0VaWprPZc3NzSpXo1+KByKP4uLiUDVNRESDzObNm7F7926vyUnHjRuHnTt34oknnsCJEyfkZampqdi0aZNWpQ46Bw8ehN1u17oM3QsqEB05cgRTp05FUlKS109iYiJKS0uVqpGIiAaxzZs346GHHsLw4cNx3XXX4aqrrgIAfPXVV3C73VizZg2efPJJJCYmalzp4DVq1Ci/3iflBdWHqLi4GAUFBTCbzV7vCyGwdu3aoAojIqLBzxOGYmJiMHnyZNx1111YvHgxXnvtNTz00EP4/PPPceONN+Jvf/sbvv/972td7qCUnJyMrKwsn8vOnz+vcjX6FVQgWrRoEX7yk5/4XFZQUBBM00RENMhZrVY8/PDDuOeee3D//ffj1ltvlW+LeTpaP/zww8jMzITFYtGy1EHtzJkzuPbaa7UuQ/eCCkTd5y/rqrGxMZimiYhokGtoaMC+fft6Hal6xYoVuO666/r8rtC7lpYWrUugfzEIIUQwDXiuBHXvIV9VVYW9e/cG03S/JElCZWUlqqqqAhoI0u12Iz4+Hk1NTYiLiwtBhURERL2rrKzEp59+yik6/BSK7++grhA98MADsFqtMJlMXj3kJUnC4cOHgy6uLzabDYcOHYIkSfx/H0REFJbmzZvX5zxvRUVFDEsqCSoQ3XzzzdiyZYvPZaEehtxsNsNsNqO6ujqk+yEiIgqVDRs2YNiwYfjv//5vrUvRvaACUV+PUK5ZsyaYpkOira0NbW1t8msOHklERFq74oortC6BEOQ4RGlpaXj77bd9LnvqqaeCaTokiouLER8fL/+kpKRoXRIREencokWLel0WHx+vYiX6FtQVomXLlkGSJDidTphMJvl9IQScTiceffTRoAtUUkFBAX7xi1/Ir91uN0MRERFpJioqqs9JXZuamlSsRt+CCkQAUFZW1uPWmRDCa9LX/lRXV2PXrl39rudrEEh/REdHIzo6OuDtiYiIlHLo0CG0t7drXQb9S1CBqKSkBLNnz/a5zJ+BGbOzs5GdnR1MKURERGHlwIED/a4TExOjQiUE+NmHyOVyeb3uLQwBPcclIiIiom/NnDkTkyZN6nOd1tZWlaohvwKRP4+4l5eX+11MIDgGERERhaPXXnsNkiRpXQb9i1+3zLZs2eI1AGNfrFZrSDtVOxwOue+RzWZDfn4+pk2bxltvREQUNjiH2eDh19QdiYmJMJlMfY4/BFy6auN0OnH27NmgCwwlTt1BRERa4SjUgdN86o6CgoIBD7gY6pGqiYiIwlVLS0uvk+J2xdCkHr/6EA3klxfIukRERHry0ksvwWq1al0GdeFXIJoyZUpI1iUiItKTuro6rUugboIemJGIiIj8ExcXh4iI/q9JjB07VoVqCAhyLjMiIiLyX1NTE9LT0/td79SpU6EvhgAwEBEREamuo6MDl112mdZlUBcMRERERBpITk7ud52oqCgVKiGAgYiIiEh1M2bMGNB6nPxVPQxEREREKvrNb36D999/X+syqBsGIiIiIhVxDs7BiYGIiIhIZQOdbmLJkiUhroQ8GIiIiIhUZDKZsHr16gGtu2PHjhBXQx4MRERERCpyOBxal0A+MBARERGppKWlResSqBcMRERERCo5e/asX2MLmUymEFZDXTEQERERqeTll1/G8OHDB7w+b6+ph5O7EhERqeTs2bMDmtSV1MffChERkYr8mcE+MjIyhJVQV7xCREREpJIVK1YgNjZ2wOt3dHSEsBrqileIiIhIUytXroTL5QIAuFwurFy5UtuCQmjz5s1+BSJSDwMRERFpavXq1cjNzcWBAweQm5s74EELw83HH3+sdQnUB94yIyIiTaWmpqKwsBDp6emora1Famqq1iWFxPHjx/3eJiMjIwSVkC+8QkRENAS88cYbOHr0KE6cOAEAaG5u7rGOy+XC73//ezQ0NPhcrhWXy4WioiLU1taiqKhIvn021LS2tmLGjBl+bbN///4QVUPdMRAREQ0BEyZMwIMPPojFixfj2LFjuPPOO72ChcvlwtKlS1FWViYvD0Uoam5uRn19vdd79fX1fe5r48aNqKiowMyZM1FRUYGNGzcqXtdgcPz4cRw9elTrMqgXDERERENAamoqtm/fDoPBgPvuuw9FRUXIzc3Fvffei3fffRdLly7F2bNn8Ytf/ALr1q3D008/jQ8//BCPPPIIgEuB6Y9//CM6OjpQV1cHt9vt1f5An3ZqaGjAPffc49VJ+p577kFDQ0Ov22zatEm+TZaamopNmzb5fwDCxC233KJ1CdQLBiIioiGiayhas2YNnnjiCXz00Uf40Y9+hNbWVjz55JOwWCx44okncO7cOdxxxx1YuHAhXC4XcnNzMWnSJGRlZWHx4sWQJElu1+Vy4a677kJ5eXmP21kulwvbtm2TX6ekpKCiosKrk3RFRQVSUlJUOgqD18iRI5Genu7XNpMmTQpNMdQDO1UTEQ0hnlC0dOlSPPzww7hw4QIAoK2tDRs3bsSePXuQl5cHg8GAPXv2YM2aNTAYDNi+fTuAS/1chBDo7OwEADksVVRUAID879TU1B7Lutagh07S/mhpaQnoFuVHH32Eu+++OwQVUXe8QkRENMSkpqaiuLgYR48exbFjx/Dkk0/i2LFjaG1tRWdnJwwGA4QQ6OjokP994sQJ5ObmYvv27dixY0ePKzypqalITU31efWne+DRSydpf+zZs0frEqgfvEJERDTEuFwuFBQUYPLkybhw4QIee+wx7NixA6WlpZg5cyYOHDgAIYR8BQdAj6s5vV3h6e/qT319vVdQ8gSoF154Qde3zQINhZz3TEVCx5qamgQA0dTUpHUpRESKcDqdYsaMGeLWW28Vf/7zn8Wtt94qzGazmDJlivxjNpvFjBkzRG1trbxubW2tyMjIEE6nUzidTpGRkeH1Xtf2e1smhBDnz58XdXV1Xu/V1dWJ8+fPq/DpB6/m5uaAtnviiScUrmRoCMX3N6MnEdEQ4Xm0XgiBDRs2YN26ddixYwc2bdoEu90OACgoKMDnn3+OtrY2REZGQggBg8GAcePGoaKiAkuXLsWSJUu8HoPPzc2Fy+Xy6jPUfZnHyJEje1wJSklJwciRI9U8FIPOhg0b8Nxzz/m9nacvF4UeAxER0RDQNQw999xzKCwsREVFBQwGAwoLC/Haa6/Jj+Q/99xzMBgMmDdvHjZs2IDt27cjNzcXHR0diImJgcFgkG/VeG57PfTQQ9i3b59XnyHPMqvVquEnDx933XWX1iVQH9iHiIhoCPj000/x7LPPIiEhAePGjcMrr7yCkSNH4pFHHkFFRQXGjRsHi8WCc+fOYdasWfjBD34Ap9OJnTt34plnnkFFRQWOHDmCmpoafPnllzAajXLbqampePnllxEZGdljv6mpqVi2bJmKnzR8JScn+70NH7tXj0EIIbQuQitutxvx8fFoampCXFyc1uUQEdEQZLfb8dJLLyEvL8/vbYuKilBYWBiCqsJbKL6/w/oKUWlpKQDI98bLysq0LIeIiKiHPXv24JtvvtG6DOpH2Aai/Px8lJSUyK+XL1+OrKws1NTUaFgVERGRt+7ToNDgFJaBSJIk2Gw2SJIk3+devnw5br75ZjgcDphMJm0LJCIi+pepU6di2LDAvm4zMjIUroZ6E7ZPmR06dAgOh0N+7QlBXeff6a6trQ1ut9vrh4iIKJQOHjwY8PQl+/fvV7YY6lVYBiKj0YjGxkaYzWb5Pc9jn31dHSouLkZ8fLz8o+dRU4mISD3XXnut1iVQP8IyEPlSXFyMsrIyr0dFuysoKEBTU5P8U19fr16BRERENGgNij5E1dXV2LVrV7/rFRQUeF0V8sjPz8eiRYtgsVj63D46OhrR0dEB10lEROSP+vp6/PSnPw14e19jP1FoDIpAlJ2djezs7IC2ra6uRlpaWr9hiIiISG1/+MMf0N7eHvBYQh0dHQpXRL0J61tmnn5DnjAkSZJXR2siIiIttbe3c8b6MBG2vyWbzQabzQaz2QyHwwGHw4Hy8nIkJiZqXRoREREAYPjw4bjuuusC3p7DyKhnUNwy85ckSZg9ezYkSUJ+fr7XskCGRiciIgqFixcvYuHChQFvz7se6gnLQOR57J6IiGgw6+zs1LoEGqCwvWVGREQ0mHnm2aTwwEBEREQUApdffjmmTZsWVBtLlixRqBrqDwMRERFRCPzf//0fbDZbUG3s2LFDoWqoPwxEREREIdDR0YGoqCity6ABYiAiIiIKke9973tal0ADxEBEREQUAtdddx1uv/32oNrgFSb1MBAREREprKWlBZ999lnQ7bS3tytQDQ0EAxEREZHC9u7dq3UJ5CcGIiIiIoWNHDkSsbGxQbczduxYBaqhgWAgIiIiUpjNZkNSUlLQ7Zw6dUqBamggGIiIiIgU1tbWhtbWVq3LID8wEBEREYXAvHnztC6B/GAQQgiti9CK2+1GfHw8mpqaEBcXp3U5REQ0RNjtdqSlpWldxpAViu9vXiEiIiJSkN1ux4svvqhIW0VFRYq0Q/1jICIiIlLQG2+8oXUJFAAGIiIiIgVduHABERH8eg03w7QugIiIaCiJj4/HbbfdpkhbMTExirRD/WOEJSIiUtAXX3wBh8OhSFt8dF89DEREREQKmz17ttYlkJ8YiIiIiBQ0YsQIRabtAC7dfiN1MBAREREp5M0338Q333yjWHtNTU2KtUV9YyAiIiJSyJEjR7QugQLEQERERKQQPnIfvjh1B6fuICIihbS0tCjWf4h6x6k7iIiIBrENGzagvr5esfY4dYd6GIiIiIgUlJKSonUJFAAGIiIiIgW0tLRoXQIFgYGIiIhIAe+88w6+853vKNrmqFGjFG2Pese5zIiIiBRw8OBBdHZ2Ktrm+fPnFW2PescrRERERAro7OzE8OHDtS6DAsQrRERERAowmUyYPHmyom2OHDlS0faod7xCREREpACHw4GbbrpJ0Tabm5sVbY96x0BEREREusdAREREFKQPPvhA6xIoSGE7dYckSaisrAQA2O12OBwObN26FUajccBtcOoOIiJSQl1dHWJjY5GcnKx1KbrAqTu6yM/PR2ZmJiwWC0pKSpCYmIicnBytyyIiIh16/vnn8eyzzyreLqfuUE/YBiKHw4Hq6mr5dVpaGg4dOqRhRUREpGcJCQlal0BBCNvH7mtqarxeHzx4EJmZmX1u09bWhra2Nvm12+0OSW1ERKQvBoMBc+fO1boMCkLYXiHqqrq6GpIkYevWrX2uV1xcjPj4ePmHE/AREZESDAYD0tLSFG83KSlJ8TbJt7AORJIkoby8HA6HAzk5Of12qC4oKEBTU5P8U19fr06hREQ0ZNXX1ys+ZYfH2bNnQ9Iu9TQobplVV1dj165d/a5XUFAAs9ksvzYajbBYLACA8vJyJCQkwOl09hqMoqOjER0drUjNREREAPDqq69qXQIpYFAEouzsbGRnZw94fUmSUFxcjIKCAjn8ZGZmQpIkWK1Wv9oiIiIKxh133IFjx46FpG3OjaaesLxl5nA4UFpaioaGBvk9SZIAwK9xiIiIiIL1/PPP4/jx4yFp+8KFCyFpl3oKy0BkNpuRl5cHk8kkv7dr1y6YzeZ+nzQjIiJS2pVXXql1CRSkQXHLLBAFBQUoLS2VX0uShLfeekvDioiISK9uv/12rUugIIVtIDIajcjLy9O6DCIi0rlJkyYhNjY2JG2vWLEiJO1ST2F5y4yIiGgwePPNN/HRRx+FrP3NmzeHrG3yxkBEREQUoA8//FDrEkghYXvLjIiISGtjxoxBS0uL1mWQAniFiIiIKEBffvklbrzxxpC1f8UVV4SsbfLGQERERBSgzs7OkIaWkydPhqxt8sZAREREFIRrr71W6xJIAQxEREREAQr1VFFRUVEhbZ++xU7VREREAdi6dStOnjyJiRMnhmwf7e3tIWubvPEKERERUQBOnTqldQmkIAYiIiKiAI0aNUrrEkghvGVGREQUgO985zv4j//4j5DuY8GCBSFtn77FK0REREQBUOOW2e7du0O+D7qEgYiIiChAoZrUldTHQEREROSnlpYWxMTEaF0GKYiBiIiIyE8vvPACOjo6Qr6fa665JuT7oEvYqZqIiMhPp0+fVmU/x48fV2U/xCtEREREAYmLi9O6BFIQrxARERH5qbCwEC0tLSHfT0QEr1uohUeaiIjIT0VFRao8YSaECPk+6BIGIiIiIj+cOXNGtX0xEKmHgYiIiMgPx44d07oECgEGIiIiIj/U1dXh+uuvV2VfM2bMUGU/xEBERETkF7vdDrvdrsq+3n//fVX2QwxEREREfjOZTFqXQApjICIiIvJDYmIi7rzzTq3LIIUxEBEREfmhoaFBtUldp02bpsp+iIGIiIhowN555x1V93fw4EFV96dnDEREREQD9OGHH2pdAoUIp+4gIiIaoLvuugvJycmq7c9gMKi2L73jFSIiIqIBev7557FlyxbV9sdApB4GIiIiIj9kZGSotq/Ozk7V9qV3DERERER+mDJlitYlUAgwEBEREQ3AmTNnEBUVpeo+1ZoihBiIiIiIBuTll19Ge3u7qvv8+9//rur+9GzIBKKsrCytSyAioiHs5MmTWpdAITQkHruvrq6G1WrVugwiIhrCsrKy0NraqnUZFCJhf4VIkiQ0NDRoXQYREQ1xe/fuRUdHh6r7vO2221Tdn56FfSCqrKzEwoULtS6DiIh0YMaMGarub+/evaruT8/C+paZ1WpFZmbmgNdva2tDW1ub/LqpqQkA4Ha7Fa+NiIiGltbWVly8eFHV74zW1lZ+R/ngOSZCCMXaNAglW1NZdXU1srOzIUkSEhIS+j0wjz/+OIqKilSqjoiIiELJbrfDZDIp0tagCETV1dXYtWtXv+sVFBTAbDYDAMrLy2GxWABgwIGo+xUiSZIwbtw41NXVIT4+PohPMLS43W6kpKSgvr4ecXFxWpczqPDY9I7Hpnc8Nr3jsekdj03vmpqacNVVV6GxsRFGo1GRNgfFLbPs7GxkZ2cPeH2bzYapU6f6vZ/o6GhER0f3eD8+Pp4nmw9xcXE8Lr3gsekdj03veGx6x2PTOx6b3kVEKNcVelAEIn81NDTAZrPJj9rb7XYAQGlpKUwmk1/hioiIiCgsA1FmZqZXZ2qbzYby8nLk5eVpWBURERGFq7B/7L66uhrFxcUAgPz8fL8GaIyOjkZhYaHP22h6xuPSOx6b3vHY9I7Hpnc8Nr3jseldKI7NoOhUTURERKSlsL9CRERERBQsBiIiIiLSPQYiIiIi0j0GIiIiItK9sHzsPhiSJKGyshLApfGLHA4Htm7d2u9Il6WlpfI6kiQN2Uf8S0tLAXw7tlNZWVmf61utVpSVlSErKwsmkwk1NTWYNm3akBwLyt9j49lmqJ83nv+mqqqqUFNT0+/6ejpn/D02gD7OGcD/zzkUz5tAftc8P3xT5PwQOmOxWITdbvd6nZmZ2ec2JSUloqSkRH5dU1MjLBZLyGrUSl5entfrgRybqqoqYTQaBQBhMplEWVlZKEvUTCDHRg/nzeHDh0VZWZkoKSkRZrN5QNvo5ZwJ5Njo4ZwRIrDPOdTOm0COAc+P3ilxfuguEGVmZnod6JKSEmE0Gvvcxmg0isbGRq/3hlqWbGxsFJmZmV6f8/DhwwKAV4DsrqqqqsexGWoCPTZ6OG88qqqq/ApEQ/2c6cqfY6OXcyaQzznUzptAjgHPj94pcX7org9RTU2N16W3gwcPeo163Z3D4YAkST5vqfkzCGQ4OHToEBwOh/zaM4OwJEkaVTR4+Hts9HTekDL0cs7o5XP2JZBjoJfjpuXn1F0foq6qq6shSRKqqqp6Xafrl2BXRqNxSAUFo9GIxsZGr/c8J5/ny783lZWVSExMRENDA+x2O0pKSkJWpxYCOTZ6OW8CNdTPmUDo5ZwJ5nMOlfMmkGPA8yP054cuA5Gno6MkScjJyem3Q7UvnoM+lBUXF6OsrKzP42M2mwF8GwzKy8uRk5PTZ8gcCgZybHzRw3nTH72eM4HSyznT3+fUw3kTyO+a58clSpwfYR2IqqursWvXrn7XKygokA8WcClpWiwWAJcOWkJCApxOp19fboP9BAz02Hjk5+dj0aJF8nHqTfcrJAsXLsTy5ct7veQ5GKh1bHwZzOdNsMdloPR4zgRjMJ8zgHLHpr/PGY7njb8C+V0P9vNDKaqcH0H1QAozjY2NIi8vz6vjld1uFwBEVVWVz208y7sDIGpqakJVqqaqqqoG3EPf13EDIA4fPqx0WYPCQI+N3s4bfztVdzfUz5mBHBu9nDOBfs6hdN4Ecgx4foT+/NBVIPL1ZJDnvb4OtNFo7PE00VDNkjU1NV4nVmNjY69PUjU2NvY4np73htLTIB7+HBsh9HXeDPRLX2/njBD+P2Wmh3PG3885FM+bQH7XPD98U+r80NVTZmazGXl5eV6X1nbt2gWz2Sw/aeZwOOQB+DwKCgq8erdXV1cHdLtksLPZbLDZbDCbzXA4HHA4HCgvL0diYiKAnsfGaDT2OJ7l5eXIzs4eMpewPfw9NoB+zhug98vZej5nPAZ6bAD9nDP9fU49nDf+HoOBbDNUaHV+GP6VvHRDkiSUl5fLrz090T0Hrby8HCUlJfJoxB6lpaXywT548GDYPt3QG0mScPXVV/vsxe85RXwdm+7H8+zZszw2XQz188bhcMh9SGw2G/Ly8rxGh9XrOQMEdmyAoX/OePT1OfVy3vh7DPrbZijR4vzQXSAiIiIi6k5Xt8yIiIiIfGEgIiIiIt1jICIiIiLdYyAiIiIi3WMgIiIiIt1jICIiIiLdYyAiIiIi3QvryV2JiLTWdUC4vLw8jashokDxChER6cby5cvlGbCVYrVacfbsWa/3ysvLsXz5clRXVyu2HyIKLQYiIlKUw+FAfn4+DAYD0tLSUFpaitLSUuTn52P58uVwOByK7ctmsyErKwtpaWn9ritJEoxGI8rKyhSd/yo7O7vH/i0WC8rKynDw4EHF9kNEocVbZkSkKJPJhJKSEthsNphMJq/bSA6HAzfffDPeeustmM3moPdlNpvloEVEFAxeISIi1ZhMJkydOhX333+/Ym0mJiYq1hYR6RevEBGRqoxGo6K3zUKptLS0R/8gAEhKSmIHaqIhhoGIiFQjSRKsViu2bt3aY1lpaSlMJhMcDgdMJhOys7PlbcrLy2EymVBTU4Ply5crcruta02e9j19izIzMwHwqTEiPWEgIqKQcTgc8pNWDocDdrvdZ/+hnJwcLFq0SA5BWVlZMJlMMJvNKC4uxvLly+WQlJaWhsOHDyvWMTo/Px/5+flyGMvPz5cD0UBYrVbU1NRAkiSvIEdE4YV9iIgoZDwBITs7G5mZmTh06BAaGhq81vGEpq5BIicnB2VlZfJyq9Xq1WbX18E6dOgQ8vPz5StTVVVVfm2fmZmJqqoq1NTUMAwRhTFeISIiVZjNZhQUFCAnJweNjY3y+1arFUaj0Svk2O12uZ+RJ6BIkgSHw4GGhoYeoSoYJSUlyMnJQVpaGsxmMw4fPqxY20QUPniFiIhUYzKZ5GDj4bnVlJmZKf+UlJSgpqYGwKWxhnJyclBZWQmTyQSTyaRoTZmZmWhsbMThw4chSRIHUyTSKQYiIlKdzWaT/202m30+dSZJEiRJwuzZs1FQUACLxQKj0SiPMq3Ek2ppaWlyO2azGcuXL1d00EYiCh8MREQUEr5ubXme5PKM4FxeXo7MzExMnTq1x5WZyspKOBwOSJLk1Qnb02bXUBUISZKQnZ3tdcXp4MGDfnWoJqKhg32IiEhRDocDZWVlmDp1KoBLj9N7ru4YjUZUVVWhpKQE1dXVchipqalBfn4+Ghoa5IEWLRYLgEuPvufn5yMrKwvApT5F+fn5WLRoEWw2G4qLi+Wnw0pKSgZcp9FoxLRp0+Qg1tDQ4HM4ACLSB4MQQmhdBBFRqEmShOLiYr9CU7D8DWlEpB3eMiMiIiLdYyAiIiIi3WMgIiLdqK6uxvLly+Un1UKlvLwcOTk5Id0HESmLfYiIiIhI93iFiIiIiHSPgYiIiIh0j4GIiIiIdI+BiIiIiHSPgYiIiIh0j4GIiIiIdI+BiIiIiHSPgYiIiIh07/8DZiralFDYijgAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "\n", "ss.EIG.plot_root_loci(ret, range(30, 51), left=-3, ymax=4, ymin=-4)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Note that the marker size increases linearly as the parameters sweep from the\n", "first to the last. We recommend sweeping the parameter in an ascending order so\n", "that the larger marker size correspond to a larger parameter." ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Cleanup" ] }, { "cell_type": "code", "execution_count": 8, "metadata": { "execution": { "iopub.execute_input": "2021-09-26T22:41:53.557089Z", "iopub.status.busy": "2021-09-26T22:41:53.556654Z", "iopub.status.idle": "2021-09-26T22:41:54.650792Z", "shell.execute_reply": "2021-09-26T22:41:54.650033Z" } }, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "\"/home/hacui/repos/andes/examples/kundur_full_out.txt\" removed.\n", "\"/home/hacui/repos/andes/examples/kundur_full_eig.txt\" removed.\n" ] } ], "source": [ "!andes misc -C" ] } ], "metadata": { "kernelspec": { "display_name": "Python 3.9.10 ('a')", "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.10" }, "vscode": { "interpreter": { "hash": "9e82087c6c724a796d11bbbca032c2fea9f2a8500ea883496f2915648d5c9903" } } }, "nbformat": 4, "nbformat_minor": 4 }