Code: Select all
for nd in self.ip3.nodes(sec):
nd.include_flux(1000.0)
Moderators: hines, wwlytton, ramcdougal
Code: Select all
for nd in self.ip3.nodes(sec):
nd.include_flux(1000.0)
Code: Select all
from neuron import h, rxd
h.load_file("stdrun.hoc")
sec = h.Section(name="sec")
sec.L = 10
sec.nseg = 11
sec.diam = 5
cyt = rxd.Region(h.allsec(), name="cyt", nrn_region="i")
ip3 = rxd.Species(cyt, name="ip3", initial=0)
# Used to rescale rate from mM/ms to molecules/ms used in include_flux
scale = 602214.129
# Specify the parameter values based using an anonymous function of the nodes
source_param = rxd.Parameter(
cyt,
name="sources",
value=lambda nd: 1.0/scale/nd.volume if nd.segment == sec(0.5) else 0)
# flux of 1000 molecules/ms restricted to sec(0.5)
my_include_flux = rxd.Rate(ip3, source_param * 1000)
h.finitialize(-70)
h.continuerun(10)
print(ip3.nodes.concentration)