// $Id: batch.hoc,v 1.124 2012/01/06 02:27:50 samn Exp $ // run this one with : // pbsb -s -ng -n 8 -i batch.hoc // @tauahpRS = (12.5,25,50,100,200,400,800,1600) /// //* setup xopen("/usr/site/nrniv/simctrl/hoc/setup.hoc") xopen("/usr/site/nrniv/simctrl/hoc/nrnoc.hoc") load_file("init.hoc") //* batch/logging related batch_flag=1 strdef logstr objref bfile bfile=new File() proc logit () { bfile.aopen("/u/samn/intfzip/data/12jan05.02_batch_tauahpRS_TEST_A_.txt") {bfile.printf($s1) bfile.printf("\n") bfile.close()} } //* load correct version of sim strdef strrcs strrcs="nqsnet.hoc,65,network.hoc,204,params.hoc,273,run.hoc,66,nload.hoc,202,basestdp.hoc,222" rcsopen(strrcs) sprint(strv,"12jan05.02_tauahpRS_%g_",tauahpRS) print "sim version: " , strv logit(strv) // write info to log //* post-sim-load modifications //* load data from sim that was already run // mytstop=tstop=htmax=10e3 // set to correct duration of sim that was run // loadminrundat(strv) //myrd(strv) //* run & save output //** main run/save // minrunsv(0,0) myrunsv(strv) // myrd(strv) // rcsopen("nqpplug.hoc",6) // loads & saves mkdrs results // savenqspec(strv) // {wgnq=mkwgnq(col) sprint(tstr,"/u/samn/intfzip/data/%s_wgnq.nqs",strv) wgnq.sv(tstr)} //* finish logging current run print strv // {pravgrates(tstop,1,tstr) logstr=strv strcat(logstr,tstr) logit(logstr)}