cigreTB880_case05ecl.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/img4_cs05.png\" width=\"400\"/> <img src=\"./img/img6_numbFlat.png\" width=\"200\"/>\n",
    "</div>\n",
    "\n",
    "- 3: Outer cable lagging phase\n",
    "- 1: Central cable\n",
    "- 2: Outer cable leading phase\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\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_air = 25 #[°C] Maximum air temperature\n",
    "theta_a = 25 #[°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": "91069ec2-7e47-43fb-bb59-a4022616c418",
   "metadata": {},
   "source": [
    "**IEC 60287-2-1, Electric Cables - Calculation of the Current Rating - Part 2-1: Thermal Resistance - Calculation of Thermal Resistance. Table 2**\n",
    "\n",
    "<div>\n",
    "<img src=\"./img/img5_2-1_Table2.png\" width=\"500\"/>\n",
    "</div>"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "b3e0f670-0b8d-4659-9067-ef006947d4bf",
   "metadata": {},
   "outputs": [],
   "source": [
    "# 4.10.1.1 Installation\n",
    "\n",
    "Z = 0.62 # From IEC 60287-2-1 Table 2\n",
    "E = 1.95\n",
    "g = 0.25"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "23db7e44-6821-48f6-9d12-dbb0aedd217e",
   "metadata": {},
   "outputs": [],
   "source": [
    "# Unfilled trough external dimensions and number of cables\n",
    "\n",
    "tr_deep = 0.360 #[m]\n",
    "tr_wide = 0.550 #[m]\n",
    "p = tr_deep*2+tr_wide #[m] Part of the trough perimeter which is effective for heat dissipation\n",
    "\n",
    "NC = 3 # Amount of cables laid into the trough"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1c97cea0-4b59-48f4-8c7e-d0072f004986",
   "metadata": {},
   "source": [
    "## 4.10 CASE #0-5: Sub-case study with cables in an unfilled trough \n",
    "\n",
    "### 4.10.7 Variant with non-neglected eddy-current losses\n",
    "\n",
    "This option was adopted by the CIGRE WG B1.56, not to neglect the eddy-current losses when the metal sheaths are bonded at both ends of an electrical section (Guidance Point 6).\n"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "1feb518f-360f-4cdc-ab3f-dc7103b9db03",
   "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": 7,
   "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": 8,
   "id": "7fab0e9f-c7a2-4273-8192-1e2b1914dd64",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "X = 5.0403313985E-05 Ω/m\n",
      "Xm = 4.3551721806E-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.3. Reactance of sheath or screen per unit length of cable for two adjacent single-core cables\n",
    "\n",
    "# 4.8.1.1 Reactance\n",
    "X_m = 2*omega*1e-7*mt.log(2) #[Ohm/m] From IEC 60287-1-1 Section 2.3.3. Mutual reactance per unit length of cable between the sheath of an outer cable and the conductors of the other two\n",
    "\n",
    "print(\"X = %.10E \\u03A9/m\\nXm = %.10E \\u03A9/m\" % (X, X_m))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "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": 10,
   "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": 11,
   "id": "1fa6cb36-8d13-4b68-bb91-8cf052a9de24",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "T3 = 0.0541996092 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",
    "# The 1,6 multiplication factor for T3 of touching cables in trefoil formation only applies to buried cables (IEC 60287-2-1 Section 4.2.4.3).\n",
    "print(\"T3 = %.10f Km/W\" % T3)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "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": "code",
   "execution_count": 13,
   "id": "1101e7b9-91f3-4dfa-b3b1-11dd4c31dccd",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "h = 3.1327826671 W/(m^2 K^(5/4))\n"
     ]
    }
   ],
   "source": [
    "# 4.10.2.2 Heat dissipation coefficient \n",
    "\n",
    "D_e = dia[5]*1e-3 #[m] External diameter in meters\n",
    "h = Z/(D_e**g) + E #[W/(m^2 K^(5/4)] From IEC 60287-2-1 Section 4.2.1.1\n",
    "print(\"h = %.10f W/(m^2 K^(5/4))\" % h)"
   ]
  },
  {
   "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": 14,
   "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.8167187646\n",
      "Convergence reached after 13 iterations for cable surface temperature\n",
      "I = 794.2632582846 A\n",
      "Wc = 24.9323178820 W/m\n",
      "Ws = 20.3626918599 W/m\n",
      "θj = 83.4465796771 °C\n",
      "θs = 85.9224258467 °C\n",
      "θc = 96.4716495593 °C\n",
      "Iteration:  2\n",
      "R = 4.0188694612E-05 Ohm/m\n",
      "Rs = 2.1125616787E-04 Ohm/m\n",
      "λ1 = 0.7914738762\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 724.0436592571 A\n",
      "Wc = 21.0684899369 W/m\n",
      "Ws = 16.6751593970 W/m\n",
      "θj = 73.9885015023 °C\n",
      "θs = 76.0550668885 °C\n",
      "θc = 84.9819794084 °C\n",
      "Iteration:  3\n",
      "R = 3.9005193467E-05 Ohm/m\n",
      "Rs = 2.0461880159E-04 Ohm/m\n",
      "λ1 = 0.8356996787\n",
      "Convergence reached after 11 iterations for cable surface temperature\n",
      "I = 781.8852505039 A\n",
      "Wc = 23.8456122512 W/m\n",
      "Ws = 19.9277704957 W/m\n",
      "θj = 81.5197693003 °C\n",
      "θs = 83.9131438813 °C\n",
      "θc = 94.0060908825 °C\n",
      "Iteration:  4\n",
      "R = 3.9934352613E-05 Ohm/m\n",
      "Rs = 2.0990460659E-04 Ohm/m\n",
      "λ1 = 0.8004756964\n",
      "Convergence reached after 11 iterations for cable surface temperature\n",
      "I = 736.0536472937 A\n",
      "Wc = 21.6354327565 W/m\n",
      "Ws = 17.3186381037 W/m\n",
      "θj = 75.5042759427 °C\n",
      "θs = 77.6364457027 °C\n",
      "θc = 86.8014013484 °C\n",
      "Iteration:  5\n",
      "R = 3.9192303798E-05 Ohm/m\n",
      "Rs = 2.0568253003E-04 Ohm/m\n",
      "λ1 = 0.8284363646\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 772.7986031539 A\n",
      "Wc = 23.4063367887 W/m\n",
      "Ws = 19.3906605582 W/m\n",
      "θj = 80.3041246527 °C\n",
      "θs = 82.6445795265 °C\n",
      "θc = 92.5530872852 °C\n",
      "Iteration:  6\n",
      "R = 3.9784558623E-05 Ohm/m\n",
      "Rs = 2.0905129555E-04 Ohm/m\n",
      "λ1 = 0.8060152205\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 743.5702708188 A\n",
      "Wc = 21.9967530693 W/m\n",
      "Ws = 17.7297177761 W/m\n",
      "θj = 76.4709477550 °C\n",
      "θs = 78.6449812923 °C\n",
      "θc = 87.9616450358 °C\n",
      "Iteration:  7\n",
      "R = 3.9311683738E-05 Ohm/m\n",
      "Rs = 2.0636093042E-04 Ohm/m\n",
      "λ1 = 0.8238491925\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 766.9675808013 A\n",
      "Wc = 23.1246761446 W/m\n",
      "Ws = 19.0512457688 W/m\n",
      "θj = 79.5300334868 °C\n",
      "θs = 81.8368263149 °C\n",
      "θc = 91.6270727995 °C\n",
      "Iteration:  8\n",
      "R = 3.9689130347E-05 Ohm/m\n",
      "Rs = 2.0850795321E-04 Ohm/m\n",
      "λ1 = 0.8095715265\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 748.3318683238 A\n",
      "Wc = 22.2259362185 W/m\n",
      "Ws = 17.9934851128 W/m\n",
      "θj = 77.0873753109 °C\n",
      "θs = 79.2881265719 °C\n",
      "θc = 88.7010177855 °C\n",
      "Iteration:  9\n",
      "R = 3.9387783543E-05 Ohm/m\n",
      "Rs = 2.0679354778E-04 Ohm/m\n",
      "λ1 = 0.8209431095\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 763.2350286182 A\n",
      "Wc = 22.9444753065 W/m\n",
      "Ws = 18.8361089038 W/m\n",
      "θj = 79.0369661644 °C\n",
      "θs = 81.3223318434 °C\n",
      "θc = 91.0369171338 °C\n",
      "Iteration: 10\n",
      "R = 3.9628328329E-05 Ohm/m\n",
      "Rs = 2.0816187395E-04 Ohm/m\n",
      "λ1 = 0.8118485805\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 751.3549588581 A\n",
      "Wc = 22.3715495709 W/m\n",
      "Ws = 18.1623107619 W/m\n",
      "θj = 77.4803694495 °C\n",
      "θs = 79.6981631816 °C\n",
      "θc = 89.1721932903 °C\n",
      "Iteration: 11\n",
      "R = 3.9436288974E-05 Ohm/m\n",
      "Rs = 2.0706936253E-04 Ohm/m\n",
      "λ1 = 0.8190980802\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 760.8495217210 A\n",
      "Wc = 22.8293519878 W/m\n",
      "Ws = 18.6994783862 W/m\n",
      "θj = 78.7228439773 °C\n",
      "θs = 80.9945646968 °C\n",
      "θc = 90.6608129880 °C\n",
      "Iteration: 12\n",
      "R = 3.9589585569E-05 Ohm/m\n",
      "Rs = 2.0794139848E-04 Ohm/m\n",
      "λ1 = 0.8133040622\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 753.2769495309 A\n",
      "Wc = 22.4641666220 W/m\n",
      "Ws = 18.2701979677 W/m\n",
      "θj = 77.7308805676 °C\n",
      "θs = 79.9595415520 °C\n",
      "θc = 89.4724589198 °C\n",
      "Iteration: 13\n",
      "R = 3.9467203951E-05 Ohm/m\n",
      "Rs = 2.0724518101E-04 Ohm/m\n",
      "λ1 = 0.8179250966\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 759.3264868379 A\n",
      "Wc = 22.7558707497 W/m\n",
      "Ws = 18.6125977814 W/m\n",
      "θj = 78.5227007118 °C\n",
      "θs = 80.7857298820 °C\n",
      "θc = 90.4211254963 °C\n",
      "Iteration: 14\n",
      "R = 3.9564897691E-05 Ohm/m\n",
      "Rs = 2.0780092389E-04 Ohm/m\n",
      "λ1 = 0.8142333876\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 754.4999566366 A\n",
      "Wc = 22.5231166110 W/m\n",
      "Ws = 18.3390735365 W/m\n",
      "θj = 77.8905522089 °C\n",
      "θs = 80.1261412886 °C\n",
      "θc = 89.6638100761 °C\n",
      "Iteration: 15\n",
      "R = 3.9486906835E-05 Ohm/m\n",
      "Rs = 2.0735724579E-04 Ohm/m\n",
      "λ1 = 0.8171787190\n",
      "Convergence reached after 10 iterations for cable surface temperature\n",
      "I = 758.3547433294 A\n",
      "Wc = 22.7089958068 W/m\n",
      "Ws = 18.5573081024 W/m\n",
      "θj = 78.3951698657 °C\n",
      "θs = 80.6526617534 °C\n",
      "θc = 90.2683759156 °C\n",
      "Iteration: 16\n",
      "R = 3.9549165461E-05 Ohm/m\n",
      "Rs = 2.0771141444E-04 Ohm/m\n",
      "λ1 = 0.8148263506\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 755.2786114701 A\n",
      "Wc = 22.5606545771 W/m\n",
      "Ws = 18.3830158355 W/m\n",
      "θj = 77.9923182814 °C\n",
      "θs = 80.2323235596 °C\n",
      "θc = 89.7857534689 °C\n",
      "Iteration: 17\n",
      "R = 3.9499463651E-05 Ohm/m\n",
      "Rs = 2.0742867024E-04 Ohm/m\n",
      "λ1 = 0.8167035303\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 757.7350042140 A\n",
      "Wc = 22.6791043450 W/m\n",
      "Ws = 18.5221045822 W/m\n",
      "θj = 78.3139036868 °C\n",
      "θs = 80.5678674519 °C\n",
      "θc = 90.1710310416 °C\n",
      "Iteration: 18\n",
      "R = 3.9539139975E-05 Ohm/m\n",
      "Rs = 2.0765437680E-04 Ohm/m\n",
      "λ1 = 0.8152045278\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 755.7745318052 A\n",
      "Wc = 22.5845647093 W/m\n",
      "Ws = 18.4110394102 W/m\n",
      "θj = 78.0571760436 °C\n",
      "θs = 80.2999961085 °C\n",
      "θc = 89.8634652005 °C\n",
      "Iteration: 19\n",
      "R = 3.9507466087E-05 Ohm/m\n",
      "Rs = 2.0747419078E-04 Ohm/m\n",
      "λ1 = 0.8164008898\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 757.3398669469 A\n",
      "Wc = 22.6600474020 W/m\n",
      "Ws = 18.4996828626 W/m\n",
      "θj = 78.2621170858 °C\n",
      "θs = 80.5138327236 °C\n",
      "θc = 90.1089948462 °C\n",
      "Iteration: 20\n",
      "R = 3.9532751076E-05 Ohm/m\n",
      "Rs = 2.0761802987E-04 Ohm/m\n",
      "λ1 = 0.8154456521\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 756.0904500135 A\n",
      "Wc = 22.5997972579 W/m\n",
      "Ws = 18.4289064114 W/m\n",
      "θj = 78.0985103388 °C\n",
      "θs = 80.3431243863 °C\n",
      "θc = 89.9129891912 °C\n",
      "Iteration: 21\n",
      "R = 3.9512565972E-05 Ohm/m\n",
      "Rs = 2.0750320139E-04 Ohm/m\n",
      "λ1 = 0.8162080993\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 757.0879765965 A\n",
      "Wc = 22.6478996617 W/m\n",
      "Ws = 18.4853991369 W/m\n",
      "θj = 78.2291156191 °C\n",
      "θs = 80.4793986818 °C\n",
      "θc = 90.0694603146 °C\n",
      "Iteration: 22\n",
      "R = 3.9528679615E-05 Ohm/m\n",
      "Rs = 2.0759486750E-04 Ohm/m\n",
      "λ1 = 0.8155993642\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 756.2917289379 A\n",
      "Wc = 22.6095026662 W/m\n",
      "Ws = 18.4402960004 W/m\n",
      "θj = 78.1248525669 °C\n",
      "θs = 80.3706099550 °C\n",
      "θc = 89.9445497842 °C\n",
      "Iteration: 23\n",
      "R = 3.9515816065E-05 Ohm/m\n",
      "Rs = 2.0752168980E-04 Ohm/m\n",
      "λ1 = 0.8160852687\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 756.9274204892 A\n",
      "Wc = 22.6401568777 W/m\n",
      "Ws = 18.4762985100 W/m\n",
      "θj = 78.2080848892 °C\n",
      "θs = 80.4574550455 °C\n",
      "θc = 90.0442657041 °C\n",
      "Iteration: 24\n",
      "R = 3.9526084978E-05 Ohm/m\n",
      "Rs = 2.0758010692E-04 Ohm/m\n",
      "λ1 = 0.8156973417\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 756.4199799200 A\n",
      "Wc = 22.6156869208 W/m\n",
      "Ws = 18.4475557011 W/m\n",
      "θj = 78.1416402308 °C\n",
      "θs = 80.3881262760 °C\n",
      "θc = 89.9646626974 °C\n",
      "Iteration: 25\n",
      "R = 3.9517887299E-05 Ohm/m\n",
      "Rs = 2.0753347231E-04 Ohm/m\n",
      "λ1 = 0.8160070038\n",
      "Convergence reached after  9 iterations for cable surface temperature\n",
      "I = 756.8250885159 A\n",
      "Wc = 22.6352220395 W/m\n",
      "Ws = 18.4704997173 W/m\n",
      "θj = 78.1946826054 °C\n",
      "θs = 80.4434710031 °C\n",
      "θc = 90.0282096639 °C\n",
      "Iteration: 26\n",
      "R = 3.9524431477E-05 Ohm/m\n",
      "Rs = 2.0757070043E-04 Ohm/m\n",
      "λ1 = 0.8157597887\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.5017035723 A\n",
      "Wc = 22.6196276946 W/m\n",
      "Ws = 18.4521827092 W/m\n",
      "θj = 78.1523387912 °C\n",
      "θs = 80.3992892069 °C\n",
      "θc = 89.9774802468 °C\n",
      "Iteration: 27\n",
      "R = 3.9519207262E-05 Ohm/m\n",
      "Rs = 2.0754098116E-04 Ohm/m\n",
      "λ1 = 0.8159571322\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.7598692004 A\n",
      "Wc = 22.6320769558 W/m\n",
      "Ws = 18.4668046096 W/m\n",
      "θj = 78.1861416704 °C\n",
      "θs = 80.4345593325 °C\n",
      "θc = 90.0179774622 °C\n",
      "Iteration: 28\n",
      "R = 3.9523377738E-05 Ohm/m\n",
      "Rs = 2.0756470592E-04 Ohm/m\n",
      "λ1 = 0.8157995882\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.5537811225 A\n",
      "Wc = 22.6221389379 W/m\n",
      "Ws = 18.4551316308 W/m\n",
      "θj = 78.1591568181 °C\n",
      "θs = 80.4064031726 °C\n",
      "θc = 89.9856486120 °C\n",
      "Iteration: 29\n",
      "R = 3.9520048450E-05 Ohm/m\n",
      "Rs = 2.0754576643E-04 Ohm/m\n",
      "λ1 = 0.8159253522\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.7183041218 A\n",
      "Wc = 22.6300725713 W/m\n",
      "Ws = 18.4644499321 W/m\n",
      "θj = 78.1806987325 °C\n",
      "θs = 80.4288801351 °C\n",
      "θc = 90.0114566809 °C\n",
      "Iteration: 30\n",
      "R = 3.9522706212E-05 Ohm/m\n",
      "Rs = 2.0756088576E-04 Ohm/m\n",
      "λ1 = 0.8158249530\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.5869677664 A\n",
      "Wc = 22.6237392491 W/m\n",
      "Ws = 18.4570110098 W/m\n",
      "θj = 78.1635018302 °C\n",
      "θs = 80.4109367826 °C\n",
      "θc = 89.9908541470 °C\n",
      "Iteration: 31\n",
      "R = 3.9520584524E-05 Ohm/m\n",
      "Rs = 2.0754881600E-04 Ohm/m\n",
      "λ1 = 0.8159051002\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.6918146685 A\n",
      "Wc = 22.6287951824 W/m\n",
      "Ws = 18.4629494002 W/m\n",
      "θj = 78.1772300686 °C\n",
      "θs = 80.4252609090 °C\n",
      "θc = 90.0073011156 °C\n",
      "Iteration: 32\n",
      "R = 3.9522278263E-05 Ohm/m\n",
      "Rs = 2.0755845125E-04 Ohm/m\n",
      "λ1 = 0.8158411180\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.6081164090 A\n",
      "Wc = 22.6247590736 W/m\n",
      "Ws = 18.4582087372 W/m\n",
      "θj = 78.1662708278 °C\n",
      "θs = 80.4138259706 °C\n",
      "θc = 89.9941715302 °C\n",
      "Iteration: 33\n",
      "R = 3.9520926153E-05 Ohm/m\n",
      "Rs = 2.0755075944E-04 Ohm/m\n",
      "λ1 = 0.8158921943\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.6749331192 A\n",
      "Wc = 22.6279811139 W/m\n",
      "Ws = 18.4619931644 W/m\n",
      "θj = 78.1750195630 °C\n",
      "θs = 80.4229544536 °C\n",
      "θc = 90.0046528561 °C\n",
      "Iteration: 34\n",
      "R = 3.9522005540E-05 Ohm/m\n",
      "Rs = 2.0755689979E-04 Ohm/m\n",
      "λ1 = 0.8158514199\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.6215937969 A\n",
      "Wc = 22.6254089785 W/m\n",
      "Ws = 18.4589720399 W/m\n",
      "θj = 78.1680354581 °C\n",
      "θs = 80.4156671962 °C\n",
      "θc = 89.9962856323 °C\n",
      "Iteration: 35\n",
      "R = 3.9521143867E-05 Ohm/m\n",
      "Rs = 2.0755199796E-04 Ohm/m\n",
      "λ1 = 0.8158839698\n",
      "Convergence reached after  8 iterations for cable surface temperature\n",
      "I = 756.6641747025 A\n",
      "Wc = 22.6274623185 W/m\n",
      "Ws = 18.4613837838 W/m\n",
      "θj = 78.1736108538 °C\n",
      "θs = 80.4214845977 °C\n",
      "θc = 90.0029651728 °C\n",
      "Iteration: 36\n",
      "R = 3.9521831739E-05 Ohm/m\n",
      "Rs = 2.0755591108E-04 Ohm/m\n",
      "λ1 = 0.8158579851\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6301825782 A\n",
      "Wc = 22.6258231463 W/m\n",
      "Ws = 18.4594584836 W/m\n",
      "θj = 78.1691600231 °C\n",
      "θs = 80.4168405740 °C\n",
      "θc = 89.9976329074 °C\n",
      "Iteration: 37\n",
      "R = 3.9521282612E-05 Ohm/m\n",
      "Rs = 2.0755278724E-04 Ohm/m\n",
      "λ1 = 0.8158787286\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6573185202 A\n",
      "Wc = 22.6271316982 W/m\n",
      "Ws = 18.4609954415 W/m\n",
      "θj = 78.1727131122 °C\n",
      "θs = 80.4205478886 °C\n",
      "θc = 90.0018896457 °C\n",
      "Iteration: 38\n",
      "R = 3.9521720979E-05 Ohm/m\n",
      "Rs = 2.0755528100E-04 Ohm/m\n",
      "λ1 = 0.8158621690\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6356560012 A\n",
      "Wc = 22.6260870858 W/m\n",
      "Ws = 18.4597684869 W/m\n",
      "θj = 78.1698766866 °C\n",
      "θs = 80.4175883449 °C\n",
      "θc = 89.9984914990 °C\n",
      "Iteration: 39\n",
      "R = 3.9521371031E-05 Ohm/m\n",
      "Rs = 2.0755329023E-04 Ohm/m\n",
      "λ1 = 0.8158753885\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6529491883 A\n",
      "Wc = 22.6269209995 W/m\n",
      "Ws = 18.4607479602 W/m\n",
      "θj = 78.1721409997 °C\n",
      "θs = 80.4199509430 °C\n",
      "θc = 90.0012042336 °C\n",
      "Iteration: 40\n",
      "R = 3.9521650394E-05 Ohm/m\n",
      "Rs = 2.0755487946E-04 Ohm/m\n",
      "λ1 = 0.8158648354\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6391440900 A\n",
      "Wc = 22.6262552885 W/m\n",
      "Ws = 18.4599660468 W/m\n",
      "θj = 78.1703334022 °C\n",
      "θs = 80.4180648847 °C\n",
      "θc = 89.9990386623 °C\n",
      "Iteration: 41\n",
      "R = 3.9521427379E-05 Ohm/m\n",
      "Rs = 2.0755361078E-04 Ohm/m\n",
      "λ1 = 0.8158732599\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6501646905 A\n",
      "Wc = 22.6267867250 W/m\n",
      "Ws = 18.4605902459 W/m\n",
      "θj = 78.1717764040 °C\n",
      "θs = 80.4195705215 °C\n",
      "θc = 90.0007674342 °C\n",
      "Iteration: 42\n",
      "R = 3.9521605411E-05 Ohm/m\n",
      "Rs = 2.0755462356E-04 Ohm/m\n",
      "λ1 = 0.8158665346\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6413669736 A\n",
      "Wc = 22.6263624806 W/m\n",
      "Ws = 18.4600919482 W/m\n",
      "θj = 78.1706244580 °C\n",
      "θs = 80.4183685741 °C\n",
      "θc = 89.9993873586 °C\n",
      "Iteration: 43\n",
      "R = 3.9521463289E-05 Ohm/m\n",
      "Rs = 2.0755381506E-04 Ohm/m\n",
      "λ1 = 0.8158719034\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6483901816 A\n",
      "Wc = 22.6267011544 W/m\n",
      "Ws = 18.4604897380 W/m\n",
      "θj = 78.1715440544 °C\n",
      "θs = 80.4193280865 °C\n",
      "θc = 90.0004890705 °C\n",
      "Iteration: 44\n",
      "R = 3.9521576745E-05 Ohm/m\n",
      "Rs = 2.0755446049E-04 Ohm/m\n",
      "λ1 = 0.8158676175\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6427835706 A\n",
      "Wc = 22.6264307918 W/m\n",
      "Ws = 18.4601721829 W/m\n",
      "θj = 78.1708099421 °C\n",
      "θs = 80.4185621094 °C\n",
      "θc = 89.9996095758 °C\n",
      "Iteration: 45\n",
      "R = 3.9521486173E-05 Ohm/m\n",
      "Rs = 2.0755394524E-04 Ohm/m\n",
      "λ1 = 0.8158710389\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6472593208 A\n",
      "Wc = 22.6266466219 W/m\n",
      "Ws = 18.4604256864 W/m\n",
      "θj = 78.1713959826 °C\n",
      "θs = 80.4191735876 °C\n",
      "θc = 90.0003116749 °C\n",
      "Iteration: 46\n",
      "R = 3.9521558476E-05 Ohm/m\n",
      "Rs = 2.0755435656E-04 Ohm/m\n",
      "λ1 = 0.8158683076\n",
      "Convergence reached after  7 iterations for cable surface temperature\n",
      "I = 756.6436863387 A\n",
      "Wc = 22.6264743252 W/m\n",
      "Ws = 18.4602233149 W/m\n",
      "θj = 78.1709281475 °C\n",
      "θs = 80.4186854456 °C\n",
      "θc = 89.9997511904 °C\n",
      "Iteration: 47\n",
      "R = 3.9521500757E-05 Ohm/m\n",
      "Rs = 2.0755402821E-04 Ohm/m\n",
      "λ1 = 0.8158704880\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6465386451 A\n",
      "Wc = 22.6266118694 W/m\n",
      "Ws = 18.4603848677 W/m\n",
      "θj = 78.1713016194 °C\n",
      "θs = 80.4190751285 °C\n",
      "θc = 90.0001986242 °C\n",
      "Iteration: 48\n",
      "R = 3.9521546834E-05 Ohm/m\n",
      "Rs = 2.0755429033E-04 Ohm/m\n",
      "λ1 = 0.8158687474\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6442616545 A\n",
      "Wc = 22.6265020681 W/m\n",
      "Ws = 18.4602559003 W/m\n",
      "θj = 78.1710034774 °C\n",
      "θs = 80.4187640453 °C\n",
      "θc = 89.9998414386 °C\n",
      "Iteration: 49\n",
      "R = 3.9521510051E-05 Ohm/m\n",
      "Rs = 2.0755408108E-04 Ohm/m\n",
      "λ1 = 0.8158701369\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6460793725 A\n",
      "Wc = 22.6265897223 W/m\n",
      "Ws = 18.4603588547 W/m\n",
      "θj = 78.1712414837 °C\n",
      "θs = 80.4190123824 °C\n",
      "θc = 90.0001265793 °C\n",
      "Iteration: 50\n",
      "R = 3.9521539415E-05 Ohm/m\n",
      "Rs = 2.0755424812E-04 Ohm/m\n",
      "λ1 = 0.8158690277\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6446282917 A\n",
      "Wc = 22.6265197481 W/m\n",
      "Ws = 18.4602766664 W/m\n",
      "θj = 78.1710514837 °C\n",
      "θs = 80.4188141353 °C\n",
      "θc = 89.9998989520 °C\n",
      "Iteration: 51\n",
      "R = 3.9521515973E-05 Ohm/m\n",
      "Rs = 2.0755411477E-04 Ohm/m\n",
      "λ1 = 0.8158699132\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6457866871 A\n",
      "Wc = 22.6265756084 W/m\n",
      "Ws = 18.4603422772 W/m\n",
      "θj = 78.1712031603 °C\n",
      "θs = 80.4189723956 °C\n",
      "θc = 90.0000806665 °C\n",
      "Iteration: 52\n",
      "R = 3.9521534687E-05 Ohm/m\n",
      "Rs = 2.0755422123E-04 Ohm/m\n",
      "λ1 = 0.8158692063\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6448619422 A\n",
      "Wc = 22.6265310153 W/m\n",
      "Ws = 18.4602899002 W/m\n",
      "θj = 78.1710820771 °C\n",
      "θs = 80.4188460567 °C\n",
      "θc = 89.9999356041 °C\n",
      "Iteration: 53\n",
      "R = 3.9521519748E-05 Ohm/m\n",
      "Rs = 2.0755413624E-04 Ohm/m\n",
      "λ1 = 0.8158697706\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6456001645 A\n",
      "Wc = 22.6265666139 W/m\n",
      "Ws = 18.4603317126 W/m\n",
      "θj = 78.1711787377 °C\n",
      "θs = 80.4189469129 °C\n",
      "θc = 90.0000514071 °C\n",
      "Iteration: 54\n",
      "R = 3.9521531674E-05 Ohm/m\n",
      "Rs = 2.0755420409E-04 Ohm/m\n",
      "λ1 = 0.8158693201\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6450108429 A\n",
      "Wc = 22.6265381956 W/m\n",
      "Ws = 18.4602983338 W/m\n",
      "θj = 78.1711015737 °C\n",
      "θs = 80.4188663996 °C\n",
      "θc = 89.9999589618 °C\n",
      "Iteration: 55\n",
      "R = 3.9521522153E-05 Ohm/m\n",
      "Rs = 2.0755414993E-04 Ohm/m\n",
      "λ1 = 0.8158696797\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6454812974 A\n",
      "Wc = 22.6265608819 W/m\n",
      "Ws = 18.4603249801 W/m\n",
      "θj = 78.1711631736 °C\n",
      "θs = 80.4189306732 °C\n",
      "θc = 90.0000327607 °C\n",
      "Iteration: 56\n",
      "R = 3.9521529753E-05 Ohm/m\n",
      "Rs = 2.0755419316E-04 Ohm/m\n",
      "λ1 = 0.8158693926\n",
      "Convergence reached after  6 iterations for cable surface temperature\n",
      "I = 756.6451057344 A\n",
      "Wc = 22.6265427714 W/m\n",
      "Ws = 18.4603037084 W/m\n",
      "θj = 78.1711139985 °C\n",
      "θs = 80.4188793637 °C\n",
      "θc = 89.9999738472 °C\n",
      "Iteration: 57\n",
      "R = 3.9521523686E-05 Ohm/m\n",
      "Rs = 2.0755415865E-04 Ohm/m\n",
      "λ1 = 0.8158696218\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6454055457 A\n",
      "Wc = 22.6265572290 W/m\n",
      "Ws = 18.4603206895 W/m\n",
      "θj = 78.1711532549 °C\n",
      "θs = 80.4189203240 °C\n",
      "θc = 90.0000208778 °C\n",
      "Iteration: 58\n",
      "R = 3.9521528530E-05 Ohm/m\n",
      "Rs = 2.0755418620E-04 Ohm/m\n",
      "λ1 = 0.8158694389\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6451662069 A\n",
      "Wc = 22.6265456875 W/m\n",
      "Ws = 18.4603071335 W/m\n",
      "θj = 78.1711219166 °C\n",
      "θs = 80.4188876254 °C\n",
      "θc = 89.9999833333 °C\n",
      "Iteration: 59\n",
      "R = 3.9521524663E-05 Ohm/m\n",
      "Rs = 2.0755416421E-04 Ohm/m\n",
      "λ1 = 0.8158695849\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6453572706 A\n",
      "Wc = 22.6265549010 W/m\n",
      "Ws = 18.4603179553 W/m\n",
      "θj = 78.1711469339 °C\n",
      "θs = 80.4189137286 °C\n",
      "θc = 90.0000133050 °C\n",
      "Iteration: 60\n",
      "R = 3.9521527750E-05 Ohm/m\n",
      "Rs = 2.0755418176E-04 Ohm/m\n",
      "λ1 = 0.8158694683\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6452047448 A\n",
      "Wc = 22.6265475459 W/m\n",
      "Ws = 18.4603093163 W/m\n",
      "θj = 78.1711269626 °C\n",
      "θs = 80.4188928905 °C\n",
      "θc = 89.9999893787 °C\n",
      "Iteration: 61\n",
      "R = 3.9521525286E-05 Ohm/m\n",
      "Rs = 2.0755416775E-04 Ohm/m\n",
      "λ1 = 0.8158695614\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6453265059 A\n",
      "Wc = 22.6265534175 W/m\n",
      "Ws = 18.4603162128 W/m\n",
      "θj = 78.1711429056 °C\n",
      "θs = 80.4189095255 °C\n",
      "θc = 90.0000084790 °C\n",
      "Iteration: 62\n",
      "R = 3.9521527253E-05 Ohm/m\n",
      "Rs = 2.0755417894E-04 Ohm/m\n",
      "λ1 = 0.8158694871\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6452293042 A\n",
      "Wc = 22.6265487302 W/m\n",
      "Ws = 18.4603107073 W/m\n",
      "θj = 78.1711301784 °C\n",
      "θs = 80.4188962458 °C\n",
      "θc = 89.9999932312 °C\n",
      "Iteration: 63\n",
      "R = 3.9521525682E-05 Ohm/m\n",
      "Rs = 2.0755417000E-04 Ohm/m\n",
      "λ1 = 0.8158695464\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6453069002 A\n",
      "Wc = 22.6265524721 W/m\n",
      "Ws = 18.4603151023 W/m\n",
      "θj = 78.1711403385 °C\n",
      "θs = 80.4189068470 °C\n",
      "θc = 90.0000054035 °C\n",
      "Iteration: 64\n",
      "R = 3.9521526936E-05 Ohm/m\n",
      "Rs = 2.0755417714E-04 Ohm/m\n",
      "λ1 = 0.8158694991\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6452449554 A\n",
      "Wc = 22.6265494850 W/m\n",
      "Ws = 18.4603115938 W/m\n",
      "θj = 78.1711322277 °C\n",
      "θs = 80.4188983841 °C\n",
      "θc = 89.9999956864 °C\n",
      "Iteration: 65\n",
      "R = 3.9521525935E-05 Ohm/m\n",
      "Rs = 2.0755417144E-04 Ohm/m\n",
      "λ1 = 0.8158695369\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6452944058 A\n",
      "Wc = 22.6265518696 W/m\n",
      "Ws = 18.4603143946 W/m\n",
      "θj = 78.1711387026 °C\n",
      "θs = 80.4189051400 °C\n",
      "θc = 90.0000034435 °C\n",
      "Iteration: 66\n",
      "R = 3.9521526734E-05 Ohm/m\n",
      "Rs = 2.0755417599E-04 Ohm/m\n",
      "λ1 = 0.8158695067\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6452549297 A\n",
      "Wc = 22.6265499659 W/m\n",
      "Ws = 18.4603121587 W/m\n",
      "θj = 78.1711335337 °C\n",
      "θs = 80.4188997468 °C\n",
      "θc = 89.9999972510 °C\n",
      "Iteration: 67\n",
      "R = 3.9521526096E-05 Ohm/m\n",
      "Rs = 2.0755417236E-04 Ohm/m\n",
      "λ1 = 0.8158695308\n",
      "Convergence reached after  5 iterations for cable surface temperature\n",
      "I = 756.6452864434 A\n",
      "Wc = 22.6265514856 W/m\n",
      "Ws = 18.4603139437 W/m\n",
      "θj = 78.1711376600 °C\n",
      "θs = 80.4189040522 °C\n",
      "θc = 90.0000021945 °C\n",
      "Iteration: 68\n",
      "R = 3.9521526606E-05 Ohm/m\n",
      "Rs = 2.0755417526E-04 Ohm/m\n",
      "λ1 = 0.8158695115\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452612860 A\n",
      "Wc = 22.6265502725 W/m\n",
      "Ws = 18.4603125188 W/m\n",
      "θj = 78.1711343660 °C\n",
      "θs = 80.4189006152 °C\n",
      "θc = 89.9999982481 °C\n",
      "Iteration: 69\n",
      "R = 3.9521526199E-05 Ohm/m\n",
      "Rs = 2.0755417294E-04 Ohm/m\n",
      "λ1 = 0.8158695269\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452813691 A\n",
      "Wc = 22.6265512409 W/m\n",
      "Ws = 18.4603136563 W/m\n",
      "θj = 78.1711369956 °C\n",
      "θs = 80.4189033589 °C\n",
      "θc = 90.0000013985 °C\n",
      "Iteration: 70\n",
      "R = 3.9521526524E-05 Ohm/m\n",
      "Rs = 2.0755417479E-04 Ohm/m\n",
      "λ1 = 0.8158695146\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452653368 A\n",
      "Wc = 22.6265504678 W/m\n",
      "Ws = 18.4603127482 W/m\n",
      "θj = 78.1711348964 °C\n",
      "θs = 80.4189011686 °C\n",
      "θc = 89.9999988836 °C\n",
      "Iteration: 71\n",
      "R = 3.9521526265E-05 Ohm/m\n",
      "Rs = 2.0755417332E-04 Ohm/m\n",
      "λ1 = 0.8158695244\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452781354 A\n",
      "Wc = 22.6265510850 W/m\n",
      "Ws = 18.4603134731 W/m\n",
      "θj = 78.1711365722 °C\n",
      "θs = 80.4189029171 °C\n",
      "θc = 90.0000008912 °C\n",
      "Iteration: 72\n",
      "R = 3.9521526471E-05 Ohm/m\n",
      "Rs = 2.0755417449E-04 Ohm/m\n",
      "λ1 = 0.8158695166\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452679183 A\n",
      "Wc = 22.6265505923 W/m\n",
      "Ws = 18.4603128944 W/m\n",
      "θj = 78.1711352344 °C\n",
      "θs = 80.4189015213 °C\n",
      "θc = 89.9999992885 °C\n",
      "Iteration: 73\n",
      "R = 3.9521526306E-05 Ohm/m\n",
      "Rs = 2.0755417355E-04 Ohm/m\n",
      "λ1 = 0.8158695229\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452760746 A\n",
      "Wc = 22.6265509856 W/m\n",
      "Ws = 18.4603133564 W/m\n",
      "θj = 78.1711363023 °C\n",
      "θs = 80.4189026356 °C\n",
      "θc = 90.0000005680 °C\n",
      "Iteration: 74\n",
      "R = 3.9521526438E-05 Ohm/m\n",
      "Rs = 2.0755417430E-04 Ohm/m\n",
      "λ1 = 0.8158695179\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452695635 A\n",
      "Wc = 22.6265506716 W/m\n",
      "Ws = 18.4603129876 W/m\n",
      "θj = 78.1711354498 °C\n",
      "θs = 80.4189017460 °C\n",
      "θc = 89.9999995466 °C\n",
      "Iteration: 75\n",
      "R = 3.9521526333E-05 Ohm/m\n",
      "Rs = 2.0755417370E-04 Ohm/m\n",
      "λ1 = 0.8158695218\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452747613 A\n",
      "Wc = 22.6265509223 W/m\n",
      "Ws = 18.4603132820 W/m\n",
      "θj = 78.1711361304 °C\n",
      "θs = 80.4189024562 °C\n",
      "θc = 90.0000003620 °C\n",
      "Iteration: 76\n",
      "R = 3.9521526417E-05 Ohm/m\n",
      "Rs = 2.0755417418E-04 Ohm/m\n",
      "λ1 = 0.8158695187\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452706119 A\n",
      "Wc = 22.6265507222 W/m\n",
      "Ws = 18.4603130470 W/m\n",
      "θj = 78.1711355871 °C\n",
      "θs = 80.4189018893 °C\n",
      "θc = 89.9999997111 °C\n",
      "Iteration: 77\n",
      "R = 3.9521526350E-05 Ohm/m\n",
      "Rs = 2.0755417380E-04 Ohm/m\n",
      "λ1 = 0.8158695212\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452739243 A\n",
      "Wc = 22.6265508819 W/m\n",
      "Ws = 18.4603132346 W/m\n",
      "θj = 78.1711360208 °C\n",
      "θs = 80.4189023418 °C\n",
      "θc = 90.0000002307 °C\n",
      "Iteration: 78\n",
      "R = 3.9521526403E-05 Ohm/m\n",
      "Rs = 2.0755417410E-04 Ohm/m\n",
      "λ1 = 0.8158695192\n",
      "Convergence reached after  4 iterations for cable surface temperature\n",
      "I = 756.6452712800 A\n",
      "Wc = 22.6265507544 W/m\n",
      "Ws = 18.4603130848 W/m\n",
      "θj = 78.1711356745 °C\n",
      "θs = 80.4189019805 °C\n",
      "θc = 89.9999998159 °C\n",
      "Iteration: 79\n",
      "R = 3.9521526361E-05 Ohm/m\n",
      "Rs = 2.0755417386E-04 Ohm/m\n",
      "λ1 = 0.8158695208\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452733910 A\n",
      "Wc = 22.6265508562 W/m\n",
      "Ws = 18.4603132044 W/m\n",
      "θj = 78.1711359509 °C\n",
      "θs = 80.4189022689 °C\n",
      "θc = 90.0000001470 °C\n",
      "Iteration: 80\n",
      "R = 3.9521526395E-05 Ohm/m\n",
      "Rs = 2.0755417406E-04 Ohm/m\n",
      "λ1 = 0.8158695195\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452717058 A\n",
      "Wc = 22.6265507749 W/m\n",
      "Ws = 18.4603131089 W/m\n",
      "θj = 78.1711357303 °C\n",
      "θs = 80.4189020387 °C\n",
      "θc = 89.9999998826 °C\n",
      "Iteration: 81\n",
      "R = 3.9521526367E-05 Ohm/m\n",
      "Rs = 2.0755417390E-04 Ohm/m\n",
      "λ1 = 0.8158695205\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452730511 A\n",
      "Wc = 22.6265508398 W/m\n",
      "Ws = 18.4603131851 W/m\n",
      "θj = 78.1711359064 °C\n",
      "θs = 80.4189022225 °C\n",
      "θc = 90.0000000937 °C\n",
      "Iteration: 82\n",
      "R = 3.9521526389E-05 Ohm/m\n",
      "Rs = 2.0755417402E-04 Ohm/m\n",
      "λ1 = 0.8158695197\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452719771 A\n",
      "Wc = 22.6265507880 W/m\n",
      "Ws = 18.4603131243 W/m\n",
      "θj = 78.1711357658 °C\n",
      "θs = 80.4189020758 °C\n",
      "θc = 89.9999999252 °C\n",
      "Iteration: 83\n",
      "R = 3.9521526372E-05 Ohm/m\n",
      "Rs = 2.0755417393E-04 Ohm/m\n",
      "λ1 = 0.8158695204\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452728345 A\n",
      "Wc = 22.6265508293 W/m\n",
      "Ws = 18.4603131729 W/m\n",
      "θj = 78.1711358781 °C\n",
      "θs = 80.4189021929 °C\n",
      "θc = 90.0000000597 °C\n",
      "Iteration: 84\n",
      "R = 3.9521526386E-05 Ohm/m\n",
      "Rs = 2.0755417400E-04 Ohm/m\n",
      "λ1 = 0.8158695199\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452721500 A\n",
      "Wc = 22.6265507963 W/m\n",
      "Ws = 18.4603131341 W/m\n",
      "θj = 78.1711357885 °C\n",
      "θs = 80.4189020994 °C\n",
      "θc = 89.9999999523 °C\n",
      "Iteration: 85\n",
      "R = 3.9521526375E-05 Ohm/m\n",
      "Rs = 2.0755417394E-04 Ohm/m\n",
      "λ1 = 0.8158695203\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452726964 A\n",
      "Wc = 22.6265508227 W/m\n",
      "Ws = 18.4603131650 W/m\n",
      "θj = 78.1711358600 °C\n",
      "θs = 80.4189021741 °C\n",
      "θc = 90.0000000380 °C\n",
      "Iteration: 86\n",
      "R = 3.9521526383E-05 Ohm/m\n",
      "Rs = 2.0755417399E-04 Ohm/m\n",
      "λ1 = 0.8158695199\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452722602 A\n",
      "Wc = 22.6265508017 W/m\n",
      "Ws = 18.4603131403 W/m\n",
      "θj = 78.1711358029 °C\n",
      "θs = 80.4189021145 °C\n",
      "θc = 89.9999999696 °C\n",
      "Iteration: 87\n",
      "R = 3.9521526376E-05 Ohm/m\n",
      "Rs = 2.0755417395E-04 Ohm/m\n",
      "λ1 = 0.8158695202\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452726084 A\n",
      "Wc = 22.6265508184 W/m\n",
      "Ws = 18.4603131601 W/m\n",
      "θj = 78.1711358485 °C\n",
      "θs = 80.4189021620 °C\n",
      "θc = 90.0000000242 °C\n",
      "Iteration: 88\n",
      "R = 3.9521526382E-05 Ohm/m\n",
      "Rs = 2.0755417398E-04 Ohm/m\n",
      "λ1 = 0.8158695200\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452723305 A\n",
      "Wc = 22.6265508050 W/m\n",
      "Ws = 18.4603131443 W/m\n",
      "θj = 78.1711358121 °C\n",
      "θs = 80.4189021241 °C\n",
      "θc = 89.9999999806 °C\n",
      "Iteration: 89\n",
      "R = 3.9521526378E-05 Ohm/m\n",
      "Rs = 2.0755417396E-04 Ohm/m\n",
      "λ1 = 0.8158695202\n",
      "Convergence reached after  3 iterations for cable surface temperature\n",
      "I = 756.6452725524 A\n",
      "Wc = 22.6265508157 W/m\n",
      "Ws = 18.4603131569 W/m\n",
      "θj = 78.1711358411 °C\n",
      "θs = 80.4189021544 °C\n",
      "θc = 90.0000000155 °C\n",
      "Iteration: 90\n",
      "R = 3.9521526381E-05 Ohm/m\n",
      "Rs = 2.0755417398E-04 Ohm/m\n",
      "λ1 = 0.8158695200\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452723752 A\n",
      "Wc = 22.6265508072 W/m\n",
      "Ws = 18.4603131468 W/m\n",
      "θj = 78.1711358179 °C\n",
      "θs = 80.4189021302 °C\n",
      "θc = 89.9999999877 °C\n",
      "Iteration: 91\n",
      "R = 3.9521526378E-05 Ohm/m\n",
      "Rs = 2.0755417396E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452725166 A\n",
      "Wc = 22.6265508140 W/m\n",
      "Ws = 18.4603131549 W/m\n",
      "θj = 78.1711358365 °C\n",
      "θs = 80.4189021495 °C\n",
      "θc = 90.0000000098 °C\n",
      "Iteration: 92\n",
      "R = 3.9521526381E-05 Ohm/m\n",
      "Rs = 2.0755417398E-04 Ohm/m\n",
      "λ1 = 0.8158695200\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724038 A\n",
      "Wc = 22.6265508086 W/m\n",
      "Ws = 18.4603131485 W/m\n",
      "θj = 78.1711358217 °C\n",
      "θs = 80.4189021341 °C\n",
      "θc = 89.9999999921 °C\n",
      "Iteration: 93\n",
      "R = 3.9521526379E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724939 A\n",
      "Wc = 22.6265508129 W/m\n",
      "Ws = 18.4603131536 W/m\n",
      "θj = 78.1711358335 °C\n",
      "θs = 80.4189021464 °C\n",
      "θc = 90.0000000063 °C\n",
      "Iteration: 94\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724219 A\n",
      "Wc = 22.6265508095 W/m\n",
      "Ws = 18.4603131495 W/m\n",
      "θj = 78.1711358241 °C\n",
      "θs = 80.4189021366 °C\n",
      "θc = 89.9999999950 °C\n",
      "Iteration: 95\n",
      "R = 3.9521526379E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724793 A\n",
      "Wc = 22.6265508122 W/m\n",
      "Ws = 18.4603131527 W/m\n",
      "θj = 78.1711358316 °C\n",
      "θs = 80.4189021444 °C\n",
      "θc = 90.0000000040 °C\n",
      "Iteration: 96\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724335 A\n",
      "Wc = 22.6265508100 W/m\n",
      "Ws = 18.4603131502 W/m\n",
      "θj = 78.1711358256 °C\n",
      "θs = 80.4189021381 °C\n",
      "θc = 89.9999999968 °C\n",
      "Iteration: 97\n",
      "R = 3.9521526379E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724701 A\n",
      "Wc = 22.6265508118 W/m\n",
      "Ws = 18.4603131522 W/m\n",
      "θj = 78.1711358304 °C\n",
      "θs = 80.4189021431 °C\n",
      "θc = 90.0000000025 °C\n",
      "Iteration: 98\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724409 A\n",
      "Wc = 22.6265508104 W/m\n",
      "Ws = 18.4603131506 W/m\n",
      "θj = 78.1711358265 °C\n",
      "θs = 80.4189021391 °C\n",
      "θc = 89.9999999980 °C\n",
      "Iteration: 99\n",
      "R = 3.9521526379E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724642 A\n",
      "Wc = 22.6265508115 W/m\n",
      "Ws = 18.4603131519 W/m\n",
      "θj = 78.1711358296 °C\n",
      "θs = 80.4189021423 °C\n",
      "θc = 90.0000000016 °C\n",
      "Iteration: 100\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724456 A\n",
      "Wc = 22.6265508106 W/m\n",
      "Ws = 18.4603131508 W/m\n",
      "θj = 78.1711358272 °C\n",
      "θs = 80.4189021398 °C\n",
      "θc = 89.9999999987 °C\n",
      "Iteration: 101\n",
      "R = 3.9521526379E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  2 iterations for cable surface temperature\n",
      "I = 756.6452724605 A\n",
      "Wc = 22.6265508113 W/m\n",
      "Ws = 18.4603131517 W/m\n",
      "θj = 78.1711358291 °C\n",
      "θs = 80.4189021418 °C\n",
      "θc = 90.0000000010 °C\n",
      "Iteration: 102\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724486 A\n",
      "Wc = 22.6265508107 W/m\n",
      "Ws = 18.4603131510 W/m\n",
      "θj = 78.1711358276 °C\n",
      "θs = 80.4189021402 °C\n",
      "θc = 89.9999999992 °C\n",
      "Iteration: 103\n",
      "R = 3.9521526379E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724581 A\n",
      "Wc = 22.6265508112 W/m\n",
      "Ws = 18.4603131515 W/m\n",
      "θj = 78.1711358288 °C\n",
      "θs = 80.4189021415 °C\n",
      "θc = 90.0000000007 °C\n",
      "Iteration: 104\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724505 A\n",
      "Wc = 22.6265508108 W/m\n",
      "Ws = 18.4603131511 W/m\n",
      "θj = 78.1711358278 °C\n",
      "θs = 80.4189021405 °C\n",
      "θc = 89.9999999995 °C\n",
      "Iteration: 105\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724565 A\n",
      "Wc = 22.6265508111 W/m\n",
      "Ws = 18.4603131515 W/m\n",
      "θj = 78.1711358286 °C\n",
      "θs = 80.4189021413 °C\n",
      "θc = 90.0000000004 °C\n",
      "Iteration: 106\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724517 A\n",
      "Wc = 22.6265508109 W/m\n",
      "Ws = 18.4603131512 W/m\n",
      "θj = 78.1711358280 °C\n",
      "θs = 80.4189021406 °C\n",
      "θc = 89.9999999997 °C\n",
      "Iteration: 107\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724556 A\n",
      "Wc = 22.6265508111 W/m\n",
      "Ws = 18.4603131514 W/m\n",
      "θj = 78.1711358285 °C\n",
      "θs = 80.4189021411 °C\n",
      "θc = 90.0000000003 °C\n",
      "Iteration: 108\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724525 A\n",
      "Wc = 22.6265508109 W/m\n",
      "Ws = 18.4603131512 W/m\n",
      "θj = 78.1711358281 °C\n",
      "θs = 80.4189021407 °C\n",
      "θc = 89.9999999998 °C\n",
      "Iteration: 109\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724550 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131514 W/m\n",
      "θj = 78.1711358284 °C\n",
      "θs = 80.4189021411 °C\n",
      "θc = 90.0000000002 °C\n",
      "Iteration: 110\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724530 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358281 °C\n",
      "θs = 80.4189021408 °C\n",
      "θc = 89.9999999999 °C\n",
      "Iteration: 111\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724546 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358283 °C\n",
      "θs = 80.4189021410 °C\n",
      "θc = 90.0000000001 °C\n",
      "Iteration: 112\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724533 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021408 °C\n",
      "θc = 89.9999999999 °C\n",
      "Iteration: 113\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724543 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358283 °C\n",
      "θs = 80.4189021410 °C\n",
      "θc = 90.0000000001 °C\n",
      "Iteration: 114\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724535 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 89.9999999999 °C\n",
      "Iteration: 115\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724541 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358283 °C\n",
      "θs = 80.4189021410 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 116\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724536 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 117\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724540 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358283 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 118\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724537 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 119\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724540 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358283 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 120\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724538 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 121\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358283 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 122\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724538 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 123\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 124\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724538 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 125\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 126\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724538 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 127\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 128\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724538 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 129\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 130\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 131\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 132\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 133\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 134\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Iteration: 135\n",
      "R = 3.9521526380E-05 Ohm/m\n",
      "Rs = 2.0755417397E-04 Ohm/m\n",
      "λ1 = 0.8158695201\n",
      "Convergence reached after  1 iterations for cable surface temperature\n",
      "I = 756.6452724539 A\n",
      "Wc = 22.6265508110 W/m\n",
      "Ws = 18.4603131513 W/m\n",
      "θj = 78.1711358282 °C\n",
      "θs = 80.4189021409 °C\n",
      "θc = 90.0000000000 °C\n",
      "Convergence reached after 135 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 = 200 # 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",
    "W_a = 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 = np.zeros((maxIte,NC)) #[Ohm/m] AC resistance of conductor\n",
    "R_s = np.zeros((maxIte,NC)) #[Ohm/m] Resistance of sheath\n",
    "lambda_1 = np.zeros((maxIte,NC)) # Loss factor\n",
    "I = np.zeros((maxIte,NC)) #[A] Current rating\n",
    "W_c = np.zeros((maxIte,NC)) #[W/m] Conductor losses \n",
    "W_s = np.zeros((maxIte,NC)) #[W/m] Sheath losses\n",
    "theta_c = np.zeros((maxIte,NC)) #[°C] Conductor operation temperature\n",
    "theta_s = np.zeros((maxIte,NC)) #[°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",
    "Delta_theta_tr = theta_air # Morello [17]\n",
    "frrt_Delta_theta_s0 = np.ones(NC)*2 # From IEC 60287-2-1 Section 4.2.1.1\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[0])\n",
    "    \n",
    "    # 4.5.1.2 Skin effect factor\n",
    "    x_s = np.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[0], y_s[0]))\n",
    "\n",
    "    # 4.5.1.3 Proximity effect factor \n",
    "    x_p = np.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[0], y_p[0]))\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][0])\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][0])\n",
    "    \n",
    "    # 4.10.3.2 Losses caused by circulating currents (4.8.2.3)\n",
    "    # For three single-core cables in flat formation, without transposition, with sheaths bonded at both ends,\n",
    "    # the resistance of the sheath or screen at its operating temperaturehas a different formula for each phase,\n",
    "    # based on the intermediate variables P and Q:\n",
    "    P = X + X_m #[Ohm/m] From IEC 60287-1-1 Section 2.3.3\n",
    "    Q = X - X_m/3 #[Ohm/m]\n",
    "    \n",
    "    lambda_1m_prime = (R_s[i][1]/R[i][1])*(Q**2/(R_s[i][1]**2+Q**2)) # From IEC 60287-1-1 Section 2.3.3\n",
    "    lambda_11_prime = (R_s[i][0]/R[i][0])*(0.75*P**2/(R_s[i][0]**2+P**2)+0.25*Q**2/(R_s[i][0]**2+Q**2)+2*R_s[i][0]*P*Q*X_m/(mt.sqrt(3)*(R_s[i][0]**2+P**2)*(R_s[i][0]**2+Q**2)))\n",
    "    lambda_12_prime = (R_s[i][2]/R[i][2])*(0.75*P**2/(R_s[i][2]**2+P**2)+0.25*Q**2/(R_s[i][2]**2+Q**2)-2*R_s[i][2]*P*Q*X_m/(mt.sqrt(3)*(R_s[i][2]**2+P**2)*(R_s[i][2]**2+Q**2)))\n",
    "    \n",
    "    lambda_1_prime = np.array([lambda_11_prime, lambda_1m_prime, lambda_12_prime]) # Array of loss factor\n",
    "    \n",
    "    # 4.10.7 Variant with non-neglected eddy-current losses \n",
    "    # 4.8.5.1 Losses caused by eddy currents (Dielectric losses taken into account) (4.6.3.1)\n",
    "    rho_s = rho_s_20*(1+alpha_s_20*(theta_s[i]-20)) #[Ohm/m]\n",
    "    \n",
    "    beta_1 = np.sqrt(4*np.pi*omega/(1e7*rho_s)) # From IEC 60287-1-1 Section 2.3.6.1\n",
    "    m = (omega/R_s[i])*1e-7\n",
    "    D_s = dia[4] # Diameter after metal sheath\n",
    "    g_s = 1 + (t_s/D_s)**1.74*(beta_1*D_s*1e-3-1.6) # From IEC 60287-1-1 Section 2.3.6.1\n",
    "    # For cables in flat formation\n",
    "    # Centre cable:\n",
    "    lambda_0_m = 6*(m[1]**2/(1+m[1]**2))*(d/(2*s))**2 # From IEC 60287-1-1 Section 2.3.6.1 (2.a)\n",
    "    Delta_1_m = (0.86*m[1]**3.08)*(d/(2*s))**(1.4*m[1]+0.7)\n",
    "    Delta_2_m = 0\n",
    "    \n",
    "    # Outer cable leading phase:\n",
    "    lambda_0_2 = 1.5*(m[2]**2/(1+m[2]**2))*(d/(2*s))**2 # From IEC 60287-1-1 Section 2.3.6.1 (2.b)\n",
    "    Delta_1_2 = (4.7*m[2]**0.7)*(d/(2*s))**(0.16*m[2]+2)\n",
    "    Delta_2_2 = (21*m[2]**3.3)*(d/(2*s))**(1.47*m[2]+5.06)\n",
    "    \n",
    "    # Outer cable lagging phase:\n",
    "    lambda_0_1 = 1.5*(m[0]**2/(1+m[0]**2))*(d/(2*s))**2 # From IEC 60287-1-1 Section 2.3.6.1 (2.c)\n",
    "    Delta_1_1 = -((0.74*(m[0]+2)*m[0]**0.5)/(2+(m[0]-0.3)**2))*(d/(2*s))**(m[0]+1)\n",
    "    Delta_2_1 = (0.92*m[0]**3.7)*(d/(2*s))**(m[0]+2)\n",
    "    \n",
    "    # Arrays\n",
    "    lambda_0 = np.array([lambda_0_1, lambda_0_m, lambda_0_2])\n",
    "    Delta_1 = np.array([Delta_1_1, Delta_1_m, Delta_1_2])\n",
    "    Delta_2 = np.array([Delta_2_1, Delta_2_m, Delta_2_2])\n",
    "    \n",
    "    lambda_1_prime2 = R_s[i]/R[i]*(g_s*lambda_0*(1+Delta_1+Delta_2)+((beta_1*t_s)**4/12e12))\n",
    "    \n",
    "    # 4.10.7.1 Reduction of the eddy currents losses by the circulating currents (4.8.6.1)\n",
    "    M = R_s[i]/(X+X_m) # From IEC 60287-1-1 Section 2.3.5\n",
    "    N = R_s[i]/(X-X_m/3) # For cables in flat formation with equidistant spacing. \n",
    "    F = (4*M**2*N**2+(M+N)**2)/(4*(M**2+1)*(N**2+1))\n",
    "    \n",
    "    lambda_1_prime2 = F*lambda_1_prime2\n",
    "    \n",
    "    lambda_1[i] = lambda_1_prime+lambda_1_prime2\n",
    "    print(\"\\u03BB1 = %.10f\" % lambda_1[i][0])\n",
    "    \n",
    "    # 4.10.3.4 Heating due to dielectric losses (4.9.3.2)\n",
    "    Delta_theta_d = W_d*(((1/(1+lambda_1[i])+lambda_2)-1/2)*T1-(n*lambda_2*T2/(1+lambda_1[i]+lambda_2))) #[K] From IEC 60287-2-1 Section 4.2.1.1\n",
    "    \n",
    "    # 4.10.3.5 KA factor (4.9.3.3)\n",
    "    K_A = (mt.pi*D_e*h/(1+lambda_1[i]+lambda_2))*(T1/n+T2*(1+lambda_1[i])+T3*(1+lambda_1[i]+lambda_2)) #[K^(-1/4)] From IEC 60287-2-1 Section 4.2.1.1\n",
    "           \n",
    "    # 4.10.3.6 Cable surface temperature\n",
    "    # Iterative method described in Section 4.9.6.4 -> Section 4.9.3.5.\n",
    "    #Delta_theta = theta_c[i] - (theta_a + Delta_theta_tr) # IEC 60287-2-1 rating is “calculated as for a cable in free air, but the ambient temperature shall be increased by Delta_theta_tr\n",
    "    Delta_theta = theta_c_max - (theta_a + theta_air) # Always keep theta_c = 90°C and Delta_theta_tr 25°C for convergence?\n",
    "    Delta_theta_ds = 0 # The cables are being protected by the trough cover and hence no direct solar radiation is considered\n",
    "    \n",
    "    Tol2 = 1e-12 # Tolerance threshold for convergence\n",
    "    maxIte2 = 20 # Maximum number of iterations\n",
    "                  \n",
    "    for j in range(maxIte2):\n",
    "        frrt_Delta_theta_s = ((Delta_theta+Delta_theta_d+Delta_theta_ds)/(1+K_A*frrt_Delta_theta_s0))**0.25 #[K^(1/4)] From IEC 60287-2-1 Section 4.2.1.2\n",
    "        error = abs(frrt_Delta_theta_s-frrt_Delta_theta_s0)\n",
    "            \n",
    "        if all(error < Tol2):\n",
    "            print(\"Convergence reached after %2d iterations for cable surface temperature\" % (j+1))\n",
    "            break\n",
    "        elif j+1 == maxIte2:\n",
    "            print(\"Convergence criteria not reached after %2d iterations for cable surface temperature\" % maxIte2)\n",
    "        else:\n",
    "            frrt_Delta_theta_s0 = frrt_Delta_theta_s # Update for next iteration\n",
    "        \n",
    "    # 4.10.3.7 External thermal resistance (4.9.3.6)\n",
    "    T4 = 1/(mt.pi*D_e*h*frrt_Delta_theta_s) #[Km/W] From IEC 60287-2-1 Section 4.2.1.1\n",
    "    \n",
    "    # 4.10.3.8 Permissible current rating (4.5.2.3)\n",
    "    #Delta_theta = theta_c[i] - theta_a\n",
    "    Delta_theta = theta_c_max - theta_a # Always keep 90°C for convergence?\n",
    "    I[i] = np.sqrt((Delta_theta-W_d*(0.5*T1+n*(T2+T3+T4))-Delta_theta_tr)/(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][0])\n",
    "    \n",
    "    # 4.10.3.9 Losses in conductor and sheath (4.5.2.4)\n",
    "    W_c[i] = R[i][0]*I[i][0]**2 #[W/m] Only consider the current in the cable 11, why?\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][0], W_s[i][0]))\n",
    "    \n",
    "    # 4.10.3.11 Heating of the air inside the trough\n",
    "    W_TOT_i = W_c[i]+W_d+W_s[i]+W_a\n",
    "    W_TOT = sum(W_TOT_i)\n",
    "    \n",
    "    Delta_theta_tr = W_TOT/(3*p)\n",
    "    \n",
    "    # 4.10.3.10 Temperature on cable components (4.5.2.5)\n",
    "    theta_j = theta_a + n*(W_c[i]+W_s[i]+W_d)*T4 + Delta_theta_tr #[°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[0], theta_s[i][0], theta_c[i][0]))\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 = np.maximum(Delta_theta_s, Delta_theta_c)\n",
    "    Delta_I = abs(I[i] - I[i-1])\n",
    "    \n",
    "    if all(Delta_theta1 < tempTol) and all(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": 15,
   "id": "5957dc03-af4f-4841-870a-9018df3a14ad",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAkMAAAGwCAYAAACq12GxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAABXtklEQVR4nO3deXxU5b0/8M+ZNQvJkD0EYsIiIhBEQWuBskgMIojLLahVlhupVQFFQCvXetXeIog/xQoutY1Bib1wW4EqVjBRQBBFCCBEUVAgYUkISxYSkplk5vz+mJwz58wSM3E2Zj7v12teJmcmk2eOveZzv8/3eR5BFEURRERERBFKE+wBEBEREQUTwxARERFFNIYhIiIiimgMQ0RERBTRGIaIiIgoojEMERERUURjGCIiIqKIpgv2AEKBzWbDqVOnEBcXB0EQgj0cIiIi6gBRFHHhwgVkZGRAo+l8fYdhCMCpU6eQmZkZ7GEQERFRJxw/fhw9evTo9M8zDAGIi4sDYL+Z8fHxQR4NERERdUR9fT0yMzPlv+OdxTAEyFNj8fHxDENERESXmJ/b4sIGaiIiIopoDENEREQU0RiGiIiIKKIxDBEREVFEYxgiIiKiiMYwRERERBGNYYiIiIgiGsMQERERRTSGISIiIopoDENEREQU0RiGiIiIKKIxDBEREVFE40GtflR70YIGcyviovQwReuDPRwiIiJyg5UhP3p+43cY8fxmvLPjWLCHQkRERB4wDPmRIAgAAKsoBnkkRERE5AnDkB9p28KQzcYwREREFKoYhvxIY89CYBYiIiIKXQxDfqTRcJqMiIgo1DEM+ZE8TcYwREREFLIYhvxIqgyxZ4iIiCh0MQz5kUauDAV5IEREROQRw5AfSQ3UVqYhIiKikMUw5EfatjQksmeIiIgoZDEM+RE3XSQiIgp9DEN+JK0ms9qCPBAiIiLyiGHIj7Rtd5fTZERERKGLYciP5GkyNlATERGFLIYhP5IaqJmFiIiIQhfDkB85ziZjGiIiIgpVDEN+pOFxHERERCGPYciPNOwZIiIiCnlBDUPZ2dkQBMHlMWvWLADA6dOnMWPGDGRkZCAmJgY33XQTDh8+rHoPs9mMOXPmIDk5GbGxsZg0aRJOnDgRjI/jwtEzxDBEREQUqoIahnbt2oXKykr5UVxcDACYPHkyRFHEbbfdhiNHjuBf//oX9u7di6ysLOTm5qKxsVF+j7lz52LdunVYvXo1tm/fjoaGBkycOBFWqzVYH0vmOKg1yAMhIiIij3TB/OUpKSmq75csWYLevXtj1KhROHz4ML788kuUlZVhwIABAIDXXnsNqamp+N///V/MnDkTdXV1KCgowKpVq5CbmwsAKCoqQmZmJkpKSjBu3Di3v9dsNsNsNsvf19fX++XzyWeTsTJEREQUskKmZ8hisaCoqAj5+fkQBEEOK1FRUfJrtFotDAYDtm/fDgAoLS1FS0sL8vLy5NdkZGRg4MCB2LFjh8fftXjxYphMJvmRmZnpl88k7UDNTReJiIhCV8iEofXr16O2thYzZswAAPTr1w9ZWVlYuHAhampqYLFYsGTJElRVVaGyshIAUFVVBYPBgISEBNV7paWloaqqyuPvWrhwIerq6uTH8ePH/fKZ2EBNREQU+kImDBUUFGD8+PHIyMgAAOj1erz33ns4dOgQEhMTERMTgy1btmD8+PHQarXtvpcoivLuz+4YjUbEx8erHv6g4aaLREREIS8kwlB5eTlKSkowc+ZM1fUhQ4Zg3759qK2tRWVlJTZu3Ihz586hZ8+eAID09HRYLBbU1NSofq66uhppaWkBG78n3HSRiIgo9IVEGCosLERqaiomTJjg9nmTyYSUlBQcPnwYu3fvxq233grAHpb0er28Cg0AKisrUVZWhmHDhgVk7O3h0noiIqLQF9TVZABgs9lQWFiI6dOnQ6dTD+cf//gHUlJScNlll+HAgQN45JFHcNttt8kN0yaTCffddx/mz5+PpKQkJCYmYsGCBcjJyZFXlwUTe4aIiIhCX9DDUElJCSoqKpCfn+/yXGVlJebNm4fTp0+jW7dumDZtGp566inVa5YtWwadTocpU6agqakJY8eOxcqVK3+yrygQ5OM4uM8QERFRyBJErvtGfX09TCYT6urqfNpMvbGsEg8U7cHQrAT888HgT9sRERGFE1/9/Q6JnqFwJa1o46aLREREoYthyI+0ApfWExERhTqGIT/StN1dG9MQERFRyGIY8iO5gZrTZERERCGLYciPuLSeiIgo9DEM+RE3XSQiIgp9DEN+pGEDNRERUchjGPIj+WwypiEiIqKQxTDkR5wmIyIiCn0MQ37ETReJiIhCH8OQH8mVIZ5NRkREFLIYhvxI7hliZYiIiChkMQz5ETddJCIiCn0MQ34kTZNZOU1GREQUshiG/IiVISIiotDHMORHWumgVoYhIiKikMUw5EcCzyYjIiIKeQxDfqRtC0MsDBEREYUuhiE/4qn1REREoY9hyI807BkiIiIKeQxDfsSzyYiIiEIfw5AfcZqMiIgo9DEM+ZFjn6EgD4SIiIg8YhjyI+lsMgCwMRERERGFJIYhP9Iq0hD7hoiIiEITw5AfSZsuAoCVYYiIiCgkMQz5kbIyxCxEREQUmhiG/EirrAyxZ4iIiCgkMQz5kaBsoGZpiIiIKCQxDPmRqoHaFsSBEBERkUcMQ36kaaeB+rUtP+DuN79Ec4s10MMiIiIiBYYhP9K0M032v19V4Isj5/DNqfoAj4qIiIiUGIb8SBAEuW/IedNFq9X+PRuriYiIgothyM+0Ho7kkKbNWtlMREREFFQMQ34mH9bqNE1mtUn/ZGWIiIgomBiG/EzTdoddpsnaKkKtDENERERBxTDkZ45pMucw1NYzZGUYIiIiCqaghqHs7Oy2JmP1Y9asWQCAhoYGzJ49Gz169EB0dDSuvPJKvP7666r3MJvNmDNnDpKTkxEbG4tJkybhxIkTwfg4bsnTZE4VIOlbVoaIiIiCK6hhaNeuXaisrJQfxcXFAIDJkycDAB599FFs3LgRRUVFOHjwIB599FHMmTMH//rXv+T3mDt3LtatW4fVq1dj+/btaGhowMSJE2G1hsb+PRqNhwbqtgvcmZqIiCi4ghqGUlJSkJ6eLj82bNiA3r17Y9SoUQCAL774AtOnT8fo0aORnZ2N+++/H1dddRV2794NAKirq0NBQQFefPFF5Obm4uqrr0ZRUREOHDiAkpISj7/XbDajvr5e9fAXaa8hl2kyeTUZwxAREVEwhUzPkMViQVFREfLz8yG0TS2NGDEC77//Pk6ePAlRFLF582YcOnQI48aNAwCUlpaipaUFeXl58vtkZGRg4MCB2LFjh8fftXjxYphMJvmRmZnpt8+l1bjvGZIaqq1cWk9ERBRUIROG1q9fj9raWsyYMUO+9sorr6B///7o0aMHDAYDbrrpJrz22msYMWIEAKCqqgoGgwEJCQmq90pLS0NVVZXH37Vw4ULU1dXJj+PHj/vlMwGQg51zz5BcGWIDNRERUVDpgj0ASUFBAcaPH4+MjAz52iuvvIIvv/wS77//PrKysvDZZ5/hoYceQrdu3ZCbm+vxvURRlEOIO0ajEUaj0afj90RaTaYsDImiKH/PfYaIiIiCKyTCUHl5OUpKSrB27Vr5WlNTE/7rv/4L69atw4QJEwAAgwYNwr59+/D//t//Q25uLtLT02GxWFBTU6OqDlVXV2PYsGEB/xzuSNNkytCj/Jo9Q0RERMEVEtNkhYWFSE1NlUMPALS0tKClpQUajXqIWq0WtrY+myFDhkCv18ur0ACgsrISZWVlIROGBDcN1K0eghEREREFXtArQzabDYWFhZg+fTp0Osdw4uPjMWrUKDz22GOIjo5GVlYWtm7dinfeeQcvvfQSAMBkMuG+++7D/PnzkZSUhMTERCxYsAA5OTntTqMFkrsGak/BiIiIiAIv6GGopKQEFRUVyM/Pd3lu9erVWLhwIe655x6cP38eWVlZWLRoER544AH5NcuWLYNOp8OUKVPQ1NSEsWPHYuXKldBqtYH8GB45Nl10XFNPmXE1GRERUTAFPQzl5eVB9LDxYHp6OgoLC9v9+aioKCxfvhzLly/3x/B+Nnf7DCnzDytDREREwRUSPUPhTKoMKQ9qVZ5gz7PJiIiIgothyM+0bo7j4GoyIiKi0MEw5Gdyz5CHBmquJiMiIgouhiE/k3YG8BSAWBkiIiIKLoYhP9O66xmyua8SERERUeAxDPmZu7PJVJUhNlATEREFFcOQn7ltoBa5zxAREVGoYBjyM/f7DLFniIiIKFQwDPmZvM+Q6GGfIYYhIiKioGIY8jPNT/UMMQwREREFFcOQn0k9Q8pFYzYP55QRERFR4DEM+ZlG46YyxFPriYiIQgbDkJ+5a6DmqfVEREShg2HIz7RuGqiVX3OfISIiouBiGPIzx6aLjmvKAMSeISIiouBiGPIzrZuzyWzsGSIiIgoZuo68KDEx0as3FQQBe/bsQVZWVqcGFU7c7jPkYZk9ERERBV6HwlBtbS1efvllmEymn3ytKIp46KGHYLVaf/bgwoG0mszmYTUZwxAREVFwdSgMAcBdd92F1NTUDr12zpw5nR5QuJEaqK2qfYYYhoiIiEJFh8KQzcvl3xcuXOjUYMKRtLRe9FANauXSeiIioqDyWQO11WrF+vXrffV2YcPdpoue+oeIiIgo8Do8TebJd999h7feegtvv/02ampqYLFYfDGusCGfTaYKQI7nuZqMiIgouDpVGWpsbMRbb72F4cOHY8CAAdizZw8WLVqEU6dO+Xp8lzypZ0h5NhkbqImIiEKHV5WhL774An/729/wf//3f7j88stxzz33YOfOnXjllVfQv39/f43xkqZpi5uejuBgZYiIiCi4OhyG+vfvj4sXL+I3v/kNdu7cKYefJ554wm+DCwfu9xmC4muGISIiomDq8DTZDz/8gJEjR2LMmDG48sor/TmmsKJ1s8+Q7SdWk+0+dh4/VHNFHhERUSB0OAwdPXoUV1xxBR588EH06NEDCxYswN69e+Wzt8g9R2XIcU3VM+R0UGtNowV3vvklZhTuCsj4iIiIIl2Hw1D37t3x5JNP4ocffsCqVatQVVWF4cOHo7W1FStXrsShQ4f8Oc5LlvvVZJ7PJjt/0QKrTUT1BXNgBkhERBThOrWa7IYbbkBRUREqKyuxYsUKfPrpp+jXrx8GDRrk6/Fd8qRNFz0d1OrcMySdaN9itak2aiQiIiL/+FmbLppMJjz00EPYvXs39uzZg9GjR/toWOHDXc9Qe5WhlrbualFkczUREVEg+GwH6sGDB+OVV17x1duFDamnytMKMpfKkOJ7i5VHdRAREflbh8LQNddcg5qamg6/6YgRI3Dy5MlODyqcaNvucMenyRwBqKWVlSEiIiJ/69A+Q/v27cPXX3+NxMTEDr3pvn37YDazARhw7ECtDECt7VSGWqysDBEREQVShzddHDt2bIcbernc3kFwE4ba22dI+X0LwxAREZHfdSgMHT161Os37tGjh9c/E460Gnc9Q46vbaI9HEmn27cqKkMMQ0RERP7XoTCUlZXl73GELWlpvbKqZnWqsFlFERrYX6gMQAxDRERE/uez1WSdkZ2dDUEQXB6zZs0CALfPCYKAF154QX4Ps9mMOXPmIDk5GbGxsZg0aRJOnDgRrI/kQiNXhtxPkzk/p/zawgZqIiIivwtqGNq1axcqKyvlR3FxMQBg8uTJAKB6rrKyEm+99RYEQcB//Md/yO8xd+5crFu3DqtXr8b27dvR0NCAiRMnwmq1BuUzOXO7A7XoeTl9C5fWExERBVSHG6j9ISUlRfX9kiVL0Lt3b4waNQoAkJ6ernr+X//6F8aMGYNevXoBAOrq6lBQUIBVq1YhNzcXAFBUVITMzEyUlJRg3LhxAfgU7ZNWkynzj0tlSNEn1MppMiIiooAKamVIyWKxoKioCPn5+W5Xo50+fRoffvgh7rvvPvlaaWkpWlpakJeXJ1/LyMjAwIEDsWPHDo+/y2w2o76+XvXwF+mjtL/RoiP0qBqoWxmGiIiI/M3rMNSrVy+cO3fO5Xptba1csemM9evXo7a2FjNmzHD7/Ntvv424uDjccccd8rWqqioYDAYkJCSoXpuWloaqqiqPv2vx4sUwmUzyIzMzs9Pj/inycRztNVCrpskcAYjTZERERP7ndRg6duyY234cs9n8s3adLigowPjx45GRkeH2+bfeegv33HMPoqKifvK9RFFsd6+jhQsXoq6uTn4cP3680+P+KW7DUDtHcKiX1rOBmoiIyN863DP0/vvvy19v2rQJJpNJ/t5qteKTTz5BdnZ2pwZRXl6OkpISrF271u3z27Ztw/fff481a9aorqenp8NisaCmpkZVHaqursawYcM8/j6j0Qij0dipsXpL3nTRw9lkzt9zaT0REVFgdTgM3XbbbQDsf9ynT5+uek6v1yM7OxsvvvhipwZRWFiI1NRUTJgwwe3zBQUFGDJkCK666irV9SFDhkCv16O4uBhTpkwBYF+BVlZWhqVLl3ZqLL6mdbOazNbOajLVQa3sGSIiIvK7DochW1tpo2fPnti1axeSk5N9MgCbzYbCwkJMnz4dOp3rcOrr6/GPf/zDbdAymUy47777MH/+fCQlJSExMRELFixATk6OvLos2KRNF23tNFBbVQ3U7BkiIiIKJK+X1nfmaI72lJSUoKKiAvn5+W6fX716NURRxN133+32+WXLlkGn02HKlCloamrC2LFjsXLlSmi1Wp+Os7M0bnuG1K9Rft/C4ziIiIgCqlP7DH3yySf45JNPUF1dLVeMJG+99ZZX75WXl9fuAbD3338/7r//fo/PR0VFYfny5Vi+fLlXvzdQHJsuOq65TpM57qGqf4jTZERERH7ndRh69tln8cc//hFDhw5Ft27deEL9T9C2rdcT21lN5mlpPVeTERER+Z/XYeiNN97AypUrMXXqVH+MJ+zIlSGb532GPC2tZ88QERGR/3m9z5DFYml32TqpSWFIOTXW3kGtqgZqTpMRERH5nddhaObMmfj73//uj7GEJXnTRUWuaXXedFHZNO1hzyEiIiLyD6+nyZqbm/Hmm2+ipKQEgwYNgl6vVz3/0ksv+Wxw4UBeWt+JyhDDEBERkf95HYb279+PwYMHAwDKyspUz7GZ2pXGzaaLrj1DHg5qZQM1ERGR33kdhjZv3uyPcYQtuWeo3U0X3U+TOTdQN7dY8c/SExjTLxXdu0b7Y7hEREQRx+ueIckPP/yATZs2oampCQDa3SsokjkOanVca/c4DuU0mVMD9cayKvxhfRle/Ph7P4yUiIgoMnkdhs6dO4exY8eib9++uPnmm1FZWQnA3lg9f/58nw/wUiftQG3taGWonR2ozzdaAADnGiw+HycREVGk8joMPfroo9Dr9aioqEBMTIx8/c4778TGjRt9Orhw4L6BWv0a9UGtns8mM7dVipparD4eJRERUeTyumfo448/xqZNm9CjRw/V9csvvxzl5eU+G1i40LrZZ8i5gdrq4TgOS6v6deZWewhqZhgiIiLyGa8rQ42NjaqKkOTs2bMwGo0+GVQ4EQTXnqF29xlqZ2m9tAkjwxAREZHveB2GRo4ciXfeeUf+XhAE2Gw2vPDCCxgzZoxPBxcOHJsuuu4zpHVzor2nYARwmoyIiMgfvJ4me+GFFzB69Gjs3r0bFosFjz/+OL755hucP38en3/+uT/GeEmTeoasbg5qNWg1aLJZVZWi9naglqbJmizcjJGIiMhXvK4M9e/fH/v378d1112HG2+8EY2Njbjjjjuwd+9e9O7d2x9jvKRp3FR/pK8NOvvt93g2mdOmi+YWW9s/WRkiIiLyFa8qQy0tLcjLy8Nf/vIXPPvss/4aU1iRG6gVxRy5MtQWhpRTY6pT6532GZJWl3GajIiIyHe8qgzp9XqUlZXx2A0vuDu1XpoyM2hdK0MtNs8N1FJlqNUm8twyIiIiH/F6mmzatGkoKCjwx1jCkqbtDlvdNFDLlSFbRxuoHRUhrigjIiLyDa8bqC0WC/72t7+huLgYQ4cORWxsrOp5nlqvpnGztF6qDBnlniHlQa2ej+NQbsLY1GJFXJTe5+MlIiKKNF6HobKyMlxzzTUAgEOHDqme4/SZK3fL56Xs464ypD6o1X0DNQA0c0UZERGRT3gVhqxWK5555hnk5OQgMTHRX2MKK+6O45CO3HDXM2Rtd2m9Igy1cpqMiIjIF7zqGdJqtRg3bhzq6ur8NZ6wI02TqQOP/Z9uK0Md2IEaAJosDENERES+4HUDdU5ODo4cOeKPsYQluWfI1tF9hjwvrVc2UHN5PRERkW94HYYWLVqEBQsWYMOGDaisrER9fb3qQWqOniHHNeUO1IDTPkPKZmqbqApRqmkyhiEiIiKf8LqB+qabbgIATJo0SdUwLYoiBEGA1co/0krSDtRW0fPSemk1mSiKaHFqmm6x2WDUaAG0H4a+PVWPdXtPYPYNl8MUzVVmREREHeV1GNq8ebM/xhG2pAZq0d2mi1IYavve6nSaPQC0WEUY2/4tWVRhSD2F9tqWH7BhfyV6pXTB3ddd5rPxExERhTuvw9CoUaP8MY6wpXXbQO28z5D9+1Z3YajVBhjtX7fXM1R7sQUAcL7R4qORExERRQavw9Bnn33W7vMjR47s9GDCkeBm00WbvOmiffpL6hlyd8SGdM1mU0+hOa8ma7S0AgDqm1t8NHIiIqLI4HUYGj16tMs1Ze8Qe4bUpAZqwB5oNBrB5aBWuTKkCDt6rYAWqyj3CVmcgpLzPkMXzfbvLzS3+vgTEBERhTevV5PV1NSoHtXV1di4cSOuvfZafPzxx/4Y4yVNkYVgFUWIoihXieTVZG0XlIe0RuntVSOpMmR26hFq9lAZYhgiIiLyjteVIZPJ5HLtxhtvhNFoxKOPPorS0lKfDCxcaJSVIVGEssCj17qvDOm1Aow6DS4A8tSY2akS5NwzdNEiVYY4TUZEROQNr8OQJykpKfj+++999XZhQysop8kAaBxTYY4dqO0JSQpFOo1GDkpyZchpA0bn1WSNZlaGiIiIOsPrMLR//37V96IoorKyEkuWLMFVV13ls4GFC42grgzB5vjeuWdICj46rSCHIYuHMKSsDLVabfLz9U2sDBEREXnD6zA0ePBgCIKg2jcHAK6//nq89dZbPhtYuNAourKsogjlXXM+m0z6p16rgV5rD00trVIY8jxNdlHxNStDRERE3vE6DB09elT1vUajQUpKCqKionw2qHCirAyJNsAqOOKQUeuhMqQRFNNkUs+Q0zllyjBkVoYhVoaIiIi84XUYysrK8sc4wpayZ8gqilBkIRj16rPJpH/qNIJcNbK0bVXgfGirsjIkrSSzf22F1SaqlvQ3mFtx0dKK1DgGViIiImdeL61/+OGH8corr7hcX7FiBebOneuLMYUVRRayryYTlXsJOe9ALfUMaeRl95ZW95Uh5aaLysoQADQ4TZX9+vUdGP3CFtSxn4iIiMiF12Hovffew/Dhw12uDxs2DP/85z+9eq/s7GwIguDymDVrlvyagwcPYtKkSTCZTIiLi8P111+PiooK+Xmz2Yw5c+YgOTkZsbGxmDRpEk6cOOHtx/IbQRDkvYZsilPoNYK9AgQ4QpA0JaZsoHbsM6QOPMrVZA1mdfhR7kJttYn4/vQFXLRYcexso68+FhERUdjwOgydO3fO7V5D8fHxOHv2rFfvtWvXLlRWVsqP4uJiAMDkyZMBAD/++CNGjBiBfv36YcuWLfj666/x1FNPqfqT5s6di3Xr1mH16tXYvn07GhoaMHHixJDaCVurOLleapLWagTotOpzy+R9hjQa6HXqMCStKpPeS3lq/UWLOgwpm6jrmlogFaPOXDD77kMRERGFCa97hvr06YONGzdi9uzZqusfffQRevXq5dV7paSkqL5fsmQJevfuLR8G++STT+Lmm2/G0qVL5dcof0ddXR0KCgqwatUq5ObmAgCKioqQmZmJkpISjBs3zqvx+Iv9uBL7ztOOypAArUZ9an2LzbG03iCtJnPagdoUrcf5RotTz5A6+CkrQzUXHQe3nmlgGCIiInLmdWVo3rx5ePzxx/H0009j69at2Lp1K/77v/8bTzzxBB599NFOD8RisaCoqAj5+fkQBAE2mw0ffvgh+vbti3HjxiE1NRW/+MUvsH79evlnSktL0dLSgry8PPlaRkYGBg4ciB07dnj8XWazGfX19aqHP0lN1DabKB/SqtUIjmky5wZqrUaxz5C6Z8gUrQfgVBkye64MSafZA6wMERERueN1GMrPz8eLL76IgoICjBkzBmPGjEFRURFef/11/Pa3v+30QNavX4/a2lrMmDEDAFBdXY2GhgYsWbIEN910Ez7++GPcfvvtuOOOO7B161YAQFVVFQwGAxISElTvlZaWhqqqKo+/a/HixTCZTPIjMzOz0+PuCLlnSBTlKTGtIDimz9quWdsqQ3rF0nqL0z5D8W1hqL3KkHJ5fa2iMnSWlSEiIiIXnTqO48EHH8SDDz6IM2fOIDo6Gl26dPnZAykoKMD48eORkZEBALC1BYNbb71VrjgNHjwYO3bswBtvvCFPpbkjimLb1JR7CxcuxLx58+Tv6+vr/RqINIrQI1WGNBrXMNReA7UUirrKlSGb/DnbqwzVsDJERETULq8rQ0opKSk+CULl5eUoKSnBzJkz5WvJycnQ6XTo37+/6rVXXnmlvJosPT0dFosFNTU1qtdUV1cjLS3N4+8zGo2Ij49XPfxJ2njRJkI+qFWrCEOtTkvr9VqNvM+QYwfqtjAUo5ffV7rW0cqQcxhqMLdi4vJtWPzRwZ/x6YiIiC5tPysM+UphYSFSU1MxYcIE+ZrBYMC1117rcvjroUOH5I0fhwwZAr1eL69CA4DKykqUlZVh2LBhgRl8B0ihRzVNpugZcqkMadw0UEvTZFGOMCTtNdTearL2Gqj3VtSg7GQ9/v5lhcvxKkRERJHCZ6fWd5bNZkNhYSGmT58OnU49nMceewx33nknRo4ciTFjxmDjxo344IMPsGXLFgCAyWTCfffdh/nz5yMpKQmJiYlYsGABcnJy5NVloUCqDCmnyZQ9Q1JFSGqg1mpcG6ilabIYgxYGrQYWqw3NbQGpsW3TxViDFo0Wq9NqMsfXZ50qQ6dqmwAAF8ytONdoQXIXo68+MhER0SUj6JWhkpISVFRUID8/3+W522+/HW+88QaWLl2KnJwc/O1vf8N7772HESNGyK9ZtmwZbrvtNkyZMgXDhw9HTEwMPvjgA2i12kB+jHa5baDWCNBp3O9ArdcKLvsMSVNiRp0GUW3HeDhXhtJN9v2X6lWryRyVoUaLFY2K/qKTtc3y19yQkYiIItXPqgw1Nzf/7ANa8/Ly2p2iyc/PdxuUJFFRUVi+fDmWL1/+s8bhT/I0mc3RH6TRwKVnqMXd0vpW9T5DBp0GUXot6ptb5RVlUs9QuikKP55pVE+TNaqP4DjbYEas0f6vXaoMAcDRs40Ymp0ofy+KIr6trMcVaXHQaYOemYmIiPzG679yNpsN//M//4Pu3bujS5cuOHLkCADgqaeeQkFBgc8HGA4cDdTqaTK5Z0jeZ8ixtN65Z0jagdqo0yLaYK96SUdySNWetHh7ML3gYdNFQN1EXVmnDkNK7+05iQmvbMeKzT904hMTERFdOrwOQ3/605+wcuVKLF26FAaDQb4uTWORq7bZMFgV02Qat6vJHEvrHafWqxuojXoNovVSGJJ6htqmyeQw5LrporHt/ZR7DZ1STpOdU4eh7YfPAABKDp7uzEcmIiK6ZHgdht555x28+eabuOeee1R9OYMGDcJ3333n08GFC6kyJIqOg1q1guvZZFIVSDlNJk2dSdNkRp0GxrYw5OgZckyTAe4rQ71T7FsgSJUhURRxUjVNdlE15oOVF+R/NjrtY1TX1IIDJ+q8uwlEREQhyuswdPLkSfTp08flus1mQ0tLi5ufIK28msxxDpl6nyFb2/PSQa2KTRed9hky6DSIbmugllaTSQ3U0jRZfZP9+yaLVf65vmnqMHSu0SL3IwH2Bmqpd8vcasWPZxrkMX19olb1eeb/39e4ZcV2bP6+utP3hIiIKFR4HYYGDBiAbdu2uVz/xz/+gauvvtongwo3gpvVZBrBsZrMJtorNcoGaoOHHaiNOi2inCpD0tL6bm2VoaYWK1qsNrkqpNMIyE6OBeDYa6iybYosMdYArUZAU4sVp+vtz/1Q3SBP2QFA6THHppYN5lZsPWQPQWu+Ou7yWfcdr0Vp+Xmv7xEREVGweL2a7Omnn8bUqVNx8uRJ2Gw2rF27Ft9//z3eeecdbNiwwR9jvOQ5VpM5Gqh1WkdlCLBXYFrlaTIBep39OZeeIZ26Z8hqE+VVZVJlCAAamlvlMNQ1xoCUOPseQmcu2K9JU2SXJcag9qIFx85dxNGzjUg3RclTZJLSCkcY+vyHs3Jo+/S7atRdbIGpbVfsyromTHnjC1hFEf9++Fe4Ij1O9T5Vdc1I7mLg6jQiIgopXv9VuuWWW7BmzRr8+9//hiAI+O///m8cPHgQH3zwAW688UZ/jPGSJ2+6KIrycRwaQR2GWm2io4Ha7UGtymkyx2oy5YGtpmi9/NyF5la5eTohRo+Utg0VpcqQtKw+o2uUXDWSVpQdrKwHAFybbT8Ad095jdzrtOX7M/Lvs1ht+HdZpfz92zvKYbHaYLWJWPRv9REfH39ThRHPf4pbX/1c1dMEAOcbLfi/XcdxvlG98o2IiCgQOrXP0Lhx4zBu3DhfjyVsqc8mcz2OA7BflxuoNRqPB7UadVpHA3WLVT6kVSPYq0ZxUTo0tdh3oZYqQwmKypC0C7UchkzRbSHsjLyi7Lsqexi6/eoeKDtZj/rmVvxwpgGXp3bBlrY+oV/2SsIXR85h3d6TuPu6y9BobsXfd5bLn+ezQ2ew5ftqjL4iFafrm/H79/aj1Sbim1P1eLBoD96acS0MOg3KzzVi2ltfofzcRaR+bMSyOwdjeJ9keRyLPjyIw6cbcPd1l2HaL7OQEGuAKIrYf6IOG/afgkYQMKpvCoZmJ8Kg08DcakXZyTqUnaxHchcj+mfEIysxBhqNgNqLFhw524jqejNS443o3jUaKV2MaLWJONdoRnW9GRarDV2j9UiINdjvpcWK+qZW1De3QKsREGvQIcaohV6rgbnFiqYWe1+WTiPAqNfKq/YsrTa0WG1otYnQaQR7U7xGAJTnB4vqL0URENsu2r+2T5/6+qSUds4wbnv+J15ARNQJXaP18j53oSY0RxVmlNNkVsVqMpfKUNv0k16xtF5eTabYgTpaEYakDRdjDToIgoC4KB2qL5hxoblVPoqja4xeMU1mhiiKqKyz9wxldI2WV7UdbWuilqbJcrqbMDizK744cg67j9VAFIHKumYYdRo8d0cObnhxC746eh4nai6i+NvTqG9uRa/kWIy+IhVvfX4Uiz48iOF9krHgH1+j5mILeqfEorKuGdt/OIsn1u7Hfw7rif9c+RXONlggCED1BTPuLdiJ+3/VCxctVry7sxxS69KykkN4Y+uPGD8wHV+fqMWPZxxbAfzlsyOIM+rQK7ULDlbWqxrDAfsRJkadRnU0iUSnEVT9UURE5B/P3Z6D3/zismAPw60OhaGEhIQO/3+L58+zedaZ6jgOUbEDteBUGbI5ltY7N1BLPUMGnQbRhrbVZC2O4zVijPaAFNd2kOuF5hbUNjoqQ9K5YxarDfXNrXLPUEbXaMS0beJ47Gwjqi+Ycb7RAo0AXJ7WBUOzE/DFkXMoLa+Rp7eu75WEnsmx+EXPRHx55DzW7z2JNbvtzdT5I3rilqsysG7vCRyubsC9f9uJnUfPw6jT4C9Th+B4TRNmvr0ba/ecxL/2nYLVJqJ/t3i8ce8QvL71R/zvVxX4y2dH5PsyfmA6xlyRire/OIZvTtVj7d6TAIAovQZ5/dOh12qw9VA1zjZY8PXxWgBAUqwBV2V2xblGC76rrMdFi1XefqCbKQqpcUacuWBGVX2zamoyJc4Ig06D2ostqGtyBKdovRZxUTrYRPvKPem9NIL9OYNOg1abCEurzTGdqdVArxWg0Qht/WCOf7+AukAk/d+WAEfVRoCjiuRyvR0i2g92P1VlEp3GRkTkK6HcLtqhMPTyyy/7eRjhTaM4nd6mmCbTaARoBPv0WavNJleGVD1DUhhS7DMUpXM0UF9UVIYAID7aHobqlZWhWD2i2v6gX2huxZkLZnmarHvXaHRta4AuP3cR35yy7x/UK6ULovRaXJNl7xsqLT8v/8yYK1IAALdf3R1fHjmP5Z/+AHOrDQkxevzHNT0QbdBibm5fPP3+N9h51B6On5xwJfqkxqFPahwW356Dx9/bD6tNxLDeSfjL1CGIi9Jj8R05GNU3GU+uK0O6KQp/mNAfv+ydBACYPLQHPv/hHEoOnkb/jHiMH5guBz+bTcT+k3UoP9eInO4m9EyOlQNGq9WGY+caYWkVkZ0cgxiD43/yrVYbqi+YEa3XwhStl/89Sc81mq3ylJiS1Sai1WaDQatx+X8SpO0JONVERHTp6FAYmj59ur/HEdbc9QxJ13Qajdx07Dio1V5VABSVIek4Dr3jOI4mxcGr0jxsXJT9nxeaW+RDWhNi7DuFp8QZcaGtKlTd1jvUrWsUEmIMMGjt4/jkoL0n6Mpu8QCAay6zh6Fj5y7ieI09DI2+IhUAcNPAbnjqX9/I1ZCp12fJY/vNLy7DO18cw49nGjH6ihRMvT5Lvh9Trs2EQadBxfmL+N2oXjDqHJt33jSwG/L6p6uCCWAPFyMuT8aIy5Nd769GwODMrhic2dXlOZ1Wgz6pcS7XpecyukZ7fM4U4/7/jbHvEeX+IGCGICKiS0+neoasVivWrVuHgwcPQhAEXHnllbj11luh07EFyR2t4mwy5aaL8j+tsE+jyPsMqVeTiaKoaKDWyPsMNbfY0Ni24aI01RUvh6FWRQO1vYKS0sWII2caUXbSXv0x6DRIijVAEARclhSDH6obsOkb+/Eb/dqWxZui9eib1gWHTjfAahORnRQjrz4zReuRe2Uq/n2gCgadBlN/mS1/Zr1Wg9fuGYJ1e0/idyN7uYSE267u7vF+OQchIiIif/I6vZSVleHWW29FVVUVrrjiCgDAoUOHkJKSgvfffx85OTk+H+SlTj6bzKY+jgOA47BWxT5Deo1G1UBtVjQEGxRhyL6arG2aTK4MOXqGHA3UjsoQAOxv21G6e9doOaRkJ8Xih+oG+eyy/m2VIQAYkpWIQ6ftO1JLVSHJ9F9mY9M3p/Gfw7Ll95dckR6HJ8b36/B9IiIiCgav25lmzpyJAQMG4MSJE9izZw/27NmD48ePY9CgQbj//vv9McZLnqa9ypDWcVir6qBWxXEcUt8Q4G41mboyFGd0VIacp8mkJuqvj9srQxldHZs09kqJVY35SlUYSpC/Ht3WLyT5Ra8klD0zjqGHiIguWV5Xhr7++mvs3r0bCQmOP5AJCQlYtGgRrr32Wp8OLlzIS+tFdQM14FwZchzHoVecWi81TwP2VUrSajKzmwbqONU0mWNpPeCoDFXVty2rNzn6ZbKTHGGoa4weafGOKs8veiZCIwAxBh2u75Xk8vmkPiEiIqJLkddh6IorrsDp06cxYMAA1fXq6mq3B7iSo6nWZlM0UCt7htC2mkxaWq8RVA3UymX1giDIq8ma2llaX3PRgvpm92FI0k3RPJydHCN/fWV6vKrHJzMxBm/NuBbx0Xp5io6IiChcdCgM1dfXy18/99xzePjhh/HMM8/g+uuvBwB8+eWX+OMf/4jnn3/eP6O8xLXlGlhFx1SYVrGaDJB2oHYsrZemyWyifQk9AHl34yiDomfIw9L6k7VN8p4yXaPVPUOS7sppsuQu8tfKKTKJc68QERFRuOhQGOratauqUiCKIqZMmSJfk/ZWueWWW2C1Wt2+RyTTKO6Tzd1qMkg9Q8ql9Y52rgvN9uqPtARdeTaZa2XI/q/0ZNsy+C5GndyMLZ1PJlEuK0+LNyJar0VTixX9urlfik5ERBSOOhSGNm/e7O9xhDXHpotQHdQKqI/qaHWztB5QhqG2ypAUhiyee4akCpQ0RQa4VoaUYUgQBAzJsu82/YueiT/r8xIREV1KOhSGRo0a5e9xhDXlPkOOylDbc4rKkPqgVkclrsGsDkNyZajVdTVZfJQj/ACOlWQAkBhrgCA4jmRQNlADwN+mD8X5RovHjQiJiIjCUYfC0P79+zFw4EBoNBrs37+/3dcOGjTIJwMLJ9I+QzZRVJ1aDzitJrM5DmoVBEHeFVoKQwanMNRiFeUztJx3oJYoK0N6rQaJMQaca7QgIUbvsgosSq9lECIioojToTA0ePBgVFVVITU1FYMHD4YgCHKfkJIgCOwZckOaErPaXMOQqmdIsbQesIciixVokKbJ2kKQUe+YQjvXYN9LSKoMdTGq/5UqK0OAfarsHKs/REREsg6FoaNHjyIlJUX+mryjPJtMniZz2YHappgms1/T6zSAxeqYJmsLSUadRp7uOt92Mr1UGdJpNYgxaOVeooQY9bSZfePFCwxDREREbToUhrKystx+TR2jbJL2uM+Q1fGcXq4M2f8prRiTKkLSXkNNLY6gFKs4jT0+Si+Hoa5uKkOA/SgOIiIi6sRxHG+//TY+/PBD+fvHH38cXbt2xbBhw1BeXu7TwYULaVcC1XEcbvcZaqsMtTVPS3sN1TutJgNcd32ONTq+V/YNOVeGRvRJhkGnwYg+rqe/ExERRSKvw9Bzzz2H6Gh7VeGLL77AihUrsHTpUiQnJ+PRRx/1+QDDgRR8rKIIq7W9fYYcmy4CkFeUOTdQA44makmMojKkCkOx6srQfwzpgW+eHYfc/mk/81MRERGFB6+P4zh+/Lh87Mb69evx61//Gvfffz+GDx+O0aNH+3p8YUEKPKIIuTIkTZNJVSDns8kAR/iRp8l0jgCkbKIGnCtDjmqQ8zQZANUeRkRERJHO67+KXbp0wblz5wAAH3/8MXJzcwEAUVFRaGpq8u3owoSgWE1mczqOQ7XPkLQDtVwZsv/raXA3TaaoDAkC5PPKgPanyYiIiEjN68rQjTfeiJkzZ+Lqq6/GoUOHMGHCBADAN998g+zsbF+PLyxIhRirTXStDLX909JqkzdD1Dk1UF8wtx+GYvRa+f0AdWXIeWk9ERERqXldGXr11Vfxy1/+EmfOnMF7772HpKQkAEBpaSnuvvtunw8wHCjPJpOO43CuDEmHsQKuDdQNZvvGisqeIeXp8TFOewvFKypDXVkZIiIiapfXlaGuXbtixYoVLtefffZZnwwoHGkUDdTSNJkUeKTVZOZWm/x6fds1vc7+mkazdGq9IwApw1Cs08oy6eR6nUZw2YSRiIiI1LyuDG3cuBHbt2+Xv3/11VcxePBg/OY3v0FNTY1PBxculJsuytNkbdc07VSG5GmyZntlyNPSeuVKMsDRM9Q1xiD3KxEREZF7Xoehxx57DPX19QCAAwcOYP78+bj55ptx5MgRzJs3z+cDDAdSz5BN2UDddk3qGWpuVYQhpwbqlrZVZqppMsXXypVkgCMMsXmaiIjop3k9h3L06FH0798fAPDee+9h4sSJeO6557Bnzx7cfPPNPh9gONAoTq13rgxJPUPmFsdRHFI1Rxl+gI5XhvqkxAEA+qbF+ewzEBERhSuvw5DBYMDFixcBACUlJZg2bRoAIDExUa4YkZpGPn8M8saKzqfWm9sqQ9IUGeBooJYYFX1CytVkzpWhnB4mlMwbie5dY3z1EYiIiMKW12FoxIgRmDdvHoYPH46vvvoKa9asAQAcOnQIPXr08PkAw4FWURmyeTi1vlmuDDkCkF6r7vdRVoaUwci5MgQAfVJZFSIiIuoIr3uGVqxYAZ1Oh3/+8594/fXX0b17dwDARx99hJtuusmr98rOzoYgCC6PWbNmAQBmzJjh8tz111+veg+z2Yw5c+YgOTkZsbGxmDRpEk6cOOHtx/IraQsgm6g4qNXp1HqpgVpZGXLeKdrTcRxcMUZERNR5Xv8Vveyyy7BhwwaX68uWLfP6l+/atQtWq6NxuKysDDfeeCMmT54sX7vppptQWFgof28wqDcRnDt3Lj744AOsXr0aSUlJmD9/PiZOnIjS0lJoterpo2BxTJOJsInOlSF7wJHDkKoy5NwzpJwmczwXYwiNz0lERHQp6lAYqq+vR3x8vPx1e6TXdURKSorq+yVLlqB3794YNWqUfM1oNCI9Pd3tz9fV1aGgoACrVq2SjwUpKipCZmYmSkpKMG7cuA6PxZ9US+udjuOQKkHSNJlyaqy9BmrVPkOsDBEREXVah6bJEhISUF1dDcC+6WJCQoLLQ7reWRaLBUVFRcjPz1ftjbNlyxakpqaib9+++O1vfyuPA7Dvet3S0oK8vDz5WkZGBgYOHIgdO3Z4/F1msxn19fWqhz9JVSCbTUTbKnnXnqEONFAbPK4mY2WIiIioszpUUvj000+RmJgIANi8ebNfBrJ+/XrU1tZixowZ8rXx48dj8uTJyMrKwtGjR/HUU0/hhhtuQGlpKYxGI6qqqmAwGFxCWFpaGqqqqjz+rsWLFwd0x2xB0TPk3ECtc2qg1rc7TaZooNYpd6BmZYiIiKizOvRXVDltpfzalwoKCjB+/HhkZGTI1+68807564EDB2Lo0KHIysrChx9+iDvuuMPje4mi2O7OywsXLlRtEFlfX4/MzMyf+Qk80yqO45AbqJ0qQ2Z3DdQ659Vkip4hZWXIyMoQERFRZ3WqpNDc3Iz9+/ejuroaNptN9dykSZO8fr/y8nKUlJRg7dq17b6uW7duyMrKwuHDhwEA6enpsFgsqKmpUVWHqqurMWzYMI/vYzQaYTQavR5nZ0mBR1Qcx6H1tJpMURlqd5pMz8oQERGRL3j9V3Tjxo2YNm0azp496/KcIAiq1WEdVVhYiNTUVEyYMKHd1507dw7Hjx9Ht27dAABDhgyBXq9HcXExpkyZAgCorKxEWVkZli5d6vU4/EWqUlltjsqQlHPk1WStrg3U7U2TRXE1GRERkU94vc/Q7NmzMXnyZFRWVsJms6kenQlCNpsNhYWFmD59OnQ6RzZraGjAggUL8MUXX+DYsWPYsmULbrnlFiQnJ+P2228HAJhMJtx3332YP38+PvnkE+zduxf33nsvcnJy5NVloUDKN1Y3+wxJeccxTdZOz5DeQ2WIq8mIiIg6zeu/otXV1Zg3bx7S0tJ8MoCSkhJUVFQgPz9fdV2r1eLAgQN45513UFtbi27dumHMmDFYs2YN4uIcuysvW7YMOp0OU6ZMQVNTE8aOHYuVK1eGzB5DgKM/SBQ97zNkbqsMSdcBdztQOz5TlJ6ryYiIiHzB6zD061//Glu2bEHv3r19MoC8vDyIbQFBKTo6Gps2bfrJn4+KisLy5cuxfPlyn4zHHzRupsk0TqvJpDPLuM8QERFRYHn9V3TFihWYPHkytm3bhpycHOj1etXzDz/8sM8GFy7a23RRWQkCfqKBWvG98nBWhiEiIqLO8/qv6N///nds2rQJ0dHR2LJli2oJuyAIDENuSBnGpjiOQ+dUGZJ4aqA2aDVyNQmwH866IK8vAJ5NRkRE9HN4/Vf0D3/4A/74xz/iiSeegEbjdf91RNIoTq33tM+QRHU2mWJazHnKDABm33C5z8dKREQUabxOMxaLBXfeeSeDkBfkniHRPlUGKHagdmqSVp9a7/ja6CYMERER0c/n9V/Y6dOnY82aNf4YS9jSKlaTuS6tV/8rcJ4akzAMERER+YfX02RWqxVLly7Fpk2bMGjQIJcG6pdeeslngwsXUluVetNF9z1Dyu8NPzFNRkRERD+f12HowIEDuPrqqwEAZWVlqufaOw8skknBRxWGPK0m87DponKPISIiIvIdr8KQ1WrFM888g5ycHPkUe/pp0pSY8mwyaXasvcqQKgzpWRkiIiLyB6/+wmq1WowbNw51dXX+Gk9Y0ihOrbfZnHeg9txAzZ4hIiIi//P6L2xOTg6OHDnij7GELSnv2ERRrgzpPIQhZTVIr/O8GzURERH5htd/YRctWoQFCxZgw4YNqKysRH19vepBrqTAY7O5W03WwWky9gwRERH5hdcN1DfddBMAYNKkSaqGaVEUIQhCp06uD3fSBos2ES7TZDqnpfXKBmplNYjTZERERP7hdRjavHmzP8YR1lQHtYrtV4b0Gvc9Q5wmIyIi8g+vw9CoUaP8MY6wpnVzHIfHfYY8Lq1nGCIiIvKHTv2F3bZtG+69914MGzYMJ0+eBACsWrUK27dv9+ngwoWqgfonVpMpj+DQagT5Z9kzRERE5B9eh6H33nsP48aNQ3R0NPbs2QOz2QwAuHDhAp577jmfDzAcaBSbLkpnk0nTZC5nk3lYXcbKEBERkX94/Rf2T3/6E9544w389a9/VR3FMWzYMOzZs8engwsXUvBplZIQOjZNBjj6htgzRERE5B9e/4X9/vvvMXLkSJfr8fHxqK2t9cWYwo6Ub1qtijDk4aBWl8qQTqoMcZqMiIjIH7wOQ926dcMPP/zgcn379u3o1auXTwYVbqTKkMVqk69ptd5VhngcBxERkX94/Rf2d7/7HR555BHs3LkTgiDg1KlTePfdd7FgwQI89NBD/hjjJU+eJlOGIU9L67XOlSH79wYtwxAREZE/eL20/vHHH0ddXR3GjBmD5uZmjBw5EkajEQsWLMDs2bP9McZLnhR4WhTTZJ4PalWHHj0rQ0RERH7ldRgC7EdyPPnkk/j2229hs9nQv39/dOnSxddjCxvSRt2WDlSGnFeXydNk7BkiIiLyiw6XGy5evIhZs2ahe/fuSE1NxcyZM5GdnY3rrruOQegnOCpDNpdrPzlNxtVkREREftXhv7BPP/00Vq5ciQkTJuCuu+5CcXExHnzwQX+OLWxIPUNtJ3FAECCf6+Z6UKv6X0nXGPv2BYkxBj+PkoiIKDJ1eJps7dq1KCgowF133QUAuPfeezF8+HBYrVZotZzCaY9GUAcereJ714Na1a99+pb++PLIefyyd5L/BkhERBTBOhyGjh8/jl/96lfy99dddx10Oh1OnTqFzMxMvwwuXDhXfzSK712nydThqE9qHPqkxvlvcERERBGuw9NkVqsVBoN6qkan06G1tdXngwo3TnnHqTLU/nEcRERE5F8drgyJoogZM2bAaDTK15qbm/HAAw8gNjZWvrZ27VrfjjAMOE+TKQOPRiNAEBz9RM7TZkRERORfHQ5D06dPd7l27733+nQw4UrTzjQZYA9H0h5Ezj1DRERE5F8dDkOFhYX+HEdY0zo3UGtcv5fCkPPSeiIiIvIvzskEgHMbkOu0mcbt10REROR//MsbAM7TYs7HjCkrRZwmIyIiCiyGoQBwmSZrp6HaeWk9ERER+Rf/8gaA87SYa6XI81J7IiIi8i+GoQBwbgNyPYJDOU3GfyVERESBxL+8AdDecRyAulLE1WRERESBxTAUAO6W0iupKkNcTUZERBRQQf3Lm52dDUEQXB6zZs1yee3vfvc7CIKAl19+WXXdbDZjzpw5SE5ORmxsLCZNmoQTJ04E6BN0jFMhqN1wxJ4hIiKiwApqGNq1axcqKyvlR3FxMQBg8uTJqtetX78eO3fuREZGhst7zJ07F+vWrcPq1auxfft2NDQ0YOLEibBarQH5DB3hMi3mYZ8hjeDaXE1ERET+FdQwlJKSgvT0dPmxYcMG9O7dG6NGjZJfc/LkScyePRvvvvsu9Hq96ufr6upQUFCAF198Ebm5ubj66qtRVFSEAwcOoKSkJNAfxyOXniEPlSE2TxMREQVeyPz1tVgsKCoqQn5+PoS28GCz2TB16lQ89thjGDBggMvPlJaWoqWlBXl5efK1jIwMDBw4EDt27PD4u8xmM+rr61UPf/rJs8namqb1rAoREREFXMiEofXr16O2thYzZsyQrz3//PPQ6XR4+OGH3f5MVVUVDAYDEhISVNfT0tJQVVXl8XctXrwYJpNJfmRmZvrkM7RHWQ1yXjDGyhAREVHwhMxf34KCAowfP17uCyotLcWf//xnrFy5Uq4UdZQoiu3+zMKFC1FXVyc/jh8//rPG3hHKoo+n1WRcVk9ERBR4IRGGysvLUVJSgpkzZ8rXtm3bhurqalx22WXQ6XTQ6XQoLy/H/PnzkZ2dDQBIT0+HxWJBTU2N6v2qq6uRlpbm8fcZjUbEx8erHv6m7Bvy1EPEZfVERESBFxJ/fQsLC5GamooJEybI16ZOnYr9+/dj37598iMjIwOPPfYYNm3aBAAYMmQI9Hq9vAoNACorK1FWVoZhw4YF/HO0RxmAXCtD9n8NPKSViIgo8HTBHoDNZkNhYSGmT58Onc4xnKSkJCQlJaleq9frkZ6ejiuuuAIAYDKZcN9992H+/PlISkpCYmIiFixYgJycHOTm5gb0c/wUVc+Qh9VkPKSViIgo8IIehkpKSlBRUYH8/PxO/fyyZcug0+kwZcoUNDU1YezYsVi5ciW0Wq2PR/rzCB3oGeKGi0RERIEX9DCUl5cHURQ79Npjx465XIuKisLy5cuxfPlyH4/Mt9SrydxXhpxDEhEREfkf52UCRBmAnPcZ4jQZERFR8PCvb4Aol/p7qgyxgZqIiCjwGIYCRFn08bjPEJfWExERBRz/+gaIpt1pMi6tJyIiChaGoQBR7TPklHl0PI6DiIgoaPjXN0A0qmky9W3X8qBWIiKioGEYChCtagdq9XM6NlATEREFDcNQgGg6sAM1p8mIiIgCj399A6S9g1q5AzUREVHwMAwFiLadg1rl1WRcWk9ERBRw/OsbIMpikKfKkJ49Q0RERAHHMBQg7Z1a3zVGDwBIiDUEdExEREQUAge1RgpNO9Nkd16biVijDuMGpAd6WERERBGPYShAlKvJnKfJ4qL0uPu6ywI9JCIiIgKnyQJG2Q7EVWNEREShg2EoQNo7m4yIiIiCh2EoQFSbLgoMQ0RERKGCYShAlMUgbjRNREQUOvhnOUCUK8g4TUZERBQ6GIYCRLW0ntNkREREIYNhKEDa22eIiIiIgodhKEA07RzHQURERMHDMBQg7R3HQURERMHDMBQgnCYjIiIKTQxDAcIwREREFJoYhgJEy00XiYiIQhLDUIAo8w/3GSIiIgodDEMBom6gDuJAiIiISIV/lgNEdVArp8mIiIhCBsNQgLCBmoiIKDQxDAWIcmqMDdREREShg2EoQFTTZKwMERERhQyGoQBRBiAdwxAREVHIYBgKEA2X1hMREYUkhqEAUfYJsWeIiIgodDAMBYjA1WREREQhiWEoQJQBiPsMERERhY6ghqHs7GwIguDymDVrFgDgmWeeQb9+/RAbG4uEhATk5uZi586dqvcwm82YM2cOkpOTERsbi0mTJuHEiRPB+DjtUu9AzTBEREQUKoIahnbt2oXKykr5UVxcDACYPHkyAKBv375YsWIFDhw4gO3btyM7Oxt5eXk4c+aM/B5z587FunXrsHr1amzfvh0NDQ2YOHEirFZrUD6TJ8piEI/jICIiCh2CKIpisAchmTt3LjZs2IDDhw+remwk9fX1MJlMKCkpwdixY1FXV4eUlBSsWrUKd955JwDg1KlTyMzMxL///W+MGzfO7e8xm80wm82q983MzERdXR3i4+P98tmWbvwOr235EQCw7qFhuPqyBL/8HiIiokgh5YKf+/c7ZGoUFosFRUVFyM/PdxuELBYL3nzzTZhMJlx11VUAgNLSUrS0tCAvL09+XUZGBgYOHIgdO3Z4/F2LFy+GyWSSH5mZmb7/QE6UfUI6TcjcdiIioogXMn+V169fj9raWsyYMUN1fcOGDejSpQuioqKwbNkyFBcXIzk5GQBQVVUFg8GAhAR1lSUtLQ1VVVUef9fChQtRV1cnP44fP+7zz+NMubcQsxAREVHo0AV7AJKCggKMHz8eGRkZqutjxozBvn37cPbsWfz1r3/FlClTsHPnTqSmpnp8L1EU3VaXJEajEUaj0Wdj7wiNqmeIDdREREShIiRqFOXl5SgpKcHMmTNdnouNjUWfPn1w/fXXo6CgADqdDgUFBQCA9PR0WCwW1NTUqH6muroaaWlpARl7R3HTRSIiotAUEmGosLAQqampmDBhwk++VhRFufl5yJAh0Ov18io0AKisrERZWRmGDRvmt/F2hnqajGGIiIgoVAR9msxms6GwsBDTp0+HTucYTmNjIxYtWoRJkyahW7duOHfuHF577TWcOHFCXnpvMplw3333Yf78+UhKSkJiYiIWLFiAnJwc5ObmBusjuaVhZYiIiCgkBT0MlZSUoKKiAvn5+arrWq0W3333Hd5++22cPXsWSUlJuPbaa7Ft2zYMGDBAft2yZcug0+kwZcoUNDU1YezYsVi5ciW0Wm2gP0q7lHsLsWeIiIgodITUPkPB4qt9Ctrzt21H8KcPDwIAPn/iBnTvGu2X30NERBQpwm6foXDHaTIiIqLQxDAUIFxaT0REFJoYhgKEB7USERGFJoahABE4TUZERBSSGIYCRMvjOIiIiEIS/ywHiGoHak6TERERhQyGoQBRzoxpOE1GREQUMhiGAoQN1ERERKGJYShAuM8QERFRaGIYChAe1EpERBSaGIYCRMo/OgYhIiKikMIwFCDS1BirQkRERKGFYShApBDEfiEiIqLQwjAUIFIDNVeSERERhRaGoQDRtt1pZiEiIqLQwjAUIAIrQ0RERCGJYShAtAxDREREIYlhKECkniEexUFERBRaGIYCRDqpnpUhIiKi0MIwFCCcJiMiIgpNDEMBIu8zxDBEREQUUnTBHkCkuLJbPPqlxyH3yrRgD4WIiIgUGIYCpItRh41zRwZ7GEREROSE02REREQU0RiGiIiIKKIxDBEREVFEYxgiIiKiiMYwRERERBGNYYiIiIgiGsMQERERRTSGISIiIopoDENEREQU0RiGiIiIKKIxDBEREVFEYxgiIiKiiMYwRERERBGNYYiIiIgimi7YAwgFoigCAOrr64M8EiIiIuoo6e+29He8sxiGAFy4cAEAkJmZGeSREBERkbcuXLgAk8nU6Z8XxJ8bp8KAzWbDqVOnEBcXB0EQftZ71dfXIzMzE8ePH0d8fLyPRnhp4r1Q4/1w4L1w4L1Q4/1w4L1w8HQvRFHEhQsXkJGRAY2m850/rAwB0Gg06NGjh0/fMz4+PuL/xyvhvVDj/XDgvXDgvVDj/XDgvXBwdy9+TkVIwgZqIiIiimgMQ0RERBTRGIZ8zGg04umnn4bRaAz2UIKO90KN98OB98KB90KN98OB98LB3/eCDdREREQU0VgZIiIioojGMEREREQRjWGIiIiIIhrDEBEREUU0hiEfe+2119CzZ09ERUVhyJAh2LZtW7CH5HeLFy/Gtddei7i4OKSmpuK2227D999/r3qNKIp45plnkJGRgejoaIwePRrffPNNkEYcGIsXL4YgCJg7d658LdLuw8mTJ3HvvfciKSkJMTExGDx4MEpLS+XnI+V+tLa24g9/+AN69uyJ6Oho9OrVC3/84x9hs9nk14Tzvfjss89wyy23ICMjA4IgYP369arnO/LZzWYz5syZg+TkZMTGxmLSpEk4ceJEAD+Fb7R3L1paWvD73/8eOTk5iI2NRUZGBqZNm4ZTp06p3iMS7oWz3/3udxAEAS+//LLquq/uBcOQD61ZswZz587Fk08+ib179+JXv/oVxo8fj4qKimAPza+2bt2KWbNm4csvv0RxcTFaW1uRl5eHxsZG+TVLly7FSy+9hBUrVmDXrl1IT0/HjTfeKJ8LF2527dqFN998E4MGDVJdj6T7UFNTg+HDh0Ov1+Ojjz7Ct99+ixdffBFdu3aVXxMp9+P555/HG2+8gRUrVuDgwYNYunQpXnjhBSxfvlx+TTjfi8bGRlx11VVYsWKF2+c78tnnzp2LdevWYfXq1di+fTsaGhowceJEWK3WQH0Mn2jvXly8eBF79uzBU089hT179mDt2rU4dOgQJk2apHpdJNwLpfXr12Pnzp3IyMhwec5n90Ikn7nuuuvEBx54QHWtX79+4hNPPBGkEQVHdXW1CEDcunWrKIqiaLPZxPT0dHHJkiXya5qbm0WTySS+8cYbwRqm31y4cEG8/PLLxeLiYnHUqFHiI488Iopi5N2H3//+9+KIESM8Ph9J92PChAlifn6+6todd9wh3nvvvaIoRta9ACCuW7dO/r4jn722tlbU6/Xi6tWr5decPHlS1Gg04saNGwM2dl9zvhfufPXVVyIAsby8XBTFyLsXJ06cELt37y6WlZWJWVlZ4rJly+TnfHkvWBnyEYvFgtLSUuTl5amu5+XlYceOHUEaVXDU1dUBABITEwEAR48eRVVVlereGI1GjBo1KizvzaxZszBhwgTk5uaqrkfafXj//fcxdOhQTJ48Gampqbj66qvx17/+VX4+ku7HiBEj8Mknn+DQoUMAgK+//hrbt2/HzTffDCCy7oWzjnz20tJStLS0qF6TkZGBgQMHhv39qaurgyAIckU1ku6FzWbD1KlT8dhjj2HAgAEuz/vyXvCgVh85e/YsrFYr0tLSVNfT0tJQVVUVpFEFniiKmDdvHkaMGIGBAwcCgPz53d2b8vLygI/Rn1avXo09e/Zg165dLs9F0n0AgCNHjuD111/HvHnz8F//9V/46quv8PDDD8NoNGLatGkRdT9+//vfo66uDv369YNWq4XVasWiRYtw9913A4i8/20odeSzV1VVwWAwICEhweU14fzf1+bmZjzxxBP4zW9+Ix9OGkn34vnnn4dOp8PDDz/s9nlf3guGIR8TBEH1vSiKLtfC2ezZs7F//35s377d5blwvzfHjx/HI488go8//hhRUVEeXxfu90Fis9kwdOhQPPfccwCAq6++Gt988w1ef/11TJs2TX5dJNyPNWvWoKioCH//+98xYMAA7Nu3D3PnzkVGRgamT58uvy4S7oUnnfns4Xx/WlpacNddd8Fms+G11177ydeH270oLS3Fn//8Z+zZs8frz9WZe8FpMh9JTk6GVqt1SaPV1dUu/x9PuJozZw7ef/99bN68GT169JCvp6enA0DY35vS0lJUV1djyJAh0Ol00Ol02Lp1K1555RXodDr5s4b7fZB069YN/fv3V1278sor5QUFkfK/CwB47LHH8MQTT+Cuu+5CTk4Opk6dikcffRSLFy8GEFn3wllHPnt6ejosFgtqamo8viactLS0YMqUKTh69CiKi4vlqhAQOfdi27ZtqK6uxmWXXSb/97S8vBzz589HdnY2AN/eC4YhHzEYDBgyZAiKi4tV14uLizFs2LAgjSowRFHE7NmzsXbtWnz66afo2bOn6vmePXsiPT1ddW8sFgu2bt0aVvdm7NixOHDgAPbt2yc/hg4dinvuuQf79u1Dr169IuI+SIYPH+6yxcKhQ4eQlZUFIHL+dwHYVwlpNOr/3Gq1WnlpfSTdC2cd+exDhgyBXq9XvaayshJlZWVhd3+kIHT48GGUlJQgKSlJ9Xyk3IupU6di//79qv+eZmRk4LHHHsOmTZsA+PheeNVuTe1avXq1qNfrxYKCAvHbb78V586dK8bGxorHjh0L9tD86sEHHxRNJpO4ZcsWsbKyUn5cvHhRfs2SJUtEk8kkrl27Vjxw4IB49913i926dRPr6+uDOHL/U64mE8XIug9fffWVqNPpxEWLFomHDx8W3333XTEmJkYsKiqSXxMp92P69Oli9+7dxQ0bNohHjx4V165dKyYnJ4uPP/64/JpwvhcXLlwQ9+7dK+7du1cEIL700kvi3r175RVSHfnsDzzwgNijRw+xpKRE3LNnj3jDDTeIV111ldja2hqsj9Up7d2LlpYWcdKkSWKPHj3Effv2qf57ajab5feIhHvhjvNqMlH03b1gGPKxV199VczKyhINBoN4zTXXyMvLwxkAt4/CwkL5NTabTXz66afF9PR00Wg0iiNHjhQPHDgQvEEHiHMYirT78MEHH4gDBw4UjUaj2K9fP/HNN99UPR8p96O+vl585JFHxMsuu0yMiooSe/XqJT755JOqP3DhfC82b97s9r8R06dPF0WxY5+9qalJnD17tpiYmChGR0eLEydOFCsqKoLwaX6e9u7F0aNHPf73dPPmzfJ7RMK9cMddGPLVvRBEURS9qyURERERhQ/2DBEREVFEYxgiIiKiiMYwRERERBGNYYiIiIgiGsMQERERRTSGISIiIopoDENEREQU0RiGiIiIKKIxDBFRRMrOzsbLL78c7GEQUQhgGCIiv5sxYwZuu+02AMDo0aMxd+7cgP3ulStXomvXri7Xd+3ahfvvvz9g4yCi0KUL9gCIiDrDYrHAYDB0+udTUlJ8OBoiupSxMkREATNjxgxs3boVf/7znyEIAgRBwLFjxwAA3377LW6++WZ06dIFaWlpmDp1Ks6ePSv/7OjRozF79mzMmzcPycnJuPHGGwEAL730EnJychAbG4vMzEw89NBDaGhoAABs2bIF//mf/4m6ujr59z3zzDMAXKfJKioqcOutt6JLly6Ij4/HlClTcPr0afn5Z555BoMHD8aqVauQnZ0Nk8mEu+66CxcuXJBf889//hM5OTmIjo5GUlIScnNz0djY6Ke7SUS+wjBERAHz5z//Gb/85S/x29/+FpWVlaisrERmZiYqKysxatQoDB48GLt378bGjRtx+vRpTJkyRfXzb7/9NnQ6HT7//HP85S9/AQBoNBq88sorKCsrw9tvv41PP/0Ujz/+OABg2LBhePnllxEfHy//vgULFriMSxRF3HbbbTh//jy2bt2K4uJi/Pjjj7jzzjtVr/vxxx+xfv16bNiwARs2bMDWrVuxZMkSAEBlZSXuvvtu5Ofn4+DBg9iyZQvuuOMO8CxsotDHaTIiChiTyQSDwYCYmBikp6fL119//XVcc801eO655+Rrb731FjIzM3Ho0CH07dsXANCnTx8sXbpU9Z7K/qOePXvif/7nf/Dggw/itddeg8FggMlkgiAIqt/nrKSkBPv378fRo0eRmZkJAFi1ahUGDBiAXbt24dprrwUA2Gw2rFy5EnFxcQCAqVOn4pNPPsGiRYtQWVmJ1tZW3HHHHcjKygIA5OTk/Iy7RUSBwsoQEQVdaWkpNm/ejC5dusiPfv36AbBXYyRDhw51+dnNmzfjxhtvRPfu3REXF4dp06bh3LlzXk1PHTx4EJmZmXIQAoD+/fuja9euOHjwoHwtOztbDkIA0K1bN1RXVwMArrrqKowdOxY5OTmYPHky/vrXv6KmpqbjN4GIgoZhiIiCzmaz4ZZbbsG+fftUj8OHD2PkyJHy62JjY1U/V15ejptvvhkDBw7Ee++9h9LSUrz66qsAgJaWlg7/flEUIQjCT17X6/Wq5wVBgM1mAwBotVoUFxfjo48+Qv/+/bF8+XJcccUVOHr0aIfHQUTBwTBERAFlMBhgtVpV16655hp88803yM7ORp8+fVQP5wCktHv3brS2tuLFF1/E9ddfj759++LUqVM/+fuc9e/fHxUVFTh+/Lh87dtvv0VdXR2uvPLKDn82QRAwfPhwPPvss9i7dy8MBgPWrVvX4Z8nouBgGCKigMrOzsbOnTtx7NgxnD17FjabDbNmzcL58+dx991346uvvsKRI0fw8ccfIz8/v90g07t3b7S2tmL58uU4cuQIVq1ahTfeeMPl9zU0NOCTTz7B2bNncfHiRZf3yc3NxaBBg3DPPfdgz549+OqrrzBt2jSMGjXK7dScOzt37sRzzz2H3bt3o6KiAmvXrsWZM2e8ClNEFBwMQ0QUUAsWLIBWq0X//v2RkpKCiooKZGRk4PPPP4fVasW4ceMwcOBAPPLIIzCZTNBoPP9navDgwXjppZfw/PPPY+DAgXj33XexePFi1WuGDRuGBx54AHfeeSdSUlJcGrABe0Vn/fr1SEhIwMiRI5Gbm4tevXphzZo1Hf5c8fHx+Oyzz3DzzTejb9+++MMf/oAXX3wR48eP7/jNIaKgEESu+yQiIqIIxsoQERERRTSGISIiIopoDENEREQU0RiGiIiIKKIxDBEREVFEYxgiIiKiiMYwRERERBGNYYiIiIgiGsMQERERRTSGISIiIopoDENEREQU0f4/cxb9ASpBO48AAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "# Permissible current convergence plot\n",
    "Iplt = I[0:i+1,0]\n",
    "\n",
    "plt.figure()\n",
    "plt.plot(list(range(1,i+2)),Iplt)\n",
    "plt.xlabel('Iterations')\n",
    "plt.ylabel('Permissible current [A]')\n",
    "plt.show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "b5519496-3d0e-462c-aa1b-7c3ac8e4d9bf",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "I3 = 756.6453 A, I1 = 855.1481 A, I2 = 780.2569 A\n",
      "The permissible current rating of one conductor is equal to   750 A\n"
     ]
    }
   ],
   "source": [
    "# Permissible current rounding\n",
    "print(\"I3 = %.4f A, I1 = %.4f A, I2 = %.4f A\" % (I[i,0], I[i,1], I[i,2]))\n",
    "\n",
    "Imin = min(I[i,:])\n",
    "\n",
    "if Imin < 200:\n",
    "    I_rating = int(Imin)\n",
    "elif Imin < 500:\n",
    "    I_rating = int(Imin/5)*5\n",
    "else:\n",
    "    I_rating = int(Imin/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
}