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