make the tics logscale

1 files changed, 19 insertions, 14 deletions

diff --git a/figs/fig-susmag.gplot b/figs/fig-susmag.gplot
index 39a3fdf..da4f8e6 100644
--- a/figs/fig-susmag.gplot
+++ b/figs/fig-susmag.gplot
@@ -34,11 +34,13 @@ set key off
 
 set size 1,1 - 0.52
 set origin 0,0.52
-set xrange [-20:1200]
-set yrange [-0.5:13]
+set xrange [-40:1200]
+set yrange [-0.4:14.4]
 set ylabel offset 1,0 '$\chi|t|^{\gamma}\times 10^3$'
 set lmargin 6
 set xtics format ''
+set mxtics 5
+set mytics 5
 set bmargin 0.2
 
 plot num using (X($2, $3)):(10**3 * $10 * t($2)**gamma):(10**3 * $11 * t($2)**gamma):(hsv2rgb(20 * t($2), 1, 1)) with yerrorbars pt 0 lc rgb variable, susfunc using (10**$1 / B):(10**(3+$2) * A * B**2) with linespoints pt 0 lw 2 lc rgb "black", susfunc using (10**$1 / B):(10**(3+$2) * A * B**2 + 10**3 * (poly(C1, lamb, 1, 10**$1) + poly(C2, lamb, 2, 10**$1) + poly(C3, lamb, 3, 10**$1))) with lines dt 2 lw 2 lc black
@@ -47,35 +49,38 @@ set bmargin -1
 set tmargin 0.2
 set size 1,0.52
 set origin 0,0
-set xrange [-20:1200]
-set yrange [1.18:1.75]
+set yrange [1.18:1.74]
 set ylabel offset 1,0 '$M|t|^{-\beta}$'
 set xlabel '$h|t|^{-\beta\delta}$'
 set xtics format '%g'
 
 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, magfunc using (10**$1 / B):(M0 + 10**($2) * Tc * A * B) with linespoints pt 0 lw 2 lc black, magfunc using (10**$1 / B):(M0 + 10**($2) * Tc * A * B + (polyint(Tc * C1, lamb, 1, 10**$1) + polyint(Tc * C2, lamb, 2, 10**$1) + polyint(Tc * C3, lamb, 3, 10**$1)) / B) smooth csplines with lines dt 2 lw 2 lc black
 
+set logscale xy
 set tmargin -1
 set lmargin -1
 set size 0.65,0.325
 set origin 0.31,0.5 + 0.29 / 2
-set xrange [-3:3]
-set yrange [-4.5:-1.5]
+set xrange [0.0015:1900]
+set yrange [0.00002:0.02]
 set ylabel offset 2.5,0 '\footnotesize$\chi|t|^\gamma$'
 set xlabel offset 0,0.5 '\footnotesize$h|t|^{-\beta\delta}$'
-set xtics format '\footnotesize$10^{%g}$' -2,2,3
-set ytics format '\footnotesize$10^{%g}$' -4,1,-2
+set mxtics 5
+set xtics format '' -2,10,1000
+set xtics add ('$\footnotesize10^{-2}$' 10**(-2), "" 0.1, '$\footnotesize10^0$' 1, "" 10, '$\footnotesize10^2$' 100, "" 1000)
+set mytics 5
+set ytics format '\footnotesize$10^{%T}$' 0.00001,10,0.01
 
-plot num using (log10(X($2, $3))):(log10($10 * t($2)**gamma)):($11 / $10):(hsv2rgb(20 * t($2), 1, 1)) with yerrorbars pt 0 lc rgb variable, susfunc using ($1-log10(B)):($2+log10(A * B**2)) with linespoints pt 0 lw 2 lc rgb "black", susfunc using ($1-log10(B)):(log10(10**$2 * A * B**2 + poly(C1, lamb, 1, 10**$1) + poly(C2, lamb, 2, 10**$1) + poly(C3, lamb, 3, 10**$1))) with lines dt 2 lw 2 lc black
+plot num using (X($2, $3)):($10 * t($2)**gamma):($11 * t($2)**gamma):(hsv2rgb(20 * t($2), 1, 1)) with yerrorbars pt 0 lc rgb variable, susfunc using (10**$1 / B):(10**$2 * A * B**2) with linespoints pt 0 lw 2 lc rgb "black", susfunc using (10**$1 / B):(10**$2 * A * B**2 + poly(C1, lamb, 1, 10**$1) + poly(C2, lamb, 2, 10**$1) + poly(C3, lamb, 3, 10**$1)) with lines dt 2 lw 2 lc black
 
+unset logscale xy
+set logscale x
 set size 0.65,0.325
 set origin 0.31,0.23 / 2
-set xrange [-3:3]
 set yrange [1.18:1.75]
 set ylabel offset 4,0 '\footnotesize$M|t|^{-\beta}$'
-set xlabel offset 0,0.5 '\footnotesize$h|t|^{-\beta\delta}$'
-set xtics format '\footnotesize$10^{%g}$' -2,2,3
-set ytics format '\footnotesize {%g}' 1.2,.2,1.6
+set ytics format '\footnotesize {%g}' 1.2,0.2,1.8
+set mytics 5
 
-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, magfunc using ($1-log10(B)):(M0 + Tc * A * B * 10**$2) with linespoints pt 0 lw 2 lc black, magfunc using ($1 - log10(B)):(M0 + 10**($2) * Tc * A * B + (polyint(Tc * C1, lamb, 1, 10**$1) + polyint(Tc * C2, lamb, 2, 10**$1) + polyint(Tc * C3, lamb, 3, 10**$1)) / B) with lines dt 2 lw 2 lc black
+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, magfunc using (10**$1 / B):(M0 + Tc * A * B * 10**$2) with linespoints pt 0 lw 2 lc black, magfunc using (10**$1 / B):(M0 + 10**($2) * Tc * A * B + (polyint(Tc * C1, lamb, 1, 10**$1) + polyint(Tc * C2, lamb, 2, 10**$1) + polyint(Tc * C3, lamb, 3, 10**$1)) / B) with lines dt 2 lw 2 lc black