//////// load needed templates//////////// load_file("nrngui.hoc") load_file("../template/load_templates.hoc") objref econ show_errs=1 debug_lev=1 econ=new ExperimentControl(show_errs,debug_lev) econ.self_define(econ) econ.morphology_dir = "../morphology/n123" // set location for morphology files econ.add_lib_dir("Terrence","../lib") // set location for library files econ.generic_dir = "../experiment" // set location for cell-setup file v_init=-70 econ.xopen_geometry_dependent("cell") // load raw cell morphology econ.xopen_geometry_dependent("cell-analysis") // load user-defined semantics on morphology cell_analysis(econ) printf("Opening cell setup\n") // load cell-setup to econ.xopen_generic("cell-setup") // specify all mechanisms, printf("Opened. Setting up cell\n") // membrane properties etc cell_setup(econ) ///////////////////Load File////////////////////// load_file("sessio.ses") cvode_active(1) ////// Spike counter////////// objref apc, v1 proc insert_APC() { apc = new APCount(0.5) apc.thresh = $1 v1 = new Vector() apc.record(v1) } ////// IClamp /////// objectvar clamp proc insert_IClamp() { clamp = new IClamp(0.5) clamp.del = $3 clamp.dur = $1 clamp.amp = $2 } ////////////////////Main/////////////////// proc main(){ del=100 tstop=1000+del+200 tstopI=1000 in=0 th=-20 access soma insert_APC(th) for u= in,15{ current=u*0.1 access soma insert_IClamp(tstopI,current,del) //////////////run//////////////////////// finitialize(v_init) fcurrent() run() ///////////////////////////////////// // print the spike number printf("\n Current: %.2f nA | # Spikes: %d\n", current, apc.n) printf("___________________\n") } // for u } main() /******** end file ******/