summaryrefslogtreecommitdiff
path: root/figs/fig-mag.gplot
diff options
context:
space:
mode:
authorJaron Kent-Dobias <jaron@kent-dobias.com>2017-06-14 14:10:53 -0400
committerJaron Kent-Dobias <jaron@kent-dobias.com>2017-06-14 14:10:53 -0400
commit1b990cf214d7e2ebf003db5f7c36565b36391633 (patch)
tree706eb2bb4a373ed4f8c4717aed9aeccda80c2711 /figs/fig-mag.gplot
parent3650d1101e092048cd4f29e444385e4a44509d63 (diff)
downloadpaper-1b990cf214d7e2ebf003db5f7c36565b36391633.tar.gz
paper-1b990cf214d7e2ebf003db5f7c36565b36391633.tar.bz2
paper-1b990cf214d7e2ebf003db5f7c36565b36391633.zip
changed the way data is processed for the plots
Diffstat (limited to 'figs/fig-mag.gplot')
-rw-r--r--figs/fig-mag.gplot15
1 files changed, 12 insertions, 3 deletions
diff --git a/figs/fig-mag.gplot b/figs/fig-mag.gplot
index 73ee9ec..df67fa9 100644
--- a/figs/fig-mag.gplot
+++ b/figs/fig-mag.gplot
@@ -2,11 +2,20 @@
set terminal pslatex rotate size 3.417,2.111
set multiplot
+Tc = 2 / log(1 + sqrt(2))
+Delta = 15. / 8
+beta = 1. / 8
+
+t(T) = abs((T - Tc) / T)
+h(T, H) = H / T
+X(T, H) = h(T, H) * t(T)**(-Delta)
+
B = 0.0861
A = 0.094
M0 = 1.21039
+
func = "figs/fig-mag_data-func.dat"
-num = "figs/fig-mag_data-num.dat"
+num = "data/data_square-2048.dat"
set samples 10000
set key off
@@ -16,7 +25,7 @@ set yrange [1.18:1.75]
set ylabel offset 1,0 '$Mt^{-\beta}$'
set xlabel '$ht^{-\Delta}$'
-plot num with yerrorbars pt 0, func using (10**$1 / B):(M0 + 10**($2) * A) with linespoints pt 0 lw 3
+plot num using (X($2, $3)):($6 * t($2)**(-beta)):($7 * t($2)**(-beta)):(hsv2rgb(20 * t($2), 1, 1)) with yerrorbars pt 0 lc rgb variable, func using (10**$1 / B):(M0 + 10**($2) * A) with linespoints pt 0 lw 2 lc rgb "black"
set size 0.6,0.6
set origin 0.38,0.2
@@ -27,7 +36,7 @@ set xlabel offset 0,0.5 '\footnotesize$ht^{-\Delta}$'
set xtics format '\footnotesize$10^{%g}$' -2,2,3
set ytics format '\footnotesize %g'
-plot num using (log10($1)):2:3 with yerrorbars pt 0, func using ($1-log10(B)):(M0 + A * 10**$2) with linespoints pt 0 lw 3
+plot num using (log10(X($2, $3))):($6 * t($2)**(-beta)):($7 * t($2)**(-beta)):(hsv2rgb(20 * t($2), 1, 1)) with yerrorbars pt 0 lc rgb variable, func using ($1-log10(B)):(M0 + A * 10**$2) with linespoints pt 0 lw 2 lc black