#!/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