User Guide: 1 Pre-processing

Introduction

At the moment this vignette has the double purpose of keeping track of my progress and testing the performance of the functions in the package. My aim is to analyse a large set of high frequency time series measurements of sunlight at seven different wavebands done with broadband sensors over one summer.

The data set will be large with seven variables and a few hundred millions of time points. During daytime measurements were acquired in bursts of 18000 time points in 15 min every 30 min.

The effect of clouds is intermittent and varies in intensity gradually. Using quantiles, mean, MAD and variance on the irradiance, and on the running differences is one possible approach. However, a recent additional approach is based on detecting change points in the running differences as a basis to quantify the duration and amplitude of cloud/shade flecks.

The data are in large text files saved with the software PC400 from Campbell Scientific, the supplier of the data logger used to collect the data set.

The first step is to read-in the data and split it into “chunks” of data from 18000 consecutive observations at the gaps between the “bursts” of observations, discarding any incomplete sets and any data measured at a lower frequency. At this stage running differences for time stamps need to be computed, and running differences for the measured variables can also be added at the same time.

The smaller values in the running differences are most likely the result of measurement “noise” and could be “cleaned” before the next steps. The differences can be used directly or discretized based on their sign.

Current state

The reading of the data from the 3.7GB text file into a tibble is the slowest step. Other computations are fast even with the functions are coded fully in R.

Set up

library(dplyr)
#> 
#> Attaching package: 'dplyr'
#> The following objects are masked from 'package:stats':
#> 
#>     filter, lag
#> The following objects are masked from 'package:base':
#> 
#>     intersect, setdiff, setequal, union
library(photobiology)
#> Loading required package: SunCalcMeeus
#> Documentation at https://docs.r4photobiology.info/
library(photobiologyInOut)
library(photobiologyFlecks)
library(ggplot2)
set_theme(theme_bw())

The data

We use data included in the package as a data frame. These data were originally saved as a CSV file by Campbell Scientific’s PC400 program. The data were acquired with a CR6 data logger.

head(three_chunks.tb)
#> # A tibble: 6 × 2
#>   time                Q_PAR
#>   <dttm>              <dbl>
#> 1 2025-07-05 09:44:35  762.
#> 2 2025-07-05 09:44:35  762.
#> 3 2025-07-05 09:44:35  762.
#> 4 2025-07-05 09:44:35  762.
#> 5 2025-07-05 09:44:35  762.
#> 6 2025-07-05 09:44:35  762.
plot(three_chunks.tb, cex = 0.1)

Split data into chunks

The time series data of sunlight have been measured in 15-min-long bursts separated by 15 min with no data acquisition. Data from each burst will be analysed separately. For this the data frame has to be split at the gaps between bursts.

# keep all qty columns, add differences
chunks.ls <-
  split_chunks(three_chunks.tb,
               time.name = "time",
               qty.name = NULL,
               chunk.min.time = 0.051,
               chunk.min.rows = 1.8e4,
               verbose = FALSE)

Some checks

length(chunks.ls) == length(unique(names(chunks.ls)))
#> [1] TRUE
nchunks_found <- length(chunks.ls)
str(chunks.ls, give.attr = FALSE)
#> List of 3
#>  $ 2025-07-05 10:00:00.05: tibble [18,000 × 4] (S3: tbl_df/tbl/data.frame)
#>   ..$ time      : POSIXct[1:18000], format: "2025-07-05 10:00:00" "2025-07-05 10:00:00" ...
#>   ..$ Q_PAR     : num [1:18000] 1396 1396 1396 1396 1396 ...
#>   ..$ time.diff : num [1:18000] 0.05 0.05 0.05 0.05 0.05 ...
#>   ..$ Q_PAR.diff: num [1:18000] 12601.57 5.08 -3.96 5.76 -2.3 ...
#>  $ 2025-07-05 10:30:00.05: tibble [18,000 × 4] (S3: tbl_df/tbl/data.frame)
#>   ..$ time      : POSIXct[1:18000], format: "2025-07-05 10:30:00" "2025-07-05 10:30:00" ...
#>   ..$ Q_PAR     : num [1:18000] 1797 1802 1805 1808 1810 ...
#>   ..$ time.diff : num [1:18000] 0.05 0.05 0.05 0.05 0.05 ...
#>   ..$ Q_PAR.diff: num [1:18000] 9414.9 19.8 95.4 64.8 56.5 ...
#>  $ 2025-07-05 11:00:00.05: tibble [18,000 × 4] (S3: tbl_df/tbl/data.frame)
#>   ..$ time      : POSIXct[1:18000], format: "2025-07-05 11:00:00" "2025-07-05 11:00:00" ...
#>   ..$ Q_PAR     : num [1:18000] 1137 1137 1137 1137 1137 ...
#>   ..$ time.diff : num [1:18000] 0.05 0.05 0.05 0.05 0.05 ...
#>   ..$ Q_PAR.diff: num [1:18000] -2573.7 -3.68 -1.02 0.18 -1.74 ...
plot(Q_PAR ~ time, data = chunks.ls[[1]], type = "l", lwd = 0.5)

plot(Q_PAR.diff ~ time, data = chunks.ls[[1]], type = "l", lwd = 0.5)

Grouping

Instead of splitting the data as in the section above, an indexing vector can be used to do computations in place with fast functions from R package ‘collapse’. Given the size of the data sets this could be worthwhile.

# generate a vector suitable for grouping
chunks.grp <-
  group_chunks(three_chunks.tb,
               time.name = "time",
               qty.name = NULL,
               chunk.min.time = 0.051,
               chunk.min.rows = 1.8e4,
               verbose = FALSE)
unique(chunks.grp)
#> [1] NA  2  3  4

Summaries

Statistics

Something seems wrong as I expect a 1 h time shift to be needed!

statistics.tb <- 
  summarize_chunks(chunks.ls, 
                   add.times = TRUE, 
                   tz = "EET",
                   time.shift = 0,
                   add.solar.times = TRUE,
                   geocode = data.frame(lon = 25.019212, lat = 60.226805),
                   verbose = FALSE)
#> Summarized 3 chunks
colnames(statistics.tb)
#> [1] "Q_PAR"      "time.diff"  "Q_PAR.diff" "parameter"  "chunk"     
#> [6] "date"       "time"       "solar.time"
length(unique(statistics.tb[["chunk"]])) == nchunks_found
#> [1] TRUE
length(unique(statistics.tb[["time"]])) == nchunks_found
#> [1] TRUE
statistics.tb
#> # A tibble: 30 × 8
#>       Q_PAR time.diff  Q_PAR.diff parameter chunk            date               
#>       <dbl>     <dbl>       <dbl> <fct>     <chr>            <dttm>             
#>  1  922.       0.0500  -269.      min       2025-07-05 10:0… 2025-07-05 00:00:00
#>  2 1006.       0.0500    -8.00    q.25      2025-07-05 10:0… 2025-07-05 00:00:00
#>  3 1135.       0.0500    -0.330   median    2025-07-05 10:0… 2025-07-05 00:00:00
#>  4 1514.       0.0500     7.76    q.75      2025-07-05 10:0… 2025-07-05 00:00:00
#>  5 2218.     900.     12602.      max       2025-07-05 10:0… 2025-07-05 00:00:00
#>  6 1265.       0.1        0.623   mean      2025-07-05 10:0… 2025-07-05 00:00:00
#>  7  317.       6.71     100.      sd        2025-07-05 10:0… 2025-07-05 00:00:00
#>  8  257.       0         11.7     mad       2025-07-05 10:0… 2025-07-05 00:00:00
#>  9    3.99     0.0149     0.00623 CVsd      2025-07-05 10:0… 2025-07-05 00:00:00
#> 10    0.226    0        -35.4     CVmad     2025-07-05 10:0… 2025-07-05 00:00:00
#> # ℹ 20 more rows
#> # ℹ 2 more variables: time <dttm>, solar.time <dbl>

Locating change points

Using diff()

Running differences after slight denoising can be used to detect changes in the sign of the local slope, the peaks and valleys of a time series. In other words, the flecks are identified based on local maxima and minima.

We use relative.threshold = 0.20 \(\mu mol\,m^{-2}\,s^{-2}\) as the minimum rate of change in irradiance per second to be considered relevant and not discarded as noise.

  chunks_dn.ls <- 
    denoise_chunks(chunks.ls, 
                   time.name = "time",
                   relative.threshold = 0.0067, 
                   add.signs = TRUE)
#> Denoised 3 chunks
plot(Q_PAR.diff ~ time, data = chunks_dn.ls[[1]], type = "l", lwd = 0.5)

In a continuous time series without noise, the derivative is equal to zero at the point where the slope changes sign. In a time series measured at discrete time points, a change in sign must be used as well as a zero value to locate maxima and minima, as the change in slope is most likely to land in-between two time points.

Once de-noised, changes in slope can be detected using sign() to discretize the data into values indicating decrease (-1) , no-change (0) and increase (+1).

Considering each of the three bursts of PAR measurements, we find many time intervals with a detectable change in irradiance.

# time steps with no change
sum(abs(chunks_dn.ls[[1]]$Q_PAR.diff) < 0.01, na.rm = TRUE)
#> [1] 11471
# increasing time steps
sum(chunks_dn.ls[[1]]$Q_PAR.diff > 0.01, na.rm = TRUE)
#> [1] 3234
# decreasing time steps
sum(chunks_dn.ls[[1]]$Q_PAR.diff > 0.01, na.rm = TRUE)
#> [1] 3234
# time steps with no change
sum(abs(chunks_dn.ls[[2]]$Q_PAR.diff) < 0.01, na.rm = TRUE)
#> [1] 13144
# increasing time steps
sum(chunks_dn.ls[[2]]$Q_PAR.diff > 0.01, na.rm = TRUE)
#> [1] 2212
# decreasing time steps
sum(chunks_dn.ls[[2]]$Q_PAR.diff > 0.01, na.rm = TRUE)
#> [1] 2212
# time steps with no change
sum(abs(chunks_dn.ls[[3]]$Q_PAR.diff) < 0.01, na.rm = TRUE)
#> [1] 13990
# increasing time steps
sum(chunks_dn.ls[[3]]$Q_PAR.diff > 0.01, na.rm = TRUE)
#> [1] 1979
# decreasing time steps
sum(chunks_dn.ls[[3]]$Q_PAR.diff > 0.01, na.rm = TRUE)
#> [1] 1979

Subsequently, changes in slope can be located using rle(). In this example all differences smaller than 1 % of the mean irradiance per second were discarded as noise.

sign_changes.ls <- list()
for (i in names(chunks_dn.ls)) {
  sign_changes.ls[[i]] <- 
    list(rle = rle(chunks_dn.ls[[i]][["Q_PAR.sign"]]))
}
length(sign_changes.ls)
#> [1] 3

Check that sum of lengths equal length of burst.

sum(sign_changes.ls[[1]][["rle"]][["lengths"]])
#> [1] 18000
sum(sign_changes.ls[[2]][["rle"]][["lengths"]])
#> [1] 18000
sum(sign_changes.ls[[3]][["rle"]][["lengths"]])
#> [1] 18000
sign_changes.ls
#> $`2025-07-05 10:00:00.05`
#> $`2025-07-05 10:00:00.05`$rle
#> Run Length Encoding
#>   lengths: int [1:2125] 1 7 1 4 1 1 10 1 63 1 ...
#>   values : num [1:2125] 1 0 -1 0 -1 1 0 -1 0 -1 ...
#> 
#> 
#> $`2025-07-05 10:30:00.05`
#> $`2025-07-05 10:30:00.05`$rle
#> Run Length Encoding
#>   lengths: int [1:1429] 12 1 33 1 9 1 1 2 1 2 ...
#>   values : num [1:1429] 1 0 1 -1 1 0 -1 1 -1 1 ...
#> 
#> 
#> $`2025-07-05 11:00:00.05`
#> $`2025-07-05 11:00:00.05`$rle
#> Run Length Encoding
#>   lengths: int [1:1073] 1 1464 1 18 1 6 1 4 1 2 ...
#>   values : num [1:1073] -1 0 1 0 1 0 1 0 1 0 ...

Durand’s et al. algorithm

# keep all qty columns, do not add differences
chunks.ls <-
  split_chunks(three_chunks.tb,
               time.name = "time",
               qty.name = NULL,
               chunk.min.time = 0.051,
               chunk.min.rows = 1.8e4,
               add.diffs = FALSE,
               verbose = FALSE)
# extract one chunk, 15 min at 0.05 s
chunk <- chunks.ls[[1]]
range(chunk$time) |> as.numeric() |> diff() # seconds
#> [1] 1799.95

# time into chunk in seconds as numeric
chunk$tod <- as_tod(chunk$time, unit.out = "seconds")
chunk$chunk.time <- chunk$tod - chunk$tod[1]
range(chunk$chunk.time)
#> [1]    0.00 1799.95

Z <- find_zeros(time = chunk$chunk.time,
                var = chunk$Q_PAR)
#> Found 3505 zeros
#> You may want to try with a different 'zero.lim' value.

sun_flecks.df <- find_flecks(time = chunk$chunk.time,
                             var = chunk$Q_PAR,
                             minTime = 0.001, 
                             minAmp = 0, 
                             minPdiff = 0.00025,
                             asymmetry = 0,
                             shadeflecks = FALSE,
                             asmMethod = "max",
                             time.digits = 6,
                             var.digits = 5,
                             verbose = FALSE)
#> Found 241 sunflecks
sun_flecks.fig <-
  ggplot(data = chunk) +
  geom_line(aes(chunk.time, Q_PAR), linewidth = 0.2, na.rm = TRUE) +
  geom_point(aes(peakTime, peak), data = sun_flecks.df, 
             size = 0.5, colour = "coral2", na.rm = TRUE) +
  geom_point(aes(baselineTime1, baseline1), data = sun_flecks.df,
             size = 0.5, colour = "chartreuse3", na.rm = TRUE) +
  geom_point(aes(baselineTime2, baseline2), data = sun_flecks.df,
             size = 0.5, colour = "chartreuse3", na.rm = TRUE) +
  geom_text(aes(peakTime, peak, label = no), data = sun_flecks.df,
            colour = "coral2",
            position = position_nudge(y = 30), size = 2.7, na.rm = TRUE,
            fontface = "bold")

sun_flecks.fig


sun_flecks.fig + xlim(0, 100)


knitr::kable(sun_flecks.df)
no peakTime peak baseline1 baseline2 baseline amplitude relAmp percDiff duration irradGain relIrrad irradTot baselineTime1 baselineTime2 symmetry z timeInterval
1 0.500000 1395.1080 1394.73 1393.15 1393.940 1.17 1.1680 0.0008 0.150 0.0831991 1.663984e+00 278.8706 0.450000 0.599999 symmetric 101 NA
2 0.649999 1394.2670 1393.15 1386.90 1390.025 4.24 4.2420 0.0031 0.950 0.5214949 1.042991e+01 1390.5462 0.599999 1.549999 symmetric 131 0.149999
3 2.949997 1383.7230 1383.30 1382.65 1382.975 0.75 0.7480 0.0005 0.200 0.0674716 1.349433e+00 345.8105 2.899997 3.099997 symmetric 591 2.299998
4 9.049991 1448.7710 1444.78 1448.23 1446.505 2.27 2.2660 0.0016 0.500 0.4343587 8.687183e+00 796.0106 8.599992 9.099991 symmetric 1811 6.099994
5 9.349991 1448.0429 1447.60 1446.07 1446.835 1.21 1.2079 0.0008 0.250 0.0897367 1.794737e+00 434.1402 9.299991 9.549991 symmetric 1871 0.300000
6 14.249986 1357.3929 1356.79 1355.86 1356.325 1.07 1.0679 0.0008 0.300 0.0869277 1.738555e+00 474.7989 14.149987 14.449986 symmetric 2851 4.899995
7 14.499986 1357.5875 1355.86 1356.59 1356.225 1.36 1.3625 0.0010 0.100 0.0681177 1.362355e+00 203.5017 14.449986 14.549986 symmetric 2901 0.250000
8 14.649986 1358.0347 1356.59 1357.69 1357.140 0.89 0.8947 0.0007 0.150 0.0276111 5.522232e-01 271.4560 14.549986 14.699986 symmetric 2931 0.150000
9 17.699983 1410.0300 1409.32 1408.76 1409.040 0.99 0.9900 0.0007 0.300 0.1197255 2.394512e+00 493.2834 17.549983 17.849983 symmetric 3541 3.049997
10 20.199981 1376.7198 1376.30 1375.54 1375.920 0.80 0.7998 0.0006 0.200 0.0743228 1.486458e+00 344.0531 20.149981 20.349981 symmetric 4041 2.499998
11 20.499980 1376.2059 1375.54 1375.47 1375.505 0.70 0.7009 0.0005 0.250 0.0817155 1.634311e+00 412.7326 20.349981 20.599980 symmetric 4101 0.299999
12 20.749980 1377.1448 1375.47 1376.56 1376.015 1.13 1.1298 0.0008 0.200 0.0837132 1.674266e+00 344.0873 20.599980 20.799980 symmetric 4151 0.250000
13 25.949975 1571.5719 1379.65 1570.94 1475.295 96.28 96.2769 0.0653 4.850 7.8083765 1.561677e+02 7236.7504 21.149980 25.999975 symmetric 5191 5.199995
14 26.199975 1572.6850 1570.94 1572.29 1571.615 1.07 1.0700 0.0007 0.300 0.1685732 3.371467e+00 550.2336 25.999975 26.299975 symmetric 5241 0.250000
15 27.699974 1582.1490 1574.10 1581.73 1577.915 4.23 4.2340 0.0027 0.700 0.1758910 3.517824e+00 1183.6111 27.049974 27.749974 symmetric 5541 1.499999
16 28.149973 1584.0859 1583.68 1583.61 1583.645 0.44 0.4409 0.0003 0.150 0.0394673 7.893464e-01 316.7679 28.049973 28.199973 symmetric 5631 0.449999
17 28.299973 1584.6589 1583.61 1582.28 1582.945 1.71 1.7139 0.0011 0.300 0.2255785 4.511573e+00 554.2561 28.199973 28.499973 symmetric 5661 0.150000
18 30.099971 1565.8758 1564.90 1564.69 1564.795 1.08 1.0808 0.0007 0.200 0.1130781 2.261563e+00 391.3111 29.999971 30.199971 symmetric 6021 1.799998
19 30.249971 1565.6095 1564.69 1564.94 1564.815 0.79 0.7945 0.0005 0.100 0.0398313 7.966276e-01 234.7615 30.199971 30.299971 symmetric 6051 0.150000
20 30.549971 1566.9129 1564.94 1566.42 1565.680 1.23 1.2329 0.0008 0.300 0.1284574 2.569149e+00 548.1147 30.299971 30.599971 symmetric 6111 0.300000
21 30.849971 1569.5229 1567.68 1569.01 1568.345 1.18 1.1779 0.0008 0.150 0.1060803 2.121608e+00 313.7750 30.749971 30.899971 symmetric 6171 0.300000
22 31.149970 1573.5104 1569.01 1572.97 1570.990 2.52 2.5204 0.0016 0.350 0.2191727 4.383458e+00 628.6140 30.899971 31.249970 symmetric 6231 0.299999
23 37.349964 1689.4083 1579.93 1688.93 1634.430 54.98 54.9783 0.0336 5.450 12.0844612 2.416895e+02 9001.4479 31.949970 37.399964 symmetric 7471 6.199994
24 37.949964 1690.5056 1689.95 1689.12 1689.535 0.97 0.9706 0.0006 0.100 0.0485279 9.705581e-01 253.4785 37.899964 37.999964 symmetric 7591 0.600000
25 38.149964 1690.3306 1689.12 1689.62 1689.370 0.96 0.9606 0.0006 0.200 0.0891919 1.783840e+00 422.4320 37.999964 38.199964 symmetric 7631 0.200000
26 38.499963 1690.3456 1689.85 1689.67 1689.760 0.59 0.5856 0.0003 0.100 0.0291395 5.827910e-01 253.4933 38.449963 38.549963 symmetric 7701 0.349999
27 38.999963 1691.0578 1690.57 1690.37 1690.470 0.59 0.5878 0.0003 0.150 0.0423968 8.479378e-01 338.1361 38.899963 39.049963 symmetric 7801 0.500000
28 39.399962 1691.6506 1691.14 1691.04 1691.090 0.56 0.5606 0.0003 0.200 0.0479610 9.592210e-01 422.8201 39.349962 39.549962 symmetric 7881 0.399999
29 53.799949 1842.8341 1651.01 1717.47 1684.240 158.59 158.5941 0.0942 0.900 -0.9374014 -1.874805e+01 1599.0890 50.549952 51.449951 symmetric 10761 14.399987
30 54.449948 1850.9821 1850.13 1850.48 1850.305 0.68 0.6771 0.0004 0.100 0.0337669 6.753394e-01 277.5795 54.399948 54.499948 symmetric 10891 0.649999
31 54.749948 1852.6326 1850.10 1850.57 1850.335 2.30 2.2976 0.0012 0.300 0.2919885 5.839776e+00 647.9087 54.599948 54.899948 symmetric 10951 0.300000
32 54.999948 1851.6934 1850.57 1849.05 1849.810 1.88 1.8834 0.0010 0.250 0.2930156 5.860317e+00 555.2347 54.899948 55.149947 symmetric 11001 0.250000
33 55.199947 1849.6014 1849.05 1750.26 1799.655 49.95 49.9464 0.0278 2.150 22.4700181 4.494008e+02 3981.7111 55.149947 57.299945 symmetric 11041 0.199999
34 62.699940 1676.5665 1675.30 1675.54 1675.420 1.15 1.1465 0.0007 0.100 0.0573898 1.147798e+00 251.3700 62.649940 62.749940 symmetric 12541 7.499993
35 63.099940 1678.3378 1676.30 1677.67 1676.985 1.35 1.3528 0.0008 0.250 0.0613269 1.226539e+00 503.1572 62.899940 63.149940 symmetric 12621 0.400000
36 63.349940 1679.9930 1677.83 1679.30 1678.565 1.43 1.4280 0.0009 0.150 0.0582623 1.165248e+00 335.7710 63.249940 63.399940 symmetric 12671 0.250000
37 70.599933 1980.2699 1911.53 1979.31 1945.420 34.85 34.8499 0.0179 1.750 0.9632243 1.926451e+01 3502.7121 68.899934 70.649933 symmetric 14121 7.249993
38 70.749933 1983.8646 1979.31 1982.44 1980.875 2.99 2.9896 0.0015 0.150 0.2102303 4.204609e+00 396.3847 70.649933 70.799932 symmetric 14151 0.150000
39 71.849931 2008.2107 1982.44 2007.52 1994.980 13.23 13.2307 0.0066 1.200 5.6636379 1.132729e+02 2499.3870 70.799932 71.999931 symmetric 14371 1.099998
40 74.749929 1936.3684 1935.56 1934.84 1935.200 1.17 1.1685 0.0006 0.150 0.0976487 1.952976e+00 387.1369 74.649929 74.799929 symmetric 14951 2.899998
41 74.849929 1936.1651 1934.84 1935.04 1934.940 1.23 1.2251 0.0006 0.100 0.0614095 1.228191e+00 290.3017 74.799929 74.899929 symmetric 14971 0.100000
42 75.199928 1933.3473 1932.16 1930.62 1931.390 1.96 1.9573 0.0010 0.150 0.0932264 1.864529e+00 386.3714 75.149928 75.299928 symmetric 15041 0.349999
43 75.399928 1932.8081 1930.62 1929.38 1930.000 2.81 2.8081 0.0015 0.300 0.2668192 5.336389e+00 675.7667 75.299928 75.599928 symmetric 15081 0.200000
44 75.649928 1930.0600 1929.38 1923.93 1926.655 3.41 3.4050 0.0018 0.400 0.4034243 8.068493e+00 867.3962 75.599928 75.999928 symmetric 15131 0.250000
45 76.049927 1925.3020 1923.93 1922.73 1923.330 1.97 1.9720 0.0010 0.100 0.0985624 1.971250e+00 288.5979 75.999928 76.099927 symmetric 15211 0.399999
46 76.149927 1923.4949 1922.73 1821.08 1871.905 51.59 51.5899 0.0276 2.150 12.7981194 2.559626e+02 4130.9904 76.099927 78.249925 symmetric 15231 0.100000
47 83.399920 1704.7368 1699.84 1704.26 1702.050 2.69 2.6868 0.0016 0.550 0.6462556 1.292512e+01 1021.8750 82.899921 83.449920 symmetric 16681 7.249993
48 83.499920 1704.7752 1704.26 1703.97 1704.115 0.66 0.6602 0.0004 0.100 0.0329580 6.591613e-01 255.6501 83.449920 83.549920 symmetric 16701 0.100000
49 83.599920 1704.5920 1703.97 1289.74 1496.855 207.74 207.7370 0.1388 8.300 71.3259778 1.426521e+03 12570.0559 83.549920 91.849912 symmetric 16721 0.100000
50 184.899824 923.7587 923.51 923.50 923.505 0.25 0.2537 0.0003 0.100 0.0126300 2.525999e-01 138.5384 184.849824 184.949824 symmetric 36981 101.299904
51 188.149821 922.8660 922.62 922.62 922.620 0.25 0.2460 0.0003 0.200 0.0271998 5.439967e-01 230.6812 187.999821 188.199821 symmetric 37631 3.249997
52 202.499807 925.0884 924.80 924.83 924.815 0.27 0.2734 0.0003 0.200 0.0269297 5.385948e-01 231.2308 202.349807 202.549807 symmetric 40501 14.349986
53 258.049754 950.6303 950.27 950.39 950.330 0.30 0.3003 0.0003 0.100 0.0149178 2.983573e-01 142.5646 257.999754 258.099754 symmetric 51611 55.549947
54 339.349676 1272.7260 1271.39 1272.31 1271.850 0.88 0.8761 0.0007 0.250 0.0847634 1.695270e+00 381.6396 339.249676 339.499676 symmetric 67871 81.299922
55 339.549676 1272.7521 1272.31 1271.84 1272.075 0.68 0.6771 0.0005 0.150 0.0561274 1.122548e+00 254.4707 339.499676 339.649676 symmetric 67911 0.200000
56 339.699676 1272.6566 1271.84 1272.18 1272.010 0.65 0.6466 0.0005 0.150 0.0586855 1.173712e+00 254.4601 339.649676 339.799676 symmetric 67941 0.150000
57 343.399673 1322.7846 1321.64 1321.90 1321.770 1.01 1.0146 0.0008 0.350 0.1222208 2.444418e+00 528.8305 343.099673 343.449672 symmetric 68681 3.699997
58 343.649672 1321.8338 1321.34 1320.87 1321.105 0.73 0.7287 0.0006 0.150 0.0390887 7.817742e-01 264.2606 343.599672 343.749672 symmetric 68731 0.249999
59 348.649668 1365.6121 1322.90 1365.23 1344.065 21.55 21.5471 0.0160 2.300 -3.4619671 -6.923941e+01 3155.0961 346.399670 348.699667 symmetric 69731 4.999996
60 352.449664 1464.7743 1456.16 1463.97 1460.065 4.71 4.7093 0.0032 0.800 1.3365654 2.673133e+01 1242.3926 351.749665 352.549664 symmetric 70491 3.799996
61 362.449654 1666.3065 1662.80 1665.35 1664.075 2.23 2.2315 0.0013 0.200 0.2021988 4.043980e+00 416.2206 362.349654 362.549654 symmetric 72491 9.999990
62 362.699654 1666.0626 1664.85 1662.53 1663.690 2.37 2.3726 0.0014 0.150 0.1549450 3.098903e+00 332.8934 362.649654 362.799654 symmetric 72541 0.250000
63 362.849654 1663.0930 1662.53 1661.54 1662.035 1.06 1.0580 0.0006 0.100 0.0527770 1.055542e+00 249.3582 362.799654 362.899654 symmetric 72571 0.150000
64 362.949654 1662.0226 1661.54 1658.53 1660.035 1.99 1.9876 0.0012 0.100 0.0994222 1.988445e+00 249.1043 362.899654 362.999654 symmetric 72591 0.100000
65 363.049654 1659.4623 1658.53 1656.09 1657.310 2.15 2.1523 0.0013 0.100 0.1076664 2.153330e+00 248.7038 362.999654 363.099654 symmetric 72611 0.100000
66 363.149654 1656.7284 1656.09 1237.13 1446.610 210.12 210.1184 0.1452 4.750 -77.1999421 -1.544000e+03 6866.5305 363.099654 367.849649 symmetric 72631 0.100000
67 422.749597 1273.9955 1235.73 1265.20 1250.465 23.53 23.5305 0.0188 4.850 32.0999111 6.419988e+02 6159.3782 418.849601 423.699596 symmetric 84551 59.599943
68 428.899591 1243.6479 1227.01 1234.28 1230.645 13.00 13.0029 0.0106 1.950 13.1284861 2.625700e+02 2474.4129 427.949592 429.899590 symmetric 85781 6.149994
69 432.349588 1360.7392 1235.05 1359.63 1297.340 63.40 63.3992 0.0489 2.350 9.0325586 1.806513e+02 3122.6545 430.049590 432.399588 symmetric 86471 3.449997
70 432.449588 1360.6157 1359.63 1330.93 1345.280 15.34 15.3357 0.0114 1.605 1.3418763 2.683755e+01 2221.1305 432.399588 434.004586 symmetric 86491 0.100000
71 435.999584 1351.9581 1350.70 1351.61 1351.155 0.80 0.8031 0.0006 0.150 0.0396131 7.922619e-01 270.2697 435.899584 436.049584 symmetric 87201 3.549996
72 436.599584 1353.1217 1352.11 1351.38 1351.745 1.38 1.3767 0.0010 0.250 0.1445404 2.890812e+00 405.6680 436.499584 436.749583 symmetric 87321 0.600000
73 441.249579 1292.6840 1292.12 1292.33 1292.225 0.46 0.4590 0.0004 0.200 0.0459497 9.189951e-01 323.1021 441.099579 441.299579 symmetric 88251 4.649995
74 441.849579 1296.9991 1293.39 1296.64 1295.015 1.98 1.9841 0.0015 0.400 0.0161521 3.230424e-01 582.7727 441.499579 441.899579 symmetric 88371 0.600000
75 448.299572 1380.7390 1380.20 1380.17 1380.185 0.55 0.5540 0.0004 0.250 0.0579765 1.159531e+00 414.1123 448.199573 448.449572 symmetric 89661 6.449993
76 456.149565 1872.1765 1871.33 1678.02 1774.675 97.50 97.5015 0.0549 2.100 44.6004750 8.920104e+02 3860.1508 456.099565 458.199563 symmetric 91231 7.849993
77 462.949558 1607.9312 1607.18 1065.00 1336.090 271.84 271.8412 0.2035 13.150 -692.0889626 -1.384179e+04 16944.2664 462.799559 475.949546 symmetric 92591 6.799993
78 499.149524 991.3258 990.89 991.05 990.970 0.36 0.3558 0.0004 0.150 0.0151307 3.026139e-01 198.2092 499.049524 499.199524 symmetric 99831 36.199966
79 516.199508 1020.7089 1020.40 1020.45 1020.425 0.28 0.2839 0.0003 0.150 0.0228912 4.578250e-01 204.1081 516.149508 516.299508 symmetric 103241 17.049984
80 529.699495 989.1727 988.87 988.78 988.825 0.35 0.3477 0.0004 0.105 0.0170750 3.414998e-01 148.3410 529.649495 529.754495 symmetric 105941 13.499987
81 557.699468 1615.9978 1615.42 1615.57 1615.495 0.50 0.5028 0.0003 0.150 0.0387117 7.742341e-01 323.1377 557.649468 557.799468 symmetric 111541 27.999973
82 558.149468 1616.9164 1615.68 1615.96 1615.820 1.10 1.0964 0.0007 0.150 0.0739744 1.479489e+00 323.2376 558.049468 558.199468 symmetric 111631 0.450000
83 561.249465 1714.0422 1712.40 1712.69 1712.545 1.50 1.4972 0.0009 0.200 0.2036269 4.072541e+00 428.3402 561.149465 561.349465 symmetric 112251 3.099997
84 580.449446 1266.0193 1265.69 1265.47 1265.580 0.44 0.4393 0.0003 0.100 0.0220445 4.408896e-01 189.8586 580.399446 580.499446 symmetric 116091 19.199981
85 580.549446 1266.4380 1265.47 1265.94 1265.705 0.73 0.7331 0.0006 0.150 0.0649588 1.299177e+00 253.2056 580.499446 580.649446 symmetric 116111 0.100000
86 580.699446 1266.3593 1265.94 1265.90 1265.920 0.44 0.4393 0.0003 0.100 0.0218442 4.368837e-01 189.9100 580.649446 580.749446 symmetric 116141 0.150000
87 594.699433 1199.7690 1199.43 1198.91 1199.170 0.60 0.5990 0.0005 0.205 0.0383461 7.669222e-01 299.8319 594.649433 594.854433 symmetric 118941 13.999987
88 605.749422 1360.7790 1360.02 1360.04 1360.030 0.75 0.7490 0.0006 0.150 0.0477232 9.544657e-01 272.0536 605.699422 605.849422 symmetric 121151 11.049989
89 638.499391 1080.4379 1076.33 1038.52 1057.425 23.01 23.0129 0.0218 3.405 20.1367512 4.027354e+02 3668.3535 637.699392 641.104389 symmetric 127701 32.749969
90 722.104311 1079.8940 1079.58 1079.50 1079.540 0.35 0.3540 0.0003 0.200 0.0363747 7.274950e-01 269.9212 722.049311 722.249311 symmetric 144422 83.604920
91 746.249288 1426.1715 1421.32 1425.69 1423.505 2.67 2.6665 0.0019 0.600 0.4471515 8.943039e+00 925.7250 745.749289 746.349288 symmetric 149251 24.144977
92 749.649285 1536.5539 1534.24 1536.11 1535.175 1.38 1.3790 0.0009 0.200 0.1244955 2.489911e+00 383.9183 749.549285 749.749285 symmetric 149931 3.399997
93 749.799285 1537.4210 1536.11 1536.70 1536.405 1.02 1.0160 0.0007 0.150 0.0683724 1.367449e+00 307.3490 749.749285 749.899285 symmetric 149961 0.150000
94 750.749284 1544.6969 1537.42 1544.27 1540.845 3.85 3.8519 0.0025 0.750 -0.0457384 -9.147684e-01 1232.6275 750.049285 750.799284 symmetric 150151 0.949999
95 750.999284 1546.2080 1544.27 1545.66 1544.965 1.24 1.2430 0.0008 0.300 0.1569068 3.138140e+00 540.8926 750.799284 751.099284 symmetric 150201 0.250000
96 751.304283 1548.6216 1545.66 1546.55 1546.105 2.52 2.5166 0.0016 0.350 0.5026857 1.005372e+01 618.9436 751.099284 751.449283 symmetric 150262 0.304999
97 751.549283 1548.4352 1546.55 1547.28 1546.915 1.52 1.5202 0.0010 0.150 0.1176765 2.353532e+00 309.5001 751.449283 751.599283 symmetric 150311 0.245000
98 751.649283 1548.0860 1547.28 1545.79 1546.535 1.55 1.5510 0.0010 0.200 0.1207685 2.415373e+00 386.7534 751.599283 751.799283 symmetric 150331 0.100000
99 752.649282 1533.6248 1532.87 1533.04 1532.955 0.67 0.6698 0.0004 0.100 0.0334120 6.682408e-01 229.9767 752.599282 752.699282 symmetric 150531 0.999999
100 757.799277 2215.3662 1644.82 2067.16 1855.990 359.38 359.3762 0.1936 1.800 -8.9127498 -1.782552e+02 3424.6704 754.199281 755.999279 symmetric 151561 5.149995
101 758.099277 2218.1729 2212.95 2214.73 2213.840 4.33 4.3329 0.0020 0.350 0.8213744 1.642750e+01 886.3579 757.899277 758.249277 symmetric 151621 0.300000
102 758.299277 2217.4570 2214.73 2214.83 2214.780 2.68 2.6770 0.0012 0.100 0.1338983 2.677968e+00 332.3504 758.249277 758.349277 symmetric 151661 0.200000
103 758.399277 2217.3937 2214.83 2215.35 2215.090 2.30 2.3037 0.0010 0.150 0.1963625 3.927255e+00 443.2139 758.349277 758.499277 symmetric 151681 0.100000
104 758.649276 2214.5860 2212.92 2213.55 2213.235 1.35 1.3510 0.0006 0.100 0.0673835 1.347672e+00 332.0528 758.599277 758.699276 symmetric 151731 0.249999
105 758.799276 2214.7150 2213.55 2212.09 2212.820 1.89 1.8950 0.0009 0.150 0.1715855 3.431713e+00 442.7355 758.699276 758.849276 symmetric 151761 0.150000
106 758.899276 2213.4347 2212.09 2212.14 2212.115 1.32 1.3197 0.0006 0.100 0.0659364 1.318728e+00 331.8830 758.849276 758.949276 symmetric 151781 0.100000
107 759.149276 2209.1088 2208.10 2204.05 2206.075 3.03 3.0338 0.0014 0.205 0.2328372 4.656748e+00 551.7632 759.099276 759.304276 symmetric 151831 0.250000
108 760.299275 2174.3042 2173.46 1849.76 2011.610 162.69 162.6942 0.0809 3.350 78.7267086 1.574536e+03 6918.2043 760.249275 763.599272 symmetric 152061 1.149999
109 765.299270 1726.8773 1726.26 1716.17 1721.215 5.66 5.6622 0.0033 0.350 0.2605282 5.210569e+00 688.7450 765.249270 765.599270 symmetric 153061 4.999995
110 765.649270 1717.9451 1716.17 1715.48 1715.825 2.12 2.1201 0.0012 0.100 0.1060725 2.121452e+00 257.4794 765.599270 765.699270 symmetric 153131 0.350000
111 765.749270 1715.9180 1715.48 1714.89 1715.185 0.73 0.7330 0.0004 0.100 0.0365187 7.303746e-01 257.3144 765.699270 765.799270 symmetric 153151 0.100000
112 774.299262 2071.7456 1839.21 1931.09 1885.150 186.60 186.5956 0.0990 2.650 2.3968735 4.793752e+01 5092.2977 768.949267 771.599264 symmetric 154861 8.549992
113 774.499261 2076.7349 2071.16 2075.69 2073.425 3.31 3.3099 0.0016 0.200 0.2307582 4.615169e+00 518.5862 774.349262 774.549261 symmetric 154901 0.199999
114 775.249261 2088.4192 2085.47 2087.59 2086.530 1.89 1.8892 0.0009 0.250 0.1613125 3.226253e+00 626.1189 775.049261 775.299261 symmetric 155051 0.750000
115 775.849260 2095.2645 2092.87 2094.10 2093.485 1.78 1.7795 0.0009 0.150 0.1087675 2.175353e+00 418.8056 775.749260 775.899260 symmetric 155171 0.599999
116 775.999260 2095.5150 2094.10 2094.82 2094.460 1.05 1.0550 0.0005 0.200 0.0917987 1.835975e+00 523.7068 775.899260 776.099260 symmetric 155201 0.150000
117 776.149260 2095.8220 2094.82 2095.22 2095.020 0.80 0.8020 0.0004 0.150 0.0541121 1.082243e+00 419.0585 776.099260 776.249260 symmetric 155231 0.150000
118 776.399260 2095.5489 2094.63 2094.28 2094.455 1.09 1.0939 0.0005 0.100 0.0547550 1.095101e+00 314.2225 776.349260 776.449260 symmetric 155281 0.250000
119 776.499259 2095.0915 2094.28 2092.05 2093.165 1.93 1.9265 0.0009 0.200 0.1707118 3.414239e+00 523.4612 776.449260 776.649259 symmetric 155301 0.099999
120 776.699259 2092.8210 2092.05 2091.40 2091.725 1.10 1.0960 0.0005 0.155 0.0758240 1.516482e+00 418.4224 776.649259 776.804259 symmetric 155341 0.200000
121 778.699257 2054.2570 2052.73 2050.95 2051.840 2.42 2.4170 0.0012 0.105 0.1161136 2.322274e+00 307.8982 778.649257 778.754257 symmetric 155741 1.999998
122 778.899257 2052.7490 2050.95 2051.31 2051.130 1.62 1.6190 0.0008 0.245 0.2036274 4.072551e+00 512.9772 778.754257 778.999257 symmetric 155781 0.200000
123 779.149257 2051.6971 2050.13 2049.25 2049.690 2.01 2.0071 0.0010 0.155 0.1632415 3.264833e+00 410.1042 779.099257 779.254257 symmetric 155831 0.250000
124 779.749256 2044.2882 2043.39 2042.36 2042.875 1.41 1.4132 0.0007 0.105 0.0684172 1.368346e+00 306.5029 779.699256 779.804256 symmetric 155951 0.599999
125 780.299256 2037.1223 2036.49 2034.47 2035.480 1.64 1.6423 0.0008 0.105 0.0759151 1.518304e+00 305.4047 780.249256 780.354256 symmetric 156061 0.550000
126 781.449255 2022.9592 2021.98 2021.22 2021.600 1.36 1.3592 0.0007 0.250 0.1589145 3.178294e+00 606.6374 781.299255 781.549255 symmetric 156291 1.149999
127 781.649255 2022.2341 2021.22 2019.58 2020.400 1.83 1.8341 0.0009 0.150 0.1801260 3.602523e+00 404.2591 781.549255 781.699255 symmetric 156331 0.200000
128 781.749254 2021.7850 2019.58 2020.63 2020.105 1.68 1.6800 0.0008 0.100 0.0841657 1.683316e+00 303.0991 781.699255 781.799254 symmetric 156351 0.099999
129 781.949254 2024.0695 2020.63 2022.56 2021.595 2.47 2.4745 0.0012 0.250 0.3317118 6.634242e+00 606.8090 781.799254 782.049254 symmetric 156391 0.200000
130 782.099254 2025.9306 2022.56 2024.82 2023.690 2.24 2.2406 0.0011 0.150 0.1889535 3.779073e+00 404.9268 782.049254 782.199254 symmetric 156421 0.150000
131 782.349254 2026.2970 2025.17 2025.66 2025.415 0.88 0.8820 0.0004 0.100 0.0439985 8.799708e-01 303.8563 782.299254 782.399254 symmetric 156471 0.250000
132 782.454254 2027.7803 2025.66 2026.82 2026.240 1.54 1.5403 0.0008 0.150 0.1448932 2.897867e+00 405.3925 782.399254 782.549254 symmetric 156492 0.105000
133 783.199253 2037.0732 2033.02 2036.24 2034.630 2.44 2.4432 0.0012 0.150 0.1151845 2.303692e+00 407.0412 783.099253 783.249253 symmetric 156641 0.744999
134 783.299253 2038.8276 2036.24 2037.14 2036.690 2.14 2.1376 0.0011 0.100 0.1067958 2.135918e+00 305.6103 783.249253 783.349253 symmetric 156661 0.100000
135 783.499253 2038.1051 2037.52 2035.83 2036.675 1.43 1.4301 0.0007 0.100 0.0713821 1.427644e+00 305.5727 783.449253 783.549253 symmetric 156701 0.200000
136 783.599253 2037.5524 2035.83 2036.50 2036.165 1.39 1.3874 0.0007 0.150 0.0899028 1.798057e+00 407.3225 783.549253 783.699253 symmetric 156721 0.100000
137 783.749253 2038.8995 2036.50 2031.66 2034.080 4.82 4.8195 0.0024 0.300 0.3300309 6.600624e+00 712.2563 783.699253 783.999252 symmetric 156751 0.150000
138 788.499248 1830.8964 1829.97 1825.31 1827.640 3.26 3.2564 0.0018 0.200 0.1164886 2.329774e+00 457.0260 788.449248 788.649248 symmetric 157701 4.749995
139 788.699248 1825.8623 1825.31 1824.09 1824.700 1.16 1.1623 0.0006 0.100 0.0582656 1.165313e+00 273.7625 788.649248 788.749248 symmetric 157741 0.200000
140 788.799248 1825.2882 1824.09 1822.57 1823.330 1.96 1.9582 0.0011 0.300 0.1771506 3.543016e+00 638.3410 788.749248 789.049248 symmetric 157761 0.100000
141 789.149247 1824.4200 1822.57 1822.73 1822.650 1.77 1.7700 0.0010 0.200 0.1323992 2.647987e+00 455.7943 789.049248 789.249247 symmetric 157831 0.349999
142 789.549247 1826.0424 1823.44 1824.36 1823.900 2.14 2.1424 0.0012 0.200 0.1793708 3.587419e+00 456.1535 789.399247 789.599247 symmetric 157911 0.400000
143 790.099247 1832.9691 1826.82 1832.14 1829.480 3.49 3.4891 0.0019 0.400 0.1866627 3.733258e+00 823.4519 789.749247 790.149246 symmetric 158021 0.550000
144 790.399246 1839.6320 1834.37 1837.84 1836.105 3.53 3.5270 0.0019 0.200 0.2672997 5.345999e+00 459.2934 790.249246 790.449246 symmetric 158081 0.299999
145 790.499246 1839.7987 1837.84 1839.07 1838.455 1.34 1.3437 0.0007 0.100 0.0669948 1.339897e+00 275.8355 790.449246 790.549246 symmetric 158101 0.100000
146 790.749246 1843.1069 1839.07 1841.86 1840.465 2.64 2.6419 0.0014 0.250 0.3866713 7.733433e+00 552.5258 790.549246 790.799246 symmetric 158151 0.250000
147 790.949246 1845.0887 1843.05 1842.63 1842.840 2.25 2.2487 0.0012 0.105 0.1112594 2.225190e+00 276.5385 790.899246 791.004246 symmetric 158191 0.200000
148 792.849244 1807.8164 1806.50 1803.80 1805.150 2.67 2.6664 0.0015 0.150 0.1886166 3.772336e+00 361.2174 792.799244 792.949244 symmetric 158571 1.899998
149 793.049244 1804.9125 1803.80 1803.11 1803.455 1.46 1.4575 0.0008 0.150 0.1181634 2.363271e+00 360.8080 792.949244 793.099244 symmetric 158611 0.200000
150 793.149244 1805.0261 1803.11 1801.77 1802.440 2.59 2.5861 0.0014 0.200 0.2520743 5.041491e+00 450.8615 793.099244 793.299243 symmetric 158631 0.100000
151 793.399243 1805.1576 1801.77 1804.22 1802.995 2.16 2.1626 0.0012 0.200 0.2357888 4.715781e+00 450.9839 793.299243 793.499243 symmetric 158681 0.249999
152 793.599243 1805.1888 1804.22 1804.00 1804.110 1.08 1.0788 0.0006 0.200 0.1065607 2.131215e+00 451.1334 793.499243 793.699243 symmetric 158721 0.200000
153 793.749243 1804.4885 1804.00 1803.64 1803.820 0.67 0.6684 0.0004 0.100 0.0334563 6.691269e-01 270.6061 793.699243 793.799243 symmetric 158751 0.150000
154 793.849243 1804.5399 1803.64 1803.19 1803.415 1.12 1.1249 0.0006 0.200 0.1098033 2.196068e+00 450.9621 793.799243 793.999243 symmetric 158771 0.100000
155 794.199243 1803.2099 1802.03 1800.11 1801.070 2.14 2.1399 0.0012 0.150 0.1708009 3.416022e+00 360.3851 794.149243 794.299242 symmetric 158841 0.350000
156 797.199240 1705.8269 1704.80 1699.97 1702.385 3.44 3.4419 0.0020 0.500 0.0832136 1.664273e+00 936.3930 797.149240 797.649239 symmetric 159441 2.999997
157 797.699239 1700.5361 1699.97 1699.04 1699.505 1.03 1.0311 0.0006 0.150 0.0915185 1.830372e+00 339.9916 797.649239 797.799239 symmetric 159541 0.499999
158 797.949239 1700.7155 1699.04 1698.28 1698.660 2.06 2.0555 0.0012 0.250 0.2868982 5.737970e+00 509.8843 797.799239 798.049239 symmetric 159591 0.250000
159 798.099239 1701.3344 1698.28 1698.91 1698.595 2.74 2.7394 0.0016 0.100 0.1369675 2.739352e+00 254.9260 798.049239 798.149239 symmetric 159621 0.150000
160 798.249239 1700.5516 1698.91 1698.86 1698.885 1.67 1.6666 0.0010 0.150 0.1447947 2.895896e+00 339.9218 798.149239 798.299239 symmetric 159651 0.150000
161 798.399239 1699.8544 1698.86 1698.37 1698.615 1.24 1.2394 0.0007 0.155 0.1034790 2.069581e+00 339.8286 798.299239 798.454239 symmetric 159681 0.150000
162 798.749238 1699.1622 1697.64 1697.99 1697.815 1.35 1.3472 0.0008 0.205 0.1079934 2.159869e+00 424.5602 798.599238 798.804238 symmetric 159751 0.349999
163 799.049238 1698.3680 1696.99 1697.10 1697.045 1.32 1.3230 0.0008 0.150 0.0897558 1.795118e+00 339.4985 798.949238 799.099238 symmetric 159811 0.300000
164 799.299238 1699.3802 1698.82 1698.15 1698.485 0.90 0.8952 0.0005 0.100 0.0448018 8.960359e-01 254.8172 799.249238 799.349238 symmetric 159861 0.250000
165 799.399238 1698.7997 1698.15 1698.25 1698.200 0.60 0.5997 0.0004 0.100 0.0298880 5.977613e-01 254.7599 799.349238 799.449238 symmetric 159881 0.100000
166 799.599237 1699.9654 1698.25 1699.16 1698.705 1.26 1.2604 0.0007 0.200 0.1519753 3.039509e+00 424.8278 799.449238 799.649237 symmetric 159921 0.199999
167 799.699237 1699.8954 1699.16 1699.17 1699.165 0.73 0.7304 0.0004 0.150 0.0537914 1.075829e+00 339.8866 799.649237 799.799237 symmetric 159941 0.100000
168 799.849237 1699.9499 1699.17 1698.20 1698.685 1.26 1.2649 0.0007 0.100 0.0632324 1.264650e+00 254.8658 799.799237 799.899237 symmetric 159971 0.150000
169 799.999237 1699.4522 1698.20 1697.85 1698.025 1.43 1.4272 0.0008 0.200 0.1373388 2.746779e+00 424.6423 799.899237 800.099237 symmetric 160001 0.150000
170 800.154237 1698.5709 1697.85 1697.22 1697.535 1.04 1.0359 0.0006 0.155 0.0997518 1.995038e+00 339.6080 800.099237 800.254237 symmetric 160032 0.155000
171 809.849228 1726.3067 1516.61 1724.84 1620.725 105.58 105.5818 0.0651 4.700 -11.1218690 -2.224376e+02 7687.3109 805.199232 809.899228 symmetric 161971 9.694991
172 809.954228 1725.9102 1724.84 1725.05 1724.945 0.97 0.9652 0.0006 0.205 0.1305126 2.610255e+00 431.3660 809.899228 810.104227 symmetric 161992 0.105000
173 810.249227 1725.3630 1724.80 1724.82 1724.810 0.55 0.5530 0.0003 0.150 0.0415444 8.308890e-01 345.0030 810.199227 810.349227 symmetric 162051 0.294999
174 810.499227 1726.1862 1724.82 1725.03 1724.925 1.26 1.2612 0.0007 0.200 0.1374278 2.748559e+00 431.3679 810.349227 810.549227 symmetric 162101 0.250000
175 810.599227 1727.7352 1725.03 1726.11 1725.570 2.17 2.1652 0.0013 0.100 0.1082480 2.164962e+00 258.9435 810.549227 810.649227 symmetric 162121 0.100000
176 810.949227 1729.5854 1726.11 1715.86 1720.985 8.60 8.6004 0.0050 0.755 3.9186013 7.837210e+01 1380.7347 810.649227 811.404226 symmetric 162191 0.350000
177 818.149220 1678.5531 1597.79 1677.77 1637.780 40.77 40.7731 0.0249 1.900 3.2912103 6.582427e+01 3196.9541 816.299222 818.199220 symmetric 163631 7.199993
178 818.399220 1678.5245 1678.05 1678.06 1678.055 0.47 0.4695 0.0003 0.100 0.0234860 4.697214e-01 251.7315 818.349220 818.449219 symmetric 163681 0.250000
179 818.499219 1678.5432 1678.06 1675.80 1676.930 1.61 1.6132 0.0010 0.205 0.1192451 2.384904e+00 419.3574 818.449219 818.654219 symmetric 163701 0.099999
180 824.349214 1925.3203 1706.06 1740.84 1723.450 201.87 201.8703 0.1171 0.345 -0.3268253 -6.536512e+00 602.0855 822.504216 822.849215 symmetric 164871 5.849995
181 826.549212 1886.1362 1885.46 1764.70 1825.080 61.06 61.0562 0.0335 1.500 6.0620012 1.212401e+02 2834.9329 826.499212 827.999210 symmetric 165311 2.199998
182 828.049210 1765.6798 1764.70 1756.94 1760.820 4.86 4.8598 0.0028 0.600 0.1075995 2.151992e+00 1144.6383 827.999210 828.599210 symmetric 165611 1.499998
183 828.649210 1757.4821 1756.94 1756.00 1756.470 1.01 1.0121 0.0006 0.150 0.0508252 1.016504e+00 351.3445 828.599210 828.749210 symmetric 165731 0.600000
184 830.499208 1935.5286 1756.00 1873.25 1814.625 120.90 120.9036 0.0666 2.750 134.0099077 2.680201e+03 5214.9560 828.749210 831.499207 symmetric 166101 1.849998
185 831.799207 1872.2450 1871.63 1867.81 1869.720 2.52 2.5250 0.0013 0.350 0.2990028 5.980062e+00 748.1860 831.749207 832.099206 symmetric 166361 1.299999
186 832.149206 1868.5270 1867.81 1868.06 1867.935 0.59 0.5920 0.0003 0.100 0.0296209 5.924191e-01 280.2195 832.099206 832.199206 symmetric 166431 0.349999
187 833.649205 1922.0046 1920.95 1921.34 1921.145 0.86 0.8596 0.0004 0.150 0.0503688 1.007378e+00 384.2791 833.549205 833.699205 symmetric 166731 1.499999
188 833.849205 1924.0857 1921.34 1920.69 1921.015 3.07 3.0707 0.0016 0.350 0.4127980 8.255968e+00 768.8180 833.699205 834.049205 symmetric 166771 0.200000
189 834.099205 1921.6218 1920.69 1920.56 1920.625 1.00 0.9968 0.0005 0.100 0.0498064 9.961280e-01 288.1434 834.049205 834.149204 symmetric 166821 0.250000
190 834.399204 1921.3433 1920.62 1920.01 1920.315 1.03 1.0283 0.0005 0.105 0.0496044 9.920893e-01 288.0989 834.349204 834.454204 symmetric 166881 0.299999
191 834.699204 1921.9203 1919.44 1909.46 1914.450 7.47 7.4703 0.0039 0.800 1.5392350 3.078473e+01 1628.8193 834.549204 835.349203 symmetric 166941 0.300000
192 836.499202 2004.8878 1909.46 1929.54 1919.500 85.39 85.3878 0.0445 2.005 87.0589020 1.741180e+03 4021.9812 835.349203 837.354201 symmetric 167301 1.799998
193 838.249201 1991.7666 1928.58 1887.83 1908.205 83.56 83.5616 0.0438 1.905 76.8787759 1.537577e+03 3797.9738 837.449201 839.354200 symmetric 167651 1.749999
194 841.349198 2027.2022 2023.86 1975.32 1999.590 27.61 27.6122 0.0138 1.250 0.0350216 7.004322e-01 2599.4954 841.299198 842.549196 symmetric 168271 3.099997
195 842.799196 1974.2240 1973.26 1971.26 1972.260 1.96 1.9640 0.0010 0.150 0.1344253 2.688508e+00 394.5861 842.749196 842.899196 symmetric 168561 1.449998
196 850.849189 1524.3273 1476.36 1523.73 1500.045 24.28 24.2823 0.0162 1.750 4.5737830 9.147575e+01 2704.6546 849.149190 850.899189 symmetric 170171 8.049993
197 850.949188 1524.1302 1523.73 1523.06 1523.395 0.74 0.7352 0.0005 0.200 0.0579707 1.159416e+00 380.9064 850.899189 851.099188 symmetric 170191 0.099999
198 851.249188 1523.9193 1523.06 1523.05 1523.055 0.86 0.8643 0.0006 0.200 0.0861730 1.723461e+00 380.8494 851.099188 851.299188 symmetric 170251 0.300000
199 851.349188 1523.4659 1523.05 1522.67 1522.860 0.61 0.6058 0.0004 0.150 0.0422844 8.456890e-01 304.6136 851.299188 851.449188 symmetric 170271 0.100000
200 852.399187 1542.6868 1540.36 1541.25 1540.805 1.88 1.8818 0.0012 0.250 0.2804676 5.609357e+00 462.5223 852.249187 852.499187 symmetric 170481 1.049999
201 852.549187 1541.7200 1541.25 1523.63 1532.440 9.28 9.2800 0.0061 0.605 0.9809838 1.961969e+01 997.1136 852.499187 853.104186 symmetric 170511 0.150000
202 854.799185 1569.1283 1521.55 1566.84 1544.195 24.93 24.9333 0.0161 1.500 3.6706230 7.341253e+01 2397.1711 853.349186 854.849185 symmetric 170961 2.249998
203 854.899185 1567.9755 1566.84 1566.09 1566.465 1.51 1.5105 0.0010 0.150 0.1373691 2.747385e+00 313.4297 854.849185 854.999185 symmetric 170981 0.100000
204 855.049185 1566.7673 1566.09 1566.17 1566.130 0.64 0.6373 0.0004 0.100 0.0318709 6.374187e-01 234.9511 854.999185 855.099185 symmetric 171011 0.150000
205 855.149184 1566.6824 1566.17 1558.36 1562.265 4.42 4.4174 0.0028 0.550 0.0693417 1.386835e+00 937.4289 855.099185 855.649184 symmetric 171031 0.099999
206 855.749184 1559.7816 1558.36 1559.31 1558.835 0.95 0.9466 0.0006 0.150 0.0564843 1.129687e+00 311.8235 855.649184 855.799184 symmetric 171151 0.600000
207 855.899184 1560.5275 1559.31 1560.13 1559.720 0.81 0.8075 0.0005 0.200 0.0990287 1.980576e+00 390.0287 855.799184 855.999184 symmetric 171181 0.150000
208 857.649182 1580.2659 1549.95 1579.40 1564.675 15.59 15.5909 0.0100 1.050 2.4452695 4.890544e+01 1723.5871 856.699183 857.749182 symmetric 171531 1.749998
209 857.799182 1579.9025 1579.40 1556.62 1568.010 11.89 11.8925 0.0076 0.755 0.1180741 2.361483e+00 1254.5849 857.749182 858.504181 symmetric 171561 0.150000
210 858.849181 1555.2396 1554.69 1553.84 1554.265 0.97 0.9746 0.0006 0.150 0.0813289 1.626580e+00 310.9340 858.799181 858.949181 symmetric 171771 1.049999
211 859.049181 1555.1323 1553.84 1554.44 1554.140 0.99 0.9923 0.0006 0.200 0.1262033 2.524069e+00 388.6607 858.949181 859.149181 symmetric 171811 0.200000
212 859.399180 1557.3869 1554.87 1555.82 1555.345 2.04 2.0419 0.0013 0.250 0.1796073 3.592148e+00 466.7835 859.249181 859.499180 symmetric 171881 0.349999
213 859.699180 1556.1443 1555.45 1553.81 1554.630 1.51 1.5143 0.0010 0.200 0.1226422 2.452847e+00 388.7792 859.649180 859.849180 symmetric 171941 0.300000
214 859.899180 1554.2090 1553.81 1550.75 1552.280 1.93 1.9290 0.0012 0.255 0.1598060 3.196123e+00 465.8512 859.849180 860.104180 symmetric 171981 0.200000
215 860.899179 1538.7753 1537.43 1538.11 1537.770 1.01 1.0053 0.0007 0.150 0.0780409 1.560819e+00 307.6313 860.849179 860.999179 symmetric 172181 0.999999
216 861.049179 1538.7699 1538.11 1537.96 1538.035 0.73 0.7349 0.0005 0.150 0.0487457 9.749159e-01 307.6550 860.999179 861.149179 symmetric 172211 0.150000
217 861.249179 1540.4204 1537.96 1540.02 1538.990 1.43 1.4304 0.0009 0.150 0.0774986 1.549973e+00 307.8750 861.149179 861.299179 symmetric 172251 0.200000
218 861.949178 1551.2470 1548.30 1550.76 1549.530 1.72 1.7170 0.0011 0.200 0.0950737 1.901476e+00 387.4772 861.799178 861.999178 symmetric 172391 0.699999
219 864.599175 1692.7021 1689.61 1691.71 1690.660 2.04 2.0421 0.0012 0.100 0.1020976 2.041954e+00 253.7009 864.549176 864.649175 symmetric 172921 2.649997
220 864.899175 1695.4936 1692.04 1693.91 1692.975 2.52 2.5186 0.0015 0.200 0.1624615 3.249233e+00 423.4062 864.799175 864.999175 symmetric 172981 0.300000
221 865.649174 1714.9563 1699.17 1713.30 1706.235 8.72 8.7213 0.0051 0.500 0.6568742 1.313750e+01 939.0835 865.199175 865.699174 symmetric 173131 0.749999
222 865.899174 1719.9895 1717.94 1719.37 1718.655 1.33 1.3345 0.0008 0.100 0.0667170 1.334342e+00 257.8647 865.849174 865.949174 symmetric 173181 0.250000
223 867.749172 1806.1346 1724.28 1804.95 1764.615 41.52 41.5196 0.0235 1.600 2.3904141 4.780833e+01 2914.0098 866.199174 867.799172 symmetric 173551 1.849998
224 868.099172 1799.8994 1799.13 1749.68 1774.405 25.49 25.4944 0.0144 1.000 1.0687871 2.137576e+01 1864.1930 868.049172 869.049171 symmetric 173621 0.350000
225 869.349171 1752.1420 1748.64 1751.56 1750.100 2.04 2.0420 0.0012 0.250 0.1865274 3.730552e+00 525.2171 869.199171 869.449171 symmetric 173871 1.249999
226 869.699171 1757.5659 1751.56 1756.86 1754.210 3.36 3.3559 0.0019 0.300 0.1133943 2.267889e+00 614.0875 869.449171 869.749171 symmetric 173941 0.350000
227 869.849170 1759.8504 1756.86 1758.56 1757.710 2.14 2.1404 0.0012 0.200 0.2511098 5.022201e+00 439.6783 869.749171 869.949170 symmetric 173971 0.149999
228 869.999170 1759.2825 1758.56 1757.55 1758.055 1.23 1.2275 0.0007 0.150 0.0602933 1.205868e+00 351.6703 869.949170 870.099170 symmetric 174001 0.150000
229 870.149170 1758.2067 1757.55 1746.90 1752.225 5.98 5.9817 0.0034 0.555 1.3281941 2.656391e+01 1052.6892 870.099170 870.654170 symmetric 174031 0.150000
230 881.399159 1624.9875 1388.65 1560.35 1474.500 150.49 150.4875 0.1021 6.900 347.4165312 6.948337e+03 10595.1790 876.149164 883.049158 symmetric 176281 11.249989
231 884.149157 1588.3262 1587.76 1567.26 1577.510 10.82 10.8162 0.0069 0.705 0.5197226 1.039446e+01 1183.7043 884.099157 884.804156 symmetric 176831 2.749998
232 886.049155 1539.0577 1537.72 1538.57 1538.145 0.91 0.9127 0.0006 0.200 0.0713116 1.426233e+00 384.6077 885.899155 886.099155 symmetric 177211 1.899998
233 886.249155 1542.3848 1538.57 1540.94 1539.755 2.63 2.6298 0.0017 0.200 0.2617347 5.234699e+00 385.2006 886.099155 886.299155 symmetric 177251 0.200000
234 886.349155 1542.5578 1540.94 1506.23 1523.585 18.97 18.9728 0.0125 1.350 -0.6924976 -1.384997e+01 2132.3259 886.299155 887.649153 symmetric 177271 0.100000
235 888.949152 1535.0426 1506.50 1534.42 1520.460 14.58 14.5826 0.0096 1.250 -2.4317080 -4.863421e+01 1974.1635 887.749153 888.999152 symmetric 177791 2.599997
236 889.149152 1536.8450 1534.42 1535.31 1534.865 1.98 1.9800 0.0013 0.300 0.3723053 7.446114e+00 537.5735 888.999152 889.299152 symmetric 177831 0.200000
237 889.349152 1536.7371 1535.31 1533.73 1534.520 2.22 2.2171 0.0014 0.250 0.1654439 3.308882e+00 460.5204 889.299152 889.549152 symmetric 177871 0.200000
238 889.599152 1534.1780 1533.73 1528.83 1531.280 2.90 2.8980 0.0019 0.500 -0.0547390 -1.094781e+00 842.1495 889.549152 890.049151 symmetric 177921 0.250000
239 890.649151 1537.9486 1528.83 1536.05 1532.440 5.51 5.5086 0.0036 0.700 1.7651013 3.530206e+01 1151.0946 890.049151 890.749151 symmetric 178131 1.049999
240 897.449144 1358.8643 1358.24 1354.38 1356.310 2.55 2.5543 0.0019 0.700 0.5902119 1.180425e+01 1017.8233 897.299144 897.999144 symmetric 179491 6.799993
241 898.149143 1354.8233 1354.38 1354.31 1354.345 0.48 0.4783 0.0004 0.250 0.0627815 1.255631e+00 406.3657 897.999144 898.249143 symmetric 179631 0.699999
shade_flecks.df <- find_flecks(time = chunk$chunk.time,
                               var = chunk$Q_PAR,
                               minTime = 0.001, 
                               minAmp = 0, 
                               minPdiff = 0.00025,
                               asymmetry = 0,
                               shadeflecks = TRUE,
                               asmMethod = "max",
                               time.digits = 6,
                               var.digits = 5,
                               verbose = FALSE)
#> Found 231 shadeflecks

shade_flecks.fig <-
  ggplot(data = chunk) +
  geom_line(aes(chunk.time, Q_PAR), linewidth = 0.2, na.rm = TRUE) +
  geom_point(aes(peakTime, peak), data = shade_flecks.df,
             size = 0.5, colour = "blue", na.rm = TRUE) +
  geom_point(aes(baselineTime1, baseline1), data = shade_flecks.df,
             size = 0.5, colour = "chartreuse3", na.rm = TRUE) +
  geom_point(aes(baselineTime2, baseline2), data = shade_flecks.df, 
             size = 0.5, colour = "chartreuse3", na.rm = TRUE) +
  geom_text(aes(peakTime, peak, label = no), data = shade_flecks.df, 
            colour = "blue", fontface = "bold",
            position = position_nudge(y = 30), size = 2.7, na.rm = TRUE)

shade_flecks.fig


shade_flecks.fig + xlim(0, 100)


knitr::kable(shade_flecks.df)
no peakTime peak baseline1 baseline2 baseline amplitude relAmp percDiff duration irradGain relIrrad irradTot baselineTime1 baselineTime2 symmetry z timeInterval
1 0.599999 1393.1500 1395.11 1394.27 1394.690 1.54 1.5400 0.0011 0.150 -0.0895977 -1.791956e+00 278.8476 0.500000 0.649999 symmetric 121 NA
2 2.899997 1383.2990 1383.70 1383.72 1383.710 0.41 0.4110 0.0003 0.100 -0.0205982 -4.119647e-01 207.5359 2.849997 2.949997 symmetric 581 2.299998
3 4.899995 1371.4480 1374.14 1372.15 1373.145 1.70 1.6970 0.0012 0.500 -0.3584355 -7.168716e+00 754.8695 4.499996 4.999995 symmetric 981 1.999998
4 9.299991 1447.6011 1448.47 1448.04 1448.255 0.65 0.6539 0.0004 0.200 -0.0250330 -5.006595e-01 362.0394 9.149991 9.349991 symmetric 1861 4.399996
5 14.149987 1356.7862 1358.99 1357.39 1358.190 1.40 1.4038 0.0010 0.250 -0.1442085 -2.884173e+00 407.3124 13.999987 14.249986 symmetric 2831 4.849996
6 14.449986 1355.8561 1357.39 1357.59 1357.490 1.63 1.6339 0.0012 0.250 -0.2563694 -5.127393e+00 406.9903 14.249986 14.499986 symmetric 2891 0.299999
7 14.549986 1356.5943 1357.59 1358.03 1357.810 1.22 1.2157 0.0009 0.150 -0.1111978 -2.223959e+00 271.4508 14.499986 14.649986 symmetric 2911 0.100000
8 14.699986 1357.6921 1358.03 1408.23 1383.130 25.44 25.4379 0.0184 2.650 -4.1943599 -8.388728e+01 3730.2632 14.649986 17.299984 symmetric 2941 0.150000
9 20.149981 1376.2953 1405.79 1376.72 1391.255 14.96 14.9597 0.0108 1.950 -4.2060751 -8.412158e+01 2778.2990 18.249983 20.199981 symmetric 4031 5.449995
10 20.349981 1375.5371 1376.72 1376.21 1376.465 0.93 0.9279 0.0007 0.300 -0.1333051 -2.666104e+00 481.6282 20.199981 20.499980 symmetric 4071 0.200000
11 20.599980 1375.4712 1376.21 1377.14 1376.675 1.20 1.2039 0.0009 0.250 -0.1402319 -2.804641e+00 412.8620 20.499980 20.749980 symmetric 4121 0.249999
12 20.799980 1376.5602 1377.14 1379.85 1378.495 1.93 1.9348 0.0014 0.350 -0.2192556 -4.385116e+00 551.1787 20.749980 21.099980 symmetric 4161 0.200000
13 25.999975 1570.9433 1571.57 1572.68 1572.125 1.18 1.1817 0.0008 0.250 -0.0712595 -1.425192e+00 471.5668 25.949975 26.199975 symmetric 5201 5.199995
14 27.749974 1581.7312 1582.15 1582.63 1582.390 0.66 0.6588 0.0004 0.150 -0.0537763 -1.075527e+00 316.4243 27.699974 27.849973 symmetric 5551 1.749999
15 28.199973 1583.6093 1584.09 1584.66 1584.375 0.77 0.7657 0.0005 0.150 -0.0337522 -6.750442e-01 316.8404 28.149973 28.299973 symmetric 5641 0.449999
16 28.499973 1582.2823 1584.66 1582.78 1583.720 1.44 1.4377 0.0009 0.250 -0.1232978 -2.465957e+00 474.9914 28.299973 28.549973 symmetric 5701 0.300000
17 29.999971 1564.8992 1567.24 1565.88 1566.560 1.66 1.6608 0.0011 0.250 -0.1689728 -3.379460e+00 469.7976 29.849972 30.099971 symmetric 6001 1.499998
18 30.199971 1564.6883 1565.88 1565.61 1565.745 1.06 1.0567 0.0007 0.150 -0.0800803 -1.601607e+00 313.0682 30.099971 30.249971 symmetric 6041 0.200000
19 30.299971 1564.9374 1565.61 1566.91 1566.260 1.32 1.3226 0.0008 0.300 -0.1165303 -2.330608e+00 548.0744 30.249971 30.549971 symmetric 6061 0.100000
20 30.599971 1566.4153 1566.91 1567.84 1567.375 0.96 0.9597 0.0006 0.100 -0.0480621 -9.612435e-01 235.0582 30.549971 30.649971 symmetric 6121 0.300000
21 30.899971 1569.0083 1569.52 1573.51 1571.515 2.51 2.5067 0.0016 0.300 -0.2427921 -4.855846e+00 549.7875 30.849971 31.149970 symmetric 6181 0.300000
22 31.249970 1572.9691 1573.51 1575.54 1574.525 1.56 1.5559 0.0010 0.200 -0.1984645 -3.969293e+00 393.4322 31.149970 31.349970 symmetric 6251 0.349999
23 37.999964 1689.1230 1690.51 1690.33 1690.420 1.30 1.2970 0.0008 0.200 -0.1279900 -2.559803e+00 422.4761 37.949964 38.149964 symmetric 7601 6.749994
24 38.299963 1689.2474 1689.74 1690.13 1689.935 0.69 0.6876 0.0004 0.150 -0.0345454 -6.909089e-01 337.9513 38.249964 38.399963 symmetric 7661 0.299999
25 39.049963 1690.3685 1691.06 1691.57 1691.315 0.95 0.9465 0.0006 0.200 -0.0904526 -1.809055e+00 422.7377 38.999963 39.199963 symmetric 7811 0.750000
26 53.849949 1842.3325 1842.83 1848.00 1845.415 3.08 3.0825 0.0017 0.300 -0.0098611 -1.972220e-01 645.8848 53.799949 54.099948 symmetric 10771 14.799986
27 54.599948 1850.1045 1850.59 1852.63 1851.610 1.51 1.5055 0.0008 0.200 -0.1621190 -3.242382e+00 462.7407 54.549948 54.749948 symmetric 10921 0.749999
28 54.899948 1850.5658 1852.63 1851.69 1852.160 1.59 1.5942 0.0009 0.250 -0.1789668 -3.579340e+00 555.4694 54.749948 54.999948 symmetric 10981 0.300000
29 55.149947 1849.0492 1851.69 1849.60 1850.645 1.60 1.5958 0.0009 0.200 -0.0306214 -6.124279e-01 462.6308 54.999948 55.199947 symmetric 11031 0.249999
30 62.749940 1675.5352 1676.57 1676.47 1676.520 0.98 0.9848 0.0006 0.150 -0.0621929 -1.243859e+00 335.2411 62.699940 62.849940 symmetric 12551 7.599993
31 63.399940 1679.2968 1679.99 1680.74 1680.365 1.07 1.0682 0.0006 0.150 -0.0789360 -1.578721e+00 335.9945 63.349940 63.499939 symmetric 12681 0.650000
32 70.649933 1979.3063 1980.27 1983.86 1982.065 2.76 2.7587 0.0014 0.150 -0.1369955 -2.739913e+00 396.2761 70.599933 70.749933 symmetric 14131 7.249993
33 70.799932 1982.4419 1983.86 2008.21 1996.035 13.59 13.5931 0.0068 1.100 1.9745191 3.949042e+01 2297.4156 70.749933 71.849931 symmetric 14161 0.149999
34 74.649929 1935.5587 2007.52 1936.37 1971.945 36.39 36.3863 0.0185 2.700 -10.7952459 -2.159051e+02 5412.0503 72.049931 74.749929 symmetric 14931 3.849997
35 74.799929 1934.8372 1936.37 1936.17 1936.270 1.43 1.4328 0.0007 0.100 -0.0714791 -1.429584e+00 290.3683 74.749929 74.849929 symmetric 14961 0.150000
36 75.149928 1932.1607 1935.08 1933.35 1934.215 2.05 2.0543 0.0011 0.200 -0.1789349 -3.578701e+00 483.3743 74.999928 75.199928 symmetric 15031 0.349999
37 75.299928 1930.6250 1933.35 1932.81 1933.080 2.46 2.4550 0.0013 0.200 -0.3263177 -6.526361e+00 482.9426 75.199928 75.399928 symmetric 15061 0.150000
38 75.599928 1929.3783 1932.81 1930.06 1931.435 2.06 2.0567 0.0011 0.250 -0.2300109 -4.600222e+00 579.1996 75.399928 75.649928 symmetric 15121 0.300000
39 75.999928 1923.9268 1930.06 1925.30 1927.680 3.75 3.7532 0.0019 0.400 -0.2632077 -5.264159e+00 867.1924 75.649928 76.049927 symmetric 15201 0.400000
40 76.099927 1922.7347 1925.30 1923.49 1924.395 1.66 1.6603 0.0009 0.100 -0.0831850 -1.663702e+00 288.5763 76.049927 76.149927 symmetric 15221 0.099999
41 82.399921 1698.1418 1923.49 1698.82 1811.155 113.01 113.0132 0.0624 6.300 -134.0097440 -2.680197e+03 11366.8288 76.149927 82.449921 symmetric 16481 6.299994
42 83.449920 1704.2638 1704.74 1704.78 1704.760 0.50 0.4962 0.0003 0.100 -0.0246123 -4.922474e-01 255.6885 83.399920 83.499920 symmetric 16691 1.049999
43 83.549920 1703.9683 1704.78 1704.59 1704.685 0.72 0.7167 0.0004 0.100 -0.0357651 -7.153031e-01 255.6665 83.499920 83.599920 symmetric 16711 0.100000
44 187.999821 922.6182 922.90 922.87 922.885 0.27 0.2668 0.0003 0.200 -0.0248217 -4.964345e-01 230.6952 187.949821 188.149821 symmetric 37601 104.449901
45 202.549807 924.8294 925.09 925.10 925.095 0.27 0.2656 0.0003 0.150 -0.0220170 -4.403398e-01 184.9964 202.499807 202.649807 symmetric 40511 14.549986
46 248.549763 948.9234 949.23 949.24 949.235 0.31 0.3116 0.0003 0.300 -0.0511195 -1.022391e+00 332.1813 248.399763 248.699763 symmetric 49711 45.999956
47 339.499676 1272.3127 1272.73 1272.75 1272.740 0.43 0.4273 0.0003 0.200 -0.0519273 -1.038546e+00 318.1325 339.349676 339.549676 symmetric 67901 90.949913
48 339.649676 1271.8355 1272.75 1272.66 1272.705 0.87 0.8695 0.0007 0.150 -0.0527379 -1.054759e+00 254.4879 339.549676 339.699676 symmetric 67931 0.150000
49 343.599672 1321.3445 1321.98 1321.83 1321.905 0.56 0.5605 0.0004 0.150 -0.0373478 -7.469568e-01 264.3442 343.499672 343.649672 symmetric 68721 3.949996
50 344.699671 1313.2488 1313.98 1313.92 1313.950 0.70 0.7012 0.0005 0.200 -0.0773336 -1.546673e+00 328.4096 344.649671 344.849671 symmetric 68941 1.099999
51 348.699667 1365.2336 1365.61 1370.48 1368.045 2.81 2.8114 0.0021 0.300 -0.1193368 -2.386738e+00 478.6969 348.649668 348.949667 symmetric 69741 3.999996
52 356.999660 1343.7689 1344.19 1345.07 1344.630 0.86 0.8611 0.0006 0.150 -0.0548275 -1.096551e+00 268.8715 356.949660 357.099659 symmetric 71401 8.299993
53 362.649654 1664.8540 1665.57 1666.06 1665.815 0.96 0.9610 0.0006 0.100 -0.0481971 -9.639434e-01 249.8243 362.599654 362.699654 symmetric 72531 5.649994
54 362.799654 1662.5341 1666.06 1663.09 1664.575 2.04 2.0409 0.0012 0.150 -0.1098496 -2.196995e+00 332.8054 362.699654 362.849654 symmetric 72561 0.150000
55 362.899654 1661.5409 1663.09 1662.02 1662.555 1.01 1.0141 0.0006 0.100 -0.0508469 -1.016940e+00 249.3326 362.849654 362.949654 symmetric 72581 0.100000
56 362.999654 1658.5274 1662.02 1659.46 1660.740 2.21 2.2126 0.0013 0.100 -0.1107522 -2.215047e+00 249.0004 362.949654 363.049654 symmetric 72601 0.100000
57 363.099654 1656.0906 1659.46 1656.73 1658.095 2.00 2.0044 0.0012 0.100 -0.1002388 -2.004778e+00 248.6138 363.049654 363.149654 symmetric 72621 0.100000
58 423.699596 1265.2037 1274.00 1265.52 1269.760 4.56 4.5563 0.0036 1.100 -1.1371840 -2.274370e+01 1459.0838 422.749597 423.849596 symmetric 84741 60.599942
59 427.949592 1227.0050 1271.77 1243.65 1257.710 30.70 30.7050 0.0244 4.150 -51.9848038 -1.039697e+03 5230.3933 424.749595 428.899591 symmetric 85591 4.249996
60 429.899590 1234.2818 1243.65 1235.08 1239.365 5.08 5.0832 0.0041 1.105 -0.5585328 -1.117067e+01 1424.7285 428.899591 430.004590 symmetric 85981 1.949998
61 432.399588 1359.6332 1360.74 1360.62 1360.680 1.05 1.0468 0.0008 0.100 -0.0522111 -1.044222e+00 204.0492 432.349588 432.449588 symmetric 86481 2.499998
62 436.049584 1351.6060 1351.96 1352.65 1352.305 0.70 0.6990 0.0005 0.350 -0.0110272 -2.205448e-01 540.9109 435.999584 436.349584 symmetric 87211 3.649996
63 436.499584 1352.1108 1352.57 1353.12 1352.845 0.73 0.7342 0.0005 0.150 -0.0423821 -8.476420e-01 270.5261 436.449584 436.599584 symmetric 87301 0.450000
64 440.699580 1291.0250 1291.50 1292.29 1291.895 0.87 0.8700 0.0007 0.500 -0.1601788 -3.203579e+00 710.3818 440.549580 441.049579 symmetric 88141 4.199996
65 441.299579 1292.3310 1292.68 1293.55 1293.115 0.78 0.7840 0.0006 0.200 -0.0523830 -1.047661e+00 323.2271 441.249579 441.449579 symmetric 88261 0.599999
66 441.899579 1296.6441 1297.00 1300.75 1298.875 2.23 2.2309 0.0017 0.655 -0.0685342 -1.370685e+00 909.1339 441.849579 442.504578 symmetric 88381 0.600000
67 444.699576 1300.4802 1302.04 1300.94 1301.490 1.01 1.0098 0.0008 0.505 -0.2191185 -4.382374e+00 715.6009 444.299576 444.804576 symmetric 88941 2.799997
68 456.099565 1871.3335 1872.98 1872.18 1872.580 1.25 1.2465 0.0007 0.100 -0.0622608 -1.245216e+00 280.8243 456.049565 456.149565 symmetric 91221 11.399989
69 460.199561 1567.6768 1750.60 1683.17 1716.885 149.21 149.2082 0.0869 0.555 -0.2203837 -4.407678e+00 1030.0907 457.599564 458.154563 symmetric 92041 4.099996
70 499.199524 991.0502 991.33 991.30 991.315 0.26 0.2648 0.0003 0.155 -0.0208830 -4.176595e-01 198.2415 499.149524 499.304524 symmetric 99841 38.999963
71 516.299508 1020.4535 1020.71 1020.79 1020.750 0.30 0.2965 0.0003 0.200 -0.0219672 -4.393437e-01 255.1654 516.199508 516.399508 symmetric 103261 17.099984
72 529.649495 988.8677 989.62 989.17 989.395 0.53 0.5273 0.0005 0.450 -0.0682845 -1.365692e+00 494.6306 529.249495 529.699495 symmetric 105931 13.349987
73 557.799468 1615.5711 1616.00 1616.05 1616.025 0.45 0.4539 0.0003 0.200 -0.0458209 -9.164198e-01 403.9598 557.699468 557.899468 symmetric 111561 28.149973
74 558.199468 1615.9574 1616.92 1617.20 1617.060 1.10 1.1026 0.0007 0.150 -0.0873303 -1.746609e+00 323.3245 558.149468 558.299468 symmetric 111641 0.400000
75 575.899451 1151.0771 1151.59 1181.21 1166.400 15.32 15.3229 0.0131 1.600 -9.5883484 -1.917672e+02 1914.9704 575.699451 577.299449 symmetric 115181 17.699983
76 580.499446 1265.4662 1266.02 1266.44 1266.230 0.76 0.7638 0.0006 0.100 -0.0381218 -7.624366e-01 189.8960 580.449446 580.549446 symmetric 116101 4.599995
77 580.649446 1265.9426 1266.44 1266.36 1266.400 0.46 0.4574 0.0004 0.150 -0.0292346 -5.846927e-01 253.2503 580.549446 580.699446 symmetric 116131 0.150000
78 580.749446 1265.9022 1266.36 1266.70 1266.530 0.63 0.6278 0.0005 0.105 -0.0304115 -6.082297e-01 189.9480 580.699446 580.804446 symmetric 116151 0.100000
79 594.649433 1199.4276 1201.59 1199.77 1200.680 1.25 1.2524 0.0010 0.500 -0.1008521 -2.017044e+00 660.2736 594.199433 594.699433 symmetric 118931 13.899987
80 597.899430 1166.4292 1198.92 1306.59 1252.755 86.33 86.3258 0.0689 6.700 -376.3099396 -7.526206e+03 8079.7884 594.899433 601.599426 symmetric 119581 3.249997
81 605.699422 1360.0178 1361.02 1360.78 1360.900 0.88 0.8822 0.0006 0.100 -0.0441571 -8.831430e-01 204.0908 605.649422 605.749422 symmetric 121141 7.799992
82 605.849422 1360.0432 1360.78 1360.94 1360.860 0.82 0.8168 0.0006 0.200 -0.0712561 -1.425123e+00 340.1436 605.749422 605.949422 symmetric 121171 0.150000
83 651.099379 1052.3797 1076.97 1052.82 1064.895 12.52 12.5153 0.0118 2.850 -5.9697537 -1.193952e+02 3082.2229 648.549381 651.399379 symmetric 130221 45.249957
84 722.049311 1079.5794 1079.86 1079.89 1079.875 0.30 0.2956 0.0003 0.155 -0.0210824 -4.216486e-01 215.9544 721.949311 722.104311 symmetric 144411 70.949932
85 722.249311 1079.5014 1079.89 1079.83 1079.860 0.36 0.3586 0.0003 0.245 -0.0449791 -8.995838e-01 269.9215 722.104311 722.349311 symmetric 144451 0.200000
86 746.349288 1425.6879 1426.17 1430.07 1428.120 2.43 2.4321 0.0017 0.450 -0.3720176 -7.440360e+00 713.6872 746.249288 746.699288 symmetric 149271 24.099977
87 749.749285 1536.1128 1536.55 1537.42 1536.985 0.87 0.8722 0.0006 0.150 -0.0845191 -1.690384e+00 307.3127 749.649285 749.799285 symmetric 149951 3.399997
88 749.899285 1536.6967 1537.42 1537.66 1537.540 0.84 0.8433 0.0005 0.200 -0.1097109 -2.194221e+00 384.2756 749.799285 749.999285 symmetric 149981 0.150000
89 750.799284 1544.2654 1544.70 1546.21 1545.455 1.19 1.1896 0.0008 0.250 -0.0824265 -1.648531e+00 463.5529 750.749284 750.999284 symmetric 150161 0.899999
90 751.099284 1545.6555 1546.21 1548.62 1547.415 1.76 1.7595 0.0011 0.305 -0.1952344 -3.904691e+00 541.3925 750.999284 751.304283 symmetric 150221 0.300000
91 751.449283 1546.5520 1548.62 1548.44 1548.530 1.98 1.9780 0.0013 0.245 -0.1655253 -3.310508e+00 386.9705 751.304283 751.549283 symmetric 150291 0.349999
92 751.599283 1547.2750 1548.44 1548.09 1548.265 0.99 0.9900 0.0006 0.100 -0.0492785 -9.855703e-01 232.1896 751.549283 751.649283 symmetric 150321 0.150000
93 752.599282 1532.8735 1534.25 1533.62 1533.935 1.06 1.0615 0.0007 0.150 -0.0529145 -1.058291e+00 306.7341 752.499282 752.649282 symmetric 150521 0.999999
94 752.699282 1533.0397 1533.62 2067.16 1800.390 267.35 267.3503 0.1485 3.350 -269.7036300 -5.394078e+03 5851.6256 752.649282 755.999279 symmetric 150541 0.100000
95 757.899277 2212.9541 2215.37 2218.17 2216.770 3.82 3.8159 0.0017 0.300 -0.3769191 -7.538388e+00 775.4917 757.799277 758.099277 symmetric 151581 5.199995
96 758.249277 2214.7328 2218.17 2217.46 2217.815 3.08 3.0822 0.0014 0.200 -0.3708933 -7.417873e+00 554.0823 758.099277 758.299277 symmetric 151651 0.350000
97 758.349277 2214.8252 2217.46 2217.39 2217.425 2.60 2.5998 0.0012 0.100 -0.1300068 -2.600138e+00 332.4835 758.299277 758.399277 symmetric 151671 0.100000
98 758.599277 2212.9240 2215.85 2214.59 2215.220 2.30 2.2960 0.0010 0.100 -0.1146339 -2.292681e+00 332.1676 758.549277 758.649276 symmetric 151721 0.250000
99 758.699276 2213.5525 2214.59 2214.71 2214.650 1.10 1.0975 0.0005 0.150 -0.0694234 -1.388468e+00 442.8602 758.649276 758.799276 symmetric 151741 0.099999
100 758.849276 2212.0906 2214.71 2213.43 2214.070 1.98 1.9794 0.0009 0.100 -0.0992100 -1.984203e+00 332.0117 758.799276 758.899276 symmetric 151771 0.150000
101 759.099276 2208.0959 2212.56 2209.11 2210.835 2.74 2.7391 0.0012 0.150 -0.2003221 -4.006446e+00 441.9658 758.999276 759.149276 symmetric 151821 0.250000
102 760.249275 2173.4627 2202.55 2174.30 2188.425 14.96 14.9623 0.0068 0.850 -0.0260366 -5.207318e-01 1969.5549 759.449276 760.299275 symmetric 152051 1.149999
103 765.249270 1726.2590 2029.61 1849.76 1939.685 213.43 213.4260 0.1100 1.500 0.7485108 1.497023e+01 3007.2569 762.099273 763.599272 symmetric 153051 4.999995
104 765.599270 1716.1664 1726.88 1717.95 1722.415 6.25 6.2486 0.0036 0.350 -0.6345536 -1.269108e+01 688.3293 765.299270 765.649270 symmetric 153121 0.350000
105 765.699270 1715.4809 1717.95 1715.92 1716.935 1.45 1.4541 0.0008 0.100 -0.0725313 -1.450627e+00 257.4670 765.649270 765.749270 symmetric 153141 0.100000
106 765.799270 1714.8943 1715.92 1715.37 1715.645 0.75 0.7507 0.0004 0.100 -0.0373800 -7.476012e-01 257.3087 765.749270 765.849270 symmetric 153161 0.100000
107 774.349262 2071.1613 2071.75 2076.73 2074.240 3.08 3.0787 0.0015 0.200 -0.1704033 -3.408069e+00 518.3892 774.299262 774.499261 symmetric 154871 8.549992
108 774.549261 2075.6861 2076.73 2085.83 2081.280 5.59 5.5939 0.0027 0.500 -0.3135638 -6.271282e+00 1144.3904 774.499261 774.999261 symmetric 154911 0.199999
109 775.299261 2087.5887 2088.42 2090.20 2089.310 1.72 1.7213 0.0008 0.150 -0.0780186 -1.560374e+00 417.7834 775.249261 775.399261 symmetric 155061 0.750000
110 775.899260 2094.0998 2095.26 2095.51 2095.385 1.29 1.2852 0.0006 0.150 -0.0954006 -1.908014e+00 418.9821 775.849260 775.999260 symmetric 155181 0.599999
111 776.099260 2094.8245 2095.51 2095.82 2095.665 0.84 0.8405 0.0004 0.150 -0.0788284 -1.576569e+00 419.0545 775.999260 776.149260 symmetric 155221 0.200000
112 776.349260 2094.6264 2095.51 2095.55 2095.530 0.90 0.9036 0.0004 0.100 -0.0452196 -9.043924e-01 314.2841 776.299260 776.399260 symmetric 155271 0.250000
113 776.449260 2094.2813 2095.55 2095.09 2095.320 1.04 1.0387 0.0005 0.100 -0.0519480 -1.038961e+00 314.2458 776.399260 776.499259 symmetric 155291 0.100000
114 776.649259 2092.0468 2095.09 2092.82 2093.955 1.91 1.9082 0.0009 0.200 -0.1003538 -2.007079e+00 523.3882 776.499259 776.699259 symmetric 155331 0.199999
115 778.649257 2052.7332 2056.10 2054.26 2055.180 2.45 2.4468 0.0012 0.200 -0.2787151 -5.574306e+00 513.5159 778.499258 778.699257 symmetric 155731 1.999998
116 778.754257 2050.9462 2054.26 2052.75 2053.505 2.56 2.5588 0.0012 0.200 -0.3000249 -6.000504e+00 513.0752 778.699257 778.899257 symmetric 155752 0.105000
117 779.099257 2050.1338 2051.58 2051.70 2051.640 1.51 1.5062 0.0007 0.100 -0.0751317 -1.502636e+00 307.6700 779.049257 779.149257 symmetric 155821 0.345000
118 779.699256 2043.3876 2044.88 2044.29 2044.585 1.20 1.1974 0.0006 0.150 -0.0807806 -1.615614e+00 408.8356 779.599257 779.749256 symmetric 155941 0.599999
119 780.249256 2036.4888 2041.49 2037.12 2039.305 2.82 2.8161 0.0014 0.350 -0.3831781 -7.663569e+00 815.3385 779.949256 780.299256 symmetric 156051 0.550000
120 781.199255 2021.5449 2024.50 2022.46 2023.480 1.94 1.9351 0.0010 0.100 -0.0967932 -1.935866e+00 303.4250 781.149255 781.249255 symmetric 156241 0.949999
121 781.549255 2021.2170 2022.96 2022.23 2022.595 1.38 1.3780 0.0007 0.200 -0.1454133 -2.908269e+00 505.5033 781.449255 781.649255 symmetric 156311 0.350000
122 781.699255 2019.5762 2022.23 2021.78 2022.005 2.43 2.4288 0.0012 0.100 -0.1216647 -2.433295e+00 303.1795 781.649255 781.749254 symmetric 156341 0.150000
123 781.799254 2020.6270 2021.78 2024.07 2022.925 2.30 2.2980 0.0011 0.200 -0.1516451 -3.032904e+00 505.5797 781.749254 781.949254 symmetric 156361 0.099999
124 782.049254 2022.5590 2024.07 2025.93 2025.000 2.44 2.4410 0.0012 0.150 -0.1813020 -3.626044e+00 404.8183 781.949254 782.099254 symmetric 156411 0.250000
125 782.399254 2025.6619 2026.30 2027.78 2027.040 1.38 1.3781 0.0007 0.105 -0.0643518 -1.287038e+00 303.9859 782.349254 782.454254 symmetric 156481 0.350000
126 782.549254 2026.8180 2027.78 2031.40 2029.590 2.77 2.7720 0.0014 0.245 -0.1884004 -3.768011e+00 507.1254 782.454254 782.699254 symmetric 156511 0.150000
127 783.249253 2036.2432 2037.07 2038.83 2037.950 1.71 1.7068 0.0008 0.100 -0.0853620 -1.707242e+00 305.6069 783.199253 783.299253 symmetric 156651 0.699999
128 783.349253 2037.1402 2038.83 2037.79 2038.310 1.17 1.1698 0.0006 0.100 -0.0584978 -1.169957e+00 305.6877 783.299253 783.399253 symmetric 156671 0.100000
129 783.549253 2035.8336 2038.11 2037.55 2037.830 2.00 1.9964 0.0010 0.100 -0.0997569 -1.995139e+00 305.5743 783.499253 783.599253 symmetric 156711 0.200000
130 783.699253 2036.4960 2037.55 2038.90 2038.225 1.73 1.7290 0.0008 0.150 -0.1690347 -3.380696e+00 407.4758 783.599253 783.749253 symmetric 156741 0.150000
131 788.449248 1829.9721 1833.60 1830.90 1832.250 2.28 2.2779 0.0012 0.150 -0.1247051 -2.494104e+00 366.3247 788.349248 788.499248 symmetric 157691 4.749995
132 788.649248 1825.3076 1830.90 1825.86 1828.380 3.07 3.0724 0.0017 0.200 -0.2738726 -5.477458e+00 456.8205 788.499248 788.699248 symmetric 157731 0.200000
133 788.749248 1824.0863 1825.86 1825.29 1825.575 1.49 1.4887 0.0008 0.100 -0.0744490 -1.488982e+00 273.7616 788.699248 788.799248 symmetric 157751 0.100000
134 789.049248 1822.5677 1825.29 1824.42 1824.855 2.29 2.2873 0.0013 0.350 -0.4344932 -8.689872e+00 729.5065 788.799248 789.149247 symmetric 157811 0.300000
135 789.249247 1822.7308 1824.42 1823.80 1824.110 1.38 1.3792 0.0008 0.200 -0.1738120 -3.476243e+00 455.8531 789.149247 789.349247 symmetric 157851 0.199999
136 789.599247 1824.3551 1826.04 1826.89 1826.465 2.11 2.1099 0.0012 0.155 -0.1610763 -3.221530e+00 365.1291 789.549247 789.704247 symmetric 157921 0.350000
137 790.149246 1832.1430 1832.97 1834.73 1833.850 1.71 1.7070 0.0009 0.105 -0.0792432 -1.584866e+00 274.9914 790.099247 790.204246 symmetric 158031 0.549999
138 790.449246 1837.8449 1839.63 1839.80 1839.715 1.87 1.8701 0.0010 0.100 -0.0935187 -1.870375e+00 275.8635 790.399246 790.499246 symmetric 158091 0.300000
139 790.549246 1839.0726 1839.80 1843.11 1841.455 2.38 2.3824 0.0013 0.250 -0.0125202 -2.504051e-01 552.4228 790.499246 790.749246 symmetric 158111 0.100000
140 790.799246 1841.8582 1843.11 1843.52 1843.315 1.46 1.4568 0.0008 0.100 -0.0726953 -1.453908e+00 276.4239 790.749246 790.849246 symmetric 158161 0.250000
141 792.799244 1806.4952 1809.06 1807.82 1808.440 1.94 1.9448 0.0011 0.150 -0.1370422 -2.740846e+00 361.5503 792.699244 792.849244 symmetric 158561 1.999998
142 792.949244 1803.7963 1807.82 1804.91 1806.365 2.57 2.5687 0.0014 0.200 -0.2348460 -4.696924e+00 451.3558 792.849244 793.049244 symmetric 158591 0.150000
143 793.099244 1803.1059 1804.91 1805.03 1804.970 1.86 1.8641 0.0010 0.100 -0.0931711 -1.863424e+00 270.6520 793.049244 793.149244 symmetric 158621 0.150000
144 793.299243 1801.7727 1805.03 1805.16 1805.095 3.32 3.3223 0.0018 0.250 -0.3603361 -7.206728e+00 541.1667 793.149244 793.399243 symmetric 158661 0.199999
145 793.499243 1804.2160 1805.16 1805.19 1805.175 0.96 0.9590 0.0005 0.200 -0.0957290 -1.914582e+00 451.1971 793.399243 793.599243 symmetric 158701 0.200000
146 793.699243 1804.0024 1805.19 1804.49 1804.840 0.84 0.8376 0.0005 0.150 -0.0653549 -1.307099e+00 360.9020 793.599243 793.749243 symmetric 158741 0.200000
147 793.799243 1803.6362 1804.49 1804.54 1804.515 0.88 0.8788 0.0005 0.100 -0.0438978 -8.779561e-01 270.6330 793.749243 793.849243 symmetric 158761 0.100000
148 794.149243 1802.0331 1803.30 1803.21 1803.255 1.22 1.2219 0.0007 0.150 -0.1145858 -2.291719e+00 360.5361 794.049243 794.199243 symmetric 158831 0.350000
149 797.149240 1704.8008 1725.62 1705.83 1715.725 10.92 10.9242 0.0064 0.800 -1.6598486 -3.319700e+01 1456.7045 796.399240 797.199240 symmetric 159431 2.999997
150 797.649239 1699.9652 1705.83 1700.54 1703.185 3.22 3.2198 0.0019 0.500 -0.5692084 -1.138418e+01 936.1797 797.199240 797.699239 symmetric 159531 0.499999
151 797.799239 1699.0384 1700.54 1700.72 1700.630 1.59 1.5916 0.0009 0.250 -0.1642746 -3.285495e+00 510.0230 797.699239 797.949239 symmetric 159561 0.150000
152 798.049239 1698.2809 1700.72 1701.33 1701.025 2.74 2.7441 0.0016 0.150 -0.1994695 -3.989393e+00 340.0052 797.949239 798.099239 symmetric 159611 0.250000
153 798.149239 1698.9092 1701.33 1700.55 1700.940 2.03 2.0308 0.0012 0.150 -0.1429534 -2.859070e+00 340.0453 798.099239 798.249239 symmetric 159631 0.100000
154 798.299239 1698.8642 1700.55 1699.85 1700.200 1.34 1.3358 0.0008 0.150 -0.1029878 -2.059758e+00 339.9373 798.249239 798.399239 symmetric 159661 0.150000
155 798.599238 1697.6356 1698.58 1699.16 1698.870 1.23 1.2344 0.0007 0.250 -0.1527820 -3.055642e+00 509.5083 798.499238 798.749238 symmetric 159721 0.299999
156 798.949238 1696.9947 1698.04 1698.37 1698.205 1.21 1.2103 0.0007 0.200 -0.1204950 -2.409903e+00 424.4302 798.849238 799.049238 symmetric 159791 0.350000
157 799.099238 1697.0959 1698.37 1699.14 1698.755 1.66 1.6591 0.0010 0.150 -0.0963066 -1.926135e+00 339.6543 799.049238 799.199238 symmetric 159821 0.150000
158 799.349238 1698.1512 1699.38 1698.80 1699.090 0.94 0.9388 0.0006 0.100 -0.0469382 -9.387642e-01 254.8163 799.299238 799.399238 symmetric 159871 0.250000
159 799.449238 1698.2526 1698.80 1699.97 1699.385 1.13 1.1324 0.0007 0.200 -0.0352854 -7.057088e-01 424.8099 799.399238 799.599237 symmetric 159891 0.100000
160 799.649237 1699.1575 1699.97 1699.90 1699.935 0.78 0.7775 0.0005 0.100 -0.0386441 -7.728825e-01 254.9507 799.599237 799.699237 symmetric 159931 0.199999
161 799.799237 1699.1741 1699.90 1699.95 1699.925 0.75 0.7509 0.0004 0.150 -0.0579509 -1.159019e+00 339.9263 799.699237 799.849237 symmetric 159961 0.150000
162 799.899237 1698.1964 1699.95 1699.45 1699.700 1.50 1.5036 0.0009 0.150 -0.0975289 -1.950579e+00 339.8424 799.849237 799.999237 symmetric 159981 0.100000
163 800.099237 1697.8467 1699.45 1698.57 1699.010 1.16 1.1633 0.0007 0.155 -0.1068678 -2.137357e+00 339.6980 799.999237 800.154237 symmetric 160021 0.200000
164 805.199232 1516.6081 1517.12 1610.46 1563.790 47.18 47.1819 0.0302 2.150 -25.9011508 -5.180235e+02 3414.4378 805.149232 807.299230 symmetric 161041 5.099995
165 809.899228 1724.8404 1726.31 1725.91 1726.110 1.27 1.2696 0.0007 0.105 -0.0651301 -1.302604e+00 258.8523 809.849228 809.954228 symmetric 161981 4.699996
166 810.349227 1724.8191 1725.36 1726.19 1725.775 0.96 0.9559 0.0006 0.250 -0.0930265 -1.860531e+00 517.6389 810.249227 810.499227 symmetric 162071 0.449999
167 810.549227 1725.0283 1726.19 1727.74 1726.965 1.94 1.9367 0.0011 0.100 -0.0966206 -1.932414e+00 258.9472 810.499227 810.599227 symmetric 162111 0.200000
168 810.649227 1726.1122 1727.74 1729.59 1728.665 2.55 2.5528 0.0015 0.350 -0.2190497 -4.380998e+00 691.2444 810.599227 810.949227 symmetric 162131 0.100000
169 814.349223 1597.1787 1601.23 1613.70 1607.465 10.29 10.2863 0.0064 1.400 -6.6656696 -1.333135e+02 2324.1549 813.749224 815.149223 symmetric 162871 3.699996
170 818.449219 1678.0575 1678.52 1678.54 1678.530 0.47 0.4725 0.0003 0.100 -0.0238151 -4.763020e-01 251.7560 818.399220 818.499219 symmetric 163691 4.099996
171 820.954217 1592.3320 1792.51 1800.50 1796.505 204.17 204.1730 0.1136 0.050 0.0000000 0.000000e+00 179.6504 823.249215 823.299215 symmetric 164192 2.504998
172 827.999210 1764.6976 1886.14 1765.68 1825.910 61.21 61.2124 0.0335 1.500 -1.2108750 -2.421752e+01 2828.9438 826.549212 828.049210 symmetric 165601 7.044993
173 828.599210 1756.9386 1765.68 1757.48 1761.580 4.64 4.6414 0.0026 0.600 -0.7490487 -1.498099e+01 1144.2775 828.049210 828.649210 symmetric 165721 0.600000
174 828.749210 1756.0013 1757.48 1804.82 1781.150 25.15 25.1487 0.0141 0.800 -9.1495772 -1.829917e+02 1504.8270 828.649210 829.449209 symmetric 165751 0.150000
175 831.749207 1871.6279 1873.33 1872.24 1872.785 1.16 1.1571 0.0006 0.250 -0.1266671 -2.533345e+00 561.7084 831.549207 831.799207 symmetric 166351 2.999997
176 832.099206 1867.8105 1872.24 1868.53 1870.385 2.57 2.5745 0.0014 0.350 -0.1227567 -2.455135e+00 748.0309 831.799207 832.149206 symmetric 166421 0.349999
177 832.199206 1868.0587 1868.53 1922.15 1895.340 27.28 27.2813 0.0144 1.105 0.4793839 9.587688e+00 2179.9790 832.149206 833.254205 symmetric 166441 0.100000
178 833.549205 1920.9525 1922.55 1922.00 1922.275 1.32 1.3225 0.0007 0.250 -0.1128484 -2.256970e+00 576.5695 833.399205 833.649205 symmetric 166711 1.349999
179 833.699205 1921.3389 1922.00 1924.09 1923.045 1.71 1.7061 0.0009 0.200 -0.1726294 -3.452592e+00 480.5882 833.649205 833.849205 symmetric 166741 0.150000
180 834.049205 1920.6909 1924.09 1921.62 1922.855 2.16 2.1641 0.0011 0.250 -0.2404220 -4.808445e+00 576.6152 833.849205 834.099205 symmetric 166811 0.350000
181 834.549204 1919.4376 1920.15 1921.92 1921.035 1.60 1.5974 0.0008 0.200 -0.1359455 -2.718912e+00 480.1227 834.499204 834.699204 symmetric 166911 0.499999
182 835.349203 1909.4604 1921.92 2004.89 1963.405 53.94 53.9447 0.0275 1.800 -38.4933599 -7.698679e+02 3593.8007 834.699204 836.499202 symmetric 167071 0.799999
183 837.449201 1928.5756 1929.66 1991.77 1960.715 32.14 32.1394 0.0164 0.850 0.0965852 1.931706e+00 1764.7364 837.399201 838.249201 symmetric 167491 2.099998
184 839.449199 1887.2367 1887.89 1889.50 1888.695 1.46 1.4583 0.0008 0.105 -0.0763654 -1.527310e+00 283.2220 839.399199 839.504199 symmetric 167891 1.999998
185 841.299198 2023.8566 2027.42 2027.20 2027.310 3.45 3.4534 0.0017 0.200 -0.3055088 -6.110183e+00 506.5213 841.149198 841.349198 symmetric 168261 1.849999
186 842.749196 1973.2640 1974.61 1974.22 1974.415 1.15 1.1510 0.0006 0.100 -0.0576444 -1.152890e+00 296.1046 842.699196 842.799196 symmetric 168551 1.449998
187 850.899189 1523.7289 1524.33 1524.13 1524.230 0.50 0.5011 0.0003 0.100 -0.0249906 -4.998122e-01 228.6091 850.849189 850.949188 symmetric 170181 8.149993
188 851.099188 1523.0616 1524.13 1523.92 1524.025 0.96 0.9634 0.0006 0.300 -0.1439722 -2.879447e+00 533.2642 850.949188 851.249188 symmetric 170221 0.199999
189 851.299188 1523.0474 1523.92 1523.47 1523.695 0.65 0.6476 0.0004 0.100 -0.0322560 -6.451209e-01 228.5214 851.249188 851.349188 symmetric 170261 0.200000
190 851.449188 1522.6689 1523.47 1523.41 1523.440 0.77 0.7711 0.0005 0.150 -0.0557628 -1.115258e+00 304.6320 851.349188 851.499188 symmetric 170291 0.150000
191 852.499187 1541.2544 1542.69 1541.72 1542.205 0.95 0.9506 0.0006 0.150 -0.0418036 -8.360732e-01 308.3986 852.399187 852.549187 symmetric 170501 1.049999
192 853.349186 1521.5502 1523.85 1569.13 1546.490 24.94 24.9398 0.0161 1.650 -5.6688312 -1.133767e+02 2623.3643 853.149186 854.799185 symmetric 170671 0.849999
193 854.849185 1566.8404 1569.13 1567.98 1568.555 1.71 1.7146 0.0011 0.100 -0.0855738 -1.711477e+00 235.1970 854.799185 854.899185 symmetric 170971 1.499999
194 854.999185 1566.0856 1567.98 1566.77 1567.375 1.29 1.2894 0.0008 0.150 -0.0479615 -9.592311e-01 313.4260 854.899185 855.049185 symmetric 171001 0.150000
195 855.099185 1566.1740 1566.77 1566.68 1566.725 0.55 0.5510 0.0004 0.100 -0.0275399 -5.507986e-01 234.9810 855.049185 855.149184 symmetric 171021 0.100000
196 855.649184 1558.3609 1566.68 1559.78 1563.230 4.87 4.8691 0.0031 0.600 -1.0395412 -2.079084e+01 1015.0603 855.149184 855.749184 symmetric 171131 0.549999
197 855.799184 1559.3120 1559.78 1560.53 1560.155 0.84 0.8430 0.0005 0.150 -0.0288797 -5.775936e-01 312.0017 855.749184 855.899184 symmetric 171161 0.150000
198 856.699183 1549.9519 1555.72 1580.27 1567.995 18.04 18.0431 0.0115 1.250 -6.8447804 -1.368957e+02 2031.5461 856.399183 857.649182 symmetric 171341 0.899999
199 857.749182 1579.3998 1580.27 1579.90 1580.085 0.69 0.6852 0.0004 0.150 -0.0475080 -9.501610e-01 315.9690 857.649182 857.799182 symmetric 171551 1.049999
200 858.799181 1554.6868 1555.39 1555.24 1555.315 0.63 0.6282 0.0004 0.100 -0.0313225 -6.264499e-01 233.2654 858.749181 858.849181 symmetric 171761 1.049999
201 858.949181 1553.8428 1555.24 1555.13 1555.185 1.34 1.3422 0.0009 0.200 -0.0893165 -1.786331e+00 388.7068 858.849181 859.049181 symmetric 171791 0.150000
202 859.149181 1554.4359 1555.13 1554.89 1555.010 0.57 0.5741 0.0004 0.155 -0.0406500 -8.130014e-01 310.9624 859.049181 859.204181 symmetric 171831 0.200000
203 859.649180 1555.4485 1556.09 1556.14 1556.115 0.67 0.6665 0.0004 0.150 -0.0442472 -8.849439e-01 311.1788 859.549180 859.699180 symmetric 171931 0.499999
204 859.849180 1553.8066 1556.14 1554.21 1555.175 1.37 1.3684 0.0009 0.200 -0.0766084 -1.532169e+00 388.7172 859.699180 859.899180 symmetric 171971 0.200000
205 860.849179 1537.4271 1551.02 1538.78 1544.900 7.47 7.4730 0.0048 0.750 -0.2164051 -4.328106e+00 1235.6998 860.149180 860.899179 symmetric 172171 0.999999
206 860.999179 1538.1084 1538.78 1538.77 1538.775 0.67 0.6666 0.0004 0.150 -0.0557919 -1.115839e+00 307.6984 860.899179 861.049179 symmetric 172201 0.150000
207 861.149179 1537.9571 1538.77 1540.42 1539.595 1.64 1.6379 0.0011 0.200 -0.1725063 -3.450129e+00 384.7259 861.049179 861.249179 symmetric 172231 0.150000
208 861.299179 1540.0210 1540.42 1548.55 1544.485 4.46 4.4640 0.0029 0.500 0.0365640 7.312810e-01 849.5023 861.249179 861.749178 symmetric 172261 0.150000
209 861.999178 1550.7634 1551.25 1574.32 1562.785 12.02 12.0216 0.0077 0.600 -1.9173247 -3.834653e+01 1013.8913 861.949178 862.549177 symmetric 172401 0.699999
210 864.649175 1691.7105 1692.70 1692.36 1692.530 0.82 0.8195 0.0005 0.105 -0.0423759 -8.475188e-01 253.8386 864.599175 864.704175 symmetric 172931 2.649997
211 864.999175 1693.9125 1695.49 1699.23 1697.360 3.45 3.4475 0.0020 0.255 -0.4733365 -9.466739e+00 508.7243 864.899175 865.154175 symmetric 173001 0.350000
212 865.699174 1713.2966 1714.96 1718.01 1716.485 3.19 3.1884 0.0019 0.155 -0.1786644 -3.573292e+00 343.1078 865.649174 865.804174 symmetric 173141 0.699999
213 865.949174 1719.3749 1719.99 1724.40 1722.195 2.82 2.8201 0.0016 0.255 -0.0579052 -1.158105e+00 516.5864 865.899174 866.154174 symmetric 173191 0.250000
214 868.049172 1799.1300 1804.19 1799.90 1802.045 2.91 2.9150 0.0016 0.150 -0.1822885 -3.645774e+00 360.2267 867.949172 868.099172 symmetric 173611 2.099998
215 869.199171 1748.6448 1750.00 1752.14 1751.070 2.43 2.4252 0.0014 0.200 -0.2154842 -4.309688e+00 437.5513 869.149171 869.349171 symmetric 173841 1.149999
216 869.449171 1751.5620 1752.14 1757.57 1754.855 3.29 3.2930 0.0019 0.350 -0.5021464 -1.004294e+01 701.4388 869.349171 869.699171 symmetric 173891 0.250000
217 869.749171 1756.8646 1757.57 1759.85 1758.710 1.85 1.8454 0.0011 0.150 -0.0801521 -1.603044e+00 351.6611 869.699171 869.849170 symmetric 173951 0.300000
218 869.949170 1758.5563 1759.85 1759.28 1759.565 1.01 1.0087 0.0006 0.150 -0.0610982 -1.221965e+00 351.8519 869.849170 869.999170 symmetric 173991 0.199999
219 870.099170 1757.5467 1759.28 1758.21 1758.745 1.20 1.1983 0.0007 0.150 -0.0958054 -1.916110e+00 351.6528 869.999170 870.149170 symmetric 174021 0.150000
220 876.149164 1388.6454 1632.86 1514.70 1573.780 185.13 185.1346 0.1176 6.850 -747.6874313 -1.495376e+04 10111.3898 871.949168 878.799162 symmetric 175231 6.049994
221 883.049158 1560.3538 1624.99 1587.94 1606.465 46.11 46.1112 0.0287 2.650 -57.3989068 -1.147979e+03 4280.0508 881.399159 884.049157 symmetric 176611 6.899994
222 885.899155 1537.7235 1567.28 1539.06 1553.170 15.45 15.4465 0.0099 1.200 -4.1453533 -8.290714e+01 1937.3149 884.849156 886.049155 symmetric 177181 2.849997
223 886.099155 1538.5704 1539.06 1542.38 1540.720 2.15 2.1496 0.0014 0.200 -0.0736493 -1.472988e+00 385.1063 886.049155 886.249155 symmetric 177221 0.200000
224 886.299155 1540.9431 1542.38 1542.56 1542.470 1.53 1.5269 0.0010 0.100 -0.0764080 -1.528161e+00 231.2941 886.249155 886.349155 symmetric 177261 0.200000
225 887.649153 1506.2289 1542.56 1506.81 1524.685 18.46 18.4561 0.0121 1.350 -3.9374706 -7.874949e+01 2130.6193 886.349155 887.699153 symmetric 177531 1.349998
226 888.999152 1534.4181 1535.04 1536.85 1535.945 1.53 1.5269 0.0010 0.200 -0.0448910 -8.978201e-01 383.9407 888.949152 889.149152 symmetric 177801 1.349999
227 889.299152 1535.3062 1536.85 1536.74 1536.795 1.49 1.4888 0.0010 0.200 -0.1335390 -2.670783e+00 384.0639 889.149152 889.349152 symmetric 177861 0.300000
228 889.549152 1533.7295 1536.74 1534.18 1535.460 1.73 1.7305 0.0011 0.250 -0.1728760 -3.457524e+00 460.4640 889.349152 889.599152 symmetric 177911 0.250000
229 890.049151 1528.8343 1534.18 1537.95 1536.065 7.23 7.2307 0.0047 1.050 -3.2183234 -6.436653e+01 1686.4497 889.599152 890.649151 symmetric 178011 0.499999
230 892.499149 1431.4406 1431.99 1439.27 1435.630 4.19 4.1894 0.0029 0.650 -0.4983086 -9.966181e+00 1004.4404 892.449149 893.099148 symmetric 178501 2.449998
231 897.999144 1354.3801 1358.86 1354.82 1356.840 2.46 2.4599 0.0018 0.700 -0.3730899 -7.461805e+00 1017.2588 897.449144 898.149143 symmetric 179601 5.499995

Using a reference

Instead of detecting sunflecks as a dynamic temporal process, or local changes in slope, especially in the case of cloud flecks an alternative is to use the estimated clear-sky irradiance as reference and to quantify the size of the deviations, both increases and decreases, in many cases using a fixed threshold. An alternative to using a clear sky simulation is to use a smoothed version of the same data as a reference, a moving median, or a spline-based smoother. However, each of these and the temporal-process approach described above are fundamentally different and describe different properties of the data.

Both observed and reference irradiances need to be available at the same time points. This in some cases can make it necessary to use interpolation to obtain suitable reference data. After this, a simple difference at each time point can be computed and compared against a threshold. As the reference varies in time, the threshold used can be expressed either as an absolute difference or as a relative difference, thus, giving two different quantifications.

Using a quantile smoother

Quantile regression, either of a linear model or a smmothing spline, differs from OLS regression and other common approaches in that it can flexibly be used to find a fitted line that describes a specific quantile. It can be imagined as a combination of a boxplot and a smoother. This approach lets us estimate local deviations from an internal reference different from the local mean or median, say, for example, to identify when irradiance is less than that in the 10% locally highest irradiances.

For vegetation shade a further possibility is to use a quantile regression directly fitted for a target shade threshold.