cigreTB880_case01.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "647901f1-6d85-4372-92dc-ff8d2354df61",
   "metadata": {
    "tags": []
   },
   "source": [
    "# CIGRE TB 880 - 2022\n",
    "## Power cable rating examples for calculation tool verification\n",
    "\n",
    "Reference IEC standards:\n",
    "- \"IEC 60287-1-1, Electric Cables - Calculation of the Current Rating - Part 1-1: Current Rating Equations (100 % load factor) and Calculation of Losses – General\", IEC,  2014-11.  \n",
    "- \"IEC 60287-2-1, Electric Cables – Calculation of the Current Rating – Part 2-1: Thermal Resistance – Calculation of Thermal Resistance, Edition 2.0\", IEC, 2015-2."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "f27d4e7f-0172-4478-a595-63ec4a5f2497",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Import libraries\n",
    "import math as mt\n",
    "#import numpy as np\n",
    "#import pandas as pd\n",
    "import matplotlib\n",
    "#matplotlib.use('TkAgg')\n",
    "import matplotlib.pyplot as plt"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3604b992-3112-4825-b45a-617695c13105",
   "metadata": {},
   "source": [
    "# 4. Introductory Case study 0\n",
    "\n",
    "<div>\n",
    "<img src=\"./img/img1_cs0.png\" width=\"150\"/>\n",
    "</div>\n",
    "\n",
    "**Input data:**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "82ce3f24-f0e2-4ac5-b9b8-057ad0175386",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 4.2.1 Grid parameters\n",
    "\n",
    "U = 132 #[kv] Voltage level\n",
    "f = 50 #[Hz] Frequency"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3e1f051e-fd4b-490a-ab4e-24c94a5ffcbe",
   "metadata": {},
   "source": [
    "<div>\n",
    "<img src=\"./img/img2_xlpe132kV.png\" width=\"300\"/> <img src=\"./img/img3_cross.png\" width=\"450\"/>\n",
    "</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "735eb1b1-ac83-46fc-b39c-da838c8f339a",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 4.2.2 Cable parameters\n",
    "\n",
    "# List of the cable composition: [Conductor, Conductor_screen, XLPE_insulation, Insulation_screen, Sheath, Oversheath]\n",
    "thk = [0, 1.5, 15.5, 1.3, 0.8, 3.5] #[mm] Thickness of the layer\n",
    "dia = [0]*6 #[mm] Diammeter after the layer. Initialize\n",
    "dia[0] = 30.3 #[mm] Conductor diameter (630 mm2)\n",
    "\n",
    "for i in range(1,len(dia)):\n",
    "    dia[i] = dia[i-1] + 2*thk[i]\n",
    "\n",
    "# Conductor\n",
    "R_0 = 28.3e-6 #[Ohm/m] DC resistance at 20°C. From IEC 60228 Table 2\n",
    "alpha_20 = 3.93e-3 #[K^-1] Temperature coefficient. From IEC 60287-1-1 Table 1 (Copper)\n",
    "k_s = 1 # Skin effect coefficient (round stranded conductor and extruded insulation system). From IEC 60287-1-1 Table 2\n",
    "k_p = 1 # Proximity effect coefficient (round stranded conductor and extruded insulation system). From IEC 60287-1-1 Table 2\n",
    "\n",
    "# Insulation layer\n",
    "epsilon = 2.5 # Relative permittivity. From IEC 60287-1-1 Table 3\n",
    "tan_delta = 0.001 # Loss factor. From IEC 60287-1-1 Table 3\n",
    "rho_i = 3.5 #[Km/W] Thermal resistivity. From CIGRE TB 640 Table 1\n",
    "rho_sc = 2.5 #[Km/W] Thermal resistivity screen (semi-conducting layer). CIGRE WG B1.56, IEC 60949\n",
    "\n",
    "# Metal sheath\n",
    "t_s = thk[4] #[mm] From Manufacturer datasheet\n",
    "rho_s_20 = 2.84e-8 #[Ohm/m] Resistivity at 20°c. From IEC 60287-1-1 Table 1\n",
    "alpha_s_20 = 4.03e-3 #[K^-1] Temperature coefficient. From IEC 60287-1-1 Table 1 (Aluminium)\n",
    "\n",
    "# Oversheath\n",
    "rho_Tos = 3.5 #[Km/W] Thermal resistivity. From IEC 60287-2-1 Table 1"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "04f470c7-a635-4761-8a4e-2f27c63076a9",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 4.2.5 Operating conditions\n",
    "\n",
    "theta_c_max = 90 #[°C] Maximum operating temperature of conductor. User\n",
    "theta_a = 20 #[°C] Ambient ground temperature. From IEC 60287-3-1 Section 3.2.1\n",
    "rho_Tsoil = 1 #[Km/W] Thermal resistivity. From IEC 60287-3-1 Section 3.2.2"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1c97cea0-4b59-48f4-8c7e-d0072f004986",
   "metadata": {},
   "source": [
    "## 4.4 CASE #0-1: Results independent of the temperature\n",
    "\n",
    "**Case #0-1 without any variation**"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "616e12f1-ee33-4de9-a0f5-016be0c5a1c6",
   "metadata": {},
   "source": [
    "### 4.4.1 Capacitance\n",
    "\n",
    "$ C = \\frac{\\epsilon}{18\\ln(\\frac{D_i}{d_c})} 10^{-9} $ (F/m). From IEC 60287-1-1 Section 2.2"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "993418d5-0a09-43e5-a8d5-1b2211d7999e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "C = 2.1107662202E-10 F/m\n"
     ]
    }
   ],
   "source": [
    "# 4.4.1 Capacitance\n",
    "\n",
    "d_c = dia[1] #[mm] diameter of conductor, including screen, if any\n",
    "D_i = dia[2] #[mm] external diameter of the insulation (excluding screen)\n",
    "\n",
    "C = (epsilon/(18*mt.log(D_i/d_c)))*1e-9 #[F/m] From IEC 60287-1-1 Section 2.2\n",
    "print(\"C = %.10E F/m\" % C)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "7fab0e9f-c7a2-4273-8192-1e2b1914dd64",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "X = 5.0403313985E-05 Ω/m\n"
     ]
    }
   ],
   "source": [
    "# 4.4.2 Reactance\n",
    "\n",
    "omega = 2*mt.pi*f #[rad/s] Angular frequency\n",
    "s = dia[5] #[mm] distance between conductor axes in the electrical section being considered\n",
    "d = D_i + 2*thk[3] + thk[4] #[mm] mean diameter of the sheath\n",
    "\n",
    "X = 2*omega*1e-7*mt.log(2*s/d) #[Ohm/m] From IEC 60287-1-1 Section 2.3.1\n",
    "print(\"X = %.10E \\u03A9/m\" % X)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "8bd37659-8cb9-4329-a1f1-d5b8d57e87be",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Rs0 = 1.6691286499E-04 Ω/m\n"
     ]
    }
   ],
   "source": [
    "# 4.4.3 Electrical resistance of the metal sheath at 20°C\n",
    "\n",
    "A_s = mt.pi*d*t_s #[mm2] cross section\n",
    "\n",
    "R_s0 = rho_s_20/(A_s*1e-6) #[Ohm/m]\n",
    "print(\"Rs0 = %.10E \\u03A9/m\" % R_s0)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "8e009965-6a4a-422c-baf0-9f2a64ed77e1",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "T1 = 0.4198714890 Km/W\n"
     ]
    }
   ],
   "source": [
    "# 4.4.4 Thermal resistance between conductor and sheath\n",
    "\n",
    "t_1 = thk[1:4] #[mm] insulation thickness\n",
    "d_c_arr = dia[0:3] #[mm] diameter\n",
    "rho_T = [rho_sc, rho_i, rho_sc] #[Km/W] Thermal resistivity\n",
    "T1_layer = [0]*len(rho_T) #[Km/W] Thermal resistance\n",
    "\n",
    "for i in range(len(T1_layer)):\n",
    "    T1_layer[i] = (rho_T[i]/(2*mt.pi))*mt.log(1+2*t_1[i]/d_c_arr[i]) #[Km/W] From IEC 60287-2-1 Section 4.1.2.1\n",
    "\n",
    "T1 = sum(T1_layer) #[Km/W]\n",
    "print(\"T1 = %.10f Km/W\" % T1)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "id": "1fa6cb36-8d13-4b68-bb91-8cf052a9de24",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "T3 = 0.0867193748 Km/W\n"
     ]
    }
   ],
   "source": [
    "# 4.4.5 Thermal resistance of the oversheath\n",
    "\n",
    "t_3 = thk[5] #[mm] oversheath thickness\n",
    "D_a = dia[4] #[mm] reference diameter\n",
    "\n",
    "T3 = (1/(2*mt.pi))*rho_Tos*mt.log(1+2*t_3/D_a) #[Km/W] From IEC 60287-2-1 Section 4.1.4.1\n",
    "T3 = 1.6*T3 #[Km/W] From IEC 60287-2-1 Section 4.2.4.3.2\n",
    "print(\"T3 = %.10f Km/W\" % T3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "2bd97694-089a-49fe-b90a-925ea00a3eb7",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "T4 = 1.5946928925 Km/W\n"
     ]
    }
   ],
   "source": [
    "# 4.4.6 External thermal resistance\n",
    "\n",
    "L = 1000 #[mm]\n",
    "u = 2*L/dia[5]\n",
    "rho_T_ex = 1 #[Km/W]\n",
    "\n",
    "T4 = 1.5/mt.pi*rho_T_ex*(mt.log(2*u)-0.630) #[Km/W] From IEC 60287-2-1 Section 4.2.4.3.2\n",
    "print(\"T4 = %.10f Km/W\" % T4)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "2dd3e5ff-c2f1-45ae-952c-0c3d6958b721",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Wd = 0.3851382172 W/m\n"
     ]
    }
   ],
   "source": [
    "# 4.4.7 Dielectric losses\n",
    "\n",
    "U_0 = U/mt.sqrt(3)*1e3 #[V]\n",
    "\n",
    "W_d = omega*C*U_0**2*tan_delta #[W/m] From IEC 60287-1-1 Section 2.2\n",
    "print(\"Wd = %.10f W/m\" % W_d)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6346190d-c4c9-48f4-a3a5-a9d4d7f0e054",
   "metadata": {
    "tags": []
   },
   "source": [
    "## 4.5 case #0-1: Results dependent on the temperature\n",
    "### 4.5.1  Results linked to the maximum operating temperature\n",
    "### 4.5.2 First iteration & 4.5.3 Next iterations"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "871b93a1-733c-4503-abc3-37f9555caa4c",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Iteration:  1\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0727239574E-04 Ohm/m\n",
      "λ1 = 0.2928142510\n",
      "I = 822.0666849587 A\n",
      "Wc = 26.7083959538 W/m\n",
      "Ws = 7.8205989571 W/m\n",
      "θj = 75.6773199466 °C\n",
      "θs = 78.7050517428 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration:  2\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0640133734E-04 Ohm/m\n",
      "λ1 = 0.2939111934\n",
      "I = 821.7745414009 A\n",
      "Wc = 26.6894162298 W/m\n",
      "Ws = 7.8443181760 W/m\n",
      "θj = 75.6848779855 °C\n",
      "θs = 78.7130207878 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration:  3\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0640669779E-04 Ohm/m\n",
      "λ1 = 0.2939044196\n",
      "I = 821.7763444974 A\n",
      "Wc = 26.6895333511 W/m\n",
      "Ws = 7.8441718079 W/m\n",
      "θj = 75.6848313459 °C\n",
      "θs = 78.7129716119 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration:  4\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0640666471E-04 Ohm/m\n",
      "λ1 = 0.2939044614\n",
      "I = 821.7763333709 A\n",
      "Wc = 26.6895326284 W/m\n",
      "Ws = 7.8441727111 W/m\n",
      "θj = 75.6848316337 °C\n",
      "θs = 78.7129719153 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration:  5\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0640666491E-04 Ohm/m\n",
      "λ1 = 0.2939044611\n",
      "I = 821.7763334396 A\n",
      "Wc = 26.6895326329 W/m\n",
      "Ws = 7.8441727055 W/m\n",
      "θj = 75.6848316319 °C\n",
      "θs = 78.7129719134 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration:  6\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0640666491E-04 Ohm/m\n",
      "λ1 = 0.2939044611\n",
      "I = 821.7763334392 A\n",
      "Wc = 26.6895326328 W/m\n",
      "Ws = 7.8441727056 W/m\n",
      "θj = 75.6848316319 °C\n",
      "θs = 78.7129719135 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration:  7\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0640666491E-04 Ohm/m\n",
      "λ1 = 0.2939044611\n",
      "I = 821.7763334392 A\n",
      "Wc = 26.6895326328 W/m\n",
      "Ws = 7.8441727056 W/m\n",
      "θj = 75.6848316319 °C\n",
      "θs = 78.7129719135 °C\n",
      "θc = 90.0000000000 °C\n",
      "Convergence reached after  7 iterations\n"
     ]
    }
   ],
   "source": [
    "# Criteria of convergence\n",
    "tempTol = 1e-11 #[K] Temperature tolerance threshold for convergence\n",
    "currTol = 1e-11 #[A] Current tolerance threshold for convergence\n",
    "maxIte = 10 # Maximum number of iterations\n",
    "\n",
    "# Other parameters\n",
    "n = 1 # Single-core cable\n",
    "d_c = dia[0] # Conductor diameter\n",
    "lambda_2 = 0 # The cable is unarmoured\n",
    "T2 = 0 # The cable is unarmoured\n",
    "# NOTE: IEC ArmoUred Cable is not the same that NEC(NFPA 70) Armored Cable (Type AC - Article 320), which is similar to Metal-Clad Cable (Type MC Article 330)\n",
    "\n",
    "# Variable Preallocation\n",
    "R = [0]*maxIte #[Ohm/m] AC resistance of conductor\n",
    "R_s = [0]*maxIte #[Ohm/m] Resistance of sheath\n",
    "lambda_1 = [0]*maxIte # Loss factor\n",
    "I = [0]*maxIte #[A] Current rating\n",
    "W_c = [0]*maxIte #[W/m] Conductor losses\n",
    "W_s = [0]*maxIte #[W/m] Sheath losses\n",
    "theta_c = [0]*maxIte #[°C] Conductor operation temperature\n",
    "theta_s = [0]*maxIte #[°C] Sheath temperature\n",
    "\n",
    "# Initial values\n",
    "theta_c[0] = theta_c_max # Maximum allowable conductor temperature\n",
    "theta_s[0] = theta_c_max-10 # Assumption to initializate!\n",
    "\n",
    "for i in range(maxIte):\n",
    "    print(\"Iteration: %2d\" % (i+1))\n",
    "    # 4.5.1  Results linked to the maximum operating temperature\n",
    "    # 4.5.1.1 DC resistance of conductor\n",
    "    R_prime = R_0*(1+alpha_20*(theta_c[i]-20)) #[Ohm/m] From IEC 60287-1-1 Section 2.1.1\n",
    "    #print(\"R' = %.10E W/m\" % R_prime)\n",
    "    \n",
    "    # 4.5.1.2 Skin effect factor\n",
    "    x_s = mt.sqrt(8*mt.pi*f/R_prime*1e-7*k_s) # From IEC 60287-1-1 Section 2.1.2\n",
    "    # For 0 < x_s <= 2,8\n",
    "    y_s = x_s**4/(192+0.8*x_s**4)\n",
    "    #print(\"xs = %.10f\\nys = %.10f\" % (x_s, y_s))\n",
    "    \n",
    "    # 4.5.1.3 Proximity effect factor \n",
    "    x_p = mt.sqrt(8*mt.pi*f/R_prime*1e-7*k_p) # From IEC 60287-1-1 Section 2.1.4.1\n",
    "    y_p = x_p**4/(192+0.8*x_p**4)*(d_c/s)**2*(0.312*(d_c/s)**2+1.18/(x_p**4/(192+0.8*x_p**4)+0.27))\n",
    "    #print(\"xp = %.10f\\nyp = %.10f\" % (x_p, y_p))\n",
    "    \n",
    "    # 4.5.1.4 AC resistance of conductor \n",
    "    R[i] = R_prime*(1+y_s+y_p) #[Ohm/m] From IEC 60287-1-1 Section 2.1\n",
    "    print(\"R = %.10E Ohm/m\" % R[i])\n",
    "    \n",
    "    # 4.5.2 First iteration & 4.5.3 Next iterations\n",
    "    # 4.5.2.1 Resistance of sheath at operating temperature\n",
    "    R_s[i] = R_s0*(1+alpha_s_20*(theta_s[i]-20)) #[Ohm/m] From IEC 60287-1-1 Section 2.3\n",
    "    print(\"Rs = %.10E Ohm/m\" % R_s[i])\n",
    "    \n",
    "    # 4.5.2.2 Losses caused by circulating currents\n",
    "    lambda_1_prime = (R_s[i]/R[i])*(1/(1+(R_s[i]/X)**2)) # From IEC 60287-1-1 Section 2.3.1\n",
    "    lambda_1_prime2 = 0\n",
    "    lambda_1[i] = lambda_1_prime + lambda_1_prime2\n",
    "    print(\"\\u03BB1 = %.10f\" % lambda_1[i])\n",
    "    \n",
    "    # 4.5.2.3 Permissible current rating\n",
    "    Delta_theta = theta_c[i] - theta_a\n",
    "    I[i] = mt.sqrt((Delta_theta-W_d*(0.5*T1+n*(T2+T3+T4)))/(R[i]*T1+n*R[i]*(1+lambda_1[i])*T2+n*R[i]*(1+lambda_1[i]+lambda_2)*(T3+T4))) #[A] From IEC 60287-1-1 Section 1.4.1.1\n",
    "    print(\"I = %.10f A\" % I[i])\n",
    "    \n",
    "    # 4.5.2.4 Losses in conductor and sheath\n",
    "    W_c[i] = R[i]*I[i]**2 #[W/m]\n",
    "    W_s[i] = lambda_1[i]*W_c[i] #[W/m]\n",
    "    print(\"Wc = %.10f W/m\\nWs = %.10f W/m\" % (W_c[i], W_s[i]))\n",
    "    \n",
    "    # 4.5.2.5 Temperature on cable components\n",
    "    theta_j = theta_a + n*(W_c[i]+W_s[i]+W_d)*T4 #[°C] Jacket. From IEC 60287-1-1 Section 1.4.1.1\n",
    "    theta_s[i] = theta_j + n*(W_c[i]+W_s[i]+W_d)*T3 #[°C] Metal sheath\n",
    "    theta_c[i] = theta_s[i] + n*(W_c[i]+W_d/2)*T1 #[°C] Conductor\n",
    "    print(\"\\u03B8j = %.10f °C\\n\\u03B8s = %.10f °C\\n\\u03B8c = %.10f °C\" % (theta_j, theta_s[i], theta_c[i]))\n",
    "    \n",
    "    # Convergence evaluation\n",
    "    Delta_theta_s = abs(theta_s[i] - theta_s[i-1])\n",
    "    Delta_theta_c = abs(theta_c[i] - theta_c[i-1])\n",
    "    Delta_theta1 = max(Delta_theta_s, Delta_theta_c)\n",
    "    Delta_I = abs(I[i] - I[i-1])\n",
    "    \n",
    "    if (Delta_theta1 < tempTol) and (Delta_I < currTol):\n",
    "        print(\"Convergence reached after %2d iterations\" % (i+1))\n",
    "        break\n",
    "    elif i+1 == maxIte:\n",
    "        print(\"Convergence criteria not reached after %2d iterations\" % maxIte)\n",
    "    else:\n",
    "        theta_s[i+1] = theta_s[i] # Update for next iteration\n",
    "        theta_c[i+1] = theta_c[i]"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "637311d6-9f48-47e4-8413-62b6267f9d5e",
   "metadata": {},
   "source": [
    "### Results"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "5957dc03-af4f-4841-870a-9018df3a14ad",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAlIAAAGwCAYAAABiu4tnAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABNL0lEQVR4nO3de1xUdf4/8NdhgOE+iMpllBkg74gCmuUltUzykpZl3io119YyS3It7OJqGpJW/tp01bWtvOWu361sXSsNr+XWmjCoaCaaCCgQgswMgtxmzu8PZHICdWYcOHN5PR+PeXydw5lz3jO73++8vp/znvcRRFEUQURERERW85C6ACIiIiJnxSBFREREZCMGKSIiIiIbMUgRERER2YhBioiIiMhGDFJERERENmKQIiIiIrKRp9QFuDKj0YjCwkIEBgZCEASpyyEiIiILiKKIiooKKJVKeHjcfM2JQaoFFRYWIjIyUuoyiIiIyAYFBQXo2LHjTfdhkGpBgYGBABr+gwgKCpK4GiIiIrKEXq9HZGSk6Xv8ZhikWlDj5bygoCAGKSIiIidjSVsOm82JiIiIbMQgRURERGQjBikiIiIiGzFIEREREdmIQYqIiIjIRgxSRERERDZikCIiIiKyEYMUERERkY0YpIiIiIhsxCBFREREZCMGKSIiIiIbMUgRERER2YhBygmJooj8sipc1F6VuhQiIiK3xiDlhFK/PIXBb+/Hhv/mSl0KERGRW2OQckLdIoIAAJp8rbSFEBERuTkGKSfUR90GAJB9UYeaeoPE1RAREbkvBiknFNXWDyH+3qitN+JkoV7qcoiIiNwWg5QTEgQBCZHBAABNXrm0xRAREbkxBiknlXjt8l4W+6SIiIgkwyDlpBJVDUFKk88VKSIiIqkwSDmp3pEKyDwEFOmqUch5UkRERJJgkHJSft6e6BYeCICrUkRERFJhkHJipst7eVppCyEiInJTDFJOrHGeFFekiIiIpMEg5cQaV6ROFupQXcfBnERERK2NQcqJRYb4ol2AN+oMIk5c1EldDhERkdthkHJigiBwDAIREZGEGKScXONgTjacExERtT4GKSfXuCKVmV8OURQlroaIiMi9SBqk6uvr8frrryM6Ohq+vr6IiYnBkiVLYDQaAQB1dXVISUlBXFwc/P39oVQqMXXqVBQWFpqOcfnyZTz//PPo2rUr/Pz8oFKp8MILL0Cnu3XP0Jo1axAdHQ0fHx/06dMH3333ndnfp0+fDkEQzB533323fT+E29SrowKeHgIuVdTgQjkHcxIREbUmSYPU8uXLsW7dOqxevRqnTp3CihUr8Pbbb2PVqlUAgKqqKmg0GixcuBAajQaff/45cnJyMHbsWNMxCgsLUVhYiHfeeQfZ2dnYsGEDdu3ahT/84Q83Pfe2bduQnJyM1157DVlZWbjnnnswcuRI5Ofnm+03YsQIFBUVmR5fffWV/T+I2+DjJUOsMggA+6SIiIhamyBKeD3owQcfRFhYGD788EPTtkcffRR+fn7YvHlzs685cuQI+vXrh7y8PKhUqmb3+de//oUnnngClZWV8PT0bHafu+66C4mJiVi7dq1pW/fu3fHwww8jLS0NQMOKlFarxRdffGHT+9Pr9VAoFNDpdAgKCrLpGJZYvOMkNnx/HtMHRGHx2NgWOw8REZE7sOb7W9IVqUGDBmHv3r3IyckBABw7dgyHDh3CqFGjbvganU4HQRAQHBx8032CgoJuGKJqa2uRmZmJpKQks+1JSUn4/vvvzbYdOHAAoaGh6NKlC55++mmUlJTc8Lw1NTXQ6/Vmj9bQ2HCemccVKSIiotbUfNJoJSkpKdDpdOjWrRtkMhkMBgNSU1MxefLkZvevrq7GggULMGXKlBsmxLKyMixduhSzZs264XlLS0thMBgQFhZmtj0sLAzFxcWm5yNHjsRjjz0GtVqN3NxcLFy4EPfddx8yMzMhl8ubHDctLQ1vvPGGJW/drhJVwQCAU0V6XK01wNdb1uo1EBERuSNJV6S2bduGLVu2YOvWrdBoNNi4cSPeeecdbNy4scm+dXV1mDRpEoxGI9asWdPs8fR6PUaPHo0ePXpg0aJFtzy/IAhmz0VRNNs2ceJEjB49Gj179sSYMWPw9ddfIycnB19++WWzx3vllVeg0+lMj4KCglvWYA8dgn0RFiRHvVHE8QvaVjknERERSbwi9dJLL2HBggWYNGkSACAuLg55eXlIS0vDtGnTTPvV1dVhwoQJyM3Nxb59+5pdjaqoqMCIESMQEBCA7du3w8vL64bnbdeuHWQymdnqEwCUlJQ0WaW6XkREBNRqNc6cOdPs3+VyebMrVS2tcTDn1yeKkZlfjrti2rZ6DURERO5I0hWpqqoqeHiYlyCTyUzjD4DfQtSZM2ewZ88etG3bNCTo9XokJSXB29sbO3bsgI+Pz03P6+3tjT59+iA9Pd1se3p6OgYMGHDD15WVlaGgoAARERGWvL1WZZpwzsGcRERErUbSFakxY8YgNTUVKpUKsbGxyMrKwsqVKzFjxgwADXOmxo8fD41Gg507d8JgMJhWkUJCQuDt7Y2KigokJSWhqqoKW7ZsMWvybt++PWSyhn6hYcOGYdy4cZgzZw4AYN68eXjyySfRt29f9O/fH+vXr0d+fj6eeeYZAMCVK1ewePFiPProo4iIiMD58+fx6quvol27dhg3blxrf1S3lKgOBgBkXRvM+fvLlkRERGR/kgapVatWYeHChZg9ezZKSkqgVCoxa9Ys/PnPfwYAXLhwATt27AAAxMfHm712//79GDp0KDIzM3H48GEAQKdOncz2yc3NRVRUFADgl19+QWlpqelvEydORFlZGZYsWYKioiL07NkTX331FdRqNYCGlbHs7Gxs2rQJWq0WERERuPfee7Ft2zYEBga2xMdxW2KVCnjLPFBWWYv8y1VQt/WXuiQiIiKXJ+kcKVfXWnOkGo1b819k5WuxckJvPJLYscXPR0RE5IqcZo4U2ZepT4oTzomIiFoFg5QLYcM5ERFR62KQciGNDec/F+tRWVMvbTFERERugEHKhUQofKFU+MAoAscKtFKXQ0RE5PIYpFxMgpp9UkRERK2FQcrF/NZwrpW2ECIiIjfAIOVi+ly3IsXJFkRERC2LQcrF9IgIgtzTA9qqOpwrrZS6HCIiIpfGIOVivD09ENdBAQDQ5LFPioiIqCUxSLmgRDX7pIiIiFoDg5QLamw4z+Iv94iIiFoUg5QLahzMefrXCuir66QthoiIyIUxSLmg0EAfdGzjC5GDOYmIiFoUg5SL4n33iIiIWh6DlIvqwwnnRERELY5BykX9NuG8HEYjB3MSERG1BAYpF9UtIhA+Xh6oqK7HL5euSF0OERGRS2KQclFeMg/07hgMgJf3iIiIWgqDlAszDeZkwzkREVGLYJByYY19UplckSIiImoRDFIuLEEVDAA4W3IFuioO5iQiIrI3BikX1i5Ajqi2fgCArAKuShEREdkbg5SL+20wJ4MUERGRvTFIubgE02BOrbSFEBERuSAGKReXeK1P6miBFgYO5iQiIrIrBikX1zUsEP7eMlypqceZkgqpyyEiInIpDFIuzlPmgd6RwQCATPZJERER2RWDlBv4reFcK20hRERELoZByg0kqoMBAFkczElERGRXDFJuICGyYUXqXGklLlfWSlwNERGR62CQcgNt/L0R094fAFeliIiI7IlByk2Y+qQYpIiIiOyGQcpNsOGciIjI/hik3ESfaxPOj13Qot5glLgaIiIi18Ag5SY6hwYgUO6JqloDfi7mYE4iIiJ7YJByEx4eAuKv3S6GDedERET2wSDlRhJUvIExERGRPTFIuZHGPin+co+IiMg+GKTcSPy1e+7llVWh9EqNtMUQERG5AAYpN6Lw9ULn0AAAgIY3MCYiIrptDFJuJpF9UkRERHbDIOVm2CdFRERkPwxSbiZRHQwAOH5BizoO5iQiIrotDFJuJqZdAIJ8PFFdZ8SpIr3U5RARETk1Bik34+EhILHx8h4bzomIiG4Lg5QbYsM5ERGRfTBIuaHGIJXJFSkiIqLbImmQqq+vx+uvv47o6Gj4+voiJiYGS5YsgdHY0ARdV1eHlJQUxMXFwd/fH0qlElOnTkVhYaHpGJcvX8bzzz+Prl27ws/PDyqVCi+88AJ0Ot0tz79mzRpER0fDx8cHffr0wXfffWf2d1EUsXjxYiiVSvj6+mLo0KE4efKkfT8ECfSOVEAQgIvaqyjRV0tdDhERkdOSNEgtX74c69atw+rVq3Hq1CmsWLECb7/9NlatWgUAqKqqgkajwcKFC6HRaPD5558jJycHY8eONR2jsLAQhYWFeOedd5CdnY0NGzZg165d+MMf/nDTc2/btg3Jycl47bXXkJWVhXvuuQcjR45Efn6+aZ8VK1Zg5cqVWL16NY4cOYLw8HAMHz4cFRUVLfOBtJJAHy90DQsEwDEIREREt0MQRVGU6uQPPvggwsLC8OGHH5q2Pfroo/Dz88PmzZubfc2RI0fQr18/5OXlQaVSNbvPv/71LzzxxBOorKyEp6dns/vcddddSExMxNq1a03bunfvjocffhhpaWkQRRFKpRLJyclISUkBANTU1CAsLAzLly/HrFmzbvn+9Ho9FAoFdDodgoKCbrl/a3p1eza2Hs7H0/dE47XRPaQuh4iIyGFY8/0t6YrUoEGDsHfvXuTk5AAAjh07hkOHDmHUqFE3fI1Op4MgCAgODr7pPkFBQTcMUbW1tcjMzERSUpLZ9qSkJHz//fcAgNzcXBQXF5vtI5fLMWTIENM+v1dTUwO9Xm/2cFRsOCciIrp9zSeNVpKSkgKdTodu3bpBJpPBYDAgNTUVkydPbnb/6upqLFiwAFOmTLlhQiwrK8PSpUtvumJUWloKg8GAsLAws+1hYWEoLi4GANP/bG6fvLy8Zo+blpaGN95444bndSSJqmAAQPZFHWrrjfD25O8OiIiIrCXpt+e2bduwZcsWbN26FRqNBhs3bsQ777yDjRs3Ntm3rq4OkyZNgtFoxJo1a5o9nl6vx+jRo9GjRw8sWrTolucXBMHsuSiKTbZZsk+jV155BTqdzvQoKCi4ZQ1SiW7njzZ+XqitN+Jk4a0b84mIiKgpSVekXnrpJSxYsACTJk0CAMTFxSEvLw9paWmYNm2aab+6ujpMmDABubm52LdvX7OrURUVFRgxYgQCAgKwfft2eHl53fC87dq1g0wmM606NSopKTGtQIWHhwNoWJmKiIhodp/fk8vlkMvlFr57aQmCgERVG+z9uQSZeeVIuHapj4iIiCwn6YpUVVUVPDzMS5DJZKbxB8BvIerMmTPYs2cP2rZt2+Q4er0eSUlJ8Pb2xo4dO+Dj43PT83p7e6NPnz5IT083256eno4BAwYAAKKjoxEeHm62T21tLQ4ePGjax9k1TjjPYp8UERGRTSRdkRozZgxSU1OhUqkQGxuLrKwsrFy5EjNmzADQMGdq/Pjx0Gg02LlzJwwGg2kVKSQkBN7e3qioqEBSUhKqqqqwZcsWsybv9u3bQyaTAQCGDRuGcePGYc6cOQCAefPm4cknn0Tfvn3Rv39/rF+/Hvn5+XjmmWcANKzYJCcnY9myZejcuTM6d+6MZcuWwc/PD1OmTGntj6pFJFzrk+IIBCIiIttIGqRWrVqFhQsXYvbs2SgpKYFSqcSsWbPw5z//GQBw4cIF7NixAwAQHx9v9tr9+/dj6NChyMzMxOHDhwEAnTp1MtsnNzcXUVFRAIBffvkFpaWlpr9NnDgRZWVlWLJkCYqKitCzZ0989dVXUKvVpn1efvllXL16FbNnz0Z5eTnuuusufPPNNwgMDLT3RyGJ3h2DIfMQUKSrRqH2KpTBvlKXRERE5FQknSPl6hx5jlSj0e9/h5OFeqyekoAHeymlLoeIiEhyTjNHiqRnmieVp5W2ECIiIifEIOXmEtXBANgnRUREZAsGKTfXRxUCADhZqEN1nUHiaoiIiJwLg5SbiwzxRbsAb9QZRJy4yMGcRERE1mCQcnOCIJiGcfLyHhERkXUYpIgN50RERDZikCL0uTbhPDO/HJyGQUREZDkGKUKvjgp4egi4VFGDC+VXpS6HiIjIaVg02TwkJMSqgwqCAI1GYzYlnByXj5cMPZRBOH5BB01+OSJD/KQuiYiIyClYFKS0Wi3ee+89KBSKW+4riiJmz54Ng4E/pXcmiao2OH5Bh6x8LR6K7yB1OURERE7B4nvtTZo0CaGhoRbt+/zzz9tcEEkjUd0GG74/z1/uERERWcGiIGU0Gq06aEVFhU3FkHQSVcEAgJ8K9bhaa4Cvt0zagoiIiJyA3ZrNDQYDvvjiC3sdjlpZh2BfhAbKUW8UcfyCVupyiIiInMJtB6mff/4ZL7/8MpRKJSZMmGCPmkgCgiD8Nk8qXyttMURERE7CpiBVWVmJjz76CAMHDkRsbCw0Gg1SU1NRWFho7/qoFTXOk2KfFBERkWUsbjYHgB9++AF///vf8X//93/o3LkzHn/8cRw+fBjvv/8+evTo0VI1UitJVAcDADR5DYM5BUGQtiAiIiIHZ/GKVI8ePTB58mSEhYXh8OHD0Gg0+NOf/sQvWxcSq1TASyagrLIW+ZerpC6HiIjI4VkcpM6ePYvBgwfj3nvvRffu3VuyJpKIj5cMPTs0zArj5T0iIqJbszhI5ebmomvXrnj22WfRsWNHzJ8/H1lZWVyRcjG8gTEREZHlLA5SHTp0wGuvvYazZ89i8+bNKC4uxsCBA1FfX48NGzYgJyenJeukVtIYpDLzuCJFRER0Kzb9au++++7Dli1bUFRUhNWrV2Pfvn3o1q0bevXqZe/6qJU1Npz/XKxHZU29tMUQERE5uNuaI6VQKDB79mxkZGRAo9Fg6NChdiqLpBKh8IVS4QOjCBzjYE4iIqKbsttk8/j4eLz//vv2OhxJKKFxnhQv7xEREd2URUEqMTER5eWWf6kOGjQIFy9etLkokhYnnBMREVnGooGcR48exbFjxxASEmLRQY8ePYqamprbKoyk03gD46x8DuYkIiK6GYsnmw8bNgyiKFq0L794nVusUgFvTw+UV9Uht7QSMe0DpC6JiIjIIVkUpHJzc60+cMeOHa1+DTkGb08P9OqgQEZeOTLzyhmkiIiIbsCiIKVWq1u6DnIwieo2yMgrhyZfi8f6RkpdDhERkUOy26/2yLVc3ydFREREzWOQomY1/nLv9K8VqKiuk7gaIiIix8QgRc0KDfJBxza+EEXgaIFW6nKIiIgcEoMU3RBvYExERHRzVgepmJgYlJWVNdmu1WoRExNjl6LIMTT2SWnYJ0VERNQsq4PU+fPnYTAYmmyvqanhNHMX00fdMIA1K78cRqNlM8SIiIjcicUDOXfs2GH69+7du6FQKEzPDQYD9u7di6ioKLsWR9LqFhEIHy8P6Kvr8culK+gcFih1SURERA7F4iD18MMPA2iYWj5t2jSzv3l5eSEqKgrvvvuuXYsjaXnJPNCrYzB+zL0MTX45gxQREdHvWHxpz2g0wmg0QqVSoaSkxPTcaDSipqYGp0+fxoMPPtiStZIE2HBORER0YxavSDWy5XYx5Lz6qK8FKTacExERNWF1kAKAvXv3Yu/evaaVqet99NFHdimMHEPCtV/unSm5Al1VHRR+XtIWRERE5ECs/tXeG2+8gaSkJOzduxelpaUoLy83e5BraRcgh7qtHwAgq4D/+RIREV3P6hWpdevWYcOGDXjyySdboh5yQImqNsgrq4ImX4uhXUOlLoeIiMhhWL0iVVtbiwEDBrRELeSgEq/1SfEGxkREROasDlIzZ87E1q1bW6IWclCNE86z8rUwcDAnERGRidWX9qqrq7F+/Xrs2bMHvXr1gpeXefPxypUr7VYcOYauYYHw85bhSk09zpRUoFt4kNQlEREROQSrg9Tx48cRHx8PADhx4oTZ3wRBsEtR5Fg8ZR7o3TEYP5wrgyZPyyBFRER0jdVBav/+/S1RBzm4Puo2DUEqvxxT7lJJXQ4REZFDsLpHqtHZs2exe/duXL16FQAgiuydcWWJ6mAAgCaPDedERESNrA5SZWVlGDZsGLp06YJRo0ahqKgIQEMT+p/+9CerjlVfX4/XX38d0dHR8PX1RUxMDJYsWWIa8llXV4eUlBTExcXB398fSqUSU6dORWFhodlx1q9fj6FDhyIoKAiCIECr1d7y3BUVFUhOToZarYavry8GDBiAI0eOmO0zffp0CIJg9rj77ruteo+uIiGy4Zd750orUV5ZK3E1REREjsHqIPXiiy/Cy8sL+fn58PPzM22fOHEidu3aZdWxli9fjnXr1mH16tU4deoUVqxYgbfffhurVq0CAFRVVUGj0WDhwoXQaDT4/PPPkZOTg7Fjx5odp6qqCiNGjMCrr75q8blnzpyJ9PR0bN68GdnZ2UhKSsL999+Pixcvmu03YsQIFBUVmR5fffWVVe/RVbTx90ZMe38AHMxJRETUyOoeqW+++Qa7d+9Gx44dzbZ37twZeXl5Vh3rhx9+wEMPPYTRo0cDAKKiovCPf/wDGRkZAACFQoH09HSz16xatQr9+vVDfn4+VKqGXp3k5GQAwIEDByw679WrV/HZZ5/h3//+NwYPHgwAWLx4Mb744gusXbsWb775pmlfuVyO8PBwi45bU1ODmpoa03O9Xm/R65xFoqoNzl2qRGZeOe7rFiZ1OURERJKzekWqsrLSbCWqUWlpKeRyuVXHGjRoEPbu3YucnBwAwLFjx3Do0CGMGjXqhq/R6XQQBAHBwcFWnet69fX1MBgM8PHxMdvu6+uLQ4cOmW07cOAAQkND0aVLFzz99NMoKSm54XHT0tKgUChMj8jISJtrdESJqms3MM7TSlsIERGRg7A6SA0ePBibNm0yPRcEAUajEW+//Tbuvfdeq46VkpKCyZMno1u3bvDy8kJCQgKSk5MxefLkZvevrq7GggULMGXKFAQF2f4T/MDAQPTv3x9Lly5FYWEhDAYDtmzZgsOHD5t6vgBg5MiR+OSTT7Bv3z68++67OHLkCO677z6zVafrvfLKK9DpdKZHQUGBzTU6osaG82MXtKg3GG++MxERkRuw+tLe22+/jaFDhyIjIwO1tbV4+eWXcfLkSVy+fBn//e9/rTrWtm3bsGXLFmzduhWxsbE4evQokpOToVQqMW3aNLN96+rqMGnSJBiNRqxZs8baspvYvHkzZsyYgQ4dOkAmkyExMRFTpkyBRqMx7TNx4kTTv3v27Im+fftCrVbjyy+/xCOPPNLkmHK53OpVOWfSOTQQgXJPVNTU4/SvFYhVKqQuiYiISFJWr0j16NEDx48fR79+/TB8+HBUVlbikUceQVZWFu644w6rjvXSSy9hwYIFmDRpEuLi4vDkk0/ixRdfRFpamtl+dXV1mDBhAnJzc5Genn5bq1GN7rjjDhw8eBBXrlxBQUEBfvzxR9TV1SE6OvqGr4mIiIBarcaZM2du+/zOSOYhIP7a7WI4BoGIiMjKFam6ujokJSXhb3/7G954443bPnlVVRU8PMyznEwmM40/aDznhAkTcObMGezfvx9t27a97fNez9/fH/7+/igvL8fu3buxYsWKG+5bVlaGgoICRERE2LUGZ5KgaoPvzpRCk6/Fk/2lroaIiEhaVgUpLy8vnDhxwm63ghkzZgxSU1OhUqkQGxuLrKwsrFy5EjNmzADQ0BQ+fvx4aDQa7Ny5EwaDAcXFxQCAkJAQeHt7AwCKi4tRXFyMs2fPAgCys7MRGBgIlUqFkJAQAMCwYcMwbtw4zJkzBwCwe/duiKKIrl274uzZs3jppZfQtWtXPPXUUwCAK1euYPHixXj00UcRERGB8+fP49VXX0W7du0wbtw4u7x/Z9R4A2NNPlekiIiIrL60N3XqVHz44Yd2OfmqVaswfvx4zJ49G927d8f8+fMxa9YsLF26FABw4cIF7NixAxcuXEB8fDwiIiJMj++//950nHXr1iEhIQFPP/00gIaG+ISEBOzYscO0zy+//ILS0lLTc51Oh+eeew7dunXD1KlTMWjQIHzzzTemmzDLZDJkZ2fjoYceQpcuXTBt2jR06dIFP/zwAwIDA+3y/p1RwrVf7uWVVaH0SvNN90RERO5CEK28t8vzzz+PTZs2oVOnTujbty/8/f3N/r5y5Uq7FujM9Ho9FAoFdDqdXfq6HMXwlQdxpuQK1j/ZB0mxls3YIiIichbWfH9b/au9EydOIDExEQBM858a2euSHzm2RFUbnCm5Ak2+lkGKiIjcmlVBymAwYPHixYiLizP1HpH7SVQHY1tGAfukiIjI7VnVIyWTyfDAAw9Ap9O1VD3kBPqoG/qkjl/Qoo6DOYmIyI1Z3WweFxeHc+fOtUQt5CRi2gUgyMcT1XVGnCpyrfsJEhERWcPqIJWamor58+dj586dKCoqgl6vN3uQ6/PwEEy/3uNgTiIicmdWN5uPGDECADB27Fiz5nJRFCEIAgwGg/2qI4eVqGqDgzmXoMnXYvpAqashIiKShtVBav/+/S1RBzmZxj4pNpwTEZE7szpIDRkypCXqICfTO1IBQQAulF9Fib4aoUE+UpdERETU6qwOUt9+++1N/z548GCbiyHnEejjha5hgfi5uAKa/HKM6Om+9x8kIiL3ZXWQGjp0aJNt1/dKsUfKfSSo2lwLUloGKSIicktW/2qvvLzc7FFSUoJdu3bhzjvvxDfffNMSNZKDMvVJ8Zd7RETkpqxekVIoFE22DR8+HHK5HC+++CIyMzPtUhg5vkRVMADg+EUdauuN8Pa0OpcTERE5Nbt987Vv3x6nT5+21+HICUS380cbPy/U1htxspDT7omIyP1YvSJ1/Phxs+eiKKKoqAhvvfUWevfubbfCyPEJQsNgzn0/l0CTrzUN6SQiInIXVgep+Ph4CIIAURTNtt9999346KOP7FYYOYc+6sYgVY4/IFrqcoiIiFqV1UEqNzfX7LmHhwfat28PHx/OEXJHCdf6pNhwTkRE7sjqIKVWq1uiDnJSvTsGw0MAinTVKNJdRYTCV+qSiIiIWo3VzeYvvPAC3n///SbbV69ejeTkZHvURE7EX+6J7hFBAABNnlbaYoiIiFqZ1UHqs88+w8CBTe9SO2DAAHz66ad2KYqcS6KK990jIiL3ZHWQKisra3aWVFBQEEpLS+1SFDmXRHUwACCTfVJERORmrA5SnTp1wq5du5ps//rrrxETE2OXosi5NK5InSzUobqOtwgiIiL3YXWz+bx58zBnzhxcunQJ9913HwBg7969ePfdd/Hee+/Zuz5yAqoQP7QL8EbplVqcLNShjzpE6pKIiIhahdVBasaMGaipqUFqaiqWLl0KAIiKisLatWsxdepUuxdIjq9xMGf6T78iM6+cQYqIiNyG1UEKAJ599lk8++yzuHTpEnx9fREQEGDvusjJJF4LUvzlHhERuRObglSj9u3b26sOcnKNNzDW5JdDFEUIgiBtQURERK3AbjctJvfWq2MwPD0ElFTU4KL2qtTlEBERtQoGKbILX28ZeigbBnNyDAIREbkLBimym8YxCFn5WmkLISIiaiW3FaSqq6vtVQe5gITr+qSIiIjcgdVBymg0YunSpejQoQMCAgJw7tw5AMDChQvx4Ycf2r1Ach591A0rUj8V6nG1loM5iYjI9VkdpN58801s2LABK1asgLe3t2l7XFwc/v73v9u1OHIuHYJ9ERooR71RxPELWqnLISIianFWB6lNmzZh/fr1ePzxxyGTyUzbe/XqhZ9//tmuxZFzEQThuhsYa6UthoiIqBVYHaQuXryITp06NdluNBpRV1dnl6LIeTXewJh9UkRE5A6sDlKxsbH47rvvmmz/17/+hYSEBLsURc6rsU8q69pgTiIiIldm9WTzRYsW4cknn8TFixdhNBrx+eef4/Tp09i0aRN27tzZEjWSE4lVKuAlE1B6pRb5l6ugbusvdUlEREQtxuoVqTFjxmDbtm346quvIAgC/vznP+PUqVP4z3/+g+HDh7dEjeREfLxkiFUqAPDyHhERuT6b7rX3wAMP4IEHHrB3LeQiElVtcLRAC02eFuMSOkpdDhERUYvhZHOyu8Y+Ka5IERGRq7NoRapNmzYQBMGiA16+fPm2CiLn1/jLvVNFelTW1MNfbtPCJxERkcOz6Bvuvffea+EyyJVEKHwRofBBka4axy5oMeCOdlKXRERE1CIsClLTpk1r6TrIxSSq2uDL7CJk5TNIERGR67LpmovBYMD27dtx6tQpCIKA7t2746GHHoKnJy/hUINEdUOQ0uSxT4qIiFyX1cnnxIkTeOihh1BcXIyuXbsCAHJyctC+fXvs2LEDcXFxdi+SnE+iKhhAQ8O5KIoW99gRERE5E6t/tTdz5kzExsbiwoUL0Gg00Gg0KCgoQK9evfDHP/6xJWokJxSrVMDb0wPlVXXILa2UuhwiIqIWYfWK1LFjx5CRkYE2bdqYtrVp0wapqam488477VocOS9vTw/EdVAgM68cmnwtYtoHSF0SERGR3Vm9ItW1a1f8+uuvTbaXlJQ0ezNjcl+cJ0VERK7OoiCl1+tNj2XLluGFF17Ap59+igsXLuDChQv49NNPkZycjOXLl1t18vr6erz++uuIjo6Gr68vYmJisGTJEhiNRgBAXV0dUlJSEBcXB39/fyiVSkydOhWFhYVmx1m/fj2GDh2KoKAgCIIArVZ7y3NXVFQgOTkZarUavr6+GDBgAI4cOWK2jyiKWLx4MZRKJXx9fTF06FCcPHnSqvfozkx9Umw4JyIiF2XRpb3g4GCzZmFRFDFhwgTTNlEUATTch89gMFh88uXLl2PdunXYuHEjYmNjkZGRgaeeegoKhQJz585FVVUVNBoNFi5ciN69e6O8vBzJyckYO3YsMjIyTMepqqrCiBEjMGLECLzyyisWnXvmzJk4ceIENm/eDKVSiS1btuD+++/HTz/9hA4dOgAAVqxYgZUrV2LDhg3o0qUL3nzzTQwfPhynT59GYGCgxe/TXSWqGlakTv9agYrqOgT6eElcERERkX0JYmMKuomDBw9afMAhQ4ZYvO+DDz6IsLAwfPjhh6Ztjz76KPz8/LB58+ZmX3PkyBH069cPeXl5UKlUZn87cOAA7r33XpSXlyM4OPiG57169SoCAwPx73//G6NHjzZtj4+Px4MPPog333wToihCqVQiOTkZKSkpAICamhqEhYVh+fLlmDVr1i3fn16vh0KhgE6nQ1BQ0C33d0WDlu/DhfKr2PKHuzCoM+dJERGR47Pm+9uiFSlrwpE1Bg0ahHXr1iEnJwddunTBsWPHcOjQoZtOUtfpdBAE4aZB6Vbq6+thMBjg4+Njtt3X1xeHDh0CAOTm5qK4uBhJSUmmv8vlcgwZMgTff/99s0GqpqYGNTU1pud6vd7mGl1FoqoNLpRfRWZeOYMUERG5HIuC1PHjx9GzZ094eHjg+PHjN923V69eFp88JSUFOp0O3bp1g0wmg8FgQGpqKiZPntzs/tXV1ViwYAGmTJlyWys8gYGB6N+/P5YuXYru3bsjLCwM//jHP3D48GF07twZAFBcXAwACAsLM3ttWFgY8vLymj1uWloa3njjDZvrckWJqmDsOFbIhnMiInJJFgWp+Ph4FBcXIzQ0FPHx8RAEAc1dERQEwaoeqW3btmHLli3YunUrYmNjcfToUSQnJ0OpVDa5LU1dXR0mTZoEo9GINWvWWHyOG9m8eTNmzJiBDh06QCaTITExEVOmTIFGo2nynq53s+GSr7zyCubNm2d6rtfrERkZedu1OrPEa7/cy8ovh9EowsODgzmJiMh1WBSkcnNz0b59e9O/7eWll17CggULMGnSJABAXFwc8vLykJaWZhak6urqMGHCBOTm5mLfvn126Te64447cPDgQVRWVkKv1yMiIgITJ05EdHQ0ACA8PBxAw8pURESE6XUlJSVNVqkayeVyyOXy267NlXSPCIKPlwf01fU4V3oFnULZpE9ERK7DoiClVqub/fftqqqqgoeH+QQGmUxmGn8A/Baizpw5g/3796Nt27Z2Oz8A+Pv7w9/fH+Xl5di9ezdWrFgBAIiOjkZ4eDjS09ORkJAAAKitrcXBgwetHvPgzrxkHujVMRg/5l5GZl45gxQREbkUqwdybty4EV9++aXp+csvv4zg4GAMGDDghr1DNzJmzBikpqbiyy+/xPnz57F9+3asXLkS48aNA9DQFD5+/HhkZGTgk08+gcFgQHFxMYqLi1FbW2s6TnFxMY4ePYqzZ88CALKzs3H06FFcvnzZtM+wYcOwevVq0/Pdu3dj165dyM3NRXp6Ou6991507doVTz31FICGS3rJyclYtmwZtm/fjhMnTmD69Onw8/PDlClTrP3Y3FrjGARNnlbaQoiIiOxNtFKXLl3EvXv3iqIoit9//73o6+sr/u1vfxPHjBkjjhs3zqpj6fV6ce7cuaJKpRJ9fHzEmJgY8bXXXhNrampEURTF3NxcEUCzj/3795uOs2jRomb3+fjjj037qNVqcdGiRabn27ZtE2NiYkRvb28xPDxcfO6550StVmtWn9FoFBctWiSGh4eLcrlcHDx4sJidnW3x+9PpdCIAUafTWfW5uJrdJ4pEdcpO8f53D0hdChER0S1Z8/1t0Ryp6/n5+eHnn3+GSqVCSkoKioqKsGnTJpw8eRJDhw7FpUuX7BjznBvnSDUovVKDvm/uAQAcW5QEhS8HcxIRkeOy5vvb6kt7AQEBKCsrAwB88803uP/++wEAPj4+uHr1qg3lkqtrFyCHuq0fgIZf7xEREbkKq4PU8OHDMXPmTMycORM5OTmmyeAnT55EVFSUvesjF2Hqk8rXSlsIERGRHVkdpP7617+if//+uHTpEj777DPTr+gyMzNvOEiTqPEGxlyRIiIiV2J1jxRZjj1SvzlZqMPo9w8hUO6Jo4uSIONgTiIiclAt2iO1a9cu0/3ogIYVqvj4eEyZMgXl5VxtoOZ1DQuEn7cMFTX1OFNSIXU5REREdmF1kHrppZdMN+PNzs7Gn/70J4waNQrnzp0zuz0K0fU8ZR7o3TEYAOdJERGR67A6SOXm5qJHjx4AgM8++wwPPvggli1bhjVr1uDrr7+2e4HkOhLVwQDAGxgTEZHLsDpIeXt7o6qqCgCwZ88eJCUlAQBCQkJMK1VEzemjbvzlHoMUERG5BovutXe9QYMGYd68eRg4cCB+/PFHbNu2DQCQk5ODjh072r1Ach0JkQ1B6tylSpRX1qKNv7fEFREREd0eq1ekVq9eDU9PT3z66adYu3YtOnToAAD4+uuvMWLECLsXSK6jjb83Ytr5AwCyCrgqRUREzs/qFSmVSoWdO3c22f7//t//s0tB5NoSVG1wrrQSmjwt7usWJnU5REREt8WiIKXX601zFG7VB+Xu85Lo5vqo2+AzzQX2SRERkUuwKEi1adMGRUVFCA0NRXBwMASh6TBFURQhCAIMBoPdiyTX0fjLvaMFWtQbjPCUWX11mYiIyGFYFKT27duHkJAQAMD+/ftbtCBybZ1DAxEg98SVmnqc/rUCsUqF1CURERHZzKIgNWTIkGb/TWQtmYeA+MhgHDpbCk2+lkGKiIicmtXN5gBQXV2N48ePo6SkBEaj0exvY8eOtUth5LoS1W1w6GwpsvLK8eTdaqnLISIispnVQWrXrl2YOnUqSktLm/yNPVJkiURVMAAgkw3nRETk5Kzu9J0zZw4ee+wxFBUVwWg0mj0YosgSjYM588qqUHqlRuJqiIiIbGd1kCopKcG8efMQFsYZQGQbhZ8XOocGAACy8rXSFkNERHQbrA5S48ePx4EDB1qgFHIniaqGVanMPF7eIyIi52V1j9Tq1avx2GOP4bvvvkNcXBy8vLzM/v7CCy/YrThyXYnqYGzLKOBgTiIicmpWB6mtW7di9+7d8PX1xYEDB8yGcwqCwCBFFmlckTp+QYs6gxFeHMxJREROyOog9frrr2PJkiVYsGABPDz45Ue2uaN9AIJ8PKGvrsfPRRWI68h5UkRE5HysTkK1tbWYOHEiQxTdFg8PAQmmPqnLEldDRERkG6vT0LRp07Bt27aWqIXcTOPlPQ1/uUdERE7K6kt7BoMBK1aswO7du9GrV68mzeYrV660W3Hk2hpvYMyGcyIiclZWB6ns7GwkJCQAAE6cOGH2t+sbz4luJT4yGIIAXCi/ihJ9NUKDfKQuiYiIyCpWBSmDwYDFixcjLi4OISEhLVUTuYlAHy90DQvEz8UV0OSXY0TPCKlLIiIisopVPVIymQwPPPAAdDpdS9VDbiaBfVJEROTErG42j4uLw7lz51qiFnJDjTcw1nDCOREROSGrg1Rqairmz5+PnTt3oqioCHq93uxBZI0+6muDOS/qUFtvlLgaIiIi61jdbD5ixAgAwNixY82ay0VRhCAIMBgM9quOXF50O3+08fNCeVUdThbqTJf6iIiInIHVQWr//v0tUQe5KUFoGMy57+cSaPK1DFJERORUrA5SQ4YMaYk6yI0lqoKvBaly/AHRUpdDRERkMZvu8/Ldd9/hiSeewIABA3Dx4kUAwObNm3Ho0CG7FkfuIfFan1QWG86JiMjJWB2kPvvsMzzwwAPw9fWFRqNBTU0NAKCiogLLli2ze4Hk+np3DIaHABTqqlGkuyp1OURERBazOki9+eabWLduHT744AOz28MMGDAAGo3GrsWRe/CXe6JbeBAAQJOnlbYYIiIiK1gdpE6fPo3Bgwc32R4UFAStVmuPmsgN8b57RETkjKwOUhERETh79myT7YcOHUJMTIxdiiL30zhPikGKiIicidVBatasWZg7dy4OHz4MQRBQWFiITz75BPPnz8fs2bNbokZyA4nXxh6cuKhDdR1nkRERkXOwevzByy+/DJ1Oh3vvvRfV1dUYPHgw5HI55s+fjzlz5rREjeQGVCF+aOvvjbLKWpws1KGPmjfFJiIix2fT+IPU1FSUlpbixx9/xP/+9z9cunQJS5cutXdt5EYaB3MCbDgnIiLnYXGQqqqqwnPPPYcOHTogNDQUM2fORFRUFPr164eAgICWrJHcBPukiIjI2VgcpBYtWoQNGzZg9OjRmDRpEtLT0/Hss8+2ZG3kZhJVwQCAzLxyiKIobTFEREQWsLhH6vPPP8eHH36ISZMmAQCeeOIJDBw4EAaDATKZrMUKJPfRq2MwPD0ElFTU4KL2Kjq28ZO6JCIiopuyeEWqoKAA99xzj+l5v3794OnpicLCwhYpjNyPr7cM3SOuDebM10pbDBERkQUsDlIGgwHe3t5m2zw9PVFfX2/3osh9mfqkeN89IiJyAhYHKVEUMX36dDzyyCOmR3V1NZ555hmzbdaor6/H66+/jujoaPj6+iImJgZLliyB0WgEANTV1SElJQVxcXHw9/eHUqnE1KlTm6yCrV+/HkOHDkVQUBAEQbBowvqtzg0A06dPhyAIZo+7777bqvdI1km41ifFhnMiInIGFvdITZs2rcm2J5544rZOvnz5cqxbtw4bN25EbGwsMjIy8NRTT0GhUGDu3LmoqqqCRqPBwoUL0bt3b5SXlyM5ORljx45FRkaG6ThVVVUYMWIERowYgVdeecUu5240YsQIfPzxx6bnv1+VI/tqHMz5U6Ee1XUG+Hix/46IiByXxUHq+jBhLz/88AMeeughjB49GgAQFRWFf/zjH6aQpFAokJ6ebvaaVatWoV+/fsjPz4dKpQIAJCcnAwAOHDhgt3M3ksvlCA8Pt+iYNTU1qKmpMT3X6/UW10MNOrbxRWigHCUVNTh+QYd+0RzMSUREjsumgZz2MmjQIOzduxc5OTkAgGPHjuHQoUMYNWrUDV+j0+kgCAKCg4Nb5dwHDhxAaGgounTpgqeffholJSU3PGZaWhoUCoXpERkZeVs1uiNBEEyrUpnskyIiIgdn9S1i7CklJQU6nQ7dunWDTCaDwWBAamoqJk+e3Oz+1dXVWLBgAaZMmYKgoKAWP/fIkSPx2GOPQa1WIzc3FwsXLsR9992HzMxMyOXyJsd85ZVXMG/ePNNzvV7PMGWDRHUwdp0sZp8UERE5PEmD1LZt27BlyxZs3boVsbGxOHr0KJKTk6FUKpv0ZNXV1WHSpEkwGo1Ys2ZNq5x74sSJpv179uyJvn37Qq1W48svv2y2sV4ulzcbsMg6jStSWfkNgzkFQZC4IiIiouZJGqReeuklLFiwwDTkMy4uDnl5eUhLSzMLUnV1dZgwYQJyc3Oxb9++216Nsubc14uIiIBarcaZM2du+/x0Yz07KOAlE1B6pRYFl69C1ZaDOYmIyDFJ2iNVVVUFDw/zEmQymdkIgsYQdebMGezZswdt27ZttXP/XllZGQoKChAREWGXGqh5Pl4yxCoVAIDM/MsSV0NERHRjkgapMWPGIDU1FV9++SXOnz+P7du3Y+XKlRg3bhyAhllP48ePR0ZGBj755BMYDAYUFxejuLgYtbW1puMUFxfj6NGjOHv2LAAgOzsbR48exeXLv30JDxs2DKtXr7b43FeuXMH8+fPxww8/4Pz58zhw4ADGjBmDdu3amfahltN4eU+Tp5W2ECIiopsRJaTX68W5c+eKKpVK9PHxEWNiYsTXXntNrKmpEUVRFHNzc0UAzT72799vOs6iRYua3efjjz827aNWq8VFixZZfO6qqioxKSlJbN++vejl5SWqVCpx2rRpYn5+vsXvT6fTiQBEnU53W5+TO/rPsYuiOmWnOOov30pdChERuRlrvr8FURRFSRKcG9Dr9VAoFNDpdHbp63InRbqr6J+2DzIPAccXJcFfLmk7HxERuRFrvr8lvbRHdCMRCl9EKHxgMIo4dkErdTlERETNYpAih/XbGASttIUQERHdAIMUOSzTDYw54ZyIiBwUgxQ5rD7qa7/cuzaYk4iIyNEwSJHDilUq4O3pgfKqOuSWVkpdDhERURMMUuSwvD09ENehYTCnhn1SRETkgBikyKElNvZJ8QbGRETkgBikyKGZ+qTYcE5ERA6IQYocWuMIhNO/VqCiuk7iaoiIiMwxSJFDCw3yQYdgX4gicKxAJ3U5REREZhikyOElXjcGgYiIyJEwSJHD68OGcyIiclAMUuTwEq9rODcaOZiTiIgcB4MUObzuEUHw8fKAvroe50qvSF0OERGRCYMUOTwvmQd6dQgGAGjytJLWQkREdD0GKXIKbDgnIiJHxCBFTqFxwnkmB3MSEZEDYZAip9C4InWm5Ap0VzmYk4iIHAODFDmFdgFyqEL8AABHC7TSFkNERHQNgxQ5jcb77vHyHhEROQoGKXIajX1SWWw4JyIiB8EgRU4j4doNjI/mazmYk4iIHAKDFDmNbuGB8POWoaKmHmdKOJiTiIikxyBFTsNT5oHeHYMBsE+KiIgcA4MUOZVEdTAADuYkIiLHwCBFTiVRxQnnRETkOBikyKk0Npyfu1SJ8spaiashIiJ3xyBFTiXE3xsx7fwBAFkFXJUiIiJpMUiR02lcldLkaaUthIiI3B6DFDkdNpwTEZGjYJAip9N4q5hjBVrUG4wSV0NERO6MQYqcTufQQATIPVFZa8DpXyukLoeIiNwYgxQ5HZmHgPjIYACAJl8raS1EROTeGKTIKZluYMwJ50REJCEGKXJKiWoO5iQiIukxSJFTSohsCFLny6pQeqVG4mqIiMhdMUiRU1L4eaFTaAAAIIt9UkREJBEGKXJajX1SvLxHRERSYZAip9U4T0rDhnMiIpIIgxQ5rcRrt4o5dkGLOg7mJCIiCTBIkdO6o30Agnw8UV1nxM9FHMxJREStj0GKnJaHh4B4FccgEBGRdBikyKn1YZAiIiIJMUiRU0tUBwMAMtlwTkREEmCQIqcWHxkMQQAulF9FSUW11OUQEZGbYZAipxbo44UuoYEAAE2eVtpiiIjI7UgapOrr6/H6668jOjoavr6+iImJwZIlS2A0NvyUva6uDikpKYiLi4O/vz+USiWmTp2KwsJCs+OsX78eQ4cORVBQEARBgFarve1zA4Aoili8eDGUSiV8fX0xdOhQnDx50q6fAd2+xvvuZbFPioiIWpmkQWr58uVYt24dVq9ejVOnTmHFihV4++23sWrVKgBAVVUVNBoNFi5cCI1Gg88//xw5OTkYO3as2XGqqqowYsQIvPrqq3Y7NwCsWLECK1euxOrVq3HkyBGEh4dj+PDhqKjgT+0dSeOEc/ZJERFRaxNEURSlOvmDDz6IsLAwfPjhh6Ztjz76KPz8/LB58+ZmX3PkyBH069cPeXl5UKlUZn87cOAA7r33XpSXlyM4OPi2zi2KIpRKJZKTk5GSkgIAqKmpQVhYGJYvX45Zs2bd8v3p9XooFArodDoEBQXdcn+yzS+XrmDYuwfh7emBE4sfgLcnr1gTEZHtrPn+lvQbZ9CgQdi7dy9ycnIAAMeOHcOhQ4cwatSoG75Gp9NBEIRbBqXbPXdubi6Ki4uRlJRkeo1cLseQIUPw/fffN3vMmpoa6PV6swe1vJh2/gj280JtvRE/FfEzJyKi1uMp5clTUlKg0+nQrVs3yGQyGAwGpKamYvLkyc3uX11djQULFmDKlCm3vcJzq3MXFxcDAMLCwsxeFxYWhry8vGaPmZaWhjfeeOO26iLrCYKARFUb7Pu5BJl55YiPDJa6JCIichOSrkht27YNW7ZswdatW6HRaLBx40a888472LhxY5N96+rqMGnSJBiNRqxZs6bVzi0IgtlzURSbbGv0yiuvQKfTmR4FBQW3XSdZprFPioM5iYioNUm6IvXSSy9hwYIFmDRpEgAgLi4OeXl5SEtLw7Rp00z71dXVYcKECcjNzcW+ffvs0m90q3OHh4cDaFiZioiIML2upKSkySpVI7lcDrlcftu1kfUab2CcxYZzIiJqRZKuSFVVVcHDw7wEmUxmNoKgMUSdOXMGe/bsQdu2bVvl3NHR0QgPD0d6errp77W1tTh48CAGDBhglxrIfnpHBsNDAAp11SjSXZW6HCIichOSrkiNGTMGqampUKlUiI2NRVZWFlauXIkZM2YAaJj1NH78eGg0GuzcuRMGg8HUuxQSEgJvb28ADatGxcXFOHv2LAAgOzsbgYGBUKlUCAkJAQAMGzYM48aNw5w5cyw6tyAISE5OxrJly9C5c2d07twZy5Ytg5+fH6ZMmdKqnxPdmr/cE93Cg/BTkR6aPC1G9/KVuiQiInIHooT0er04d+5cUaVSiT4+PmJMTIz42muviTU1NaIoimJubq4IoNnH/v37TcdZtGhRs/t8/PHHpn3UarW4aNEii88tiqJoNBrFRYsWieHh4aJcLhcHDx4sZmdnW/z+dDqdCEDU6XQ2f0Zkude2HxfVKTvFJf85KXUpRETkxKz5/pZ0jpSr4xyp1vW55gLm/d8xJKiCsX32QKnLISIiJ+U0c6SI7KnPtVvFnLyoR3WdQeJqiIjIHTBIkctQhfihrb83ag1GnCzUSV0OERG5AQYpchmCICDh2hgETZ5W2mKIiMgtMEiRS0lUBwPgYE4iImodDFLkUvo0rkjll4O/oyAiopbGIEUupVfHYHh6CPhVX4OLWg7mJCKilsUgRS7F11uG7hENP1XV5GulLYaIiFwegxS5HNMNjHnfPSIiamEMUuRyEq/Nk8piwzkREbUwBilyOYnXGs5PFnIwJxERtSwGKXI5Hdv4on2gHPVGEccvcDAnERG1HAYpcjmCIPzWJ8XLe0RE1IIYpMglNd53jw3nRETUkhikyCUlcjAnERG1AgYpckk9OyjgJRNQeqUWBZc5mJOIiFoGgxS5JB8vGXooFQDYJ0VERC2HQYpc1vX33SMiImoJDFLkshLVwQCATDacExFRC/GUugCiltLYcP5zcQWqauvh583/ugOAKIqoN4qorTc2PAwN/7Pm2nNjKzTnC0KLnwICWuEkaJ33QkQ31tbfG6FBPpKdn98s5LKUwb4ID/JBsb4axwp06H9HW0nqEEXRFFaaCy41v9ve8G+D6d81v//b757XmB3P0PRY155ffx7+kJGIXMXsoXfg5RHdJDs/gxS5tD7qNvgyuwj7T5egjb+XWbiouUk4MQse1wWbmhuFk2Zef/2/HZnMQ4C3zAPeng0PWQsvsYho+RTXWkGReZRIev5yaaMMgxS5tARVML7MLsL6b89h/bfnpC4HAOAlMw8u3p4e157L4O3pAbns99s9IPc037/pPjKz59fvL//dsRpe/9v+Mg9emyIishWDFLm0UXER2PjDeZRX1v0WMrzMQ0WTsGL2N5lZIGl+n9+Flute0+T4Mg94MLgQEbkMBilyacpgX3z38n1Sl0FERC6K4w+IiIiIbMQgRURERGQjBikiIiIiGzFIEREREdmIQYqIiIjIRgxSRERERDZikCIiIiKyEYMUERERkY0YpIiIiIhsxCBFREREZCMGKSIiIiIbMUgRERER2YhBioiIiMhGDFJERERENvKUugBXJooiAECv10tcCREREVmq8Xu78Xv8ZhikWlBFRQUAIDIyUuJKiIiIyFoVFRVQKBQ33UcQLYlbZBOj0YjCwkIEBgZCEAS7Hluv1yMyMhIFBQUICgqy67FdDT8ry/Gzshw/K8vxs7IcPyvLteRnJYoiKioqoFQq4eFx8y4orki1IA8PD3Ts2LFFzxEUFMT/ZbMQPyvL8bOyHD8ry/Gzshw/K8u11Gd1q5WoRmw2JyIiIrIRgxQRERGRjRiknJRcLseiRYsgl8ulLsXh8bOyHD8ry/Gzshw/K8vxs7Kco3xWbDYnIiIishFXpIiIiIhsxCBFREREZCMGKSIiIiIbMUgRERER2YhBysl8++23GDNmDJRKJQRBwBdffCF1SQ4pLS0Nd955JwIDAxEaGoqHH34Yp0+flrosh7V27Vr06tXLNNiuf//++Prrr6Uuy+GlpaVBEAQkJydLXYpDWrx4MQRBMHuEh4dLXZbDunjxIp544gm0bdsWfn5+iI+PR2ZmptRlOZyoqKgm/70SBAHPPfecJPUwSDmZyspK9O7dG6tXr5a6FId28OBBPPfcc/jf//6H9PR01NfXIykpCZWVlVKX5pA6duyIt956CxkZGcjIyMB9992Hhx56CCdPnpS6NId15MgRrF+/Hr169ZK6FIcWGxuLoqIi0yM7O1vqkhxSeXk5Bg4cCC8vL3z99df46aef8O677yI4OFjq0hzOkSNHzP47lZ6eDgB47LHHJKmHt4hxMiNHjsTIkSOlLsPh7dq1y+z5xx9/jNDQUGRmZmLw4MESVeW4xowZY/Y8NTUVa9euxf/+9z/ExsZKVJXjunLlCh5//HF88MEHePPNN6Uux6F5enpyFcoCy5cvR2RkJD7++GPTtqioKOkKcmDt27c3e/7WW2/hjjvuwJAhQySphytS5BZ0Oh0AICQkROJKHJ/BYMA///lPVFZWon///lKX45Cee+45jB49Gvfff7/UpTi8M2fOQKlUIjo6GpMmTcK5c+ekLskh7dixA3379sVjjz2G0NBQJCQk4IMPPpC6LIdXW1uLLVu2YMaMGRAEQZIaGKTI5YmiiHnz5mHQoEHo2bOn1OU4rOzsbAQEBEAul+OZZ57B9u3b0aNHD6nLcjj//Oc/odFokJaWJnUpDu+uu+7Cpk2bsHv3bnzwwQcoLi7GgAEDUFZWJnVpDufcuXNYu3YtOnfujN27d+OZZ57BCy+8gE2bNkldmkP74osvoNVqMX36dMlq4KU9cnlz5szB8ePHcejQIalLcWhdu3bF0aNHodVq8dlnn2HatGk4ePAgw9R1CgoKMHfuXHzzzTfw8fGRuhyHd30bQlxcHPr374877rgDGzduxLx58ySszPEYjUb07dsXy5YtAwAkJCTg5MmTWLt2LaZOnSpxdY7rww8/xMiRI6FUKiWrgStS5NKef/557NixA/v370fHjh2lLseheXt7o1OnTujbty/S0tLQu3dv/OUvf5G6LIeSmZmJkpIS9OnTB56envD09MTBgwfx/vvvw9PTEwaDQeoSHZq/vz/i4uJw5swZqUtxOBEREU3+n5bu3bsjPz9fooocX15eHvbs2YOZM2dKWgdXpMgliaKI559/Htu3b8eBAwcQHR0tdUlORxRF1NTUSF2GQxk2bFiTX5099dRT6NatG1JSUiCTySSqzDnU1NTg1KlTuOeee6QuxeEMHDiwyYiWnJwcqNVqiSpyfI0/Iho9erSkdTBIOZkrV67g7Nmzpue5ubk4evQoQkJCoFKpJKzMsTz33HPYunUr/v3vfyMwMBDFxcUAAIVCAV9fX4mrczyvvvoqRo4cicjISFRUVOCf//wnDhw40OTXj+4uMDCwSZ+dv78/2rZty/67ZsyfPx9jxoyBSqVCSUkJ3nzzTej1ekybNk3q0hzOiy++iAEDBmDZsmWYMGECfvzxR6xfvx7r16+XujSHZDQa8fHHH2PatGnw9JQ4yojkVPbv3y8CaPKYNm2a1KU5lOY+IwDixx9/LHVpDmnGjBmiWq0Wvb29xfbt24vDhg0Tv/nmG6nLcgpDhgwR586dK3UZDmnixIliRESE6OXlJSqVSvGRRx4RT548KXVZDus///mP2LNnT1Eul4vdunUT169fL3VJDmv37t0iAPH06dNSlyIKoiiK0kQ4IiIiIufGZnMiIiIiGzFIEREREdmIQYqIiIjIRgxSRERERDZikCIiIiKyEYMUERERkY0YpIiIiIhsxCBFREREZCMGKSIiO4qKisJ7770ndRlE1EoYpIjIaU2fPh0PP/wwAGDo0KFITk5utXNv2LABwcHBTbYfOXIEf/zjH1utDiKSFm9aTER0ndraWnh7e9v8+vbt29uxGiJydFyRIiKnN336dBw8eBB/+ctfIAgCBEHA+fPnAQA//fQTRo0ahYCAAISFheHJJ59EaWmp6bVDhw7FnDlzMG/ePLRr1w7Dhw8HAKxcuRJxcXHw9/dHZGQkZs+ejStXrgAADhw4gKeeego6nc50vsWLFwNoemkvPz8fDz30EAICAhAUFIQJEybg119/Nf198eLFiI+Px+bNmxEVFQWFQoFJkyahoqLCtM+nn36KuLg4+Pr6om3btrj//vtRWVnZQp8mEVmDQYqInN5f/vIX9O/fH08//TSKiopQVFSEyMhIFBUVYciQIYiPj0dGRgZ27dqFX3/9FRMmTDB7/caNG+Hp6Yn//ve/+Nvf/gYA8PDwwPvvv48TJ05g48aN2LdvH15++WUAwIABA/Dee+8hKCjIdL758+c3qUsURTz88MO4fPkyDh48iPT0dPzyyy+YOHGi2X6//PILvvjiC+zcuRM7d+7EwYMH8dZbbwEAioqKMHnyZMyYMQOnTp3CgQMH8Mgjj4D3mydyDLy0R0ROT6FQwNvbG35+fggPDzdtX7t2LRITE7Fs2TLTto8++giRkZHIyclBly5dAACdOnXCihUrzI55fb9VdHQ0li5dimeffRZr1qyBt7c3FAoFBEEwO9/v7dmzB8ePH0dubi4iIyMBAJs3b0ZsbCyOHDmCO++8EwBgNBqxYcMGBAYGAgCefPJJ7N27F6mpqSgqKkJ9fT0eeeQRqNVqAEBcXNxtfFpEZE9ckSIil5WZmYn9+/cjICDA9OjWrRuAhlWgRn379m3y2v3792P48OHo0KEDAgMDMXXqVJSVlVl1Se3UqVOIjIw0hSgA6NGjB4KDg3Hq1CnTtqioKFOIAoCIiAiUlJQAAHr37o1hw4YhLi4Ojz32GD744AOUl5db/iEQUYtikCIil2U0GjFmzBgcPXrU7HHmzBkMHjzYtJ+/v7/Z6/Ly8jBq1Cj07NkTn332GTIzM/HXv/4VAFBXV2fx+UVRhCAIt9zu5eVl9ndBEGA0GgEAMpkM6enp+Prrr9GjRw+sWrUKXbt2RW5ursV1EFHLYZAiIpfg7e0Ng8Fgti0xMREnT55EVFQUOnXqZPb4fXi6XkZGBurr6/Huu+/i7rvvRpcuXVBYWHjL8/1ejx49kJ+fj4KCAtO2n376CTqdDt27d7f4vQmCgIEDB+KNN95AVlYWvL29sX37dotfT0Qth0GKiFxCVFQUDh8+jPPnz6O0tBRGoxHPPfccLl++jMmTJ+PHH3/EuXPn8M0332DGjBk3DUF33HEH6uvrsWrVKpw7dw6bN2/GunXrmpzvypUr2Lt3L0pLS1FVVdXkOPfffz969eqFxx9/HBqNBj/++COmTp2KIUOGNHs5sTmHDx/GsmXLkJGRgfz8fHz++ee4dOmSVUGMiFoOgxQRuYT58+dDJpOhR48eaN++PfLz86FUKvHf//4XBoMBDzzwAHr27Im5c+dCoVDAw+PG/+cvPj4eK1euxPLly9GzZ0988sknSEtLM9tnwIABeOaZZzBx4kS0b9++SbM60LCS9MUXX6BNmzYYPHgw7r//fsTExGDbtm0Wv6+goCB8++23GDVqFLp06YLXX38d7777LkaOHGn5h0NELUYQ+RtaIiIiIptwRYqIiIjIRgxSRERERDZikCIiIiKyEYMUERERkY0YpIiIiIhsxCBFREREZCMGKSIiIiIbMUgRERER2YhBioiIiMhGDFJERERENmKQIiIiIrLR/wdCShJHKWwuGgAAAABJRU5ErkJggg==",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Permissible current convergence plot\n",
    "I = I[:i+1]\n",
    "\n",
    "plt.figure()\n",
    "plt.plot(list(range(1,i+2)),I)\n",
    "plt.xlabel('Iterations')\n",
    "plt.ylabel('Permissible current [A]')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "b5519496-3d0e-462c-aa1b-7c3ac8e4d9bf",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "The permissible current rating of one conductor is equal to   820 A\n"
     ]
    }
   ],
   "source": [
    "# Permissible current rounding\n",
    "\n",
    "if I[i] < 200:\n",
    "    I_rating = int(I[i])\n",
    "elif I[i] < 500:\n",
    "    I_rating = int(I[i]/5)*5\n",
    "else:\n",
    "    I_rating = int(I[i]/10)*10\n",
    "        \n",
    "print(\"The permissible current rating of one conductor is equal to %5d A\" % I_rating)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.10.9"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}