Brought in cantilever-data generation scripts.
authorW. Trevor King <wking@drexel.edu>
Mon, 15 Mar 2010 11:18:13 +0000 (07:18 -0400)
committerW. Trevor King <wking@drexel.edu>
Mon, 15 Mar 2010 11:18:13 +0000 (07:18 -0400)
46 files changed:
tex/src/cantilever/simulations.tex
tex/src/figures/Makefile
tex/src/figures/cantilever-data/Makefile
tex/src/figures/cantilever-data/README [new file with mode: 0644]
tex/src/figures/cantilever-data/avg_data.py [new file with mode: 0755]
tex/src/figures/cantilever-data/data [new file with mode: 0644]
tex/src/figures/cantilever-data/extract_f_v_k_data.sh [new file with mode: 0755]
tex/src/figures/cantilever-data/fit_data.py [new file with mode: 0755]
tex/src/figures/cantilever-data/get_loading_rates.py [new file with mode: 0755]
tex/src/figures/cantilever-data/make_plots.sh [new file with mode: 0755]
tex/src/figures/cantilever-data/makeplotsplits.sh [new file with mode: 0755]
tex/src/figures/cantilever-data/plot_avg.gp [new file with mode: 0644]
tex/src/figures/cantilever-data/plot_raw.gp [new file with mode: 0644]
tex/src/figures/cantilever-data/plot_splits.pdf [deleted file]
tex/src/figures/cantilever-data/plot_splits_nolandscape.pdf [deleted file]
tex/src/figures/cantilever-data/plot_splits_rotated.pdf [deleted file]
tex/src/figures/cantilever-sim/Makefile [new file with mode: 0644]
tex/src/figures/cantilever-sim/i-dep.asy [moved from tex/src/figures/cantilever-data/i-dep.asy with 100% similarity]
tex/src/figures/cantilever-sim/i-dep.d/i_dep_0.1_1_1e-6 [moved from tex/src/figures/cantilever-data/i-dep.d/i_dep_0.1_1_1e-6 with 100% similarity]
tex/src/figures/cantilever-sim/i-dep.d/i_dep_0.1_30_1e-6 [moved from tex/src/figures/cantilever-data/i-dep.d/i_dep_0.1_30_1e-6 with 100% similarity]
tex/src/figures/cantilever-sim/i-dep.d/i_dep_127_30_1e-6 [moved from tex/src/figures/cantilever-data/i-dep.d/i_dep_127_30_1e-6 with 100% similarity]
tex/src/figures/cantilever-sim/i-dep.d/i_dep_127_8_1e-6 [moved from tex/src/figures/cantilever-data/i-dep.d/i_dep_127_8_1e-6 with 100% similarity]
tex/src/figures/cantilever-sim/i-dep.d/i_dep_27_30_1e-6 [moved from tex/src/figures/cantilever-data/i-dep.d/i_dep_27_30_1e-6 with 100% similarity]
tex/src/figures/cantilever-sim/i-dep.d/i_dep_27_8_1e-6 [moved from tex/src/figures/cantilever-data/i-dep.d/i_dep_27_8_1e-6 with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.asy [moved from tex/src/figures/cantilever-data/sim-loading-rate.asy with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.d/loading_rate_dep_0.1_1 [moved from tex/src/figures/cantilever-data/sim-loading-rate.d/loading_rate_dep_0.1_1 with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.d/loading_rate_dep_0.1_30 [moved from tex/src/figures/cantilever-data/sim-loading-rate.d/loading_rate_dep_0.1_30 with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.d/loading_rate_dep_127_30 [moved from tex/src/figures/cantilever-data/sim-loading-rate.d/loading_rate_dep_127_30 with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.d/loading_rate_dep_127_8 [moved from tex/src/figures/cantilever-data/sim-loading-rate.d/loading_rate_dep_127_8 with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.d/loading_rate_dep_27_30 [moved from tex/src/figures/cantilever-data/sim-loading-rate.d/loading_rate_dep_27_30 with 100% similarity]
tex/src/figures/cantilever-sim/loading-rate.d/loading_rate_dep_27_8 [moved from tex/src/figures/cantilever-data/sim-loading-rate.d/loading_rate_dep_27_8 with 100% similarity]
tex/src/figures/cantilever-sim/stiffness.d/stiffness_data-127_8 [moved from tex/src/figures/cantilever-data/stiffness.d/stiffness_data-127_8 with 100% similarity]
tex/src/figures/cantilever-sim/stiffness.d/stiffness_data-27_8 [moved from tex/src/figures/cantilever-data/stiffness.d/stiffness_data-27_8 with 100% similarity]
tex/src/figures/cantilever-sim/stiffness.d/stiffness_force-127_8 [moved from tex/src/figures/cantilever-data/stiffness.d/stiffness_force-127_8 with 100% similarity]
tex/src/figures/cantilever-sim/stiffness.d/stiffness_force-27_8 [moved from tex/src/figures/cantilever-data/stiffness.d/stiffness_force-27_8 with 100% similarity]
tex/src/figures/cantilever-sim/stiffness.d/stiffness_stiffness-127_8 [moved from tex/src/figures/cantilever-data/stiffness.d/stiffness_stiffness-127_8 with 100% similarity]
tex/src/figures/cantilever-sim/stiffness.d/stiffness_stiffness-27_8 [moved from tex/src/figures/cantilever-data/stiffness.d/stiffness_stiffness-27_8 with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.asy [moved from tex/src/figures/cantilever-data/v-dep.asy with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.d/v_dep_0.1_1 [moved from tex/src/figures/cantilever-data/v-dep.d/v_dep_0.1_1 with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.d/v_dep_0.1_30 [moved from tex/src/figures/cantilever-data/v-dep.d/v_dep_0.1_30 with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.d/v_dep_127_30 [moved from tex/src/figures/cantilever-data/v-dep.d/v_dep_127_30 with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.d/v_dep_127_8 [moved from tex/src/figures/cantilever-data/v-dep.d/v_dep_127_8 with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.d/v_dep_27_30 [moved from tex/src/figures/cantilever-data/v-dep.d/v_dep_27_30 with 100% similarity]
tex/src/figures/cantilever-sim/v-dep.d/v_dep_27_8 [moved from tex/src/figures/cantilever-data/v-dep.d/v_dep_27_8 with 100% similarity]
tex/src/figures/cantilever-sim/wtk_graph.asy [moved from tex/src/figures/cantilever-data/wtk_graph.asy with 100% similarity]
tex/src/root.bib

index 7b14616dbce6b9530341e5c4577489ff2994cddc..fd9ff150e377e1eaef2a622e4028732f0bb77140 100644 (file)
@@ -9,19 +9,19 @@ following results
 (\cref{fig:cant:sim:v-dep,fig:cant:sim:load-dep,fig:cant:sim:i-dep}).
 
 \begin{figure}
-  \asyfig{figures/cantilever-data/v-dep}
+  \asyfig{figures/cantilever-sim/v-dep}
   \caption{Unfolding force velocity dependence for different
     cantilevers.\label{fig:cant:sim:v-dep}}
 \end{figure}
 
 \begin{figure}
-  \asyfig{figures/cantilever-data/loading-rate-sim}
+  \asyfig{figures/cantilever-sim/loading-rate}
   \caption{Unfolding force loading rate dependence simulations for
     different cantilevers.\label{fig:cant:sim:load-dep}}
 \end{figure}
 
 \begin{figure}
-  \asyfig{figures/cantilever-data/i-dep}
+  \asyfig{figures/cantilever-sim/i-dep}
   \caption{Unfolding force peak index dependence simulations for
     different cantilevers.\label{fig:cant:sim:i-dep}}
 \end{figure}
index 742e0bff645a513078872f6574717ee80009a4ca..96fdf01863b52b0771d75584c0181497440876fa 100644 (file)
@@ -1,5 +1,6 @@
-SUBDIRS = biotin-streptavidin cantilever-data contour expt-sawtooth fit-space i27 \
-       kappa-sawteeth order-dep schematic sim-hist sim-sawtooth v-dep \
+SUBDIRS = biotin-streptavidin cantilever-data cantilever-sim contour \
+       expt-sawtooth fit-space i27 kappa-sawteeth order-dep schematic \
+       sim-hist sim-sawtooth v-dep \
        # calibration-concept-map
 
 all :
index 6f968f86e0c140b6c27b241c2c117e02f2ac4da8..02e8e7678596eda2ee2a34ac048be71f1f279d05 100644 (file)
@@ -3,7 +3,7 @@ ASYPROCESS = ../asy/asyprocess \
   --pretex='\documentclass{drexel-thesis} \input{packages}' \
   --
 
-FIGS = loading-rate sim-loading-rate v-dep i-dep 
+FIGS = plot-splits loading-rate 
 
 all : $(FIGS:%=%_.tex)
 
diff --git a/tex/src/figures/cantilever-data/README b/tex/src/figures/cantilever-data/README
new file mode 100644 (file)
index 0000000..5c49101
--- /dev/null
@@ -0,0 +1,5 @@
+./extract_f_v_k_data.sh > data
+./avg_data.py
+./make_plots.sh
+./get_loading_rates
+./makeplotsplits.sh
diff --git a/tex/src/figures/cantilever-data/avg_data.py b/tex/src/figures/cantilever-data/avg_data.py
new file mode 100755 (executable)
index 0000000..b1bb276
--- /dev/null
@@ -0,0 +1,66 @@
+#!/usr/bin/python
+
+import numpy, sys
+
+DATAFILE='./data'
+KFILE='./spring_constants'
+AVGFILE='./averaged_data'
+
+KCUTS=[45, 90]
+
+FofV = {}
+Ks = []
+for i in range(len(KCUTS)+1):
+    Ks.append([])
+
+for line in file(DATAFILE, 'r'):
+    if line[0] == '#':
+        continue # ignore comments
+    fields = line.strip().split('\t')
+    force = float(fields[0])
+    rate = float(fields[1])
+    spring = float(fields[2])
+    
+    ispring = 0
+    while ispring < len(KCUTS) and spring > KCUTS[ispring]:
+        ispring += 1
+    Ks[ispring].append(spring)
+    if rate not in FofV.keys():
+        FofV[rate]=[]
+        for i in range(len(Ks)):
+            FofV[rate].append([])
+    FofV[rate][ispring].append(force)
+
+avgK = [0]*len(Ks)
+stdK = [0]*len(Ks)
+numK = [0]*len(Ks)
+kf = file(KFILE, 'w')
+for i in range(len(Ks)):
+    K = numpy.array(Ks[i])
+    avgK[i] = K.mean()
+    stdK[i] = K.std()
+    numK[i] = len(K)
+    
+    if i == 1:
+        continue # poor calibration bumps for the older cantilevers
+    print >> kf, "K(%d) = %g +/- %g, (# %d)" % (i, avgK[i], stdK[i], numK[i])
+kf.close()
+
+
+Vs = FofV.keys()
+Vs.sort()
+af = file(AVGFILE, 'w')
+print >> af, '#'+'\t'.join(['Pulling speed (nm/s)','Spring constant (pN/nm)',
+                            'Mean force (pN)','Std. force (pN)','Events (#)'])
+for V in Vs:
+    for i,k in enumerate(avgK):
+        if i == 1:
+            continue # poor calibration bumps for the older cantilevers
+        F = numpy.array(FofV[V][i])
+        if len(F) == 0:
+            continue
+        outs= [V, k, F.mean(), F.std(), len(F)]
+        souts = ["%.2f" % (x) for x in outs]
+        souts[-1] = "%d" % outs[-1] # special treatment for the integer
+        print >> af, '\t'.join(souts)
+af.close()
diff --git a/tex/src/figures/cantilever-data/data b/tex/src/figures/cantilever-data/data
new file mode 100644 (file)
index 0000000..6bc2313
--- /dev/null
@@ -0,0 +1,433 @@
+#Unfolding force (pN)  Unfolding rate (nm/s)   Spring constant (pN/nm) Datafile
+231.2  1357.2088083    14.1236 34.7852617567   /home/wking/rsrch/data/scaled_unfold/20090114/20090114122724_unfold
+232.8  1357.2088083    14.1236 28.9221197049   /home/wking/rsrch/data/scaled_unfold/20090114/20090114122724_unfold
+249.9  1357.2088083    14.1236 31.2565188551   /home/wking/rsrch/data/scaled_unfold/20090114/20090114122724_unfold
+271    3237.39401435   14.7206 139.466934138   /home/wking/rsrch/data/scaled_unfold/20090114/20090114125600_unfold
+264.6  3237.39401435   14.7206 109.650535266   /home/wking/rsrch/data/scaled_unfold/20090114/20090114125600_unfold
+246.3  2096.14400083   14.7206 144.445283097   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130514_unfold
+256.2  2096.14400083   14.7206 180.184538311   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130514_unfold
+215.3  3237.39401435   14.7206 361.616911403   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+246.7  3237.39401435   14.7206 222.570838487   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+234.3  3237.39401435   14.7206 153.549598101   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+250.6  3237.39401435   14.7206 168.991967549   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+268.8  3237.39401435   14.7206 162.258187999   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+279.3  3237.39401435   14.7206 153.290606579   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+186.1  3237.39401435   14.7206 57.1723782934   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130522_unfold
+229.3  878.76393444    14.7206 149.829250822   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130808_unfold
+193    878.76393444    14.7206 51.8646474106   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130808_unfold
+243.6  2096.14400083   14.7206 121.806927888   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130823_unfold
+193.9  2096.14400083   14.7206 67.6844897868   /home/wking/rsrch/data/scaled_unfold/20090114/20090114130823_unfold
+199    368.403149864   14.7206 12.5809675679   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131207_unfold
+253.4  368.403149864   14.7206 17.8712367999   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131207_unfold
+228.6  368.403149864   14.7206 10.4773855821   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131207_unfold
+283.4  368.403149864   14.7206 12.4557104969   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131207_unfold
+208.1  2096.14400083   14.7206 318.404273726   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131556_unfold
+255.8  2096.14400083   14.7206 190.476605355   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131556_unfold
+173    238.533230447   14.7206 6.00149607805   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131823_unfold
+171.5  238.533230447   14.7206 6.74571975704   /home/wking/rsrch/data/scaled_unfold/20090114/20090114131823_unfold
+267.4  878.76393444    14.7206 241.835834758   /home/wking/rsrch/data/scaled_unfold/20090114/20090114133826_unfold
+174.6  878.76393444    14.7206 60.7928889846   /home/wking/rsrch/data/scaled_unfold/20090114/20090114133826_unfold
+283.3  878.76393444    14.7206 101.057852461   /home/wking/rsrch/data/scaled_unfold/20090114/20090114133826_unfold
+221.9  2096.14400083   13.7398 225.125865689   /home/wking/rsrch/data/scaled_unfold/20090114/20090114142553_unfold
+200.2  2096.14400083   13.7398 128.032475571   /home/wking/rsrch/data/scaled_unfold/20090114/20090114142553_unfold
+129.2  238.533230447   13.7398 6.09452403792   /home/wking/rsrch/data/scaled_unfold/20090114/20090114142643_unfold
+151.9  238.533230447   13.7398 4.63947133219   /home/wking/rsrch/data/scaled_unfold/20090114/20090114142643_unfold
+100.2  238.533230447   13.7398 1.87367852516   /home/wking/rsrch/data/scaled_unfold/20090114/20090114142643_unfold
+221.2  3237.39401435   13.7398 541.616018601   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143217_unfold
+258.1  3237.39401435   13.7398 392.695893941   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143217_unfold
+209.4  3237.39401435   13.7398 123.862694989   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143217_unfold
+181.6  3237.39401435   13.7398 88.0571171903   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143217_unfold
+199.1  2096.14400083   13.7398 154.779273021   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143721_unfold
+165.4  2096.14400083   13.7398 67.8731427469   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143721_unfold
+214.2  2096.14400083   13.7398 123.483843089   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143721_unfold
+187    154.445210495   13.7398 6.10521917087   /home/wking/rsrch/data/scaled_unfold/20090114/20090114143800_unfold
+181.9  1357.2088083    13.7398 111.073968871   /home/wking/rsrch/data/scaled_unfold/20090114/20090114144025_unfold
+150.3  1357.2088083    13.7398 65.403892472    /home/wking/rsrch/data/scaled_unfold/20090114/20090114144025_unfold
+178.1  1357.2088083    13.7398 44.6793139692   /home/wking/rsrch/data/scaled_unfold/20090114/20090114144025_unfold
+195.4  1357.2088083    13.7398 29.0171243215   /home/wking/rsrch/data/scaled_unfold/20090114/20090114144025_unfold
+179    1357.2088083    13.7398 25.664818565    /home/wking/rsrch/data/scaled_unfold/20090114/20090114144025_unfold
+221.8  100.0   13.7398 9.425   /home/wking/rsrch/data/scaled_unfold/20090114/20090114144407_unfold
+315.1  5000.0  13.7398 2566.5  /home/wking/rsrch/data/scaled_unfold/20090114/20090114144713_unfold
+311.1  5000.0  13.7398 698.5   /home/wking/rsrch/data/scaled_unfold/20090114/20090114144713_unfold
+248.1  5000.0  13.7398 266.35  /home/wking/rsrch/data/scaled_unfold/20090114/20090114144713_unfold
+225.7  5000.0  13.7398 242.95  /home/wking/rsrch/data/scaled_unfold/20090114/20090114144713_unfold
+249.4  5000.0  13.7398 254.0   /home/wking/rsrch/data/scaled_unfold/20090114/20090114144713_unfold
+174.1  2096.14400083   13.7398 183.831828873   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145142_unfold
+186    2096.14400083   13.7398 84.7890248336   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145142_unfold
+237.7  2096.14400083   13.7398 130.296311092   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145142_unfold
+272.1  2096.14400083   13.7398 130.820347092   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145142_unfold
+223.3  2096.14400083   13.7398 62.6327827448   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145142_unfold
+240.8  2096.14400083   13.7398 65.0223869057   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145142_unfold
+215.4  2096.14400083   13.7398 239.170030495   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145808_unfold
+169.2  238.533230447   13.7398 8.40114037634   /home/wking/rsrch/data/scaled_unfold/20090114/20090114145858_unfold
+224.5  2096.14400083   13.7398 426.146075369   /home/wking/rsrch/data/scaled_unfold/20090114/20090114150120_unfold
+210.6  5000.0  13.7398 293.5   /home/wking/rsrch/data/scaled_unfold/20090114/20090114150312_unfold
+242.9  5000.0  13.7398 235.2   /home/wking/rsrch/data/scaled_unfold/20090114/20090114150312_unfold
+219.9  878.76393444    13.7398 58.8244577714   /home/wking/rsrch/data/scaled_unfold/20090114/20090114150730_unfold
+208.4  878.76393444    13.7398 27.3031954431   /home/wking/rsrch/data/scaled_unfold/20090114/20090114150730_unfold
+219.2  238.533230447   13.7398 11.876569544    /home/wking/rsrch/data/scaled_unfold/20090114/20090114150942_unfold
+202.1  238.533230447   13.7398 7.71177934035   /home/wking/rsrch/data/scaled_unfold/20090114/20090114150942_unfold
+197.5  2096.14400083   13.7398 139.833766295   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151028_unfold
+217.1  2096.14400083   13.7398 155.911190782   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151028_unfold
+199.4  2096.14400083   13.7398 87.8493950748   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151028_unfold
+227.8  2096.14400083   13.7398 52.6341758608   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151028_unfold
+264.2  2096.14400083   13.7398 56.8893481825   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151028_unfold
+228.2  5000.0  13.7398 691.5   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151218_unfold
+177.9  154.445210495   13.7398 10.8420537767   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+194.6  154.445210495   13.7398 7.97246176575   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+197.2  154.445210495   13.7398 4.9607801611    /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+222.7  154.445210495   13.7398 5.5847388115    /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+203.5  154.445210495   13.7398 3.55069538928   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+184.4  154.445210495   13.7398 6.5839993234    /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+205.3  154.445210495   13.7398 2.92828119099   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+228    154.445210495   13.7398 2.79545830996   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151241_unfold
+194.1  3237.39401435   13.7398 109.618161326   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151523_unfold
+216.4  3237.39401435   13.7398 128.58929025    /home/wking/rsrch/data/scaled_unfold/20090114/20090114151523_unfold
+264.8  3237.39401435   13.7398 170.966777898   /home/wking/rsrch/data/scaled_unfold/20090114/20090114151523_unfold
+178.9  5000.0  13.7398 296.55  /home/wking/rsrch/data/scaled_unfold/20090114/20090114151530_unfold
+234.7  5000.0  13.7398 295.55  /home/wking/rsrch/data/scaled_unfold/20090114/20090114151530_unfold
+206.7  1357.2088083    13.7398 79.4102873736   /home/wking/rsrch/data/scaled_unfold/20090114/20090114152437_unfold
+204.8  1357.2088083    13.7398 86.4270569125   /home/wking/rsrch/data/scaled_unfold/20090114/20090114152437_unfold
+176.1  238.533230447   28.8414 12.7162065151   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031036_unfold
+108.5  238.533230447   28.8414 3.90478898242   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031036_unfold
+213.2  238.533230447   28.8414 8.95930813559   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031036_unfold
+162.1  238.533230447   28.8414 4.61800334145   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031036_unfold
+233.3  238.533230447   28.8414 8.41783770247   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031036_unfold
+206.1  568.981020276   28.8414 33.1431444311   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031057_unfold
+220    568.981020276   28.8414 23.8573741802   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031057_unfold
+181.6  568.981020276   28.8414 13.974173858    /home/wking/rsrch/data/scaled_unfold/20090203/20090203031057_unfold
+206.2  568.981020276   28.8414 13.4051928377   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031057_unfold
+206.6  568.981020276   28.8414 17.8147957448   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031057_unfold
+219.7  238.533230447   28.8414 28.266187808    /home/wking/rsrch/data/scaled_unfold/20090203/20090203031725_unfold
+208.7  238.533230447   28.8414 14.736582977    /home/wking/rsrch/data/scaled_unfold/20090203/20090203031725_unfold
+267.4  238.533230447   28.8414 16.4420955747   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031725_unfold
+177.6  368.403149864   28.8414 65.2441978409   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031911_unfold
+236.3  368.403149864   28.8414 79.7592819456   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031911_unfold
+374.5  1357.2088083    28.8414 305.643423629   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031937_unfold
+256    1357.2088083    28.8414 119.529379747   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031937_unfold
+287.4  1357.2088083    28.8414 147.528597462   /home/wking/rsrch/data/scaled_unfold/20090203/20090203031937_unfold
+221.7  1357.2088083    28.8414 130.292045597   /home/wking/rsrch/data/scaled_unfold/20090203/20090203032640_unfold
+236.1  1357.2088083    28.8414 103.962194716   /home/wking/rsrch/data/scaled_unfold/20090203/20090203032640_unfold
+273    945.741609003   28.8414 97.8842565318   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033033_unfold
+244.1  945.741609003   28.8414 85.2586060516   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033033_unfold
+283.7  945.741609003   28.8414 55.4110008715   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033033_unfold
+321    945.741609003   28.8414 185.459929525   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033715_unfold
+250    945.741609003   28.8414 71.6682991302   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033715_unfold
+225.5  211.474252688   28.8414 12.546767412    /home/wking/rsrch/data/scaled_unfold/20090203/20090203033800_unfold
+169.1  211.474252688   28.8414 5.98895083612   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033800_unfold
+221    211.474252688   28.8414 7.68285960016   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033800_unfold
+284.7  211.474252688   28.8414 8.20520100429   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033800_unfold
+230.7  211.474252688   28.8414 14.2174140082   /home/wking/rsrch/data/scaled_unfold/20090203/20090203033915_unfold
+226.3  100.0   28.8414 19.4    /home/wking/rsrch/data/scaled_unfold/20090203/20090203035105_unfold
+221.2  100.0   28.8414 8.418   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035105_unfold
+213.7  100.0   28.8414 6.458   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035105_unfold
+235.7  100.0   28.8414 5.068   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035105_unfold
+220.8  100.0   28.8414 3.153   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035105_unfold
+236.4  100.0   28.8414 3.118   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035105_unfold
+289.4  945.741609003   28.8414 280.601535391   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035154_unfold
+270.7  945.741609003   28.8414 108.476562553   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035154_unfold
+274    945.741609003   28.8414 82.98882619     /home/wking/rsrch/data/scaled_unfold/20090203/20090203035154_unfold
+414.4  447.2135955     28.8414 205.002712177   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035328_unfold
+238    945.741609003   28.8414 197.281699638   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035911_unfold
+295    945.741609003   28.8414 150.656638314   /home/wking/rsrch/data/scaled_unfold/20090203/20090203035911_unfold
+198.4  100.0   28.8414 15.1    /home/wking/rsrch/data/scaled_unfold/20090203/20090203041905_unfold
+207.7  100.0   28.8414 9.036   /home/wking/rsrch/data/scaled_unfold/20090203/20090203041905_unfold
+212.1  100.0   28.8414 4.502   /home/wking/rsrch/data/scaled_unfold/20090203/20090203041905_unfold
+181.8  100.0   28.8414 2.738   /home/wking/rsrch/data/scaled_unfold/20090203/20090203041905_unfold
+250.7  100.0   28.8414 3.672   /home/wking/rsrch/data/scaled_unfold/20090203/20090203041905_unfold
+224.4  447.2135955     28.8414 178.93015956    /home/wking/rsrch/data/scaled_unfold/20090203/20090203042138_unfold
+217.4  447.2135955     28.8414 56.0358635161   /home/wking/rsrch/data/scaled_unfold/20090203/20090203042138_unfold
+210.9  447.2135955     28.8414 83.0475646844   /home/wking/rsrch/data/scaled_unfold/20090203/20090203042634_unfold
+282.4  447.2135955     28.8414 108.31513283    /home/wking/rsrch/data/scaled_unfold/20090203/20090203042634_unfold
+268.7  447.2135955     28.8414 47.0468702466   /home/wking/rsrch/data/scaled_unfold/20090203/20090203042634_unfold
+264.6  447.2135955     28.8414 33.089333931    /home/wking/rsrch/data/scaled_unfold/20090203/20090203042634_unfold
+310.6  447.2135955     28.8414 153.438984616   /home/wking/rsrch/data/scaled_unfold/20090203/20090203043209_unfold
+327.6  447.2135955     28.8414 92.9309851449   /home/wking/rsrch/data/scaled_unfold/20090203/20090203043209_unfold
+412.5  447.2135955     28.8414 83.5394996394   /home/wking/rsrch/data/scaled_unfold/20090203/20090203043209_unfold
+228.9  447.2135955     28.8414 119.048259122   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044027_unfold
+233.7  447.2135955     28.8414 57.0644547858   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044027_unfold
+242.1  945.741609003   28.8414 199.267757017   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044312_unfold
+234.9  945.741609003   28.8414 74.4676942929   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044312_unfold
+307.3  945.741609003   28.8414 90.8006518804   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044312_unfold
+168    211.474252688   28.8414 42.7389464682   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044536_unfold
+193.7  211.474252688   28.8414 24.4887184613   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044536_unfold
+253.6  211.474252688   28.8414 18.6562585721   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044536_unfold
+234.6  211.474252688   28.8414 58.5149257188   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044652_unfold
+265.5  211.474252688   28.8414 29.7544273532   /home/wking/rsrch/data/scaled_unfold/20090203/20090203044652_unfold
+203.9  945.741609003   28.8414 116.893662873   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045402_unfold
+262.8  945.741609003   28.8414 147.0628202     /home/wking/rsrch/data/scaled_unfold/20090203/20090203045402_unfold
+235.1  945.741609003   28.8414 58.200938618    /home/wking/rsrch/data/scaled_unfold/20090203/20090203045402_unfold
+302.9  945.741609003   28.8414 60.6314945532   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045402_unfold
+253.2  945.741609003   28.8414 40.8087504285   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045402_unfold
+258.7  945.741609003   28.8414 32.2687036992   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045402_unfold
+204.2  945.741609003   28.8414 111.881232345   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045518_unfold
+259.5  945.741609003   28.8414 71.9898512773   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045518_unfold
+268    945.741609003   28.8414 55.2880544623   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045518_unfold
+255.2  945.741609003   28.8414 45.6698622988   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045518_unfold
+291.1  945.741609003   28.8414 42.1422460972   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045518_unfold
+239.1  211.474252688   28.8414 17.579854626    /home/wking/rsrch/data/scaled_unfold/20090203/20090203045618_unfold
+252.7  211.474252688   28.8414 13.7585148799   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045618_unfold
+303.5  211.474252688   28.8414 13.7923507603   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045618_unfold
+212.5  945.741609003   28.8414 54.7111520808   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+230.8  945.741609003   28.8414 35.6733734916   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+200.4  945.741609003   28.8414 28.4006205184   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+217.7  945.741609003   28.8414 25.1756416317   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+263.9  945.741609003   28.8414 30.2637314881   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+244.5  945.741609003   28.8414 26.2443296498   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+277.2  945.741609003   28.8414 29.204500886    /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+203.9  945.741609003   28.8414 13.864571988    /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+233.2  945.741609003   28.8414 15.5763643003   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+275.6  945.741609003   28.8414 19.0188637571   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045647_unfold
+194.3  447.2135955     28.8414 98.9683686841   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045749_unfold
+247.7  447.2135955     28.8414 63.1912810442   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045749_unfold
+258.5  447.2135955     28.8414 41.13917865     /home/wking/rsrch/data/scaled_unfold/20090203/20090203045906_unfold
+260.7  447.2135955     28.8414 33.8809019951   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045906_unfold
+331    447.2135955     28.8414 34.3191713187   /home/wking/rsrch/data/scaled_unfold/20090203/20090203045906_unfold
+231.5  211.474252688   28.8414 13.5026310341   /home/wking/rsrch/data/scaled_unfold/20090203/20090203050203_unfold
+231.6  211.474252688   28.8414 9.74261882134   /home/wking/rsrch/data/scaled_unfold/20090203/20090203050203_unfold
+255.3  211.474252688   28.8414 10.0386827751   /home/wking/rsrch/data/scaled_unfold/20090203/20090203050203_unfold
+243.3  211.474252688   28.8414 6.9786503387    /home/wking/rsrch/data/scaled_unfold/20090203/20090203050203_unfold
+227    945.741609003   28.8414 253.364177052   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051002_unfold
+259.4  945.741609003   28.8414 110.651768253   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051002_unfold
+255.2  945.741609003   28.8414 62.4378610264   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051002_unfold
+257.3  2000.0  28.8414 1154.4  /home/wking/rsrch/data/scaled_unfold/20090203/20090203051658_unfold
+246.8  2000.0  28.8414 369.0   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051658_unfold
+298.4  2000.0  28.8414 332.4   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051658_unfold
+193.7  447.2135955     28.8414 12.8305580549   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+243.1  447.2135955     28.8414 13.5326833998   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+180    447.2135955     28.8414 6.90050577856   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+177.6  447.2135955     28.8414 5.88980305274   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+222.7  447.2135955     28.8414 7.81729364934   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+241.5  447.2135955     28.8414 8.20636947743   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+222.7  447.2135955     28.8414 6.29676742464   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+261.9  447.2135955     28.8414 8.44339268304   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+288.3  447.2135955     28.8414 8.13481530215   /home/wking/rsrch/data/scaled_unfold/20090203/20090203051942_unfold
+262.8  447.2135955     28.8414 55.6780926397   /home/wking/rsrch/data/scaled_unfold/20090203/20090203052936_unfold
+271.9  447.2135955     28.8414 48.5226751117   /home/wking/rsrch/data/scaled_unfold/20090203/20090203052936_unfold
+232    211.474252688   28.8414 18.3559651333   /home/wking/rsrch/data/scaled_unfold/20090203/20090203053331_unfold
+195.4  100.0   28.4299 9.235   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055154_unfold
+235.5  100.0   28.4299 15.22   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055154_unfold
+196    211.474252688   28.4299 27.1955888957   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055432_unfold
+195.3  211.474252688   28.4299 12.5298494718   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055432_unfold
+199.3  211.474252688   28.4299 9.04263904494   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055432_unfold
+202.1  211.474252688   28.4299 6.31462118526   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055432_unfold
+187.7  447.2135955     28.4299 129.691942695   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055458_unfold
+164.7  447.2135955     28.4299 31.8505522715   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055458_unfold
+213.1  100.0   28.4299 9.703   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055810_unfold
+217.2  100.0   28.4299 6.617   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055810_unfold
+230.8  100.0   28.4299 4.691   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055810_unfold
+212.5  100.0   28.4299 3.34    /home/wking/rsrch/data/scaled_unfold/20090203/20090203055810_unfold
+258    945.741609003   28.4299 264.618502199   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055930_unfold
+273.1  945.741609003   28.4299 146.117078591   /home/wking/rsrch/data/scaled_unfold/20090203/20090203055930_unfold
+212.5  211.474252688   28.4299 46.418598465    /home/wking/rsrch/data/scaled_unfold/20090203/20090203061239_unfold
+267.5  211.474252688   28.4299 29.8601644795   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061239_unfold
+270.5  945.741609003   28.4299 118.028552804   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061309_unfold
+277    945.741609003   28.4299 103.369557864   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061309_unfold
+263.6  945.741609003   28.4299 63.0904227366   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061309_unfold
+116.7  100.0   28.4299 1.764   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061335_unfold
+223.7  100.0   28.4299 4.17    /home/wking/rsrch/data/scaled_unfold/20090203/20090203061335_unfold
+225.4  100.0   28.4299 2.992   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061335_unfold
+204.6  100.0   28.4299 2.322   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061335_unfold
+236.2  100.0   28.4299 2.555   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061335_unfold
+272.5  447.2135955     28.4299 149.011570021   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061926_unfold
+293.8  447.2135955     28.4299 42.0067730253   /home/wking/rsrch/data/scaled_unfold/20090203/20090203061926_unfold
+206.4  2000.0  28.4299 402.8   /home/wking/rsrch/data/scaled_unfold/20090203/20090203063254_unfold
+266.6  2000.0  28.4299 257.4   /home/wking/rsrch/data/scaled_unfold/20090203/20090203063254_unfold
+287.1  2000.0  28.4299 171.62  /home/wking/rsrch/data/scaled_unfold/20090203/20090203063254_unfold
+238.2  2000.0  28.4299 101.3   /home/wking/rsrch/data/scaled_unfold/20090203/20090203063254_unfold
+118.5  212.936037324   24.221  3.66462920235   /home/wking/rsrch/data/scaled_unfold/20090212/20090212005005_unfold
+93.12  212.936037324   24.221  1.8604221581    /home/wking/rsrch/data/scaled_unfold/20090212/20090212005005_unfold
+129.7  212.936037324   24.221  2.55310308751   /home/wking/rsrch/data/scaled_unfold/20090212/20090212005005_unfold
+123.5  212.936037324   24.221  2.0258734591    /home/wking/rsrch/data/scaled_unfold/20090212/20090212005005_unfold
+141.1  212.936037324   24.221  2.12403697231   /home/wking/rsrch/data/scaled_unfold/20090212/20090212005005_unfold
+107.3  100.0   131.914 1.338   /home/wking/rsrch/data/scaled_unfold/20090212/20090212014843_unfold
+140.8  100.0   131.914 1.565   /home/wking/rsrch/data/scaled_unfold/20090212/20090212014843_unfold
+329.2  3000.0  131.914 332.1   /home/wking/rsrch/data/scaled_unfold/20090212/20090212015334_unfold
+230.2  3000.0  131.914 109.92  /home/wking/rsrch/data/scaled_unfold/20090212/20090212015334_unfold
+323.9  3000.0  131.914 165.6   /home/wking/rsrch/data/scaled_unfold/20090212/20090212015334_unfold
+391.1  3000.0  131.914 247.83  /home/wking/rsrch/data/scaled_unfold/20090212/20090212015334_unfold
+179.3  1408.87378093   131.914 36.5180084017   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025044_unfold
+269.9  1408.87378093   131.914 54.9883436697   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025044_unfold
+273.9  1408.87378093   131.914 47.3663365149   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025044_unfold
+264.3  1408.87378093   131.914 42.0267048851   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025044_unfold
+259.6  1408.87378093   131.914 35.0809571452   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025044_unfold
+279.9  310.723250595   131.914 56.6137762584   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025949_unfold
+317    310.723250595   131.914 39.7415037511   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025949_unfold
+246.8  310.723250595   131.914 15.1974741866   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025949_unfold
+278.5  310.723250595   131.914 15.4771251121   /home/wking/rsrch/data/scaled_unfold/20090212/20090212025949_unfold
+216.9  145.923280296   131.914 15.3073521031   /home/wking/rsrch/data/scaled_unfold/20090212/20090212030956_unfold
+217    145.923280296   131.914 8.24466533672   /home/wking/rsrch/data/scaled_unfold/20090212/20090212030956_unfold
+256.3  145.923280296   131.914 9.5141978753    /home/wking/rsrch/data/scaled_unfold/20090212/20090212030956_unfold
+348.6  1408.87378093   131.914 367.575169445   /home/wking/rsrch/data/scaled_unfold/20090212/20090212032408_unfold
+325.3  1408.87378093   131.914 182.026492496   /home/wking/rsrch/data/scaled_unfold/20090212/20090212032408_unfold
+295.8  1408.87378093   131.914 109.610380156   /home/wking/rsrch/data/scaled_unfold/20090212/20090212032408_unfold
+110.3  965.489384606   62.1096 28.4336623766   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044850_unfold
+141.3  3000.0  62.1096 196.23  /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+150.9  3000.0  62.1096 108.9   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+186.9  3000.0  62.1096 94.02   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+209.3  3000.0  62.1096 97.32   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+201.6  3000.0  62.1096 82.38   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+201.9  3000.0  62.1096 62.67   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+209.6  3000.0  62.1096 51.18   /home/wking/rsrch/data/scaled_unfold/20090212/20090212044927_unfold
+115.6  310.723250595   62.1096 25.1561543682   /home/wking/rsrch/data/scaled_unfold/20090212/20090212045258_unfold
+148.9  310.723250595   62.1096 13.597249446    /home/wking/rsrch/data/scaled_unfold/20090212/20090212045258_unfold
+154.5  310.723250595   62.1096 9.10419124243   /home/wking/rsrch/data/scaled_unfold/20090212/20090212045258_unfold
+159.3  2055.87483636   62.1096 92.7816313649   /home/wking/rsrch/data/scaled_unfold/20090212/20090212045406_unfold
+145.7  2055.87483636   62.1096 60.668866421    /home/wking/rsrch/data/scaled_unfold/20090212/20090212045406_unfold
+148.9  145.923280296   72.4529 10.0949725309   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051429_unfold
+184.5  145.923280296   72.4529 7.42019880305   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051429_unfold
+160.3  2055.87483636   72.4529 90.8285502704   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051610_unfold
+167    2055.87483636   72.4529 69.6530394559   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051610_unfold
+204.8  2055.87483636   72.4529 68.4195145541   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051610_unfold
+229.2  2055.87483636   72.4529 68.1728095737   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051610_unfold
+247.5  2055.87483636   72.4529 71.1127105897   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051610_unfold
+261    2055.87483636   72.4529 71.7294730406   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051610_unfold
+148.5  965.489384606   72.4529 74.3909570839   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051810_unfold
+179.4  965.489384606   72.4529 51.1033531272   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051810_unfold
+218.6  965.489384606   72.4529 51.0261139764   /home/wking/rsrch/data/scaled_unfold/20090212/20090212051810_unfold
+222.9  661.641776861   72.4529 69.4062223927   /home/wking/rsrch/data/scaled_unfold/20090212/20090212053527_unfold
+194.1  661.641776861   72.4529 30.4289053178   /home/wking/rsrch/data/scaled_unfold/20090212/20090212053527_unfold
+237.3  661.641776861   72.4529 30.9052873972   /home/wking/rsrch/data/scaled_unfold/20090212/20090212053527_unfold
+217.8  310.723250595   72.4529 41.3261923291   /home/wking/rsrch/data/scaled_unfold/20090212/20090212053725_unfold
+183.8  310.723250595   72.4529 14.5884566154   /home/wking/rsrch/data/scaled_unfold/20090212/20090212053725_unfold
+196.5  2055.87483636   72.4529 237.042368632   /home/wking/rsrch/data/scaled_unfold/20090212/20090212054838_unfold
+147.8  2055.87483636   72.4529 78.9455937162   /home/wking/rsrch/data/scaled_unfold/20090212/20090212054838_unfold
+149.7  2055.87483636   72.4529 53.1032470232   /home/wking/rsrch/data/scaled_unfold/20090212/20090212054838_unfold
+164.9  2055.87483636   72.4529 48.8681448603   /home/wking/rsrch/data/scaled_unfold/20090212/20090212054838_unfold
+134.3  2000.0  72.4529 82.58   /home/wking/rsrch/data/scaled_unfold/20090212/20090212061558_unfold
+157    2000.0  72.4529 65.52   /home/wking/rsrch/data/scaled_unfold/20090212/20090212061558_unfold
+194.3  2000.0  72.4529 79.1    /home/wking/rsrch/data/scaled_unfold/20090212/20090212061558_unfold
+208    2000.0  72.4529 70.06   /home/wking/rsrch/data/scaled_unfold/20090212/20090212061558_unfold
+217.1  2000.0  72.4529 56.02   /home/wking/rsrch/data/scaled_unfold/20090212/20090212061558_unfold
+173.8  447.2135955     72.4529 31.0455677996   /home/wking/rsrch/data/scaled_unfold/20090212/20090212062622_unfold
+214.3  447.2135955     72.4529 27.2934457334   /home/wking/rsrch/data/scaled_unfold/20090212/20090212062622_unfold
+212.3  447.2135955     72.4529 18.6756397481   /home/wking/rsrch/data/scaled_unfold/20090212/20090212062622_unfold
+151.4  447.2135955     72.4529 8.39867132349   /home/wking/rsrch/data/scaled_unfold/20090212/20090212062622_unfold
+155.1  2000.0  72.4529 75.02   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065445_unfold
+167.4  2000.0  72.4529 63.58   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065445_unfold
+123.7  2000.0  72.4529 56.14   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065905_unfold
+138.6  2000.0  72.4529 47.76   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065905_unfold
+199.8  2000.0  72.4529 65.58   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065905_unfold
+228.1  2000.0  72.4529 65.64   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065905_unfold
+218.4  2000.0  72.4529 51.84   /home/wking/rsrch/data/scaled_unfold/20090212/20090212065905_unfold
+142.2  100.0   72.4529 4.783   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072150_unfold
+217.8  100.0   72.4529 6.317   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072150_unfold
+214.5  100.0   72.4529 4.249   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072150_unfold
+134.7  945.741609003   72.4529 50.5971760817   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+163.1  945.741609003   72.4529 35.1532156066   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+205    945.741609003   72.4529 38.028270098    /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+188.9  945.741609003   72.4529 26.8590616957   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+161.7  945.741609003   72.4529 16.9949767138   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+157    945.741609003   72.4529 12.9282877951   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+142    945.741609003   72.4529 9.40823752636   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+219.2  945.741609003   72.4529 16.8152858081   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+214    945.741609003   72.4529 14.7441116844   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+156.1  945.741609003   72.4529 7.76643009313   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+151.7  945.741609003   72.4529 6.68450169243   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072249_unfold
+171.7  211.474252688   72.4529 27.0687043441   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072650_unfold
+120.8  211.474252688   72.4529 5.72249327774   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072650_unfold
+200.8  945.741609003   72.4529 120.014610182   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072725_unfold
+110.7  945.741609003   72.4529 22.4235335495   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072725_unfold
+213.2  2000.0  72.4529 223.2   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072911_unfold
+158.7  2000.0  72.4529 66.92   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072911_unfold
+147.3  211.474252688   72.4529 9.67283231795   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072955_unfold
+176.1  211.474252688   72.4529 8.2686432801    /home/wking/rsrch/data/scaled_unfold/20090212/20090212072955_unfold
+224.2  211.474252688   72.4529 8.97919676913   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072955_unfold
+214.1  211.474252688   72.4529 7.05266632714   /home/wking/rsrch/data/scaled_unfold/20090212/20090212072955_unfold
+181.5  2000.0  72.4529 90.96   /home/wking/rsrch/data/scaled_unfold/20090212/20090212073046_unfold
+218.4  2000.0  72.4529 71.94   /home/wking/rsrch/data/scaled_unfold/20090212/20090212073046_unfold
+202.3  2000.0  72.4529 58.36   /home/wking/rsrch/data/scaled_unfold/20090212/20090212073046_unfold
+146.3  945.741609003   28.5347 32.4105649405   /home/wking/rsrch/data/scaled_unfold/20090214/20090214220914_unfold
+167.6  945.741609003   28.5347 24.0218368687   /home/wking/rsrch/data/scaled_unfold/20090214/20090214220914_unfold
+157.2  945.741609003   28.5347 17.2124972839   /home/wking/rsrch/data/scaled_unfold/20090214/20090214220914_unfold
+173.1  945.741609003   28.5347 16.0114054404   /home/wking/rsrch/data/scaled_unfold/20090214/20090214220914_unfold
+150.1  945.741609003   28.5347 18.9526618444   /home/wking/rsrch/data/scaled_unfold/20090214/20090214221045_unfold
+166    945.741609003   28.5347 16.6166800702   /home/wking/rsrch/data/scaled_unfold/20090214/20090214221045_unfold
+159.6  945.741609003   28.5347 13.4200734318   /home/wking/rsrch/data/scaled_unfold/20090214/20090214221045_unfold
+143.4  945.741609003   28.5347 10.0248610554   /home/wking/rsrch/data/scaled_unfold/20090214/20090214221045_unfold
+215.8  2000.0  28.7156 156.58  /home/wking/rsrch/data/scaled_unfold/20090219/20090219174814_unfold
+184.2  2000.0  28.7156 58.26   /home/wking/rsrch/data/scaled_unfold/20090219/20090219174814_unfold
+150.6  945.741609003   28.7156 64.5846944788   /home/wking/rsrch/data/scaled_unfold/20090219/20090219182043_unfold
+174.1  945.741609003   28.7156 44.752492938    /home/wking/rsrch/data/scaled_unfold/20090219/20090219182043_unfold
+161.7  945.741609003   28.7156 20.8252302302   /home/wking/rsrch/data/scaled_unfold/20090219/20090219182043_unfold
+159.5  945.741609003   28.7156 16.2478408427   /home/wking/rsrch/data/scaled_unfold/20090219/20090219182043_unfold
+271.2  100.0   123.857 11.75   /home/wking/rsrch/data/scaled_unfold/20090219/20090219190020_unfold
+273.2  100.0   123.857 7.565   /home/wking/rsrch/data/scaled_unfold/20090219/20090219190020_unfold
+215.8  100.0   123.857 3.825   /home/wking/rsrch/data/scaled_unfold/20090219/20090219190020_unfold
+321.9  945.741609003   123.857 178.461441619   /home/wking/rsrch/data/scaled_unfold/20090219/20090219190927_unfold
+302.7  945.741609003   123.857 107.814543426   /home/wking/rsrch/data/scaled_unfold/20090219/20090219190927_unfold
+274.7  211.474252688   35.1535 14.5346253872   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+264.2  211.474252688   35.1535 9.52902982612   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+153.6  211.474252688   35.1535 4.1618132929    /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+263.3  211.474252688   35.1535 8.07620171015   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+255.7  211.474252688   35.1535 5.88321370978   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+246.6  211.474252688   35.1535 5.30377425742   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+292.4  211.474252688   35.1535 5.94242650053   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+289.4  211.474252688   35.1535 6.28290004736   /home/wking/rsrch/data/scaled_unfold/20090222/20090222211752_unfold
+230.8  211.474252688   35.1535 20.8640497702   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+256    211.474252688   35.1535 19.0897807901   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+207.5  211.474252688   35.1535 10.3706973518   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+215.8  211.474252688   35.1535 7.34661553838   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+225.7  211.474252688   35.1535 6.23003148419   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+298.8  211.474252688   35.1535 7.82243260693   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+324.5  211.474252688   35.1535 8.23480739967   /home/wking/rsrch/data/scaled_unfold/20090222/20090222213222_unfold
+201.5  945.741609003   35.1535 248.162598202   /home/wking/rsrch/data/scaled_unfold/20090222/20090222214527_unfold
+262.9  945.741609003   35.1535 134.579030961   /home/wking/rsrch/data/scaled_unfold/20090222/20090222214527_unfold
+237.2  945.741609003   35.1535 627.12126093    /home/wking/rsrch/data/scaled_unfold/20090222/20090222220126_unfold
+299.1  945.741609003   35.1535 259.133200867   /home/wking/rsrch/data/scaled_unfold/20090222/20090222220126_unfold
+319.5  447.2135955     132.727 248.650759098   /home/wking/rsrch/data/scaled_unfold/20090222/20090222224649_unfold
+378.4  447.2135955     132.727 110.103987212   /home/wking/rsrch/data/scaled_unfold/20090222/20090222224649_unfold
+258    945.741609003   132.727 82.5726998821   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230350_unfold
+276.8  945.741609003   132.727 55.7041807703   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230350_unfold
+366.8  945.741609003   132.727 69.1620838664   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230350_unfold
+338.7  945.741609003   132.727 46.1427331033   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230350_unfold
+316.9  2000.0  132.727 217.6   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230533_unfold
+386    2000.0  132.727 230.8   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230533_unfold
+284.8  2000.0  132.727 94.88   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230533_unfold
+326.7  945.741609003   132.727 61.9460753897   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230646_unfold
+347.7  945.741609003   132.727 49.263680413    /home/wking/rsrch/data/scaled_unfold/20090222/20090222230646_unfold
+247.4  945.741609003   132.727 30.8973783661   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230646_unfold
+444    945.741609003   132.727 65.1899691086   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230646_unfold
+247.8  447.2135955     132.727 32.897032085    /home/wking/rsrch/data/scaled_unfold/20090222/20090222230924_unfold
+388.7  447.2135955     132.727 45.9735576174   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230924_unfold
+505.5  447.2135955     132.727 59.6135722802   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230924_unfold
+355    447.2135955     132.727 13.6802638863   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230924_unfold
+342.6  447.2135955     132.727 10.9925101774   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230924_unfold
+336.5  447.2135955     132.727 10.0444173549   /home/wking/rsrch/data/scaled_unfold/20090222/20090222230924_unfold
+190.9  211.474252688   132.727 8.47165856268   /home/wking/rsrch/data/scaled_unfold/20090222/20090222231202_unfold
+315.9  211.474252688   132.727 11.2800366384   /home/wking/rsrch/data/scaled_unfold/20090222/20090222231202_unfold
+357.2  211.474252688   132.727 14.0736115164   /home/wking/rsrch/data/scaled_unfold/20090222/20090222231202_unfold
+500.8  211.474252688   132.727 20.5827890141   /home/wking/rsrch/data/scaled_unfold/20090222/20090222231202_unfold
+377.5  447.2135955     132.727 101.741092976   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232706_unfold
+396    447.2135955     132.727 73.9691286957   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232706_unfold
+259.9  211.474252688   132.727 9.31544083091   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232816_unfold
+357.7  211.474252688   132.727 12.390276465    /home/wking/rsrch/data/scaled_unfold/20090222/20090222232816_unfold
+410.8  211.474252688   132.727 14.4225440333   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232816_unfold
+410.4  211.474252688   132.727 11.5782153347   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232816_unfold
+252.7  447.2135955     132.727 35.1062672468   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232834_unfold
+282.3  447.2135955     132.727 25.7326702851   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232834_unfold
+269.3  100.0   132.727 5.993   /home/wking/rsrch/data/scaled_unfold/20090222/20090222232921_unfold
+202.3  945.741609003   132.727 57.1227931838   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+314.4  945.741609003   132.727 70.0794532271   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+303.4  945.741609003   132.727 52.5926908767   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+339    945.741609003   132.727 55.5150324485   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+397.2  945.741609003   132.727 62.2297978724   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+402.9  945.741609003   132.727 57.7280678135   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+397.9  945.741609003   132.727 48.2990239718   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+380.4  945.741609003   132.727 36.2502758731   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+396.3  945.741609003   132.727 35.6355438272   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+374.6  945.741609003   132.727 29.8570625962   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+443.3  945.741609003   132.727 36.5718280201   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233019_unfold
+343.5  945.741609003   132.727 146.589949395   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233613_unfold
+400.8  945.741609003   132.727 146.211652752   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233613_unfold
+387.4  2000.0  132.727 294.6   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233757_unfold
+393.2  2000.0  132.727 255.6   /home/wking/rsrch/data/scaled_unfold/20090222/20090222233757_unfold
+304.9  2000.0  132.727 227.6   /home/wking/rsrch/data/scaled_unfold/20090222/20090222235356_unfold
+397.7  2000.0  132.727 224.4   /home/wking/rsrch/data/scaled_unfold/20090222/20090222235356_unfold
+248.1  447.2135955     132.727 45.9735576174   /home/wking/rsrch/data/scaled_unfold/20090222/20090223000513_unfold
+313.3  447.2135955     132.727 39.1669666939   /home/wking/rsrch/data/scaled_unfold/20090222/20090223000513_unfold
+373.6  945.741609003   132.727 76.6239851614   /home/wking/rsrch/data/scaled_unfold/20090222/20090223001009_unfold
+449.6  945.741609003   132.727 117.933978643   /home/wking/rsrch/data/scaled_unfold/20090222/20090223001009_unfold
+346.7  945.741609003   132.727 49.5190306474   /home/wking/rsrch/data/scaled_unfold/20090222/20090223001009_unfold
+332.2  945.741609003   132.727 40.439911201    /home/wking/rsrch/data/scaled_unfold/20090222/20090223001009_unfold
+396.2  945.741609003   132.727 43.8256661612   /home/wking/rsrch/data/scaled_unfold/20090222/20090223001009_unfold
diff --git a/tex/src/figures/cantilever-data/extract_f_v_k_data.sh b/tex/src/figures/cantilever-data/extract_f_v_k_data.sh
new file mode 100755 (executable)
index 0000000..976d4b7
--- /dev/null
@@ -0,0 +1,92 @@
+#!/bin/bash
+#
+# Extract v,f,k data.  Ignore index, temperature, etc.
+
+set -e # exit immediately on failed command
+set -o pipefail # pipes fail if any stage fails
+
+BASE_DATA_DIR="$HOME/rsrch/data"
+SCALED_DATA_DIR="scaled_unfold"
+CALIBRATION_DATA_DIR="calibration"
+
+EXTRACT_UNFOLDS="$HOME/src/lab/analysis/sawmodel/extract_unfolds.py"
+SCALE_UNFOLD="$HOME/src/lab/analysis/sawtooth_worming-0.1/scale_unfold.py"
+
+DATA_DAYS[8]=20090114
+#DATA_DAYS[9]=20090128   # No good sawteeth
+DATA_DAYS[10]=20090203
+#DATA_DAYS[11]=20090211  # No good sawteeth
+DATA_DAYS[12]=20090212
+DATA_DAYS[13]=20090214
+DATA_DAYS[14]=20090219
+DATA_DAYS[15]=20090222
+
+echo -e "#Unfolding force (pN)\tUnfolding rate (nm/s)\tSpring constant (pN/nm)\tLoading rate (pN/s)\tDatafile"
+
+function NotEmpty ()
+{
+    NAME=$1
+    VALUE=$2
+    if [ -z "$VALUE" ]; then
+       echo "Error parsing"
+       echo "$UNFOLD_LINE"
+       echo "$NAME not defined.  File format must have changed."
+       exit 1
+    fi
+}
+
+FIELD_INDEX=""
+function FieldIndex ()
+{
+    LINE="$1"
+    TITLE_REGEXP="$2"
+    FIELD_INDEX=`echo "$LINE" | sed 's/\t/\n/g' | grep -n "$TITLE_REGEXP" | cut -d: -f1`
+}
+
+for DAY in ${DATA_DAYS[*]}
+do
+    echo "---Processing day $DAY---" >&2
+    DATAFILE_FIELD=""
+    FORCE_FIELD=""
+    LAST_DATAFILE="" # Caching for efficiency
+    while read UNFOLD_LINE
+    do
+       echo "Parsing '$UNFOLD_LINE'" >&2
+       if [ "${UNFOLD_LINE:0:1}" == "#" ] ; then
+           # reading the header line
+           FieldIndex "$UNFOLD_LINE" "Datafile"
+           DATAFILE_FIELD="$FIELD_INDEX"
+           FieldIndex "$UNFOLD_LINE" "Force (pN)"
+           FORCE_FIELD="$FIELD_INDEX"
+           FieldIndex "$UNFOLD_LINE" "Stiffness (pN/nm)"
+           STIFFNESS_FIELD="$FIELD_INDEX"
+           NotEmpty 'Datafile field' "$DATAFILE_FIELD"
+           NotEmpty 'Force field' "$FORCE_FIELD"
+           NotEmpty 'Stiffness field' "$STIFFNESS_FIELD"
+           continue
+       fi
+       NotEmpty "$DAY extraction header" "$DATAFILE_FIELD"
+       DATAFILE=`echo "$UNFOLD_LINE" | cut -f $DATAFILE_FIELD`
+        NotEmpty 'Datafile' "$DATAFILE"
+       DATAFILE="$BASE_DATA_DIR/$SCALED_DATA_DIR/$DATAFILE"
+       UNFOLDING_FORCE=`echo "$UNFOLD_LINE" | cut -f $FORCE_FIELD`
+       NotEmpty 'Unfolding force' "$UNFOLDING_FORCE"
+       STIFFNESS=`echo "$UNFOLD_LINE" | cut -f $STIFFNESS_FIELD`
+       NotEmpty 'Stiffness' "$STIFFNESS"
+       if [ "$DATAFILE" != "$LAST_DATAFILE" ]; then
+           UNFOLDING_RATE=`sed -n 's/^Unfold rate (nm\/s):\t//p' "${DATAFILE}_param"`
+           NotEmpty 'Unfolding rate' "$UNFOLDING_RATE"
+           SPRING_CONSTANT_STRING=`$SCALE_UNFOLD --print-k "$DATAFILE"`
+           NotEmpty 'Spring constant string' "$SPRING_CONSTANT_STRING"
+           SPRING_CONSTANT_N_PER_M=`echo "$SPRING_CONSTANT_STRING" | sed -n 's/ N\/m$//p'`
+           NotEmpty 'Spring constant (N/m)' "$SPRING_CONSTANT_N_PER_M"
+           SPRING_CONSTANT=`echo "print ($SPRING_CONSTANT_N_PER_M * 1e3)" | python`
+           NotEmpty 'Spring constant (pN/nm)' "$SPRING_CONSTANT"
+       fi
+       LOADING_RATE=`python -c "print $STIFFNESS * $UNFOLDING_RATE"` # pN/s
+       echo -e "$UNFOLDING_FORCE\t$UNFOLDING_RATE\t$SPRING_CONSTANT\t$LOADING_RATE\t$DATAFILE"
+       LAST_DATAFILE="$DATAFILE"
+    done < <(cd $BASE_DATA_DIR/$SCALED_DATA_DIR/ && $EXTRACT_UNFOLDS ${DAY}_select_fit.xml)
+done
+
+exit 0
diff --git a/tex/src/figures/cantilever-data/fit_data.py b/tex/src/figures/cantilever-data/fit_data.py
new file mode 100755 (executable)
index 0000000..f87a6e0
--- /dev/null
@@ -0,0 +1,45 @@
+#!/usr/bin/python
+
+import sys
+from scipy.stats import linregress
+from math import log10
+
+fname = sys.argv[1]
+
+WEIGHTHACK=True
+LOGX=True
+
+x = []
+y = []
+w = [] # weights
+
+print >> sys.stderr, "Fitting data from", fname
+for line in file(fname, 'r'):
+    fields = [float(a) for a in line.strip().split()]
+    if WEIGHTHACK==True:
+        w = int(fields[3])
+        x.extend([fields[0]]*w)
+        y.extend([fields[1]]*w)
+    else:
+        x.append(fields[0])
+        y.append(fields[1])
+        w.append(fields[3])
+
+if LOGX == True:
+    slope,intercept,r,two_tailed_prob,stderr_of_the_estimate = linregress([log10(xi) for xi in x],y)
+else:
+    slope,intercept,r,two_tailed_prob,stderr_of_the_estimate = linregress(x,y)
+#print slope,intercept,r,two_tailed_prob,stderr_of_the_estimate
+
+def f(x):
+    if LOGX == True:
+        x = log10(x)
+    val = intercept+slope*x
+    return val
+
+x.sort()
+xmin = x[0]
+xmax = x[-1]
+print '#fitted line'
+print '%g\t%g' % (xmin, f(xmin))
+print '%g\t%g' % (xmax, f(xmax))
diff --git a/tex/src/figures/cantilever-data/get_loading_rates.py b/tex/src/figures/cantilever-data/get_loading_rates.py
new file mode 100755 (executable)
index 0000000..8010910
--- /dev/null
@@ -0,0 +1,40 @@
+#!/usr/bin/python
+
+DATA="./data"
+UNFOLD_FORCE_COL=0
+UNFOLD_RATE_COL=1
+SPRING_CONST_COL=2
+LOAD_RATE_COL=3
+
+from avg_data import KCUTS
+
+Ks = [24.33, None, 131.98]   # HACK!
+OFILEs = [None, None, None]
+
+for i,K in enumerate(Ks):
+    if K == None:
+        continue
+    OFILEs[i] = file('loading_rate_%.2f' % K, 'w')
+
+for line in file(DATA, 'r'):
+    if line[0] == '#':
+        continue
+    fields = line.strip().split()
+    F = float(fields[UNFOLD_FORCE_COL])
+    V = float(fields[UNFOLD_RATE_COL])
+    K = float(fields[SPRING_CONST_COL])
+    L = float(fields[LOAD_RATE_COL])
+    
+    ispring = 0
+    while ispring < len(KCUTS) and K > KCUTS[ispring]:
+        ispring += 1
+    if ispring == 1:
+        continue # drop middle spring constants
+    
+    print >> OFILEs[ispring], "%g\t%g" % (L, F)
+    
+for i,K in enumerate(Ks):
+    if K == None:
+        continue
+    OFILEs[i].close()
+    OFILEs[i] = None
diff --git a/tex/src/figures/cantilever-data/make_plots.sh b/tex/src/figures/cantilever-data/make_plots.sh
new file mode 100755 (executable)
index 0000000..43baa96
--- /dev/null
@@ -0,0 +1,41 @@
+#!/bin/bash
+
+KFILE='./spring_constants'
+AVGFILE='./averaged_data'
+GPSPLIT='./plot_splits.gp'
+
+PLOTS=""
+while read LINE
+do
+    Kprecise=`echo "$LINE" | sed 's/.*= \([0-9.]*\) +.*/\1/'`
+    K=`python -c "print '%.2f' % $Kprecise"`
+    FILE="v_dep_$K"
+    echo "clearing $FILE"
+    > "$FILE"
+    rm "$FILE.fit"
+    PLOTS="$PLOTS, '$FILE' using 1:2:(sqrt(\$4)) with points pt 6 pointsize variable t '$K (pN/nm)'"
+    PLOTS="$PLOTS, '$FILE.fit' using 1:2 with lines notitle"
+done < <(tac "$KFILE")
+PLOTS="${PLOTS:2}" # remove leading ' ,'
+
+while read DATA
+do
+    V=`echo "$DATA" | awk '{print $1}'`
+    K=`echo "$DATA" | awk '{print $2}'`
+    F=`echo "$DATA" | awk '{print $3}'`
+    Fstd=`echo "$DATA" | awk '{print $4}'`
+    N=`echo "$DATA" | awk '{print $5}'`
+    FILE="v_dep_$K"
+    echo "adding $V, $F, $Fstd, $N to $FILE"
+    echo -e "$V\t$F\t$Fstd\t$N" >> "$FILE"
+done < <(cat "$AVGFILE" | grep -v '^#')
+
+for FILE in v_dep_*
+do
+    ./fit_data.py "$FILE" > "$FILE.fit"
+done
+
+echo "set xlabel 'Pulling speed (nm/s)'
+set ylabel 'Unfolding force (pN)'
+set logscale x;
+plot $PLOTS" > "$GPSPLIT"
diff --git a/tex/src/figures/cantilever-data/makeplotsplits.sh b/tex/src/figures/cantilever-data/makeplotsplits.sh
new file mode 100755 (executable)
index 0000000..38d148e
--- /dev/null
@@ -0,0 +1,6 @@
+#!/bin/bash
+
+gnuplot plot_splits.gp 
+epstopdf plot_splits.eps 
+pdftk plot_splits.pdf cat 1E output plot_splits_rotated.pdf
+cp plot_splits_rotated.pdf ../../posters/2009_biophysical_society_poster/plot_splits.pdf
diff --git a/tex/src/figures/cantilever-data/plot_avg.gp b/tex/src/figures/cantilever-data/plot_avg.gp
new file mode 100644 (file)
index 0000000..f9ce6a6
--- /dev/null
@@ -0,0 +1,10 @@
+
+# see: show palette rgbformulae
+set palette rgbformulae 3,13,-3
+set xlabel 'Pulling speed (nm/s)'
+set ylabel 'Unfolding force (pN)'
+set zlabel 'Spring constant (pN/nm)'
+set view map
+set logscale x
+set style line 1 lt 1 lw 1 pt 2 ps 3
+splot 'averaged_data' using 1:3:2:(sqrt($5)) with points pt 6 pointsize variable lt palette z
diff --git a/tex/src/figures/cantilever-data/plot_raw.gp b/tex/src/figures/cantilever-data/plot_raw.gp
new file mode 100644 (file)
index 0000000..d8772ac
--- /dev/null
@@ -0,0 +1,8 @@
+# see: show palette rgbformulae
+set palette rgbformulae 3,13,-3
+set xlabel 'Pulling speed (nm/s)'
+set ylabel 'Unfolding force (pN)'
+set zlabel 'Spring constant (pN/nm)'
+set view map
+set logscale x
+splot 'data' using 2:1:3 with points palette z
diff --git a/tex/src/figures/cantilever-data/plot_splits.pdf b/tex/src/figures/cantilever-data/plot_splits.pdf
deleted file mode 100644 (file)
index a9702d0..0000000
Binary files a/tex/src/figures/cantilever-data/plot_splits.pdf and /dev/null differ
diff --git a/tex/src/figures/cantilever-data/plot_splits_nolandscape.pdf b/tex/src/figures/cantilever-data/plot_splits_nolandscape.pdf
deleted file mode 100644 (file)
index 0b1f403..0000000
Binary files a/tex/src/figures/cantilever-data/plot_splits_nolandscape.pdf and /dev/null differ
diff --git a/tex/src/figures/cantilever-data/plot_splits_rotated.pdf b/tex/src/figures/cantilever-data/plot_splits_rotated.pdf
deleted file mode 100644 (file)
index 28a2e30..0000000
Binary files a/tex/src/figures/cantilever-data/plot_splits_rotated.pdf and /dev/null differ
diff --git a/tex/src/figures/cantilever-sim/Makefile b/tex/src/figures/cantilever-sim/Makefile
new file mode 100644 (file)
index 0000000..7a28e49
--- /dev/null
@@ -0,0 +1,14 @@
+ASYPROCESS = ../asy/asyprocess \
+  --texinputs=../.. \
+  --pretex='\documentclass{drexel-thesis} \input{packages}' \
+  --
+
+FIGS = loading-rate v-dep i-dep 
+
+all : $(FIGS:%=%_.tex)
+
+clean :
+       rm -f $(FIGS:%=%_*) $(FIGS:%=%-*)
+
+%_.tex : %.asy
+       $(ASYPROCESS) $(patsubst %.asy, %, $<)
index 47583a27ab4361b42765f817c96abd055cdedb48..f874527fc0842617219bda56ee5b67516a596b59 100644 (file)
@@ -42,6 +42,7 @@
 %  Listed in alphebetical order.
 @String{AIP = "AIP"}
 @String{APS = "American Physical Society"}
+@String{SPRINGER = "Springer Science + Business Media, LLC"}
 
 %  Define some journal name shortcuts (useful for consistency)
 %  Listed in alphebetical order.
 @String{SNg = "Ng, Sean P."}
 @String{SNie = "Nie, S."}
 @String{RNome = "Nome, Rene A."}
+@String{ANoy = "Noy, Aleksandr"}
 @String{AOberhauser = "Oberhauser, Andres F."}
 @String{FOesterhelt = "Oesterhelt, Filipp"}
 @String{TOhashi = "Ohashi, Tomoo"}
   doi =          "10.1110/ps.0224602",
   URL =          "http://www.proteinscience.org/cgi/content/abstract/11/12/2759",
   eprint =       "http://www.proteinscience.org/cgi/reprint/11/12/2759.pdf",
-  note =         "Introduces unfolding-order effects on average unfolding force.",
+  note =         "Introduces unfolding-order \emph{scaffold effect} on average unfolding force.",
   project =      "sawtooth simulation",
 }
 
@@ -6002,6 +6004,7 @@ doi = {10.1063/1.439715}
   doi =          "10.1038/nsb968",
   url = "http://www.nature.com/nsmb/journal/v10/n9/abs/nsb968.html",
   eprint = "http://www.nature.com/nsmb/journal/v10/n9/pdf/nsb968.pdf",
+  note = "Another scaffold effect paper.  TODO: details",
 }
 
 @article{friddle08,
@@ -7755,7 +7758,7 @@ doi = "DOI: 10.1016/j.ijbiomac.2009.12.001",
 url = "http://www.sciencedirect.com/science/article/B6T7J-4XWMND2-1/2/7ef768562b4157fc201d450553e5de5e",
 }
 
-@Article {metropolis87,
+@Article{metropolis87,
   author = "Nicholas Metropolis",
   title = "The Beginning of the Monte Carlo Method",
   journal = "Los Alamos Science",
@@ -7765,3 +7768,12 @@ url = "http://www.sciencedirect.com/science/article/B6T7J-4XWMND2-1/2/7ef768562b
   publisher = "Los Alamos National Laboratory",
   url = "http://library.lanl.gov/cgi-bin/getfile?15-12.pdf",
 }
+
+@Book{noy08,
+  editor = ANoy,
+  title = "Handbook of Molecular Force Spectroscopy",
+  publisher = SPRINGER,
+  year = 2008,
+  ISBN="978-0-387-49987-1",
+  note="The first book about force spectroscopy.  Discusses the scaffold effect in section 8.4.1.",
+}