diff --git a/analyses/all b/analyses/all
deleted file mode 100644
index ceae1b0..0000000
--- a/analyses/all
+++ /dev/null
@@ -1,651 +0,0 @@
-==> g(𝛌＝4775)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--9.745e-01 2.756e-01
--1.384e+00 2.128e-01
--2.748e+00 2.454e-01
--3.305e+00 2.335e-01
--3.314e+00 1.777e-01
--3.389e+00 1.613e-01
--4.198e+00 1.693e-01
--4.465e+00 1.915e-01
--4.700e-01 1.237e+00
--7.487e-01 1.302e+00
--2.046e+00 1.348e+00
--2.428e+00 1.335e+00
--2.953e+00 1.392e+00
--3.086e+00 1.396e+00
--3.761e+00 1.443e+00
--4.290e+00 1.408e+00
-9.862e-02 3.790e-01
-3.899e-01 4.121e-01
-8.650e-01 4.290e-01
-5.516e-01 4.096e-01
-2.228e-01 4.090e-01
-9.508e-01 4.215e-01
--2.872e-01 3.988e-01
-9.059e-02 3.866e-01
-***************************************
-
-==> HST(𝛌＝1158)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
-2.402e-01 2.838e-01
-1.329e-01 2.104e-01
--8.990e-01 2.184e-01
--1.229e+00 1.729e-01
--1.825e+00 1.514e-01
--2.236e+00 1.569e-01
--2.758e+00 1.438e-01
--4.749e+00 3.931e+00
-1.518e-01 3.556e-02
-1.794e-02 4.241e-03
--1.040e+00 2.554e-01
--1.377e+00 3.397e-01
--1.989e+00 4.898e-01
--2.375e+00 5.854e-01
--2.886e+00 7.127e-01
--3.737e+00 9.211e-01
--3.034e-02 7.945e-03
--6.855e-02 1.742e-02
--6.829e-02 1.598e-02
--4.907e-02 1.126e-02
--5.681e-02 1.297e-02
--1.644e-01 3.816e-02
--3.273e-01 8.072e-02
-1.333e-01 3.447e-02
-***************************************
-
-==> HST(𝛌＝1367)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
-6.534e-02 2.733e-01
--9.689e-02 2.184e-01
--1.175e+00 2.494e-01
--1.525e+00 1.446e-01
--2.167e+00 1.617e-01
--2.492e+00 1.545e-01
--3.258e+00 1.095e-01
--1.115e+01 5.296e+03
-6.378e-02 1.554e-02
--9.676e-02 2.400e-02
--1.181e+00 2.892e-01
--1.520e+00 3.738e-01
--2.180e+00 5.369e-01
--2.532e+00 6.061e-01
--3.129e+00 7.691e-01
--3.794e+00 9.654e-01
--8.535e-05 2.174e-05
--5.126e-04 1.233e-04
--4.066e-04 9.595e-05
-1.080e-03 2.497e-04
--2.682e-03 6.077e-04
--5.462e-03 1.334e-03
--2.209e-02 5.376e-03
-1.564e-01 3.945e-02
-***************************************
-
-==> HST(𝛌＝1479)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
-1.224e-01 2.666e-01
--7.607e-02 2.200e-01
--1.174e+00 2.513e-01
--1.500e+00 1.446e-01
--2.075e+00 1.576e-01
--2.554e+00 1.536e-01
--3.337e+00 1.052e-01
--1.152e+01 9.495e+03
-9.340e-02 2.202e-02
--8.842e-02 2.148e-02
--1.173e+00 2.885e-01
--1.512e+00 3.787e-01
--2.113e+00 5.217e-01
--2.547e+00 6.151e-01
--3.178e+00 7.800e-01
--7.000e+00 1.520e+00
--6.035e-02 1.487e-02
-2.044e-02 5.014e-03
-1.258e-01 3.039e-02
-5.764e-02 1.497e-02
-2.580e-02 6.105e-03
--4.737e-02 1.144e-02
-9.579e-02 2.320e-02
-6.411e-01 1.538e-01
-***************************************
-
-==> HST(𝛌＝1746)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
-3.934e-02 2.650e-01
--1.638e-01 2.223e-01
--1.339e+00 2.641e-01
--1.655e+00 1.414e-01
--2.232e+00 1.613e-01
--2.793e+00 1.845e-01
--4.025e+00 1.005e-01
--1.072e+01 3.151e+03
-5.043e-02 1.188e-02
--1.304e-01 3.177e-02
--1.254e+00 3.099e-01
--1.587e+00 4.030e-01
--2.207e+00 5.484e-01
--2.693e+00 6.488e-01
--3.209e+00 8.109e-01
--3.820e+00 9.624e-01
-3.040e-02 8.141e-03
-2.438e-02 6.215e-03
-1.450e-01 3.324e-02
-1.458e-01 3.170e-02
-1.772e-01 3.882e-02
-2.076e-01 4.966e-02
--1.390e-01 3.462e-02
-3.929e-01 1.097e-01
-***************************************
-
-==> i(𝛌＝7647)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.510e+00 2.756e-01
--1.872e+00 2.087e-01
--3.441e+00 2.422e-01
--3.722e+00 1.931e-01
--4.794e+00 3.059e-01
--4.336e+00 1.760e-01
--4.644e+00 1.360e-01
--5.061e+00 1.844e-01
--7.402e-01 1.117e+00
--9.967e-01 1.174e+00
--2.393e+00 1.214e+00
--2.723e+00 1.271e+00
--3.701e+00 1.349e+00
--3.742e+00 1.343e+00
--4.317e+00 1.393e+00
--4.979e+00 1.451e+00
-2.638e-01 3.565e-01
-7.095e-01 3.780e-01
-1.297e+00 4.108e-01
-1.035e+00 3.960e-01
-9.267e-01 3.966e-01
-8.982e-01 4.049e-01
-7.318e-01 3.899e-01
-9.092e-01 3.982e-01
-***************************************
-
-==> I(𝛌＝8560)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.497e+00 3.256e-01
--1.786e+00 2.808e-01
--3.555e+00 5.022e-01
--4.203e+00 1.412e-01
--3.826e+00 2.505e-01
--4.096e+00 5.759e-01
--4.304e+00 1.379e-01
--1.109e+01 5.792e+03
--7.542e-01 1.301e+00
--9.532e-01 1.415e+00
--2.397e+00 1.461e+00
--3.214e+00 1.450e+00
--3.008e+00 1.544e+00
--3.489e+00 1.514e+00
--4.353e+00 1.499e+00
--4.061e+00 1.520e+00
-1.675e-01 4.363e-01
-7.660e-01 4.491e-01
-1.781e+00 4.834e-01
-1.971e+00 4.573e-01
-6.572e-01 4.531e-01
-8.520e-01 4.811e-01
-2.405e+00 4.975e-01
--1.573e-01 4.341e-01
-***************************************
-
-==> plots <==
-
-==> r(𝛌＝6175)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.372e+00 2.752e-01
--1.786e+00 2.125e-01
--3.214e+00 2.522e-01
--3.558e+00 2.087e-01
--4.416e+00 3.027e-01
--4.366e+00 2.391e-01
--4.559e+00 1.732e-01
--4.586e+00 1.687e-01
--6.693e-01 4.559e-01
--5.998e+00 1.480e+00
--2.390e+00 9.034e-01
--2.927e+00 1.017e+00
--3.411e+00 1.073e+00
--3.731e+00 1.204e+00
--3.930e+00 1.152e+00
--7.000e+00 1.390e+00
-2.905e-01 1.874e-01
-1.835e+00 4.578e-01
-8.760e-01 3.014e-01
-1.476e+00 3.951e-01
--3.280e-01 2.003e-01
-1.439e+00 4.270e-01
-9.464e-02 1.784e-01
-1.532e+00 4.331e-01
-***************************************
-
-==> R(𝛌＝6439)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.287e+00 2.795e-01
--1.580e+00 2.127e-01
--2.879e+00 2.577e-01
--3.569e+00 2.889e-01
--3.661e+00 2.296e-01
--4.471e+00 5.134e-01
--3.998e+00 1.564e-01
--5.001e+00 4.661e-01
--6.583e-01 7.106e-01
--6.900e+00 1.424e+00
--2.060e+00 9.613e-01
--2.727e+00 1.047e+00
--2.954e+00 1.131e+00
--3.654e+00 1.178e+00
--4.406e+00 1.288e+00
--5.203e+00 1.386e+00
-1.448e-01 2.645e-01
--1.004e-01 2.434e-01
-9.754e-01 3.634e-01
-1.140e+00 3.905e-01
-9.284e-01 3.632e-01
-1.116e+00 3.672e-01
-2.092e+00 4.941e-01
-8.634e-01 3.384e-01
-***************************************
-
-==> Swiftb(𝛌＝4392)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--5.521e-01 3.134e-01
--1.135e+00 2.562e-01
--2.691e+00 4.640e-01
--3.626e+00 1.546e-01
--4.510e+00 2.215e-01
--1.036e+01 5.116e-01
--1.076e+01 1.445e-01
--1.213e+01 8.139e+04
--2.515e-01 3.686e-01
--6.156e-01 4.904e-01
--1.924e+00 8.403e-01
--2.417e+00 8.016e-01
--3.053e+00 1.013e+00
--6.570e+00 1.397e+00
--3.450e+00 1.038e+00
--6.333e+00 1.336e+00
-7.408e-02 1.220e-01
-3.089e-01 2.006e-01
-5.057e-01 2.378e-01
-4.831e-01 1.961e-01
-5.779e-01 2.557e-01
-6.277e-01 2.966e-01
-4.823e-01 2.573e-01
-6.110e-01 2.531e-01
-***************************************
-
-==> Swiftu(𝛌＝3465)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--3.178e-01 3.178e-01
--7.537e-01 2.541e-01
--2.402e+00 4.048e-01
--2.624e+00 1.540e-01
--3.413e+00 1.944e-01
--3.834e+00 2.633e-01
--1.236e+01 1.174e-01
--1.071e+01 3.318e+03
--1.376e-01 3.796e-01
--4.256e-01 4.494e-01
--1.794e+00 6.810e-01
--2.123e+00 6.992e-01
--2.817e+00 9.699e-01
--3.029e+00 9.065e-01
--4.139e+00 1.273e+00
--7.000e+00 1.331e+00
-7.840e-02 1.316e-01
-3.388e-01 1.716e-01
-6.616e-01 2.532e-01
-1.349e-01 1.501e-01
-9.909e-01 3.749e-01
-8.538e-01 3.311e-01
-1.397e+00 4.674e-01
-6.412e-01 2.404e-01
-***************************************
-
-==> Swiftv(𝛌＝5468)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--8.802e-01 3.314e-01
--1.441e+00 2.854e-01
--2.883e+00 4.166e-01
--4.483e+00 1.547e-01
--1.045e+01 2.071e-01
--1.041e+01 3.271e-01
--1.042e+01 1.274e-01
--1.078e+01 5.879e+03
--4.161e-01 9.142e-01
--7.660e-01 1.130e+00
--2.020e+00 1.179e+00
--3.201e+00 1.143e+00
--2.699e+00 1.335e+00
--2.859e+00 1.254e+00
--3.653e+00 1.312e+00
--5.501e+00 1.434e+00
-1.900e-02 2.266e-01
-3.560e-01 3.285e-01
-1.007e+00 3.907e-01
--1.270e+00 3.724e-01
-1.858e+00 4.560e-01
-9.302e-01 4.464e-01
--1.502e+00 4.728e-01
-2.653e+00 4.849e-01
-***************************************
-
-==> tables <==
-
-==> u(𝛌＝3471)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--4.289e-01 2.733e-01
--9.017e-01 2.118e-01
--2.581e+00 3.058e-01
--2.751e+00 2.567e-01
--2.764e+00 1.891e-01
--3.105e+00 2.139e-01
--3.298e+00 1.431e-01
--3.870e+00 2.222e-01
--1.925e-01 3.059e-01
--5.093e-01 3.513e-01
--1.923e+00 6.007e-01
--2.204e+00 6.982e-01
--2.568e+00 8.240e-01
--2.768e+00 8.080e-01
--3.425e+00 9.196e-01
--5.982e+00 1.396e+00
-1.646e-01 1.104e-01
-4.376e-01 1.640e-01
-7.202e-01 1.945e-01
-9.556e-01 2.502e-01
-4.856e-01 1.911e-01
-2.985e-01 1.448e-01
--6.801e-01 2.362e-01
--2.416e+00 4.887e-01
-***************************************
-
-==> UVM2(𝛌＝2246)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.003e-01 3.133e-01
--5.235e-01 2.446e-01
--1.863e+00 3.017e-01
--2.159e+00 1.943e-01
--2.718e+00 2.380e-01
--3.504e+00 4.807e-01
--1.139e+01 1.494e+04
--4.119e+00 8.683e-01
--2.431e-02 1.310e+00
--3.121e-01 1.446e+00
--1.526e+00 1.472e+00
--1.877e+00 1.484e+00
--2.455e+00 1.493e+00
--2.945e+00 1.508e+00
--3.356e+00 1.510e+00
--3.602e+00 1.508e+00
--1.862e-04 4.058e-01
-1.380e-01 4.538e-01
-2.475e-01 4.630e-01
-3.102e-02 4.416e-01
-6.040e-01 4.608e-01
--1.488e-01 4.709e-01
-1.171e+00 4.637e-01
-1.617e+00 4.645e-01
-***************************************
-
-==> UVW1(𝛌＝2600)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--2.417e-01 3.190e-01
--7.027e-01 2.509e-01
--1.998e+00 3.947e-01
--2.387e+00 1.542e-01
--2.912e+00 1.882e-01
--4.026e+00 2.472e-01
--1.069e+01 1.123e-01
--9.716e+00 2.578e+02
--9.867e-02 6.982e-01
--4.150e-01 8.299e-01
--1.587e+00 8.040e-01
--1.963e+00 9.829e-01
--2.539e+00 1.024e+00
--3.274e+00 1.148e+00
--3.460e+00 1.229e+00
--6.742e+00 1.428e+00
-6.734e-02 2.545e-01
-1.353e-01 2.890e-01
-3.585e-01 2.644e-01
-1.587e-01 2.261e-01
-2.959e-01 2.424e-01
-1.111e-01 2.882e-01
-2.364e-01 2.837e-01
-7.463e-01 3.484e-01
-***************************************
-
-==> UVW2(𝛌＝1928)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.648e-01 2.888e-01
--4.974e-01 2.247e-01
--1.686e+00 2.872e-01
--1.897e+00 1.430e-01
--2.600e+00 1.717e-01
--2.801e+00 2.023e-01
--3.814e+00 1.056e-01
--9.888e+00 3.109e+02
--6.101e-02 2.539e-01
--3.020e-01 3.388e-01
--1.433e+00 4.484e-01
--1.756e+00 5.656e-01
--2.485e+00 7.064e-01
--2.799e+00 8.287e-01
--3.317e+00 9.293e-01
--7.000e+00 1.393e+00
-4.622e-02 9.341e-02
-1.473e-01 1.171e-01
-2.243e-01 1.191e-01
-1.606e-01 1.081e-01
-2.365e-01 1.249e-01
-1.484e-01 1.182e-01
-8.110e-01 2.388e-01
-9.689e-01 2.720e-01
-***************************************
-
-==> V(𝛌＝5404)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.044e+00 2.743e-01
--1.457e+00 2.109e-01
--2.848e+00 2.322e-01
--2.957e+00 1.713e-01
--4.232e+00 2.257e-01
--4.052e+00 1.742e-01
--4.515e+00 1.593e-01
--5.250e+00 3.713e-01
--5.026e-01 1.143e+00
--7.878e-01 1.224e+00
--2.118e+00 1.280e+00
--2.276e+00 1.263e+00
--3.421e+00 1.357e+00
--3.378e+00 1.364e+00
--4.224e+00 1.433e+00
--4.930e+00 1.494e+00
-9.142e-02 3.750e-01
-4.814e-01 3.731e-01
-8.550e-01 3.884e-01
-3.929e-01 3.841e-01
-3.975e-01 3.691e-01
-5.311e-01 3.951e-01
--1.203e-01 3.770e-01
--2.905e+00 5.243e-01
-***************************************
-
-==> z(𝛌＝9157)≻HST(𝛌＝1367)_Δt＝0.1 <==
-# 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603 
-6.533e-02 2.816e-01
--9.694e-02 2.103e-01
--1.175e+00 2.185e-01
--1.525e+00 1.729e-01
--2.166e+00 1.536e-01
--2.492e+00 1.538e-01
--3.258e+00 1.803e-01
--9.328e+00 7.607e+01
--1.852e+00 2.783e-01
--2.156e+00 2.086e-01
--3.798e+00 2.699e-01
--4.146e+00 2.393e-01
--4.577e+00 2.541e-01
--4.338e+00 1.832e-01
--4.672e+00 1.459e-01
--5.230e+00 2.311e-01
--9.339e-01 1.264e+00
--1.131e+00 1.351e+00
--2.536e+00 1.393e+00
--3.031e+00 1.413e+00
--4.060e+00 1.498e+00
--3.801e+00 1.463e+00
--4.139e+00 1.444e+00
--5.151e+00 1.548e+00
-2.445e-01 3.942e-01
-8.131e-01 4.160e-01
-1.621e+00 4.545e-01
-1.460e+00 4.611e-01
-2.909e-01 4.349e-01
-6.991e-01 4.403e-01
--5.831e-01 4.255e-01
--2.921e-01 4.034e-01
-***************************************
diff --git a/pres/log b/pres/log
index a95a218..b012e58 100644
--- a/pres/log
+++ b/pres/log
@@ -11,18 +11,24 @@ Given calibration data in the form of light curves produced by a known transfer
 
 Did things, forgot to keep track
 
-
+Week 3-4ish:
 A script has been created to read an output file from psdlag and then plot the source and reprocessed power spectral densities, the cross-correlation power spectral density, and the phase difference between the signals.
 
 scripts/power_lag_plot.pl
 
 
+2016-
+
 
 2016-07-11
 
 	created gnuplot templates for plotting the timelags and power spectrum densities
 		- still need to modify this for points to be evenly spaced in log space
 
+2016-07-14
+	Ran through results with simple error checking and determined they made sense.
+	Project can continue.
+
 2016-07-15
 	Determined frequency binning for PSD and timelags:
 
@@ -34,3 +40,8 @@ scripts/power_lag_plot.pl
 		0.16658029	0.25819945      
 		0.25819945	0.40020915      
 		0.40020915	0.62032418      
+
+
+
+2016-07-18
+	Ran first MCMC run, need to analyze results
diff --git a/scripts/analyze_lightcurves.sh b/scripts/analyze_lightcurves.sh
index 5aa6c30..331c020 100755
--- a/scripts/analyze_lightcurves.sh
+++ b/scripts/analyze_lightcurves.sh
@@ -11,20 +11,85 @@ error_type="2"
 echo Using lightcurves in $lc_dir.
 
 for echo_curve in $lc_dir/*
-do 
+do
+    # Determine band and inputs for band
     echo_band=$(basename $echo_curve|sed 's@\(.*\)\.lc@\1@')
     echo -n $(date "+%R")\: Running psdlag for $echo_band against $ref_band.
     newfilename="analyses/${echo_band}≻${ref_band}_Δt＝${timestep}"
+    case $echo_band in
+        "g(𝛌＝4775)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -9.745e-01 -1.384e+00 -2.748e+00 -3.305e+00 -3.314e+00 -3.389e+00 -4.198e+00 -4.465e+00 -4.700e-01 -7.487e-01 -2.046e+00 -2.428e+00 -2.953e+00 -3.086e+00 -3.761e+00 -4.290e+00 9.862e-02 3.899e-01 8.650e-01 5.516e-01 2.228e-01 9.508e-01 -2.872e-01 9.059e-02"
+        ;;
+        "HST(𝛌＝1158)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 2.402e-01 1.329e-01 -8.990e-01 -1.229e+00 -1.825e+00 -2.236e+00 -2.758e+00 -4.749e+00 1.518e-01 1.794e-02 -1.040e+00 -1.377e+00 -1.989e+00 -2.375e+00 -2.886e+00 -3.737e+00 -3.034e-02 -6.855e-02 -6.829e-02 -4.907e-02 -5.681e-02 -1.644e-01 -3.273e-01 1.333e-01"
+        ;;
+        "HST(𝛌＝1367)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 6.534e-02 -9.689e-02 -1.175e+00 -1.525e+00 -2.167e+00 -2.492e+00 -3.258e+00 -1.115e+01 6.378e-02 -9.676e-02 -1.181e+00 -1.520e+00 -2.180e+00 -2.532e+00 -3.129e+00 -3.794e+00 -8.535e-05 -5.126e-04 -4.066e-04 1.080e-03 -2.682e-03 -5.462e-03 -2.209e-02 1.564e-01"
+        ;;
+        "HST(𝛌＝1479)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 1.224e-01 -7.607e-02 -1.174e+00 -1.500e+00 -2.075e+00 -2.554e+00 -3.337e+00 -1.152e+01 9.340e-02 -8.842e-02 -1.173e+00 -1.512e+00 -2.113e+00 -2.547e+00 -3.178e+00 -7.000e+00 -6.035e-02 2.044e-02 1.258e-01 5.764e-02 2.580e-02 -4.737e-02 9.579e-02 6.411e-01"
+        ;;
+        "HST(𝛌＝1746)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 3.934e-02 -1.638e-01 -1.339e+00 -1.655e+00 -2.232e+00 -2.793e+00 -4.025e+00 -1.072e+01 5.043e-02 -1.304e-01 -1.254e+00 -1.587e+00 -2.207e+00 -2.693e+00 -3.209e+00 -3.820e+00 3.040e-02 2.438e-02 1.450e-01 1.458e-01 1.772e-01 2.076e-01 -1.390e-01 3.929e-01"
+        ;;
+        "i(𝛌＝7647)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.510e+00 -1.872e+00 -3.441e+00 -3.722e+00 -4.794e+00 -4.336e+00 -4.644e+00 -5.061e+00 -7.402e-01 -9.967e-01 -2.393e+00 -2.723e+00 -3.701e+00 -3.742e+00 -4.317e+00 -4.979e+00 2.638e-01 7.095e-01 1.297e+00 1.035e+00 9.267e-01 8.982e-01 7.318e-01 9.092e-01"
+        ;;
+        "I(𝛌＝8560)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.497e+00 -1.786e+00 -3.555e+00 -4.203e+00 -3.826e+00 -4.096e+00 -4.304e+00 -1.109e+01 -7.542e-01 -9.532e-01 -2.397e+00 -3.214e+00 -3.008e+00 -3.489e+00 -4.353e+00 -4.061e+00 1.675e-01 7.660e-01 1.781e+00 1.971e+00 6.572e-01 8.520e-01 2.405e+00 -1.573e-01"
+        ;;
+        "r(𝛌＝6175)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.372e+00 -1.786e+00 -3.214e+00 -3.558e+00 -4.416e+00 -4.366e+00 -4.559e+00 -4.586e+00 -6.693e-01 -5.998e+00 -2.390e+00 -2.927e+00 -3.411e+00 -3.731e+00 -3.930e+00 -7.000e+00 2.905e-01 1.835e+00 8.760e-01 1.476e+00 -3.280e-01 1.439e+00 9.464e-02 1.532e+00"
+        ;;
+        "R(𝛌＝6439)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.287e+00 -1.580e+00 -2.879e+00 -3.569e+00 -3.661e+00 -4.471e+00 -3.998e+00 -5.001e+00 -6.583e-01 -6.900e+00 -2.060e+00 -2.727e+00 -2.954e+00 -3.654e+00 -4.406e+00 -5.203e+00 1.448e-01 -1.004e-01 9.754e-01 1.140e+00 9.284e-01 1.116e+00 2.092e+00 8.634e-01"
+        ;;
+        "Swiftb(𝛌＝4392)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -5.521e-01 -1.135e+00 -2.691e+00 -3.626e+00 -4.510e+00 -1.036e+01 -1.076e+01 -1.213e+01 -2.515e-01 -6.156e-01 -1.924e+00 -2.417e+00 -3.053e+00 -6.570e+00 -3.450e+00 -6.333e+00 7.408e-02 3.089e-01 5.057e-01 4.831e-01 5.779e-01 6.277e-01 4.823e-01 6.110e-01"
+        ;;
+        "Swiftu(𝛌＝3465)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -3.178e-01 -7.537e-01 -2.402e+00 -2.624e+00 -3.413e+00 -3.834e+00 -1.236e+01 -1.071e+01 -1.376e-01 -4.256e-01 -1.794e+00 -2.123e+00 -2.817e+00 -3.029e+00 -4.139e+00 -7.000e+00 7.840e-02 3.388e-01 6.616e-01 1.349e-01 9.909e-01 8.538e-01 1.397e+00 6.412e-01"
+        ;;
+        "Swiftv(𝛌＝5468)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -8.802e-01 -1.441e+00 -2.883e+00 -4.483e+00 -1.045e+01 -1.041e+01 -1.042e+01 -1.078e+01 -4.161e-01 -7.660e-01 -2.020e+00 -3.201e+00 -2.699e+00 -2.859e+00 -3.653e+00 -5.501e+00 1.900e-02 3.560e-01 1.007e+00 -1.270e+00 1.858e+00 9.302e-01 -1.502e+00 2.653e+00"
+        ;;
+        "u(𝛌＝3471)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -4.289e-01 -9.017e-01 -2.581e+00 -2.751e+00 -2.764e+00 -3.105e+00 -3.298e+00 -3.870e+00 -1.925e-01 -5.093e-01 -1.923e+00 -2.204e+00 -2.568e+00 -2.768e+00 -3.425e+00 -5.982e+00 1.646e-01 4.376e-01 7.202e-01 9.556e-01 4.856e-01 2.985e-01 -6.801e-01 -2.416e+00"
+        ;;
+        "UVM2(𝛌＝2246)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.003e-01 -5.235e-01 -1.863e+00 -2.159e+00 -2.718e+00 -3.504e+00 -1.139e+01 -4.119e+00 -2.431e-02 -3.121e-01 -1.526e+00 -1.877e+00 -2.455e+00 -2.945e+00 -3.356e+00 -3.602e+00 -1.862e-04 1.380e-01 2.475e-01 3.102e-02 6.040e-01 -1.488e-01 1.171e+00 1.617e+00"
+        ;;
+        "UVW1(𝛌＝2600)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -2.417e-01 -7.027e-01 -1.998e+00 -2.387e+00 -2.912e+00 -4.026e+00 -1.069e+01 -9.716e+00 -9.867e-02 -4.150e-01 -1.587e+00 -1.963e+00 -2.539e+00 -3.274e+00 -3.460e+00 -6.742e+00 6.734e-02 1.353e-01 3.585e-01 1.587e-01 2.959e-01 1.111e-01 2.364e-01 7.463e-01"
+        ;;
+        "UVW2(𝛌＝1928)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.648e-01 -4.974e-01 -1.686e+00 -1.897e+00 -2.600e+00 -2.801e+00 -3.814e+00 -9.888e+00 -6.101e-02 -3.020e-01 -1.433e+00 -1.756e+00 -2.485e+00 -2.799e+00 -3.317e+00 -7.000e+00 4.622e-02 1.473e-01 2.243e-01 1.606e-01 2.365e-01 1.484e-01 8.110e-01 9.689e-01"
+        ;;
+        "V(𝛌＝5404)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.044e+00 -1.457e+00 -2.848e+00 -2.957e+00 -4.232e+00 -4.052e+00 -4.515e+00 -5.250e+00 -5.026e-01 -7.878e-01 -2.118e+00 -2.276e+00 -3.421e+00 -3.378e+00 -4.224e+00 -4.930e+00 9.142e-02 4.814e-01 8.550e-01 3.929e-01 3.975e-01 5.311e-01 -1.203e-01 -2.905e+00"
+        ;;
+        "z(𝛌＝9157)")
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -1.852e+00 -2.156e+00 -3.798e+00 -4.146e+00 -4.577e+00 -4.338e+00 -4.672e+00 -5.230e+00 -9.339e-01 -1.131e+00 -2.536e+00 -3.031e+00 -4.060e+00 -3.801e+00 -4.139e+00 -5.151e+00 2.445e-01 8.131e-01 1.621e+00 1.460e+00 2.909e-01 6.991e-01 -5.831e-01 -2.921e-01"
+        ;;
+        *)
+            echo "Unknown band: using default parameters."
+            initial_params="6.533e-02 -9.694e-02 -1.175e+00 -1.525e+00 -2.166e+00 -2.492e+00 -3.258e+00 -9.328e+00 -9.745e-01 -1.384e+00 -2.748e+00 -3.305e+00 -3.314e+00 -3.389e+00 -4.198e+00 -4.465e+00 -4.700e-01 -7.487e-01 -2.046e+00 -2.428e+00 -2.953e+00 -3.086e+00 -3.761e+00 -4.290e+00 9.862e-02 3.899e-01 8.650e-01 5.516e-01 2.228e-01 9.508e-01 -2.872e-01 9.059e-02"
+        ;;
+    esac
+
+    # Write temporary input file
     echo "2" > tmp.psdlagargs
     echo "$ref_curve 0" >> tmp.psdlagargs
     echo "$echo_curve 0" >> tmp.psdlagargs
     echo "8 0.005 0.019 0.0425 0.07 0.11 0.17 0.24 0.4 0.603" >> tmp.psdlagargs
     echo "0" >> tmp.psdlagargs
-    echo ".9999 .9999 .9999 0.49 0.42 0.30 0.36 1.2 -0.60 -0.89 -0.72 2.39 -0.74 -1.5 -1 -3.1 -1.3 -2.3 -2 -1.9 -1.6 -2.6 -2 -1.3 -1.8 -2.6 -1.2 0 -2 -2.7 -1.3 0.63" >> tmp.psdlagargs
+    echo "$initial_params" >> tmp.psdlagargs
     echo "0:0 0" >> tmp.psdlagargs
     echo $error_type >> tmp.psdlagargs
     echo "0" >> tmp.psdlagargs
     echo "100 50 50 mcmc.dat" >> tmp.psdlagargs
+
+    # Run psdlag with inputs
     time bin/psdlag tmp.psdlagargs >> $newfilename
     echo ""
-done
\ No newline at end of file
+done