Polishing up fits and referee responses.

1 files changed, 14 insertions, 14 deletions

diff --git a/fig-stiffnesses.gplot b/fig-stiffnesses.gplot
index 52762dd..93b1f37 100644
--- a/fig-stiffnesses.gplot
+++ b/fig-stiffnesses.gplot
@@ -9,11 +9,11 @@ cc5 = "#000000"
 
 Tc = 17.26
 
-a1 = 72.9479505025133
-b1 = 0.01227234014115862
-c1 = 2675.3172792640935
-d1 = 8.418493827841177
-e1 = 19.456906054458745
+a1 = 71.13597475161484
+b1 = 0.010425748328804992
+c1 = 1665.064389018724
+d1 = 6.283065722894796
+e1 = 14.579664883988315
 
 C10(T) = a1 - b1 * T
 C1(T) = C10(T) / (1 + 1 / (1 / d1 + abs(T - Tc) / c1) / C10(T))
@@ -63,16 +63,16 @@ set xrange [0:300]
 
 set title '(a)' offset 12,-2.7
 set ylabel '$C_{\mathrm{B_{2\mathrm g}}} / \mathrm{GPa}$' offset 3
-set yrange [140:145]
-set ytics 141,1,144 offset 0.5 mirror
+set yrange [140.2:144.7]
+set ytics 139,1,146 offset 0.5 mirror
 plot "data/c66.dat" using 1:(100 * $2) with lines lw 6 lc rgb cc3, \
      C20(x) dt 3 lw 4 lc rgb cc4
 
 set title '(b)' offset 12,-2.7
 unset ylabel
 set y2label '$C_{\mathrm{B_{1\mathrm g}}} / \mathrm{GPa}$' offset -4 rotate by -90
-set yrange [65.05:65.7]
-set y2tics 62.1,0.1,65.6 offset -0.5 mirror
+set yrange [65.05:65.68]
+set y2tics 62.1,0.1,66 offset -0.5 mirror
 set samples 1000
 plot "data/c11mc12.dat" using 1:(100 * $2) with points pt 6 ps 0.4 lc rgb cc3, \
      x > Tc ? C12(x) : 1 / 0 lw 3 lc rgb cc5, \
@@ -92,7 +92,7 @@ set ytics 63,1,74
 unset y2label
 set ylabel '$\textcolor{mathc3}{C_{\mathrm{B_{1\mathrm g}}}},\textcolor{mathc4}{C_{\mathrm{B_{1\mathrm g}}}^0} / \mathrm{GPa}$' offset 3
 
-set yrange [64.5:73.5]
+set yrange [64.5:71.5]
 plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 6 lc rgb cc3, \
      C10(x) dt 3 lw 4 lc rgb cc4
 
@@ -102,8 +102,8 @@ set y2label '$\Big[C^0_{\mathrm{B_{1\mathrm g}}}(C^0_{\mathrm{B_{1\mathrm g}}}/C
 set title '(d)' offset 12,-2.7
 set format y2 '\scriptsize$%0.2f$'
 set format y '\scriptsize$%0.2f$'
-set yrange [0.11:0.25]
-set y2tics 0.05,0.025,0.25 offset -0.7 mirror
+set yrange [0.14:0.38]
+set y2tics 0.05,0.05,0.5 offset -0.7 mirror
 plot "data/c11mc12.dat" using 1:(1 / (C10($1)*(C10($1) / (100 * $2) - 1))) with lines lw 6 lc rgb cc3, \
      x > Tc ? 1/(C10(x) * (C10(x) / C12(x) - 1)) : 1 / 0 lw 3 lc cc5, \
      x <= Tc ? 1/(C10(x) * (C10(x) / C12(x) - 1)) : 1 / 0 dt (2,1) lw 3 lc cc5, \
@@ -127,8 +127,8 @@ plot "data/c11mc12.dat" using 1:(100 * $2) with points pt 2 ps 0.15 lc rgb cc3,
 
 unset margins
 set origin 0.5,0.26
-set yrange [0.1155:0.123]
-set y2tics 0.105,0.0025,0.131
+set yrange [0.154:0.166]
+set y2tics 0.155,0.005,0.167
 set format y2 '\scriptsize$%0.3f$'
 
 plot "data/c11mc12.dat" using 1:(1 / (C10($1)*(C10($1) / (100 * $2) - 1))) with points pt 2 ps 0.15 lc rgb cc3, \