if __name__ == '__main__':
sim = Simulation('test_ss_soma')
soma = MyCompartment('soma', sim.model)
soma.length = 20e-6
soma.diameter = 2e-6 * 7.5
soma.setSpecificCm(9e-3)
soma.setSpecificRm(5.0)
soma.setSpecificRa(1.0)
soma.Em = -65e-3
soma.initVm = -65e-3
channel_lib = init_channels()
gk = {}
for channel, density in channel_density.items():
chan = channel_lib[channel]
new_chan = moose.HHChannel(chan, chan.name, soma)
chan = soma.insertChannel(new_chan, density)
chan.X = 0.0
gk[channel] = moose.Table(channel, sim.data)
gk[channel].stepMode = 3
gk[channel].connect('inputRequest', chan, 'Gk')
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:
import config
from cell import *
from capool import CaPool
from compartment import MyCompartment
from simulation import Simulation
sim = Simulation()
cell = moose.Cell('cell', sim.model)
soma = MyCompartment('soma', cell)
soma.diameter = 16e-6
soma.length = 15e-6
soma.Em = -70e-3
soma.initVm = -65e-3
soma.setSpecificRm(5.0)
soma.setSpecificRa(2.5)
soma.setSpecificCm(9e-3)
soma.insertChannel('NaF', specificGbar=1875.0, Ek=50e-3, shift=-3.5e-3)
soma.insertChannel('CaL', specificGbar=10.0, Ek=125e-3)
soma.insertCaPool(2600000.0, 100e-3)
vmTable = soma.insertRecorder('Vm', 'Vm', sim.data)
caTable = moose.Table('ca', sim.data)
caTable.stepMode = 3
soma.ca_pool.connect('Ca', caTable, 'inputRequest')
sim.schedule()
sim.run(200e-3)
sim.dump_data('data')
from pylab import *
mus_Ca = array(caTable)
mus_data = loadtxt('data/Vm.plot')
nrn_data = loadtxt('../nrn/mydata/Vm.plot')
nrn_Vm = nrn_data[:,1]
if __name__ == '__main__':
sim = Simulation('test_ss_soma')
soma = MyCompartment('soma', sim.model)
soma.length = 20e-6
soma.diameter = 2e-6 * 7.5
soma.setSpecificCm(9e-3)
soma.setSpecificRm(5.0)
soma.setSpecificRa(1.0)
soma.Em = -65e-3
soma.initVm = -65e-3
channel_lib = init_channels()
gk = {}
for channel, density in channel_density.items():
chan = channel_lib[channel]
new_chan = moose.HHChannel(chan, chan.name, soma)
chan = soma.insertChannel(new_chan, density)
chan.X = 0.0
gk[channel] = moose.Table(channel, sim.data)
gk[channel].stepMode = 3
gk[channel].connect('inputRequest', chan, 'Gk')
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:
