set encoding iso_8859_1
set xla 'Air wavelength [Å]'
set yla 'Wavelength difference O-C [Å]'
set term post eps enh color solid 20
# 
# gamma prime
#
unset label
set label 99 "O-C wavelength differences for TiO {/Symbol g \242} system\nO: Davis et al. 1986 (from Plez)\nC: Plez 2008" at graph 0.03,0.95 front
set key right
#set xra [5800:7000]
#set yra [-6:2]
#
set output 'O-C-H.ps'
set label  '0-0' at 6250, 0.2 front
set label  '1-0' at 5900, 0.3 front
set label  '0-1' at 6620, 0.4 front
set label  '1' at 5820, 0.10 front
set label  '1' at 6130, 0.10 front
set label  '1' at 6530, 0.20 front
set label  '1' at 6050,-1.20 front
set label  '1' at 6430,-2.50 front
set label  '1' at 6900,-4.00 front
set label  '2' at 5880,-0.07 front
set label  '2' at 6190,-0.07 front
set label  '2' at 6600,-0.07 front
set label  '3' at 6050,-0.50 front
set label  '3' at 6420,-1.10 front
set label  '3' at 6780,-1.00 front
set label  '3' at 6050,-2.50 front
set label  '3' at 6400,-5.00 front
set label  '3' at 6770,-5.00 front
plot \
"< awk '$7==\"H\" && $6~\"P\" && ($2==0 && $3==0)' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.5, \
"< awk '$7==\"H\" && $6~\"Q\" && ($2==0 && $3==0)' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.5, \
"< awk '$7==\"H\" && $6~\"R\" && ($2==0 && $3==0)' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.5, \
"< awk '$7==\"H\" && $6~\"P\" && ($2==1 && $3==0)' TiO-Phillips-O-C.dat" u 8:9 t ''  w p pt 7 lt 7 ps 0.5, \
"< awk '$7==\"H\" && $6~\"Q\" && ($2==1 && $3==0)' TiO-Phillips-O-C.dat" u 8:9 t ''  w p pt 7 lt 1 ps 0.5, \
"< awk '$7==\"H\" && $6~\"R\" && ($2==1 && $3==0)' TiO-Phillips-O-C.dat" u 8:9 t ''  w p pt 7 lt 3 ps 0.5, \
"< awk '$7==\"H\" && $6~\"P\" && ($2==0 && $3==1)' TiO-Phillips-O-C.dat" u 8:9 t ''  w p pt 7 lt 7 ps 0.5, \
"< awk '$7==\"H\" && $6~\"Q\" && ($2==0 && $3==1)' TiO-Phillips-O-C.dat" u 8:9 t ''  w p pt 7 lt 1 ps 0.5, \
"< awk '$7==\"H\" && $6~\"R\" && ($2==0 && $3==1)' TiO-Phillips-O-C.dat" u 8:9 t ''  w p pt 7 lt 3 ps 0.5
# 
# gamma
#
unset label
set label 99 "O-C wavelength differences for TiO {/Symbol g} system\nO: Davis et al. 1986 (from Plez)\nC: Plez 2008" at graph 0.03,0.15 front
set key left
set xra [6250:8800]
set yra [-0.7:0.7]
#
set output 'O-C-G0.ps'
set label 1 '5-3' at 6500,-0.2 front
set label 2 '0-0' at 7000,-0.2 front
set label 3 '0-1' at 7600,-0.2 front
set label 4 '0-2' at 8100,-0.2 front
set label 11 '1' at 7050,-0.13 front
set label 12 '2' at 7000,-0.05 front
set label 13 '3' at 6950, 0.05 front
set label 21 '1' at 7500, 0.50 front
set label 22 '2' at 7500, 0.18 front
set label 23 '3' at 7450,-0.23 front
plot \
"< awk '$7==\"G\" && $6~\"P\" && ($2==0 || $2==5)' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"G\" && $6~\"Q\" && ($2==0 || $2==5)' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"G\" && $6~\"R\" && ($2==0 || $2==5)' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.3
#
set output 'O-C-G1.ps'
set label 1 '1-0' at 6600,-0.15 front
set label 2 '1-1' at 7100,-0.15 front
set label 3 '1-2' at 7600,-0.15 front
set label 4 '1-3' at 8200,-0.15 front
set label 11 '1' at 7150,-0.10 front
set label 12 '2' at 7100,-0.02 front
set label 13 '3' at 7050, 0.07 front
set label 21 '1' at 7600, 0.50 front
set label 22 '2' at 7580, 0.18 front
set label 23 '3' at 7450,-0.35 front
plot \
"< awk '$7==\"G\" && $6~\"P\" && ($2==1)' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"G\" && $6~\"Q\" && ($2==1)' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"G\" && $6~\"R\" && ($2==1)' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.3
#
set output 'O-C-G2.ps'
set label 1 '2-0' at 6400,-0.15 front
set label 2 '2-1' at 6700,-0.15 front
set label 3 '2-3' at 7600,-0.15 front
set label 4 '2-4' at 8250,-0.15 front
unset label 11
unset label 12
unset label 13
unset label 21
unset label 22
unset label 23
plot \
"< awk '$7==\"G\" && $6~\"P\" && ($2==2)' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"G\" && $6~\"Q\" && ($2==2)' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"G\" && $6~\"R\" && ($2==2)' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.3
#
set output 'O-C-G3.ps'
set label 1 '3-1' at 6400,-0.20 front
set label 2 '3-2' at 6750,-0.20 front
set label 3 '3-4' at 7700,-0.20 front
set label 4 '3-5' at 8350,-0.20 front
plot \
"< awk '$7==\"G\" && $6~\"P\" && ($2==3)' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"G\" && $6~\"Q\" && ($2==3)' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"G\" && $6~\"R\" && ($2==3)' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.3
#
set output 'O-C-G4.ps'
set label 1 '4-2' at 6400,-0.25 front
set label 2 '4-3' at 6800,-0.25 front
set label 3 '4-5' at 7800,-0.25 front
unset label 4
plot \
"< awk '$7==\"G\" && $6~\"P\" && ($2==4)' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"G\" && $6~\"Q\" && ($2==4)' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"G\" && $6~\"R\" && ($2==4)' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.3
# 
# beta
#
set output 'O-C-B.ps'
unset label
set label "O-C wavelength differences for TiO {/Symbol b} system\nO: Davis et al. 1986 (from Plez)\nC: Plez 2008" at graph 0.03,0.95
set key bottom
set xra [5580:5740]
set yra [-2:2]
set label '0-0' at 5595,-0.2 front
set label '1-1' at 5625,-0.7 front
set label '2-2' at 5660,-1.2 front
plot \
"< awk '$7==\"B\" && $6==\"P\"' TiO-Phillips-O-C.dat" u 8:9 t 'P' w p pt 7 lt 7 ps 0.5, \
"< awk '$7==\"B\" && $6==\"Q\"' TiO-Phillips-O-C.dat" u 8:9 t 'Q' w p pt 7 lt 1 ps 0.5, \
"< awk '$7==\"B\" && $6==\"R\"' TiO-Phillips-O-C.dat" u 8:9 t 'R' w p pt 7 lt 3 ps 0.5
#
# alpha
#
set output 'O-C-A.ps'
unset label
# title
set label "O-C wavelength differences for TiO {/Symbol a} system\nO: Davis et al. 1986 (from Plez and Kurucz)\nC: Plez 2008" at graph 0.30,0.95
#set key horizontal samplen -2 graph 0.95,0.75
unset key
# key
set label 'P1' tc lt 7 point pt 7 lt 7 offset 1,-0.3 at graph 0.30,0.75
set label 'P2' tc lt 3 point pt 7 lt 3 offset 1,-0.3 at graph 0.40,0.75
set label 'P3' tc lt 9 point pt 7 lt 9 offset 1,-0.3 at graph 0.50,0.75
set label 'R1' tc lt 1 point pt 7 lt 1 offset 1,-0.3 at graph 0.60,0.75
set label 'R2' tc lt 5 point pt 7 lt 5 offset 1,-0.3 at graph 0.70,0.75
set label 'R3' tc lt 8 point pt 7 lt 8 offset 1,-0.3 at graph 0.80,0.75
set xra [4350:6400]
#set xra [4460:4600]
set yra [-3.0:4.0]
#set yra [-1.0:1.0]
# Deltav=-3
set label '0-3' at 6100,-0.7 front
# Deltav=-2
set label '0-2' at 5750,-0.5 front
set label '1-3' at 5800,-0.7 front
set label '2-4' at 5850, 0.8 front
set label '3-5' at 5900, 0.6 front
# Deltav=-1
set label '0-1' at 5450,-0.5 front
set label '1-2' at 5500,-0.7 front
set label '2-3' at 5550, 0.7 front
# Deltav=0
set label '0-0' at 5150, 0.2 front
# Deltav=1
set label '1-0' at 4950,-0.5 front
set label '2-1' at 5000, 0.6 front
# Deltav=2
set label '2-0' at 4750, 0.5 front
set label '3-1' at 4800, 0.3 front
set label '4-2' at 4850, 0.1 front
# Deltav=3
set label '3-0' at 4580, 0.7 front
set label '4-1' at 4590, 0.45 front
set label '5-2' at 4600, 0.2 front
set label '6-3' at 4710, 2.0 front
# Deltav=4
set label '4-0' at 4390, 0.60 front
set label '5-1' at 4400, 0.30 front
set label '6-2' at 4480, 2.1 front
set label '7-3' at 4540, 1.8 front
#
plot \
"< awk '$7==\"A\" && $6==\"P1\"' TiO-Phillips-O-C.dat" u 8:9 t 'P1' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"A\" && $6==\"P2\"' TiO-Phillips-O-C.dat" u 8:9 t 'P2' w p pt 7 lt 3 ps 0.3, \
"< awk '$7==\"A\" && $6==\"P3\"' TiO-Phillips-O-C.dat" u 8:9 t 'P3' w p pt 7 lt 9 ps 0.3, \
"< awk '$7==\"A\" && $6==\"R1\"' TiO-Phillips-O-C.dat" u 8:9 t 'R1' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"A\" && $6==\"R2\"' TiO-Phillips-O-C.dat" u 8:9 t 'R2' w p pt 7 lt 5 ps 0.3, \
"< awk '$7==\"A\" && $6==\"R3\"' TiO-Phillips-O-C.dat" u 8:9 t 'R3' w p pt 7 lt 8 ps 0.3
#quit
#
# all systems
#
set output 'O-C-all.ps'
unset label
# title
set label "O-C wavelength differences for TiO\nO: Davis et al. 1986 (from Plez and Kurucz)\nC: Plez 2008" at graph 0.30,0.95
#set key bottom
unset key
# key
set label '{/Symbol a}' tc lt 1 point pt 7 lt 1 offset 1,-0.3 at graph 0.80,0.40
set label '{/Symbol b}' tc lt 3 point pt 7 lt 3 offset 1,-0.3 at graph 0.80,0.30
set label '{/Symbol g}' tc lt 7 point pt 7 lt 7 offset 1,-0.3 at graph 0.80,0.20
set label '{/Symbol g \242}' tc lt 9 point pt 7 lt 9 offset 1,-0.3 at graph 0.80,0.10
set xra [4300:9000]
set yra [-5.0:4.0]
plot \
"< awk '$7==\"A\"' TiO-Phillips-O-C.dat" u 8:9 t '{/Symbol a}' w p pt 7 lt 1 ps 0.3, \
"< awk '$7==\"B\"' TiO-Phillips-O-C.dat" u 8:9 t '{/Symbol b}' w p pt 7 lt 3 ps 0.3, \
"< awk '$7==\"G\"' TiO-Phillips-O-C.dat" u 8:9 t '{/Symbol g}' w p pt 7 lt 7 ps 0.3, \
"< awk '$7==\"H\"' TiO-Phillips-O-C.dat" u 8:9 t '{/Symbol g \242}' w p pt 7 lt 9 ps 0.3

# && ($2==4 && $3==2)    && $3==0   || $2==5
# && ($2==4 && $3==2)    && $3==0   || $2==5
# && ($2==4 && $3==2)    && $3==0   || $2==5
# && ($2==4 && $3==2)
# && ($2==4 && $3==2)
# && ($2==4 && $3==2)
quit