Simulation NREL 5MW wind turbine

[1]:
%config InlineBackend.figure_format = 'svg'
from IPython.display import Image
url = 'https://raw.githubusercontent.com/ImperialCollegeLondon/sharpy/dev_doc/docs/source/content/example_notebooks/images/turbulence_no_legend.png'
Image(url=url, width=800)
[1]:

In this notebook 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.

OpenFAST: https://openfast.readthedocs.io

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

Load the required packages:

[2]:
# Required packages
%matplotlib inline
import numpy as np
import os
import matplotlib.pyplot as plt

# Required SHARPy modules
import sharpy.sharpy_main
import sharpy.utils.algebra as algebra
import sharpy.utils.generate_cases as gc
import sharpy.cases.templates.template_wt as template_wt

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

[3]:
of_rR = np.array([0.20158356, 0.3127131, 0.40794048, 0.5984148, 0.6936519, 0.85238045, 0.899999, 0.95555407, 0.98729974, 1.0])
of_cNdrR = np.array([0.08621394, 0.14687876, 0.19345148, 0.2942731, 0.36003628, 0.43748564, 0.44762507, 0.38839236, 0.29782477, 0.0])
of_cTdrR = np.array([0.048268348, 0.051957503, 0.05304592, 0.052862607, 0.056001827, 0.0536646, 0.050112925, 0.038993906, 0.023664437, 0.0])

of_ct = 0.69787693
of_cp = 0.48813498

Create SHARPy case

Definition of parameters

[4]:
# Mathematical constants
deg2rad = np.pi/180.

# Case
case = 'rotor'
route = './'

# Geometry discretization
chord_panels = np.array([8], dtype=int)
revs_in_wake = 5

# Operation
rotation_velocity = 12.1*2*np.pi/60
pitch_deg = 0. #degrees

# Wind
WSP = 12.
air_density = 1.225

# Simulation
dphi = 4.*deg2rad
revs_to_simulate = 5

Computation of associated parameters

[5]:
dt = dphi/rotation_velocity
time_steps = int(revs_to_simulate*2.*np.pi/dphi)
mstar = int(revs_in_wake*2.*np.pi/dphi)

Generation of the rotor geometry based on the information in the excel file

[6]:
op_params = {'rotation_velocity': rotation_velocity,
             'pitch_deg': pitch_deg,
             'wsp': WSP,
             'dt': dt}

geom_params = {'chord_panels':chord_panels,
               'tol_remove_points': 1e-8,
               'n_points_camber': 100,
               'm_distribution': 'uniform'}

excel_description = {'excel_file_name': route + '/source/type02_db_NREL5MW_v02.xlsx',
                     'excel_sheet_parameters': 'parameters',
                     'excel_sheet_structural_blade': 'structural_blade',
                     'excel_sheet_discretization_blade': 'discretization_blade',
                     'excel_sheet_aero_blade': 'aero_blade',
                     'excel_sheet_airfoil_info': 'airfoil_info',
                     'excel_sheet_airfoil_chord': 'airfoil_coord'}

options = {'camber_effect_on_twist': False,
           'user_defined_m_distribution_type': None,
           'include_polars': False}

rotor = template_wt.rotor_from_excel_type03(op_params,
                                            geom_params,
                                            excel_description,
                                            options)

WARNING: The poisson cofficient is assumed equal to 0.3
WARNING: Cross-section area is used as shear area
WARNING: Using perpendicular axis theorem to compute the inertia around xB
WARNING: Replacing node 29 by node 0
WARNING: Replacing node 58 by node 0

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.

[7]:
steady_simulation = False
[8]:
SimInfo = gc.SimulationInformation()
SimInfo.set_default_values()

if steady_simulation:
    SimInfo.solvers['SHARPy']['flow'] = ['BeamLoader',
                            'AerogridLoader',
                            'StaticCoupledRBM',
                            'BeamPlot',
                            'AerogridPlot',
                            'SaveData']
else:
    SimInfo.solvers['SHARPy']['flow'] = ['BeamLoader',
                            'AerogridLoader',
                            'StaticCoupledRBM',
                            'DynamicCoupled']

SimInfo.solvers['SHARPy']['case'] = case
SimInfo.solvers['SHARPy']['route'] = route
SimInfo.solvers['SHARPy']['write_log'] = True
SimInfo.set_variable_all_dicts('dt', dt)
SimInfo.set_variable_all_dicts('rho', air_density)

SimInfo.solvers['SteadyVelocityField']['u_inf'] = WSP
SimInfo.solvers['SteadyVelocityField']['u_inf_direction'] = np.array([0., 0., 1.])

SimInfo.solvers['BeamLoader']['unsteady'] = 'on'

SimInfo.solvers['AerogridLoader']['unsteady'] = 'on'
SimInfo.solvers['AerogridLoader']['mstar'] = mstar
SimInfo.solvers['AerogridLoader']['freestream_dir'] = np.array([0.,0.,0.])
SimInfo.solvers['AerogridLoader']['wake_shape_generator'] = 'HelicoidalWake'
SimInfo.solvers['AerogridLoader']['wake_shape_generator_input'] = {'u_inf': WSP,
                                                                   'u_inf_direction': SimInfo.solvers['SteadyVelocityField']['u_inf_direction'],
                                                                   'rotation_velocity': rotation_velocity*np.array([0., 0., 1.]),
                                                                   'dt': dt,
                                                                   'dphi1': dphi,
                                                                   'ndphi1': mstar,
                                                                   'r': 1.,
                                                                   'dphimax': 10*deg2rad}

SimInfo.solvers['StaticCoupledRBM']['structural_solver'] = 'RigidDynamicPrescribedStep'
SimInfo.solvers['StaticCoupledRBM']['structural_solver_settings'] = SimInfo.solvers['RigidDynamicPrescribedStep']
SimInfo.solvers['StaticCoupledRBM']['aero_solver'] = 'StaticUvlm'
SimInfo.solvers['StaticCoupledRBM']['aero_solver_settings'] = SimInfo.solvers['StaticUvlm']

SimInfo.solvers['StaticCoupledRBM']['tolerance'] = 1e-8
SimInfo.solvers['StaticCoupledRBM']['n_load_steps'] = 0
SimInfo.solvers['StaticCoupledRBM']['relaxation_factor'] = 0.

SimInfo.solvers['StaticUvlm']['num_cores'] = 8
SimInfo.solvers['StaticUvlm']['velocity_field_generator'] = 'SteadyVelocityField'
SimInfo.solvers['StaticUvlm']['velocity_field_input'] = SimInfo.solvers['SteadyVelocityField']

SimInfo.solvers['SaveData']['compress_float'] = True

# Only used for steady_simulation = False
SimInfo.solvers['StepUvlm']['convection_scheme'] = 3
SimInfo.solvers['StepUvlm']['num_cores'] = 8
SimInfo.solvers['StepUvlm']['velocity_field_generator'] = 'SteadyVelocityField'
SimInfo.solvers['StepUvlm']['velocity_field_input'] = SimInfo.solvers['SteadyVelocityField']

SimInfo.solvers['DynamicCoupled']['structural_solver'] = 'RigidDynamicPrescribedStep'
SimInfo.solvers['DynamicCoupled']['structural_solver_settings'] = SimInfo.solvers['RigidDynamicPrescribedStep']
SimInfo.solvers['DynamicCoupled']['aero_solver'] = 'StepUvlm'
SimInfo.solvers['DynamicCoupled']['aero_solver_settings'] = SimInfo.solvers['StepUvlm']
SimInfo.solvers['DynamicCoupled']['postprocessors'] = ['BeamPlot', 'AerogridPlot', 'Cleanup', 'SaveData']
SimInfo.solvers['DynamicCoupled']['postprocessors_settings'] = {'BeamPlot': SimInfo.solvers['BeamPlot'],
                                                             'AerogridPlot': SimInfo.solvers['AerogridPlot'],
                                                             'Cleanup': SimInfo.solvers['Cleanup'],
                                                             'SaveData': SimInfo.solvers['SaveData']}
SimInfo.solvers['DynamicCoupled']['minimum_steps'] = 0
SimInfo.solvers['DynamicCoupled']['include_unsteady_force_contribution'] = True
SimInfo.solvers['DynamicCoupled']['relaxation_factor'] = 0.
SimInfo.solvers['DynamicCoupled']['final_relaxation_factor'] = 0.
SimInfo.solvers['DynamicCoupled']['dynamic_relaxation'] = False
SimInfo.solvers['DynamicCoupled']['relaxation_steps'] = 0

# Define dynamic simulation (used regardless the value of "steady_simulation" variable)
SimInfo.define_num_steps(time_steps)
SimInfo.with_forced_vel = True
SimInfo.for_vel = np.zeros((time_steps,6), dtype=float)
SimInfo.for_vel[:,5] = rotation_velocity
SimInfo.for_acc = np.zeros((time_steps,6), dtype=float)
SimInfo.with_dynamic_forces = True
SimInfo.dynamic_forces = np.zeros((time_steps,rotor.StructuralInformation.num_node,6), dtype=float)

Generate simulation files

[9]:
gc.clean_test_files(SimInfo.solvers['SHARPy']['route'], SimInfo.solvers['SHARPy']['case'])
rotor.generate_h5_files(SimInfo.solvers['SHARPy']['route'], SimInfo.solvers['SHARPy']['case'])
SimInfo.generate_solver_file()
SimInfo.generate_dyn_file(time_steps)

Run SHARPy case

[10]:
sharpy_output = sharpy.sharpy_main.main(['', SimInfo.solvers['SHARPy']['route'] + SimInfo.solvers['SHARPy']['case'] + '.sharpy'])
--------------------------------------------------------------------------------
            ######  ##     ##    ###    ########  ########  ##    ##
           ##    ## ##     ##   ## ##   ##     ## ##     ##  ##  ##
           ##       ##     ##  ##   ##  ##     ## ##     ##   ####
            ######  ######### ##     ## ########  ########     ##
                 ## ##     ## ######### ##   ##   ##           ##
           ##    ## ##     ## ##     ## ##    ##  ##           ##
            ######  ##     ## ##     ## ##     ## ##           ##
--------------------------------------------------------------------------------
Aeroelastics Lab, Aeronautics Department.
    Copyright (c), Imperial College London.
    All rights reserved.
    License available at https://github.com/imperialcollegelondon/sharpy
Running SHARPy from /home/am3717/code/sharpy/docs/source/content/example_notebooks
SHARPy being run is in /home/am3717/code/sharpy
The branch being run is dev_setting_error
The version and commit hash are: v1.2.1-357-gbb7c4b4-bb7c4b4
SHARPy output folder set
     ./output//rotor/
Generating an instance of BeamLoader
Generating an instance of AerogridLoader
Variable shear_direction has no assigned value in the settings file.
    will default to the value: [1. 0. 0.]
Variable shear_exp has no assigned value in the settings file.
    will default to the value: 0.0
Variable h_ref has no assigned value in the settings file.
    will default to the value: 1.0
Variable h_corr has no assigned value in the settings file.
    will default to the value: 1.0
The aerodynamic grid contains 3 surfaces
  Surface 0, M=8, N=26
     Wake 0, M=450, N=26
  Surface 1, M=8, N=26
     Wake 1, M=450, N=26
  Surface 2, M=8, N=26
     Wake 2, M=450, N=26
  In total: 624 bound panels
  In total: 35100 wake panels
  Total number of panels = 35724
Generating an instance of StaticCoupledRBM
Generating an instance of RigidDynamicPrescribedStep
Generating an instance of StaticUvlm
i_step: 0, i_iter: 0
Resultant forces and moments: (array([2.54814836e-08, 5.82946313e-09, 2.09459948e+04]), array([1.53411975e-07, 3.24440634e-08, 5.96558361e+06]))
Generating an instance of DynamicCoupled
Generating an instance of RigidDynamicPrescribedStep
Generating an instance of StepUvlm
Generating an instance of BeamPlot
Generating an instance of AerogridPlot
Generating an instance of Cleanup
Generating an instance of SaveData



|=======|========|======|==============|==============|==============|==============|==============|
|  ts   |   t    | iter | struc ratio  |  iter time   | residual vel |  FoR_vel(x)  |  FoR_vel(z)  |
|=======|========|======|==============|==============|==============|==============|==============|
|   1   | 0.0551 |  1   |   0.000251   |  56.804234   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   2   | 0.1102 |  1   |   0.000177   |  60.232129   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   3   | 0.1653 |  1   |   0.000234   |  59.442611   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   4   | 0.2204 |  1   |   0.000290   |  64.383652   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   5   | 0.2755 |  1   |   0.000275   |  52.774436   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   6   | 0.3306 |  1   |   0.000293   |  48.828457   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   7   | 0.3857 |  1   |   0.000316   |  49.762675   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   8   | 0.4408 |  1   |   0.000227   |  49.627842   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|   9   | 0.4959 |  1   |   0.000193   |  58.254770   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  10   | 0.5510 |  1   |   0.000233   |  49.347327   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  11   | 0.6061 |  1   |   0.000222   |  51.950339   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  12   | 0.6612 |  1   |   0.000195   |  58.489715   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  13   | 0.7163 |  1   |   0.000190   |  55.497720   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  14   | 0.7713 |  1   |   0.000287   |  49.864559   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  15   | 0.8264 |  1   |   0.000198   |  56.915311   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  16   | 0.8815 |  1   |   0.000273   |  52.475129   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  17   | 0.9366 |  1   |   0.000212   |  53.472742   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  18   | 0.9917 |  1   |   0.000220   |  51.679432   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  19   | 1.0468 |  1   |   0.000230   |  49.468328   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  20   | 1.1019 |  1   |   0.000233   |  48.806374   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  21   | 1.1570 |  1   |   0.000214   |  51.507987   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  22   | 1.2121 |  1   |   0.000152   |  75.114590   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  23   | 1.2672 |  1   |   0.000178   |  63.764692   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  24   | 1.3223 |  1   |   0.000232   |  48.927093   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  25   | 1.3774 |  1   |   0.000205   |  55.065928   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  26   | 1.4325 |  1   |   0.000209   |  53.747631   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  27   | 1.4876 |  1   |   0.000215   |  65.487096   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  28   | 1.5427 |  1   |   0.000366   |  52.159906   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  29   | 1.5978 |  1   |   0.000286   |  50.405248   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  30   | 1.6529 |  1   |   0.000243   |  53.903698   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  31   | 1.7080 |  1   |   0.000189   |  55.762230   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  32   | 1.7631 |  1   |   0.000230   |  49.679477   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  33   | 1.8182 |  1   |   0.000267   |  49.140923   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  34   | 1.8733 |  1   |   0.000217   |  51.686471   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  35   | 1.9284 |  1   |   0.000227   |  49.924635   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  36   | 1.9835 |  1   |   0.000290   |  49.295199   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  37   | 2.0386 |  1   |   0.000231   |  49.096186   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  38   | 2.0937 |  1   |   0.000267   |  53.824351   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  39   | 2.1488 |  1   |   0.000274   |  50.208116   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  40   | 2.2039 |  1   |   0.000164   |  68.475332   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  41   | 2.2590 |  1   |   0.000219   |  51.476537   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  42   | 2.3140 |  1   |   0.000258   |  50.780835   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  43   | 2.3691 |  1   |   0.000205   |  54.508414   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  44   | 2.4242 |  1   |   0.000229   |  49.643120   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  45   | 2.4793 |  1   |   0.000423   |  49.510440   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  46   | 2.5344 |  1   |   0.000229   |  49.386177   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  47   | 2.5895 |  1   |   0.000224   |  51.096846   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  48   | 2.6446 |  1   |   0.000429   |  67.054611   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  49   | 2.6997 |  1   |   0.000674   |  49.034094   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  50   | 2.7548 |  1   |   0.000473   |  67.770922   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  51   | 2.8099 |  1   |   0.000623   |  52.057162   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  52   | 2.8650 |  1   |   0.000197   |  57.338530   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  53   | 2.9201 |  1   |   0.000211   |  53.323214   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  54   | 2.9752 |  1   |   0.000228   |  49.446313   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  55   | 3.0303 |  1   |   0.000226   |  50.531894   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  56   | 3.0854 |  1   |   0.000216   |  49.277498   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  57   | 3.1405 |  1   |   0.000229   |  49.378688   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  58   | 3.1956 |  1   |   0.000456   |  49.096131   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  59   | 3.2507 |  1   |   0.000638   |  49.407950   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  60   | 3.3058 |  1   |   0.000644   |  49.355487   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  61   | 3.3609 |  1   |   0.000429   |  52.381709   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  62   | 3.4160 |  1   |   0.000227   |  50.244867   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  63   | 3.4711 |  1   |   0.000226   |  50.090238   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  64   | 3.5262 |  1   |   0.000221   |  50.890983   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  65   | 3.5813 |  1   |   0.000292   |  52.582257   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  66   | 3.6364 |  1   |   0.000229   |  48.970770   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  67   | 3.6915 |  1   |   0.000211   |  53.921150   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  68   | 3.7466 |  1   |   0.000276   |  49.332785   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  69   | 3.8017 |  1   |   0.000405   |  49.861912   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  70   | 3.8567 |  1   |   0.000678   |  50.109032   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  71   | 3.9118 |  1   |   0.000627   |  50.052867   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  72   | 3.9669 |  1   |   0.000231   |  49.125818   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  73   | 4.0220 |  1   |   0.000215   |  52.475354   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  74   | 4.0771 |  1   |   0.000302   |  51.044905   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  75   | 4.1322 |  1   |   0.000232   |  49.154350   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  76   | 4.1873 |  1   |   0.000263   |  49.587372   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  77   | 4.2424 |  1   |   0.000382   |  48.959923   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  78   | 4.2975 |  1   |   0.000207   |  54.498279   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  79   | 4.3526 |  1   |   0.000209   |  55.866552   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  80   | 4.4077 |  1   |   0.000222   |  50.659905   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  81   | 4.4628 |  1   |   0.000232   |  49.229321   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  82   | 4.5179 |  1   |   0.000230   |  49.197802   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  83   | 4.5730 |  1   |   0.000244   |  58.532026   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  84   | 4.6281 |  1   |   0.000229   |  49.010943   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  85   | 4.6832 |  1   |   0.000236   |  49.181471   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  86   | 4.7383 |  1   |   0.000234   |  50.404876   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  87   | 4.7934 |  1   |   0.000223   |  50.863604   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  88   | 4.8485 |  1   |   0.000232   |  49.386036   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  89   | 4.9036 |  1   |   0.000292   |  48.966426   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  90   | 4.9587 |  1   |   0.000204   |  55.176201   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  91   | 5.0138 |  1   |   0.000381   |  49.363449   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  92   | 5.0689 |  1   |   0.000254   |  51.436284   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  93   | 5.1240 |  1   |   0.000427   |  54.203413   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  94   | 5.1791 |  1   |   0.000278   |  51.088371   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  95   | 5.2342 |  1   |   0.000221   |  51.267494   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  96   | 5.2893 |  1   |   0.000220   |  51.246126   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  97   | 5.3444 |  1   |   0.000261   |  50.044367   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  98   | 5.3994 |  1   |   0.000237   |  48.928531   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  99   | 5.4545 |  1   |   0.000230   |  48.869230   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  100  | 5.5096 |  1   |   0.000229   |  49.124024   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  101  | 5.5647 |  1   |   0.000230   |  49.170615   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  102  | 5.6198 |  1   |   0.000227   |  51.072559   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  103  | 5.6749 |  1   |   0.000230   |  49.106752   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  104  | 5.7300 |  1   |   0.000230   |  49.157374   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  105  | 5.7851 |  1   |   0.000230   |  48.729726   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  106  | 5.8402 |  1   |   0.000230   |  49.307981   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  107  | 5.8953 |  1   |   0.000232   |  49.303297   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  108  | 5.9504 |  1   |   0.000228   |  49.388923   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  109  | 6.0055 |  1   |   0.000251   |  49.505395   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  110  | 6.0606 |  1   |   0.000294   |  49.671154   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  111  | 6.1157 |  1   |   0.000637   |  49.106910   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  112  | 6.1708 |  1   |   0.000668   |  49.600747   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  113  | 6.2259 |  1   |   0.000694   |  48.788273   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  114  | 6.2810 |  1   |   0.000605   |  49.259812   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  115  | 6.3361 |  1   |   0.000235   |  49.096078   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  116  | 6.3912 |  1   |   0.000235   |  48.627974   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  117  | 6.4463 |  1   |   0.000231   |  48.687230   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  118  | 6.5014 |  1   |   0.000261   |  50.182136   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  119  | 6.5565 |  1   |   0.000220   |  51.480792   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  120  | 6.6116 |  1   |   0.000212   |  52.759361   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  121  | 6.6667 |  1   |   0.000234   |  61.395446   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  122  | 6.7218 |  1   |   0.000196   |  57.312678   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  123  | 6.7769 |  1   |   0.000196   |  57.288595   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  124  | 6.8320 |  1   |   0.000219   |  52.597533   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  125  | 6.8871 |  1   |   0.000279   |  51.215139   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  126  | 6.9421 |  1   |   0.000232   |  48.835388   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  127  | 6.9972 |  1   |   0.000228   |  49.442145   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  128  | 7.0523 |  1   |   0.000232   |  49.943434   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  129  | 7.1074 |  1   |   0.000265   |  48.943739   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  130  | 7.1625 |  1   |   0.000233   |  48.552538   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  131  | 7.2176 |  1   |   0.000235   |  48.576140   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  132  | 7.2727 |  1   |   0.000227   |  50.729522   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  133  | 7.3278 |  1   |   0.000199   |  57.042671   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  134  | 7.3829 |  1   |   0.000229   |  49.647710   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  135  | 7.4380 |  1   |   0.000377   |  53.191903   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  136  | 7.4931 |  1   |   0.000666   |  50.832518   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  137  | 7.5482 |  1   |   0.000560   |  52.129724   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  138  | 7.6033 |  1   |   0.000593   |  48.825715   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  139  | 7.6584 |  1   |   0.000231   |  49.220850   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  140  | 7.7135 |  1   |   0.000347   |  64.642874   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  141  | 7.7686 |  1   |   0.000615   |  54.791703   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  142  | 7.8237 |  1   |   0.000557   |  55.332506   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  143  | 7.8788 |  1   |   0.000649   |  51.958248   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  144  | 7.9339 |  1   |   0.000431   |  52.224418   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  145  | 7.9890 |  1   |   0.000532   |  60.207598   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  146  | 8.0441 |  1   |   0.000267   |  48.557877   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  147  | 8.0992 |  1   |   0.000283   |  50.720948   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  148  | 8.1543 |  1   |   0.000264   |  54.107664   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  149  | 8.2094 |  1   |   0.000225   |  49.998135   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  150  | 8.2645 |  1   |   0.000297   |  48.723537   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  151  | 8.3196 |  1   |   0.000229   |  48.991915   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  152  | 8.3747 |  1   |   0.000223   |  50.525580   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  153  | 8.4298 |  1   |   0.000372   |  48.880427   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  154  | 8.4848 |  1   |   0.000232   |  48.851891   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  155  | 8.5399 |  1   |   0.000231   |  48.731096   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  156  | 8.5950 |  1   |   0.000209   |  53.669843   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  157  | 8.6501 |  1   |   0.000203   |  55.525421   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  158  | 8.7052 |  1   |   0.000216   |  48.727216   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  159  | 8.7603 |  1   |   0.000240   |  52.981239   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  160  | 8.8154 |  1   |   0.000294   |  50.516618   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  161  | 8.8705 |  1   |   0.000198   |  56.957355   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  162  | 8.9256 |  1   |   0.000268   |  48.780843   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  163  | 8.9807 |  1   |   0.000228   |  49.216988   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  164  | 9.0358 |  1   |   0.000228   |  49.020854   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  165  | 9.0909 |  1   |   0.000236   |  48.516192   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  166  | 9.1460 |  1   |   0.000275   |  48.936249   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  167  | 9.2011 |  1   |   0.000229   |  48.984558   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  168  | 9.2562 |  1   |   0.000228   |  50.011495   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  169  | 9.3113 |  1   |   0.000230   |  49.079509   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  170  | 9.3664 |  1   |   0.000232   |  48.940832   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  171  | 9.4215 |  1   |   0.000266   |  48.539668   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  172  | 9.4766 |  1   |   0.000329   |  48.707876   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  173  | 9.5317 |  1   |   0.000344   |  49.214713   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  174  | 9.5868 |  1   |   0.000221   |  50.762513   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  175  | 9.6419 |  1   |   0.000336   |  55.302099   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  176  | 9.6970 |  1   |   0.000217   |  48.736059   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  177  | 9.7521 |  1   |   0.000231   |  48.750098   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  178  | 9.8072 |  1   |   0.000229   |  48.794680   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  179  | 9.8623 |  1   |   0.000225   |  50.519301   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  180  | 9.9174 |  1   |   0.000225   |  49.930789   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  181  | 9.9725 |  1   |   0.000265   |  50.077232   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  182  |10.0275 |  1   |   0.000285   |  50.327260   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  183  |10.0826 |  1   |   0.000319   |  50.215326   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  184  |10.1377 |  1   |   0.000288   |  48.350528   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  185  |10.1928 |  1   |   0.000233   |  48.372314   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  186  |10.2479 |  1   |   0.000232   |  49.104808   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  187  |10.3030 |  1   |   0.000292   |  48.972235   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  188  |10.3581 |  1   |   0.000655   |  48.972110   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  189  |10.4132 |  1   |   0.000233   |  48.457135   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  190  |10.4683 |  1   |   0.000229   |  49.120724   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  191  |10.5234 |  1   |   0.000222   |  51.106254   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  192  |10.5785 |  1   |   0.000232   |  48.794933   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  193  |10.6336 |  1   |   0.000292   |  49.537670   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  194  |10.6887 |  1   |   0.000230   |  49.236260   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  195  |10.7438 |  1   |   0.000603   |  52.529460   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  196  |10.7989 |  1   |   0.000697   |  48.383116   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  197  |10.8540 |  1   |   0.000347   |  48.288813   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  198  |10.9091 |  1   |   0.000227   |  49.550065   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  199  |10.9642 |  1   |   0.000216   |  48.675859   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  200  |11.0193 |  1   |   0.000231   |  49.152487   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  201  |11.0744 |  1   |   0.000235   |  48.778197   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  202  |11.1295 |  1   |   0.000373   |  57.791454   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  203  |11.1846 |  1   |   0.000582   |  57.753718   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  204  |11.2397 |  1   |   0.000451   |  48.675944   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  205  |11.2948 |  1   |   0.000235   |  48.274112   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  206  |11.3499 |  1   |   0.000232   |  48.917871   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  207  |11.4050 |  1   |   0.000293   |  49.482251   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  208  |11.4601 |  1   |   0.000297   |  48.304309   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  209  |11.5152 |  1   |   0.000263   |  49.162448   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  210  |11.5702 |  1   |   0.000227   |  50.845835   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  211  |11.6253 |  1   |   0.000232   |  48.809305   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  212  |11.6804 |  1   |   0.000293   |  48.950330   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  213  |11.7355 |  1   |   0.000298   |  48.481146   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  214  |11.7906 |  1   |   0.000234   |  48.222179   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  215  |11.8457 |  1   |   0.000225   |  48.228569   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  216  |11.9008 |  1   |   0.000218   |  48.173668   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  217  |11.9559 |  1   |   0.000230   |  48.496389   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  218  |12.0110 |  1   |   0.000212   |  53.240908   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  219  |12.0661 |  1   |   0.000224   |  51.372105   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  220  |12.1212 |  1   |   0.000419   |  49.490550   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  221  |12.1763 |  1   |   0.000459   |  49.025501   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  222  |12.2314 |  1   |   0.000210   |  53.157816   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  223  |12.2865 |  1   |   0.000233   |  48.099241   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  224  |12.3416 |  1   |   0.000241   |  59.450060   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  225  |12.3967 |  1   |   0.000224   |  50.153783   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  226  |12.4518 |  1   |   0.000169   |  67.550282   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  227  |12.5069 |  1   |   0.000635   |  51.062920   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  228  |12.5620 |  1   |   0.000233   |  48.498621   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  229  |12.6171 |  1   |   0.000647   |  48.837410   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  230  |12.6722 |  1   |   0.000685   |  48.524438   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  231  |12.7273 |  1   |   0.000645   |  47.972726   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  232  |12.7824 |  1   |   0.000698   |  48.332650   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  233  |12.8375 |  1   |   0.000291   |  49.344210   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  234  |12.8926 |  1   |   0.000231   |  48.621207   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  235  |12.9477 |  1   |   0.000211   |  53.099352   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  236  |13.0028 |  1   |   0.000231   |  48.779174   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  237  |13.0579 |  1   |   0.000221   |  52.508748   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  238  |13.1129 |  1   |   0.000197   |  57.071891   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  239  |13.1680 |  1   |   0.000232   |  48.635785   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  240  |13.2231 |  1   |   0.000228   |  49.653808   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  241  |13.2782 |  1   |   0.000232   |  49.880283   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  242  |13.3333 |  1   |   0.000219   |  50.989864   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  243  |13.3884 |  1   |   0.000218   |  51.562733   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  244  |13.4435 |  1   |   0.000175   |  60.525012   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  245  |13.4986 |  1   |   0.000208   |  53.714735   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  246  |13.5537 |  1   |   0.000188   |  60.430187   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  247  |13.6088 |  1   |   0.000580   |  55.369422   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  248  |13.6639 |  1   |   0.000376   |  52.945929   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  249  |13.7190 |  1   |   0.000221   |  50.621281   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  250  |13.7741 |  1   |   0.000222   |  51.930545   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  251  |13.8292 |  1   |   0.000203   |  55.750479   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  252  |13.8843 |  1   |   0.000228   |  49.736165   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  253  |13.9394 |  1   |   0.000239   |  48.733683   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  254  |13.9945 |  1   |   0.000228   |  48.860423   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  255  |14.0496 |  1   |   0.000234   |  48.158656   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  256  |14.1047 |  1   |   0.000222   |  48.995898   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  257  |14.1598 |  1   |   0.000181   |  61.699569   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  258  |14.2149 |  1   |   0.000217   |  51.853006   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  259  |14.2700 |  1   |   0.000231   |  48.370357   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  260  |14.3251 |  1   |   0.000231   |  48.536412   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  261  |14.3802 |  1   |   0.000231   |  48.706123   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  262  |14.4353 |  1   |   0.000321   |  48.889459   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  263  |14.4904 |  1   |   0.000232   |  48.413389   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  264  |14.5455 |  1   |   0.000183   |  61.637297   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  265  |14.6006 |  1   |   0.000266   |  51.996729   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  266  |14.6556 |  1   |   0.000440   |  51.124606   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  267  |14.7107 |  1   |   0.000206   |  54.482224   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  268  |14.7658 |  1   |   0.000227   |  49.396567   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  269  |14.8209 |  1   |   0.000216   |  48.399076   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  270  |14.8760 |  1   |   0.000225   |  49.932379   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  271  |14.9311 |  1   |   0.000229   |  48.837313   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  272  |14.9862 |  1   |   0.000234   |  48.303012   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  273  |15.0413 |  1   |   0.000297   |  48.557886   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  274  |15.0964 |  1   |   0.000232   |  48.462573   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  275  |15.1515 |  1   |   0.000235   |  48.074680   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  276  |15.2066 |  1   |   0.000237   |  47.978437   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  277  |15.2617 |  1   |   0.000218   |  48.799704   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  278  |15.3168 |  1   |   0.000180   |  61.907709   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  279  |15.3719 |  1   |   0.000223   |  50.256748   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  280  |15.4270 |  1   |   0.000229   |  49.243923   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  281  |15.4821 |  1   |   0.000453   |  49.924487   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  282  |15.5372 |  1   |   0.000350   |  49.501563   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  283  |15.5923 |  1   |   0.000234   |  48.333380   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  284  |15.6474 |  1   |   0.000223   |  50.163273   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  285  |15.7025 |  1   |   0.000225   |  50.068699   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  286  |15.7576 |  1   |   0.000229   |  48.458516   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  287  |15.8127 |  1   |   0.000231   |  49.273332   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  288  |15.8678 |  1   |   0.000232   |  48.548550   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  289  |15.9229 |  1   |   0.000224   |  51.276163   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  290  |15.9780 |  1   |   0.000233   |  48.570122   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  291  |16.0331 |  1   |   0.000192   |  50.266936   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  292  |16.0882 |  1   |   0.000295   |  48.692444   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  293  |16.1433 |  1   |   0.000202   |  48.893389   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  294  |16.1983 |  1   |   0.000443   |  50.813372   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  295  |16.2534 |  1   |   0.000529   |  60.601432   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  296  |16.3085 |  1   |   0.000571   |  50.993108   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  297  |16.3636 |  1   |   0.000228   |  49.065097   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  298  |16.4187 |  1   |   0.000234   |  48.226984   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  299  |16.4738 |  1   |   0.000233   |  48.534926   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  300  |16.5289 |  1   |   0.000627   |  48.788802   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  301  |16.5840 |  1   |   0.000555   |  60.530895   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  302  |16.6391 |  1   |   0.000628   |  48.944436   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  303  |16.6942 |  1   |   0.000427   |  49.049230   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  304  |16.7493 |  1   |   0.000236   |  47.892993   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  305  |16.8044 |  1   |   0.000325   |  48.441506   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  306  |16.8595 |  1   |   0.000236   |  48.047404   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  307  |16.9146 |  1   |   0.000232   |  48.215284   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  308  |16.9697 |  1   |   0.000226   |  50.124858   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  309  |17.0248 |  1   |   0.000213   |  52.395805   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  310  |17.0799 |  1   |   0.000222   |  51.302776   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  311  |17.1350 |  1   |   0.000439   |  51.171773   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  312  |17.1901 |  1   |   0.000232   |  49.032143   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  313  |17.2452 |  1   |   0.000233   |  48.628400   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  314  |17.3003 |  1   |   0.000237   |  47.940311   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  315  |17.3554 |  1   |   0.000216   |  51.865635   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  316  |17.4105 |  1   |   0.000373   |  54.905208   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  317  |17.4656 |  1   |   0.000310   |  48.191338   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  318  |17.5207 |  1   |   0.000409   |  48.607705   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  319  |17.5758 |  1   |   0.000350   |  53.126807   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  320  |17.6309 |  1   |   0.000236   |  47.827493   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  321  |17.6860 |  1   |   0.000233   |  48.142512   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  322  |17.7410 |  1   |   0.000233   |  48.315800   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  323  |17.7961 |  1   |   0.000203   |  55.821999   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  324  |17.8512 |  1   |   0.000227   |  49.786573   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  325  |17.9063 |  1   |   0.000424   |  53.387859   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  326  |17.9614 |  1   |   0.000447   |  50.685506   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  327  |18.0165 |  1   |   0.000235   |  48.078798   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  328  |18.0716 |  1   |   0.000232   |  48.322922   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  329  |18.1267 |  1   |   0.000223   |  48.098878   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  330  |18.1818 |  1   |   0.000201   |  48.656743   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  331  |18.2369 |  1   |   0.000273   |  48.314217   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  332  |18.2920 |  1   |   0.000224   |  48.475068   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  333  |18.3471 |  1   |   0.000230   |  48.642246   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  334  |18.4022 |  1   |   0.000194   |  57.654265   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  335  |18.4573 |  1   |   0.000225   |  50.092508   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  336  |18.5124 |  1   |   0.000291   |  49.505815   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  337  |18.5675 |  1   |   0.000228   |  49.799143   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  338  |18.6226 |  1   |   0.000272   |  48.344175   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  339  |18.6777 |  1   |   0.000219   |  48.669889   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  340  |18.7328 |  1   |   0.000221   |  48.424671   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  341  |18.7879 |  1   |   0.000204   |  55.293948   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  342  |18.8430 |  1   |   0.000361   |  48.517776   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  343  |18.8981 |  1   |   0.000218   |  52.256841   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  344  |18.9532 |  1   |   0.000232   |  48.441100   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  345  |19.0083 |  1   |   0.000234   |  48.389977   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  346  |19.0634 |  1   |   0.000230   |  48.982534   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  347  |19.1185 |  1   |   0.000231   |  48.746604   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  348  |19.1736 |  1   |   0.000229   |  56.961390   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  349  |19.2287 |  1   |   0.000188   |  59.830517   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  350  |19.2837 |  1   |   0.000296   |  51.906979   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  351  |19.3388 |  1   |   0.000217   |  51.754677   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  352  |19.3939 |  1   |   0.000244   |  48.090983   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  353  |19.4490 |  1   |   0.000234   |  48.219884   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  354  |19.5041 |  1   |   0.000211   |  53.274204   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  355  |19.5592 |  1   |   0.000236   |  47.988238   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  356  |19.6143 |  1   |   0.000236   |  48.149312   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  357  |19.6694 |  1   |   0.000231   |  48.684465   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  358  |19.7245 |  1   |   0.000245   |  48.573752   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  359  |19.7796 |  1   |   0.000255   |  48.387298   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  360  |19.8347 |  1   |   0.000233   |  48.496933   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  361  |19.8898 |  1   |   0.000298   |  48.227055   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  362  |19.9449 |  1   |   0.000277   |  50.033512   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  363  |20.0000 |  1   |   0.000230   |  48.865905   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  364  |20.0551 |  1   |   0.000233   |  48.282476   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  365  |20.1102 |  1   |   0.000231   |  48.869913   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  366  |20.1653 |  1   |   0.000240   |  48.197653   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  367  |20.2204 |  1   |   0.000236   |  47.995288   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  368  |20.2755 |  1   |   0.000218   |  48.416211   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  369  |20.3306 |  1   |   0.000229   |  49.650859   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  370  |20.3857 |  1   |   0.000281   |  50.925766   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  371  |20.4408 |  1   |   0.000293   |  49.072721   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  372  |20.4959 |  1   |   0.000237   |  47.848434   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  373  |20.5510 |  1   |   0.000461   |  48.577324   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  374  |20.6061 |  1   |   0.000288   |  48.149040   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  375  |20.6612 |  1   |   0.000232   |  48.140545   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  376  |20.7163 |  1   |   0.000627   |  48.528365   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  377  |20.7713 |  1   |   0.000649   |  48.815492   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  378  |20.8264 |  1   |   0.000220   |  48.691541   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  379  |20.8815 |  1   |   0.000234   |  48.163874   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  380  |20.9366 |  1   |   0.000397   |  47.970319   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  381  |20.9917 |  1   |   0.000266   |  49.183593   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  382  |21.0468 |  1   |   0.000228   |  49.448159   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  383  |21.1019 |  1   |   0.000230   |  48.558385   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  384  |21.1570 |  1   |   0.000279   |  48.382636   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  385  |21.2121 |  1   |   0.000237   |  47.482272   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  386  |21.2672 |  1   |   0.000233   |  48.183362   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  387  |21.3223 |  1   |   0.000232   |  48.062299   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  388  |21.3774 |  1   |   0.000235   |  48.362948   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  389  |21.4325 |  1   |   0.000233   |  48.422455   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  390  |21.4876 |  1   |   0.000236   |  48.077169   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  391  |21.5427 |  1   |   0.000218   |  48.317787   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  392  |21.5978 |  1   |   0.000231   |  48.382455   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  393  |21.6529 |  1   |   0.000233   |  48.310206   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  394  |21.7080 |  1   |   0.000219   |  49.225511   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  395  |21.7631 |  1   |   0.000238   |  47.725182   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  396  |21.8182 |  1   |   0.000403   |  48.100903   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  397  |21.8733 |  1   |   0.000227   |  49.886770   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  398  |21.9284 |  1   |   0.000226   |  49.658201   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  399  |21.9835 |  1   |   0.000236   |  48.305113   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  400  |22.0386 |  1   |   0.000213   |  53.014337   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  401  |22.0937 |  1   |   0.000323   |  48.752367   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  402  |22.1488 |  1   |   0.000231   |  48.655110   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  403  |22.2039 |  1   |   0.000302   |  49.894845   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  404  |22.2590 |  1   |   0.000192   |  54.227928   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  405  |22.3140 |  1   |   0.000213   |  53.815790   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  406  |22.3691 |  1   |   0.000231   |  49.235416   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  407  |22.4242 |  1   |   0.000232   |  48.387503   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  408  |22.4793 |  1   |   0.000220   |  52.815450   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  409  |22.5344 |  1   |   0.000232   |  48.253041   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  410  |22.5895 |  1   |   0.000218   |  52.584628   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  411  |22.6446 |  1   |   0.000222   |  50.874607   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  412  |22.6997 |  1   |   0.000234   |  48.125687   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  413  |22.7548 |  1   |   0.000314   |  65.912785   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  414  |22.8099 |  1   |   0.000294   |  48.712801   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  415  |22.8650 |  1   |   0.000234   |  48.264049   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  416  |22.9201 |  1   |   0.000219   |  52.101943   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  417  |22.9752 |  1   |   0.000233   |  48.181907   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  418  |23.0303 |  1   |   0.000239   |  47.737343   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  419  |23.0854 |  1   |   0.000189   |  55.077763   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  420  |23.1405 |  1   |   0.000252   |  48.412007   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  421  |23.1956 |  1   |   0.000207   |  52.025445   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  422  |23.2507 |  1   |   0.000202   |  55.306192   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  423  |23.3058 |  1   |   0.000235   |  48.069254   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  424  |23.3609 |  1   |   0.000268   |  48.786326   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  425  |23.4160 |  1   |   0.000234   |  48.662434   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  426  |23.4711 |  1   |   0.000430   |  48.229208   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  427  |23.5262 |  1   |   0.000660   |  48.355178   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  428  |23.5813 |  1   |   0.000475   |  47.940589   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  429  |23.6364 |  1   |   0.000464   |  48.627356   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  430  |23.6915 |  1   |   0.000231   |  48.519762   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  431  |23.7466 |  1   |   0.000233   |  48.359547   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  432  |23.8017 |  1   |   0.000273   |  48.313746   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  433  |23.8567 |  1   |   0.000237   |  47.977213   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  434  |23.9118 |  1   |   0.000226   |  49.786822   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  435  |23.9669 |  1   |   0.000231   |  48.641675   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  436  |24.0220 |  1   |   0.000263   |  50.056180   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  437  |24.0771 |  1   |   0.000239   |  48.169318   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  438  |24.1322 |  1   |   0.000237   |  47.903935   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  439  |24.1873 |  1   |   0.000235   |  48.071315   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  440  |24.2424 |  1   |   0.000298   |  47.997039   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  441  |24.2975 |  1   |   0.000204   |  51.404477   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  442  |24.3526 |  1   |   0.000234   |  48.367934   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  443  |24.4077 |  1   |   0.000297   |  48.626223   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  444  |24.4628 |  1   |   0.000254   |  48.214660   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  445  |24.5179 |  1   |   0.000234   |  47.731147   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  446  |24.5730 |  1   |   0.000268   |  48.701918   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  447  |24.6281 |  1   |   0.000233   |  48.425494   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  448  |24.6832 |  1   |   0.000233   |  48.324203   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  449  |24.7383 |  1   |   0.000228   |  49.378323   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
|  450  |24.7934 |  1   |   0.000229   |  48.777402   |   0.000000   | 0.000000e+00 | 0.000000e+00 |
...Finished
FINISHED - Elapsed time = 23245.3692466 seconds
FINISHED - CPU process time = 92349.9106726 seconds

Postprocessing

Read the structural and aerodynamic information of the last time step

[11]:
tstep = sharpy_output.structure.timestep_info[-1]
astep = sharpy_output.aero.timestep_info[-1]

Separate the structure into blades

[12]:
# Define beams
ielem = 0
nblades = np.max(sharpy_output.structure.beam_number) + 1
nodes_blade = []
first_node = 0
for iblade in range(nblades):
    nodes_blade.append(np.zeros((sharpy_output.structure.num_node,), dtype=bool))
    while sharpy_output.structure.beam_number[ielem] <= iblade:
        ielem += 1
        if ielem == sharpy_output.structure.num_elem:
            break
    nodes_blade[iblade][first_node:sharpy_output.structure.connectivities[ielem-1,1]+1] = True
    first_node = sharpy_output.structure.connectivities[ielem-1,1]+1

Compute the radial position of the nodes and initialise the rest of the variables

[13]:
r = []
c = []
dr = []
forces = []
CN_drR = []
CTan_drR = []
CP_drR = []
nodes_num = []
for iblade in range(nblades):
    forces.append(tstep.steady_applied_forces[nodes_blade[iblade]].copy())

    nodes_num.append(np.arange(0, sharpy_output.structure.num_node, 1)[nodes_blade[iblade]])

    r.append(np.linalg.norm(tstep.pos[nodes_blade[iblade], :], axis=1))
    dr.append(np.zeros(np.sum(nodes_blade[iblade])))
    dr[iblade][0] = 0.5*(r[iblade][1]-r[iblade][0])
    dr[iblade][-1] = 0.5 * (r[iblade][-1] - r[iblade][-2])
    for inode in range(1,len(r[iblade]) - 1):
        dr[iblade][inode] = 0.5*(r[iblade][inode+1] - r[iblade][inode-1])

    CN_drR.append(np.zeros(len(r[iblade])))
    c.append(np.zeros(len(r[iblade])))
    CTan_drR.append(np.zeros(len(r[iblade])))
    CP_drR.append(np.zeros(len(r[iblade])))

Transform the loads computed by SHARPy into out-of-plane and in-plane components

[14]:
rho = sharpy_output.settings['StaticCoupledRBM']['aero_solver_settings']['rho']
uinf = sharpy_output.settings['StaticCoupledRBM']['aero_solver_settings']['velocity_field_input']['u_inf']
R = np.max(r[0])
Cp = 0
Ct = 0

global_force_factor = 0.5 * rho * uinf** 2 * np.pi * R**2
global_power_factor = global_force_factor*uinf
for iblade in range(nblades):
    for inode in range(len(r[iblade])):
        forces[iblade][inode, 0] *= 0. # Discard the spanwise component

        node_global_index = nodes_num[iblade][inode]
        ielem = sharpy_output.structure.node_master_elem[node_global_index, 0]
        inode_in_elem = sharpy_output.structure.node_master_elem[node_global_index, 1]
        CAB = algebra.crv2rotation(tstep.psi[ielem, inode_in_elem, :])

        c[iblade][inode] = sharpy_output.aero.aero_dict['chord'][ielem,inode_in_elem]

        forces_AFoR = np.dot(CAB, forces[iblade][inode, 0:3])

        CN_drR[iblade][inode] = forces_AFoR[2]/dr[iblade][inode]*R / global_force_factor
        CTan_drR[iblade][inode] = np.linalg.norm(forces_AFoR[0:2])/dr[iblade][inode]*R  / global_force_factor
        CP_drR[iblade][inode] = np.linalg.norm(forces_AFoR[0:2])/dr[iblade][inode]*R  * r[iblade][inode]*rotation_velocity / global_power_factor

    Cp += np.sum(CP_drR[iblade]*dr[iblade]/R)
    Ct += np.sum(CN_drR[iblade]*dr[iblade]/R)

Results

Plot of the loads along the blade

[15]:
fig, list_plots = plt.subplots(1, 2, figsize=(12, 3))

list_plots[0].grid()
list_plots[0].set_xlabel("r/R [-]")
list_plots[0].set_ylabel("CN/d(r/R) [-]")
list_plots[0].plot(r[0]/R, CN_drR[0], '-', label='SHARPy')
list_plots[0].plot(of_rR, of_cNdrR, '-', label='OpenFAST')
list_plots[0].legend()

list_plots[1].grid()
list_plots[1].set_xlabel("r/R [-]")
list_plots[1].set_ylabel("CT/d(r/R) [-]")
list_plots[1].plot(r[0]/R, CTan_drR[0], '-', label='SHARPy')
list_plots[1].plot(of_rR, of_cTdrR, '-', label='OpenFAST')
list_plots[1].legend()

plt.show()
../../_images/content_example_notebooks_wind_turbine_34_0.svg

Print the rotor thrust and power coefficients

[16]:
print("      OpenFAST SHARPy")
print("Cp[-] %.2f       %.2f" % (of_cp, Cp))
print("Ct[-] %.2f       %.2f" % (of_ct, Ct))
      OpenFAST SHARPy
Cp[-] 0.49       0.52
Ct[-] 0.70       0.70