docs/source/content/example_notebooks/wind_turbine.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Simulation NREL 5MW wind turbine"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/html": [
       "<img src=\"https://raw.githubusercontent.com/ImperialCollegeLondon/sharpy/dev_doc/docs/source/content/example_notebooks/images/turbulence_no_legend.png\" width=\"800\"/>"
      ],
      "text/plain": [
       "<IPython.core.display.Image object>"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "%config InlineBackend.figure_format = 'svg'\n",
    "from IPython.display import Image\n",
    "url = 'https://raw.githubusercontent.com/ImperialCollegeLondon/sharpy/dev_doc/docs/source/content/example_notebooks/images/turbulence_no_legend.png'\n",
    "Image(url=url, width=800)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "In this notebook:\n",
    "\n",
    "The blade loads on the NREL-5MW reference wind turbine computed with SHARPy and OpenFAST will be compared. However, zero-drag airfoils have been used."
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "OpenFAST: _https://openfast.readthedocs.io_"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "NREL-5MW: Jonkman, J.; Butterfield, S.; Musial, W. and Scott, G.. _Definition of a 5-MW Reference Wind Turbine for Offshore System Development_, Technical Report, NREL 2009"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Load the required packages:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Required packages\n",
    "%matplotlib inline\n",
    "import numpy as np\n",
    "import os\n",
    "import matplotlib.pyplot as plt\n",
    "\n",
    "# Required SHARPy modules\n",
    "import sharpy.sharpy_main\n",
    "import sharpy.utils.algebra as algebra\n",
    "import sharpy.utils.generate_cases as gc\n",
    "import sharpy.cases.templates.template_wt as template_wt"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "These are the results from the OpenFAST simulation for comparison: out-of-plane `of_cNdrR` and in-plane `of_cTdrR` coefficients along the blade and thrust `of_ct` and power `of_cp` rotor coefficients"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "of_rR = np.array([0.20158356, 0.3127131, 0.40794048, 0.5984148, 0.6936519, 0.85238045, 0.899999, 0.95555407, 0.98729974, 1.0])\n",
    "of_cNdrR = np.array([0.08621394, 0.14687876, 0.19345148, 0.2942731, 0.36003628, 0.43748564, 0.44762507, 0.38839236, 0.29782477, 0.0])\n",
    "of_cTdrR = np.array([0.048268348, 0.051957503, 0.05304592, 0.052862607, 0.056001827, 0.0536646, 0.050112925, 0.038993906, 0.023664437, 0.0])\n",
    "\n",
    "of_ct = 0.69787693\n",
    "of_cp = 0.48813498"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Create SHARPy case"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "We define our parameters:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Mathematical constants\n",
    "deg2rad = np.pi/180.\n",
    "\n",
    "# Case\n",
    "case = 'rotor'\n",
    "route = './'\n",
    "\n",
    "# Geometry discretization\n",
    "chord_panels = np.array([8], dtype=int)\n",
    "revs_in_wake = 5\n",
    "\n",
    "# Operation\n",
    "rotation_velocity = 12.1*2*np.pi/60\n",
    "pitch_deg = 0. #degrees\n",
    "\n",
    "# Wind\n",
    "WSP = 12.\n",
    "air_density = 1.225\n",
    "\n",
    "# Simulation\n",
    "dphi = 4.*deg2rad\n",
    "revs_to_simulate = 5"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Computation of associated parameters"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "dt = dphi/rotation_velocity\n",
    "time_steps = int(revs_to_simulate*2.*np.pi/dphi)\n",
    "mstar = int(revs_in_wake*2.*np.pi/dphi)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Generation of the rotor geometry based on information in the excel file:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[33mWARNING: The poisson cofficient is assumed equal to 0.3\u001b[0m\n",
      "\u001b[33mWARNING: Cross-section area is used as shear area\u001b[0m\n",
      "\u001b[33mWARNING: Using perpendicular axis theorem to compute the inertia around xB\u001b[0m\n",
      "\u001b[33mWARNING: Replacing node 29 by node 0\u001b[0m\n",
      "\u001b[33mWARNING: Replacing node 58 by node 0\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "op_params = {}                                                                                                                                                                                     \n",
    "op_params['rotation_velocity'] = rotation_velocity                                                                                                                                                 \n",
    "op_params['pitch_deg'] = pitch_deg                                                                                                                                                                 \n",
    "op_params['wsp'] = WSP                                                                                                                                                                             \n",
    "op_params['dt'] = dt                                                                                                                                                                               \n",
    "                                                                                                                                                                                                   \n",
    "geom_params = {}                                                                                                                                                                                   \n",
    "geom_params['chord_panels'] = chord_panels                                                                                                                                                         \n",
    "geom_params['tol_remove_points'] = 1e-8                                                                                                                                                            \n",
    "geom_params['n_points_camber'] = 100                                                                                                                                                               \n",
    "geom_params['h5_cross_sec_prop'] = None                                                                                                                                                            \n",
    "geom_params['m_distribution'] = 'uniform'                                                                                                                                                          \n",
    "                                                                                                                                                                                                   \n",
    "options = {}                                                                                                                                                                                       \n",
    "options['camber_effect_on_twist'] = False                                                                                                                                                          \n",
    "options['user_defined_m_distribution_type'] = None                                                                                                                                                 \n",
    "options['include_polars'] = False                                                                                                                                                                  \n",
    "options['separate_blades'] = False                                                                                                                                                                 \n",
    "                                                                                                                                                                                                   \n",
    "excel_description = {}                                                                                                                                                                             \n",
    "excel_description['excel_file_name'] = 'source/type04_db_nrel5mw_oc3_v06.xlsx'                                                                          \n",
    "excel_description['excel_sheet_parameters'] = 'parameters'                                                                                                                                         \n",
    "excel_description['excel_sheet_structural_blade'] = 'structural_blade'                                                                                                                             \n",
    "excel_description['excel_sheet_discretization_blade'] = 'discretization_blade'                                                                                                                     \n",
    "excel_description['excel_sheet_aero_blade'] = 'aero_blade'                                                                                                                                         \n",
    "excel_description['excel_sheet_airfoil_info'] = 'airfoil_info'                                                                                                                                     \n",
    "excel_description['excel_sheet_airfoil_chord'] = 'airfoil_coord'                                                                                                                                   \n",
    "                                                                                                                                                                                                   \n",
    "rotor, hub_nodes = template_wt.rotor_from_excel_type03(op_params,                                                                                                                                  \n",
    "                                            geom_params,                                                                                                                                           \n",
    "                                            excel_description,                                                                                                                                     \n",
    "                                            options) "
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Define simulation details. The steady simulation is faster than the dynamic simulation. However, the dynamic simulation includes wake self-induction and provides more accurate results."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [],
   "source": [
    "steady_simulation = False"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [],
   "source": [
    "SimInfo = gc.SimulationInformation()\n",
    "SimInfo.set_default_values()\n",
    "\n",
    "if steady_simulation:\n",
    "    SimInfo.solvers['SHARPy']['flow'] = ['BeamLoader',\n",
    "                            'AerogridLoader',\n",
    "                            'StaticCoupled',\n",
    "                            'BeamPlot',\n",
    "                            'AerogridPlot',            \n",
    "                            'SaveData'] \n",
    "else:\n",
    "    SimInfo.solvers['SHARPy']['flow'] = ['BeamLoader',\n",
    "                            'AerogridLoader',\n",
    "                            'StaticCoupled',\n",
    "                            'DynamicCoupled']\n",
    "                                         \n",
    "SimInfo.solvers['SHARPy']['case'] = case\n",
    "SimInfo.solvers['SHARPy']['route'] = route\n",
    "SimInfo.solvers['SHARPy']['write_log'] = True\n",
    "SimInfo.set_variable_all_dicts('dt', dt)\n",
    "SimInfo.set_variable_all_dicts('rho', air_density)\n",
    "\n",
    "SimInfo.solvers['SteadyVelocityField']['u_inf'] = WSP\n",
    "SimInfo.solvers['SteadyVelocityField']['u_inf_direction'] = np.array([0., 0., 1.])\n",
    "\n",
    "SimInfo.solvers['BeamLoader']['unsteady'] = 'on'\n",
    "\n",
    "SimInfo.solvers['AerogridLoader']['unsteady'] = 'on'\n",
    "SimInfo.solvers['AerogridLoader']['mstar'] = mstar\n",
    "SimInfo.solvers['AerogridLoader']['freestream_dir'] = np.array([0.,0.,0.])\n",
    "SimInfo.solvers['AerogridLoader']['wake_shape_generator'] = 'HelicoidalWake'\n",
    "SimInfo.solvers['AerogridLoader']['wake_shape_generator_input'] = {'u_inf': WSP,\n",
    "                                                                   'u_inf_direction': SimInfo.solvers['SteadyVelocityField']['u_inf_direction'],\n",
    "                                                                   'rotation_velocity': rotation_velocity*np.array([0., 0., 1.]),\n",
    "                                                                   'dt': dt,\n",
    "                                                                   'dphi1': dphi,\n",
    "                                                                   'ndphi1': mstar,\n",
    "                                                                   'r': 1.,\n",
    "                                                                   'dphimax': 10*deg2rad}\n",
    "    \n",
    "SimInfo.solvers['StaticCoupled']['structural_solver'] = 'RigidDynamicPrescribedStep'\n",
    "SimInfo.solvers['StaticCoupled']['structural_solver_settings'] = SimInfo.solvers['RigidDynamicPrescribedStep']\n",
    "SimInfo.solvers['StaticCoupled']['aero_solver'] = 'StaticUvlm'\n",
    "SimInfo.solvers['StaticCoupled']['aero_solver_settings'] = SimInfo.solvers['StaticUvlm']\n",
    "\n",
    "SimInfo.solvers['StaticCoupled']['tolerance'] = 1e-8\n",
    "SimInfo.solvers['StaticCoupled']['n_load_steps'] = 0\n",
    "SimInfo.solvers['StaticCoupled']['relaxation_factor'] = 0.\n",
    "\n",
    "SimInfo.solvers['StaticUvlm']['num_cores'] = 8\n",
    "SimInfo.solvers['StaticUvlm']['velocity_field_generator'] = 'SteadyVelocityField'\n",
    "SimInfo.solvers['StaticUvlm']['velocity_field_input'] = SimInfo.solvers['SteadyVelocityField']\n",
    "\n",
    "SimInfo.solvers['SaveData']['compress_float'] = True\n",
    "                                         \n",
    "# Only used for steady_simulation = False\n",
    "SimInfo.solvers['StepUvlm']['convection_scheme'] = 3\n",
    "SimInfo.solvers['StepUvlm']['num_cores'] = 8\n",
    "SimInfo.solvers['StepUvlm']['velocity_field_generator'] = 'SteadyVelocityField'\n",
    "SimInfo.solvers['StepUvlm']['velocity_field_input'] = SimInfo.solvers['SteadyVelocityField']\n",
    "\n",
    "SimInfo.solvers['DynamicCoupled']['structural_solver'] = 'RigidDynamicPrescribedStep'\n",
    "SimInfo.solvers['DynamicCoupled']['structural_solver_settings'] = SimInfo.solvers['RigidDynamicPrescribedStep']\n",
    "SimInfo.solvers['DynamicCoupled']['aero_solver'] = 'StepUvlm'\n",
    "SimInfo.solvers['DynamicCoupled']['aero_solver_settings'] = SimInfo.solvers['StepUvlm']\n",
    "SimInfo.solvers['DynamicCoupled']['postprocessors'] = ['BeamPlot', 'AerogridPlot', 'Cleanup', 'SaveData']\n",
    "SimInfo.solvers['DynamicCoupled']['postprocessors_settings'] = {'BeamPlot': SimInfo.solvers['BeamPlot'],\n",
    "                                                             'AerogridPlot': SimInfo.solvers['AerogridPlot'],\n",
    "                                                             'Cleanup': SimInfo.solvers['Cleanup'],\n",
    "                                                             'SaveData': SimInfo.solvers['SaveData']}\n",
    "SimInfo.solvers['DynamicCoupled']['minimum_steps'] = 0\n",
    "SimInfo.solvers['DynamicCoupled']['include_unsteady_force_contribution'] = True\n",
    "SimInfo.solvers['DynamicCoupled']['relaxation_factor'] = 0.\n",
    "SimInfo.solvers['DynamicCoupled']['final_relaxation_factor'] = 0.\n",
    "SimInfo.solvers['DynamicCoupled']['dynamic_relaxation'] = False\n",
    "SimInfo.solvers['DynamicCoupled']['relaxation_steps'] = 0\n",
    "\n",
    "# Define dynamic simulation (used regardless the value of \"steady_simulation\" variable)\n",
    "SimInfo.define_num_steps(time_steps)\n",
    "SimInfo.with_forced_vel = True\n",
    "SimInfo.for_vel = np.zeros((time_steps,6), dtype=float)\n",
    "SimInfo.for_vel[:,5] = rotation_velocity\n",
    "SimInfo.for_acc = np.zeros((time_steps,6), dtype=float)\n",
    "SimInfo.with_dynamic_forces = True\n",
    "SimInfo.dynamic_forces = np.zeros((time_steps,rotor.StructuralInformation.num_node,6), dtype=float)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Generate simulation files"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {},
   "outputs": [],
   "source": [
    "gc.clean_test_files(SimInfo.solvers['SHARPy']['route'], SimInfo.solvers['SHARPy']['case'])\n",
    "rotor.generate_h5_files(SimInfo.solvers['SHARPy']['route'], SimInfo.solvers['SHARPy']['case'])\n",
    "SimInfo.generate_solver_file()\n",
    "SimInfo.generate_dyn_file(time_steps)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Run SHARPy case"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "--------------------------------------------------------------------------------\u001b[0m\n",
      "            ######  ##     ##    ###    ########  ########  ##    ##\u001b[0m\n",
      "           ##    ## ##     ##   ## ##   ##     ## ##     ##  ##  ##\u001b[0m\n",
      "           ##       ##     ##  ##   ##  ##     ## ##     ##   ####\u001b[0m\n",
      "            ######  ######### ##     ## ########  ########     ##\u001b[0m\n",
      "                 ## ##     ## ######### ##   ##   ##           ##\u001b[0m\n",
      "           ##    ## ##     ## ##     ## ##    ##  ##           ##\u001b[0m\n",
      "            ######  ##     ## ##     ## ##     ## ##           ##\u001b[0m\n",
      "--------------------------------------------------------------------------------\u001b[0m\n",
      "Aeroelastics Lab, Aeronautics Department.\u001b[0m\n",
      "    Copyright (c), Imperial College London.\u001b[0m\n",
      "    All rights reserved.\u001b[0m\n",
      "    License available at https://github.com/imperialcollegelondon/sharpy\u001b[0m\n",
      "\u001b[36mRunning SHARPy from /home/arturo/code/sharpy/docs/source/content/example_notebooks\u001b[0m\n",
      "\u001b[36mSHARPy being run is in /home/arturo/code/sharpy\u001b[0m\n",
      "\u001b[36mThe branch being run is dev_blade_pitch_v2\u001b[0m\n",
      "\u001b[36mThe version and commit hash are: v1.2.1-546-g05602d1f-05602d1f\u001b[0m\n",
      "SHARPy output folder set\u001b[0m\n",
      "\u001b[34m\t./output//rotor/\u001b[0m\n",
      "\u001b[36mGenerating an instance of BeamLoader\u001b[0m\n",
      "\u001b[36mGenerating an instance of AerogridLoader\u001b[0m\n",
      "Variable shear_direction has no assigned value in the settings file.\u001b[0m\n",
      "\u001b[34m    will default to the value: [1. 0. 0.]\u001b[0m\n",
      "Variable shear_exp has no assigned value in the settings file.\u001b[0m\n",
      "\u001b[34m    will default to the value: 0.0\u001b[0m\n",
      "Variable h_ref has no assigned value in the settings file.\u001b[0m\n",
      "\u001b[34m    will default to the value: 1.0\u001b[0m\n",
      "Variable h_corr has no assigned value in the settings file.\u001b[0m\n",
      "\u001b[34m    will default to the value: 1.0\u001b[0m\n",
      "\u001b[34mThe aerodynamic grid contains 3 surfaces\u001b[0m\n",
      "\u001b[34m  Surface 0, M=8, N=26\u001b[0m\n",
      "     Wake 0, M=450, N=26\u001b[0m\n",
      "\u001b[34m  Surface 1, M=8, N=26\u001b[0m\n",
      "     Wake 1, M=450, N=26\u001b[0m\n",
      "\u001b[34m  Surface 2, M=8, N=26\u001b[0m\n",
      "     Wake 2, M=450, N=26\u001b[0m\n",
      "  In total: 624 bound panels\u001b[0m\n",
      "  In total: 35100 wake panels\u001b[0m\n",
      "  Total number of panels = 35724\u001b[0m\n",
      "\u001b[36mGenerating an instance of StaticCoupled\u001b[0m\n",
      "\u001b[36mGenerating an instance of RigidDynamicPrescribedStep\u001b[0m\n",
      "\u001b[36mGenerating an instance of StaticUvlm\u001b[0m\n",
      "\u001b[0m\n",
      "\u001b[0m\n",
      "\u001b[0m\n",
      "|=====|=====|============|==========|==========|==========|==========|==========|==========|\u001b[0m\n",
      "|iter |step | log10(res) |    Fx    |    Fy    |    Fz    |    Mx    |    My    |    Mz    |\u001b[0m\n",
      "|=====|=====|============|==========|==========|==========|==========|==========|==========|\u001b[0m\n",
      "|  0  |  0  |  0.00000   | -0.0000  |  0.0000  |21927.8250|  0.0000  | -0.0000  |5962997.5361|\u001b[0m\n",
      "|  1  |  0  |    -inf    | -0.0000  |  0.0000  |21927.8250|  0.0000  | -0.0000  |5962997.5361|\u001b[0m\n",
      "\u001b[36mGenerating an instance of DynamicCoupled\u001b[0m\n",
      "\u001b[36mGenerating an instance of RigidDynamicPrescribedStep\u001b[0m\n",
      "\u001b[36mGenerating an instance of StepUvlm\u001b[0m\n",
      "\u001b[36mGenerating an instance of BeamPlot\u001b[0m\n",
      "\u001b[36mGenerating an instance of AerogridPlot\u001b[0m\n",
      "\u001b[36mGenerating an instance of Cleanup\u001b[0m\n",
      "\u001b[36mGenerating an instance of SaveData\u001b[0m\n",
      "\u001b[0m\n",
      "\u001b[0m\n",
      "\u001b[0m\n",
      "|=======|========|======|==============|==============|==============|==============|==============|\u001b[0m\n",
      "|  ts   |   t    | iter | struc ratio  |  iter time   | residual vel |  FoR_vel(x)  |  FoR_vel(z)  |\u001b[0m\n",
      "|=======|========|======|==============|==============|==============|==============|==============|\u001b[0m\n",
      "|   1   | 0.0551 |  1   |   0.001471   |  11.072773   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   2   | 0.1102 |  1   |   0.001529   |  11.094774   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   3   | 0.1653 |  1   |   0.001487   |  11.017530   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   4   | 0.2204 |  1   |   0.001487   |  11.053551   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   5   | 0.2755 |  1   |   0.001495   |  10.972424   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   6   | 0.3306 |  1   |   0.001496   |  10.951425   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   7   | 0.3857 |  1   |   0.001492   |  11.018625   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   8   | 0.4408 |  1   |   0.001476   |  11.044770   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|   9   | 0.4959 |  1   |   0.001478   |  11.122489   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  10   | 0.5510 |  1   |   0.001500   |  10.997606   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  11   | 0.6061 |  1   |   0.001481   |  11.076450   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  12   | 0.6612 |  1   |   0.001513   |  10.976762   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  13   | 0.7163 |  1   |   0.001489   |  11.056770   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  14   | 0.7713 |  1   |   0.001473   |  11.147662   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  15   | 0.8264 |  1   |   0.001485   |  11.047932   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  16   | 0.8815 |  1   |   0.001500   |  10.984375   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  17   | 0.9366 |  1   |   0.001492   |  10.999334   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  18   | 0.9917 |  1   |   0.001492   |  11.036468   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  19   | 1.0468 |  1   |   0.001490   |  10.989860   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  20   | 1.1019 |  1   |   0.001487   |  10.981977   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  21   | 1.1570 |  1   |   0.001480   |  11.026050   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  22   | 1.2121 |  1   |   0.001501   |  10.964416   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  23   | 1.2672 |  1   |   0.001483   |  11.061231   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  24   | 1.3223 |  1   |   0.001502   |  11.026285   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  25   | 1.3774 |  1   |   0.001497   |  10.937541   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  26   | 1.4325 |  1   |   0.001464   |  11.216671   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  27   | 1.4876 |  1   |   0.001489   |  11.077631   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  28   | 1.5427 |  1   |   0.001487   |  11.029765   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  29   | 1.5978 |  1   |   0.001484   |  11.057365   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  30   | 1.6529 |  1   |   0.001499   |  10.988523   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  31   | 1.7080 |  1   |   0.001480   |  11.082585   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  32   | 1.7631 |  1   |   0.001493   |  11.026814   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  33   | 1.8182 |  1   |   0.001500   |  10.964227   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  34   | 1.8733 |  1   |   0.001496   |  11.070894   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  35   | 1.9284 |  1   |   0.001499   |  10.991359   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  36   | 1.9835 |  1   |   0.001493   |  10.956692   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  37   | 2.0386 |  1   |   0.001489   |  11.022953   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  38   | 2.0937 |  1   |   0.001493   |  11.002933   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  39   | 2.1488 |  1   |   0.001494   |  10.971854   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  40   | 2.2039 |  1   |   0.001499   |  10.956611   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  41   | 2.2590 |  1   |   0.001484   |  11.020658   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  42   | 2.3140 |  1   |   0.001484   |  11.016274   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  43   | 2.3691 |  1   |   0.001490   |  11.171444   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  44   | 2.4242 |  1   |   0.001496   |  10.975332   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  45   | 2.4793 |  1   |   0.001479   |  11.094840   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  46   | 2.5344 |  1   |   0.001484   |  11.060059   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  47   | 2.5895 |  1   |   0.001329   |  12.309849   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  48   | 2.6446 |  1   |   0.001455   |  11.269876   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  49   | 2.6997 |  1   |   0.001489   |  11.076657   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  50   | 2.7548 |  1   |   0.001496   |  10.998456   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  51   | 2.8099 |  1   |   0.001471   |  11.215534   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  52   | 2.8650 |  1   |   0.001494   |  10.981076   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  53   | 2.9201 |  1   |   0.001476   |  11.089418   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  54   | 2.9752 |  1   |   0.001482   |  10.973016   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  55   | 3.0303 |  1   |   0.001469   |  11.115966   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  56   | 3.0854 |  1   |   0.001458   |  11.172810   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  57   | 3.1405 |  1   |   0.001496   |  10.984580   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  58   | 3.1956 |  1   |   0.001451   |  11.196615   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  59   | 3.2507 |  1   |   0.001498   |  10.961545   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  60   | 3.3058 |  1   |   0.001492   |  11.053819   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  61   | 3.3609 |  1   |   0.001492   |  11.066293   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  62   | 3.4160 |  1   |   0.001489   |  11.025009   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  63   | 3.4711 |  1   |   0.001485   |  10.999920   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  64   | 3.5262 |  1   |   0.001505   |  10.973439   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  65   | 3.5813 |  1   |   0.001495   |  10.981763   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  66   | 3.6364 |  1   |   0.001481   |  10.997240   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  67   | 3.6915 |  1   |   0.001489   |  11.030007   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  68   | 3.7466 |  1   |   0.001487   |  10.977346   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  69   | 3.8017 |  1   |   0.001490   |  11.086072   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  70   | 3.8567 |  1   |   0.001478   |  11.027897   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  71   | 3.9118 |  1   |   0.001492   |  10.997195   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  72   | 3.9669 |  1   |   0.001476   |  11.129629   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  73   | 4.0220 |  1   |   0.001491   |  10.996616   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  74   | 4.0771 |  1   |   0.001513   |  10.996785   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  75   | 4.1322 |  1   |   0.001487   |  11.010984   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  76   | 4.1873 |  1   |   0.001504   |  10.960333   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  77   | 4.2424 |  1   |   0.001483   |  11.062480   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  78   | 4.2975 |  1   |   0.001502   |  10.999189   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  79   | 4.3526 |  1   |   0.001504   |  10.989805   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  80   | 4.4077 |  1   |   0.001499   |  10.986292   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  81   | 4.4628 |  1   |   0.001471   |  11.064248   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  82   | 4.5179 |  1   |   0.001478   |  11.030763   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  83   | 4.5730 |  1   |   0.001423   |  11.418318   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  84   | 4.6281 |  1   |   0.001513   |  10.934981   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  85   | 4.6832 |  1   |   0.001477   |  10.947891   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  86   | 4.7383 |  1   |   0.001449   |  11.337449   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  87   | 4.7934 |  1   |   0.001446   |  11.356113   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  88   | 4.8485 |  1   |   0.001493   |  11.006986   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  89   | 4.9036 |  1   |   0.001508   |  10.995290   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  90   | 4.9587 |  1   |   0.001519   |  10.961908   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  91   | 5.0138 |  1   |   0.001511   |  11.004637   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  92   | 5.0689 |  1   |   0.001493   |  11.060139   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  93   | 5.1240 |  1   |   0.001475   |  11.108975   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  94   | 5.1791 |  1   |   0.001501   |  10.946505   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  95   | 5.2342 |  1   |   0.001504   |  10.936386   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  96   | 5.2893 |  1   |   0.001486   |  10.994701   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  97   | 5.3444 |  1   |   0.001483   |  11.023811   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  98   | 5.3994 |  1   |   0.001494   |  10.979126   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  99   | 5.4545 |  1   |   0.001481   |  11.163518   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  100  | 5.5096 |  1   |   0.001479   |  11.047313   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  101  | 5.5647 |  1   |   0.001473   |  11.060119   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  102  | 5.6198 |  1   |   0.001486   |  11.053105   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  103  | 5.6749 |  1   |   0.001486   |  11.092742   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  104  | 5.7300 |  1   |   0.001477   |  11.137425   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  105  | 5.7851 |  1   |   0.001491   |  11.010371   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  106  | 5.8402 |  1   |   0.001496   |  10.990369   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  107  | 5.8953 |  1   |   0.001473   |  11.243003   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  108  | 5.9504 |  1   |   0.001388   |  11.929604   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  109  | 6.0055 |  1   |   0.001458   |  11.222341   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  110  | 6.0606 |  1   |   0.001499   |  11.021925   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  111  | 6.1157 |  1   |   0.001489   |  10.988726   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  112  | 6.1708 |  1   |   0.001464   |  11.127002   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  113  | 6.2259 |  1   |   0.001490   |  10.982098   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  114  | 6.2810 |  1   |   0.001490   |  11.058744   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  115  | 6.3361 |  1   |   0.001501   |  11.036917   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  116  | 6.3912 |  1   |   0.001487   |  11.093127   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  117  | 6.4463 |  1   |   0.001511   |  10.986041   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  118  | 6.5014 |  1   |   0.001473   |  11.067499   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  119  | 6.5565 |  1   |   0.001495   |  10.978576   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  120  | 6.6116 |  1   |   0.001497   |  11.027401   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  121  | 6.6667 |  1   |   0.001487   |  11.056618   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  122  | 6.7218 |  1   |   0.001492   |  10.946566   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  123  | 6.7769 |  1   |   0.001492   |  10.977270   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  124  | 6.8320 |  1   |   0.001466   |  11.212596   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  125  | 6.8871 |  1   |   0.001493   |  10.996106   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  126  | 6.9421 |  1   |   0.001480   |  11.018607   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  127  | 6.9972 |  1   |   0.001509   |  10.976320   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  128  | 7.0523 |  1   |   0.001329   |  12.417034   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  129  | 7.1074 |  1   |   0.001510   |  10.938081   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  130  | 7.1625 |  1   |   0.001496   |  10.964072   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  131  | 7.2176 |  1   |   0.001497   |  11.006588   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  132  | 7.2727 |  1   |   0.001493   |  11.018465   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  133  | 7.3278 |  1   |   0.001493   |  11.030553   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  134  | 7.3829 |  1   |   0.001490   |  11.034491   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  135  | 7.4380 |  1   |   0.001493   |  10.942718   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  136  | 7.4931 |  1   |   0.001485   |  11.063141   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  137  | 7.5482 |  1   |   0.001490   |  10.940797   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  138  | 7.6033 |  1   |   0.001496   |  11.094267   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  139  | 7.6584 |  1   |   0.001476   |  11.092305   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  140  | 7.7135 |  1   |   0.001489   |  10.999871   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  141  | 7.7686 |  1   |   0.001486   |  10.969563   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  142  | 7.8237 |  1   |   0.001485   |  11.045270   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  143  | 7.8788 |  1   |   0.001493   |  10.999192   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  144  | 7.9339 |  1   |   0.001498   |  11.051950   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  145  | 7.9890 |  1   |   0.001489   |  11.134410   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  146  | 8.0441 |  1   |   0.001474   |  11.101941   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  147  | 8.0992 |  1   |   0.001480   |  11.063681   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  148  | 8.1543 |  1   |   0.001486   |  10.946801   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  149  | 8.2094 |  1   |   0.001484   |  11.076774   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  150  | 8.2645 |  1   |   0.001480   |  11.143664   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  151  | 8.3196 |  1   |   0.001509   |  10.947056   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  152  | 8.3747 |  1   |   0.001447   |  11.226169   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  153  | 8.4298 |  1   |   0.001479   |  11.063754   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  154  | 8.4848 |  1   |   0.001484   |  10.961506   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  155  | 8.5399 |  1   |   0.001497   |  11.050401   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  156  | 8.5950 |  1   |   0.001503   |  11.018204   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  157  | 8.6501 |  1   |   0.001472   |  11.145519   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  158  | 8.7052 |  1   |   0.001486   |  11.037257   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  159  | 8.7603 |  1   |   0.001492   |  11.036405   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  160  | 8.8154 |  1   |   0.001480   |  11.015381   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  161  | 8.8705 |  1   |   0.001498   |  10.979052   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  162  | 8.9256 |  1   |   0.001496   |  11.051611   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  163  | 8.9807 |  1   |   0.001486   |  11.052557   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  164  | 9.0358 |  1   |   0.001491   |  11.034617   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  165  | 9.0909 |  1   |   0.001488   |  11.049273   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  166  | 9.1460 |  1   |   0.001491   |  11.085268   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  167  | 9.2011 |  1   |   0.001503   |  10.958939   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  168  | 9.2562 |  1   |   0.001497   |  11.062020   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  169  | 9.3113 |  1   |   0.001518   |  10.996845   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  170  | 9.3664 |  1   |   0.001470   |  11.209611   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  171  | 9.4215 |  1   |   0.001477   |  11.023668   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  172  | 9.4766 |  1   |   0.001481   |  11.090508   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  173  | 9.5317 |  1   |   0.001512   |  10.958598   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  174  | 9.5868 |  1   |   0.001472   |  11.093911   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  175  | 9.6419 |  1   |   0.001495   |  11.022321   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  176  | 9.6970 |  1   |   0.001497   |  10.963686   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  177  | 9.7521 |  1   |   0.001495   |  11.014677   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  178  | 9.8072 |  1   |   0.001494   |  10.971079   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  179  | 9.8623 |  1   |   0.001478   |  11.075295   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  180  | 9.9174 |  1   |   0.001492   |  10.967837   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  181  | 9.9725 |  1   |   0.001502   |  11.002055   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  182  |10.0275 |  1   |   0.001478   |  11.256881   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  183  |10.0826 |  1   |   0.001503   |  10.965836   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  184  |10.1377 |  1   |   0.001492   |  11.003513   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  185  |10.1928 |  1   |   0.001499   |  11.033890   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  186  |10.2479 |  1   |   0.001477   |  11.086396   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  187  |10.3030 |  1   |   0.001494   |  11.022237   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  188  |10.3581 |  1   |   0.001488   |  11.065405   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  189  |10.4132 |  1   |   0.001479   |  11.172747   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  190  |10.4683 |  1   |   0.001488   |  11.135500   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  191  |10.5234 |  1   |   0.001485   |  10.995766   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  192  |10.5785 |  1   |   0.001494   |  10.953639   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  193  |10.6336 |  1   |   0.001492   |  11.006385   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  194  |10.6887 |  1   |   0.001384   |  11.898180   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  195  |10.7438 |  1   |   0.001505   |  10.971906   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  196  |10.7989 |  1   |   0.001506   |  10.982670   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  197  |10.8540 |  1   |   0.001324   |  12.325860   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  198  |10.9091 |  1   |   0.001502   |  11.026861   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  199  |10.9642 |  1   |   0.001505   |  10.987199   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  200  |11.0193 |  1   |   0.001488   |  11.061289   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  201  |11.0744 |  1   |   0.001489   |  11.049199   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  202  |11.1295 |  1   |   0.001501   |  10.973437   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  203  |11.1846 |  1   |   0.001500   |  10.937448   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  204  |11.2397 |  1   |   0.001499   |  10.940185   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  205  |11.2948 |  1   |   0.001506   |  10.956871   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  206  |11.3499 |  1   |   0.001496   |  10.959225   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  207  |11.4050 |  1   |   0.001380   |  11.911324   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  208  |11.4601 |  1   |   0.001505   |  10.983367   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  209  |11.5152 |  1   |   0.001483   |  11.116040   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  210  |11.5702 |  1   |   0.001499   |  10.946843   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  211  |11.6253 |  1   |   0.001495   |  11.007702   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  212  |11.6804 |  1   |   0.001484   |  11.036537   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  213  |11.7355 |  1   |   0.001501   |  10.962850   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  214  |11.7906 |  1   |   0.001458   |  11.231765   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  215  |11.8457 |  1   |   0.001500   |  10.974388   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  216  |11.9008 |  1   |   0.001502   |  10.971881   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  217  |11.9559 |  1   |   0.001484   |  11.094227   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  218  |12.0110 |  1   |   0.001465   |  11.285769   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  219  |12.0661 |  1   |   0.001472   |  11.138884   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  220  |12.1212 |  1   |   0.001462   |  11.085258   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  221  |12.1763 |  1   |   0.001506   |  10.962956   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  222  |12.2314 |  1   |   0.001481   |  11.063845   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  223  |12.2865 |  1   |   0.001490   |  10.940562   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  224  |12.3416 |  1   |   0.001487   |  11.074080   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  225  |12.3967 |  1   |   0.001500   |  10.966088   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  226  |12.4518 |  1   |   0.001503   |  10.940730   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  227  |12.5069 |  1   |   0.001500   |  10.983767   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  228  |12.5620 |  1   |   0.001495   |  11.087621   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  229  |12.6171 |  1   |   0.001491   |  11.030768   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  230  |12.6722 |  1   |   0.001501   |  10.986019   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  231  |12.7273 |  1   |   0.001471   |  11.210606   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  232  |12.7824 |  1   |   0.001496   |  10.989540   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  233  |12.8375 |  1   |   0.001479   |  10.991125   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  234  |12.8926 |  1   |   0.001484   |  10.997315   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  235  |12.9477 |  1   |   0.001495   |  10.981242   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  236  |13.0028 |  1   |   0.001494   |  11.034151   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  237  |13.0579 |  1   |   0.001521   |  11.393696   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  238  |13.1129 |  1   |   0.001483   |  11.187673   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  239  |13.1680 |  1   |   0.001476   |  11.165030   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  240  |13.2231 |  1   |   0.001461   |  11.121420   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  241  |13.2782 |  1   |   0.001497   |  11.077922   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  242  |13.3333 |  1   |   0.001472   |  11.274707   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  243  |13.3884 |  1   |   0.001477   |  11.087600   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  244  |13.4435 |  1   |   0.001469   |  11.061153   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  245  |13.4986 |  1   |   0.001449   |  11.292755   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  246  |13.5537 |  1   |   0.001495   |  11.001245   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  247  |13.6088 |  1   |   0.001493   |  11.024726   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  248  |13.6639 |  1   |   0.001480   |  10.968597   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  249  |13.7190 |  1   |   0.001478   |  11.008988   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  250  |13.7741 |  1   |   0.001464   |  11.120604   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  251  |13.8292 |  1   |   0.001494   |  10.976703   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  252  |13.8843 |  1   |   0.001464   |  11.173856   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  253  |13.9394 |  1   |   0.001478   |  11.182188   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  254  |13.9945 |  1   |   0.001474   |  11.138903   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  255  |14.0496 |  1   |   0.001498   |  10.998538   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  256  |14.1047 |  1   |   0.001490   |  10.959735   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  257  |14.1598 |  1   |   0.001497   |  11.015843   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  258  |14.2149 |  1   |   0.001442   |  11.296236   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  259  |14.2700 |  1   |   0.001470   |  11.035043   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  260  |14.3251 |  1   |   0.001499   |  11.025067   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  261  |14.3802 |  1   |   0.001483   |  11.032883   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  262  |14.4353 |  1   |   0.001450   |  11.225693   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  263  |14.4904 |  1   |   0.001490   |  11.049877   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  264  |14.5455 |  1   |   0.001500   |  10.961913   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  265  |14.6006 |  1   |   0.001468   |  11.227395   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  266  |14.6556 |  1   |   0.001500   |  10.950821   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  267  |14.7107 |  1   |   0.001496   |  11.034819   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  268  |14.7658 |  1   |   0.001478   |  11.039118   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  269  |14.8209 |  1   |   0.001480   |  11.054862   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  270  |14.8760 |  1   |   0.001499   |  10.966170   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  271  |14.9311 |  1   |   0.001492   |  11.013549   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  272  |14.9862 |  1   |   0.001487   |  11.063011   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  273  |15.0413 |  1   |   0.001494   |  11.025787   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  274  |15.0964 |  1   |   0.001492   |  11.056323   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  275  |15.1515 |  1   |   0.001502   |  11.017958   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  276  |15.2066 |  1   |   0.001498   |  10.993303   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  277  |15.2617 |  1   |   0.001492   |  11.045624   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  278  |15.3168 |  1   |   0.001488   |  11.001706   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  279  |15.3719 |  1   |   0.001501   |  11.019076   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  280  |15.4270 |  1   |   0.001482   |  10.980229   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  281  |15.4821 |  1   |   0.001500   |  10.990641   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  282  |15.5372 |  1   |   0.001482   |  11.039660   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  283  |15.5923 |  1   |   0.001506   |  10.960431   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  284  |15.6474 |  1   |   0.001483   |  11.035762   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  285  |15.7025 |  1   |   0.001483   |  10.972536   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  286  |15.7576 |  1   |   0.001497   |  11.078446   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  287  |15.8127 |  1   |   0.001494   |  11.003910   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  288  |15.8678 |  1   |   0.001524   |  10.954572   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  289  |15.9229 |  1   |   0.001468   |  11.170064   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  290  |15.9780 |  1   |   0.001466   |  11.078832   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  291  |16.0331 |  1   |   0.001500   |  10.965774   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  292  |16.0882 |  1   |   0.001492   |  10.988034   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  293  |16.1433 |  1   |   0.001498   |  11.034559   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  294  |16.1983 |  1   |   0.001485   |  10.967131   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  295  |16.2534 |  1   |   0.001480   |  11.102506   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  296  |16.3085 |  1   |   0.001504   |  10.998471   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  297  |16.3636 |  1   |   0.001499   |  10.959420   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  298  |16.4187 |  1   |   0.001483   |  11.099791   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  299  |16.4738 |  1   |   0.001493   |  10.974991   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  300  |16.5289 |  1   |   0.001487   |  10.989717   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  301  |16.5840 |  1   |   0.001484   |  11.069229   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  302  |16.6391 |  1   |   0.001503   |  10.985806   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  303  |16.6942 |  1   |   0.001496   |  11.041107   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  304  |16.7493 |  1   |   0.001487   |  10.991144   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  305  |16.8044 |  1   |   0.001492   |  11.005035   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  306  |16.8595 |  1   |   0.001488   |  11.087163   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  307  |16.9146 |  1   |   0.001491   |  11.002862   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  308  |16.9697 |  1   |   0.001519   |  11.018252   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  309  |17.0248 |  1   |   0.001489   |  11.030986   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  310  |17.0799 |  1   |   0.001509   |  11.001859   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  311  |17.1350 |  1   |   0.001492   |  11.019300   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  312  |17.1901 |  1   |   0.001488   |  11.036436   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  313  |17.2452 |  1   |   0.001493   |  11.046922   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  314  |17.3003 |  1   |   0.001485   |  10.976898   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  315  |17.3554 |  1   |   0.001497   |  10.941875   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  316  |17.4105 |  1   |   0.001497   |  11.056621   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  317  |17.4656 |  1   |   0.001471   |  11.117604   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  318  |17.5207 |  1   |   0.001469   |  11.168160   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  319  |17.5758 |  1   |   0.001476   |  11.020628   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  320  |17.6309 |  1   |   0.001490   |  11.020211   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  321  |17.6860 |  1   |   0.001488   |  11.010496   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  322  |17.7410 |  1   |   0.001475   |  11.125369   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  323  |17.7961 |  1   |   0.001494   |  10.955690   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  324  |17.8512 |  1   |   0.001493   |  11.078510   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  325  |17.9063 |  1   |   0.001362   |  12.022863   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  326  |17.9614 |  1   |   0.001489   |  11.021870   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  327  |18.0165 |  1   |   0.001463   |  11.139775   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  328  |18.0716 |  1   |   0.001475   |  11.015395   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  329  |18.1267 |  1   |   0.001475   |  10.995628   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  330  |18.1818 |  1   |   0.001490   |  11.037015   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  331  |18.2369 |  1   |   0.001489   |  10.994272   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  332  |18.2920 |  1   |   0.001485   |  11.023925   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  333  |18.3471 |  1   |   0.001465   |  11.131615   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  334  |18.4022 |  1   |   0.001481   |  10.996215   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  335  |18.4573 |  1   |   0.001493   |  11.042847   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  336  |18.5124 |  1   |   0.001493   |  11.153712   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  337  |18.5675 |  1   |   0.001497   |  10.981005   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  338  |18.6226 |  1   |   0.001490   |  10.993124   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  339  |18.6777 |  1   |   0.001498   |  10.970379   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  340  |18.7328 |  1   |   0.001492   |  10.992345   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  341  |18.7879 |  1   |   0.001492   |  10.998224   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  342  |18.8430 |  1   |   0.001523   |  10.982035   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  343  |18.8981 |  1   |   0.001485   |  11.054384   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  344  |18.9532 |  1   |   0.001486   |  11.027520   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  345  |19.0083 |  1   |   0.001482   |  11.016985   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  346  |19.0634 |  1   |   0.001484   |  11.169659   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  347  |19.1185 |  1   |   0.001476   |  11.164135   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  348  |19.1736 |  1   |   0.001499   |  10.963815   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  349  |19.2287 |  1   |   0.001481   |  11.041115   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  350  |19.2837 |  1   |   0.001478   |  11.025983   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  351  |19.3388 |  1   |   0.001504   |  10.984140   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  352  |19.3939 |  1   |   0.001478   |  11.033161   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  353  |19.4490 |  1   |   0.001490   |  11.040692   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  354  |19.5041 |  1   |   0.001479   |  11.018249   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  355  |19.5592 |  1   |   0.001471   |  10.995813   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  356  |19.6143 |  1   |   0.001478   |  11.009744   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  357  |19.6694 |  1   |   0.001489   |  11.012060   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  358  |19.7245 |  1   |   0.001487   |  11.043458   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  359  |19.7796 |  1   |   0.001492   |  10.996924   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  360  |19.8347 |  1   |   0.001486   |  10.997914   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  361  |19.8898 |  1   |   0.001462   |  11.193126   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  362  |19.9449 |  1   |   0.001476   |  11.131790   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  363  |20.0000 |  1   |   0.001456   |  11.255347   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  364  |20.0551 |  1   |   0.001443   |  11.262457   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  365  |20.1102 |  1   |   0.001442   |  11.282331   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  366  |20.1653 |  1   |   0.001460   |  11.133270   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  367  |20.2204 |  1   |   0.001461   |  11.137590   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  368  |20.2755 |  1   |   0.001458   |  11.279287   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  369  |20.3306 |  1   |   0.001473   |  11.155934   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  370  |20.3857 |  1   |   0.001480   |  11.140440   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  371  |20.4408 |  1   |   0.001464   |  11.196868   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  372  |20.4959 |  1   |   0.001473   |  11.196430   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  373  |20.5510 |  1   |   0.001467   |  11.166742   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  374  |20.6061 |  1   |   0.001488   |  11.095840   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  375  |20.6612 |  1   |   0.001456   |  11.119288   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  376  |20.7163 |  1   |   0.001514   |  10.997535   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  377  |20.7713 |  1   |   0.001500   |  10.973636   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  378  |20.8264 |  1   |   0.001485   |  11.041631   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  379  |20.8815 |  1   |   0.001474   |  11.054704   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  380  |20.9366 |  1   |   0.001505   |  10.956344   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  381  |20.9917 |  1   |   0.001485   |  11.037928   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  382  |21.0468 |  1   |   0.001494   |  11.004072   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  383  |21.1019 |  1   |   0.001460   |  11.191711   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  384  |21.1570 |  1   |   0.001464   |  11.213192   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  385  |21.2121 |  1   |   0.001297   |  12.715977   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  386  |21.2672 |  1   |   0.001473   |  11.214244   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  387  |21.3223 |  1   |   0.001488   |  11.175076   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  388  |21.3774 |  1   |   0.001477   |  11.140157   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  389  |21.4325 |  1   |   0.001476   |  11.217773   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  390  |21.4876 |  1   |   0.001487   |  11.070076   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  391  |21.5427 |  1   |   0.001476   |  11.079977   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  392  |21.5978 |  1   |   0.001497   |  10.963420   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  393  |21.6529 |  1   |   0.001492   |  10.965270   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  394  |21.7080 |  1   |   0.001499   |  11.031488   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  395  |21.7631 |  1   |   0.001379   |  11.928026   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  396  |21.8182 |  1   |   0.001474   |  11.124487   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  397  |21.8733 |  1   |   0.001482   |  11.037149   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  398  |21.9284 |  1   |   0.001491   |  11.062619   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  399  |21.9835 |  1   |   0.001483   |  11.227740   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  400  |22.0386 |  1   |   0.001460   |  11.253708   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  401  |22.0937 |  1   |   0.001474   |  11.210907   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  402  |22.1488 |  1   |   0.001483   |  11.119041   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  403  |22.2039 |  1   |   0.001490   |  11.014906   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  404  |22.2590 |  1   |   0.001499   |  10.956471   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  405  |22.3140 |  1   |   0.001468   |  11.109550   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  406  |22.3691 |  1   |   0.001490   |  10.982853   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  407  |22.4242 |  1   |   0.001478   |  11.193154   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  408  |22.4793 |  1   |   0.001498   |  11.031459   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  409  |22.5344 |  1   |   0.001508   |  10.978516   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  410  |22.5895 |  1   |   0.001481   |  11.028180   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  411  |22.6446 |  1   |   0.001492   |  11.037275   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  412  |22.6997 |  1   |   0.001483   |  10.967960   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  413  |22.7548 |  1   |   0.001493   |  10.990911   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  414  |22.8099 |  1   |   0.001484   |  11.067482   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  415  |22.8650 |  1   |   0.001486   |  10.987234   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  416  |22.9201 |  1   |   0.001453   |  11.105183   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  417  |22.9752 |  1   |   0.001484   |  11.019325   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  418  |23.0303 |  1   |   0.001489   |  11.033711   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  419  |23.0854 |  1   |   0.001486   |  10.961752   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  420  |23.1405 |  1   |   0.001488   |  11.020402   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  421  |23.1956 |  1   |   0.001485   |  10.985927   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  422  |23.2507 |  1   |   0.001484   |  11.122986   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  423  |23.3058 |  1   |   0.001488   |  11.010117   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  424  |23.3609 |  1   |   0.001494   |  11.027276   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  425  |23.4160 |  1   |   0.001498   |  11.002853   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  426  |23.4711 |  1   |   0.001504   |  10.948009   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  427  |23.5262 |  1   |   0.001537   |  10.965587   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  428  |23.5813 |  1   |   0.001480   |  11.023513   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  429  |23.6364 |  1   |   0.001495   |  11.016430   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  430  |23.6915 |  1   |   0.001498   |  10.990407   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  431  |23.7466 |  1   |   0.001482   |  11.029547   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  432  |23.8017 |  1   |   0.001489   |  11.028503   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  433  |23.8567 |  1   |   0.001494   |  10.994521   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  434  |23.9118 |  1   |   0.001497   |  11.033860   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  435  |23.9669 |  1   |   0.001472   |  11.085127   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  436  |24.0220 |  1   |   0.001519   |  10.940756   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  437  |24.0771 |  1   |   0.001477   |  11.058848   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  438  |24.1322 |  1   |   0.001479   |  11.067911   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  439  |24.1873 |  1   |   0.001488   |  11.029004   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  440  |24.2424 |  1   |   0.001497   |  10.977903   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  441  |24.2975 |  1   |   0.001486   |  11.038936   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  442  |24.3526 |  1   |   0.001492   |  11.042946   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  443  |24.4077 |  1   |   0.001492   |  10.970813   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  444  |24.4628 |  1   |   0.001486   |  11.023437   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  445  |24.5179 |  1   |   0.001504   |  10.981979   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  446  |24.5730 |  1   |   0.001487   |  10.943550   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  447  |24.6281 |  1   |   0.001504   |  10.990585   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  448  |24.6832 |  1   |   0.001488   |  11.020756   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  449  |24.7383 |  1   |   0.001502   |  11.020719   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "|  450  |24.7934 |  1   |   0.001490   |  10.944653   |   0.000000   | 0.000000e+00 | 0.000000e+00 |\u001b[0m\n",
      "\u001b[34m...Finished\u001b[0m\n",
      "\u001b[36mFINISHED - Elapsed time = 5192.1956986 seconds\u001b[0m\n",
      "\u001b[36mFINISHED - CPU process time = 38238.0321676 seconds\u001b[0m\n"
     ]
    }
   ],
   "source": [
    "sharpy_output = sharpy.sharpy_main.main(['', SimInfo.solvers['SHARPy']['route'] + SimInfo.solvers['SHARPy']['case'] + '.sharpy'])"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Postprocessing"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "This reads the structural and aerodynamic information of the last time step."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {},
   "outputs": [],
   "source": [
    "tstep = sharpy_output.structure.timestep_info[-1]\n",
    "astep = sharpy_output.aero.timestep_info[-1]"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Now we separate the structure into blades:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {},
   "outputs": [],
   "source": [
    "# Define beams\n",
    "ielem = 0\n",
    "nblades = np.max(sharpy_output.structure.beam_number) + 1\n",
    "nodes_blade = []\n",
    "first_node = 0\n",
    "for iblade in range(nblades):\n",
    "    nodes_blade.append(np.zeros((sharpy_output.structure.num_node,), dtype=bool))\n",
    "    while sharpy_output.structure.beam_number[ielem] <= iblade:\n",
    "        ielem += 1\n",
    "        if ielem == sharpy_output.structure.num_elem:\n",
    "            break\n",
    "    nodes_blade[iblade][first_node:sharpy_output.structure.connectivities[ielem-1,1]+1] = True\n",
    "    first_node = sharpy_output.structure.connectivities[ielem-1,1]+1"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Compute the radial position of the nodes and initialise the rest of the variables."
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {},
   "outputs": [],
   "source": [
    "r = []\n",
    "c = []\n",
    "dr = []\n",
    "forces = []\n",
    "CN_drR = []\n",
    "CTan_drR = []\n",
    "CP_drR = []\n",
    "nodes_num = []\n",
    "for iblade in range(nblades):\n",
    "    forces.append(tstep.steady_applied_forces[nodes_blade[iblade]].copy())\n",
    "\n",
    "    nodes_num.append(np.arange(0, sharpy_output.structure.num_node, 1)[nodes_blade[iblade]])\n",
    "\n",
    "    r.append(np.linalg.norm(tstep.pos[nodes_blade[iblade], :], axis=1))\n",
    "    dr.append(np.zeros(np.sum(nodes_blade[iblade])))\n",
    "    dr[iblade][0] = 0.5*(r[iblade][1]-r[iblade][0])\n",
    "    dr[iblade][-1] = 0.5 * (r[iblade][-1] - r[iblade][-2])\n",
    "    for inode in range(1,len(r[iblade]) - 1):\n",
    "        dr[iblade][inode] = 0.5*(r[iblade][inode+1] - r[iblade][inode-1])\n",
    "\n",
    "    CN_drR.append(np.zeros(len(r[iblade])))\n",
    "    c.append(np.zeros(len(r[iblade])))\n",
    "    CTan_drR.append(np.zeros(len(r[iblade])))\n",
    "    CP_drR.append(np.zeros(len(r[iblade])))"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Transform the loads computed by SHARPy into out-of-plane and in-plane components:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [],
   "source": [
    "rho = sharpy_output.settings['StaticCoupled']['aero_solver_settings']['rho']\n",
    "uinf = sharpy_output.settings['StaticCoupled']['aero_solver_settings']['velocity_field_input']['u_inf']\n",
    "R = np.max(r[0])\n",
    "Cp = 0\n",
    "Ct = 0\n",
    "\n",
    "global_force_factor = 0.5 * rho * uinf** 2 * np.pi * R**2\n",
    "global_power_factor = global_force_factor*uinf\n",
    "for iblade in range(nblades):\n",
    "    for inode in range(len(r[iblade])):\n",
    "        forces[iblade][inode, 0] *= 0. # Discard the spanwise component\n",
    "\n",
    "        node_global_index = nodes_num[iblade][inode]\n",
    "        ielem = sharpy_output.structure.node_master_elem[node_global_index, 0]\n",
    "        inode_in_elem = sharpy_output.structure.node_master_elem[node_global_index, 1]\n",
    "        CAB = algebra.crv2rotation(tstep.psi[ielem, inode_in_elem, :])\n",
    "\n",
    "        c[iblade][inode] = sharpy_output.aero.data_dict['chord'][ielem,inode_in_elem]\n",
    "\n",
    "        forces_AFoR = np.dot(CAB, forces[iblade][inode, 0:3])\n",
    "\n",
    "        CN_drR[iblade][inode] = forces_AFoR[2]/dr[iblade][inode]*R / global_force_factor\n",
    "        CTan_drR[iblade][inode] = np.linalg.norm(forces_AFoR[0:2])/dr[iblade][inode]*R  / global_force_factor\n",
    "        CP_drR[iblade][inode] = np.linalg.norm(forces_AFoR[0:2])/dr[iblade][inode]*R  * r[iblade][inode]*rotation_velocity / global_power_factor\n",
    "\n",
    "    Cp += np.sum(CP_drR[iblade]*dr[iblade]/R)\n",
    "    Ct += np.sum(CN_drR[iblade]*dr[iblade]/R)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Results"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Plot of the loads along the blade:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\" standalone=\"no\"?>\n",
       "<!DOCTYPE svg PUBLIC \"-//W3C//DTD SVG 1.1//EN\"\n",
       "  \"http://www.w3.org/Graphics/SVG/1.1/DTD/svg11.dtd\">\n",
       "<!-- Created with matplotlib (https://matplotlib.org/) -->\n",
       "<svg height=\"207.83625pt\" version=\"1.1\" viewBox=\"0 0 720.58125 207.83625\" width=\"720.58125pt\" xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\">\n",
       " <defs>\n",
       "  <style type=\"text/css\">\n",
       "*{stroke-linecap:butt;stroke-linejoin:round;}\n",
       "  </style>\n",
       " </defs>\n",
       " <g id=\"figure_1\">\n",
       "  <g id=\"patch_1\">\n",
       "   <path d=\"M 0 207.83625 \n",
       "L 720.58125 207.83625 \n",
       "L 720.58125 0 \n",
       "L 0 0 \n",
       "z\n",
       "\" style=\"fill:none;\"/>\n",
       "  </g>\n",
       "  <g id=\"axes_1\">\n",
       "   <g id=\"patch_2\">\n",
       "    <path d=\"M 43.78125 170.28 \n",
       "L 348.144886 170.28 \n",
       "L 348.144886 7.2 \n",
       "L 43.78125 7.2 \n",
       "z\n",
       "\" style=\"fill:#ffffff;\"/>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_1\">\n",
       "    <g id=\"xtick_1\">\n",
       "     <g id=\"line2d_1\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 57.615961 170.28 \n",
       "L 57.615961 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_2\">\n",
       "      <defs>\n",
       "       <path d=\"M 0 0 \n",
       "L 0 3.5 \n",
       "\" id=\"m6e3d9edeb4\" style=\"stroke:#000000;stroke-width:0.8;\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"57.615961\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_1\">\n",
       "      <!-- 0.0 -->\n",
       "      <defs>\n",
       "       <path d=\"M 31.78125 66.40625 \n",
       "Q 24.171875 66.40625 20.328125 58.90625 \n",
       "Q 16.5 51.421875 16.5 36.375 \n",
       "Q 16.5 21.390625 20.328125 13.890625 \n",
       "Q 24.171875 6.390625 31.78125 6.390625 \n",
       "Q 39.453125 6.390625 43.28125 13.890625 \n",
       "Q 47.125 21.390625 47.125 36.375 \n",
       "Q 47.125 51.421875 43.28125 58.90625 \n",
       "Q 39.453125 66.40625 31.78125 66.40625 \n",
       "z\n",
       "M 31.78125 74.21875 \n",
       "Q 44.046875 74.21875 50.515625 64.515625 \n",
       "Q 56.984375 54.828125 56.984375 36.375 \n",
       "Q 56.984375 17.96875 50.515625 8.265625 \n",
       "Q 44.046875 -1.421875 31.78125 -1.421875 \n",
       "Q 19.53125 -1.421875 13.0625 8.265625 \n",
       "Q 6.59375 17.96875 6.59375 36.375 \n",
       "Q 6.59375 54.828125 13.0625 64.515625 \n",
       "Q 19.53125 74.21875 31.78125 74.21875 \n",
       "z\n",
       "\" id=\"DejaVuSans-48\"/>\n",
       "       <path d=\"M 10.6875 12.40625 \n",
       "L 21 12.40625 \n",
       "L 21 0 \n",
       "L 10.6875 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-46\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(49.664398 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_2\">\n",
       "     <g id=\"line2d_3\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 112.954804 170.28 \n",
       "L 112.954804 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_4\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"112.954804\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_2\">\n",
       "      <!-- 0.2 -->\n",
       "      <defs>\n",
       "       <path d=\"M 19.1875 8.296875 \n",
       "L 53.609375 8.296875 \n",
       "L 53.609375 0 \n",
       "L 7.328125 0 \n",
       "L 7.328125 8.296875 \n",
       "Q 12.9375 14.109375 22.625 23.890625 \n",
       "Q 32.328125 33.6875 34.8125 36.53125 \n",
       "Q 39.546875 41.84375 41.421875 45.53125 \n",
       "Q 43.3125 49.21875 43.3125 52.78125 \n",
       "Q 43.3125 58.59375 39.234375 62.25 \n",
       "Q 35.15625 65.921875 28.609375 65.921875 \n",
       "Q 23.96875 65.921875 18.8125 64.3125 \n",
       "Q 13.671875 62.703125 7.8125 59.421875 \n",
       "L 7.8125 69.390625 \n",
       "Q 13.765625 71.78125 18.9375 73 \n",
       "Q 24.125 74.21875 28.421875 74.21875 \n",
       "Q 39.75 74.21875 46.484375 68.546875 \n",
       "Q 53.21875 62.890625 53.21875 53.421875 \n",
       "Q 53.21875 48.921875 51.53125 44.890625 \n",
       "Q 49.859375 40.875 45.40625 35.40625 \n",
       "Q 44.1875 33.984375 37.640625 27.21875 \n",
       "Q 31.109375 20.453125 19.1875 8.296875 \n",
       "z\n",
       "\" id=\"DejaVuSans-50\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(105.003241 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-50\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_3\">\n",
       "     <g id=\"line2d_5\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 168.293647 170.28 \n",
       "L 168.293647 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_6\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"168.293647\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_3\">\n",
       "      <!-- 0.4 -->\n",
       "      <defs>\n",
       "       <path d=\"M 37.796875 64.3125 \n",
       "L 12.890625 25.390625 \n",
       "L 37.796875 25.390625 \n",
       "z\n",
       "M 35.203125 72.90625 \n",
       "L 47.609375 72.90625 \n",
       "L 47.609375 25.390625 \n",
       "L 58.015625 25.390625 \n",
       "L 58.015625 17.1875 \n",
       "L 47.609375 17.1875 \n",
       "L 47.609375 0 \n",
       "L 37.796875 0 \n",
       "L 37.796875 17.1875 \n",
       "L 4.890625 17.1875 \n",
       "L 4.890625 26.703125 \n",
       "z\n",
       "\" id=\"DejaVuSans-52\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(160.342084 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-52\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_4\">\n",
       "     <g id=\"line2d_7\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 223.63249 170.28 \n",
       "L 223.63249 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_8\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"223.63249\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_4\">\n",
       "      <!-- 0.6 -->\n",
       "      <defs>\n",
       "       <path d=\"M 33.015625 40.375 \n",
       "Q 26.375 40.375 22.484375 35.828125 \n",
       "Q 18.609375 31.296875 18.609375 23.390625 \n",
       "Q 18.609375 15.53125 22.484375 10.953125 \n",
       "Q 26.375 6.390625 33.015625 6.390625 \n",
       "Q 39.65625 6.390625 43.53125 10.953125 \n",
       "Q 47.40625 15.53125 47.40625 23.390625 \n",
       "Q 47.40625 31.296875 43.53125 35.828125 \n",
       "Q 39.65625 40.375 33.015625 40.375 \n",
       "z\n",
       "M 52.59375 71.296875 \n",
       "L 52.59375 62.3125 \n",
       "Q 48.875 64.0625 45.09375 64.984375 \n",
       "Q 41.3125 65.921875 37.59375 65.921875 \n",
       "Q 27.828125 65.921875 22.671875 59.328125 \n",
       "Q 17.53125 52.734375 16.796875 39.40625 \n",
       "Q 19.671875 43.65625 24.015625 45.921875 \n",
       "Q 28.375 48.1875 33.59375 48.1875 \n",
       "Q 44.578125 48.1875 50.953125 41.515625 \n",
       "Q 57.328125 34.859375 57.328125 23.390625 \n",
       "Q 57.328125 12.15625 50.6875 5.359375 \n",
       "Q 44.046875 -1.421875 33.015625 -1.421875 \n",
       "Q 20.359375 -1.421875 13.671875 8.265625 \n",
       "Q 6.984375 17.96875 6.984375 36.375 \n",
       "Q 6.984375 53.65625 15.1875 63.9375 \n",
       "Q 23.390625 74.21875 37.203125 74.21875 \n",
       "Q 40.921875 74.21875 44.703125 73.484375 \n",
       "Q 48.484375 72.75 52.59375 71.296875 \n",
       "z\n",
       "\" id=\"DejaVuSans-54\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(215.680927 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-54\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_5\">\n",
       "     <g id=\"line2d_9\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 278.971333 170.28 \n",
       "L 278.971333 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_10\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"278.971333\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_5\">\n",
       "      <!-- 0.8 -->\n",
       "      <defs>\n",
       "       <path d=\"M 31.78125 34.625 \n",
       "Q 24.75 34.625 20.71875 30.859375 \n",
       "Q 16.703125 27.09375 16.703125 20.515625 \n",
       "Q 16.703125 13.921875 20.71875 10.15625 \n",
       "Q 24.75 6.390625 31.78125 6.390625 \n",
       "Q 38.8125 6.390625 42.859375 10.171875 \n",
       "Q 46.921875 13.96875 46.921875 20.515625 \n",
       "Q 46.921875 27.09375 42.890625 30.859375 \n",
       "Q 38.875 34.625 31.78125 34.625 \n",
       "z\n",
       "M 21.921875 38.8125 \n",
       "Q 15.578125 40.375 12.03125 44.71875 \n",
       "Q 8.5 49.078125 8.5 55.328125 \n",
       "Q 8.5 64.0625 14.71875 69.140625 \n",
       "Q 20.953125 74.21875 31.78125 74.21875 \n",
       "Q 42.671875 74.21875 48.875 69.140625 \n",
       "Q 55.078125 64.0625 55.078125 55.328125 \n",
       "Q 55.078125 49.078125 51.53125 44.71875 \n",
       "Q 48 40.375 41.703125 38.8125 \n",
       "Q 48.828125 37.15625 52.796875 32.3125 \n",
       "Q 56.78125 27.484375 56.78125 20.515625 \n",
       "Q 56.78125 9.90625 50.3125 4.234375 \n",
       "Q 43.84375 -1.421875 31.78125 -1.421875 \n",
       "Q 19.734375 -1.421875 13.25 4.234375 \n",
       "Q 6.78125 9.90625 6.78125 20.515625 \n",
       "Q 6.78125 27.484375 10.78125 32.3125 \n",
       "Q 14.796875 37.15625 21.921875 38.8125 \n",
       "z\n",
       "M 18.3125 54.390625 \n",
       "Q 18.3125 48.734375 21.84375 45.5625 \n",
       "Q 25.390625 42.390625 31.78125 42.390625 \n",
       "Q 38.140625 42.390625 41.71875 45.5625 \n",
       "Q 45.3125 48.734375 45.3125 54.390625 \n",
       "Q 45.3125 60.0625 41.71875 63.234375 \n",
       "Q 38.140625 66.40625 31.78125 66.40625 \n",
       "Q 25.390625 66.40625 21.84375 63.234375 \n",
       "Q 18.3125 60.0625 18.3125 54.390625 \n",
       "z\n",
       "\" id=\"DejaVuSans-56\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(271.01977 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-56\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_6\">\n",
       "     <g id=\"line2d_11\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 334.310176 170.28 \n",
       "L 334.310176 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_12\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"334.310176\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_6\">\n",
       "      <!-- 1.0 -->\n",
       "      <defs>\n",
       "       <path d=\"M 12.40625 8.296875 \n",
       "L 28.515625 8.296875 \n",
       "L 28.515625 63.921875 \n",
       "L 10.984375 60.40625 \n",
       "L 10.984375 69.390625 \n",
       "L 28.421875 72.90625 \n",
       "L 38.28125 72.90625 \n",
       "L 38.28125 8.296875 \n",
       "L 54.390625 8.296875 \n",
       "L 54.390625 0 \n",
       "L 12.40625 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-49\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(326.358613 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-49\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_7\">\n",
       "     <!-- r/R [-] -->\n",
       "     <defs>\n",
       "      <path d=\"M 41.109375 46.296875 \n",
       "Q 39.59375 47.171875 37.8125 47.578125 \n",
       "Q 36.03125 48 33.890625 48 \n",
       "Q 26.265625 48 22.1875 43.046875 \n",
       "Q 18.109375 38.09375 18.109375 28.8125 \n",
       "L 18.109375 0 \n",
       "L 9.078125 0 \n",
       "L 9.078125 54.6875 \n",
       "L 18.109375 54.6875 \n",
       "L 18.109375 46.1875 \n",
       "Q 20.953125 51.171875 25.484375 53.578125 \n",
       "Q 30.03125 56 36.53125 56 \n",
       "Q 37.453125 56 38.578125 55.875 \n",
       "Q 39.703125 55.765625 41.0625 55.515625 \n",
       "z\n",
       "\" id=\"DejaVuSans-114\"/>\n",
       "      <path d=\"M 25.390625 72.90625 \n",
       "L 33.6875 72.90625 \n",
       "L 8.296875 -9.28125 \n",
       "L 0 -9.28125 \n",
       "z\n",
       "\" id=\"DejaVuSans-47\"/>\n",
       "      <path d=\"M 44.390625 34.1875 \n",
       "Q 47.5625 33.109375 50.5625 29.59375 \n",
       "Q 53.5625 26.078125 56.59375 19.921875 \n",
       "L 66.609375 0 \n",
       "L 56 0 \n",
       "L 46.6875 18.703125 \n",
       "Q 43.0625 26.03125 39.671875 28.421875 \n",
       "Q 36.28125 30.8125 30.421875 30.8125 \n",
       "L 19.671875 30.8125 \n",
       "L 19.671875 0 \n",
       "L 9.8125 0 \n",
       "L 9.8125 72.90625 \n",
       "L 32.078125 72.90625 \n",
       "Q 44.578125 72.90625 50.734375 67.671875 \n",
       "Q 56.890625 62.453125 56.890625 51.90625 \n",
       "Q 56.890625 45.015625 53.6875 40.46875 \n",
       "Q 50.484375 35.9375 44.390625 34.1875 \n",
       "z\n",
       "M 19.671875 64.796875 \n",
       "L 19.671875 38.921875 \n",
       "L 32.078125 38.921875 \n",
       "Q 39.203125 38.921875 42.84375 42.21875 \n",
       "Q 46.484375 45.515625 46.484375 51.90625 \n",
       "Q 46.484375 58.296875 42.84375 61.546875 \n",
       "Q 39.203125 64.796875 32.078125 64.796875 \n",
       "z\n",
       "\" id=\"DejaVuSans-82\"/>\n",
       "      <path id=\"DejaVuSans-32\"/>\n",
       "      <path d=\"M 8.59375 75.984375 \n",
       "L 29.296875 75.984375 \n",
       "L 29.296875 69 \n",
       "L 17.578125 69 \n",
       "L 17.578125 -6.203125 \n",
       "L 29.296875 -6.203125 \n",
       "L 29.296875 -13.1875 \n",
       "L 8.59375 -13.1875 \n",
       "z\n",
       "\" id=\"DejaVuSans-91\"/>\n",
       "      <path d=\"M 4.890625 31.390625 \n",
       "L 31.203125 31.390625 \n",
       "L 31.203125 23.390625 \n",
       "L 4.890625 23.390625 \n",
       "z\n",
       "\" id=\"DejaVuSans-45\"/>\n",
       "      <path d=\"M 30.421875 75.984375 \n",
       "L 30.421875 -13.1875 \n",
       "L 9.71875 -13.1875 \n",
       "L 9.71875 -6.203125 \n",
       "L 21.390625 -6.203125 \n",
       "L 21.390625 69 \n",
       "L 9.71875 69 \n",
       "L 9.71875 75.984375 \n",
       "z\n",
       "\" id=\"DejaVuSans-93\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(181.454474 198.556563)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-114\"/>\n",
       "      <use x=\"41.113281\" xlink:href=\"#DejaVuSans-47\"/>\n",
       "      <use x=\"74.804688\" xlink:href=\"#DejaVuSans-82\"/>\n",
       "      <use x=\"144.287109\" xlink:href=\"#DejaVuSans-32\"/>\n",
       "      <use x=\"176.074219\" xlink:href=\"#DejaVuSans-91\"/>\n",
       "      <use x=\"215.087891\" xlink:href=\"#DejaVuSans-45\"/>\n",
       "      <use x=\"251.171875\" xlink:href=\"#DejaVuSans-93\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_2\">\n",
       "    <g id=\"ytick_1\">\n",
       "     <g id=\"line2d_13\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 43.78125 162.867273 \n",
       "L 348.144886 162.867273 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_14\">\n",
       "      <defs>\n",
       "       <path d=\"M 0 0 \n",
       "L -3.5 0 \n",
       "\" id=\"m60ac5f1668\" style=\"stroke:#000000;stroke-width:0.8;\"/>\n",
       "      </defs>\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"43.78125\" xlink:href=\"#m60ac5f1668\" y=\"162.867273\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_8\">\n",
       "      <!-- 0.0 -->\n",
       "      <g transform=\"translate(20.878125 166.666491)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_2\">\n",
       "     <g id=\"line2d_15\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 43.78125 129.747022 \n",
       "L 348.144886 129.747022 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_16\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"43.78125\" xlink:href=\"#m60ac5f1668\" y=\"129.747022\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_9\">\n",
       "      <!-- 0.1 -->\n",
       "      <g transform=\"translate(20.878125 133.546241)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-49\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_3\">\n",
       "     <g id=\"line2d_17\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 43.78125 96.626771 \n",
       "L 348.144886 96.626771 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_18\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"43.78125\" xlink:href=\"#m60ac5f1668\" y=\"96.626771\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_10\">\n",
       "      <!-- 0.2 -->\n",
       "      <g transform=\"translate(20.878125 100.42599)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-50\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_4\">\n",
       "     <g id=\"line2d_19\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 43.78125 63.506521 \n",
       "L 348.144886 63.506521 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_20\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"43.78125\" xlink:href=\"#m60ac5f1668\" y=\"63.506521\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_11\">\n",
       "      <!-- 0.3 -->\n",
       "      <defs>\n",
       "       <path d=\"M 40.578125 39.3125 \n",
       "Q 47.65625 37.796875 51.625 33 \n",
       "Q 55.609375 28.21875 55.609375 21.1875 \n",
       "Q 55.609375 10.40625 48.1875 4.484375 \n",
       "Q 40.765625 -1.421875 27.09375 -1.421875 \n",
       "Q 22.515625 -1.421875 17.65625 -0.515625 \n",
       "Q 12.796875 0.390625 7.625 2.203125 \n",
       "L 7.625 11.71875 \n",
       "Q 11.71875 9.328125 16.59375 8.109375 \n",
       "Q 21.484375 6.890625 26.8125 6.890625 \n",
       "Q 36.078125 6.890625 40.9375 10.546875 \n",
       "Q 45.796875 14.203125 45.796875 21.1875 \n",
       "Q 45.796875 27.640625 41.28125 31.265625 \n",
       "Q 36.765625 34.90625 28.71875 34.90625 \n",
       "L 20.21875 34.90625 \n",
       "L 20.21875 43.015625 \n",
       "L 29.109375 43.015625 \n",
       "Q 36.375 43.015625 40.234375 45.921875 \n",
       "Q 44.09375 48.828125 44.09375 54.296875 \n",
       "Q 44.09375 59.90625 40.109375 62.90625 \n",
       "Q 36.140625 65.921875 28.71875 65.921875 \n",
       "Q 24.65625 65.921875 20.015625 65.03125 \n",
       "Q 15.375 64.15625 9.8125 62.3125 \n",
       "L 9.8125 71.09375 \n",
       "Q 15.4375 72.65625 20.34375 73.4375 \n",
       "Q 25.25 74.21875 29.59375 74.21875 \n",
       "Q 40.828125 74.21875 47.359375 69.109375 \n",
       "Q 53.90625 64.015625 53.90625 55.328125 \n",
       "Q 53.90625 49.265625 50.4375 45.09375 \n",
       "Q 46.96875 40.921875 40.578125 39.3125 \n",
       "z\n",
       "\" id=\"DejaVuSans-51\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(20.878125 67.305739)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-51\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_5\">\n",
       "     <g id=\"line2d_21\">\n",
       "      <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 43.78125 30.38627 \n",
       "L 348.144886 30.38627 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_22\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"43.78125\" xlink:href=\"#m60ac5f1668\" y=\"30.38627\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_12\">\n",
       "      <!-- 0.4 -->\n",
       "      <g transform=\"translate(20.878125 34.185489)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-52\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_13\">\n",
       "     <!-- CN/d(r/R) [-] -->\n",
       "     <defs>\n",
       "      <path d=\"M 64.40625 67.28125 \n",
       "L 64.40625 56.890625 \n",
       "Q 59.421875 61.53125 53.78125 63.8125 \n",
       "Q 48.140625 66.109375 41.796875 66.109375 \n",
       "Q 29.296875 66.109375 22.65625 58.46875 \n",
       "Q 16.015625 50.828125 16.015625 36.375 \n",
       "Q 16.015625 21.96875 22.65625 14.328125 \n",
       "Q 29.296875 6.6875 41.796875 6.6875 \n",
       "Q 48.140625 6.6875 53.78125 8.984375 \n",
       "Q 59.421875 11.28125 64.40625 15.921875 \n",
       "L 64.40625 5.609375 \n",
       "Q 59.234375 2.09375 53.4375 0.328125 \n",
       "Q 47.65625 -1.421875 41.21875 -1.421875 \n",
       "Q 24.65625 -1.421875 15.125 8.703125 \n",
       "Q 5.609375 18.84375 5.609375 36.375 \n",
       "Q 5.609375 53.953125 15.125 64.078125 \n",
       "Q 24.65625 74.21875 41.21875 74.21875 \n",
       "Q 47.75 74.21875 53.53125 72.484375 \n",
       "Q 59.328125 70.75 64.40625 67.28125 \n",
       "z\n",
       "\" id=\"DejaVuSans-67\"/>\n",
       "      <path d=\"M 9.8125 72.90625 \n",
       "L 23.09375 72.90625 \n",
       "L 55.421875 11.921875 \n",
       "L 55.421875 72.90625 \n",
       "L 64.984375 72.90625 \n",
       "L 64.984375 0 \n",
       "L 51.703125 0 \n",
       "L 19.390625 60.984375 \n",
       "L 19.390625 0 \n",
       "L 9.8125 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-78\"/>\n",
       "      <path d=\"M 45.40625 46.390625 \n",
       "L 45.40625 75.984375 \n",
       "L 54.390625 75.984375 \n",
       "L 54.390625 0 \n",
       "L 45.40625 0 \n",
       "L 45.40625 8.203125 \n",
       "Q 42.578125 3.328125 38.25 0.953125 \n",
       "Q 33.9375 -1.421875 27.875 -1.421875 \n",
       "Q 17.96875 -1.421875 11.734375 6.484375 \n",
       "Q 5.515625 14.40625 5.515625 27.296875 \n",
       "Q 5.515625 40.1875 11.734375 48.09375 \n",
       "Q 17.96875 56 27.875 56 \n",
       "Q 33.9375 56 38.25 53.625 \n",
       "Q 42.578125 51.265625 45.40625 46.390625 \n",
       "z\n",
       "M 14.796875 27.296875 \n",
       "Q 14.796875 17.390625 18.875 11.75 \n",
       "Q 22.953125 6.109375 30.078125 6.109375 \n",
       "Q 37.203125 6.109375 41.296875 11.75 \n",
       "Q 45.40625 17.390625 45.40625 27.296875 \n",
       "Q 45.40625 37.203125 41.296875 42.84375 \n",
       "Q 37.203125 48.484375 30.078125 48.484375 \n",
       "Q 22.953125 48.484375 18.875 42.84375 \n",
       "Q 14.796875 37.203125 14.796875 27.296875 \n",
       "z\n",
       "\" id=\"DejaVuSans-100\"/>\n",
       "      <path d=\"M 31 75.875 \n",
       "Q 24.46875 64.65625 21.28125 53.65625 \n",
       "Q 18.109375 42.671875 18.109375 31.390625 \n",
       "Q 18.109375 20.125 21.3125 9.0625 \n",
       "Q 24.515625 -2 31 -13.1875 \n",
       "L 23.1875 -13.1875 \n",
       "Q 15.875 -1.703125 12.234375 9.375 \n",
       "Q 8.59375 20.453125 8.59375 31.390625 \n",
       "Q 8.59375 42.28125 12.203125 53.3125 \n",
       "Q 15.828125 64.359375 23.1875 75.875 \n",
       "z\n",
       "\" id=\"DejaVuSans-40\"/>\n",
       "      <path d=\"M 8.015625 75.875 \n",
       "L 15.828125 75.875 \n",
       "Q 23.140625 64.359375 26.78125 53.3125 \n",
       "Q 30.421875 42.28125 30.421875 31.390625 \n",
       "Q 30.421875 20.453125 26.78125 9.375 \n",
       "Q 23.140625 -1.703125 15.828125 -13.1875 \n",
       "L 8.015625 -13.1875 \n",
       "Q 14.5 -2 17.703125 9.0625 \n",
       "Q 20.90625 20.125 20.90625 31.390625 \n",
       "Q 20.90625 42.671875 17.703125 53.65625 \n",
       "Q 14.5 64.65625 8.015625 75.875 \n",
       "z\n",
       "\" id=\"DejaVuSans-41\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(14.798438 119.240781)rotate(-90)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-67\"/>\n",
       "      <use x=\"69.824219\" xlink:href=\"#DejaVuSans-78\"/>\n",
       "      <use x=\"144.628906\" xlink:href=\"#DejaVuSans-47\"/>\n",
       "      <use x=\"178.320312\" xlink:href=\"#DejaVuSans-100\"/>\n",
       "      <use x=\"241.796875\" xlink:href=\"#DejaVuSans-40\"/>\n",
       "      <use x=\"280.810547\" xlink:href=\"#DejaVuSans-114\"/>\n",
       "      <use x=\"321.923828\" xlink:href=\"#DejaVuSans-47\"/>\n",
       "      <use x=\"355.615234\" xlink:href=\"#DejaVuSans-82\"/>\n",
       "      <use x=\"425.097656\" xlink:href=\"#DejaVuSans-41\"/>\n",
       "      <use x=\"464.111328\" xlink:href=\"#DejaVuSans-32\"/>\n",
       "      <use x=\"495.898438\" xlink:href=\"#DejaVuSans-91\"/>\n",
       "      <use x=\"534.912109\" xlink:href=\"#DejaVuSans-45\"/>\n",
       "      <use x=\"570.996094\" xlink:href=\"#DejaVuSans-93\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"line2d_23\">\n",
       "    <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 57.615961 162.867273 \n",
       "L 73.130329 162.867273 \n",
       "L 88.595884 151.198309 \n",
       "L 103.963964 136.634632 \n",
       "L 119.186216 128.573349 \n",
       "L 126.789514 118.267999 \n",
       "L 149.002266 112.315897 \n",
       "L 163.502253 100.781432 \n",
       "L 182.128357 93.633305 \n",
       "L 191.458178 83.798905 \n",
       "L 204.826158 76.740354 \n",
       "L 217.730961 68.28959 \n",
       "L 230.131982 59.797073 \n",
       "L 241.990204 51.604997 \n",
       "L 253.268316 43.766327 \n",
       "L 265.136622 37.507221 \n",
       "L 273.944209 31.387125 \n",
       "L 283.276935 26.902704 \n",
       "L 291.899648 23.217086 \n",
       "L 299.785218 21.080717 \n",
       "L 306.908834 21.163383 \n",
       "L 313.248083 24.028153 \n",
       "L 317.708523 32.386809 \n",
       "L 323.496225 42.364374 \n",
       "L 327.372874 57.246164 \n",
       "L 330.400768 71.982645 \n",
       "L 332.57038 89.376094 \n",
       "L 333.874884 109.07552 \n",
       "L 334.310176 133.385298 \n",
       "\" style=\"fill:none;stroke:#1f77b4;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_24\">\n",
       "    <path clip-path=\"url(#p19e32ab1d1)\" d=\"M 113.392966 134.313 \n",
       "L 144.141866 114.220659 \n",
       "L 170.490732 98.795658 \n",
       "L 223.193874 65.403284 \n",
       "L 249.545429 43.622354 \n",
       "L 293.4647 17.970932 \n",
       "L 306.640477 14.612727 \n",
       "L 322.012244 34.230749 \n",
       "L 330.796087 64.226962 \n",
       "L 334.310176 162.867273 \n",
       "\" style=\"fill:none;stroke:#ff7f0e;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_3\">\n",
       "    <path d=\"M 43.78125 170.28 \n",
       "L 43.78125 7.2 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_4\">\n",
       "    <path d=\"M 348.144886 170.28 \n",
       "L 348.144886 7.2 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_5\">\n",
       "    <path d=\"M 43.78125 170.28 \n",
       "L 348.144886 170.28 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_6\">\n",
       "    <path d=\"M 43.78125 7.2 \n",
       "L 348.144886 7.2 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"legend_1\">\n",
       "    <g id=\"patch_7\">\n",
       "     <path d=\"M 50.78125 44.55625 \n",
       "L 134.526563 44.55625 \n",
       "Q 136.526563 44.55625 136.526563 42.55625 \n",
       "L 136.526563 14.2 \n",
       "Q 136.526563 12.2 134.526563 12.2 \n",
       "L 50.78125 12.2 \n",
       "Q 48.78125 12.2 48.78125 14.2 \n",
       "L 48.78125 42.55625 \n",
       "Q 48.78125 44.55625 50.78125 44.55625 \n",
       "z\n",
       "\" style=\"fill:#ffffff;opacity:0.8;stroke:#cccccc;stroke-linejoin:miter;\"/>\n",
       "    </g>\n",
       "    <g id=\"line2d_25\">\n",
       "     <path d=\"M 52.78125 20.298437 \n",
       "L 72.78125 20.298437 \n",
       "\" style=\"fill:none;stroke:#1f77b4;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "    </g>\n",
       "    <g id=\"line2d_26\"/>\n",
       "    <g id=\"text_14\">\n",
       "     <!-- SHARPy -->\n",
       "     <defs>\n",
       "      <path d=\"M 53.515625 70.515625 \n",
       "L 53.515625 60.890625 \n",
       "Q 47.90625 63.578125 42.921875 64.890625 \n",
       "Q 37.9375 66.21875 33.296875 66.21875 \n",
       "Q 25.25 66.21875 20.875 63.09375 \n",
       "Q 16.5 59.96875 16.5 54.203125 \n",
       "Q 16.5 49.359375 19.40625 46.890625 \n",
       "Q 22.3125 44.4375 30.421875 42.921875 \n",
       "L 36.375 41.703125 \n",
       "Q 47.40625 39.59375 52.65625 34.296875 \n",
       "Q 57.90625 29 57.90625 20.125 \n",
       "Q 57.90625 9.515625 50.796875 4.046875 \n",
       "Q 43.703125 -1.421875 29.984375 -1.421875 \n",
       "Q 24.8125 -1.421875 18.96875 -0.25 \n",
       "Q 13.140625 0.921875 6.890625 3.21875 \n",
       "L 6.890625 13.375 \n",
       "Q 12.890625 10.015625 18.65625 8.296875 \n",
       "Q 24.421875 6.59375 29.984375 6.59375 \n",
       "Q 38.421875 6.59375 43.015625 9.90625 \n",
       "Q 47.609375 13.234375 47.609375 19.390625 \n",
       "Q 47.609375 24.75 44.3125 27.78125 \n",
       "Q 41.015625 30.8125 33.5 32.328125 \n",
       "L 27.484375 33.5 \n",
       "Q 16.453125 35.6875 11.515625 40.375 \n",
       "Q 6.59375 45.0625 6.59375 53.421875 \n",
       "Q 6.59375 63.09375 13.40625 68.65625 \n",
       "Q 20.21875 74.21875 32.171875 74.21875 \n",
       "Q 37.3125 74.21875 42.625 73.28125 \n",
       "Q 47.953125 72.359375 53.515625 70.515625 \n",
       "z\n",
       "\" id=\"DejaVuSans-83\"/>\n",
       "      <path d=\"M 9.8125 72.90625 \n",
       "L 19.671875 72.90625 \n",
       "L 19.671875 43.015625 \n",
       "L 55.515625 43.015625 \n",
       "L 55.515625 72.90625 \n",
       "L 65.375 72.90625 \n",
       "L 65.375 0 \n",
       "L 55.515625 0 \n",
       "L 55.515625 34.71875 \n",
       "L 19.671875 34.71875 \n",
       "L 19.671875 0 \n",
       "L 9.8125 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-72\"/>\n",
       "      <path d=\"M 34.1875 63.1875 \n",
       "L 20.796875 26.90625 \n",
       "L 47.609375 26.90625 \n",
       "z\n",
       "M 28.609375 72.90625 \n",
       "L 39.796875 72.90625 \n",
       "L 67.578125 0 \n",
       "L 57.328125 0 \n",
       "L 50.6875 18.703125 \n",
       "L 17.828125 18.703125 \n",
       "L 11.1875 0 \n",
       "L 0.78125 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-65\"/>\n",
       "      <path d=\"M 19.671875 64.796875 \n",
       "L 19.671875 37.40625 \n",
       "L 32.078125 37.40625 \n",
       "Q 38.96875 37.40625 42.71875 40.96875 \n",
       "Q 46.484375 44.53125 46.484375 51.125 \n",
       "Q 46.484375 57.671875 42.71875 61.234375 \n",
       "Q 38.96875 64.796875 32.078125 64.796875 \n",
       "z\n",
       "M 9.8125 72.90625 \n",
       "L 32.078125 72.90625 \n",
       "Q 44.34375 72.90625 50.609375 67.359375 \n",
       "Q 56.890625 61.8125 56.890625 51.125 \n",
       "Q 56.890625 40.328125 50.609375 34.8125 \n",
       "Q 44.34375 29.296875 32.078125 29.296875 \n",
       "L 19.671875 29.296875 \n",
       "L 19.671875 0 \n",
       "L 9.8125 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-80\"/>\n",
       "      <path d=\"M 32.171875 -5.078125 \n",
       "Q 28.375 -14.84375 24.75 -17.8125 \n",
       "Q 21.140625 -20.796875 15.09375 -20.796875 \n",
       "L 7.90625 -20.796875 \n",
       "L 7.90625 -13.28125 \n",
       "L 13.1875 -13.28125 \n",
       "Q 16.890625 -13.28125 18.9375 -11.515625 \n",
       "Q 21 -9.765625 23.484375 -3.21875 \n",
       "L 25.09375 0.875 \n",
       "L 2.984375 54.6875 \n",
       "L 12.5 54.6875 \n",
       "L 29.59375 11.921875 \n",
       "L 46.6875 54.6875 \n",
       "L 56.203125 54.6875 \n",
       "z\n",
       "\" id=\"DejaVuSans-121\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(80.78125 23.798437)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-83\"/>\n",
       "      <use x=\"63.476562\" xlink:href=\"#DejaVuSans-72\"/>\n",
       "      <use x=\"138.671875\" xlink:href=\"#DejaVuSans-65\"/>\n",
       "      <use x=\"207.080078\" xlink:href=\"#DejaVuSans-82\"/>\n",
       "      <use x=\"276.5625\" xlink:href=\"#DejaVuSans-80\"/>\n",
       "      <use x=\"336.865234\" xlink:href=\"#DejaVuSans-121\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"line2d_27\">\n",
       "     <path d=\"M 52.78125 34.976562 \n",
       "L 72.78125 34.976562 \n",
       "\" style=\"fill:none;stroke:#ff7f0e;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "    </g>\n",
       "    <g id=\"line2d_28\"/>\n",
       "    <g id=\"text_15\">\n",
       "     <!-- OpenFAST -->\n",
       "     <defs>\n",
       "      <path d=\"M 39.40625 66.21875 \n",
       "Q 28.65625 66.21875 22.328125 58.203125 \n",
       "Q 16.015625 50.203125 16.015625 36.375 \n",
       "Q 16.015625 22.609375 22.328125 14.59375 \n",
       "Q 28.65625 6.59375 39.40625 6.59375 \n",
       "Q 50.140625 6.59375 56.421875 14.59375 \n",
       "Q 62.703125 22.609375 62.703125 36.375 \n",
       "Q 62.703125 50.203125 56.421875 58.203125 \n",
       "Q 50.140625 66.21875 39.40625 66.21875 \n",
       "z\n",
       "M 39.40625 74.21875 \n",
       "Q 54.734375 74.21875 63.90625 63.9375 \n",
       "Q 73.09375 53.65625 73.09375 36.375 \n",
       "Q 73.09375 19.140625 63.90625 8.859375 \n",
       "Q 54.734375 -1.421875 39.40625 -1.421875 \n",
       "Q 24.03125 -1.421875 14.8125 8.828125 \n",
       "Q 5.609375 19.09375 5.609375 36.375 \n",
       "Q 5.609375 53.65625 14.8125 63.9375 \n",
       "Q 24.03125 74.21875 39.40625 74.21875 \n",
       "z\n",
       "\" id=\"DejaVuSans-79\"/>\n",
       "      <path d=\"M 18.109375 8.203125 \n",
       "L 18.109375 -20.796875 \n",
       "L 9.078125 -20.796875 \n",
       "L 9.078125 54.6875 \n",
       "L 18.109375 54.6875 \n",
       "L 18.109375 46.390625 \n",
       "Q 20.953125 51.265625 25.265625 53.625 \n",
       "Q 29.59375 56 35.59375 56 \n",
       "Q 45.5625 56 51.78125 48.09375 \n",
       "Q 58.015625 40.1875 58.015625 27.296875 \n",
       "Q 58.015625 14.40625 51.78125 6.484375 \n",
       "Q 45.5625 -1.421875 35.59375 -1.421875 \n",
       "Q 29.59375 -1.421875 25.265625 0.953125 \n",
       "Q 20.953125 3.328125 18.109375 8.203125 \n",
       "z\n",
       "M 48.6875 27.296875 \n",
       "Q 48.6875 37.203125 44.609375 42.84375 \n",
       "Q 40.53125 48.484375 33.40625 48.484375 \n",
       "Q 26.265625 48.484375 22.1875 42.84375 \n",
       "Q 18.109375 37.203125 18.109375 27.296875 \n",
       "Q 18.109375 17.390625 22.1875 11.75 \n",
       "Q 26.265625 6.109375 33.40625 6.109375 \n",
       "Q 40.53125 6.109375 44.609375 11.75 \n",
       "Q 48.6875 17.390625 48.6875 27.296875 \n",
       "z\n",
       "\" id=\"DejaVuSans-112\"/>\n",
       "      <path d=\"M 56.203125 29.59375 \n",
       "L 56.203125 25.203125 \n",
       "L 14.890625 25.203125 \n",
       "Q 15.484375 15.921875 20.484375 11.0625 \n",
       "Q 25.484375 6.203125 34.421875 6.203125 \n",
       "Q 39.59375 6.203125 44.453125 7.46875 \n",
       "Q 49.3125 8.734375 54.109375 11.28125 \n",
       "L 54.109375 2.78125 \n",
       "Q 49.265625 0.734375 44.1875 -0.34375 \n",
       "Q 39.109375 -1.421875 33.890625 -1.421875 \n",
       "Q 20.796875 -1.421875 13.15625 6.1875 \n",
       "Q 5.515625 13.8125 5.515625 26.8125 \n",
       "Q 5.515625 40.234375 12.765625 48.109375 \n",
       "Q 20.015625 56 32.328125 56 \n",
       "Q 43.359375 56 49.78125 48.890625 \n",
       "Q 56.203125 41.796875 56.203125 29.59375 \n",
       "z\n",
       "M 47.21875 32.234375 \n",
       "Q 47.125 39.59375 43.09375 43.984375 \n",
       "Q 39.0625 48.390625 32.421875 48.390625 \n",
       "Q 24.90625 48.390625 20.390625 44.140625 \n",
       "Q 15.875 39.890625 15.1875 32.171875 \n",
       "z\n",
       "\" id=\"DejaVuSans-101\"/>\n",
       "      <path d=\"M 54.890625 33.015625 \n",
       "L 54.890625 0 \n",
       "L 45.90625 0 \n",
       "L 45.90625 32.71875 \n",
       "Q 45.90625 40.484375 42.875 44.328125 \n",
       "Q 39.84375 48.1875 33.796875 48.1875 \n",
       "Q 26.515625 48.1875 22.3125 43.546875 \n",
       "Q 18.109375 38.921875 18.109375 30.90625 \n",
       "L 18.109375 0 \n",
       "L 9.078125 0 \n",
       "L 9.078125 54.6875 \n",
       "L 18.109375 54.6875 \n",
       "L 18.109375 46.1875 \n",
       "Q 21.34375 51.125 25.703125 53.5625 \n",
       "Q 30.078125 56 35.796875 56 \n",
       "Q 45.21875 56 50.046875 50.171875 \n",
       "Q 54.890625 44.34375 54.890625 33.015625 \n",
       "z\n",
       "\" id=\"DejaVuSans-110\"/>\n",
       "      <path d=\"M 9.8125 72.90625 \n",
       "L 51.703125 72.90625 \n",
       "L 51.703125 64.59375 \n",
       "L 19.671875 64.59375 \n",
       "L 19.671875 43.109375 \n",
       "L 48.578125 43.109375 \n",
       "L 48.578125 34.8125 \n",
       "L 19.671875 34.8125 \n",
       "L 19.671875 0 \n",
       "L 9.8125 0 \n",
       "z\n",
       "\" id=\"DejaVuSans-70\"/>\n",
       "      <path d=\"M -0.296875 72.90625 \n",
       "L 61.375 72.90625 \n",
       "L 61.375 64.59375 \n",
       "L 35.5 64.59375 \n",
       "L 35.5 0 \n",
       "L 25.59375 0 \n",
       "L 25.59375 64.59375 \n",
       "L -0.296875 64.59375 \n",
       "z\n",
       "\" id=\"DejaVuSans-84\"/>\n",
       "     </defs>\n",
       "     <g transform=\"translate(80.78125 38.476562)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-79\"/>\n",
       "      <use x=\"78.710938\" xlink:href=\"#DejaVuSans-112\"/>\n",
       "      <use x=\"142.1875\" xlink:href=\"#DejaVuSans-101\"/>\n",
       "      <use x=\"203.710938\" xlink:href=\"#DejaVuSans-110\"/>\n",
       "      <use x=\"267.089844\" xlink:href=\"#DejaVuSans-70\"/>\n",
       "      <use x=\"324.46875\" xlink:href=\"#DejaVuSans-65\"/>\n",
       "      <use x=\"392.876953\" xlink:href=\"#DejaVuSans-83\"/>\n",
       "      <use x=\"456.353516\" xlink:href=\"#DejaVuSans-84\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       "  <g id=\"axes_2\">\n",
       "   <g id=\"patch_8\">\n",
       "    <path d=\"M 409.017614 170.28 \n",
       "L 713.38125 170.28 \n",
       "L 713.38125 7.2 \n",
       "L 409.017614 7.2 \n",
       "z\n",
       "\" style=\"fill:#ffffff;\"/>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_3\">\n",
       "    <g id=\"xtick_7\">\n",
       "     <g id=\"line2d_29\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 422.852324 170.28 \n",
       "L 422.852324 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_30\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"422.852324\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_16\">\n",
       "      <!-- 0.0 -->\n",
       "      <g transform=\"translate(414.900762 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_8\">\n",
       "     <g id=\"line2d_31\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 478.191167 170.28 \n",
       "L 478.191167 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_32\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"478.191167\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_17\">\n",
       "      <!-- 0.2 -->\n",
       "      <g transform=\"translate(470.239605 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-50\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_9\">\n",
       "     <g id=\"line2d_33\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 533.53001 170.28 \n",
       "L 533.53001 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_34\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"533.53001\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_18\">\n",
       "      <!-- 0.4 -->\n",
       "      <g transform=\"translate(525.578448 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-52\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_10\">\n",
       "     <g id=\"line2d_35\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 588.868853 170.28 \n",
       "L 588.868853 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_36\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"588.868853\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_19\">\n",
       "      <!-- 0.6 -->\n",
       "      <g transform=\"translate(580.917291 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-54\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_11\">\n",
       "     <g id=\"line2d_37\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 644.207696 170.28 \n",
       "L 644.207696 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_38\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"644.207696\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_20\">\n",
       "      <!-- 0.8 -->\n",
       "      <g transform=\"translate(636.256134 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-56\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"xtick_12\">\n",
       "     <g id=\"line2d_39\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 699.546539 170.28 \n",
       "L 699.546539 7.2 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_40\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"699.546539\" xlink:href=\"#m6e3d9edeb4\" y=\"170.28\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_21\">\n",
       "      <!-- 1.0 -->\n",
       "      <g transform=\"translate(691.594977 184.878438)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-49\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_22\">\n",
       "     <!-- r/R [-] -->\n",
       "     <g transform=\"translate(546.690838 198.556563)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-114\"/>\n",
       "      <use x=\"41.113281\" xlink:href=\"#DejaVuSans-47\"/>\n",
       "      <use x=\"74.804688\" xlink:href=\"#DejaVuSans-82\"/>\n",
       "      <use x=\"144.287109\" xlink:href=\"#DejaVuSans-32\"/>\n",
       "      <use x=\"176.074219\" xlink:href=\"#DejaVuSans-91\"/>\n",
       "      <use x=\"215.087891\" xlink:href=\"#DejaVuSans-45\"/>\n",
       "      <use x=\"251.171875\" xlink:href=\"#DejaVuSans-93\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"matplotlib.axis_4\">\n",
       "    <g id=\"ytick_6\">\n",
       "     <g id=\"line2d_41\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 162.867273 \n",
       "L 713.38125 162.867273 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_42\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"162.867273\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_23\">\n",
       "      <!-- 0.00 -->\n",
       "      <g transform=\"translate(379.751989 166.666491)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_7\">\n",
       "     <g id=\"line2d_43\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 138.65226 \n",
       "L 713.38125 138.65226 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_44\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"138.65226\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_24\">\n",
       "      <!-- 0.01 -->\n",
       "      <g transform=\"translate(379.751989 142.451479)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-49\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_8\">\n",
       "     <g id=\"line2d_45\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 114.437248 \n",
       "L 713.38125 114.437248 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_46\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"114.437248\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_25\">\n",
       "      <!-- 0.02 -->\n",
       "      <g transform=\"translate(379.751989 118.236467)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-50\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_9\">\n",
       "     <g id=\"line2d_47\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 90.222236 \n",
       "L 713.38125 90.222236 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_48\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"90.222236\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_26\">\n",
       "      <!-- 0.03 -->\n",
       "      <g transform=\"translate(379.751989 94.021454)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-51\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_10\">\n",
       "     <g id=\"line2d_49\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 66.007223 \n",
       "L 713.38125 66.007223 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_50\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"66.007223\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_27\">\n",
       "      <!-- 0.04 -->\n",
       "      <g transform=\"translate(379.751989 69.806442)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-52\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_11\">\n",
       "     <g id=\"line2d_51\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 41.792211 \n",
       "L 713.38125 41.792211 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_52\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"41.792211\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_28\">\n",
       "      <!-- 0.05 -->\n",
       "      <defs>\n",
       "       <path d=\"M 10.796875 72.90625 \n",
       "L 49.515625 72.90625 \n",
       "L 49.515625 64.59375 \n",
       "L 19.828125 64.59375 \n",
       "L 19.828125 46.734375 \n",
       "Q 21.96875 47.46875 24.109375 47.828125 \n",
       "Q 26.265625 48.1875 28.421875 48.1875 \n",
       "Q 40.625 48.1875 47.75 41.5 \n",
       "Q 54.890625 34.8125 54.890625 23.390625 \n",
       "Q 54.890625 11.625 47.5625 5.09375 \n",
       "Q 40.234375 -1.421875 26.90625 -1.421875 \n",
       "Q 22.3125 -1.421875 17.546875 -0.640625 \n",
       "Q 12.796875 0.140625 7.71875 1.703125 \n",
       "L 7.71875 11.625 \n",
       "Q 12.109375 9.234375 16.796875 8.0625 \n",
       "Q 21.484375 6.890625 26.703125 6.890625 \n",
       "Q 35.15625 6.890625 40.078125 11.328125 \n",
       "Q 45.015625 15.765625 45.015625 23.390625 \n",
       "Q 45.015625 31 40.078125 35.4375 \n",
       "Q 35.15625 39.890625 26.703125 39.890625 \n",
       "Q 22.75 39.890625 18.8125 39.015625 \n",
       "Q 14.890625 38.140625 10.796875 36.28125 \n",
       "z\n",
       "\" id=\"DejaVuSans-53\"/>\n",
       "      </defs>\n",
       "      <g transform=\"translate(379.751989 45.59143)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-53\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"ytick_12\">\n",
       "     <g id=\"line2d_53\">\n",
       "      <path clip-path=\"url(#p8917949baa)\" d=\"M 409.017614 17.577199 \n",
       "L 713.38125 17.577199 \n",
       "\" style=\"fill:none;stroke:#b0b0b0;stroke-linecap:square;stroke-width:0.8;\"/>\n",
       "     </g>\n",
       "     <g id=\"line2d_54\">\n",
       "      <g>\n",
       "       <use style=\"stroke:#000000;stroke-width:0.8;\" x=\"409.017614\" xlink:href=\"#m60ac5f1668\" y=\"17.577199\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "     <g id=\"text_29\">\n",
       "      <!-- 0.06 -->\n",
       "      <g transform=\"translate(379.751989 21.376417)scale(0.1 -0.1)\">\n",
       "       <use xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"63.623047\" xlink:href=\"#DejaVuSans-46\"/>\n",
       "       <use x=\"95.410156\" xlink:href=\"#DejaVuSans-48\"/>\n",
       "       <use x=\"159.033203\" xlink:href=\"#DejaVuSans-54\"/>\n",
       "      </g>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"text_30\">\n",
       "     <!-- CT/d(r/R) [-] -->\n",
       "     <g transform=\"translate(373.672301 118.554063)rotate(-90)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-67\"/>\n",
       "      <use x=\"69.824219\" xlink:href=\"#DejaVuSans-84\"/>\n",
       "      <use x=\"130.908203\" xlink:href=\"#DejaVuSans-47\"/>\n",
       "      <use x=\"164.599609\" xlink:href=\"#DejaVuSans-100\"/>\n",
       "      <use x=\"228.076172\" xlink:href=\"#DejaVuSans-40\"/>\n",
       "      <use x=\"267.089844\" xlink:href=\"#DejaVuSans-114\"/>\n",
       "      <use x=\"308.203125\" xlink:href=\"#DejaVuSans-47\"/>\n",
       "      <use x=\"341.894531\" xlink:href=\"#DejaVuSans-82\"/>\n",
       "      <use x=\"411.376953\" xlink:href=\"#DejaVuSans-41\"/>\n",
       "      <use x=\"450.390625\" xlink:href=\"#DejaVuSans-32\"/>\n",
       "      <use x=\"482.177734\" xlink:href=\"#DejaVuSans-91\"/>\n",
       "      <use x=\"521.191406\" xlink:href=\"#DejaVuSans-45\"/>\n",
       "      <use x=\"557.275391\" xlink:href=\"#DejaVuSans-93\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "   <g id=\"line2d_55\">\n",
       "    <path clip-path=\"url(#p8917949baa)\" d=\"M 422.852324 162.867273 \n",
       "L 438.366693 162.867273 \n",
       "L 453.832247 116.564768 \n",
       "L 469.200327 51.646525 \n",
       "L 484.422579 36.028121 \n",
       "L 492.025878 30.676052 \n",
       "L 514.23863 29.475378 \n",
       "L 528.738616 28.128245 \n",
       "L 547.364721 27.272783 \n",
       "L 556.694542 26.129473 \n",
       "L 570.062522 25.488263 \n",
       "L 582.967324 24.655189 \n",
       "L 595.368346 23.863476 \n",
       "L 607.226568 23.459217 \n",
       "L 618.50468 23.385531 \n",
       "L 630.372986 23.643781 \n",
       "L 639.180573 24.599816 \n",
       "L 648.513299 26.230545 \n",
       "L 657.136012 29.062176 \n",
       "L 665.021582 33.411016 \n",
       "L 672.145198 39.692453 \n",
       "L 678.484446 47.388722 \n",
       "L 682.944886 57.442048 \n",
       "L 688.732588 66.136341 \n",
       "L 692.609238 81.04617 \n",
       "L 695.637132 99.399803 \n",
       "L 697.806744 119.163547 \n",
       "L 699.111248 141.56949 \n",
       "L 699.546539 14.612727 \n",
       "\" style=\"fill:none;stroke:#1f77b4;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "   </g>\n",
       "   <g id=\"line2d_56\">\n",
       "    <path clip-path=\"url(#p8917949baa)\" d=\"M 478.629329 45.985408 \n",
       "L 509.37823 37.052115 \n",
       "L 535.727095 34.416512 \n",
       "L 588.430238 34.860405 \n",
       "L 614.781792 27.25878 \n",
       "L 658.701064 32.918378 \n",
       "L 671.876841 41.518763 \n",
       "L 687.248608 68.443481 \n",
       "L 696.032451 105.563809 \n",
       "L 699.546539 162.867273 \n",
       "\" style=\"fill:none;stroke:#ff7f0e;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_9\">\n",
       "    <path d=\"M 409.017614 170.28 \n",
       "L 409.017614 7.2 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_10\">\n",
       "    <path d=\"M 713.38125 170.28 \n",
       "L 713.38125 7.2 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_11\">\n",
       "    <path d=\"M 409.017614 170.28 \n",
       "L 713.38125 170.28 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"patch_12\">\n",
       "    <path d=\"M 409.017614 7.2 \n",
       "L 713.38125 7.2 \n",
       "\" style=\"fill:none;stroke:#000000;stroke-linecap:square;stroke-linejoin:miter;stroke-width:0.8;\"/>\n",
       "   </g>\n",
       "   <g id=\"legend_2\">\n",
       "    <g id=\"patch_13\">\n",
       "     <path d=\"M 519.326776 165.28 \n",
       "L 603.072088 165.28 \n",
       "Q 605.072088 165.28 605.072088 163.28 \n",
       "L 605.072088 134.92375 \n",
       "Q 605.072088 132.92375 603.072088 132.92375 \n",
       "L 519.326776 132.92375 \n",
       "Q 517.326776 132.92375 517.326776 134.92375 \n",
       "L 517.326776 163.28 \n",
       "Q 517.326776 165.28 519.326776 165.28 \n",
       "z\n",
       "\" style=\"fill:#ffffff;opacity:0.8;stroke:#cccccc;stroke-linejoin:miter;\"/>\n",
       "    </g>\n",
       "    <g id=\"line2d_57\">\n",
       "     <path d=\"M 521.326776 141.022188 \n",
       "L 541.326776 141.022188 \n",
       "\" style=\"fill:none;stroke:#1f77b4;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "    </g>\n",
       "    <g id=\"line2d_58\"/>\n",
       "    <g id=\"text_31\">\n",
       "     <!-- SHARPy -->\n",
       "     <g transform=\"translate(549.326776 144.522188)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-83\"/>\n",
       "      <use x=\"63.476562\" xlink:href=\"#DejaVuSans-72\"/>\n",
       "      <use x=\"138.671875\" xlink:href=\"#DejaVuSans-65\"/>\n",
       "      <use x=\"207.080078\" xlink:href=\"#DejaVuSans-82\"/>\n",
       "      <use x=\"276.5625\" xlink:href=\"#DejaVuSans-80\"/>\n",
       "      <use x=\"336.865234\" xlink:href=\"#DejaVuSans-121\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "    <g id=\"line2d_59\">\n",
       "     <path d=\"M 521.326776 155.700312 \n",
       "L 541.326776 155.700312 \n",
       "\" style=\"fill:none;stroke:#ff7f0e;stroke-linecap:square;stroke-width:1.5;\"/>\n",
       "    </g>\n",
       "    <g id=\"line2d_60\"/>\n",
       "    <g id=\"text_32\">\n",
       "     <!-- OpenFAST -->\n",
       "     <g transform=\"translate(549.326776 159.200312)scale(0.1 -0.1)\">\n",
       "      <use xlink:href=\"#DejaVuSans-79\"/>\n",
       "      <use x=\"78.710938\" xlink:href=\"#DejaVuSans-112\"/>\n",
       "      <use x=\"142.1875\" xlink:href=\"#DejaVuSans-101\"/>\n",
       "      <use x=\"203.710938\" xlink:href=\"#DejaVuSans-110\"/>\n",
       "      <use x=\"267.089844\" xlink:href=\"#DejaVuSans-70\"/>\n",
       "      <use x=\"324.46875\" xlink:href=\"#DejaVuSans-65\"/>\n",
       "      <use x=\"392.876953\" xlink:href=\"#DejaVuSans-83\"/>\n",
       "      <use x=\"456.353516\" xlink:href=\"#DejaVuSans-84\"/>\n",
       "     </g>\n",
       "    </g>\n",
       "   </g>\n",
       "  </g>\n",
       " </g>\n",
       " <defs>\n",
       "  <clipPath id=\"p19e32ab1d1\">\n",
       "   <rect height=\"163.08\" width=\"304.363636\" x=\"43.78125\" y=\"7.2\"/>\n",
       "  </clipPath>\n",
       "  <clipPath id=\"p8917949baa\">\n",
       "   <rect height=\"163.08\" width=\"304.363636\" x=\"409.017614\" y=\"7.2\"/>\n",
       "  </clipPath>\n",
       " </defs>\n",
       "</svg>\n"
      ],
      "text/plain": [
       "<Figure size 864x216 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "fig, list_plots = plt.subplots(1, 2, figsize=(12, 3))\n",
    "\n",
    "list_plots[0].grid()\n",
    "list_plots[0].set_xlabel(\"r/R [-]\")\n",
    "list_plots[0].set_ylabel(\"CN/d(r/R) [-]\")\n",
    "list_plots[0].plot(r[0]/R, CN_drR[0], '-', label='SHARPy')\n",
    "list_plots[0].plot(of_rR, of_cNdrR, '-', label='OpenFAST')\n",
    "list_plots[0].legend()\n",
    "\n",
    "list_plots[1].grid()\n",
    "list_plots[1].set_xlabel(\"r/R [-]\")\n",
    "list_plots[1].set_ylabel(\"CT/d(r/R) [-]\")\n",
    "list_plots[1].plot(r[0]/R, CTan_drR[0], '-', label='SHARPy')\n",
    "list_plots[1].plot(of_rR, of_cTdrR, '-', label='OpenFAST')\n",
    "list_plots[1].legend()\n",
    "\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "Print the rotor thrust and power coefficients:"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "      OpenFAST SHARPy\n",
      "Cp[-] 0.49       0.52\n",
      "Ct[-] 0.70       0.70\n"
     ]
    }
   ],
   "source": [
    "print(\"      OpenFAST SHARPy\")\n",
    "print(\"Cp[-] %.2f       %.2f\" % (of_cp, Cp))\n",
    "print(\"Ct[-] %.2f       %.2f\" % (of_ct, Ct))"
   ]
  }
 ],
 "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": 4
}