The literal answer to this question is probably not the best way to accomplish your goal. Please describe what you are trying to accomplish. As near as I can tell, your goal is to deliver a series of current pulses to some location in a model cell. Did I guess correctly?how to make Iclamps in shape of arrays
Read about the File and Vector classes in the Programmer's Reference. Relevant Vector class methods include scanf, scantil, vread, fread. Relevant File class methods include scanvar, scanstr, gets, vread. Also check the NEURON Forum's "Hot tips" section for relevant items.how neuron can read data from a file
So your goal is to deliver a series of current pulses to a location in a model cell. Go to the "Hot tips" area of the Forum and look for the item "How to change parameters during a simulation". That will show youmahsa wrote:yes you guessed correctly
Code: Select all
load_file("nrngui.hoc")
create soma
soma {
nseg = 1
diam = 18.8
L = 18.8
Ra = 123.0
insert hh
gnabar_hh=0.25
gl_hh = .0001666
el_hh = -60.0
}
soma psection()
objectvar stim1,stim2,stim3,stim4,stim5,stim6,stim7,stim8,stim9,stim10,stim11
stim1 = new IClamp(0.5)
stim1.del = 0
stim1.dur = 10
stim1.amp =-4
stim2 = new IClamp(0.5)
stim2.del = 10
stim2.dur = 10
stim2.amp = 5
stim3 = new IClamp(0.5)
stim3.del = 20
stim3.dur = 10
stim3.amp =-2
stim4 = new IClamp(0.5)
stim4.del = 30
stim4.dur = 10
stim4.amp = 5
stim5 = new IClamp(0.5)
stim5.del = 40
stim5.dur = 10
stim5.amp = 0
stim6 = new IClamp(0.5)
stim6.del = 50
stim6.dur = 10
stim6.amp = 5
stim7 = new IClamp(0.5)
stim7.del = 60
stim7.dur = 10
stim7.amp = 3
stim8 = new IClamp(0.5)
stim8.del = 70
stim8.dur = 30
stim8.amp = 5
stim9 = new IClamp(0.5)
stim9.del = 100
stim9.dur = 10
stim9.amp = 7
stim10 = new IClamp(0.5)
stim10.del = 110
stim10.dur = 10
stim10.amp = 5
stim11 = new IClamp(0.5)
stim11.del = 120
stim11.dur = 10
stim11.amp = 0
tstop = 500
objref rect, recv
rect = new Vector()
recv = new Vector()
recv.record(&soma.v(0.5))
rect.record(&t)
objref savv, savt
savv = new File()
savt = new File()
savv.wopen("Nv.dat")
savt.wopen("Nt.dat")
recv.printf(savv)
rect.printf(savt)
savv.close()
savt.close()
Code: Select all
proc myrun() { local ii
run()
print " "
for ii=0,rect.size()-1 print ii, " ", rect.x[ii], " ", recv.x[ii]
}Code: Select all
objref savv, savt
savv = new File()
savt = new File()
proc myrun() { local ii
run()
// print " "
// for ii=0,rect.size()-1 print ii, " ", rect.x[ii], " ", recv.x[ii]
// statements that open and write the output files
savv.wopen("Nv.dat")
savt.wopen("Nt.dat")
recv.printf(savv)
rect.printf(savt)
savv.close()
savt.close()
}Code: Select all
tvec.x[0] = 0 // start of step 0
ivec.x[0] = 0
tvec.x[1] = 10 // end of step 0
ivec.x[1] = 0
tvec.x[2] = 10 // start of step 1
ivec.x[2] = -4
tvec.x[3] = 20 // end of step 1
ivec.x[3] = -4
tvec.x[4] = 20 // start of step 2
ivec.x[4] = 5
tvec.x[5] = 30 // end of step 2
ivec.x[5] = 5
tvec.x[6] = 30 // start of remainder of simulation
ivec.x[6] = 0
tvec.x[7] = 1e9 // after returning to 0, stay there forever
ivec.x[7] = 0Code: Select all
objref ivec, tvec
ivec = new Vector()
tvec = new Vector()
// creates a step that applies no stimulus at all,
// so we can be sure the model is properly initialized
// just one argument: duration of this first step
proc step0() {
// throw out anything that might have been in ivec and tvec
ivec = new Vector(2) // this contains a pair of 0s, so we don't have to do anything to it
tvec = new Vector(2)
tvec.x[1] = $1
}
// arguments are step amplitude and duration
proc addstep() {
ivec.append($1)
ivec.append($1)
// last element of tvec contains the start time
tvec.append(tvec.x[tvec.size()-1])
// step ends $2 later
tvec.append(tvec.x[tvec.size()-1] + $2)
}
proc set_up_stim() {
step0 // initialize ivec and tvec, and fill them with the "no stimulus" step
// add three steps:
addstep(-4, 10) // -4 nA x 10 ms
addstep(5, 10) // 5 nA x 10 ms
addstep(0, 1e9) // 0 nA forever
}