// This script shows NMDA oscillations when
// Na channels are blocked
// Part of the resulting trace is shown in
// Fig 8a.
// This oscillation is, however, very
// sensitive using these particular parameter
// values.
// The oscillations can be more stable by
// various means, e.g. changing the KCaNMDA
// activation function.
float PI = 3.141592654
int i
float dt = 1e-4
float sim_time = 10
setclock 0 {dt}
setclock 1 1e-4
setmethod 11
// setmethod 0
float NMDA=0.7
float AMPA=0.0
include library
include finish
make_cell cell.p neuron
include xout_fig8
create spikegen /neuron/soma/spike
setfield /neuron/soma/spike thresh -0.035 abs_refract 0.010 output_amp 1
create spikehistory spike.history
setfield spike.history ident_toggle 0 \ // index specification
filename "spikes.dat" \
initialize 1 leave_open 1 flush 1
addmsg /neuron/soma/spike spike.history SPIKESAVE
reset
str comp
foreach comp ({el /neuron/##[OBJECT=compartment]})
setfield {comp} Em -0.070
setfield {comp} initVm -0.070
setfield {comp}/fshNa Gbar 0
setfield {comp}/fshNa-is Gbar 0
setfield {comp}/fshNa-dend Gbar 0
end
// setfield neuron chanmode 0
// if hsolve is used
call /neuron SETUP
reset
step 10 -t
|