geo_MI_lock_fin4.gnu

reset
set xrange[0:200]
set yrange[-1:1]
set y2range[-400:0]
set xlabel "phi [deg] (MFE)"
set ylabel "error signal [mW]"
set y2label "feedback [deg]"
set mxtics 2
set mytics 2
set my2tics 2
set ytics .5
set y2tics 100
set zero 0.0
set title "geo_MI_lock                Mon Jan 15 13:45:42 2007", 2
set nolog y
set nolog y2
set term x11 
set size ratio .5 
set key below 
set grid xtics ytics 
y(gray)=abs(1-2*gray)
f7(x)=sqrt(y(x))
f5(x)=y(x)**3
f6(x)=y(x)**4
f15(x)=sin(2*pi*y(x))
f9(x)=sin(0.5*pi*y(x))
flip76(x)=0.5*((sgn(x-0.5)+1)*f7(x)+(1-sgn(x-0.5))*f6(x))
flip59(x)=0.5*((sgn(x-0.5)+1)*f5(x)+(1-sgn(x-0.5))*f9(x))
set palette model RGB functions flip76(gray), flip59(gray), f15(gray)
unset surface
unset hidden3d
set pm3d map
set size ratio 1
unset title
plot\
'geo_MI_lock.out' using ($1):($2*1000) axes x1y1 title " error signal [mW] " w l lt 1 lw 2, \
'geo_MI_lock.out' using ($1):($3) axes x1y2 title " feedback signal [deg] " w l lt 2 lw 2
set terminal pdf  
set key 100,-.7 box
set output "geo_MI_lock4.pdf"
replot