Code: Select all
{
load_file("nrngui.hoc")
}
create soma
soma{
nseg =1
diam = 18.8
L = 18.8
}
Rm = 2.62e-5
global_ra = 100
Vleak = -50
soma {
insert pas g_pas=1/(Rm) Ra=global_ra e_pas=Vleak
}
objref i_i
soma {
i_i = new SEClamp(0.5)
}
tstop = 1000
objref tvec
vecsize = tstop
tvec = new Vector( vecsize , 0 )
for i=0, vecsize-2 {
tvec.x[i+1] = tvec.x[i] + 1
}
objref condvec_i
vecsize = tvec.size()
condvec_i = new Vector( vecsize , 0 )
for i=0, vecsize-1 {
t = tvec.x[i]
condvec_i.x[i] = 1e2
}
i_i.amp1 = -70
i_i.dur1 = 100
condvec_i.play(&i_i.rs, tvec, 1)
objref vm_trace, time_trace
vm_trace = new Vector()
time_trace = new Vector()
objref se_amp, se_rs,se_vc,se_i
se_amp = new Vector()
se_rs = new Vector()
se_vc = new Vector()
se_i = new Vector()
objref f1, ftime1,f_seamp,f_sers,f_sevc,f_sei
f1 = new File()
f1.wopen("./results/vm_trace_sample_soma.dat")
ftime1 = new File()
ftime1.wopen("./results/time_trace_sample_soma.dat")
f_seamp = new File()
f_seamp.wopen("./results/seamp_trace_sample_soma.dat")
f_sers = new File()
f_sers.wopen("./results/sers_trace_sample_soma.dat")
f_sevc = new File()
f_sevc.wopen("./results/sevc_trace_sample_soma.dat")
f_sei = new File()
f_sei.wopen("./results/sei_trace_sample_soma.dat")
vm_trace.record( &soma.v(0.1) )
time_trace.record( &t )
se_amp.record( &i_i.amp1 )
se_rs.record( &i_i.rs )
se_vc.record( &i_i.vc )
se_i.record( &i_i.i )
v_init = -65
{
tstop = 1000
dt = 0.1
steps_per_ms = 10
init()
run()
}
vm_trace.printf(f1)
time_trace.printf(ftime1)
se_amp.printf(f_seamp)
se_rs.printf(f_sers)
se_vc.printf(f_sevc)
se_i.printf(f_sei)
ftime1.close()
f1.close()
f_seamp.close()
f_sers.close()
f_sei.close()
f_sevc.close()
quit()