In [1]:
from pydelft import grd, dep
import matplotlib.pyplot as plt
import numpy as np
import sys 
import os
%matplotlib inline

In [2]:
gfname = '../test_data/natori_10_v1.0.grd'
dfname = '../test_data/natori_10_v1.0.dep'

dep.py


In [3]:
grd_fname = None

f = open(dfname, 'r')
lines = f.readlines()

if '\t' in lines[0]:
    delimiter = '\t'
else:
    delimiter = 'space'
print('Delimiter = %s' % delimiter)


z = []

if delimiter != 'space':
    for line in lines:
        for s in line.split('%s' % delimiter):
            z.append(float(s))
else:
    for line in lines:
        for s in line.split():
            z.append(float(s))

# Reshape depth
if grd_fname:
    grid = grd.grd()
    grid.read_grd(grd_fname)
    z = np.ma.masked_equal(z, -999.)
    z = np.ma.compressed(z)
    z = np.reshape(z, (grid.m, grid.n))
else:
    if delimiter != 'space':
        z = np.reshape(z, (np.size(lines),
            np.size(line.split('%s' % delimiter))))
        z = z[:-1, :-1] # get rid of -999. Nans
        z = np.reshape(z, (np.size(lines)-1, np.size(line.split())-1))

        z = z[:-1, :-1] # get rid of -999. Nans
    else:
        a = np.where(np.array(z) == -999.)
        cols = np.max(np.ediff1d(a))
        rows = np.size(np.array_split(a[0],np.where(np.diff(a[0])!=1)[0]+1))+1
        z = np.reshape(z, (rows, cols))


Delimiter = space

Test functionality


In [4]:
grid = grd.grd()
depth = dep.dep()

grid.read_grd(gfname)
depth.read_dep(dfname)

In [6]:
depth.depth.shape


Out[6]:
(326, 975)

In [7]:
grid.x.shape


Out[7]:
(325, 974)

In [9]:
depth.depth[0]


Out[9]:
array([ -2.82500000e+00,  -2.95000000e+00,  -3.00000000e+00,
        -3.07500000e+00,  -2.97500000e+00,  -2.87500000e+00,
        -2.85000000e+00,  -3.00000000e+00,  -2.90000000e+00,
        -2.90000000e+00,  -2.85000000e+00,  -2.72500000e+00,
        -2.92500000e+00,  -2.97500000e+00,  -3.02500000e+00,
        -2.95000000e+00,  -2.80000000e+00,  -2.77500000e+00,
        -2.75000000e+00,  -2.62500000e+00,  -2.42500000e+00,
        -2.47500000e+00,  -2.52500000e+00,  -2.55000000e+00,
        -2.50000000e+00,  -2.22500000e+00,  -2.00000000e+00,
        -2.02500000e+00,  -1.92500000e+00,  -1.90000000e+00,
        -2.50000000e+00,  -2.55000000e+00,  -2.62500000e+00,
        -3.05000000e+00,  -2.45000000e+00,  -1.22500000e+00,
        -1.12500000e+00,  -1.12500000e+00,  -1.10000000e+00,
        -1.10000000e+00,  -1.12500000e+00,  -1.00000000e+00,
        -1.00000000e+00,  -1.00000000e+00,  -1.00000000e+00,
        -1.00000000e+00,  -1.00000000e+00,  -1.10000000e+00,
        -7.75000000e-01,  -8.00000000e-01,  -8.00000000e-01,
        -8.00000000e-01,  -7.50000000e-01,  -8.00000000e-01,
        -7.50000000e-01,  -6.50000000e-01,  -6.75000000e-01,
        -7.00000000e-01,  -7.00000000e-01,  -7.00000000e-01,
        -7.00000000e-01,  -7.00000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.50000000e-01,  -7.50000000e-01,
        -7.25000000e-01,  -7.75000000e-01,  -7.00000000e-01,
        -6.00000000e-01,  -7.00000000e-01,  -7.75000000e-01,
        -7.25000000e-01,  -7.75000000e-01,  -6.75000000e-01,
        -5.25000000e-01,  -1.25000000e+00,  -1.45000000e+00,
        -6.87500000e+00,  -7.92500000e+00,  -8.07500000e+00,
        -7.52500000e+00,  -5.55000000e+00,  -1.65000000e+00,
        -6.00000000e-01,  -6.00000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.75000000e-01,  -6.00000000e-01,
        -6.50000000e-01,  -6.75000000e-01,  -5.50000000e-01,
        -5.50000000e-01,  -6.00000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -6.75000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -6.50000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -5.75000000e-01,
        -5.50000000e-01,  -5.50000000e-01,  -5.75000000e-01,
        -6.25000000e-01,  -5.25000000e-01,  -5.50000000e-01,
        -5.75000000e-01,  -5.25000000e-01,  -5.50000000e-01,
        -5.00000000e-01,  -5.50000000e-01,  -5.00000000e-01,
        -5.25000000e-01,  -5.00000000e-01,  -5.00000000e-01,
        -5.00000000e-01,  -5.00000000e-01,  -5.25000000e-01,
        -5.00000000e-01,  -5.00000000e-01,  -5.25000000e-01,
        -6.00000000e-01,  -5.75000000e-01,  -6.75000000e-01,
        -1.10000000e+00,  -5.75000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.25000000e-01,  -5.50000000e-01,
        -6.25000000e-01,  -6.00000000e-01,  -5.75000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -6.25000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -5.75000000e-01,
        -5.50000000e-01,  -5.75000000e-01,  -5.50000000e-01,
        -7.00000000e-01,  -6.50000000e-01,  -6.75000000e-01,
        -7.00000000e-01,  -6.75000000e-01,  -6.50000000e-01,
        -6.75000000e-01,  -7.25000000e-01,  -6.00000000e-01,
        -7.00000000e-01,  -6.50000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -1.15000000e+00,
        -1.20000000e+00,  -1.20000000e+00,  -1.22500000e+00,
        -1.17500000e+00,  -1.07500000e+00,  -1.42500000e+00,
        -1.20000000e+00,  -1.27500000e+00,  -1.27500000e+00,
        -1.15000000e+00,  -1.20000000e+00,  -1.05000000e+00,
        -1.02500000e+00,  -1.00000000e+00,  -5.00000000e-01,
        -5.25000000e-01,  -4.50000000e-01,  -5.50000000e-01,
        -6.25000000e-01,  -6.00000000e-01,  -5.75000000e-01,
        -5.50000000e-01,  -5.00000000e-01,  -5.25000000e-01,
        -5.25000000e-01,  -5.50000000e-01,  -5.25000000e-01,
        -5.50000000e-01,  -5.25000000e-01,  -5.75000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -6.25000000e-01,
        -6.00000000e-01,  -6.25000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -7.50000000e-01,  -8.50000000e-01,
        -7.25000000e-01,  -6.00000000e-01,  -5.00000000e-01,
        -5.25000000e-01,  -5.00000000e-01,  -5.50000000e-01,
        -5.75000000e-01,  -6.00000000e-01,  -6.00000000e-01,
        -5.75000000e-01,  -6.00000000e-01,  -5.00000000e-01,
        -5.00000000e-01,  -5.50000000e-01,  -6.00000000e-01,
        -6.00000000e-01,  -6.00000000e-01,  -5.75000000e-01,
        -5.25000000e-01,  -5.75000000e-01,   1.00000000e+00,
        -7.50000000e-01,  -5.00000000e-01,  -5.00000000e-01,
        -5.00000000e-01,  -4.50000000e-01,  -4.00000000e-01,
        -4.00000000e-01,  -5.00000000e-01,  -5.00000000e-01,
        -4.75000000e-01,  -5.25000000e-01,  -7.00000000e-01,
        -6.25000000e-01,  -6.00000000e-01,  -4.50000000e-01,
        -4.00000000e-01,  -1.85000000e+00,  -1.92500000e+00,
        -1.87500000e+00,  -4.75000000e-01,  -2.00000000e-01,
        -2.25000000e-01,  -2.25000000e-01,  -2.00000000e-01,
        -3.00000000e-01,  -2.75000000e-01,  -3.00000000e-01,
        -3.00000000e-01,  -3.00000000e-01,  -2.75000000e-01,
        -3.00000000e-01,  -3.00000000e-01,  -3.75000000e-01,
        -3.75000000e-01,  -3.75000000e-01,  -4.75000000e-01,
        -5.00000000e-01,  -1.17500000e+00,  -2.10000000e+00,
        -2.25000000e+00,  -2.50000000e+00,  -2.50000000e+00,
        -2.45000000e+00,  -2.31666670e+00,  -2.13333330e+00,
        -2.03333330e+00,  -2.05000000e+00,  -1.91666670e+00,
        -1.91666670e+00,  -2.20000000e+00,  -2.26666670e+00,
        -2.18333330e+00,  -2.23333330e+00,  -2.18333330e+00,
        -2.10000000e+00,  -2.10000000e+00,  -2.05000000e+00,
        -2.00000000e+00,  -2.01666670e+00,  -1.93333330e+00,
        -1.80000000e+00,  -2.15000000e+00,  -2.30000000e+00,
        -2.35000000e+00,  -2.36666670e+00,  -2.36666670e+00,
        -2.30000000e+00,  -2.26666670e+00,  -2.30000000e+00,
        -1.96666670e+00,  -9.83333330e-01,  -1.06666670e+00,
        -1.25000000e+00,  -5.33333330e-01,   1.00000000e+00,
        -1.23333330e+00,  -1.63333330e+00,  -1.61666670e+00,
        -1.73333330e+00,  -1.71666670e+00,  -1.71666670e+00,
        -1.73333330e+00,  -1.71666670e+00,  -1.75000000e+00,
        -1.76666670e+00,  -1.53333330e+00,  -1.55000000e+00,
        -1.66666670e+00,  -1.63333330e+00,  -1.66666670e+00,
        -1.71666670e+00,  -1.85000000e+00,  -1.81666670e+00,
        -1.95000000e+00,  -1.80000000e+00,  -1.71666670e+00,
        -1.58333330e+00,  -1.56666670e+00,  -1.43333330e+00,
        -1.50000000e+00,  -1.50000000e+00,  -1.60000000e+00,
        -1.60000000e+00,  -1.70000000e+00,  -1.66666670e+00,
        -1.60000000e+00,  -1.46666670e+00,  -1.50000000e+00,
        -1.65000000e+00,  -1.93333330e+00,  -1.90000000e+00,
        -1.98333330e+00,  -1.53333330e+00,  -1.60000000e+00,
        -1.71666670e+00,  -1.85000000e+00,  -1.80000000e+00,
        -1.60000000e+00,  -1.43333330e+00,  -1.43333330e+00,
        -1.45000000e+00,  -1.63333330e+00,  -1.68333330e+00,
        -1.56666670e+00,  -1.61666670e+00,  -1.63333330e+00,
        -1.70000000e+00,  -1.85000000e+00,  -1.85000000e+00,
        -1.78333330e+00,  -1.68333330e+00,  -1.66666670e+00,
        -1.53333330e+00,  -1.51666670e+00,  -1.70000000e+00,
        -1.65000000e+00,  -1.66666670e+00,  -1.70000000e+00,
        -1.80000000e+00,  -1.80000000e+00,  -1.80000000e+00,
        -1.80000000e+00,  -1.76666670e+00,  -1.80000000e+00,
        -1.83333330e+00,  -1.81666670e+00,  -1.78333330e+00,
        -1.73333330e+00,  -1.70000000e+00,  -1.75000000e+00,
        -1.66666670e+00,  -1.61666670e+00,  -1.53333330e+00,
        -1.65000000e+00,  -1.51666670e+00,  -1.50000000e+00,
        -1.43333330e+00,  -1.28333330e+00,  -1.66666670e+00,
        -1.75000000e+00,  -1.85000000e+00,  -1.88333330e+00,
        -1.90000000e+00,  -1.90000000e+00,  -1.90000000e+00,
        -1.91666670e+00,  -1.91666670e+00,  -1.88333330e+00,
        -1.91666670e+00,  -1.91666670e+00,  -1.90000000e+00,
        -1.90000000e+00,  -1.90000000e+00,  -1.86666670e+00,
        -1.81666670e+00,  -1.73333330e+00,  -1.80000000e+00,
        -1.90000000e+00,  -1.86666670e+00,  -1.80000000e+00,
        -1.80000000e+00,  -1.83333330e+00,  -1.61666670e+00,
        -1.38333330e+00,  -1.40000000e+00,  -1.50000000e+00,
        -1.76666670e+00,  -1.70000000e+00,  -1.70000000e+00,
        -1.60000000e+00,  -1.58333330e+00,  -1.60000000e+00,
        -1.66666670e+00,  -1.70000000e+00,  -1.70000000e+00,
        -1.63333330e+00,  -1.58333330e+00,  -1.68333330e+00,
        -1.66666670e+00,  -1.63333330e+00,  -1.70000000e+00,
        -1.55000000e+00,  -1.70000000e+00,  -1.73333330e+00,
        -1.68333330e+00,  -1.73333330e+00,  -1.83333330e+00,
        -1.85000000e+00,  -1.96666670e+00,  -1.95000000e+00,
        -1.86666670e+00,  -1.90000000e+00,  -2.16666670e+00,
        -2.55000000e+00,  -2.78333330e+00,  -2.83333330e+00,
        -2.00000000e+00,  -2.23333330e+00,  -2.26666670e+00,
        -2.11666670e+00,  -2.25000000e+00,  -2.33333330e+00,
        -8.00000000e-01,   5.33333330e-01,   2.66666670e-01,
        -1.15000000e+00,  -2.46666670e+00,  -2.62000000e+00,
        -2.82500000e+00,  -3.10000000e+00,  -3.32500000e+00,
        -3.15000000e+00,  -3.10000000e+00,  -3.07500000e+00,
        -2.97500000e+00,  -2.95000000e+00,  -2.62500000e+00,
        -2.10000000e+00,  -2.07500000e+00,  -2.17500000e+00,
        -2.05000000e+00,  -2.00000000e+00,  -2.00000000e+00,
        -2.00000000e+00,  -1.92500000e+00,  -1.90000000e+00,
        -1.87500000e+00,  -1.95000000e+00,  -1.90000000e+00,
        -1.97500000e+00,  -1.87500000e+00,  -1.82500000e+00,
        -1.80000000e+00,  -1.90000000e+00,  -1.97500000e+00,
        -2.05000000e+00,  -2.02500000e+00,  -1.77500000e+00,
        -1.70000000e+00,  -1.85000000e+00,  -1.80000000e+00,
        -1.85000000e+00,  -1.85000000e+00,  -1.77500000e+00,
        -1.70000000e+00,  -1.72500000e+00,  -1.75000000e+00,
        -1.80000000e+00,  -1.80000000e+00,  -1.80000000e+00,
        -1.82500000e+00,  -1.75000000e+00,  -1.77500000e+00,
        -1.80000000e+00,  -1.55000000e+00,  -1.12500000e+00,
        -1.85000000e+00,  -1.90000000e+00,  -1.85000000e+00,
        -1.90000000e+00,  -1.90000000e+00,  -1.90000000e+00,
        -1.90000000e+00,  -1.85000000e+00,  -1.85000000e+00,
        -1.90000000e+00,  -1.87500000e+00,  -1.85000000e+00,
        -1.85000000e+00,  -1.90000000e+00,  -2.10000000e+00,
        -2.15000000e+00,  -1.95000000e+00,  -1.75000000e+00,
        -1.72500000e+00,  -2.27500000e+00,  -2.37500000e+00,
        -2.47500000e+00,  -2.50000000e+00,  -2.35000000e+00,
        -2.10000000e+00,  -2.05000000e+00,  -1.85000000e+00,
        -2.10000000e+00,  -2.07500000e+00,  -2.20000000e+00,
        -2.45000000e+00,  -2.72500000e+00,  -3.07500000e+00,
        -3.07500000e+00,  -3.17500000e+00,  -3.30000000e+00,
        -3.35000000e+00,  -3.47500000e+00,  -3.30000000e+00,
         4.45000000e-01,   1.94000000e+00,   1.94000000e+00,
         1.94000000e+00,   1.94000000e+00,   1.35500000e+00,
        -2.25000000e+00,  -3.55000000e+00,  -2.60000000e+00,
        -1.65000000e+00,  -1.47500000e+00,  -1.50000000e+00,
        -1.55000000e+00,  -1.50000000e+00,  -1.30000000e+00,
        -1.10000000e+00,  -8.75000000e-01,  -8.25000000e-01,
        -8.50000000e-01,  -1.02500000e+00,  -1.20000000e+00,
        -1.32500000e+00,  -1.25000000e+00,  -1.07500000e+00,
        -1.07500000e+00,  -1.00000000e+00,  -1.07500000e+00,
        -1.30000000e+00,  -1.90000000e+00,  -1.87500000e+00,
        -1.85000000e+00,  -1.82500000e+00,  -1.75000000e+00,
        -1.80000000e+00,  -1.67500000e+00,  -1.67500000e+00,
        -1.65000000e+00,  -1.60000000e+00,  -1.65000000e+00,
        -1.72500000e+00,  -1.65000000e+00,  -1.65000000e+00,
        -1.67500000e+00,  -1.67500000e+00,  -1.75000000e+00,
        -1.65000000e+00,  -1.80000000e+00,  -1.67500000e+00,
        -1.50000000e+00,  -1.55000000e+00,  -1.45000000e+00,
        -1.70000000e+00,  -1.82500000e+00,  -1.77500000e+00,
        -1.72500000e+00,  -1.75000000e+00,  -1.57500000e+00,
        -1.45000000e+00,  -1.60000000e+00,  -1.70000000e+00,
        -1.77500000e+00,  -1.90000000e+00,  -2.00000000e+00,
        -2.10000000e+00,  -2.30000000e+00,  -2.60000000e+00,
        -2.87500000e+00,  -3.15000000e+00,  -3.20000000e+00,
        -3.27500000e+00,  -3.35000000e+00,  -3.60000000e+00,
        -3.65000000e+00,  -3.55000000e+00,  -3.57500000e+00,
        -3.65000000e+00,  -3.75000000e+00,  -3.72500000e+00,
        -4.05000000e+00,  -4.10000000e+00,  -4.12500000e+00,
        -4.10000000e+00,  -4.15000000e+00,  -4.37500000e+00,
        -4.60000000e+00,  -5.35000000e+00,  -5.40000000e+00,
        -5.05000000e+00,  -4.42500000e+00,  -4.23333330e+00,
        -3.92500000e+00,  -3.50000000e+00,  -3.40000000e+00,
        -3.32500000e+00,  -3.05000000e+00,  -2.87500000e+00,
        -2.95000000e+00,  -3.12500000e+00,  -3.45000000e+00,
        -3.40000000e+00,  -3.10000000e+00,  -2.85000000e+00,
        -3.05000000e+00,  -3.37500000e+00,  -4.10000000e+00,
        -4.32500000e+00,  -4.35000000e+00,  -4.65000000e+00,
        -4.52500000e+00,  -3.15000000e+00,  -3.60000000e+00,
        -3.45000000e+00,  -3.45000000e+00,  -3.37500000e+00,
        -3.45000000e+00,  -3.50000000e+00,  -3.27500000e+00,
        -3.20000000e+00,  -3.05000000e+00,  -2.00000000e+00,
        -5.00000000e-01,  -7.50000000e-02,  -8.50000000e-02,
         1.45000000e-01,   3.77500000e-01,   6.92500000e-01,
         9.30000000e-01,   1.17250000e+00,   1.48500000e+00,
         1.63500000e+00,   1.92500000e+00,   2.14000000e+00,
         2.28750000e+00,   2.44000000e+00,   2.51750000e+00,
         2.67000000e+00,   2.82500000e+00,   2.90250000e+00,
         3.05750000e+00,   3.17500000e+00,   3.29000000e+00,
         3.44500000e+00,   3.52500000e+00,   3.68500000e+00,
         3.81000000e+00,   3.93750000e+00,   4.11000000e+00,
         4.24500000e+00,   4.39000000e+00,   4.58500000e+00,
         4.68500000e+00,   4.89000000e+00,   5.04500000e+00,
         5.20750000e+00,   5.42500000e+00,   5.59000000e+00,
         5.76250000e+00,   5.99250000e+00,   6.11000000e+00,
         6.35000000e+00,   6.53500000e+00,   6.72250000e+00,
         6.97500000e+00,   7.17000000e+00,   7.37000000e+00,
         7.64000000e+00,   7.78000000e+00,   8.06000000e+00,
         8.27500000e+00,   8.49500000e+00,   8.79750000e+00,
         8.95000000e+00,   9.25750000e+00,   9.57500000e+00,
         9.73500000e+00,   1.00550000e+01,   1.02950000e+01,
         1.05450000e+01,   1.08725000e+01,   1.10350000e+01,
         1.13600000e+01,   1.16000000e+01,   1.18400000e+01,
         1.21575000e+01,   1.23950000e+01,   1.26250000e+01,
         1.29300000e+01,   1.30800000e+01,   1.33675000e+01,
         1.35700000e+01,   1.37625000e+01,   1.39950000e+01,
         1.41600000e+01,   1.43100000e+01,   1.45100000e+01,
         1.46050000e+01,   1.48000000e+01,   1.49500000e+01,
         1.51075000e+01,   1.53175000e+01,   1.54800000e+01,
         1.56450000e+01,   1.58650000e+01,   1.59700000e+01,
         1.61850000e+01,   1.63400000e+01,   1.64900000e+01,
         1.66775000e+01,   1.68050000e+01,   1.69250000e+01,
         1.70700000e+01,   1.71400000e+01,   1.72700000e+01,
         1.73550000e+01,   1.74450000e+01,   1.75500000e+01,
         1.76000000e+01,   1.76975000e+01,   1.77875000e+01,
         1.78325000e+01,   1.79150000e+01,   1.79800000e+01,
         1.80400000e+01,   1.81150000e+01,   1.81525000e+01,
         1.82200000e+01,   1.82700000e+01,   1.83200000e+01,
         1.83800000e+01,   1.84250000e+01,   1.84775000e+01,
         1.85375000e+01,   1.85700000e+01,   1.86300000e+01,
         1.86700000e+01,   1.87100000e+01,   1.87700000e+01,
         1.88150000e+01,   1.88500000e+01,   1.89000000e+01,
         1.89250000e+01,   1.89650000e+01,   1.89950000e+01,
         1.90175000e+01,   1.90500000e+01,   1.90700000e+01,
         1.90900000e+01,   1.91175000e+01,   1.91300000e+01,
         1.91550000e+01,   1.91700000e+01,   1.91900000e+01,
         1.92150000e+01,   1.92250000e+01,   1.92500000e+01,
         1.92700000e+01,   1.92800000e+01,   1.93050000e+01,
         1.93200000e+01,   1.93375000e+01,   1.93600000e+01,
         1.93700000e+01,   1.93950000e+01,   1.94100000e+01,
         1.94250000e+01,   1.94500000e+01,   1.94650000e+01,
         1.94850000e+01,   1.95050000e+01,   1.95150000e+01,
         1.95375000e+01,   1.95550000e+01,   1.95700000e+01,
         1.95900000e+01,   1.96050000e+01,   1.96250000e+01,
         1.96450000e+01,   1.96550000e+01,   1.96750000e+01,
         1.96900000e+01,   1.97050000e+01,   1.97250000e+01,
         1.97400000e+01,   1.97550000e+01,   1.97775000e+01,
         1.97900000e+01,   1.98100000e+01,   1.98250000e+01,
         1.98400000e+01,   1.98600000e+01,   1.98700000e+01,
         1.98900000e+01,   1.99100000e+01,   1.99200000e+01,
         1.99400000e+01,   1.99550000e+01,   1.99700000e+01,
         1.99850000e+01,   1.99950000e+01,   2.00100000e+01,
         2.00250000e+01,   2.00350000e+01,   2.00525000e+01,
         2.00650000e+01,   2.00750000e+01,   2.00875000e+01,
         2.00950000e+01,   2.01075000e+01,   2.01200000e+01,
         2.01280000e+01,   2.01400000e+01,   2.01500000e+01,
         2.01600000e+01,   2.01700000e+01,   2.01766670e+01,
         2.01900000e+01,   2.02033330e+01,   2.02116670e+01,
         2.02266670e+01,   2.02366670e+01,   2.02466670e+01,
         2.02600000e+01,   2.02666670e+01,   2.02800000e+01,
         2.02900000e+01,   2.03000000e+01,   2.03133330e+01,
         2.03233330e+01,   2.03316670e+01,   2.03450000e+01,
         2.03500000e+01,   2.03650000e+01,   2.03766670e+01,
         2.03866670e+01,   2.04000000e+01,   2.04066670e+01,
         2.04200000e+01,   2.04333330e+01,   2.04400000e+01,
         2.04533330e+01,   2.04633330e+01,   2.04733330e+01,
         2.04850000e+01,   2.04916670e+01,   2.05033330e+01,
         2.05133330e+01,   2.05233330e+01,   2.05350000e+01,
         2.05466670e+01,   2.05550000e+01,   2.05666670e+01,
         2.05733330e+01,   2.05850000e+01,   2.05933330e+01,
         2.06033330e+01,   2.06150000e+01,   2.06266670e+01,
         2.06366670e+01,   2.06500000e+01,   2.06566670e+01,
         2.06700000e+01,   2.06766670e+01,   2.06900000e+01,
         2.07000000e+01,   2.07133330e+01,   2.07200000e+01,
         2.07333330e+01,   2.07400000e+01,   2.07516670e+01,
         2.07600000e+01,   2.07700000e+01,   2.07833330e+01,
         2.07900000e+01,   2.08033330e+01,   2.08166670e+01,
         2.08233330e+01,   2.08350000e+01,   2.08433330e+01,
         2.08533330e+01,   2.08683330e+01,   2.08750000e+01,
         2.08900000e+01,   2.08966670e+01,   2.09066670e+01,
         2.09200000e+01,   2.09300000e+01,   2.09400000e+01,
         2.09500000e+01,   2.09566670e+01,   2.09700000e+01,
         2.09800000e+01,   2.09900000e+01,   2.10016670e+01,
         2.10100000e+01,   2.10200000e+01,   2.10333330e+01,
         2.10400000e+01,   2.10516670e+01,   2.10600000e+01,
         2.10700000e+01,   2.10816670e+01,   2.10900000e+01,
         2.11000000e+01,   2.11116670e+01,   2.11183330e+01,
         2.11300000e+01,   2.11400000e+01,   2.11500000e+01,
         2.11633330e+01,   2.11700000e+01,   2.11800000e+01,
         2.11933330e+01,   2.12000000e+01,   2.12150000e+01,
         2.12233330e+01,   2.12333330e+01,   2.12466670e+01,
         2.12516670e+01,   2.12633330e+01,   2.12760000e+01,
         2.12816670e+01,   2.12966670e+01,   2.13033330e+01,
         2.13125000e+01,   2.13250000e+01,   2.13300000e+01,
         2.13425000e+01,   2.13500000e+01,   2.13625000e+01,
         2.13750000e+01,   2.13850000e+01,   2.13950000e+01,
         2.14075000e+01,   2.14150000e+01,  -9.99000000e+02])

In [ ]: