set term postscript eps enhanced "Times-Roman" 26 color
set out 'figSI07b.eps'
set size 1,1
set log y
set format y "10^{%T}"
set pointsize 2.6
set key left
t_0=1e3
set xlabel 'Z (=P/{/Symbol r})'
set ylabel '{/Symbol t}_{/Symbol a}/{/Symbol t}_0'
set size square
#t0_vft=exp(-4.07)
#A0_vft=2308.0
#Z0_vft=37.7
#t2_vft=exp(-2.11)
#A2_vft=3008.0
#Z2_vft=49.8

t30=1800
t32=3000
Z0_par=19.84
t0_par=8.25e4/t30
A0_par=0.22

A0_par=0.22
t0_par=8.25e4/t30
Z0_par=19.84
t2_par=2.37e4/t32
A2_par=0.039
Z2_par=21.55

#t0_par=12e4/t30
#A0_par=0.34
#Z0_par=19.81
#A2_par=0.10
#Z2_par=21.6
#t2_par=0.92e4/t30
#unset key
set ytics 1e-5
set xtics 5
p [18:35][:1e24]"tau3D_N2000_M7_Rmax1.614Rmin0.724SWP0.00.dat" u 5:($2/t30) pt 6 lw 3 lc rgb '#A0522D' t'std',\
"tau3D_N2000_M7_Rmax1.614Rmin0.724SWP0.20.dat" u 5:($2/t32) pt 8 lc 3 lw 3 t'SWAP',\
t0_par*exp(A0_par*(x-Z0_par)**2.0) lt 1 lw 5 lc rgb '#A0522D' t'',\
t2_par*exp(A2_par*(x-Z2_par)**2.0) lt 1 lw 5 lc 3 t'',\
"line3d.dat" w l lc -1 lt 1 lw 2 t''
#t2_vft*exp(A2_vft/(Z2_vft-x)**2.0) lt 1 lw 5 lc 3 t''