#!/bin/sh
# Plotting trajectories for all parameter sets, locked to the stimulus
# TC (VPM and POm) and RE populations
pra=tlp
prg=plp
fla=d1
flb=d2
grfl=ptrint.gr
tlp.ex $fla
plp.ex $fla
stm.ex $fla
cat > ptrint.xx.fsed <<EOF
s/scan=u/scan=e/
s/INPUT ff 2.0 5.0 8.0 11.0/INPUT ff parmin=1.0 parmax=11.5 npar=21/
EOF
sed -f ptrint.xx.fsed $pra.n.$fla > $pra.n.$flb
tlp.ex $flb
plp.ex $flb
awk '{if ($2 == "1") print $1, 1000.0 * $7}' $prg.res.$flb > ptrint.xx.ult.in
awk '{if ($2 == "1") print $1, 1000.0 * $6}' $prg.res.$flb > ptrint.xx.ult.t5
awk '{if ($2 == "2") print $1, 1000.0 * $7}' $prg.res.$flb > ptrint.xx.upt.in
awk '{if ($2 == "2") print $1, 1000.0 * $6}' $prg.res.$flb > ptrint.xx.upt.t5
awk '{if ($1 != "#") print $1, $2}' delay_exp > delay.exp.vpm.xx
awk '{if ($1 != "#") print $1, $3}' delay_exp > delay.exp.pom.xx
awk '{print $1, $2 * 40.0}' intg.res.a3 > intg.res.a3.xx
awk '{print $1, $2 * 40.0}' intg.res.a4 > intg.res.a4.xx
intnor.ex ptrint ult 40.0
intnor.ex ptrint upt 40.0
xmgrace \
-graph 0 plp.xx.$fla.1.ult plp.xx.$fla.2.ult \
plp.xx.$fla.3.ult plp.xx.$fla.4.ult \
-graph 1 plp.xx.$fla.1.upt plp.xx.$fla.2.upt \
plp.xx.$fla.3.upt plp.xx.$fla.4.upt \
-graph 2 stm.flt.$fla \
-graph 3 stm.fpt.$fla \
-graph 4 intn.res.ult intg.res.a3.xx \
ptrint.xx.ult.t5 delay.exp.vpm.xx \
-graph 5 intn.res.upt intg.res.a4.xx \
ptrint.xx.upt.t5 delay.exp.pom.xx \
-graph 6 plp.xx.$fla.1.ura plp.xx.$fla.2.ura \
plp.xx.$fla.3.ura plp.xx.$fla.4.ura \
-graph 7 plp.xx.$fla.1.urb plp.xx.$fla.2.urb \
plp.xx.$fla.3.urb plp.xx.$fla.4.urb \
-hdevice EPS -p $grfl -printfile ptrint.$fla.eps
/bin/rm $pra.out.$fla $pra.tmp.$fla $pra.col.$fla
/bin/rm $pra.out.$flb $pra.tmp.$flb $pra.col.$flb
/bin/rm plp.xx.$fla.*.*
/bin/rm stm.flt.$fla stm.fpt.$fla stm.ffr.$fla
/bin/rm ptrint.xx.fsed
/bin/rm plp.xx.$flb.*.*
/bin/rm ptrint.xx.ult.in ptrint.xx.ult.t5 ptrint.xx.upt.in ptrint.xx.upt.t5
/bin/rm delay.exp.vpm.xx delay.exp.pom.xx
/bin/rm intg.res.a3.xx intg.res.a4.xx
/bin/rm intn.out.ult intn.res.ult intn.out.upt intn.res.upt
| 
