gk[channel].connect('inputRequest', chan, 'Gk')
chandict[channel] = chan
print chan.name, chan.Gbar
if channel.startswith('K'):
chan.Ek = EK
elif channel.startswith('Na'):
chan.Ek = ENa
elif channel.startswith('Ca'):
chan.Ek = ECa
elif channel.startswith('AR'):
chan.Ek = EAR
# chan.X = 0.25
else:
print 'Error: unknown channel', channel
cad = moose.CaConc('CaPool', soma)
cad.B = 2.6e7 / soma.sarea()
cad.tau = 50e-3
cad.connect('concSrc', chandict['KAHP_SLOWER'], 'concen')
cad.connect('concSrc', chandict['KC_FAST'], 'concen')
chandict['CaL'].connect('IkSrc', cad, 'current')
ca_table = moose.Table('Ca', sim.data)
ca_table.stepMode = 3
ca_table.connect('inputRequest', cad, 'Ca')
vm_table = soma.insertRecorder('Vm', 'Vm', sim.data)
soma.insertPulseGen('pulsegen',
sim.model,
firstLevel=3e-10,
firstDelay=100e-3,
firstWidth=50e-3)
print 'Rm:', soma.Rm
print 'Cm:', soma.Cm,
firstLevel=3e-10,
firstDelay=20.0e-3,
firstWidth=1e3)
print 'Rm:', soma.Rm
print 'Cm:', soma.Cm,
print 'Em:', soma.Em
gk_naf2_table = moose.Table('Gk_NaF2', sim.data)
gk_naf2_table.stepMode = 3
chan = moose.HHChannel(soma.path + '/NaF2')
print chan.Gbar, chan.Ek
chan.connect('Gk', gk_naf2_table, 'inputRequest')
sim.schedule()
config.context.useClock(0, sim.model.path + '/##')
soma.useClock(1, 'init')
for channel in soma.channels:
channel.useClock(0)
sim.run(50e-3)
for i in range(len(gk_naf2_table)):
gk_naf2_table[i] = gk_naf2_table[i] / soma.sarea()
sim.dump_data('data')
pylab.subplot(211)
pylab.plot(vm_table)
pylab.subplot(212)
for key, value in gk.items():
pylab.plot(value, label=key)
pylab.legend()
pylab.show()
#
# test_ss_soma.py ends here
soma.insertPulseGen('pulsegen', sim.model, firstLevel=3e-10, firstDelay=20.0e-3, firstWidth=1e3)
print 'Rm:', soma.Rm
print 'Cm:', soma.Cm,
print 'Em:', soma.Em
gk_naf2_table = moose.Table('Gk_NaF2', sim.data)
gk_naf2_table.stepMode = 3
chan = moose.HHChannel(soma.path + '/NaF2')
print chan.Gbar, chan.Ek
chan.connect('Gk', gk_naf2_table, 'inputRequest')
sim.schedule()
config.context.useClock(0, sim.model.path + '/##')
soma.useClock(1, 'init')
for channel in soma.channels:
channel.useClock(0)
sim.run(50e-3)
for i in range(len(gk_naf2_table)):
gk_naf2_table[i] = gk_naf2_table[i] / soma.sarea()
sim.dump_data('data')
pylab.subplot(211)
pylab.plot(vm_table)
pylab.subplot(212)
for key, value in gk.items():
pylab.plot(value, label=key)
pylab.legend()
pylab.show()
#
# test_ss_soma.py ends here
