Full Inversion Method for a single chromosome using OpenMiChroM

This tutorial enables performing Full Inversion Method using Adam (First-order optimization)

The first step is to import the OpenMiChroM modules.

To install OpenMM and OpenMiChroM, follow the instalation guide

The inputs and apps used in this tutorial can be downloaded here

Adam optimization is available in OpenMichroM versions >= 1.0.5

from OpenMiChroM.ChromDynamics import MiChroM # OpenMiChroM simulation module
from OpenMiChroM.Optimization import AdamTraining # Adam optimization module
from OpenMiChroM.CndbTools import cndbTools # Analysis' tools module

# modules to load and plot
import matplotlib.pyplot as plt
import matplotlib as mpl
from sklearn.preprocessing import normalize
import numpy as np
import pandas as pd
import h5py

To create a model trained by the Full Inversion Method using Adam optimization 3 files are required:

1. The experimental Hi-C matrix in text format (dense file)

2. The sequence file

3. Initial force field

In the next steps we will create/extract theses files.

1. Extract the experimental Hi-C matrix

A Hi-C file is required for the analysis and training of the Full inversion optmization. The file format chosen here is a matrix .txt file (we will call it the dense file).

For this tutorial, we will use chromosome 10 from GM12878 cell line in 100 kb resolution.

To extract it from the .hic file we can use juicer_tools with this command:

java -jar juicer_tools_1.22.01.jar dump observed Balanced -d https://hicfiles.s3.amazonaws.com/hiseq/gm12878/in-situ/combined.hic 10 10 BP 100000 input/chr10_100k.dense

%%bash
java -jar apps/juicer_tools_1.22.01.jar dump observed Balanced -d https://hicfiles.s3.amazonaws.com/hiseq/gm12878/in-situ/combined.hic 10 10 BP 100000 input/chr10_100k.dense

WARN [2022-09-14T13:46:43,872]  [Globals.java:138] [main]  Development mode is enabled
INFO [2022-09-14T13:46:47,853]  [DirectoryManager.java:179] [main]  IGV Directory: /Users/mm146/igv
INFO [2022-09-14T13:46:48,195]  [HttpUtils.java:937] [main]  Range-byte request succeeded

This command accesses part of the .hic file from the web and extracts the chromosome 10 in .dense format to the folder “input/”.

You can get more information about it at the JuicerTools documentation.

Visualize the .dense file for inspection

filename = 'input/chr10_100k.dense'
hic_file = np.loadtxt(filename)

r=np.triu(hic_file, k=1)
r[np.isnan(r)]= 0.0
r = normalize(r, axis=1, norm='max')
rd = np.transpose(r)
r=r+rd + np.diag(np.ones(len(r)))
print("number of beads: ", len(r))
plt.matshow(r,norm=mpl.colors.LogNorm(vmin=0.0001, vmax=r.max()),cmap="Reds")
plt.colorbar()

number of beads:  1356

<matplotlib.colorbar.Colorbar at 0x7f81cf0f90d0>

The Hi-C map has a resolution of 100 kb per bead, so this chromosome 10 model has a polymer chain with a total of 1356 beads.

2. Create the sequence file

The second required input is the sequence file.

The sequence file is a 2 colunms text file that labels each bead of your polymer chain.

In OpenMiChroM, each bead must be assigned a chromatin type. In the full inversion method each pair of interaction is trained individually and therefore each bead is unique.

To create the sequence file for chr10 (1356 beads) use the following bash command:

%%bash
rm input/seq_chr10
for i in {1..1356}; do echo "$i Bead$i" >> input/seq_chr10; done

This command creates a file inside the folder “input” named seq_chr10.

The first column is the index for each bead and the second column is the flavor (or type) for each bead. In this tutorial, we chose the names as beadX, where X is the index. You can choose any name for the beads as long as each name is unique. Additionally, you will need to have all names in the force field file (next step).

%%bash
head input/seq_chr10

1 Bead1
2 Bead2
3 Bead3
4 Bead4
5 Bead5
6 Bead6
7 Bead7
8 Bead8
9 Bead9
10 Bead10

3. Create the initial force field file

Last but not least, we will create the third required file, the initial force field. The force field is represented by a \(N \times N\) matrix, where it pair \((i,j)\) is the interaction between beads \(i\) and \(j\).

This file format is a csv-file with 1 header line, delimited by comma (“,”).

To be simple, as a tutorial must be, the initial force field will have all entries equal to \(0\).

(In the end of this tutorial, we present a smarter method to fix certain interactions and to solve the centromere problem)

#read the sequence file
seq = np.loadtxt("input/seq_chr10", dtype=str)[:,1]
pol_size = len(seq)
#create a NxN matrix fill with zeros, where N is the chromossome lenght
data = np.zeros((pol_size,pol_size))
#transform it in a Pandas dataframe using the sequence names as header
lamb = pd.DataFrame(data, columns=seq)
#save to a csv file format in the  "input" folder
lamb.to_csv("input/lambda_1", index=False)
lamb

	Bead1	Bead2	Bead3	Bead4	Bead5	Bead6	Bead7	Bead8	Bead9	Bead10	...	Bead1347	Bead1348	Bead1349	Bead1350	Bead1351	Bead1352	Bead1353	Bead1354	Bead1355	Bead1356
0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
2	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
3	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
4	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
1351	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1352	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1353	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1354	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1355	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0

1356 rows × 1356 columns

With all required files inside the “input” folder, we can start to make iterations to converge the force field. The pipeline for this optimization is the following:

|b5f296c0a3244465a3f8445a13c989ed|

The next step is to carry out simulations to obtain a contact map from an ensemble of simulated structures. Multiple simulations are required in order to get a converged contact map in each iteration.

Here we will show how to run a single chromosome simulation and run one iteration of the Adam optimization procedure.

To initiate MiChroM, we need to setup some variables:

time_step=0.01 (the time step used for integration, default is 0.01) temperature=1 (set the temperature of your simulation) name=‘opt_chr10_100K’ (the simulation name)

sim = MiChroM(name='opt_chr10_100K',temperature=1.0, time_step=0.01)

    ***************************************************************************************
     **** **** *** *** *** *** *** *** OpenMiChroM-1.0.5 *** *** *** *** *** *** **** ****

         OpenMiChroM is a Python library for performing chromatin dynamics simulations.
                            OpenMiChroM uses the OpenMM Python API,
                employing the MiChroM (Minimal Chromatin Model) energy function.
      The chromatin dynamics simulations generate an ensemble of 3D chromosomal structures
      that are consistent with experimental Hi-C maps, also allows simulations of a single
                 or multiple chromosome chain using High-Performance Computing
                            in different platforms (GPUs and CPUs).
         OpenMiChroM documentation is available at https://open-michrom.readthedocs.io

         OpenMiChroM is described in: Oliveira Junior, A. B & Contessoto, V, G et. al.
      A Scalable Computational Approach for Simulating Complexes of Multiple Chromosomes.
                  Journal of Molecular Biology. doi:10.1016/j.jmb.2020.10.034.
                                              and
                                 Oliveira Junior, A. B. et al.
     Chromosome Modeling on Downsampled Hi-C Maps Enhances the Compartmentalization Signal.
                        J. Phys. Chem. B, doi:10.1021/acs.jpcb.1c04174.

                    Copyright (c) 2022, The OpenMiChroM development team at
                                        Rice University
    ***************************************************************************************

Now we need to setup the platform that we will use. The options are:

platform=“cuda” (remember that you need to install CUDA in your system) GPU=“0” (optional) (if you have more than one GPU device, you can set which gpu you want [“0”, “1”,…,“n”]) platform=“cpu” platform=“opencl”

sim.setup(platform="opencl")

Set the folder name where the output will be saved

sim.saveFolder('iteration_1')

The next step is to setup your chromosome sequence and the initial configuration

mychro = sim.createSpringSpiral(ChromSeq="input/seq_chr10")

Load the initial structure into the “sim” object

sim.loadStructure(mychro, center=True)

Now it is time to include the force field in the simulation object “sim”.

Let’s separate the forces in two sets:

Homopolymer Potentials

sim.addFENEBonds(kfb=30.0)
sim.addAngles(ka=2.0)
sim.addRepulsiveSoftCore(Ecut=4.0)

Chromosome Potential

In this tutorial, we will use the CustomTypes potential. Here you need to insert the force field file that contains a matrix of interactions for each bead. For each iteration of the optimization of the parameters, a new force field file will be created.

To check that, you can look on the documentation https://open-michrom.readthedocs.io/en/latest/OpenMiChroM.html#OpenMiChroM.ChromDynamics.MiChroM.addCustomTypes

sim.addCustomTypes(mu=3.22, rc = 1.78, TypesTable='input/lambda_1')

Note: these values for \(mu\) and \(r_c\) were calculated for human GM12878 cells and can be changed for other species.

The last potential to be added is the spherical restraint in order to collapse the initial structure

sim.addFlatBottomHarmonic(kr=5*10**-3, n_rad=8.0)

Now we will run a short simulation in order to get a collapsed structure.

There are two variables that control the chromosome simulation steps:

block: The number of time steps performed in each cycle (n_blocks) n_blocks: The number of blocks that will be perfomed.

In this example, to perfom the collapse we will run \(5\times10^2 \times 10^3 = 5\times10^5\) time steps

block = 5*10**2
n_blocks = 10**3

We can save the radius of gyration of each block to observe the convergence into the collapsed state (the time required here depends on the size of your chromosome)

rg = []

for _ in range(n_blocks):
    sim.runSimBlock(block, increment=False)
    rg.append(sim.chromRG())

#save the collapsed structure in pdb format for inspection
sim.saveStructure(mode = 'pdb')

Number of exceptions: 1355
adding force  FENEBond 0
adding force  AngleForce 1
Add exclusions for RepulsiveSoftCore force
adding force  RepulsiveSoftCore 2
Add exclusions for CustomTypes force
adding force  CustomTypes 3
adding force  FlatBottomHarmonic 4
Positions...
 loaded!
potential energy is 31.227323
bl=0 pos[1]=[103.8 -6.9 1.7] dr=1.22 t=0.0ps kin=1.51 pot=31.48 Rg=73.005 SPS=2934
bl=0 pos[1]=[103.3 -6.6 4.0] dr=1.86 t=0.0ps kin=1.61 pot=31.24 Rg=71.826 SPS=14103
bl=0 pos[1]=[101.9 -6.3 5.2] dr=1.91 t=0.0ps kin=1.72 pot=30.91 Rg=70.422 SPS=14505
bl=0 pos[1]=[100.7 -6.0 4.7] dr=1.96 t=0.0ps kin=1.85 pot=30.57 Rg=68.921 SPS=13752
bl=0 pos[1]=[99.3 -4.6 4.8] dr=1.98 t=0.0ps kin=1.89 pot=30.12 Rg=67.367 SPS=14440
bl=0 pos[1]=[97.6 -2.7 6.1] dr=1.97 t=0.0ps kin=1.93 pot=29.76 Rg=65.807 SPS=15452
bl=0 pos[1]=[94.9 -0.6 7.3] dr=1.95 t=0.0ps kin=1.88 pot=29.37 Rg=64.273 SPS=15480
bl=0 pos[1]=[92.4 -0.5 7.3] dr=1.90 t=0.0ps kin=1.90 pot=28.95 Rg=62.808 SPS=17235
bl=0 pos[1]=[90.7 -0.5 8.5] dr=1.88 t=0.0ps kin=1.85 pot=28.62 Rg=61.385 SPS=17314
bl=0 pos[1]=[87.8 -0.5 9.4] dr=1.84 t=0.0ps kin=1.86 pot=28.22 Rg=60.033 SPS=17292
bl=0 pos[1]=[84.3 0.1 9.3] dr=1.85 t=0.0ps kin=1.84 pot=27.90 Rg=58.697 SPS=16504
bl=0 pos[1]=[81.1 -0.8 10.5] dr=1.82 t=0.0ps kin=1.83 pot=27.56 Rg=57.382 SPS=17570
bl=0 pos[1]=[80.2 -0.7 11.4] dr=1.84 t=0.0ps kin=1.83 pot=27.28 Rg=56.065 SPS=17480
bl=0 pos[1]=[76.9 -0.5 12.9] dr=1.81 t=0.0ps kin=1.75 pot=26.99 Rg=54.792 SPS=17399
bl=0 pos[1]=[74.3 -1.3 13.2] dr=1.76 t=0.0ps kin=1.81 pot=26.73 Rg=53.565 SPS=17461
bl=0 pos[1]=[72.6 -2.9 13.8] dr=1.75 t=0.0ps kin=1.77 pot=26.44 Rg=52.331 SPS=17060
bl=0 pos[1]=[71.5 -3.3 13.5] dr=1.75 t=0.0ps kin=1.74 pot=26.15 Rg=51.107 SPS=17400
bl=0 pos[1]=[68.8 -3.8 11.8] dr=1.68 t=0.0ps kin=1.72 pot=25.95 Rg=49.951 SPS=17750
bl=0 pos[1]=[66.9 -4.5 9.5] dr=1.65 t=0.0ps kin=1.72 pot=25.77 Rg=48.878 SPS=17742
bl=0 pos[1]=[64.2 -3.6 8.1] dr=1.65 t=0.0ps kin=1.75 pot=25.55 Rg=47.822 SPS=16747
bl=0 pos[1]=[61.7 -4.1 7.5] dr=1.66 t=0.0ps kin=1.74 pot=25.36 Rg=46.750 SPS=17638
bl=0 pos[1]=[60.6 -6.5 6.7] dr=1.62 t=0.0ps kin=1.71 pot=25.15 Rg=45.736 SPS=17850
bl=0 pos[1]=[60.4 -5.4 5.8] dr=1.59 t=0.0ps kin=1.73 pot=24.94 Rg=44.822 SPS=18043
bl=0 pos[1]=[58.9 -5.2 5.9] dr=1.58 t=0.0ps kin=1.68 pot=24.76 Rg=43.940 SPS=17832
bl=0 pos[1]=[56.6 -6.0 7.6] dr=1.57 t=0.0ps kin=1.69 pot=24.60 Rg=43.033 SPS=17849
bl=0 pos[1]=[55.1 -6.4 9.9] dr=1.57 t=0.0ps kin=1.61 pot=24.47 Rg=42.106 SPS=17763
bl=0 pos[1]=[52.5 -6.8 9.0] dr=1.55 t=0.0ps kin=1.62 pot=24.29 Rg=41.198 SPS=17837
bl=0 pos[1]=[51.0 -6.2 8.0] dr=1.54 t=0.0ps kin=1.67 pot=24.18 Rg=40.336 SPS=17422
bl=0 pos[1]=[48.7 -5.3 7.5] dr=1.56 t=0.0ps kin=1.63 pot=24.06 Rg=39.490 SPS=17423
bl=0 pos[1]=[46.6 -4.3 7.5] dr=1.54 t=0.0ps kin=1.62 pot=23.91 Rg=38.646 SPS=15926
bl=0 pos[1]=[44.4 -3.6 7.2] dr=1.51 t=0.0ps kin=1.61 pot=23.80 Rg=37.834 SPS=16716
bl=0 pos[1]=[42.9 -2.7 8.8] dr=1.51 t=0.0ps kin=1.61 pot=23.65 Rg=36.994 SPS=17139
bl=0 pos[1]=[43.0 -3.7 9.1] dr=1.50 t=0.0ps kin=1.56 pot=23.51 Rg=36.191 SPS=17881
bl=0 pos[1]=[43.1 -3.6 9.2] dr=1.49 t=0.0ps kin=1.60 pot=23.43 Rg=35.413 SPS=16878
bl=0 pos[1]=[43.1 -3.9 10.7] dr=1.47 t=0.0ps kin=1.63 pot=23.29 Rg=34.688 SPS=16727
bl=0 pos[1]=[44.2 -3.5 10.3] dr=1.43 t=0.0ps kin=1.60 pot=23.21 Rg=34.001 SPS=17193
bl=0 pos[1]=[43.0 -3.0 10.9] dr=1.44 t=0.0ps kin=1.53 pot=23.17 Rg=33.324 SPS=17711
bl=0 pos[1]=[41.3 -2.6 11.7] dr=1.41 t=0.0ps kin=1.58 pot=23.07 Rg=32.712 SPS=17386
bl=0 pos[1]=[40.7 -3.8 12.6] dr=1.41 t=0.0ps kin=1.58 pot=23.00 Rg=32.084 SPS=16220
bl=0 pos[1]=[40.2 -5.3 12.9] dr=1.41 t=0.0ps kin=1.53 pot=22.89 Rg=31.436 SPS=16462
bl=0 pos[1]=[39.1 -6.7 13.5] dr=1.41 t=0.0ps kin=1.57 pot=22.81 Rg=30.814 SPS=17130
bl=0 pos[1]=[38.7 -5.6 13.3] dr=1.39 t=0.0ps kin=1.55 pot=22.76 Rg=30.207 SPS=17235
bl=0 pos[1]=[38.4 -4.0 12.6] dr=1.40 t=0.0ps kin=1.56 pot=22.69 Rg=29.637 SPS=16268
bl=0 pos[1]=[36.8 -2.1 13.1] dr=1.38 t=0.0ps kin=1.55 pot=22.68 Rg=29.117 SPS=18108
bl=0 pos[1]=[36.0 0.1 13.2] dr=1.40 t=0.0ps kin=1.57 pot=22.60 Rg=28.609 SPS=17989
bl=0 pos[1]=[33.4 1.6 12.9] dr=1.39 t=0.0ps kin=1.56 pot=22.54 Rg=28.120 SPS=18060
bl=0 pos[1]=[30.6 2.2 14.2] dr=1.37 t=0.0ps kin=1.48 pot=22.50 Rg=27.635 SPS=17033
bl=0 pos[1]=[29.9 2.8 14.4] dr=1.39 t=0.0ps kin=1.55 pot=22.43 Rg=27.117 SPS=16332
bl=0 pos[1]=[29.5 2.4 14.0] dr=1.36 t=0.0ps kin=1.51 pot=22.43 Rg=26.613 SPS=17331
bl=0 pos[1]=[30.8 1.9 13.2] dr=1.32 t=0.0ps kin=1.55 pot=22.38 Rg=26.133 SPS=17662
bl=0 pos[1]=[32.0 1.7 13.8] dr=1.36 t=0.0ps kin=1.58 pot=22.33 Rg=25.700 SPS=17185
bl=0 pos[1]=[31.7 2.6 13.9] dr=1.31 t=0.0ps kin=1.52 pot=22.30 Rg=25.307 SPS=15488
bl=0 pos[1]=[30.2 2.1 15.3] dr=1.34 t=0.0ps kin=1.54 pot=22.25 Rg=24.873 SPS=16077
bl=0 pos[1]=[29.4 1.7 13.5] dr=1.33 t=0.0ps kin=1.52 pot=22.22 Rg=24.431 SPS=15575
bl=0 pos[1]=[31.1 2.0 12.2] dr=1.32 t=0.0ps kin=1.52 pot=22.19 Rg=24.015 SPS=15193
bl=0 pos[1]=[31.7 1.6 12.6] dr=1.30 t=0.0ps kin=1.53 pot=22.19 Rg=23.595 SPS=15574
bl=0 pos[1]=[31.6 1.1 14.4] dr=1.34 t=0.0ps kin=1.50 pot=22.12 Rg=23.179 SPS=16172
bl=0 pos[1]=[29.2 0.6 15.1] dr=1.31 t=0.0ps kin=1.47 pot=22.09 Rg=22.806 SPS=15418
bl=0 pos[1]=[27.2 1.9 12.8] dr=1.28 t=0.0ps kin=1.51 pot=22.07 Rg=22.433 SPS=16059
bl=0 pos[1]=[24.8 3.5 12.5] dr=1.33 t=0.0ps kin=1.55 pot=22.01 Rg=22.041 SPS=15369
bl=0 pos[1]=[23.3 4.5 12.7] dr=1.35 t=0.0ps kin=1.56 pot=21.98 Rg=21.686 SPS=15285
bl=0 pos[1]=[23.1 4.3 12.9] dr=1.34 t=0.0ps kin=1.53 pot=22.04 Rg=21.364 SPS=17384
bl=0 pos[1]=[21.8 4.5 12.6] dr=1.33 t=0.0ps kin=1.55 pot=22.04 Rg=21.049 SPS=16052
bl=0 pos[1]=[20.6 4.8 12.1] dr=1.33 t=0.0ps kin=1.56 pot=21.99 Rg=20.733 SPS=17021
bl=0 pos[1]=[20.1 3.9 10.2] dr=1.31 t=0.0ps kin=1.54 pot=21.91 Rg=20.416 SPS=15706
bl=0 pos[1]=[19.4 1.9 9.6] dr=1.31 t=0.0ps kin=1.47 pot=21.93 Rg=20.113 SPS=17062
bl=0 pos[1]=[19.3 3.0 9.5] dr=1.27 t=0.0ps kin=1.47 pot=21.89 Rg=19.833 SPS=16363
bl=0 pos[1]=[17.4 4.3 9.7] dr=1.28 t=0.0ps kin=1.51 pot=21.91 Rg=19.548 SPS=14951
bl=0 pos[1]=[15.7 4.1 8.9] dr=1.28 t=0.0ps kin=1.53 pot=21.84 Rg=19.320 SPS=16253
bl=0 pos[1]=[15.5 3.3 7.9] dr=1.29 t=0.0ps kin=1.48 pot=21.84 Rg=19.095 SPS=16024
bl=0 pos[1]=[14.2 2.0 8.6] dr=1.25 t=0.0ps kin=1.50 pot=21.78 Rg=18.884 SPS=15265
bl=0 pos[1]=[13.8 1.9 8.8] dr=1.24 t=0.0ps kin=1.45 pot=21.84 Rg=18.707 SPS=16479
bl=0 pos[1]=[14.4 1.5 7.5] dr=1.29 t=0.0ps kin=1.48 pot=21.83 Rg=18.521 SPS=15760
bl=0 pos[1]=[14.1 -0.5 6.7] dr=1.32 t=0.0ps kin=1.50 pot=21.77 Rg=18.294 SPS=15834
bl=0 pos[1]=[14.2 -3.6 7.3] dr=1.34 t=0.0ps kin=1.53 pot=21.77 Rg=18.027 SPS=15507
bl=0 pos[1]=[13.5 -5.0 7.9] dr=1.31 t=0.0ps kin=1.51 pot=21.76 Rg=17.776 SPS=14977
bl=0 pos[1]=[11.7 -6.0 9.5] dr=1.30 t=0.0ps kin=1.52 pot=21.74 Rg=17.512 SPS=15046
bl=0 pos[1]=[12.5 -6.3 9.7] dr=1.28 t=0.0ps kin=1.46 pot=21.76 Rg=17.308 SPS=15968
bl=0 pos[1]=[10.5 -5.7 8.9] dr=1.28 t=0.0ps kin=1.50 pot=21.71 Rg=17.141 SPS=15638
bl=0 pos[1]=[11.2 -5.6 10.2] dr=1.26 t=0.0ps kin=1.49 pot=21.73 Rg=16.963 SPS=14641
bl=0 pos[1]=[9.6 -5.5 10.3] dr=1.26 t=0.0ps kin=1.50 pot=21.72 Rg=16.799 SPS=15726
bl=0 pos[1]=[7.6 -3.7 10.2] dr=1.29 t=0.0ps kin=1.50 pot=21.73 Rg=16.645 SPS=15978
bl=0 pos[1]=[8.0 -3.3 10.3] dr=1.25 t=0.0ps kin=1.48 pot=21.75 Rg=16.468 SPS=14907
bl=0 pos[1]=[8.0 -3.3 9.3] dr=1.26 t=0.0ps kin=1.48 pot=21.76 Rg=16.306 SPS=15538
bl=0 pos[1]=[8.2 -2.9 9.0] dr=1.27 t=0.0ps kin=1.49 pot=21.70 Rg=16.155 SPS=15881
bl=0 pos[1]=[8.3 -2.9 8.4] dr=1.26 t=0.0ps kin=1.45 pot=21.74 Rg=16.036 SPS=15195
bl=0 pos[1]=[9.6 -2.1 8.9] dr=1.28 t=0.0ps kin=1.54 pot=21.76 Rg=15.930 SPS=15667
bl=0 pos[1]=[9.7 -4.0 8.6] dr=1.30 t=0.0ps kin=1.54 pot=21.75 Rg=15.824 SPS=15813
bl=0 pos[1]=[11.9 -1.6 10.9] dr=1.28 t=0.0ps kin=1.53 pot=21.77 Rg=15.702 SPS=15025
bl=0 pos[1]=[10.9 -0.3 9.8] dr=1.26 t=0.0ps kin=1.51 pot=21.74 Rg=15.585 SPS=15269
bl=0 pos[1]=[10.6 -1.0 9.6] dr=1.30 t=0.0ps kin=1.50 pot=21.71 Rg=15.463 SPS=14977
bl=0 pos[1]=[10.5 -1.5 9.9] dr=1.29 t=0.0ps kin=1.54 pot=21.66 Rg=15.330 SPS=15329
bl=0 pos[1]=[9.4 0.1 11.2] dr=1.28 t=0.0ps kin=1.51 pot=21.69 Rg=15.175 SPS=15634
bl=0 pos[1]=[9.3 0.6 13.1] dr=1.25 t=0.0ps kin=1.50 pot=21.70 Rg=15.047 SPS=14532
bl=0 pos[1]=[10.9 0.2 12.8] dr=1.25 t=0.0ps kin=1.52 pot=21.68 Rg=14.976 SPS=15172
bl=0 pos[1]=[9.0 -2.2 12.4] dr=1.26 t=0.0ps kin=1.55 pot=21.71 Rg=14.872 SPS=15249
bl=0 pos[1]=[8.2 -1.7 12.3] dr=1.26 t=0.0ps kin=1.51 pot=21.69 Rg=14.751 SPS=14576
bl=0 pos[1]=[8.6 -2.2 12.1] dr=1.27 t=0.0ps kin=1.50 pot=21.70 Rg=14.626 SPS=15352
bl=0 pos[1]=[9.8 -2.1 13.5] dr=1.28 t=0.0ps kin=1.50 pot=21.67 Rg=14.539 SPS=14347
bl=0 pos[1]=[11.4 -1.8 12.7] dr=1.29 t=0.0ps kin=1.52 pot=21.67 Rg=14.481 SPS=15101
bl=0 pos[1]=[11.1 -1.3 12.1] dr=1.27 t=0.0ps kin=1.55 pot=21.69 Rg=14.404 SPS=15225
bl=0 pos[1]=[10.1 -0.9 11.6] dr=1.28 t=0.0ps kin=1.56 pot=21.66 Rg=14.341 SPS=14536
bl=0 pos[1]=[7.9 -1.6 10.1] dr=1.25 t=0.0ps kin=1.49 pot=21.72 Rg=14.295 SPS=15357
bl=0 pos[1]=[8.2 -3.5 8.5] dr=1.25 t=0.0ps kin=1.52 pot=21.67 Rg=14.239 SPS=14538
bl=0 pos[1]=[8.3 -2.7 5.7] dr=1.25 t=0.0ps kin=1.49 pot=21.68 Rg=14.166 SPS=14819
bl=0 pos[1]=[9.5 -2.0 4.7] dr=1.25 t=0.0ps kin=1.44 pot=21.62 Rg=14.074 SPS=14742
bl=0 pos[1]=[8.7 -1.1 7.0] dr=1.23 t=0.0ps kin=1.50 pot=21.60 Rg=13.981 SPS=14529
bl=0 pos[1]=[9.7 -2.3 6.4] dr=1.25 t=0.0ps kin=1.47 pot=21.61 Rg=13.903 SPS=15026
bl=0 pos[1]=[10.0 -2.8 6.1] dr=1.27 t=0.0ps kin=1.46 pot=21.66 Rg=13.817 SPS=14228
bl=0 pos[1]=[10.7 -4.1 5.2] dr=1.28 t=0.0ps kin=1.46 pot=21.60 Rg=13.703 SPS=14853
bl=0 pos[1]=[10.4 -4.1 5.6] dr=1.23 t=0.0ps kin=1.45 pot=21.65 Rg=13.621 SPS=14234
bl=0 pos[1]=[9.4 -4.5 6.4] dr=1.25 t=0.0ps kin=1.48 pot=21.63 Rg=13.559 SPS=14643
bl=0 pos[1]=[9.4 -3.9 6.1] dr=1.26 t=0.0ps kin=1.49 pot=21.58 Rg=13.454 SPS=14587
bl=0 pos[1]=[9.4 -3.2 5.7] dr=1.27 t=0.0ps kin=1.49 pot=21.68 Rg=13.361 SPS=14174
bl=0 pos[1]=[8.6 -1.6 5.6] dr=1.26 t=0.0ps kin=1.48 pot=21.64 Rg=13.247 SPS=14783
bl=0 pos[1]=[8.9 -2.4 5.0] dr=1.23 t=0.0ps kin=1.48 pot=21.65 Rg=13.137 SPS=13895
bl=0 pos[1]=[8.6 -2.9 5.4] dr=1.29 t=0.0ps kin=1.50 pot=21.63 Rg=13.033 SPS=15268
bl=0 pos[1]=[8.1 -2.7 5.5] dr=1.23 t=0.0ps kin=1.47 pot=21.59 Rg=13.018 SPS=14657
bl=0 pos[1]=[9.0 -3.4 5.6] dr=1.24 t=0.0ps kin=1.50 pot=21.57 Rg=13.020 SPS=15318
bl=0 pos[1]=[10.3 -3.4 5.2] dr=1.26 t=0.0ps kin=1.47 pot=21.62 Rg=12.993 SPS=15619
bl=0 pos[1]=[12.2 -3.1 4.5] dr=1.24 t=0.0ps kin=1.44 pot=21.62 Rg=12.958 SPS=15009
bl=0 pos[1]=[13.1 -4.2 3.2] dr=1.27 t=0.0ps kin=1.47 pot=21.60 Rg=12.973 SPS=16050
bl=0 pos[1]=[15.7 -4.1 2.5] dr=1.22 t=0.0ps kin=1.53 pot=21.61 Rg=12.973 SPS=15933
bl=0 pos[1]=[16.6 -2.7 2.8] dr=1.27 t=0.0ps kin=1.60 pot=21.63 Rg=12.967 SPS=14931
bl=0 pos[1]=[16.6 -2.4 3.0] dr=1.31 t=0.0ps kin=1.58 pot=21.60 Rg=12.976 SPS=15430
bl=0 pos[1]=[16.8 -3.2 3.1] dr=1.27 t=0.0ps kin=1.52 pot=21.59 Rg=12.996 SPS=15985
bl=0 pos[1]=[15.8 -3.3 2.8] dr=1.25 t=0.0ps kin=1.51 pot=21.67 Rg=13.008 SPS=14636
bl=0 pos[1]=[14.6 -2.7 3.2] dr=1.26 t=0.0ps kin=1.53 pot=21.64 Rg=13.001 SPS=14895
bl=0 pos[1]=[15.8 -3.8 3.4] dr=1.24 t=0.0ps kin=1.53 pot=21.60 Rg=13.052 SPS=14783
bl=0 pos[1]=[13.8 -4.0 3.2] dr=1.30 t=0.0ps kin=1.53 pot=21.61 Rg=13.043 SPS=15287
bl=0 pos[1]=[13.0 -4.9 2.2] dr=1.28 t=0.0ps kin=1.54 pot=21.62 Rg=12.968 SPS=16326
bl=0 pos[1]=[14.2 -5.7 2.7] dr=1.31 t=0.0ps kin=1.49 pot=21.63 Rg=12.896 SPS=14765
bl=0 pos[1]=[14.1 -6.7 2.3] dr=1.32 t=0.0ps kin=1.49 pot=21.58 Rg=12.902 SPS=14747
bl=0 pos[1]=[13.1 -7.2 3.1] dr=1.28 t=0.0ps kin=1.47 pot=21.60 Rg=12.843 SPS=15693
bl=0 pos[1]=[13.5 -7.8 3.4] dr=1.25 t=0.0ps kin=1.44 pot=21.59 Rg=12.788 SPS=14657
bl=0 pos[1]=[13.6 -7.7 4.0] dr=1.28 t=0.0ps kin=1.55 pot=21.58 Rg=12.710 SPS=15304
bl=0 pos[1]=[12.0 -6.8 4.8] dr=1.31 t=0.0ps kin=1.53 pot=21.59 Rg=12.667 SPS=14126
bl=0 pos[1]=[11.2 -4.6 3.3] dr=1.30 t=0.0ps kin=1.53 pot=21.62 Rg=12.625 SPS=14725
bl=0 pos[1]=[10.0 -2.3 0.8] dr=1.28 t=0.0ps kin=1.52 pot=21.61 Rg=12.537 SPS=14998
bl=0 pos[1]=[10.1 -1.0 0.0] dr=1.24 t=0.0ps kin=1.49 pot=21.61 Rg=12.500 SPS=14985
bl=0 pos[1]=[9.8 -0.2 -0.1] dr=1.23 t=0.0ps kin=1.51 pot=21.62 Rg=12.467 SPS=15193
bl=0 pos[1]=[8.9 -0.9 0.8] dr=1.28 t=0.0ps kin=1.52 pot=21.60 Rg=12.458 SPS=15583
bl=0 pos[1]=[8.7 -1.6 1.2] dr=1.26 t=0.0ps kin=1.49 pot=21.64 Rg=12.488 SPS=15535
bl=0 pos[1]=[8.2 -2.2 1.7] dr=1.27 t=0.0ps kin=1.49 pot=21.63 Rg=12.516 SPS=15794
bl=0 pos[1]=[8.8 -3.2 1.5] dr=1.27 t=0.0ps kin=1.49 pot=21.61 Rg=12.568 SPS=14251
bl=0 pos[1]=[9.5 -4.4 1.4] dr=1.25 t=0.0ps kin=1.49 pot=21.62 Rg=12.583 SPS=15538
bl=0 pos[1]=[9.6 -3.4 2.9] dr=1.29 t=0.0ps kin=1.52 pot=21.58 Rg=12.574 SPS=15765
bl=0 pos[1]=[9.8 -2.2 3.3] dr=1.26 t=0.0ps kin=1.55 pot=21.60 Rg=12.525 SPS=14635
bl=0 pos[1]=[11.3 -1.9 3.2] dr=1.27 t=0.0ps kin=1.50 pot=21.65 Rg=12.525 SPS=15586
bl=0 pos[1]=[12.3 -1.0 3.6] dr=1.27 t=0.0ps kin=1.53 pot=21.63 Rg=12.543 SPS=15903
bl=0 pos[1]=[13.3 0.2 4.5] dr=1.26 t=0.0ps kin=1.51 pot=21.64 Rg=12.557 SPS=14694
bl=0 pos[1]=[15.1 -0.3 3.5] dr=1.23 t=0.0ps kin=1.50 pot=21.64 Rg=12.550 SPS=15851
bl=0 pos[1]=[15.3 -1.0 4.7] dr=1.24 t=0.0ps kin=1.49 pot=21.59 Rg=12.550 SPS=15151
bl=0 pos[1]=[15.5 -0.2 4.7] dr=1.27 t=0.0ps kin=1.50 pot=21.60 Rg=12.526 SPS=14689
bl=0 pos[1]=[17.1 -0.6 3.6] dr=1.26 t=0.0ps kin=1.50 pot=21.65 Rg=12.487 SPS=15510
bl=0 pos[1]=[17.7 -0.3 2.6] dr=1.26 t=0.0ps kin=1.53 pot=21.58 Rg=12.496 SPS=14952
bl=0 pos[1]=[19.3 -0.2 3.1] dr=1.26 t=0.0ps kin=1.55 pot=21.62 Rg=12.544 SPS=15678
bl=0 pos[1]=[20.2 -1.5 4.5] dr=1.30 t=0.0ps kin=1.53 pot=21.60 Rg=12.615 SPS=15646
bl=0 pos[1]=[19.3 -3.0 4.8] dr=1.25 t=0.0ps kin=1.48 pot=21.63 Rg=12.657 SPS=14483
bl=0 pos[1]=[18.8 -3.5 6.4] dr=1.23 t=0.0ps kin=1.46 pot=21.61 Rg=12.653 SPS=15484
bl=0 pos[1]=[17.9 -4.0 7.2] dr=1.25 t=0.0ps kin=1.47 pot=21.62 Rg=12.643 SPS=15727
bl=0 pos[1]=[16.7 -5.1 7.0] dr=1.26 t=0.0ps kin=1.51 pot=21.60 Rg=12.612 SPS=15453
bl=0 pos[1]=[16.7 -4.6 7.1] dr=1.24 t=0.0ps kin=1.43 pot=21.61 Rg=12.543 SPS=15341
bl=0 pos[1]=[16.9 -4.8 8.0] dr=1.23 t=0.0ps kin=1.48 pot=21.58 Rg=12.479 SPS=15349
bl=0 pos[1]=[16.0 -4.9 9.2] dr=1.25 t=0.0ps kin=1.48 pot=21.59 Rg=12.412 SPS=15213
bl=0 pos[1]=[14.3 -5.5 10.1] dr=1.25 t=0.0ps kin=1.49 pot=21.59 Rg=12.357 SPS=15405
bl=0 pos[1]=[13.5 -5.7 9.7] dr=1.25 t=0.0ps kin=1.50 pot=21.58 Rg=12.295 SPS=15740
bl=0 pos[1]=[13.7 -5.6 10.6] dr=1.25 t=0.0ps kin=1.47 pot=21.57 Rg=12.219 SPS=15298
bl=0 pos[1]=[13.1 -4.0 10.2] dr=1.25 t=0.0ps kin=1.46 pot=21.62 Rg=12.166 SPS=15613
bl=0 pos[1]=[11.7 -2.8 12.3] dr=1.28 t=0.0ps kin=1.50 pot=21.60 Rg=12.104 SPS=15885
bl=0 pos[1]=[11.6 -2.6 12.2] dr=1.27 t=0.0ps kin=1.49 pot=21.60 Rg=12.051 SPS=14927
bl=0 pos[1]=[11.7 -3.1 12.0] dr=1.27 t=0.0ps kin=1.51 pot=21.62 Rg=12.054 SPS=15853
bl=0 pos[1]=[11.7 -3.4 9.4] dr=1.27 t=0.0ps kin=1.48 pot=21.64 Rg=12.082 SPS=14691
bl=0 pos[1]=[12.6 -3.5 10.1] dr=1.27 t=0.0ps kin=1.46 pot=21.58 Rg=12.105 SPS=15485
bl=0 pos[1]=[14.8 -3.0 12.1] dr=1.29 t=0.0ps kin=1.46 pot=21.62 Rg=12.081 SPS=15465
bl=0 pos[1]=[15.5 -3.2 13.2] dr=1.26 t=0.0ps kin=1.48 pot=21.57 Rg=12.013 SPS=14774
bl=0 pos[1]=[16.1 -3.7 12.6] dr=1.27 t=0.0ps kin=1.50 pot=21.60 Rg=11.947 SPS=15412
bl=0 pos[1]=[16.1 -2.8 11.8] dr=1.28 t=0.0ps kin=1.54 pot=21.56 Rg=11.882 SPS=15319
bl=0 pos[1]=[17.6 -2.6 12.5] dr=1.27 t=0.0ps kin=1.46 pot=21.59 Rg=11.782 SPS=14661
bl=0 pos[1]=[17.3 -1.3 11.4] dr=1.26 t=0.0ps kin=1.56 pot=21.57 Rg=11.716 SPS=15855
bl=0 pos[1]=[17.0 -2.1 10.9] dr=1.29 t=0.0ps kin=1.56 pot=21.58 Rg=11.703 SPS=15822
bl=0 pos[1]=[16.1 -3.4 10.3] dr=1.27 t=0.0ps kin=1.52 pot=21.63 Rg=11.698 SPS=15103
bl=0 pos[1]=[14.1 -3.1 9.9] dr=1.25 t=0.0ps kin=1.47 pot=21.63 Rg=11.687 SPS=16030
bl=0 pos[1]=[12.9 -1.7 11.0] dr=1.23 t=0.0ps kin=1.48 pot=21.59 Rg=11.679 SPS=15538
bl=0 pos[1]=[12.9 -1.6 11.4] dr=1.23 t=0.0ps kin=1.48 pot=21.60 Rg=11.692 SPS=14640
bl=0 pos[1]=[13.0 -1.2 11.9] dr=1.27 t=0.0ps kin=1.51 pot=21.60 Rg=11.712 SPS=15223
bl=0 pos[1]=[13.9 0.8 13.1] dr=1.28 t=0.0ps kin=1.50 pot=21.54 Rg=11.732 SPS=15002
bl=0 pos[1]=[13.6 3.4 12.0] dr=1.25 t=0.0ps kin=1.52 pot=21.61 Rg=11.722 SPS=15281
bl=0 pos[1]=[15.0 5.1 10.7] dr=1.28 t=0.0ps kin=1.48 pot=21.62 Rg=11.697 SPS=15488
bl=0 pos[1]=[15.7 5.7 10.3] dr=1.25 t=0.0ps kin=1.55 pot=21.57 Rg=11.701 SPS=15135
bl=0 pos[1]=[15.6 5.2 8.9] dr=1.30 t=0.0ps kin=1.52 pot=21.60 Rg=11.727 SPS=15764
bl=0 pos[1]=[15.0 6.4 7.1] dr=1.28 t=0.0ps kin=1.51 pot=21.59 Rg=11.743 SPS=15725
bl=0 pos[1]=[13.3 5.4 7.0] dr=1.25 t=0.0ps kin=1.49 pot=21.62 Rg=11.696 SPS=14909
bl=0 pos[1]=[12.4 4.5 6.4] dr=1.30 t=0.0ps kin=1.49 pot=21.64 Rg=11.649 SPS=15487
bl=0 pos[1]=[13.3 4.1 6.5] dr=1.26 t=0.0ps kin=1.55 pot=21.62 Rg=11.658 SPS=15225
bl=0 pos[1]=[13.7 1.9 6.5] dr=1.27 t=0.0ps kin=1.56 pot=21.57 Rg=11.669 SPS=14727
bl=0 pos[1]=[13.4 3.0 8.3] dr=1.29 t=0.0ps kin=1.52 pot=21.63 Rg=11.682 SPS=15424
bl=0 pos[1]=[12.9 3.4 7.3] dr=1.27 t=0.0ps kin=1.46 pot=21.60 Rg=11.671 SPS=15406
bl=0 pos[1]=[13.3 2.7 6.8] dr=1.23 t=0.0ps kin=1.51 pot=21.60 Rg=11.657 SPS=15134
bl=0 pos[1]=[13.1 2.2 6.2] dr=1.26 t=0.0ps kin=1.44 pot=21.57 Rg=11.658 SPS=15449
bl=0 pos[1]=[13.6 3.6 6.1] dr=1.26 t=0.0ps kin=1.49 pot=21.59 Rg=11.668 SPS=14817
bl=0 pos[1]=[13.2 4.2 5.8] dr=1.28 t=0.0ps kin=1.55 pot=21.60 Rg=11.695 SPS=15689
bl=0 pos[1]=[11.9 4.5 5.6] dr=1.30 t=0.0ps kin=1.61 pot=21.63 Rg=11.666 SPS=15229
bl=0 pos[1]=[10.6 5.4 5.4] dr=1.31 t=0.0ps kin=1.51 pot=21.60 Rg=11.658 SPS=15060
bl=0 pos[1]=[10.0 3.6 4.8] dr=1.26 t=0.0ps kin=1.52 pot=21.56 Rg=11.637 SPS=15574
bl=0 pos[1]=[9.1 2.6 3.9] dr=1.25 t=0.0ps kin=1.45 pot=21.61 Rg=11.582 SPS=15346
bl=0 pos[1]=[9.8 4.2 3.4] dr=1.25 t=0.0ps kin=1.50 pot=21.58 Rg=11.529 SPS=14922
bl=0 pos[1]=[10.8 4.3 2.7] dr=1.26 t=0.0ps kin=1.53 pot=21.57 Rg=11.521 SPS=15497
bl=0 pos[1]=[12.3 3.7 1.4] dr=1.26 t=0.0ps kin=1.48 pot=21.59 Rg=11.569 SPS=14988
bl=0 pos[1]=[12.3 4.0 0.7] dr=1.27 t=0.0ps kin=1.46 pot=21.59 Rg=11.590 SPS=15306
bl=0 pos[1]=[12.3 5.6 1.2] dr=1.26 t=0.0ps kin=1.48 pot=21.62 Rg=11.660 SPS=14902
bl=0 pos[1]=[11.6 6.7 0.9] dr=1.26 t=0.0ps kin=1.53 pot=21.63 Rg=11.731 SPS=13760
bl=0 pos[1]=[12.2 6.7 0.1] dr=1.26 t=0.0ps kin=1.48 pot=21.62 Rg=11.757 SPS=14693
bl=0 pos[1]=[10.9 7.7 0.8] dr=1.27 t=0.0ps kin=1.51 pot=21.56 Rg=11.752 SPS=14069
bl=0 pos[1]=[8.9 8.3 0.6] dr=1.26 t=0.0ps kin=1.45 pot=21.62 Rg=11.764 SPS=15348
bl=0 pos[1]=[8.1 6.2 1.9] dr=1.25 t=0.0ps kin=1.54 pot=21.53 Rg=11.751 SPS=13892
bl=0 pos[1]=[8.6 6.3 2.6] dr=1.32 t=0.0ps kin=1.55 pot=21.59 Rg=11.687 SPS=15426
bl=0 pos[1]=[6.3 5.8 3.0] dr=1.30 t=0.0ps kin=1.55 pot=21.60 Rg=11.611 SPS=14976
bl=0 pos[1]=[5.8 5.6 4.3] dr=1.29 t=0.0ps kin=1.51 pot=21.61 Rg=11.547 SPS=14934
bl=0 pos[1]=[4.3 7.6 4.5] dr=1.23 t=0.0ps kin=1.47 pot=21.62 Rg=11.528 SPS=15963
bl=0 pos[1]=[5.0 8.9 3.1] dr=1.25 t=0.0ps kin=1.52 pot=21.56 Rg=11.524 SPS=15907
bl=0 pos[1]=[5.5 10.2 2.4] dr=1.29 t=0.0ps kin=1.50 pot=21.59 Rg=11.519 SPS=15570
bl=0 pos[1]=[8.3 9.0 3.5] dr=1.27 t=0.0ps kin=1.50 pot=21.58 Rg=11.496 SPS=15071
bl=0 pos[1]=[8.6 7.2 3.7] dr=1.24 t=0.0ps kin=1.45 pot=21.60 Rg=11.524 SPS=15630
bl=0 pos[1]=[8.4 7.2 4.8] dr=1.24 t=0.0ps kin=1.53 pot=21.58 Rg=11.551 SPS=14633
bl=0 pos[1]=[7.9 8.7 5.0] dr=1.29 t=0.0ps kin=1.51 pot=21.61 Rg=11.574 SPS=15573
bl=0 pos[1]=[8.0 9.1 5.0] dr=1.23 t=0.0ps kin=1.49 pot=21.62 Rg=11.603 SPS=14067
bl=0 pos[1]=[8.4 9.1 4.9] dr=1.25 t=0.0ps kin=1.55 pot=21.61 Rg=11.556 SPS=14884
bl=0 pos[1]=[7.2 8.9 4.0] dr=1.27 t=0.0ps kin=1.55 pot=21.58 Rg=11.530 SPS=15804
bl=0 pos[1]=[7.9 9.4 4.8] dr=1.28 t=0.0ps kin=1.50 pot=21.63 Rg=11.508 SPS=14850
bl=0 pos[1]=[7.4 9.3 6.3] dr=1.28 t=0.0ps kin=1.54 pot=21.62 Rg=11.462 SPS=15185
bl=0 pos[1]=[7.7 9.9 6.9] dr=1.26 t=0.0ps kin=1.50 pot=21.55 Rg=11.455 SPS=15104
bl=0 pos[1]=[7.1 10.0 7.0] dr=1.23 t=0.0ps kin=1.48 pot=21.53 Rg=11.477 SPS=15178
bl=0 pos[1]=[8.0 10.6 7.4] dr=1.21 t=0.0ps kin=1.45 pot=21.58 Rg=11.493 SPS=14900
bl=0 pos[1]=[10.0 8.8 5.8] dr=1.24 t=0.0ps kin=1.49 pot=21.52 Rg=11.521 SPS=14086
bl=0 pos[1]=[8.8 7.7 3.3] dr=1.22 t=0.0ps kin=1.51 pot=21.60 Rg=11.447 SPS=16072
bl=0 pos[1]=[7.7 6.9 3.9] dr=1.24 t=0.0ps kin=1.47 pot=21.57 Rg=11.406 SPS=15798
bl=0 pos[1]=[8.2 8.7 4.4] dr=1.23 t=0.0ps kin=1.44 pot=21.61 Rg=11.421 SPS=14803
bl=0 pos[1]=[10.3 9.9 4.0] dr=1.24 t=0.0ps kin=1.48 pot=21.61 Rg=11.436 SPS=16213
bl=0 pos[1]=[11.0 9.8 3.7] dr=1.24 t=0.0ps kin=1.50 pot=21.56 Rg=11.466 SPS=16050
bl=0 pos[1]=[11.4 9.3 4.0] dr=1.21 t=0.0ps kin=1.46 pot=21.59 Rg=11.544 SPS=15034
bl=0 pos[1]=[13.0 9.9 4.4] dr=1.25 t=0.0ps kin=1.45 pot=21.57 Rg=11.574 SPS=15253
bl=0 pos[1]=[12.6 8.4 5.2] dr=1.26 t=0.0ps kin=1.47 pot=21.61 Rg=11.581 SPS=14861
bl=0 pos[1]=[12.1 7.4 4.3] dr=1.27 t=0.0ps kin=1.50 pot=21.58 Rg=11.566 SPS=15985
bl=0 pos[1]=[11.2 5.6 4.6] dr=1.27 t=0.0ps kin=1.49 pot=21.57 Rg=11.522 SPS=15361
bl=0 pos[1]=[11.7 6.8 3.1] dr=1.26 t=0.0ps kin=1.40 pot=21.56 Rg=11.464 SPS=14883
bl=0 pos[1]=[11.5 6.4 3.6] dr=1.29 t=0.0ps kin=1.47 pot=21.54 Rg=11.401 SPS=15308
bl=0 pos[1]=[10.4 6.8 3.0] dr=1.27 t=0.0ps kin=1.51 pot=21.55 Rg=11.339 SPS=15368
bl=0 pos[1]=[8.7 5.6 2.4] dr=1.27 t=0.0ps kin=1.47 pot=21.56 Rg=11.304 SPS=15015
bl=0 pos[1]=[7.0 6.6 2.8] dr=1.25 t=0.0ps kin=1.47 pot=21.57 Rg=11.287 SPS=15242
bl=0 pos[1]=[4.6 5.8 1.4] dr=1.25 t=0.0ps kin=1.42 pot=21.59 Rg=11.275 SPS=12877
bl=0 pos[1]=[3.6 3.5 2.2] dr=1.24 t=0.0ps kin=1.46 pot=21.55 Rg=11.277 SPS=15548
bl=0 pos[1]=[4.1 2.7 2.5] dr=1.21 t=0.0ps kin=1.45 pot=21.56 Rg=11.260 SPS=15089
bl=0 pos[1]=[4.1 2.5 1.7] dr=1.22 t=0.0ps kin=1.48 pot=21.55 Rg=11.280 SPS=15123
bl=0 pos[1]=[3.6 2.8 3.6] dr=1.25 t=0.0ps kin=1.45 pot=21.56 Rg=11.303 SPS=15307
bl=0 pos[1]=[2.1 3.5 2.8] dr=1.25 t=0.0ps kin=1.45 pot=21.59 Rg=11.310 SPS=15295
bl=0 pos[1]=[1.6 3.1 1.9] dr=1.23 t=0.0ps kin=1.48 pot=21.60 Rg=11.319 SPS=15658
bl=0 pos[1]=[0.1 3.9 2.8] dr=1.25 t=0.0ps kin=1.54 pot=21.59 Rg=11.377 SPS=15241
bl=0 pos[1]=[-0.4 4.4 4.5] dr=1.28 t=0.0ps kin=1.57 pot=21.62 Rg=11.431 SPS=14641
bl=0 pos[1]=[-1.5 4.0 4.2] dr=1.26 t=0.0ps kin=1.49 pot=21.65 Rg=11.480 SPS=15910
bl=0 pos[1]=[-2.0 3.2 4.6] dr=1.27 t=0.0ps kin=1.49 pot=21.59 Rg=11.509 SPS=15875
bl=0 pos[1]=[-0.8 3.3 4.5] dr=1.25 t=0.0ps kin=1.55 pot=21.59 Rg=11.535 SPS=15391
bl=0 pos[1]=[0.9 3.1 5.4] dr=1.27 t=0.0ps kin=1.51 pot=21.61 Rg=11.544 SPS=15960
bl=0 pos[1]=[2.1 3.7 4.8] dr=1.27 t=0.0ps kin=1.50 pot=21.61 Rg=11.526 SPS=14958
bl=0 pos[1]=[1.6 1.2 5.1] dr=1.28 t=0.0ps kin=1.49 pot=21.57 Rg=11.512 SPS=14556
bl=0 pos[1]=[0.0 0.8 5.3] dr=1.26 t=0.0ps kin=1.50 pot=21.58 Rg=11.501 SPS=15425
bl=0 pos[1]=[-0.4 0.5 5.2] dr=1.25 t=0.0ps kin=1.51 pot=21.62 Rg=11.532 SPS=14282
bl=0 pos[1]=[-0.3 -0.2 6.0] dr=1.25 t=0.0ps kin=1.54 pot=21.59 Rg=11.567 SPS=15298
bl=0 pos[1]=[-0.6 -0.1 5.0] dr=1.28 t=0.0ps kin=1.53 pot=21.62 Rg=11.582 SPS=15590
bl=0 pos[1]=[-0.4 -0.2 4.8] dr=1.29 t=0.0ps kin=1.50 pot=21.65 Rg=11.576 SPS=15143
bl=0 pos[1]=[-0.2 -0.9 5.5] dr=1.29 t=0.0ps kin=1.54 pot=21.60 Rg=11.549 SPS=15769
bl=0 pos[1]=[0.2 -0.7 4.5] dr=1.26 t=0.0ps kin=1.53 pot=21.59 Rg=11.504 SPS=15588
bl=0 pos[1]=[0.7 -0.4 4.0] dr=1.26 t=0.0ps kin=1.52 pot=21.57 Rg=11.499 SPS=14497
bl=0 pos[1]=[0.9 1.6 4.1] dr=1.26 t=0.0ps kin=1.49 pot=21.56 Rg=11.502 SPS=16029
bl=0 pos[1]=[1.9 3.2 5.4] dr=1.28 t=0.0ps kin=1.49 pot=21.59 Rg=11.498 SPS=15561
bl=0 pos[1]=[1.3 2.4 5.8] dr=1.25 t=0.0ps kin=1.48 pot=21.61 Rg=11.487 SPS=15137
bl=0 pos[1]=[3.4 2.2 5.7] dr=1.26 t=0.0ps kin=1.47 pot=21.57 Rg=11.499 SPS=15750
bl=0 pos[1]=[4.3 2.5 6.0] dr=1.26 t=0.0ps kin=1.57 pot=21.54 Rg=11.502 SPS=15969
bl=0 pos[1]=[3.2 3.6 5.7] dr=1.27 t=0.0ps kin=1.54 pot=21.58 Rg=11.456 SPS=14584
bl=0 pos[1]=[2.9 3.7 4.8] dr=1.31 t=0.0ps kin=1.52 pot=21.59 Rg=11.421 SPS=15340
bl=0 pos[1]=[3.4 3.8 6.2] dr=1.28 t=0.0ps kin=1.53 pot=21.58 Rg=11.392 SPS=15004
bl=0 pos[1]=[5.7 4.8 7.8] dr=1.26 t=0.0ps kin=1.48 pot=21.58 Rg=11.341 SPS=15638
bl=0 pos[1]=[6.3 7.6 7.8] dr=1.28 t=0.0ps kin=1.53 pot=21.57 Rg=11.351 SPS=15443
bl=0 pos[1]=[5.3 9.9 6.7] dr=1.25 t=0.0ps kin=1.51 pot=21.59 Rg=11.356 SPS=14598
bl=0 pos[1]=[6.6 9.1 8.0] dr=1.28 t=0.0ps kin=1.54 pot=21.60 Rg=11.377 SPS=15460
bl=0 pos[1]=[7.5 7.8 6.3] dr=1.28 t=0.0ps kin=1.53 pot=21.61 Rg=11.405 SPS=15348
bl=0 pos[1]=[6.3 7.6 5.8] dr=1.26 t=0.0ps kin=1.52 pot=21.60 Rg=11.384 SPS=14498
bl=0 pos[1]=[5.3 8.1 6.3] dr=1.24 t=0.0ps kin=1.57 pot=21.59 Rg=11.346 SPS=15382
bl=0 pos[1]=[6.0 8.5 6.7] dr=1.28 t=0.0ps kin=1.53 pot=21.59 Rg=11.366 SPS=14421
bl=0 pos[1]=[8.4 5.7 7.3] dr=1.30 t=0.0ps kin=1.54 pot=21.59 Rg=11.415 SPS=15602
bl=0 pos[1]=[9.3 6.4 5.1] dr=1.30 t=0.0ps kin=1.55 pot=21.57 Rg=11.450 SPS=15288
bl=0 pos[1]=[9.5 7.0 4.3] dr=1.26 t=0.0ps kin=1.51 pot=21.63 Rg=11.464 SPS=13880
bl=0 pos[1]=[9.8 6.0 3.0] dr=1.28 t=0.0ps kin=1.60 pot=21.57 Rg=11.445 SPS=15141
bl=0 pos[1]=[9.3 6.8 1.5] dr=1.27 t=0.0ps kin=1.47 pot=21.59 Rg=11.437 SPS=14515
bl=0 pos[1]=[8.7 5.4 0.7] dr=1.25 t=0.0ps kin=1.58 pot=21.58 Rg=11.436 SPS=15241
bl=0 pos[1]=[8.0 6.8 1.1] dr=1.29 t=0.0ps kin=1.53 pot=21.58 Rg=11.440 SPS=15042
bl=0 pos[1]=[8.5 6.3 0.8] dr=1.26 t=0.0ps kin=1.50 pot=21.55 Rg=11.446 SPS=14376
bl=0 pos[1]=[9.1 6.4 -0.9] dr=1.25 t=0.0ps kin=1.41 pot=21.55 Rg=11.507 SPS=15507
bl=0 pos[1]=[6.6 4.6 -0.6] dr=1.24 t=0.0ps kin=1.44 pot=21.53 Rg=11.574 SPS=13916
bl=0 pos[1]=[5.4 4.6 -0.8] dr=1.25 t=0.0ps kin=1.44 pot=21.57 Rg=11.607 SPS=15513
bl=0 pos[1]=[4.9 4.7 -0.6] dr=1.27 t=0.0ps kin=1.49 pot=21.55 Rg=11.622 SPS=15355
bl=0 pos[1]=[5.0 4.8 -0.1] dr=1.27 t=0.0ps kin=1.52 pot=21.54 Rg=11.591 SPS=14799
bl=0 pos[1]=[4.9 4.7 1.1] dr=1.25 t=0.0ps kin=1.46 pot=21.59 Rg=11.495 SPS=15607
bl=0 pos[1]=[5.0 5.6 0.9] dr=1.26 t=0.0ps kin=1.46 pot=21.63 Rg=11.421 SPS=15649
bl=0 pos[1]=[5.7 5.7 -0.0] dr=1.23 t=0.0ps kin=1.49 pot=21.60 Rg=11.419 SPS=14767
bl=0 pos[1]=[6.7 4.5 1.6] dr=1.25 t=0.0ps kin=1.51 pot=21.65 Rg=11.421 SPS=15406
bl=0 pos[1]=[8.2 4.3 1.6] dr=1.26 t=0.0ps kin=1.51 pot=21.62 Rg=11.403 SPS=14683
bl=0 pos[1]=[7.8 4.7 0.6] dr=1.24 t=0.0ps kin=1.53 pot=21.60 Rg=11.366 SPS=14926
bl=0 pos[1]=[7.8 6.2 -0.5] dr=1.27 t=0.0ps kin=1.52 pot=21.64 Rg=11.371 SPS=15150
bl=0 pos[1]=[7.2 6.8 -0.8] dr=1.28 t=0.0ps kin=1.50 pot=21.59 Rg=11.383 SPS=14266
bl=0 pos[1]=[6.7 7.1 -0.8] dr=1.24 t=0.0ps kin=1.52 pot=21.54 Rg=11.416 SPS=15577
bl=0 pos[1]=[8.1 6.0 -0.5] dr=1.21 t=0.0ps kin=1.47 pot=21.62 Rg=11.391 SPS=14314
bl=0 pos[1]=[8.5 6.0 0.9] dr=1.28 t=0.0ps kin=1.49 pot=21.58 Rg=11.401 SPS=15167
bl=0 pos[1]=[8.8 6.2 -0.3] dr=1.24 t=0.0ps kin=1.44 pot=21.56 Rg=11.427 SPS=15749
bl=0 pos[1]=[11.1 7.2 2.1] dr=1.23 t=0.0ps kin=1.50 pot=21.54 Rg=11.467 SPS=15318
bl=0 pos[1]=[10.1 7.1 2.2] dr=1.25 t=0.0ps kin=1.47 pot=21.61 Rg=11.549 SPS=15424
bl=0 pos[1]=[11.7 7.4 2.8] dr=1.30 t=0.0ps kin=1.52 pot=21.60 Rg=11.650 SPS=15118
bl=0 pos[1]=[13.9 9.1 2.5] dr=1.28 t=0.0ps kin=1.56 pot=21.60 Rg=11.735 SPS=14509
bl=0 pos[1]=[12.4 8.2 0.8] dr=1.26 t=0.0ps kin=1.52 pot=21.60 Rg=11.769 SPS=15139
bl=0 pos[1]=[11.0 7.8 1.0] dr=1.27 t=0.0ps kin=1.54 pot=21.62 Rg=11.790 SPS=14834
bl=0 pos[1]=[11.1 8.1 1.2] dr=1.26 t=0.0ps kin=1.48 pot=21.58 Rg=11.816 SPS=15916
bl=0 pos[1]=[11.2 8.4 1.8] dr=1.25 t=0.0ps kin=1.52 pot=21.59 Rg=11.844 SPS=15685
bl=0 pos[1]=[10.8 9.6 2.5] dr=1.28 t=0.0ps kin=1.54 pot=21.59 Rg=11.891 SPS=14799
bl=0 pos[1]=[10.5 11.2 1.9] dr=1.25 t=0.0ps kin=1.48 pot=21.61 Rg=11.924 SPS=16088
bl=0 pos[1]=[9.7 9.3 1.3] dr=1.24 t=0.0ps kin=1.46 pot=21.59 Rg=11.955 SPS=15311
bl=0 pos[1]=[9.9 10.0 0.0] dr=1.26 t=0.0ps kin=1.51 pot=21.57 Rg=11.974 SPS=15034
bl=0 pos[1]=[10.9 9.9 -1.3] dr=1.28 t=0.0ps kin=1.48 pot=21.58 Rg=11.985 SPS=15112
bl=0 pos[1]=[10.8 8.9 -0.7] dr=1.27 t=0.0ps kin=1.46 pot=21.55 Rg=11.957 SPS=15244
bl=0 pos[1]=[10.5 7.3 1.5] dr=1.26 t=0.0ps kin=1.48 pot=21.57 Rg=11.921 SPS=15311
bl=0 pos[1]=[10.8 7.2 3.2] dr=1.25 t=0.0ps kin=1.50 pot=21.55 Rg=11.897 SPS=15243
bl=0 pos[1]=[9.7 7.0 3.7] dr=1.27 t=0.0ps kin=1.46 pot=21.60 Rg=11.872 SPS=14756
bl=0 pos[1]=[8.2 5.2 3.1] dr=1.25 t=0.0ps kin=1.51 pot=21.56 Rg=11.863 SPS=16067
bl=0 pos[1]=[8.6 5.9 2.5] dr=1.23 t=0.0ps kin=1.52 pot=21.58 Rg=11.883 SPS=15459
bl=0 pos[1]=[8.3 7.8 1.3] dr=1.24 t=0.0ps kin=1.47 pot=21.56 Rg=11.897 SPS=15477
bl=0 pos[1]=[7.6 10.2 1.1] dr=1.26 t=0.0ps kin=1.50 pot=21.56 Rg=11.899 SPS=15593
bl=0 pos[1]=[7.2 11.3 -1.1] dr=1.25 t=0.0ps kin=1.53 pot=21.57 Rg=11.921 SPS=15478
bl=0 pos[1]=[7.1 11.7 -2.2] dr=1.26 t=0.0ps kin=1.46 pot=21.60 Rg=11.905 SPS=15304
bl=0 pos[1]=[6.5 12.9 -1.4] dr=1.28 t=0.0ps kin=1.50 pot=21.56 Rg=11.928 SPS=15518
bl=0 pos[1]=[6.5 12.6 1.4] dr=1.29 t=0.0ps kin=1.50 pot=21.56 Rg=11.940 SPS=14940
bl=0 pos[1]=[4.9 11.8 1.3] dr=1.26 t=0.0ps kin=1.46 pot=21.58 Rg=11.929 SPS=14451
bl=0 pos[1]=[4.2 12.6 2.8] dr=1.23 t=0.0ps kin=1.47 pot=21.57 Rg=11.928 SPS=15003
bl=0 pos[1]=[4.9 14.6 3.1] dr=1.23 t=0.0ps kin=1.47 pot=21.61 Rg=11.941 SPS=14478
bl=0 pos[1]=[4.8 14.3 2.7] dr=1.24 t=0.0ps kin=1.48 pot=21.54 Rg=11.949 SPS=14911
bl=0 pos[1]=[3.4 14.2 1.6] dr=1.21 t=0.0ps kin=1.46 pot=21.59 Rg=11.951 SPS=14889
bl=0 pos[1]=[1.3 14.7 2.8] dr=1.23 t=0.0ps kin=1.49 pot=21.59 Rg=11.904 SPS=14624
bl=0 pos[1]=[1.9 15.2 4.3] dr=1.27 t=0.0ps kin=1.52 pot=21.62 Rg=11.927 SPS=15737
bl=0 pos[1]=[2.1 14.3 6.6] dr=1.28 t=0.0ps kin=1.52 pot=21.61 Rg=11.933 SPS=14514
bl=0 pos[1]=[3.6 15.4 7.1] dr=1.27 t=0.0ps kin=1.47 pot=21.61 Rg=11.940 SPS=14650
bl=0 pos[1]=[5.0 17.9 6.5] dr=1.27 t=0.0ps kin=1.49 pot=21.59 Rg=11.911 SPS=16222
bl=0 pos[1]=[5.4 18.5 4.6] dr=1.24 t=0.0ps kin=1.54 pot=21.53 Rg=11.840 SPS=14735
bl=0 pos[1]=[5.9 19.0 2.6] dr=1.25 t=0.0ps kin=1.49 pot=21.56 Rg=11.729 SPS=15246
bl=0 pos[1]=[5.8 17.7 4.1] dr=1.25 t=0.0ps kin=1.47 pot=21.57 Rg=11.650 SPS=15573
bl=0 pos[1]=[7.0 18.3 4.9] dr=1.26 t=0.0ps kin=1.51 pot=21.53 Rg=11.623 SPS=14427
bl=0 pos[1]=[7.0 18.3 3.2] dr=1.26 t=0.0ps kin=1.47 pot=21.59 Rg=11.616 SPS=15326
bl=0 pos[1]=[4.7 18.3 3.0] dr=1.22 t=0.0ps kin=1.48 pot=21.59 Rg=11.624 SPS=14654
bl=0 pos[1]=[3.0 18.9 1.5] dr=1.25 t=0.0ps kin=1.53 pot=21.60 Rg=11.624 SPS=15614
bl=0 pos[1]=[3.3 19.1 1.0] dr=1.26 t=0.0ps kin=1.51 pot=21.57 Rg=11.643 SPS=16030
bl=0 pos[1]=[3.7 18.2 0.3] dr=1.27 t=0.0ps kin=1.52 pot=21.59 Rg=11.633 SPS=14643
bl=0 pos[1]=[4.2 17.3 0.6] dr=1.24 t=0.0ps kin=1.44 pot=21.62 Rg=11.594 SPS=15126
bl=0 pos[1]=[3.2 17.5 1.5] dr=1.22 t=0.0ps kin=1.45 pot=21.57 Rg=11.590 SPS=13579
bl=0 pos[1]=[2.3 16.5 1.5] dr=1.23 t=0.0ps kin=1.46 pot=21.58 Rg=11.643 SPS=14483
bl=0 pos[1]=[2.3 16.2 2.3] dr=1.21 t=0.0ps kin=1.48 pot=21.61 Rg=11.668 SPS=14565
bl=0 pos[1]=[0.9 15.4 2.4] dr=1.27 t=0.0ps kin=1.49 pot=21.56 Rg=11.688 SPS=15540
bl=0 pos[1]=[-0.9 14.0 2.6] dr=1.27 t=0.0ps kin=1.47 pot=21.62 Rg=11.678 SPS=15674
bl=0 pos[1]=[-1.3 12.1 3.2] dr=1.24 t=0.0ps kin=1.50 pot=21.60 Rg=11.627 SPS=14832
bl=0 pos[1]=[-1.5 11.1 2.7] dr=1.21 t=0.0ps kin=1.50 pot=21.60 Rg=11.632 SPS=15716
bl=0 pos[1]=[-1.4 11.8 4.0] dr=1.25 t=0.0ps kin=1.53 pot=21.56 Rg=11.715 SPS=16122
bl=0 pos[1]=[0.3 12.8 4.5] dr=1.25 t=0.0ps kin=1.57 pot=21.54 Rg=11.768 SPS=14995
bl=0 pos[1]=[2.5 13.6 4.2] dr=1.30 t=0.0ps kin=1.51 pot=21.57 Rg=11.756 SPS=13949
bl=0 pos[1]=[2.8 13.5 4.5] dr=1.26 t=0.0ps kin=1.49 pot=21.60 Rg=11.756 SPS=14822
bl=0 pos[1]=[2.5 16.1 4.5] dr=1.30 t=0.0ps kin=1.52 pot=21.56 Rg=11.768 SPS=15616
bl=0 pos[1]=[2.8 17.0 4.0] dr=1.26 t=0.0ps kin=1.49 pot=21.57 Rg=11.789 SPS=15202
bl=0 pos[1]=[5.0 15.6 5.8] dr=1.26 t=0.0ps kin=1.49 pot=21.62 Rg=11.818 SPS=15381
bl=0 pos[1]=[6.6 14.4 7.4] dr=1.30 t=0.0ps kin=1.57 pot=21.58 Rg=11.791 SPS=15474
bl=0 pos[1]=[6.1 14.6 9.0] dr=1.26 t=0.0ps kin=1.50 pot=21.58 Rg=11.746 SPS=15067
bl=0 pos[1]=[5.7 13.3 9.1] dr=1.23 t=0.0ps kin=1.48 pot=21.56 Rg=11.727 SPS=14807
bl=0 pos[1]=[5.2 11.7 8.8] dr=1.26 t=0.0ps kin=1.48 pot=21.56 Rg=11.721 SPS=15734
bl=0 pos[1]=[5.9 11.0 7.3] dr=1.23 t=0.0ps kin=1.48 pot=21.56 Rg=11.698 SPS=15824
bl=0 pos[1]=[6.4 11.7 6.7] dr=1.27 t=0.0ps kin=1.49 pot=21.57 Rg=11.685 SPS=14496
bl=0 pos[1]=[5.1 10.6 7.1] dr=1.28 t=0.0ps kin=1.45 pot=21.58 Rg=11.718 SPS=15761
bl=0 pos[1]=[5.9 11.5 7.7] dr=1.28 t=0.0ps kin=1.50 pot=21.65 Rg=11.758 SPS=14670
bl=0 pos[1]=[5.2 12.4 9.5] dr=1.28 t=0.0ps kin=1.58 pot=21.62 Rg=11.786 SPS=15269
bl=0 pos[1]=[4.3 15.6 9.0] dr=1.28 t=0.0ps kin=1.58 pot=21.62 Rg=11.825 SPS=14865
bl=0 pos[1]=[3.1 16.4 8.7] dr=1.27 t=0.0ps kin=1.50 pot=21.63 Rg=11.858 SPS=14738
bl=0 pos[1]=[-0.2 16.2 7.4] dr=1.25 t=0.0ps kin=1.51 pot=21.57 Rg=11.884 SPS=15555
bl=0 pos[1]=[-1.8 15.8 7.6] dr=1.25 t=0.0ps kin=1.53 pot=21.59 Rg=11.942 SPS=15647
bl=0 pos[1]=[-2.7 16.2 7.5] dr=1.25 t=0.0ps kin=1.49 pot=21.61 Rg=12.043 SPS=14928
bl=0 pos[1]=[-1.4 15.3 5.2] dr=1.26 t=0.0ps kin=1.50 pot=21.60 Rg=12.078 SPS=15003
bl=0 pos[1]=[-0.2 16.6 3.1] dr=1.26 t=0.0ps kin=1.54 pot=21.57 Rg=12.096 SPS=14607
bl=0 pos[1]=[-1.1 18.1 1.9] dr=1.26 t=0.0ps kin=1.51 pot=21.62 Rg=12.087 SPS=16222
bl=0 pos[1]=[-1.2 18.5 3.5] dr=1.28 t=0.0ps kin=1.55 pot=21.56 Rg=12.080 SPS=15733
bl=0 pos[1]=[-1.1 19.5 4.9] dr=1.24 t=0.0ps kin=1.49 pot=21.60 Rg=12.058 SPS=14348
bl=0 pos[1]=[0.4 19.3 5.7] dr=1.23 t=0.0ps kin=1.47 pot=21.64 Rg=12.019 SPS=15288
bl=0 pos[1]=[1.3 19.4 5.9] dr=1.21 t=0.0ps kin=1.51 pot=21.58 Rg=12.012 SPS=15142
bl=0 pos[1]=[3.2 19.3 6.5] dr=1.24 t=0.0ps kin=1.47 pot=21.56 Rg=12.011 SPS=14530
bl=0 pos[1]=[4.3 16.5 6.2] dr=1.26 t=0.0ps kin=1.46 pot=21.57 Rg=12.000 SPS=15494
bl=0 pos[1]=[4.4 16.8 7.6] dr=1.27 t=0.0ps kin=1.48 pot=21.60 Rg=11.985 SPS=14728
bl=0 pos[1]=[6.4 17.2 7.9] dr=1.25 t=0.0ps kin=1.58 pot=21.59 Rg=11.979 SPS=14530
bl=0 pos[1]=[6.3 18.8 7.0] dr=1.24 t=0.0ps kin=1.47 pot=21.62 Rg=11.980 SPS=15446
bl=0 pos[1]=[6.3 18.3 6.1] dr=1.25 t=0.0ps kin=1.50 pot=21.60 Rg=11.987 SPS=15151
bl=0 pos[1]=[5.4 16.6 5.5] dr=1.26 t=0.0ps kin=1.49 pot=21.58 Rg=12.003 SPS=15394
bl=0 pos[1]=[5.7 16.3 6.6] dr=1.25 t=0.0ps kin=1.53 pot=21.55 Rg=12.009 SPS=15505
bl=0 pos[1]=[4.9 17.0 6.5] dr=1.24 t=0.0ps kin=1.52 pot=21.57 Rg=11.978 SPS=15430
bl=0 pos[1]=[6.1 17.0 8.5] dr=1.25 t=0.0ps kin=1.50 pot=21.58 Rg=11.906 SPS=15184
bl=0 pos[1]=[5.2 16.0 10.7] dr=1.25 t=0.0ps kin=1.51 pot=21.60 Rg=11.851 SPS=14792
bl=0 pos[1]=[5.8 14.2 10.1] dr=1.25 t=0.0ps kin=1.53 pot=21.64 Rg=11.769 SPS=16055
bl=0 pos[1]=[6.8 14.2 9.2] dr=1.27 t=0.0ps kin=1.52 pot=21.61 Rg=11.732 SPS=15562
bl=0 pos[1]=[8.1 13.6 7.7] dr=1.29 t=0.0ps kin=1.49 pot=21.62 Rg=11.713 SPS=15357
bl=0 pos[1]=[8.7 11.6 6.4] dr=1.28 t=0.0ps kin=1.51 pot=21.64 Rg=11.724 SPS=15244
bl=0 pos[1]=[8.0 10.5 7.2] dr=1.29 t=0.0ps kin=1.52 pot=21.59 Rg=11.728 SPS=15944
bl=0 pos[1]=[8.4 10.5 7.2] dr=1.26 t=0.0ps kin=1.45 pot=21.59 Rg=11.685 SPS=15161
bl=0 pos[1]=[8.5 11.2 6.0] dr=1.26 t=0.0ps kin=1.49 pot=21.58 Rg=11.627 SPS=15916
bl=0 pos[1]=[7.5 12.1 5.7] dr=1.24 t=0.0ps kin=1.49 pot=21.57 Rg=11.596 SPS=15792
bl=0 pos[1]=[9.2 13.8 4.4] dr=1.27 t=0.0ps kin=1.53 pot=21.56 Rg=11.644 SPS=15226
bl=0 pos[1]=[10.3 14.3 5.7] dr=1.30 t=0.0ps kin=1.51 pot=21.57 Rg=11.698 SPS=16051
bl=0 pos[1]=[12.0 14.2 5.9] dr=1.27 t=0.0ps kin=1.53 pot=21.55 Rg=11.748 SPS=15785
bl=0 pos[1]=[12.6 14.4 6.3] dr=1.29 t=0.0ps kin=1.54 pot=21.55 Rg=11.754 SPS=14919
bl=0 pos[1]=[12.7 13.9 5.8] dr=1.26 t=0.0ps kin=1.49 pot=21.59 Rg=11.698 SPS=14968
bl=0 pos[1]=[11.5 13.4 5.2] dr=1.23 t=0.0ps kin=1.45 pot=21.58 Rg=11.655 SPS=15456
bl=0 pos[1]=[10.8 14.4 4.5] dr=1.24 t=0.0ps kin=1.44 pot=21.56 Rg=11.647 SPS=15596
bl=0 pos[1]=[10.5 12.2 3.7] dr=1.25 t=0.0ps kin=1.50 pot=21.55 Rg=11.645 SPS=15459
bl=0 pos[1]=[9.1 11.7 1.9] dr=1.24 t=0.0ps kin=1.53 pot=21.53 Rg=11.595 SPS=15617
bl=0 pos[1]=[7.7 12.9 -0.3] dr=1.26 t=0.0ps kin=1.55 pot=21.60 Rg=11.539 SPS=14885
bl=0 pos[1]=[8.1 13.5 -0.5] dr=1.28 t=0.0ps kin=1.51 pot=21.65 Rg=11.524 SPS=15627
bl=0 pos[1]=[7.3 12.0 -1.2] dr=1.28 t=0.0ps kin=1.56 pot=21.58 Rg=11.514 SPS=13846
bl=0 pos[1]=[6.5 11.7 -1.3] dr=1.29 t=0.0ps kin=1.50 pot=21.57 Rg=11.507 SPS=15385
bl=0 pos[1]=[5.5 10.5 -0.8] dr=1.28 t=0.0ps kin=1.47 pot=21.59 Rg=11.480 SPS=15372
bl=0 pos[1]=[6.0 10.3 0.8] dr=1.27 t=0.0ps kin=1.53 pot=21.56 Rg=11.428 SPS=14014
bl=0 pos[1]=[6.8 11.4 2.5] dr=1.26 t=0.0ps kin=1.50 pot=21.60 Rg=11.403 SPS=14677
bl=0 pos[1]=[6.6 10.9 0.3] dr=1.27 t=0.0ps kin=1.52 pot=21.61 Rg=11.429 SPS=14508
bl=0 pos[1]=[7.6 12.9 -0.2] dr=1.25 t=0.0ps kin=1.48 pot=21.61 Rg=11.502 SPS=15046
bl=0 pos[1]=[9.2 13.5 -1.6] dr=1.27 t=0.0ps kin=1.51 pot=21.57 Rg=11.571 SPS=14682
bl=0 pos[1]=[8.7 13.0 -2.2] dr=1.29 t=0.0ps kin=1.51 pot=21.59 Rg=11.597 SPS=14980
bl=0 pos[1]=[8.8 12.2 -1.9] dr=1.27 t=0.0ps kin=1.46 pot=21.62 Rg=11.592 SPS=14717
bl=0 pos[1]=[9.1 13.0 -3.1] dr=1.26 t=0.0ps kin=1.51 pot=21.58 Rg=11.599 SPS=13925
bl=0 pos[1]=[10.4 13.0 -2.6] dr=1.27 t=0.0ps kin=1.53 pot=21.59 Rg=11.589 SPS=14274
bl=0 pos[1]=[10.0 13.1 -3.1] dr=1.28 t=0.0ps kin=1.52 pot=21.61 Rg=11.542 SPS=14183
bl=0 pos[1]=[9.0 12.2 -3.2] dr=1.26 t=0.0ps kin=1.51 pot=21.59 Rg=11.479 SPS=15679
bl=0 pos[1]=[7.2 12.6 -4.5] dr=1.26 t=0.0ps kin=1.51 pot=21.61 Rg=11.404 SPS=15416
bl=0 pos[1]=[6.2 14.6 -4.4] dr=1.28 t=0.0ps kin=1.52 pot=21.63 Rg=11.348 SPS=16134
bl=0 pos[1]=[6.0 14.5 -4.8] dr=1.27 t=0.0ps kin=1.52 pot=21.56 Rg=11.367 SPS=14773
bl=0 pos[1]=[7.0 13.6 -3.6] dr=1.27 t=0.0ps kin=1.50 pot=21.54 Rg=11.342 SPS=14545
bl=0 pos[1]=[7.7 13.2 -2.9] dr=1.27 t=0.0ps kin=1.51 pot=21.50 Rg=11.309 SPS=13732
bl=0 pos[1]=[8.2 13.7 -4.0] dr=1.26 t=0.0ps kin=1.45 pot=21.57 Rg=11.276 SPS=13899
bl=0 pos[1]=[8.9 14.5 -3.2] dr=1.24 t=0.0ps kin=1.50 pot=21.54 Rg=11.265 SPS=15524
bl=0 pos[1]=[7.9 13.2 -2.2] dr=1.24 t=0.0ps kin=1.47 pot=21.63 Rg=11.261 SPS=14919
bl=0 pos[1]=[9.7 12.0 -1.6] dr=1.26 t=0.0ps kin=1.49 pot=21.59 Rg=11.252 SPS=14731
bl=0 pos[1]=[11.3 12.0 -1.3] dr=1.26 t=0.0ps kin=1.50 pot=21.58 Rg=11.293 SPS=15438
bl=0 pos[1]=[11.5 10.7 -0.3] dr=1.25 t=0.0ps kin=1.51 pot=21.60 Rg=11.355 SPS=14726
bl=0 pos[1]=[12.6 11.8 -0.0] dr=1.27 t=0.0ps kin=1.51 pot=21.61 Rg=11.453 SPS=15509
bl=0 pos[1]=[12.1 13.2 -0.6] dr=1.26 t=0.0ps kin=1.48 pot=21.60 Rg=11.555 SPS=15026
bl=0 pos[1]=[11.3 14.9 0.7] dr=1.24 t=0.0ps kin=1.49 pot=21.56 Rg=11.594 SPS=14180
bl=0 pos[1]=[9.8 14.1 2.2] dr=1.25 t=0.0ps kin=1.50 pot=21.60 Rg=11.619 SPS=15133
bl=0 pos[1]=[10.3 13.5 2.7] dr=1.26 t=0.0ps kin=1.48 pot=21.61 Rg=11.687 SPS=14096
bl=0 pos[1]=[10.3 14.4 3.0] dr=1.23 t=0.0ps kin=1.49 pot=21.58 Rg=11.754 SPS=14581
bl=0 pos[1]=[10.0 14.1 3.9] dr=1.30 t=0.0ps kin=1.52 pot=21.56 Rg=11.824 SPS=15088
bl=0 pos[1]=[9.8 13.6 4.2] dr=1.29 t=0.0ps kin=1.49 pot=21.58 Rg=11.925 SPS=15267
bl=0 pos[1]=[9.2 12.9 4.7] dr=1.28 t=0.0ps kin=1.50 pot=21.58 Rg=12.006 SPS=15473
bl=0 pos[1]=[8.4 13.5 5.3] dr=1.23 t=0.0ps kin=1.47 pot=21.64 Rg=12.017 SPS=14429
bl=0 pos[1]=[7.9 13.6 5.1] dr=1.25 t=0.0ps kin=1.51 pot=21.61 Rg=12.007 SPS=15493
bl=0 pos[1]=[7.7 14.4 5.1] dr=1.22 t=0.0ps kin=1.50 pot=21.54 Rg=12.020 SPS=15450
bl=0 pos[1]=[8.6 13.9 5.7] dr=1.27 t=0.0ps kin=1.50 pot=21.55 Rg=12.052 SPS=14219
bl=0 pos[1]=[9.0 12.8 5.4] dr=1.24 t=0.0ps kin=1.43 pot=21.58 Rg=12.040 SPS=15895
bl=0 pos[1]=[7.7 11.5 5.5] dr=1.25 t=0.0ps kin=1.46 pot=21.56 Rg=12.008 SPS=15294
bl=0 pos[1]=[6.9 10.5 5.1] dr=1.27 t=0.0ps kin=1.48 pot=21.59 Rg=12.040 SPS=14041
bl=0 pos[1]=[5.2 11.6 4.3] dr=1.26 t=0.0ps kin=1.45 pot=21.62 Rg=12.074 SPS=15221
bl=0 pos[1]=[4.8 12.7 4.2] dr=1.23 t=0.0ps kin=1.48 pot=21.62 Rg=12.095 SPS=13564
bl=0 pos[1]=[6.1 14.3 4.9] dr=1.26 t=0.0ps kin=1.57 pot=21.57 Rg=12.159 SPS=15436
bl=0 pos[1]=[6.6 12.5 5.3] dr=1.25 t=0.0ps kin=1.48 pot=21.59 Rg=12.197 SPS=14542
bl=0 pos[1]=[7.1 13.1 4.7] dr=1.22 t=0.0ps kin=1.51 pot=21.59 Rg=12.208 SPS=14816
bl=0 pos[1]=[7.1 13.6 3.9] dr=1.29 t=0.0ps kin=1.52 pot=21.58 Rg=12.203 SPS=15181
bl=0 pos[1]=[6.3 14.7 4.1] dr=1.29 t=0.0ps kin=1.51 pot=21.58 Rg=12.224 SPS=14865
bl=0 pos[1]=[4.8 15.8 3.9] dr=1.27 t=0.0ps kin=1.53 pot=21.56 Rg=12.261 SPS=15235
bl=0 pos[1]=[4.9 16.1 1.9] dr=1.31 t=0.0ps kin=1.50 pot=21.59 Rg=12.343 SPS=14562
bl=0 pos[1]=[4.0 15.8 0.8] dr=1.32 t=0.0ps kin=1.52 pot=21.60 Rg=12.371 SPS=15078
bl=0 pos[1]=[4.6 14.2 -1.1] dr=1.28 t=0.0ps kin=1.57 pot=21.57 Rg=12.331 SPS=15550
bl=0 pos[1]=[6.8 14.2 -2.1] dr=1.28 t=0.0ps kin=1.50 pot=21.61 Rg=12.252 SPS=14613
bl=0 pos[1]=[6.7 15.1 -2.8] dr=1.26 t=0.0ps kin=1.48 pot=21.66 Rg=12.146 SPS=14280
bl=0 pos[1]=[6.8 15.4 -3.2] dr=1.27 t=0.0ps kin=1.49 pot=21.63 Rg=12.076 SPS=13838
bl=0 pos[1]=[7.0 15.6 -2.4] dr=1.24 t=0.0ps kin=1.51 pot=21.63 Rg=12.019 SPS=15198
bl=0 pos[1]=[8.0 15.0 -1.7] dr=1.24 t=0.0ps kin=1.47 pot=21.63 Rg=12.013 SPS=15690
bl=0 pos[1]=[8.9 16.4 -2.3] dr=1.24 t=0.0ps kin=1.54 pot=21.60 Rg=11.993 SPS=15286
bl=0 pos[1]=[9.5 14.9 -1.9] dr=1.28 t=0.0ps kin=1.53 pot=21.63 Rg=11.987 SPS=15493
bl=0 pos[1]=[10.3 13.4 -0.5] dr=1.31 t=0.0ps kin=1.56 pot=21.64 Rg=11.962 SPS=13843
bl=0 pos[1]=[9.8 13.0 -0.7] dr=1.28 t=0.0ps kin=1.53 pot=21.63 Rg=11.964 SPS=14480
bl=0 pos[1]=[9.6 12.9 -0.2] dr=1.26 t=0.0ps kin=1.50 pot=21.62 Rg=11.974 SPS=15336
bl=0 pos[1]=[9.3 13.5 -0.0] dr=1.24 t=0.0ps kin=1.51 pot=21.62 Rg=11.959 SPS=15035
bl=0 pos[1]=[7.9 12.3 -0.3] dr=1.25 t=0.0ps kin=1.52 pot=21.60 Rg=11.874 SPS=15215
bl=0 pos[1]=[9.3 10.4 0.2] dr=1.26 t=0.0ps kin=1.51 pot=21.58 Rg=11.855 SPS=14117
bl=0 pos[1]=[11.6 12.2 -1.0] dr=1.27 t=0.0ps kin=1.51 pot=21.59 Rg=11.853 SPS=15139
bl=0 pos[1]=[12.8 12.8 -3.5] dr=1.26 t=0.0ps kin=1.49 pot=21.59 Rg=11.843 SPS=15094
bl=0 pos[1]=[12.0 11.2 -5.3] dr=1.25 t=0.0ps kin=1.51 pot=21.54 Rg=11.866 SPS=14067
bl=0 pos[1]=[12.7 10.1 -5.1] dr=1.27 t=0.0ps kin=1.50 pot=21.58 Rg=11.896 SPS=15537
bl=0 pos[1]=[14.2 10.6 -3.1] dr=1.23 t=0.0ps kin=1.51 pot=21.58 Rg=11.919 SPS=15163
bl=0 pos[1]=[14.8 10.2 -2.2] dr=1.23 t=0.0ps kin=1.46 pot=21.60 Rg=11.928 SPS=15143
bl=0 pos[1]=[15.8 10.2 -4.5] dr=1.26 t=0.0ps kin=1.49 pot=21.61 Rg=11.903 SPS=15407
bl=0 pos[1]=[14.8 9.7 -2.9] dr=1.24 t=0.0ps kin=1.54 pot=21.60 Rg=11.861 SPS=14401
bl=0 pos[1]=[14.6 9.7 -2.6] dr=1.22 t=0.0ps kin=1.48 pot=21.63 Rg=11.843 SPS=15047
bl=0 pos[1]=[13.2 9.4 -3.2] dr=1.22 t=0.0ps kin=1.53 pot=21.63 Rg=11.825 SPS=13932
bl=0 pos[1]=[13.2 9.2 -3.4] dr=1.28 t=0.0ps kin=1.50 pot=21.63 Rg=11.796 SPS=14650
bl=0 pos[1]=[13.2 8.5 -3.3] dr=1.25 t=0.0ps kin=1.53 pot=21.57 Rg=11.761 SPS=15033
bl=0 pos[1]=[13.8 8.7 -3.5] dr=1.26 t=0.0ps kin=1.56 pot=21.64 Rg=11.664 SPS=14740
bl=0 pos[1]=[14.9 8.4 -3.5] dr=1.29 t=0.0ps kin=1.58 pot=21.61 Rg=11.591 SPS=15200
bl=0 pos[1]=[15.5 9.3 -4.0] dr=1.29 t=0.0ps kin=1.49 pot=21.59 Rg=11.538 SPS=14298
bl=0 pos[1]=[16.5 10.2 -3.8] dr=1.22 t=0.0ps kin=1.46 pot=21.58 Rg=11.537 SPS=15419
bl=0 pos[1]=[16.2 9.0 -3.8] dr=1.27 t=0.0ps kin=1.49 pot=21.54 Rg=11.516 SPS=15181
bl=0 pos[1]=[15.1 8.0 -5.6] dr=1.26 t=0.0ps kin=1.47 pot=21.59 Rg=11.468 SPS=14796
bl=0 pos[1]=[13.1 7.6 -6.1] dr=1.25 t=0.0ps kin=1.48 pot=21.61 Rg=11.479 SPS=14794
bl=0 pos[1]=[13.1 8.3 -5.3] dr=1.23 t=0.0ps kin=1.48 pot=21.59 Rg=11.508 SPS=14382
bl=0 pos[1]=[13.0 9.2 -5.6] dr=1.24 t=0.0ps kin=1.54 pot=21.56 Rg=11.515 SPS=15385
bl=0 pos[1]=[14.1 10.2 -6.4] dr=1.26 t=0.0ps kin=1.53 pot=21.53 Rg=11.566 SPS=15152
bl=0 pos[1]=[14.7 10.4 -7.1] dr=1.24 t=0.0ps kin=1.42 pot=21.62 Rg=11.614 SPS=14194
bl=0 pos[1]=[14.2 10.1 -6.1] dr=1.25 t=0.0ps kin=1.46 pot=21.62 Rg=11.613 SPS=15414
bl=0 pos[1]=[14.6 9.0 -3.5] dr=1.23 t=0.0ps kin=1.49 pot=21.58 Rg=11.600 SPS=14618
bl=0 pos[1]=[14.3 8.3 -2.3] dr=1.23 t=0.0ps kin=1.53 pot=21.59 Rg=11.579 SPS=15365
bl=0 pos[1]=[15.0 9.1 -1.5] dr=1.24 t=0.0ps kin=1.50 pot=21.59 Rg=11.530 SPS=15023
bl=0 pos[1]=[16.4 9.6 -1.4] dr=1.24 t=0.0ps kin=1.51 pot=21.63 Rg=11.517 SPS=14392
bl=0 pos[1]=[17.1 9.0 -1.8] dr=1.28 t=0.0ps kin=1.53 pot=21.63 Rg=11.553 SPS=15165
bl=0 pos[1]=[17.4 10.1 -2.8] dr=1.29 t=0.0ps kin=1.53 pot=21.67 Rg=11.569 SPS=14889
bl=0 pos[1]=[17.3 8.4 -3.9] dr=1.27 t=0.0ps kin=1.52 pot=21.60 Rg=11.567 SPS=15775
bl=0 pos[1]=[16.4 7.5 -3.3] dr=1.25 t=0.0ps kin=1.57 pot=21.60 Rg=11.534 SPS=16002
bl=0 pos[1]=[16.1 7.1 -2.8] dr=1.28 t=0.0ps kin=1.58 pot=21.61 Rg=11.528 SPS=14810
bl=0 pos[1]=[15.8 6.5 -2.3] dr=1.27 t=0.0ps kin=1.56 pot=21.58 Rg=11.512 SPS=16148
bl=0 pos[1]=[16.7 7.0 -1.7] dr=1.27 t=0.0ps kin=1.51 pot=21.62 Rg=11.469 SPS=14941
bl=0 pos[1]=[16.8 8.3 -2.3] dr=1.27 t=0.0ps kin=1.54 pot=21.58 Rg=11.495 SPS=14530
bl=0 pos[1]=[16.7 8.1 0.1] dr=1.26 t=0.0ps kin=1.51 pot=21.61 Rg=11.548 SPS=15569
bl=0 pos[1]=[15.3 8.6 2.5] dr=1.31 t=0.0ps kin=1.54 pot=21.62 Rg=11.606 SPS=15542
bl=0 pos[1]=[13.8 7.5 4.1] dr=1.28 t=0.0ps kin=1.50 pot=21.62 Rg=11.611 SPS=15047
bl=0 pos[1]=[12.5 7.4 4.8] dr=1.26 t=0.0ps kin=1.54 pot=21.63 Rg=11.615 SPS=15173
bl=0 pos[1]=[10.5 7.3 6.9] dr=1.26 t=0.0ps kin=1.49 pot=21.58 Rg=11.616 SPS=15766
bl=0 pos[1]=[9.5 6.5 6.8] dr=1.28 t=0.0ps kin=1.52 pot=21.57 Rg=11.648 SPS=14900
bl=0 pos[1]=[7.7 8.1 5.4] dr=1.27 t=0.0ps kin=1.52 pot=21.57 Rg=11.702 SPS=15392
bl=0 pos[1]=[7.4 10.8 6.7] dr=1.26 t=0.0ps kin=1.46 pot=21.62 Rg=11.756 SPS=15789
bl=0 pos[1]=[6.7 13.2 5.5] dr=1.25 t=0.0ps kin=1.52 pot=21.56 Rg=11.782 SPS=15868
bl=0 pos[1]=[7.1 14.8 4.4] dr=1.24 t=0.0ps kin=1.51 pot=21.59 Rg=11.711 SPS=15578
bl=0 pos[1]=[7.7 15.2 4.9] dr=1.27 t=0.0ps kin=1.45 pot=21.65 Rg=11.679 SPS=15016
bl=0 pos[1]=[6.9 15.1 5.7] dr=1.23 t=0.0ps kin=1.51 pot=21.63 Rg=11.663 SPS=15924
bl=0 pos[1]=[7.0 16.5 6.2] dr=1.27 t=0.0ps kin=1.50 pot=21.61 Rg=11.624 SPS=15485
bl=0 pos[1]=[7.8 15.9 5.2] dr=1.29 t=0.0ps kin=1.49 pot=21.59 Rg=11.553 SPS=15240
bl=0 pos[1]=[7.1 15.6 7.5] dr=1.24 t=0.0ps kin=1.53 pot=21.61 Rg=11.502 SPS=15795
bl=0 pos[1]=[6.8 14.5 7.4] dr=1.24 t=0.0ps kin=1.50 pot=21.67 Rg=11.530 SPS=15946
bl=0 pos[1]=[6.5 14.4 7.2] dr=1.26 t=0.0ps kin=1.55 pot=21.62 Rg=11.534 SPS=14095
bl=0 pos[1]=[6.7 14.8 6.0] dr=1.28 t=0.0ps kin=1.52 pot=21.62 Rg=11.528 SPS=15763
bl=0 pos[1]=[7.2 14.4 7.0] dr=1.24 t=0.0ps kin=1.52 pot=21.58 Rg=11.538 SPS=14915
bl=0 pos[1]=[6.5 12.6 7.2] dr=1.28 t=0.0ps kin=1.51 pot=21.63 Rg=11.556 SPS=14796
bl=0 pos[1]=[5.8 11.9 9.8] dr=1.26 t=0.0ps kin=1.50 pot=21.64 Rg=11.574 SPS=15816
bl=0 pos[1]=[4.6 11.5 9.3] dr=1.30 t=0.0ps kin=1.52 pot=21.60 Rg=11.539 SPS=15681
bl=0 pos[1]=[5.6 11.7 7.9] dr=1.28 t=0.0ps kin=1.48 pot=21.63 Rg=11.493 SPS=15412
bl=0 pos[1]=[6.3 10.5 7.6] dr=1.27 t=0.0ps kin=1.55 pot=21.58 Rg=11.470 SPS=15985
bl=0 pos[1]=[5.7 10.3 8.3] dr=1.27 t=0.0ps kin=1.55 pot=21.58 Rg=11.474 SPS=15454
bl=0 pos[1]=[5.9 10.6 7.7] dr=1.26 t=0.0ps kin=1.51 pot=21.59 Rg=11.447 SPS=15298
bl=0 pos[1]=[5.8 12.6 6.0] dr=1.26 t=0.0ps kin=1.49 pot=21.57 Rg=11.393 SPS=14544
bl=0 pos[1]=[4.2 12.8 7.0] dr=1.30 t=0.0ps kin=1.53 pot=21.54 Rg=11.342 SPS=14763
bl=0 pos[1]=[3.6 14.0 7.6] dr=1.28 t=0.0ps kin=1.56 pot=21.56 Rg=11.299 SPS=15452
bl=0 pos[1]=[4.1 15.4 6.7] dr=1.23 t=0.0ps kin=1.52 pot=21.51 Rg=11.274 SPS=15141
bl=0 pos[1]=[5.1 16.5 6.7] dr=1.24 t=0.0ps kin=1.52 pot=21.59 Rg=11.243 SPS=14290
bl=0 pos[1]=[5.0 17.1 6.6] dr=1.27 t=0.0ps kin=1.52 pot=21.53 Rg=11.273 SPS=14552
bl=0 pos[1]=[4.5 16.9 7.0] dr=1.23 t=0.0ps kin=1.41 pot=21.61 Rg=11.266 SPS=14291
bl=0 pos[1]=[3.7 15.2 7.0] dr=1.22 t=0.0ps kin=1.45 pot=21.60 Rg=11.295 SPS=14788
bl=0 pos[1]=[2.9 14.3 5.9] dr=1.24 t=0.0ps kin=1.53 pot=21.58 Rg=11.359 SPS=14975
bl=0 pos[1]=[2.4 14.7 5.7] dr=1.23 t=0.0ps kin=1.49 pot=21.61 Rg=11.396 SPS=15083
bl=0 pos[1]=[2.5 14.1 5.4] dr=1.22 t=0.0ps kin=1.51 pot=21.63 Rg=11.381 SPS=15412
bl=0 pos[1]=[3.0 14.1 5.5] dr=1.25 t=0.0ps kin=1.54 pot=21.58 Rg=11.362 SPS=14164
bl=0 pos[1]=[4.7 13.3 6.3] dr=1.24 t=0.0ps kin=1.48 pot=21.60 Rg=11.344 SPS=15530
bl=0 pos[1]=[7.0 15.1 5.9] dr=1.27 t=0.0ps kin=1.52 pot=21.56 Rg=11.321 SPS=15475
bl=0 pos[1]=[9.4 13.7 4.7] dr=1.29 t=0.0ps kin=1.50 pot=21.53 Rg=11.257 SPS=14924
bl=0 pos[1]=[10.0 14.0 4.4] dr=1.27 t=0.0ps kin=1.48 pot=21.57 Rg=11.220 SPS=15032
bl=0 pos[1]=[9.1 12.9 4.2] dr=1.28 t=0.0ps kin=1.54 pot=21.59 Rg=11.169 SPS=14250
bl=0 pos[1]=[8.0 12.7 1.9] dr=1.29 t=0.0ps kin=1.54 pot=21.58 Rg=11.142 SPS=15110
bl=0 pos[1]=[7.9 13.8 2.8] dr=1.29 t=0.0ps kin=1.50 pot=21.59 Rg=11.149 SPS=15232
bl=0 pos[1]=[9.0 13.2 4.5] dr=1.29 t=0.0ps kin=1.56 pot=21.56 Rg=11.201 SPS=14520
bl=0 pos[1]=[10.8 11.3 3.1] dr=1.27 t=0.0ps kin=1.52 pot=21.56 Rg=11.253 SPS=15080
bl=0 pos[1]=[10.8 10.2 3.0] dr=1.27 t=0.0ps kin=1.52 pot=21.58 Rg=11.294 SPS=14424
bl=0 pos[1]=[10.9 8.9 4.1] dr=1.31 t=0.0ps kin=1.48 pot=21.59 Rg=11.300 SPS=14721
bl=0 pos[1]=[10.4 8.3 5.9] dr=1.22 t=0.0ps kin=1.50 pot=21.57 Rg=11.270 SPS=13729
bl=0 pos[1]=[9.4 9.1 7.3] dr=1.22 t=0.0ps kin=1.52 pot=21.55 Rg=11.234 SPS=14649
bl=0 pos[1]=[10.1 10.0 7.6] dr=1.29 t=0.0ps kin=1.52 pot=21.56 Rg=11.169 SPS=15054
bl=0 pos[1]=[10.6 10.5 6.3] dr=1.29 t=0.0ps kin=1.54 pot=21.55 Rg=11.110 SPS=13729
bl=0 pos[1]=[9.6 10.5 6.4] dr=1.27 t=0.0ps kin=1.49 pot=21.56 Rg=11.100 SPS=14563
bl=0 pos[1]=[9.1 9.7 5.5] dr=1.23 t=0.0ps kin=1.50 pot=21.54 Rg=11.140 SPS=13892
bl=0 pos[1]=[8.5 8.6 4.4] dr=1.24 t=0.0ps kin=1.47 pot=21.53 Rg=11.189 SPS=15372
bl=0 pos[1]=[8.4 7.3 2.7] dr=1.26 t=0.0ps kin=1.47 pot=21.59 Rg=11.219 SPS=14535
bl=0 pos[1]=[6.9 6.1 1.7] dr=1.24 t=0.0ps kin=1.44 pot=21.56 Rg=11.298 SPS=15139
bl=0 pos[1]=[6.2 6.2 2.3] dr=1.20 t=0.0ps kin=1.43 pot=21.58 Rg=11.353 SPS=15267
bl=0 pos[1]=[6.9 5.7 1.9] dr=1.21 t=0.0ps kin=1.53 pot=21.56 Rg=11.377 SPS=14895
bl=0 pos[1]=[7.2 5.4 1.5] dr=1.28 t=0.0ps kin=1.50 pot=21.55 Rg=11.399 SPS=15297
bl=0 pos[1]=[8.7 5.7 1.4] dr=1.22 t=0.0ps kin=1.44 pot=21.58 Rg=11.387 SPS=15510
bl=0 pos[1]=[9.3 8.0 2.9] dr=1.21 t=0.0ps kin=1.52 pot=21.56 Rg=11.358 SPS=15264
bl=0 pos[1]=[8.0 7.3 3.0] dr=1.25 t=0.0ps kin=1.48 pot=21.55 Rg=11.328 SPS=15202
bl=0 pos[1]=[7.3 7.8 4.1] dr=1.25 t=0.0ps kin=1.45 pot=21.60 Rg=11.318 SPS=14731
bl=0 pos[1]=[5.9 7.8 2.6] dr=1.25 t=0.0ps kin=1.51 pot=21.61 Rg=11.329 SPS=15523
bl=0 pos[1]=[4.7 8.4 4.0] dr=1.29 t=0.0ps kin=1.55 pot=21.57 Rg=11.348 SPS=14674
bl=0 pos[1]=[3.7 10.7 5.5] dr=1.27 t=0.0ps kin=1.47 pot=21.56 Rg=11.320 SPS=14519
bl=0 pos[1]=[5.0 11.6 4.6] dr=1.27 t=0.0ps kin=1.49 pot=21.57 Rg=11.254 SPS=14949
bl=0 pos[1]=[4.0 12.1 5.5] dr=1.28 t=0.0ps kin=1.45 pot=21.59 Rg=11.223 SPS=14205
bl=0 pos[1]=[2.5 12.2 5.4] dr=1.29 t=0.0ps kin=1.55 pot=21.55 Rg=11.221 SPS=15084
bl=0 pos[1]=[2.5 11.6 4.8] dr=1.26 t=0.0ps kin=1.48 pot=21.57 Rg=11.180 SPS=15772
bl=0 pos[1]=[1.6 12.1 4.1] dr=1.28 t=0.0ps kin=1.56 pot=21.57 Rg=11.139 SPS=14671
bl=0 pos[1]=[2.5 11.6 4.0] dr=1.28 t=0.0ps kin=1.52 pot=21.61 Rg=11.111 SPS=14798
bl=0 pos[1]=[4.7 13.0 3.0] dr=1.24 t=0.0ps kin=1.51 pot=21.58 Rg=11.084 SPS=14057
bl=0 pos[1]=[6.0 15.3 4.5] dr=1.25 t=0.0ps kin=1.53 pot=21.57 Rg=11.027 SPS=14868
bl=0 pos[1]=[6.6 14.7 7.8] dr=1.24 t=0.0ps kin=1.47 pot=21.62 Rg=10.956 SPS=14636
bl=0 pos[1]=[6.0 14.0 8.8] dr=1.25 t=0.0ps kin=1.51 pot=21.58 Rg=10.903 SPS=15513
bl=0 pos[1]=[5.0 13.6 8.2] dr=1.26 t=0.0ps kin=1.50 pot=21.56 Rg=10.861 SPS=14851
bl=0 pos[1]=[3.5 13.0 8.2] dr=1.24 t=0.0ps kin=1.53 pot=21.59 Rg=10.886 SPS=14377
bl=0 pos[1]=[3.6 11.5 7.5] dr=1.26 t=0.0ps kin=1.52 pot=21.57 Rg=10.917 SPS=12716
bl=0 pos[1]=[3.0 10.0 5.8] dr=1.27 t=0.0ps kin=1.51 pot=21.58 Rg=10.914 SPS=12583
bl=0 pos[1]=[2.1 11.0 5.4] dr=1.27 t=0.0ps kin=1.54 pot=21.54 Rg=10.925 SPS=14218
bl=0 pos[1]=[2.5 10.0 3.3] dr=1.25 t=0.0ps kin=1.51 pot=21.56 Rg=10.973 SPS=14327
bl=0 pos[1]=[2.8 10.2 3.8] dr=1.23 t=0.0ps kin=1.57 pot=21.59 Rg=11.021 SPS=14723
bl=0 pos[1]=[3.2 10.0 3.0] dr=1.27 t=0.0ps kin=1.53 pot=21.63 Rg=11.071 SPS=14693
bl=0 pos[1]=[3.5 10.6 4.3] dr=1.23 t=0.0ps kin=1.57 pot=21.62 Rg=11.169 SPS=14080
bl=0 pos[1]=[3.9 12.5 4.1] dr=1.25 t=0.0ps kin=1.49 pot=21.61 Rg=11.255 SPS=13813
bl=0 pos[1]=[4.7 12.8 6.0] dr=1.25 t=0.0ps kin=1.51 pot=21.63 Rg=11.350 SPS=15194
bl=0 pos[1]=[5.1 14.3 5.9] dr=1.27 t=0.0ps kin=1.52 pot=21.60 Rg=11.433 SPS=15519
bl=0 pos[1]=[6.3 14.8 5.4] dr=1.25 t=0.0ps kin=1.51 pot=21.67 Rg=11.472 SPS=14805
bl=0 pos[1]=[5.9 12.3 3.6] dr=1.26 t=0.0ps kin=1.55 pot=21.58 Rg=11.491 SPS=14567
bl=0 pos[1]=[5.1 12.0 4.3] dr=1.24 t=0.0ps kin=1.47 pot=21.60 Rg=11.509 SPS=14593
bl=0 pos[1]=[7.4 11.7 4.0] dr=1.22 t=0.0ps kin=1.55 pot=21.57 Rg=11.526 SPS=15199
bl=0 pos[1]=[8.9 11.7 3.6] dr=1.24 t=0.0ps kin=1.56 pot=21.56 Rg=11.519 SPS=15570
bl=0 pos[1]=[8.8 11.7 3.4] dr=1.26 t=0.0ps kin=1.53 pot=21.61 Rg=11.507 SPS=14521
bl=0 pos[1]=[8.4 11.6 3.3] dr=1.24 t=0.0ps kin=1.54 pot=21.58 Rg=11.446 SPS=14533
bl=0 pos[1]=[9.8 13.0 3.5] dr=1.21 t=0.0ps kin=1.46 pot=21.56 Rg=11.423 SPS=14008
bl=0 pos[1]=[9.7 13.2 3.7] dr=1.22 t=0.0ps kin=1.48 pot=21.58 Rg=11.414 SPS=14704
bl=0 pos[1]=[9.8 13.5 2.5] dr=1.25 t=0.0ps kin=1.53 pot=21.57 Rg=11.467 SPS=13919
bl=0 pos[1]=[10.3 12.3 3.4] dr=1.28 t=0.0ps kin=1.50 pot=21.58 Rg=11.477 SPS=15753
bl=0 pos[1]=[11.9 13.2 3.0] dr=1.30 t=0.0ps kin=1.47 pot=21.63 Rg=11.499 SPS=15615
bl=0 pos[1]=[11.2 12.8 2.4] dr=1.26 t=0.0ps kin=1.51 pot=21.60 Rg=11.559 SPS=15018
bl=0 pos[1]=[9.0 13.3 2.3] dr=1.26 t=0.0ps kin=1.56 pot=21.57 Rg=11.618 SPS=15626
bl=0 pos[1]=[11.6 12.7 1.7] dr=1.28 t=0.0ps kin=1.54 pot=21.63 Rg=11.612 SPS=15235
bl=0 pos[1]=[10.2 11.8 3.4] dr=1.27 t=0.0ps kin=1.53 pot=21.64 Rg=11.646 SPS=15630
bl=0 pos[1]=[10.0 13.1 4.4] dr=1.24 t=0.0ps kin=1.52 pot=21.67 Rg=11.681 SPS=15562
bl=0 pos[1]=[10.0 14.4 5.3] dr=1.27 t=0.0ps kin=1.52 pot=21.62 Rg=11.703 SPS=15510
bl=0 pos[1]=[10.4 16.3 4.6] dr=1.26 t=0.0ps kin=1.55 pot=21.61 Rg=11.707 SPS=15191
bl=0 pos[1]=[9.4 16.7 2.9] dr=1.30 t=0.0ps kin=1.53 pot=21.56 Rg=11.679 SPS=15225
bl=0 pos[1]=[8.7 17.3 3.1] dr=1.28 t=0.0ps kin=1.49 pot=21.57 Rg=11.697 SPS=14355
bl=0 pos[1]=[9.0 16.5 4.7] dr=1.29 t=0.0ps kin=1.46 pot=21.55 Rg=11.700 SPS=15479
bl=0 pos[1]=[9.3 17.0 5.1] dr=1.27 t=0.0ps kin=1.42 pot=21.52 Rg=11.673 SPS=15146
bl=0 pos[1]=[8.2 18.2 4.6] dr=1.27 t=0.0ps kin=1.46 pot=21.53 Rg=11.650 SPS=14949
bl=0 pos[1]=[7.4 16.1 3.2] dr=1.23 t=0.0ps kin=1.44 pot=21.56 Rg=11.630 SPS=15015
bl=0 pos[1]=[9.5 15.9 4.1] dr=1.22 t=0.0ps kin=1.43 pot=21.57 Rg=11.649 SPS=14963
bl=0 pos[1]=[8.4 16.8 6.1] dr=1.22 t=0.0ps kin=1.47 pot=21.58 Rg=11.638 SPS=15988
bl=0 pos[1]=[8.8 15.5 5.6] dr=1.23 t=0.0ps kin=1.48 pot=21.61 Rg=11.576 SPS=14920
bl=0 pos[1]=[9.1 14.9 6.1] dr=1.29 t=0.0ps kin=1.45 pot=21.58 Rg=11.494 SPS=14556
bl=0 pos[1]=[6.9 14.8 5.3] dr=1.27 t=0.0ps kin=1.48 pot=21.55 Rg=11.450 SPS=15414
bl=0 pos[1]=[6.0 14.1 4.1] dr=1.26 t=0.0ps kin=1.50 pot=21.60 Rg=11.388 SPS=14576
bl=0 pos[1]=[7.5 15.8 3.5] dr=1.25 t=0.0ps kin=1.46 pot=21.60 Rg=11.359 SPS=14722
bl=0 pos[1]=[7.6 16.9 5.0] dr=1.22 t=0.0ps kin=1.45 pot=21.61 Rg=11.368 SPS=14909
bl=0 pos[1]=[8.2 17.5 4.4] dr=1.25 t=0.0ps kin=1.50 pot=21.57 Rg=11.354 SPS=14674
bl=0 pos[1]=[7.2 18.8 4.2] dr=1.28 t=0.0ps kin=1.50 pot=21.55 Rg=11.375 SPS=15420
bl=0 pos[1]=[7.5 19.5 2.7] dr=1.24 t=0.0ps kin=1.50 pot=21.58 Rg=11.387 SPS=15534
bl=0 pos[1]=[6.9 19.5 1.4] dr=1.25 t=0.0ps kin=1.47 pot=21.55 Rg=11.379 SPS=14480
bl=0 pos[1]=[5.2 19.6 1.1] dr=1.27 t=0.0ps kin=1.52 pot=21.57 Rg=11.374 SPS=15908
bl=0 pos[1]=[4.1 19.3 1.4] dr=1.29 t=0.0ps kin=1.53 pot=21.59 Rg=11.356 SPS=14851
bl=0 pos[1]=[3.2 19.6 0.6] dr=1.26 t=0.0ps kin=1.50 pot=21.61 Rg=11.374 SPS=15336
bl=0 pos[1]=[3.4 18.3 0.3] dr=1.30 t=0.0ps kin=1.52 pot=21.58 Rg=11.390 SPS=14647
bl=0 pos[1]=[2.7 17.1 0.6] dr=1.30 t=0.0ps kin=1.54 pot=21.63 Rg=11.451 SPS=14328
bl=0 pos[1]=[2.8 16.7 -0.6] dr=1.26 t=0.0ps kin=1.49 pot=21.63 Rg=11.474 SPS=14921
bl=0 pos[1]=[1.6 17.1 -0.7] dr=1.23 t=0.0ps kin=1.52 pot=21.63 Rg=11.430 SPS=14421
bl=0 pos[1]=[3.9 16.0 -1.8] dr=1.28 t=0.0ps kin=1.51 pot=21.61 Rg=11.425 SPS=15171
bl=0 pos[1]=[2.9 17.4 -1.0] dr=1.25 t=0.0ps kin=1.50 pot=21.62 Rg=11.427 SPS=15637
bl=0 pos[1]=[2.4 17.3 -1.5] dr=1.26 t=0.0ps kin=1.54 pot=21.57 Rg=11.437 SPS=14576
bl=0 pos[1]=[2.1 16.8 -0.6] dr=1.30 t=0.0ps kin=1.52 pot=21.60 Rg=11.467 SPS=15411
bl=0 pos[1]=[2.9 16.8 -1.2] dr=1.26 t=0.0ps kin=1.48 pot=21.60 Rg=11.527 SPS=14686
bl=0 pos[1]=[4.4 16.7 -0.8] dr=1.26 t=0.0ps kin=1.48 pot=21.60 Rg=11.520 SPS=15479
bl=0 pos[1]=[6.1 15.3 0.4] dr=1.26 t=0.0ps kin=1.46 pot=21.55 Rg=11.467 SPS=14861
bl=0 pos[1]=[5.7 15.8 -0.2] dr=1.24 t=0.0ps kin=1.56 pot=21.54 Rg=11.455 SPS=13809
bl=0 pos[1]=[6.5 15.2 -1.1] dr=1.27 t=0.0ps kin=1.49 pot=21.58 Rg=11.511 SPS=15363
bl=0 pos[1]=[6.9 14.0 -1.2] dr=1.22 t=0.0ps kin=1.49 pot=21.61 Rg=11.543 SPS=14826
bl=0 pos[1]=[8.5 11.7 -2.7] dr=1.25 t=0.0ps kin=1.53 pot=21.62 Rg=11.552 SPS=15360
bl=0 pos[1]=[8.3 10.4 -2.8] dr=1.27 t=0.0ps kin=1.54 pot=21.56 Rg=11.548 SPS=15171
bl=0 pos[1]=[7.9 10.4 -4.1] dr=1.26 t=0.0ps kin=1.52 pot=21.58 Rg=11.512 SPS=14282
bl=0 pos[1]=[9.7 10.1 -4.3] dr=1.26 t=0.0ps kin=1.55 pot=21.60 Rg=11.494 SPS=15022
bl=0 pos[1]=[9.0 11.2 -5.0] dr=1.25 t=0.0ps kin=1.50 pot=21.61 Rg=11.499 SPS=14180
bl=0 pos[1]=[9.4 12.3 -5.1] dr=1.24 t=0.0ps kin=1.49 pot=21.58 Rg=11.561 SPS=14573
bl=0 pos[1]=[9.0 14.0 -4.2] dr=1.24 t=0.0ps kin=1.51 pot=21.61 Rg=11.578 SPS=14898
bl=0 pos[1]=[10.6 13.2 -2.2] dr=1.27 t=0.0ps kin=1.52 pot=21.64 Rg=11.520 SPS=14854
bl=0 pos[1]=[7.4 12.1 -0.8] dr=1.23 t=0.0ps kin=1.47 pot=21.61 Rg=11.465 SPS=15288
bl=0 pos[1]=[5.5 10.9 -1.2] dr=1.27 t=0.0ps kin=1.54 pot=21.58 Rg=11.408 SPS=14207
bl=0 pos[1]=[6.3 12.0 -1.1] dr=1.26 t=0.0ps kin=1.48 pot=21.61 Rg=11.358 SPS=13869
bl=0 pos[1]=[5.6 11.3 -1.3] dr=1.27 t=0.0ps kin=1.51 pot=21.60 Rg=11.285 SPS=14367
bl=0 pos[1]=[4.9 10.9 -0.8] dr=1.25 t=0.0ps kin=1.46 pot=21.62 Rg=11.221 SPS=15048
bl=0 pos[1]=[5.2 12.7 0.2] dr=1.27 t=0.0ps kin=1.58 pot=21.57 Rg=11.192 SPS=15402
bl=0 pos[1]=[4.9 14.6 -0.0] dr=1.26 t=0.0ps kin=1.55 pot=21.58 Rg=11.166 SPS=14803
bl=0 pos[1]=[3.5 16.9 -0.7] dr=1.30 t=0.0ps kin=1.52 pot=21.61 Rg=11.177 SPS=15296
bl=0 pos[1]=[2.2 17.7 -1.0] dr=1.27 t=0.0ps kin=1.54 pot=21.62 Rg=11.206 SPS=14397
bl=0 pos[1]=[2.2 18.0 -0.8] dr=1.23 t=0.0ps kin=1.53 pot=21.57 Rg=11.229 SPS=14957
bl=0 pos[1]=[1.8 17.2 -1.6] dr=1.23 t=0.0ps kin=1.50 pot=21.60 Rg=11.248 SPS=15029
bl=0 pos[1]=[2.2 16.2 -1.8] dr=1.23 t=0.0ps kin=1.47 pot=21.61 Rg=11.246 SPS=14561
bl=0 pos[1]=[1.9 15.1 -3.0] dr=1.24 t=0.0ps kin=1.47 pot=21.58 Rg=11.278 SPS=14217
bl=0 pos[1]=[2.1 16.6 -3.9] dr=1.23 t=0.0ps kin=1.52 pot=21.57 Rg=11.280 SPS=14562
bl=0 pos[1]=[4.4 17.7 -4.0] dr=1.29 t=0.0ps kin=1.54 pot=21.57 Rg=11.298 SPS=15042
bl=0 pos[1]=[3.8 16.1 -2.8] dr=1.29 t=0.0ps kin=1.48 pot=21.59 Rg=11.335 SPS=15223
bl=0 pos[1]=[4.2 16.1 -3.2] dr=1.24 t=0.0ps kin=1.45 pot=21.54 Rg=11.359 SPS=14188
bl=0 pos[1]=[3.1 14.6 -2.2] dr=1.29 t=0.0ps kin=1.52 pot=21.53 Rg=11.327 SPS=14593
bl=0 pos[1]=[1.1 14.7 -2.6] dr=1.26 t=0.0ps kin=1.52 pot=21.58 Rg=11.278 SPS=13957
bl=0 pos[1]=[-0.1 15.0 -4.5] dr=1.28 t=0.0ps kin=1.57 pot=21.53 Rg=11.318 SPS=14739
bl=0 pos[1]=[1.2 13.7 -6.3] dr=1.29 t=0.0ps kin=1.52 pot=21.57 Rg=11.369 SPS=14175
bl=0 pos[1]=[4.1 14.1 -6.2] dr=1.26 t=0.0ps kin=1.51 pot=21.55 Rg=11.403 SPS=15117
bl=0 pos[1]=[5.4 12.8 -6.0] dr=1.24 t=0.0ps kin=1.45 pot=21.57 Rg=11.451 SPS=14449
bl=0 pos[1]=[6.5 12.4 -7.1] dr=1.26 t=0.0ps kin=1.50 pot=21.51 Rg=11.451 SPS=14637
bl=0 pos[1]=[5.2 14.5 -7.4] dr=1.30 t=0.0ps kin=1.48 pot=21.57 Rg=11.448 SPS=15339
bl=0 pos[1]=[4.1 14.8 -6.0] dr=1.29 t=0.0ps kin=1.50 pot=21.54 Rg=11.420 SPS=14371
bl=0 pos[1]=[4.9 14.2 -5.1] dr=1.28 t=0.0ps kin=1.52 pot=21.58 Rg=11.393 SPS=14675
bl=0 pos[1]=[4.7 13.0 -5.1] dr=1.27 t=0.0ps kin=1.48 pot=21.62 Rg=11.342 SPS=13986
bl=0 pos[1]=[3.8 13.2 -6.4] dr=1.27 t=0.0ps kin=1.49 pot=21.56 Rg=11.343 SPS=14722
bl=0 pos[1]=[4.1 14.0 -7.0] dr=1.26 t=0.0ps kin=1.52 pot=21.59 Rg=11.372 SPS=13439
bl=0 pos[1]=[3.7 12.9 -8.8] dr=1.24 t=0.0ps kin=1.46 pot=21.56 Rg=11.389 SPS=14546
bl=0 pos[1]=[4.1 12.0 -9.3] dr=1.23 t=0.0ps kin=1.48 pot=21.58 Rg=11.372 SPS=13843
bl=0 pos[1]=[4.6 12.2 -8.9] dr=1.25 t=0.0ps kin=1.46 pot=21.62 Rg=11.355 SPS=15084
bl=0 pos[1]=[4.0 13.2 -7.9] dr=1.23 t=0.0ps kin=1.53 pot=21.59 Rg=11.432 SPS=15213
bl=0 pos[1]=[2.7 11.9 -7.9] dr=1.26 t=0.0ps kin=1.54 pot=21.55 Rg=11.481 SPS=14578
bl=0 pos[1]=[0.7 9.9 -6.6] dr=1.25 t=0.0ps kin=1.45 pot=21.55 Rg=11.528 SPS=15128
bl=0 pos[1]=[0.4 11.7 -9.3] dr=1.22 t=0.0ps kin=1.43 pot=21.56 Rg=11.512 SPS=14552
bl=0 pos[1]=[1.0 12.7 -10.3] dr=1.22 t=0.0ps kin=1.45 pot=21.58 Rg=11.487 SPS=15391
bl=0 pos[1]=[3.0 13.3 -9.1] dr=1.23 t=0.0ps kin=1.47 pot=21.63 Rg=11.491 SPS=15838
bl=0 pos[1]=[4.4 14.5 -7.6] dr=1.25 t=0.0ps kin=1.51 pot=21.59 Rg=11.486 SPS=14833
bl=0 pos[1]=[3.8 13.0 -8.4] dr=1.27 t=0.0ps kin=1.50 pot=21.57 Rg=11.471 SPS=14912
bl=0 pos[1]=[4.1 11.2 -10.1] dr=1.26 t=0.0ps kin=1.48 pot=21.63 Rg=11.423 SPS=13591
bl=0 pos[1]=[4.0 10.7 -11.3] dr=1.24 t=0.0ps kin=1.47 pot=21.60 Rg=11.428 SPS=14599
bl=0 pos[1]=[2.7 10.2 -11.8] dr=1.25 t=0.0ps kin=1.49 pot=21.59 Rg=11.416 SPS=15032
bl=0 pos[1]=[0.7 9.5 -11.9] dr=1.27 t=0.0ps kin=1.54 pot=21.64 Rg=11.419 SPS=14903
bl=0 pos[1]=[0.7 9.0 -10.7] dr=1.22 t=0.0ps kin=1.48 pot=21.63 Rg=11.432 SPS=15718
bl=0 pos[1]=[2.2 8.4 -10.2] dr=1.21 t=0.0ps kin=1.53 pot=21.61 Rg=11.457 SPS=15631
bl=0 pos[1]=[1.8 8.6 -8.8] dr=1.22 t=0.0ps kin=1.49 pot=21.59 Rg=11.491 SPS=15106
bl=0 pos[1]=[1.1 10.6 -8.8] dr=1.28 t=0.0ps kin=1.51 pot=21.57 Rg=11.516 SPS=15416
bl=0 pos[1]=[2.0 9.9 -8.2] dr=1.29 t=0.0ps kin=1.51 pot=21.59 Rg=11.515 SPS=14927
bl=0 pos[1]=[2.5 10.1 -8.0] dr=1.29 t=0.0ps kin=1.52 pot=21.57 Rg=11.513 SPS=15332
bl=0 pos[1]=[4.4 10.6 -6.9] dr=1.27 t=0.0ps kin=1.51 pot=21.59 Rg=11.522 SPS=15445
bl=0 pos[1]=[6.9 12.0 -7.3] dr=1.27 t=0.0ps kin=1.48 pot=21.61 Rg=11.553 SPS=14794
bl=0 pos[1]=[6.1 12.5 -7.3] dr=1.21 t=0.0ps kin=1.48 pot=21.63 Rg=11.570 SPS=15371
bl=0 pos[1]=[5.9 11.9 -5.9] dr=1.24 t=0.0ps kin=1.48 pot=21.59 Rg=11.572 SPS=15672
bl=0 pos[1]=[5.7 14.1 -4.6] dr=1.22 t=0.0ps kin=1.51 pot=21.59 Rg=11.555 SPS=14351
bl=0 pos[1]=[5.0 13.0 -5.4] dr=1.28 t=0.0ps kin=1.53 pot=21.59 Rg=11.530 SPS=14981
bl=0 pos[1]=[5.6 10.7 -5.3] dr=1.29 t=0.0ps kin=1.54 pot=21.57 Rg=11.488 SPS=11741
bl=0 pos[1]=[6.0 9.3 -3.7] dr=1.26 t=0.0ps kin=1.51 pot=21.60 Rg=11.435 SPS=15038
bl=0 pos[1]=[5.1 10.3 -3.5] dr=1.27 t=0.0ps kin=1.60 pot=21.61 Rg=11.431 SPS=12183
bl=0 pos[1]=[4.1 11.3 -3.9] dr=1.28 t=0.0ps kin=1.54 pot=21.63 Rg=11.376 SPS=13150
bl=0 pos[1]=[3.6 11.9 -3.7] dr=1.29 t=0.0ps kin=1.56 pot=21.58 Rg=11.259 SPS=14974
bl=0 pos[1]=[3.8 13.1 -3.9] dr=1.28 t=0.0ps kin=1.53 pot=21.58 Rg=11.116 SPS=15216
bl=0 pos[1]=[4.2 13.0 -2.3] dr=1.29 t=0.0ps kin=1.51 pot=21.63 Rg=11.020 SPS=13723
bl=0 pos[1]=[4.6 14.1 -2.5] dr=1.26 t=0.0ps kin=1.56 pot=21.57 Rg=11.033 SPS=14242
bl=0 pos[1]=[4.0 12.8 -3.6] dr=1.26 t=0.0ps kin=1.51 pot=21.58 Rg=11.030 SPS=14244
bl=0 pos[1]=[4.5 10.3 -2.5] dr=1.27 t=0.0ps kin=1.48 pot=21.61 Rg=11.017 SPS=14515
bl=0 pos[1]=[4.4 10.5 -2.6] dr=1.26 t=0.0ps kin=1.50 pot=21.55 Rg=11.017 SPS=13881
bl=0 pos[1]=[4.5 10.4 -3.4] dr=1.25 t=0.0ps kin=1.49 pot=21.56 Rg=11.006 SPS=14760
bl=0 pos[1]=[5.1 11.8 -2.7] dr=1.25 t=0.0ps kin=1.49 pot=21.54 Rg=11.022 SPS=13977
bl=0 pos[1]=[4.8 9.6 -0.7] dr=1.26 t=0.0ps kin=1.50 pot=21.57 Rg=11.030 SPS=14837
bl=0 pos[1]=[5.5 8.9 -1.1] dr=1.26 t=0.0ps kin=1.47 pot=21.58 Rg=11.040 SPS=15141
bl=0 pos[1]=[6.1 9.1 -1.4] dr=1.24 t=0.0ps kin=1.52 pot=21.59 Rg=11.071 SPS=15265
bl=0 pos[1]=[6.1 10.7 -0.7] dr=1.22 t=0.0ps kin=1.49 pot=21.60 Rg=11.084 SPS=15123
bl=0 pos[1]=[6.4 10.7 -0.7] dr=1.22 t=0.0ps kin=1.44 pot=21.60 Rg=11.095 SPS=13893
bl=0 pos[1]=[6.6 10.5 -2.1] dr=1.25 t=0.0ps kin=1.45 pot=21.55 Rg=11.096 SPS=15398
bl=0 pos[1]=[6.3 10.2 -2.1] dr=1.22 t=0.0ps kin=1.48 pot=21.54 Rg=11.052 SPS=13781
bl=0 pos[1]=[6.7 10.0 -4.4] dr=1.23 t=0.0ps kin=1.47 pot=21.54 Rg=11.002 SPS=14235
bl=0 pos[1]=[7.8 12.6 -6.7] dr=1.27 t=0.0ps kin=1.44 pot=21.60 Rg=10.955 SPS=14924
bl=0 pos[1]=[5.8 13.1 -8.5] dr=1.29 t=0.0ps kin=1.52 pot=21.54 Rg=10.917 SPS=13963
bl=0 pos[1]=[5.4 13.7 -10.1] dr=1.29 t=0.0ps kin=1.52 pot=21.61 Rg=10.894 SPS=14865
bl=0 pos[1]=[6.0 14.3 -10.6] dr=1.29 t=0.0ps kin=1.50 pot=21.57 Rg=10.919 SPS=13846
bl=0 pos[1]=[5.9 14.7 -9.9] dr=1.25 t=0.0ps kin=1.48 pot=21.56 Rg=10.967 SPS=14738
bl=0 pos[1]=[6.2 13.8 -7.8] dr=1.24 t=0.0ps kin=1.41 pot=21.61 Rg=11.023 SPS=13638
bl=0 pos[1]=[7.0 12.0 -6.9] dr=1.26 t=0.0ps kin=1.44 pot=21.56 Rg=11.072 SPS=14812
bl=0 pos[1]=[7.1 12.3 -5.7] dr=1.25 t=0.0ps kin=1.43 pot=21.60 Rg=11.114 SPS=14703
bl=0 pos[1]=[7.0 11.3 -5.6] dr=1.27 t=0.0ps kin=1.53 pot=21.59 Rg=11.139 SPS=14840
bl=0 pos[1]=[7.4 9.8 -4.6] dr=1.26 t=0.0ps kin=1.47 pot=21.58 Rg=11.152 SPS=15275
bl=0 pos[1]=[7.5 9.9 -4.6] dr=1.25 t=0.0ps kin=1.52 pot=21.56 Rg=11.200 SPS=13919
bl=0 pos[1]=[7.9 10.2 -4.9] dr=1.24 t=0.0ps kin=1.44 pot=21.60 Rg=11.313 SPS=15044
bl=0 pos[1]=[8.0 10.5 -4.8] dr=1.26 t=0.0ps kin=1.51 pot=21.58 Rg=11.413 SPS=13203
bl=0 pos[1]=[8.0 12.0 -5.3] dr=1.29 t=0.0ps kin=1.51 pot=21.61 Rg=11.497 SPS=15229
bl=0 pos[1]=[7.2 14.6 -5.3] dr=1.26 t=0.0ps kin=1.54 pot=21.56 Rg=11.532 SPS=14474
bl=0 pos[1]=[6.0 17.7 -6.2] dr=1.30 t=0.0ps kin=1.48 pot=21.60 Rg=11.544 SPS=15012
bl=0 pos[1]=[4.0 18.3 -5.8] dr=1.28 t=0.0ps kin=1.48 pot=21.57 Rg=11.538 SPS=15584
bl=0 pos[1]=[2.3 18.0 -5.8] dr=1.27 t=0.0ps kin=1.49 pot=21.54 Rg=11.532 SPS=14679
bl=0 pos[1]=[2.6 16.8 -4.1] dr=1.25 t=0.0ps kin=1.47 pot=21.57 Rg=11.500 SPS=15210
bl=0 pos[1]=[4.0 18.6 -3.1] dr=1.24 t=0.0ps kin=1.51 pot=21.52 Rg=11.475 SPS=14006
bl=0 pos[1]=[4.8 20.2 -0.8] dr=1.28 t=0.0ps kin=1.52 pot=21.55 Rg=11.460 SPS=15301
bl=0 pos[1]=[4.2 19.7 -0.2] dr=1.25 t=0.0ps kin=1.49 pot=21.55 Rg=11.446 SPS=15187
bl=0 pos[1]=[3.9 19.0 0.1] dr=1.24 t=0.0ps kin=1.49 pot=21.60 Rg=11.412 SPS=14216
bl=0 pos[1]=[3.9 16.6 1.6] dr=1.25 t=0.0ps kin=1.52 pot=21.59 Rg=11.403 SPS=14868
bl=0 pos[1]=[2.4 16.5 1.7] dr=1.24 t=0.0ps kin=1.54 pot=21.58 Rg=11.404 SPS=14819
bl=0 pos[1]=[1.1 16.2 1.8] dr=1.26 t=0.0ps kin=1.49 pot=21.59 Rg=11.413 SPS=15821
bl=0 pos[1]=[0.8 16.2 0.3] dr=1.20 t=0.0ps kin=1.46 pot=21.63 Rg=11.405 SPS=15954
bl=0 pos[1]=[1.8 15.3 -0.2] dr=1.24 t=0.0ps kin=1.48 pot=21.61 Rg=11.387 SPS=14595
bl=0 pos[1]=[2.8 15.4 1.1] dr=1.22 t=0.0ps kin=1.42 pot=21.59 Rg=11.380 SPS=14904
bl=0 pos[1]=[5.2 14.5 0.8] dr=1.22 t=0.0ps kin=1.43 pot=21.56 Rg=11.348 SPS=14873
bl=0 pos[1]=[6.4 15.2 0.6] dr=1.25 t=0.0ps kin=1.46 pot=21.56 Rg=11.305 SPS=15001
bl=0 pos[1]=[5.9 15.4 1.5] dr=1.27 t=0.0ps kin=1.45 pot=21.61 Rg=11.284 SPS=14970
bl=0 pos[1]=[5.4 15.1 1.9] dr=1.25 t=0.0ps kin=1.47 pot=21.60 Rg=11.226 SPS=14288
bl=0 pos[1]=[4.1 14.8 3.0] dr=1.23 t=0.0ps kin=1.47 pot=21.60 Rg=11.239 SPS=14838
bl=0 pos[1]=[3.4 14.4 3.6] dr=1.22 t=0.0ps kin=1.52 pot=21.54 Rg=11.308 SPS=15626
bl=0 pos[1]=[2.8 14.5 2.6] dr=1.22 t=0.0ps kin=1.49 pot=21.56 Rg=11.373 SPS=13935
bl=0 pos[1]=[3.9 16.7 2.9] dr=1.23 t=0.0ps kin=1.51 pot=21.57 Rg=11.393 SPS=14753
bl=0 pos[1]=[4.5 17.1 1.1] dr=1.22 t=0.0ps kin=1.49 pot=21.61 Rg=11.420 SPS=13444
bl=0 pos[1]=[3.6 16.2 1.7] dr=1.27 t=0.0ps kin=1.49 pot=21.59 Rg=11.492 SPS=15691
bl=0 pos[1]=[3.9 16.6 1.6] dr=1.27 t=0.0ps kin=1.54 pot=21.55 Rg=11.569 SPS=14227
bl=0 pos[1]=[4.9 16.1 2.8] dr=1.24 t=0.0ps kin=1.45 pot=21.60 Rg=11.610 SPS=14708
bl=0 pos[1]=[4.2 17.8 2.2] dr=1.26 t=0.0ps kin=1.45 pot=21.63 Rg=11.638 SPS=14888
bl=0 pos[1]=[4.2 18.6 2.4] dr=1.25 t=0.0ps kin=1.50 pot=21.56 Rg=11.664 SPS=13730
bl=0 pos[1]=[2.7 18.2 2.4] dr=1.23 t=0.0ps kin=1.44 pot=21.62 Rg=11.688 SPS=14603
bl=0 pos[1]=[0.8 18.2 2.3] dr=1.25 t=0.0ps kin=1.49 pot=21.57 Rg=11.689 SPS=14354
bl=0 pos[1]=[0.6 17.3 2.3] dr=1.27 t=0.0ps kin=1.49 pot=21.59 Rg=11.690 SPS=15746
bl=0 pos[1]=[-0.7 15.7 2.4] dr=1.27 t=0.0ps kin=1.49 pot=21.59 Rg=11.687 SPS=14272
bl=0 pos[1]=[-0.4 15.2 1.7] dr=1.25 t=0.0ps kin=1.55 pot=21.60 Rg=11.682 SPS=15147
bl=0 pos[1]=[0.8 15.5 1.2] dr=1.25 t=0.0ps kin=1.50 pot=21.63 Rg=11.688 SPS=15091
bl=0 pos[1]=[1.6 14.4 0.9] dr=1.23 t=0.0ps kin=1.45 pot=21.59 Rg=11.709 SPS=14505
bl=0 pos[1]=[0.4 16.6 0.2] dr=1.24 t=0.0ps kin=1.51 pot=21.57 Rg=11.750 SPS=14911
bl=0 pos[1]=[0.2 17.9 -0.7] dr=1.25 t=0.0ps kin=1.48 pot=21.58 Rg=11.736 SPS=13680
bl=0 pos[1]=[-1.1 17.6 -0.7] dr=1.25 t=0.0ps kin=1.51 pot=21.58 Rg=11.694 SPS=14813
bl=0 pos[1]=[1.7 16.4 0.2] dr=1.27 t=0.0ps kin=1.55 pot=21.57 Rg=11.654 SPS=14507
bl=0 pos[1]=[2.7 14.1 1.7] dr=1.30 t=0.0ps kin=1.48 pot=21.60 Rg=11.636 SPS=15112
bl=0 pos[1]=[3.2 13.7 1.6] dr=1.27 t=0.0ps kin=1.52 pot=21.60 Rg=11.583 SPS=15066
bl=0 pos[1]=[3.4 14.1 3.0] dr=1.29 t=0.0ps kin=1.52 pot=21.58 Rg=11.501 SPS=14712
bl=0 pos[1]=[2.8 14.1 2.5] dr=1.27 t=0.0ps kin=1.49 pot=21.60 Rg=11.393 SPS=14858
bl=0 pos[1]=[3.4 13.8 2.7] dr=1.28 t=0.0ps kin=1.46 pot=21.62 Rg=11.295 SPS=14184
bl=0 pos[1]=[4.5 13.4 1.8] dr=1.28 t=0.0ps kin=1.49 pot=21.61 Rg=11.212 SPS=15126
bl=0 pos[1]=[4.7 12.4 2.5] dr=1.29 t=0.0ps kin=1.54 pot=21.63 Rg=11.160 SPS=14904
bl=0 pos[1]=[5.3 12.4 4.7] dr=1.29 t=0.0ps kin=1.57 pot=21.61 Rg=11.158 SPS=14567
bl=0 pos[1]=[5.0 12.2 5.1] dr=1.27 t=0.0ps kin=1.51 pot=21.59 Rg=11.132 SPS=14287
bl=0 pos[1]=[3.5 12.4 4.0] dr=1.25 t=0.0ps kin=1.50 pot=21.61 Rg=11.094 SPS=14565
bl=0 pos[1]=[4.1 11.2 4.1] dr=1.27 t=0.0ps kin=1.52 pot=21.58 Rg=11.104 SPS=14916
bl=0 pos[1]=[5.1 11.0 4.2] dr=1.24 t=0.0ps kin=1.49 pot=21.60 Rg=11.150 SPS=13832
bl=0 pos[1]=[4.8 12.2 3.2] dr=1.25 t=0.0ps kin=1.56 pot=21.58 Rg=11.132 SPS=15024
bl=0 pos[1]=[4.7 13.0 3.5] dr=1.25 t=0.0ps kin=1.51 pot=21.63 Rg=11.079 SPS=14950
bl=0 pos[1]=[4.5 12.7 3.0] dr=1.25 t=0.0ps kin=1.52 pot=21.56 Rg=11.079 SPS=15128
bl=0 pos[1]=[4.3 12.4 1.8] dr=1.25 t=0.0ps kin=1.54 pot=21.58 Rg=11.084 SPS=14920
bl=0 pos[1]=[4.4 12.1 2.4] dr=1.25 t=0.0ps kin=1.51 pot=21.59 Rg=11.029 SPS=14472
bl=0 pos[1]=[5.1 11.9 3.3] dr=1.23 t=0.0ps kin=1.47 pot=21.61 Rg=11.022 SPS=15143
bl=0 pos[1]=[5.4 10.7 2.8] dr=1.23 t=0.0ps kin=1.51 pot=21.56 Rg=11.074 SPS=12343
bl=0 pos[1]=[5.8 11.1 1.6] dr=1.25 t=0.0ps kin=1.49 pot=21.56 Rg=11.119 SPS=14593
bl=0 pos[1]=[6.7 10.7 1.3] dr=1.26 t=0.0ps kin=1.54 pot=21.57 Rg=11.147 SPS=14153
bl=0 pos[1]=[5.8 11.1 1.3] dr=1.28 t=0.0ps kin=1.48 pot=21.58 Rg=11.147 SPS=14381
bl=0 pos[1]=[6.0 12.0 0.3] dr=1.27 t=0.0ps kin=1.50 pot=21.55 Rg=11.190 SPS=15121
bl=0 pos[1]=[6.8 13.3 -0.6] dr=1.25 t=0.0ps kin=1.50 pot=21.54 Rg=11.264 SPS=14265
bl=0 pos[1]=[6.6 15.5 -1.2] dr=1.30 t=0.0ps kin=1.48 pot=21.59 Rg=11.285 SPS=14862
bl=0 pos[1]=[5.2 16.3 -1.1] dr=1.28 t=0.0ps kin=1.46 pot=21.58 Rg=11.314 SPS=14936
bl=0 pos[1]=[4.4 16.6 -0.2] dr=1.27 t=0.0ps kin=1.48 pot=21.61 Rg=11.346 SPS=15462
bl=0 pos[1]=[3.6 15.7 -1.3] dr=1.25 t=0.0ps kin=1.55 pot=21.59 Rg=11.394 SPS=14998
bl=0 pos[1]=[4.4 14.2 -3.1] dr=1.24 t=0.0ps kin=1.53 pot=21.61 Rg=11.385 SPS=14996
bl=0 pos[1]=[4.1 14.5 -3.8] dr=1.25 t=0.0ps kin=1.52 pot=21.63 Rg=11.348 SPS=14541
bl=0 pos[1]=[5.3 15.4 -3.1] dr=1.24 t=0.0ps kin=1.49 pot=21.58 Rg=11.309 SPS=13760
bl=0 pos[1]=[3.4 16.3 -3.2] dr=1.23 t=0.0ps kin=1.46 pot=21.59 Rg=11.273 SPS=14258
bl=0 pos[1]=[4.1 15.2 -4.6] dr=1.23 t=0.0ps kin=1.49 pot=21.58 Rg=11.240 SPS=13516
bl=0 pos[1]=[3.7 12.6 -4.3] dr=1.21 t=0.0ps kin=1.46 pot=21.56 Rg=11.167 SPS=14826
bl=0 pos[1]=[1.8 12.1 -3.3] dr=1.24 t=0.0ps kin=1.49 pot=21.55 Rg=11.140 SPS=13501
bl=0 pos[1]=[1.6 13.5 -3.2] dr=1.24 t=0.0ps kin=1.50 pot=21.54 Rg=11.133 SPS=14859
bl=0 pos[1]=[2.0 14.5 -2.4] dr=1.24 t=0.0ps kin=1.50 pot=21.58 Rg=11.132 SPS=14580
bl=0 pos[1]=[2.4 15.4 -2.2] dr=1.24 t=0.0ps kin=1.56 pot=21.58 Rg=11.124 SPS=15186
bl=0 pos[1]=[1.0 14.7 -2.6] dr=1.26 t=0.0ps kin=1.49 pot=21.61 Rg=11.122 SPS=14980
bl=0 pos[1]=[1.4 14.8 -5.0] dr=1.28 t=0.0ps kin=1.45 pot=21.58 Rg=11.073 SPS=14327
bl=0 pos[1]=[1.2 16.0 -5.2] dr=1.23 t=0.0ps kin=1.45 pot=21.58 Rg=11.072 SPS=15121
bl=0 pos[1]=[2.4 16.5 -5.4] dr=1.27 t=0.0ps kin=1.51 pot=21.53 Rg=11.111 SPS=14535
bl=0 pos[1]=[4.1 15.5 -6.5] dr=1.25 t=0.0ps kin=1.47 pot=21.57 Rg=11.136 SPS=14729
bl=0 pos[1]=[5.2 13.6 -6.5] dr=1.24 t=0.0ps kin=1.48 pot=21.57 Rg=11.180 SPS=14774
bl=0 pos[1]=[5.1 13.5 -3.6] dr=1.25 t=0.0ps kin=1.49 pot=21.58 Rg=11.204 SPS=14328
bl=0 pos[1]=[4.6 13.9 -3.2] dr=1.26 t=0.0ps kin=1.53 pot=21.60 Rg=11.242 SPS=14581
bl=0 pos[1]=[4.6 13.9 -4.8] dr=1.28 t=0.0ps kin=1.57 pot=21.64 Rg=11.317 SPS=14719
bl=0 pos[1]=[5.7 13.9 -6.5] dr=1.25 t=0.0ps kin=1.50 pot=21.61 Rg=11.424 SPS=15486
bl=0 pos[1]=[4.0 13.2 -5.1] dr=1.26 t=0.0ps kin=1.49 pot=21.63 Rg=11.441 SPS=15547
bl=0 pos[1]=[3.5 14.1 -4.4] dr=1.27 t=0.0ps kin=1.51 pot=21.61 Rg=11.467 SPS=14295
bl=0 pos[1]=[2.5 15.2 -5.5] dr=1.24 t=0.0ps kin=1.49 pot=21.62 Rg=11.504 SPS=15586
bl=0 pos[1]=[3.5 15.9 -5.6] dr=1.25 t=0.0ps kin=1.54 pot=21.58 Rg=11.545 SPS=14878
bl=0 pos[1]=[4.0 16.3 -5.0] dr=1.29 t=0.0ps kin=1.53 pot=21.57 Rg=11.630 SPS=12155
bl=0 pos[1]=[4.3 14.7 -5.9] dr=1.31 t=0.0ps kin=1.53 pot=21.58 Rg=11.642 SPS=14441
bl=0 pos[1]=[4.5 13.4 -6.3] dr=1.26 t=0.0ps kin=1.53 pot=21.58 Rg=11.608 SPS=15325
bl=0 pos[1]=[5.0 12.1 -8.3] dr=1.31 t=0.0ps kin=1.56 pot=21.59 Rg=11.611 SPS=15282
bl=0 pos[1]=[3.4 13.4 -10.5] dr=1.27 t=0.0ps kin=1.54 pot=21.59 Rg=11.626 SPS=14968
bl=0 pos[1]=[4.7 14.7 -10.7] dr=1.24 t=0.0ps kin=1.42 pot=21.61 Rg=11.628 SPS=15159
bl=0 pos[1]=[6.6 14.4 -10.9] dr=1.19 t=0.0ps kin=1.46 pot=21.58 Rg=11.609 SPS=13759
bl=0 pos[1]=[6.9 12.5 -9.9] dr=1.23 t=0.0ps kin=1.47 pot=21.55 Rg=11.616 SPS=15099
bl=0 pos[1]=[7.6 11.8 -10.4] dr=1.27 t=0.0ps kin=1.46 pot=21.61 Rg=11.640 SPS=16083
bl=0 pos[1]=[8.6 9.9 -11.4] dr=1.28 t=0.0ps kin=1.54 pot=21.56 Rg=11.650 SPS=14773
bl=0 pos[1]=[9.5 8.6 -10.5] dr=1.27 t=0.0ps kin=1.52 pot=21.66 Rg=11.604 SPS=15338
bl=0 pos[1]=[8.7 8.5 -8.9] dr=1.27 t=0.0ps kin=1.55 pot=21.65 Rg=11.590 SPS=14406
bl=0 pos[1]=[11.4 8.8 -8.4] dr=1.27 t=0.0ps kin=1.53 pot=21.63 Rg=11.566 SPS=15060
bl=0 pos[1]=[12.2 9.2 -10.3] dr=1.27 t=0.0ps kin=1.53 pot=21.63 Rg=11.558 SPS=14607
bl=0 pos[1]=[12.1 9.5 -12.0] dr=1.27 t=0.0ps kin=1.52 pot=21.68 Rg=11.550 SPS=14066
bl=0 pos[1]=[12.5 8.4 -11.7] dr=1.27 t=0.0ps kin=1.50 pot=21.65 Rg=11.574 SPS=15142
bl=0 pos[1]=[11.5 8.5 -10.8] dr=1.22 t=0.0ps kin=1.50 pot=21.63 Rg=11.593 SPS=14627
bl=0 pos[1]=[10.6 8.9 -8.6] dr=1.26 t=0.0ps kin=1.53 pot=21.60 Rg=11.580 SPS=14580
bl=0 pos[1]=[11.0 8.8 -6.8] dr=1.26 t=0.0ps kin=1.50 pot=21.62 Rg=11.565 SPS=14483
bl=0 pos[1]=[11.0 7.9 -6.7] dr=1.27 t=0.0ps kin=1.48 pot=21.64 Rg=11.547 SPS=15337
bl=0 pos[1]=[10.8 7.2 -7.0] dr=1.24 t=0.0ps kin=1.50 pot=21.64 Rg=11.501 SPS=14749
bl=0 pos[1]=[10.7 6.6 -7.6] dr=1.28 t=0.0ps kin=1.50 pot=21.59 Rg=11.430 SPS=14549
bl=0 pos[1]=[8.7 5.6 -9.2] dr=1.27 t=0.0ps kin=1.50 pot=21.60 Rg=11.393 SPS=14934
bl=0 pos[1]=[8.5 5.7 -9.5] dr=1.29 t=0.0ps kin=1.49 pot=21.62 Rg=11.392 SPS=14225
bl=0 pos[1]=[9.7 5.2 -8.9] dr=1.26 t=0.0ps kin=1.46 pot=21.60 Rg=11.398 SPS=15503
bl=0 pos[1]=[7.9 4.5 -8.9] dr=1.27 t=0.0ps kin=1.49 pot=21.61 Rg=11.396 SPS=14900
bl=0 pos[1]=[6.0 5.5 -9.7] dr=1.25 t=0.0ps kin=1.53 pot=21.59 Rg=11.412 SPS=14630
bl=0 pos[1]=[4.8 5.6 -10.7] dr=1.26 t=0.0ps kin=1.50 pot=21.61 Rg=11.418 SPS=15524
bl=0 pos[1]=[5.4 5.0 -12.0] dr=1.27 t=0.0ps kin=1.51 pot=21.59 Rg=11.426 SPS=11907
bl=0 pos[1]=[6.8 4.3 -13.2] dr=1.28 t=0.0ps kin=1.51 pot=21.57 Rg=11.415 SPS=15206
bl=0 pos[1]=[7.1 3.8 -13.4] dr=1.29 t=0.0ps kin=1.50 pot=21.56 Rg=11.388 SPS=13823
bl=0 pos[1]=[7.5 4.5 -11.3] dr=1.25 t=0.0ps kin=1.50 pot=21.55 Rg=11.322 SPS=14900
bl=0 pos[1]=[6.7 5.7 -11.7] dr=1.26 t=0.0ps kin=1.50 pot=21.59 Rg=11.285 SPS=14703
bl=0 pos[1]=[7.2 4.2 -12.7] dr=1.27 t=0.0ps kin=1.54 pot=21.61 Rg=11.256 SPS=14417
bl=0 pos[1]=[7.6 6.0 -13.8] dr=1.27 t=0.0ps kin=1.56 pot=21.58 Rg=11.213 SPS=15203
bl=0 pos[1]=[7.0 6.6 -14.1] dr=1.26 t=0.0ps kin=1.52 pot=21.57 Rg=11.157 SPS=14053
bl=0 pos[1]=[8.3 7.4 -11.9] dr=1.29 t=0.0ps kin=1.54 pot=21.55 Rg=11.122 SPS=14993
bl=0 pos[1]=[8.5 7.1 -11.4] dr=1.28 t=0.0ps kin=1.48 pot=21.57 Rg=11.126 SPS=14203
bl=0 pos[1]=[7.8 7.1 -11.6] dr=1.28 t=0.0ps kin=1.56 pot=21.58 Rg=11.173 SPS=14913
bl=0 pos[1]=[9.5 7.1 -11.9] dr=1.29 t=0.0ps kin=1.52 pot=21.58 Rg=11.224 SPS=14177
bl=0 pos[1]=[9.6 8.1 -11.1] dr=1.23 t=0.0ps kin=1.45 pot=21.61 Rg=11.242 SPS=14714
bl=0 pos[1]=[10.5 8.2 -10.3] dr=1.25 t=0.0ps kin=1.47 pot=21.60 Rg=11.281 SPS=14948
bl=0 pos[1]=[10.0 8.4 -9.0] dr=1.24 t=0.0ps kin=1.51 pot=21.57 Rg=11.337 SPS=13961
bl=0 pos[1]=[8.4 9.3 -7.7] dr=1.22 t=0.0ps kin=1.49 pot=21.58 Rg=11.400 SPS=14744
bl=0 pos[1]=[7.8 10.0 -7.8] dr=1.22 t=0.0ps kin=1.52 pot=21.54 Rg=11.425 SPS=14162
bl=0 pos[1]=[7.3 10.8 -7.0] dr=1.24 t=0.0ps kin=1.48 pot=21.59 Rg=11.428 SPS=14925
bl=0 pos[1]=[8.4 9.1 -6.6] dr=1.26 t=0.0ps kin=1.54 pot=21.55 Rg=11.445 SPS=14424
bl=0 pos[1]=[6.5 7.4 -8.5] dr=1.27 t=0.0ps kin=1.50 pot=21.58 Rg=11.461 SPS=15230
bl=0 pos[1]=[4.3 8.0 -8.5] dr=1.25 t=0.0ps kin=1.46 pot=21.62 Rg=11.491 SPS=15375
bl=0 pos[1]=[5.0 6.4 -8.1] dr=1.23 t=0.0ps kin=1.51 pot=21.61 Rg=11.560 SPS=15371
bl=0 pos[1]=[4.4 7.2 -8.0] dr=1.24 t=0.0ps kin=1.51 pot=21.56 Rg=11.603 SPS=14887
bl=0 pos[1]=[4.3 8.2 -8.4] dr=1.26 t=0.0ps kin=1.54 pot=21.60 Rg=11.607 SPS=15105
bl=0 pos[1]=[5.7 8.7 -8.6] dr=1.24 t=0.0ps kin=1.48 pot=21.59 Rg=11.607 SPS=14981
bl=0 pos[1]=[4.2 7.0 -8.3] dr=1.24 t=0.0ps kin=1.56 pot=21.60 Rg=11.580 SPS=15146
bl=0 pos[1]=[4.4 5.3 -6.8] dr=1.30 t=0.0ps kin=1.58 pot=21.59 Rg=11.521 SPS=13328
bl=0 pos[1]=[4.6 5.1 -6.5] dr=1.28 t=0.0ps kin=1.49 pot=21.58 Rg=11.535 SPS=16195
bl=0 pos[1]=[3.6 4.8 -5.3] dr=1.25 t=0.0ps kin=1.48 pot=21.63 Rg=11.554 SPS=15638
bl=0 pos[1]=[3.4 5.6 -4.8] dr=1.24 t=0.0ps kin=1.50 pot=21.58 Rg=11.580 SPS=14816
bl=0 pos[1]=[4.4 6.1 -3.8] dr=1.24 t=0.0ps kin=1.46 pot=21.58 Rg=11.602 SPS=14917
bl=0 pos[1]=[4.4 7.4 -4.3] dr=1.26 t=0.0ps kin=1.48 pot=21.60 Rg=11.552 SPS=12700
bl=0 pos[1]=[4.2 7.8 -4.4] dr=1.21 t=0.0ps kin=1.52 pot=21.61 Rg=11.576 SPS=15436
bl=0 pos[1]=[6.2 9.8 -5.8] dr=1.25 t=0.0ps kin=1.56 pot=21.55 Rg=11.631 SPS=14398
bl=0 pos[1]=[7.2 9.8 -5.7] dr=1.29 t=0.0ps kin=1.50 pot=21.60 Rg=11.655 SPS=15138
bl=0 pos[1]=[7.8 10.7 -4.1] dr=1.28 t=0.0ps kin=1.52 pot=21.59 Rg=11.660 SPS=15052
bl=0 pos[1]=[8.6 9.2 -4.5] dr=1.29 t=0.0ps kin=1.50 pot=21.61 Rg=11.693 SPS=14961
bl=0 pos[1]=[7.2 10.2 -4.2] dr=1.28 t=0.0ps kin=1.55 pot=21.54 Rg=11.746 SPS=15520
bl=0 pos[1]=[8.0 12.6 -5.2] dr=1.29 t=0.0ps kin=1.50 pot=21.59 Rg=11.766 SPS=14401
bl=0 pos[1]=[8.2 12.9 -6.5] dr=1.24 t=0.0ps kin=1.48 pot=21.59 Rg=11.787 SPS=15644
bl=0 pos[1]=[9.7 12.9 -7.9] dr=1.20 t=0.0ps kin=1.50 pot=21.58 Rg=11.788 SPS=15606
bl=0 pos[1]=[12.1 12.0 -8.2] dr=1.23 t=0.0ps kin=1.48 pot=21.59 Rg=11.747 SPS=14639
bl=0 pos[1]=[13.4 12.0 -8.6] dr=1.20 t=0.0ps kin=1.50 pot=21.56 Rg=11.703 SPS=14765
bl=0 pos[1]=[13.7 12.2 -8.6] dr=1.22 t=0.0ps kin=1.45 pot=21.63 Rg=11.680 SPS=14261
bl=0 pos[1]=[14.6 12.7 -8.6] dr=1.22 t=0.0ps kin=1.50 pot=21.59 Rg=11.652 SPS=15184
bl=0 pos[1]=[13.0 13.0 -8.5] dr=1.24 t=0.0ps kin=1.56 pot=21.60 Rg=11.633 SPS=15003
bl=0 pos[1]=[12.3 12.3 -8.9] dr=1.25 t=0.0ps kin=1.52 pot=21.59 Rg=11.578 SPS=14386
bl=0 pos[1]=[11.0 12.9 -8.4] dr=1.26 t=0.0ps kin=1.49 pot=21.59 Rg=11.540 SPS=15033
bl=0 pos[1]=[10.1 12.9 -9.0] dr=1.23 t=0.0ps kin=1.51 pot=21.60 Rg=11.467 SPS=14169
bl=0 pos[1]=[10.8 12.3 -9.9] dr=1.29 t=0.0ps kin=1.52 pot=21.59 Rg=11.401 SPS=13994
bl=0 pos[1]=[10.8 12.1 -10.5] dr=1.29 t=0.0ps kin=1.51 pot=21.62 Rg=11.337 SPS=14157
bl=0 pos[1]=[10.3 13.5 -10.0] dr=1.26 t=0.0ps kin=1.50 pot=21.59 Rg=11.319 SPS=14781
bl=0 pos[1]=[11.9 13.7 -10.6] dr=1.24 t=0.0ps kin=1.55 pot=21.54 Rg=11.295 SPS=14839
bl=0 pos[1]=[11.4 12.6 -10.3] dr=1.26 t=0.0ps kin=1.49 pot=21.60 Rg=11.277 SPS=14477
bl=0 pos[1]=[10.3 12.8 -9.5] dr=1.26 t=0.0ps kin=1.50 pot=21.61 Rg=11.286 SPS=14677
bl=0 pos[1]=[11.8 13.1 -8.9] dr=1.28 t=0.0ps kin=1.57 pot=21.53 Rg=11.312 SPS=13968
bl=0 pos[1]=[11.2 13.8 -6.9] dr=1.27 t=0.0ps kin=1.50 pot=21.59 Rg=11.309 SPS=14679
bl=0 pos[1]=[10.9 14.4 -6.6] dr=1.26 t=0.0ps kin=1.53 pot=21.55 Rg=11.277 SPS=14271
bl=0 pos[1]=[9.6 13.4 -6.9] dr=1.26 t=0.0ps kin=1.47 pot=21.56 Rg=11.218 SPS=14968
bl=0 pos[1]=[9.0 12.7 -9.2] dr=1.27 t=0.0ps kin=1.50 pot=21.55 Rg=11.185 SPS=14328
bl=0 pos[1]=[7.5 13.7 -10.1] dr=1.26 t=0.0ps kin=1.52 pot=21.57 Rg=11.177 SPS=12810
bl=0 pos[1]=[4.7 14.9 -9.6] dr=1.25 t=0.0ps kin=1.47 pot=21.55 Rg=11.184 SPS=14137
bl=0 pos[1]=[4.2 15.6 -7.9] dr=1.22 t=0.0ps kin=1.46 pot=21.56 Rg=11.219 SPS=14580
bl=0 pos[1]=[5.6 16.8 -6.3] dr=1.21 t=0.0ps kin=1.52 pot=21.56 Rg=11.321 SPS=14409
bl=0 pos[1]=[6.3 17.9 -7.7] dr=1.26 t=0.0ps kin=1.49 pot=21.56 Rg=11.400 SPS=15007
bl=0 pos[1]=[4.6 16.6 -8.6] dr=1.22 t=0.0ps kin=1.47 pot=21.57 Rg=11.409 SPS=14983
bl=0 pos[1]=[5.0 15.6 -8.3] dr=1.22 t=0.0ps kin=1.46 pot=21.58 Rg=11.411 SPS=14280
bl=0 pos[1]=[5.9 15.1 -7.9] dr=1.20 t=0.0ps kin=1.45 pot=21.61 Rg=11.401 SPS=15341
bl=0 pos[1]=[6.9 17.6 -7.4] dr=1.25 t=0.0ps kin=1.49 pot=21.57 Rg=11.396 SPS=15287
bl=0 pos[1]=[6.2 17.8 -6.1] dr=1.24 t=0.0ps kin=1.53 pot=21.58 Rg=11.349 SPS=14955
bl=0 pos[1]=[6.8 17.2 -5.6] dr=1.26 t=0.0ps kin=1.48 pot=21.58 Rg=11.303 SPS=14678
bl=0 pos[1]=[7.8 17.2 -5.7] dr=1.24 t=0.0ps kin=1.46 pot=21.61 Rg=11.282 SPS=13784
bl=0 pos[1]=[9.7 16.6 -4.1] dr=1.25 t=0.0ps kin=1.49 pot=21.61 Rg=11.204 SPS=15248
bl=0 pos[1]=[9.8 16.7 -2.8] dr=1.26 t=0.0ps kin=1.55 pot=21.57 Rg=11.155 SPS=14220
bl=0 pos[1]=[9.7 17.2 -3.2] dr=1.25 t=0.0ps kin=1.47 pot=21.60 Rg=11.103 SPS=14423
bl=0 pos[1]=[8.6 17.9 -2.3] dr=1.27 t=0.0ps kin=1.48 pot=21.60 Rg=11.086 SPS=14795
bl=0 pos[1]=[6.9 18.8 -2.6] dr=1.26 t=0.0ps kin=1.51 pot=21.55 Rg=11.090 SPS=14605
bl=0 pos[1]=[5.6 18.1 -2.7] dr=1.24 t=0.0ps kin=1.46 pot=21.62 Rg=11.104 SPS=14324
bl=0 pos[1]=[5.2 18.4 -2.3] dr=1.25 t=0.0ps kin=1.49 pot=21.59 Rg=11.139 SPS=13517
bl=0 pos[1]=[5.6 17.0 -2.8] dr=1.26 t=0.0ps kin=1.51 pot=21.55 Rg=11.183 SPS=15119
bl=0 pos[1]=[4.6 15.4 -1.6] dr=1.23 t=0.0ps kin=1.51 pot=21.57 Rg=11.195 SPS=14360
bl=0 pos[1]=[3.8 15.3 -1.6] dr=1.24 t=0.0ps kin=1.48 pot=21.55 Rg=11.167 SPS=15966
bl=0 pos[1]=[4.7 14.6 -1.9] dr=1.26 t=0.0ps kin=1.45 pot=21.56 Rg=11.132 SPS=15034
bl=0 pos[1]=[5.0 14.4 -2.6] dr=1.25 t=0.0ps kin=1.47 pot=21.58 Rg=11.110 SPS=14345
bl=0 pos[1]=[6.2 15.3 -2.2] dr=1.22 t=0.0ps kin=1.50 pot=21.59 Rg=11.096 SPS=14920
bl=0 pos[1]=[5.3 16.2 -1.3] dr=1.24 t=0.0ps kin=1.53 pot=21.61 Rg=11.095 SPS=14239
bl=0 pos[1]=[6.3 16.6 -1.6] dr=1.28 t=0.0ps kin=1.54 pot=21.62 Rg=11.075 SPS=14562
bl=0 pos[1]=[6.6 18.5 -1.6] dr=1.26 t=0.0ps kin=1.53 pot=21.61 Rg=11.034 SPS=13436
bl=0 pos[1]=[6.6 19.7 -2.9] dr=1.22 t=0.0ps kin=1.50 pot=21.66 Rg=11.010 SPS=14621
bl=0 pos[1]=[7.1 18.9 -5.9] dr=1.25 t=0.0ps kin=1.46 pot=21.61 Rg=10.979 SPS=13823
bl=0 pos[1]=[6.4 16.5 -6.9] dr=1.24 t=0.0ps kin=1.49 pot=21.56 Rg=10.956 SPS=14835
bl=0 pos[1]=[6.4 15.2 -6.5] dr=1.26 t=0.0ps kin=1.50 pot=21.57 Rg=10.915 SPS=14853
bl=0 pos[1]=[4.0 14.4 -4.8] dr=1.26 t=0.0ps kin=1.51 pot=21.60 Rg=10.898 SPS=13865
bl=0 pos[1]=[2.0 13.9 -4.0] dr=1.25 t=0.0ps kin=1.46 pot=21.61 Rg=10.884 SPS=11310
bl=0 pos[1]=[2.1 13.7 -4.3] dr=1.22 t=0.0ps kin=1.48 pot=21.52 Rg=10.878 SPS=14068
bl=0 pos[1]=[4.3 12.2 -4.6] dr=1.26 t=0.0ps kin=1.47 pot=21.62 Rg=10.836 SPS=14063
bl=0 pos[1]=[6.1 12.2 -3.7] dr=1.23 t=0.0ps kin=1.53 pot=21.57 Rg=10.876 SPS=14536
bl=0 pos[1]=[7.3 12.0 -2.6] dr=1.24 t=0.0ps kin=1.46 pot=21.60 Rg=10.955 SPS=14182
bl=0 pos[1]=[8.0 11.1 -2.3] dr=1.22 t=0.0ps kin=1.52 pot=21.58 Rg=11.024 SPS=14548
bl=0 pos[1]=[8.4 11.3 -2.7] dr=1.25 t=0.0ps kin=1.49 pot=21.60 Rg=11.124 SPS=14879
bl=0 pos[1]=[8.9 11.3 -2.2] dr=1.23 t=0.0ps kin=1.53 pot=21.60 Rg=11.177 SPS=14002
bl=0 pos[1]=[9.6 10.5 -3.2] dr=1.22 t=0.0ps kin=1.52 pot=21.60 Rg=11.176 SPS=14717
bl=0 pos[1]=[9.3 10.0 -2.8] dr=1.25 t=0.0ps kin=1.50 pot=21.60 Rg=11.188 SPS=14164
bl=0 pos[1]=[10.1 10.2 -2.6] dr=1.25 t=0.0ps kin=1.51 pot=21.56 Rg=11.209 SPS=14505
bl=0 pos[1]=[8.2 10.0 -3.0] dr=1.26 t=0.0ps kin=1.54 pot=21.59 Rg=11.234 SPS=14060
bl=0 pos[1]=[6.4 11.3 -4.3] dr=1.29 t=0.0ps kin=1.48 pot=21.60 Rg=11.244 SPS=14579
bl=0 pos[1]=[5.9 10.6 -5.0] dr=1.26 t=0.0ps kin=1.54 pot=21.57 Rg=11.286 SPS=14090
bl=0 pos[1]=[5.1 10.1 -6.1] dr=1.27 t=0.0ps kin=1.48 pot=21.59 Rg=11.328 SPS=14923
bl=0 pos[1]=[4.9 9.1 -8.0] dr=1.25 t=0.0ps kin=1.52 pot=21.60 Rg=11.321 SPS=15368
bl=0 pos[1]=[5.7 9.3 -8.6] dr=1.25 t=0.0ps kin=1.51 pot=21.62 Rg=11.362 SPS=13698
bl=0 pos[1]=[5.3 9.6 -9.0] dr=1.26 t=0.0ps kin=1.56 pot=21.59 Rg=11.408 SPS=14979

Some details about the output for each performed block:

bl=0 is the number of each simulated block. In this case we set increment=False, so the number of steps is not being accounted. pos[1]=[X,Y,Z] is the spatial position for the first bead. dr=1.26 shows the average position displacement in each block (in units of sigma). t=0 is the total time (?). In this case we set increment=False, so the time is not being increased. kin=1.5 is the kinect energy of the system. pot=19.90 is the potential energy of the system. RG=7.654 is the radius of gyration in the end of each block. SPS=12312 is the number of steps per second of each block. A measure of how fast the computations are being performed.

Check the convergence of the radius of gyration:

plt.plot(rg)
plt.xlabel('block number')
plt.ylabel(r'radius of gyration ($\sigma$)')

Text(0, 0.5, 'radius of gyration ($\\sigma$)')

The next step is to remove the restraining force in order to run the sampling simulation

sim.removeFlatBottomHarmonic()

Initiate the Adam optimization object. In this tutorial, this object is named “opt”

opt = AdamTraining(mu=3.22, rc = 1.78, eta=0.01, it=1)

Let’s take a look on Adam initialization parameters:

mu and rc are calculated to fit the experimental Hi-C map. These specific values are for GM12878 (see Di Pierro et al., 2016). eta is the learning rate (default value is 0.01, works well to train the pairwise interactions based on the experimental Hi-C maps) it is the iteration step.

Load the contact matrix file to the object opt

opt.getPars(HiC="input/chr10_100k.dense")

The optimization step uses the contact map from the simulation to update the values in the force field. Therefore, the next step is to perform a long simulation with enough statistics to obtain a converged simulated contact map.

In order to get a good inversion calculation, it is important to have around \(1\times10^5\) frames from a set of different simulations. For example, \(20\) replicas of \(5000\) saved frames.

This can take some time, so in this tutorial we will use just 1 replica of \(5000\) frames saved every \(1000\) steps.

\(block = 1\times10^3\) \(n\_blocks = 5000\)

block = 1000
n_blocks = 50

for _ in range(n_blocks):
    sim.runSimBlock(block, increment=True) #perform 1 block of the simulation
    state = sim.getPositions() #get the positions for each block
    opt.probCalc(state) #calculate and store the contact map for each block

bl=1 pos[1]=[5.5 10.3 -7.9] dr=1.88 t=0.0ps kin=1.46 pot=21.58 Rg=11.463 SPS=11324
bl=2 pos[1]=[6.7 11.7 -9.4] dr=1.88 t=0.0ps kin=1.50 pot=21.53 Rg=11.610 SPS=11625
bl=3 pos[1]=[6.7 11.2 -9.9] dr=1.89 t=0.0ps kin=1.50 pot=21.51 Rg=11.857 SPS=13000
bl=4 pos[1]=[6.1 8.8 -10.1] dr=1.93 t=0.0ps kin=1.52 pot=21.51 Rg=12.107 SPS=15198
bl=5 pos[1]=[8.1 8.0 -8.6] dr=1.95 t=0.0ps kin=1.54 pot=21.53 Rg=12.210 SPS=15751
bl=6 pos[1]=[8.0 8.4 -4.9] dr=1.92 t=0.0ps kin=1.52 pot=21.52 Rg=12.351 SPS=15347
bl=7 pos[1]=[4.6 8.0 -7.0] dr=1.89 t=0.0ps kin=1.48 pot=21.56 Rg=12.410 SPS=15001
bl=8 pos[1]=[5.7 11.8 -4.5] dr=1.91 t=0.0ps kin=1.49 pot=21.59 Rg=12.485 SPS=14968
bl=9 pos[1]=[6.4 13.6 -5.2] dr=1.92 t=0.0ps kin=1.52 pot=21.53 Rg=12.653 SPS=14891
bl=10 pos[1]=[4.5 10.4 -4.5] dr=1.91 t=0.0ps kin=1.55 pot=21.49 Rg=12.824 SPS=15702
bl=11 pos[1]=[6.1 9.3 -7.0] dr=1.90 t=0.0ps kin=1.52 pot=21.52 Rg=12.994 SPS=15009
bl=12 pos[1]=[5.9 13.4 -9.1] dr=1.90 t=0.0ps kin=1.53 pot=21.58 Rg=13.158 SPS=14747
bl=13 pos[1]=[5.5 14.2 -9.4] dr=1.93 t=0.0ps kin=1.53 pot=21.49 Rg=13.266 SPS=14975
bl=14 pos[1]=[5.3 10.1 -9.4] dr=1.91 t=0.0ps kin=1.48 pot=21.55 Rg=13.334 SPS=14952
bl=15 pos[1]=[3.8 8.2 -8.7] dr=1.91 t=0.0ps kin=1.53 pot=21.58 Rg=13.469 SPS=14917
bl=16 pos[1]=[3.3 7.3 -8.3] dr=1.83 t=0.0ps kin=1.52 pot=21.49 Rg=13.521 SPS=15231
bl=17 pos[1]=[3.3 6.4 -9.3] dr=1.94 t=0.0ps kin=1.44 pot=21.55 Rg=13.493 SPS=14712
bl=18 pos[1]=[4.3 6.2 -11.1] dr=1.92 t=0.0ps kin=1.49 pot=21.52 Rg=13.548 SPS=15783
bl=19 pos[1]=[7.4 6.9 -12.6] dr=1.87 t=0.0ps kin=1.44 pot=21.56 Rg=13.681 SPS=15239
bl=20 pos[1]=[7.0 8.8 -12.4] dr=1.85 t=0.0ps kin=1.48 pot=21.52 Rg=13.740 SPS=14830
bl=21 pos[1]=[5.1 9.8 -11.0] dr=1.94 t=0.0ps kin=1.50 pot=21.55 Rg=13.804 SPS=15532
bl=22 pos[1]=[4.0 8.0 -10.9] dr=1.89 t=0.0ps kin=1.52 pot=21.52 Rg=13.964 SPS=15286
bl=23 pos[1]=[3.8 10.5 -9.6] dr=1.87 t=0.0ps kin=1.52 pot=21.57 Rg=14.019 SPS=15235
bl=24 pos[1]=[3.5 11.2 -10.9] dr=1.90 t=0.0ps kin=1.43 pot=21.54 Rg=14.053 SPS=16519
bl=25 pos[1]=[6.2 10.7 -11.7] dr=1.93 t=0.0ps kin=1.53 pot=21.56 Rg=14.008 SPS=16787
bl=26 pos[1]=[6.6 7.7 -12.1] dr=1.97 t=0.0ps kin=1.51 pot=21.53 Rg=14.083 SPS=16116
bl=27 pos[1]=[4.4 8.6 -10.6] dr=1.96 t=0.0ps kin=1.49 pot=21.58 Rg=14.146 SPS=16359
bl=28 pos[1]=[4.8 12.2 -10.4] dr=1.95 t=0.0ps kin=1.46 pot=21.55 Rg=14.174 SPS=16202
bl=29 pos[1]=[5.2 12.6 -9.8] dr=1.97 t=0.0ps kin=1.51 pot=21.54 Rg=14.196 SPS=16520
bl=30 pos[1]=[3.7 12.9 -9.9] dr=1.91 t=0.0ps kin=1.45 pot=21.52 Rg=14.144 SPS=15429
bl=31 pos[1]=[5.8 9.6 -7.4] dr=1.88 t=0.0ps kin=1.49 pot=21.54 Rg=14.221 SPS=15646
bl=32 pos[1]=[6.3 10.4 -6.3] dr=1.88 t=0.0ps kin=1.48 pot=21.52 Rg=14.266 SPS=15488
bl=33 pos[1]=[5.6 6.7 -7.0] dr=1.91 t=0.0ps kin=1.47 pot=21.50 Rg=14.302 SPS=16323
bl=34 pos[1]=[6.1 6.5 -8.3] dr=1.94 t=0.0ps kin=1.54 pot=21.54 Rg=14.444 SPS=15755
bl=35 pos[1]=[7.9 7.1 -7.6] dr=1.90 t=0.0ps kin=1.53 pot=21.51 Rg=14.456 SPS=15281
bl=36 pos[1]=[6.2 5.0 -10.2] dr=1.96 t=0.0ps kin=1.52 pot=21.54 Rg=14.557 SPS=15919
bl=37 pos[1]=[5.0 7.4 -10.4] dr=1.96 t=0.0ps kin=1.49 pot=21.52 Rg=14.650 SPS=15861
bl=38 pos[1]=[7.3 9.5 -9.1] dr=1.83 t=0.0ps kin=1.54 pot=21.56 Rg=14.604 SPS=16106
bl=39 pos[1]=[8.0 8.5 -10.2] dr=1.98 t=0.0ps kin=1.51 pot=21.56 Rg=14.597 SPS=16383
bl=40 pos[1]=[6.7 7.6 -8.3] dr=1.91 t=0.0ps kin=1.47 pot=21.53 Rg=14.613 SPS=16118
bl=41 pos[1]=[6.8 7.1 -7.4] dr=1.94 t=0.0ps kin=1.53 pot=21.53 Rg=14.552 SPS=15405
bl=42 pos[1]=[7.9 7.4 -5.8] dr=1.96 t=0.0ps kin=1.54 pot=21.54 Rg=14.560 SPS=15734
bl=43 pos[1]=[8.6 6.8 -6.3] dr=1.94 t=0.0ps kin=1.56 pot=21.56 Rg=14.638 SPS=16091
bl=44 pos[1]=[10.7 8.1 -7.1] dr=1.92 t=0.0ps kin=1.52 pot=21.52 Rg=14.659 SPS=15797
bl=45 pos[1]=[13.2 7.3 -4.7] dr=1.94 t=0.0ps kin=1.50 pot=21.53 Rg=14.690 SPS=15806
bl=46 pos[1]=[12.8 5.8 -6.6] dr=1.91 t=0.0ps kin=1.50 pot=21.56 Rg=14.725 SPS=15625
bl=47 pos[1]=[10.2 4.3 -7.1] dr=1.93 t=0.0ps kin=1.46 pot=21.50 Rg=14.795 SPS=16080
bl=48 pos[1]=[9.6 8.2 -9.8] dr=1.89 t=0.0ps kin=1.46 pot=21.47 Rg=14.775 SPS=15822
bl=49 pos[1]=[7.1 5.9 -9.0] dr=1.94 t=0.0ps kin=1.51 pot=21.48 Rg=14.797 SPS=15511

Statistics for the simulation opt_chr10_100K, number of particles: 1356,  number of chains: 1

Statistics for particle position
     mean position is:  [-1.23957781  1.7779258   2.78054735]   Rg =  14.796688
     median bond size is  0.9680976445921962
     three shortest/longest (<10)/ bonds are  [0.86936723 0.88400501 0.88756599]    [1.09114228 1.09780597 1.10176906]
     95 percentile of distance to center is:    21.431837046446862
     density of closest 95% monomers is:    0.031240393096327784
     density of the core monomers is:    0.07548378809911002
     min/median/mean/max coordinates are:
     x: -21.97, -1.27, -1.24, 20.78
     y: -18.23, 1.69, 1.78, 28.29
     z: -12.99, 2.12, 2.78, 20.77

Statistics for velocities:
     mean kinetic energy is:  1.512121087969534 should be: 1.5
     fastest particles are (in kT):  [7.25928234 7.28251669 7.57800841 7.84269129 8.89079614]

Statistics for the system:
     Forces are:  ['FENEBond', 'AngleForce', 'RepulsiveSoftCore', 'CustomTypes']
     Number of exceptions:   1355

Potential Energy Ep =  21.483418602507374
bl=50 pos[1]=[5.9 7.4 -7.8] dr=1.97 t=0.0ps kin=1.50 pot=21.44 Rg=14.789 SPS=16144

In the end of each simulation replica we need to save some important values required to calculate the Adam step.

We save these values using the H5 compressed files because of the efficiency writing/reading them.

Note: attention to this step. We have two files for each replica for each iteration step. Be organized!

replica=1

with h5py.File(sim.folder + "/Pi_" + str(replica)+".h5", 'w') as hf:
    hf.create_dataset("Pi",  data=opt.Pi)

np.savetxt(sim.folder + "/Nframes_" + str(replica), [opt.NFrames])

The first part of the optimization is finished. Inside the output folder, for each iteration, there are these 2 files used in next step:

Nframes_1 Pi_1.h5

%%bash
ls iteration_1/

Nframes_1
opt_chr10_100K_0_block0.pdb
Pi_1.h5
platform_info.dat
probdist_1

The second part is the inversion. It is quite simple, just feed the optmization object with the files from all replicas and make the inversion to get the updated force field file.

opt2 = AdamTraining(mu=3.22, rc = 1.78, eta=0.01, it=1)
opt2.getPars(HiC="input/chr10_100k.dense")

Nreplicas = 1
path = "iteration_1"
for rep in range(1,Nreplicas+1):
    print("reading files from replica {}".format(rep))
    with h5py.File(path + "/Pi_" + str(rep)+".h5", 'r') as hf:
            opt2.Pi += hf['Pi'][:]
            opt2.NFrames += int(np.loadtxt(path + "/Nframes_" + str(rep)))

reading files from replica 1

With the contact probabilities from with all replicas, we calculate the inversion and get the updated parameters for the force field.

To calculate the new FF, we need to load the old one.

lamb_new = opt2.getLamb(Lambdas="input/lambda_1")

lamb_new

	Bead1	Bead2	Bead3	Bead4	Bead5	Bead6	Bead7	Bead8	Bead9	Bead10	...	Bead1347	Bead1348	Bead1349	Bead1350	Bead1351	Bead1352	Bead1353	Bead1354	Bead1355	Bead1356
0	-0.00999	-0.01000	0.01000	-0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000
1	-0.01000	-0.00999	-0.01000	0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000
2	0.01000	-0.01000	-0.00999	-0.01000	0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000
3	-0.01000	0.01000	-0.01000	-0.00999	-0.01000	0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000
4	-0.01000	-0.01000	0.01000	-0.01000	-0.00999	-0.01	0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
1351	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	0.01	-0.01	-0.00999	-0.01000	0.01000	-0.01000	-0.01000
1352	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	0.01	-0.01000	-0.00999	-0.01000	0.01000	-0.01000
1353	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	0.01000	-0.01000	-0.00999	-0.01000	0.01000
1354	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	0.01000	-0.01000	-0.00999	-0.01000
1355	-0.01000	-0.01000	-0.01000	-0.01000	-0.01000	-0.01	-0.01	-0.01	-0.01	-0.01	...	-0.01	-0.01	-0.01	-0.01	-0.01	-0.01000	-0.01000	0.01000	-0.01000	-0.00999

1356 rows × 1356 columns

Save the new force field to a file. To keep the organization, we use the same name followed by the number of the new iteration:

lamb_new.to_csv("input/lambda_2", index=False)

We can plot the parameters of the force field to check the changes between the 2 steps:

ff_old = pd.read_csv("input/lambda_1")
ff_new = pd.read_csv("input/lambda_2")

fig, axes = plt.subplots(nrows=1, ncols=2, figsize=(8,4))

f1 = axes[0].matshow(ff_old.values, vmin=-0.1, vmax=0.1, cmap="bwr")
axes[0].set_title('lambda_1', loc='center')
fig.colorbar(f1, orientation='horizontal', ax=axes[0], shrink=0.9, fraction=0.046, pad=0.04)

f2 = axes[1].matshow(ff_new.values, vmin=-0.1, vmax=0.1, cmap="bwr")
axes[1].set_title('lambda_2', loc='center')
fig.colorbar(f2, orientation='horizontal', ax=axes[1], shrink=0.9, fraction=0.046, pad=0.04)

<matplotlib.colorbar.Colorbar at 0x7f76c41f9a60>

We can also check the error between experimental and simulation Hi-C maps:

print(opt2.error)

0.9527989958833853

Save/plot the simulation Hi-C and compare to the experimental one

#Save
np.savetxt("iteration_1/probdist_1", opt2.phi_sim)

#Plot
fig, axes = plt.subplots(nrows=1, ncols=2, figsize=(8,4))

f1 = axes[0].matshow(opt2.phi_exp,norm=mpl.colors.LogNorm(vmin=0.0001, vmax=opt2.phi_exp.max()),cmap="Reds")
axes[0].set_title('Experimental Hi-C', loc='center')
fig.colorbar(f1, orientation='horizontal', ax=axes[0], shrink=0.9, fraction=0.046, pad=0.04)

f2 = axes[1].matshow(opt2.phi_sim,norm=mpl.colors.LogNorm(vmin=0.0001, vmax=opt2.phi_exp.max()),cmap="Reds")
axes[1].set_title('Simulation Hi-C', loc='center')
fig.colorbar(f2, orientation='horizontal', ax=axes[1], shrink=0.9, fraction=0.046, pad=0.04)

<matplotlib.colorbar.Colorbar at 0x7f76c4103d30>

As one can see, the simulation Hi-C map is not converged. This is due to two main reasons: 1 - We just performed a single simulation instead of multiple replicas; 2 - The force field remains almost with zero interactions due to the number of optimization steps.

Redo these steps using the new lambda file (lambda_2) as input for customTypes potential in the next iteration.

Fixing values on the force field

If we inspect again the experimental Hi-C map you can observe that some regions have a white band. These gaps arise from difficulty of sequencing those regions, often highly repetitive regions. The centromere of the chromosomes is one of them.

filename = 'input/chr10_100k.dense'
hic_file = np.loadtxt(filename)

r=np.triu(hic_file, k=1)
r[np.isnan(r)]= 0.0
r = normalize(r, axis=1, norm='max')
rd = np.transpose(r)
r=r+rd + np.diag(np.ones(len(r)))
plt.matshow(r,norm=mpl.colors.LogNorm(vmin=0.0001, vmax=r.max()),cmap="Reds")
plt.colorbar()

<matplotlib.colorbar.Colorbar at 0x7f81b07e71f0>

So, when we perform the Adam optimization, the parameters regarding the interactions of these regions will not be trained because of the lack of data.

One way to solve this problem is to define some interactions initially and then keep these fixed over iterations.

Find the index position of each white band

#get the first diagonal (first neighbor)
diag1 = np.diag(r, k=1)

#find the index where the first neighbor is equal to zero
centromere_index = np.where(diag1 == 0)[0]
print(centromere_index)

[391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408
 409 410 411 412 413 414 415 416 417 418 419 420 421 422 463 464 466 467
 468 477 478 480 481 482 490]

Now we can create the new initial force field parameters’ matrix with a specific value for these interactions. Here we chose \(-0.3\) for centromere-centromere interactions and \(-0.2\) for centromere-other beads.

#read the sequence file
seq = np.loadtxt("input/seq_chr10", dtype=str)[:,1]
pol_size = len(seq)
#create a NxN matrix fill with zeros, where N is the chromossome lenght
data = np.zeros((pol_size,pol_size))

#Set values for centromere regions and NA regions
cent = [391,422]
NAs = [463, 464, 466, 467, 468, 477, 478, 480, 481, 482, 490]
data[cent[0]:cent[1],:] = -0.2
data[:,cent[0]:cent[1]] = -0.2
data[cent[0]:cent[1],cent[0]:cent[1]] = -0.3
for na in NAs:
    data[na,:] = -0.2
    data[:,na] = -0.2

#transform it in a Pandas dataframe using the sequence names as header
lamb = pd.DataFrame(data, columns=seq)
#save to a csv file format in the  "input" folder
lamb.to_csv("input/lambda_fixed_1", index=False)
lamb

	Bead1	Bead2	Bead3	Bead4	Bead5	Bead6	Bead7	Bead8	Bead9	Bead10	...	Bead1347	Bead1348	Bead1349	Bead1350	Bead1351	Bead1352	Bead1353	Bead1354	Bead1355	Bead1356
0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
2	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
3	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
4	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
1351	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1352	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1353	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1354	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
1355	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	...	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0

1356 rows × 1356 columns

plt.matshow(lamb.values, vmin=-0.4, vmax=0.4, cmap="bwr")
plt.colorbar()

<matplotlib.colorbar.Colorbar at 0x7f81b0a783d0>

Save the information about the index of the gaps.

When calling the Adam optimization step, we can define fixed points. We provide a list of pair of interactions that will remain unchanged throughout the optimization procedure.

import itertools
#list of all beads
N = np.arange(opt2.expHiC.shape[0])

#use the itertools to create all pair of iteractions
fix_points = list(itertools.product(centromere_index,N))

#Use the following line to update the parameters maintaining the desired interactions fixed
#lamb_new = opt2.getLamb(Lambdas="input/lambda_1", fixedPoints=fix_points)