#!icc -O3 -xHost -I /usr/include/x86_64-linux-gnu/c++/4.7 NaP_KA.cc -o bart
!g++ -O2 NaP-EL-3cell.cc -o bart
start=20
fin=80
Y1=0
Y2=0
aa1=0
aa2=12
x1=0
x2=1.2
set print "tmp1.do"
print sprintf("./bart -dt .1 -f 100 -l %g",fin)
set print
!chmod +x tmp1.do
!./tmp1.do >dat

set term post eps enhanced solid
set out "tmp.eps"
set grid
unset key
set lmargin 10
set multi lay 6,1
set ylabel "TONIC"
set ytics 1
plot [start:fin*100][0:1]"dat" u 0:($1<-50?0:($1>0?1:($1+50)/50)) w l
set ylabel "HIGH\nFREQ"
plot [start:fin*100][0:1]"dat" u 0:($2<-50?0:($2>0?1:($2+50)/50)) w l
set ylabel "LOW\nFREQ"
plot [start:fin*100][0:1]"dat" u 0:($3<-50?0:($3>0?1:($3+50)/50)) w l
set ylabel "SILENT"
plot [start:fin*100][0:1]"dat" u 0:($4<-50?0:($4>0?1:($4+50)/50)) w l
set ylabel "SUMMATION"
set bmargin at screen 0.1
set ytics 100
set xlabel "TIME"
#plot [][0:1]"dat" u 0:(($1+$2+$3+$4)<-200?0:(($1+$2+$3+$4)>0?1:($1+$2+$3+$4+200)/200)) w l
plot [start:fin*100][-250:-50]"dat" u 0:($1+$2+$3+$4) w l lc 1 lw 3
#plot [2000:fin*100][0:1]"dat" u 0:($5<-50?0:($5>0?1:($5+50)/50)) w l lc 1 lw 3
#plot [2000:fin*100][]"dat" u 0:5 w l lc 1 lw 3

unset multi