set term pdf enhanced color font "helvetica,20" linewidth 1 size 10,6
set output "tau_O2O2_CIA.pdf"

set view map
set size square
set samples 50, 50
set isosamples 400,400
unset surface
set style data lines

set contour base
set cntrparam bspline
set cntrparam levels incremental 1e-112,1e-112,1e-111
#set cntrparam levels auto
set table 'cont.dat'
splot "./tau_O2O2_cia.txt" using 1:($2/1000):3
unset table
load '<./contour_norotate.sh cont.dat 0 15 0 > ./labels.txt; ./contour_norotate.sh cont.dat 1 15 0 > ./contours.txt'

reset
unset key
unset surface
set contour base
set title "N2-N2 CIA opt. thickness" 
set xlabel "{/Symbol n} (cm^{-1})"
set ylabel "Altitude (km)"
#load "./labels.txt"
p './contours.txt' w l lt -1 lw 0.5

quit
