cc1 = "#5e81b5"
cc2 = "#e19c24"
cc3 = "#8fb032"
cc4 = "#eb6235"

set terminal epslatex size 8.68cm, 8.2cm standalone
set output "fig-stiffnesses.tex"

set multiplot  layout 3, 2  margins 0.13, 0.9, 0.1, 0.998  spacing 0.01, 0.01

set nokey
set xrange [0:315]
set arrow 1 from 17.26,graph 0 to 17.26,graph 1 nohead lw 4 lc rgb cc2

lam(T) = 71.1212 - 0.0104105 * T

set format x ""
set xtics 0, 100, 300 offset 0,0.5
set mxtics 2
set ylabel '\tiny $C / \mathrm{GPa}$' offset 4.5
set format y '\tiny $%.1f$'
set format y2 '\tiny $%.1f$'

set yrange [140:145]
set ytics 141,1,144 offset 0.5
set title '\tiny (e) $C_{\mathrm{B_{2\mathrm g}}}$' offset -1,-6
plot "data/c66.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3, \
     144.345 - 0.019492 * x**2 / (120.462 + x) dt 3 lw 4 lc rgb cc4

set ylabel ''

set yrange [65.05:65.7]
set y2tics 62.1,0.1,65.6 offset -0.5 mirror
set title '\tiny (c) $C_{\mathrm{B_{1\mathrm g}}}$' offset -1,-6
plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3

set format x '\tiny $%.0f$'
set ylabel '\tiny $C / \mathrm{GPa}$'
set xlabel '\tiny $T / \mathrm K$' offset 0,1
unset y2tics

set yrange [65:71.5]
set ytics 63,1,72
set title '\tiny (c) $C_{\mathrm{B_{1\mathrm g}}}$' offset -1,-6
plot "data/c11mc12.dat" using 1:(100 * $2) with lines lw 3 lc rgb cc3, \
     lam(x) dt 3 lw 4 lc rgb cc4

set ylabel ''
set y2label '\tiny $\log[(C^0 - C) / \mathrm{GPa}]$' offset -5 rotate by -90
set xlabel '\tiny $\log(T / \mathrm K+D_\perp q_*^4/a)$' offset 0,1

set yrange [1:1.85]
set format x '\tiny $%0.1f$'
set xrange [5.5:6.35]
set y2tics 1,0.1,1.9 offset -0.7 mirror
set xtics 5.6,0.2,6.4
set title '\tiny (c) $C_{\mathrm{B_{1\mathrm g}}}$' offset -1,-6
plot "data/c11mc12.dat" using (log($1 + 1 / 0.00378087)):(log(lam($1) - 100 * $2)) with lines lw 3 lc rgb cc3, \
     7.3 - x