reset
#### Terminal settings
set term svg enhanced size 600,600 font "Helvetica, 16";set output 'fbrt.svg'
#### Mathematical definitions
## Fireball temperature
p0= 2.8526650801 ; q0= 69.18513566
p1=-1.0863702704 ; q1=115.00394909
p2= 4.7729619842 ; q2= 79.77893973
p3=10.026967289 ; q3= 27.50798341
p4= 7.1862268373 ; q4= 4.74096116
p5= 2.2724484357 ; q5= 0.32749001
p6= 0.27143255706
g1(x)=-1.103046*x+1.032988
g2(x)=(((((p6*x+p5)*x+p4)*x+p3)*x+p2)*x+p1)*x+p0
g3(x)=-0.547586*x+3.505928
h1(x)=((((q5*x+q4)*x+q3)*x+q2)*x+q1)*x+q0
g(x) = x<=-2.507 ? g1(x) : x<=-0.686 ? g2(x) : g3(x)
h(x) = x<=-2.22638 ? h1(x) : NaN
fg(x)= 10.**g(log10(x))
fh(x)= 10.**h(log10(x))
## Fireball diameter
rm4=13.; rm1=215.; re=236.; pr=log10(rm1/rm4)/3.; pk=-3.
d(x)=2*((rm4*(10000.*x)**pr)**pk+re**pk)**(1./pk)
#### Format definitions
set tmargin 0.2; set bmargin 2.8; set lmargin 9.6; set rmargin 1
lw1=3.; lw2=1.
unset key
set logscale xy
set xrange [1e-4:10]
set yrange [20:3e5]
set xtics ("10^{-4} " 0.0001,"" 0.0002,"" 0.0003,"" 0.0004,"" 0.0005,\
"" 0.0006,"" 0.0007,"" 0.0008,"" 0.0009,\
"10^{-3} " 0.001,"" 0.002,"" 0.003,"" 0.004,"" 0.005,\
"" 0.006,"" 0.007,"" 0.008,"" 0.009,\
0.01,"" 0.02,"" 0.03,"" 0.04,"" 0.05,"" 0.06,"" 0.07,"" 0.08,"" 0.09,\
0.1,"" 0.2,"" 0.3,"" 0.4,"" 0.5,"" 0.6,"" 0.7,"" 0.8,"" 0.9,\
1,"" 2,"" 3,"" 4,"" 5,"" 6,"" 7,"" 8,"" 9,10)
set ytics (10,"20" 20,"" 30,"" 40,"" 50,"" 60,"" 70,"" 80,"" 90,\
100,"" 200,"" 300,"" 400,"" 500,"" 600,"" 700,"" 800,"" 900,\
1000,"" 2000,"" 3000,"" 4000,"" 5000,"" 6000,"" 7000,"" 8000,"" 9000,\
"10000" 10000,"" 20000,"" 30000,"" 40000,"" 50000,\
"" 60000,"" 70000,"" 80000,"" 90000,\
"100000" 100000,"" 200000,"" 300000,"" 400000,"" 500000,\
"" 600000,"" 700000,"" 800000,"" 900000,\
"10^6 " 1e6)
set grid xtics ytics mxtics lt 1 lw lw2 lc rgb '#A0A0A0'
#### Labels
## Use oblique fonts if available
#set xlabel '{/Oblique t} / s'
set xlabel 't / s'
#set ylabel offset 1,0 '{/Oblique D} / m &{123456789012345} {/Oblique T}_{eff} / K'
set ylabel offset 0.5,0 'D / m &{123456789012345} T_{eff} / K'
#### Plotting
plot fg(x) w l lt 1 lw lw1 lc rgb '#700',\
fh(x) w l lt 1 lw lw1 lc rgb '#F80',\
d(x) w l lt 1 lw lw1 lc rgb '#00F'
reset