k=lambda V: 1.0e3 * 3. * beta_cap(V * 1.0e3) * Qt)
CaPm2m1 = smodel.VDepSReac('CaPm2m1',
ssys,
slhs=[CaP_m2],
srhs=[CaP_m1],
k=lambda V: 1.0e3 * 2. * beta_cap(V * 1.0e3) * Qt)
CaPm1m0 = smodel.VDepSReac('CaPm1m0',
ssys,
slhs=[CaP_m1],
srhs=[CaP_m0],
k=lambda V: 1.0e3 * 1. * beta_cap(V * 1.0e3) * Qt)
if cyl160:
OC_CaP = smodel.GHKcurr('OC_CaP',
ssys,
CaP_m3,
Ca,
virtual_oconc=Ca_oconc,
computeflux=True)
else:
OC_CaP = smodel.GHKcurr('OC_CaP', ssys, CaP_m3, Ca, computeflux=True)
OC_CaP.setP(CaP_P)
######## CaT channel ##########
CaTchan = smodel.Chan('CaTchan', mdl)
CaT_m0h0 = smodel.ChanState('CaT_m0h0', mdl, CaTchan)
CaT_m0h1 = smodel.ChanState('CaT_m0h1', mdl, CaTchan)
CaT_m1h0 = smodel.ChanState('CaT_m1h0', mdl, CaTchan)
CaT_m1h1 = smodel.ChanState('CaT_m1h1', mdl, CaTchan)
k=lambda V: 1.0e3 * 3. * beta_cap(V * 1.0e3) * Qt)
CaPm2m1 = smodel.VDepSReac('CaPm2m1',
ssys_stoch,
slhs=[CaP_m2],
srhs=[CaP_m1],
k=lambda V: 1.0e3 * 2. * beta_cap(V * 1.0e3) * Qt)
CaPm1m0 = smodel.VDepSReac('CaPm1m0',
ssys_stoch,
slhs=[CaP_m1],
srhs=[CaP_m0],
k=lambda V: 1.0e3 * 1. * beta_cap(V * 1.0e3) * Qt)
if cyl160:
OC_CaP = smodel.GHKcurr('OC_CaP',
ssys_stoch,
CaP_m3,
Ca_stoch,
virtual_oconc=Ca_oconc,
computeflux=True)
else:
OC_CaP = smodel.GHKcurr('OC_CaP',
ssys_stoch,
CaP_m3,
Ca_stoch,
computeflux=True)
OC_CaP.setP(CaP_P)
######## CaT channel ##########
CaT_m0h0 = smodel.Spec('CaT_m0h0', mdl_det)
CaT_m0h1 = smodel.Spec('CaT_m0h1', mdl_det)
CaT_m1h0 = smodel.Spec('CaT_m1h0', mdl_det)
k=lambda V: 1.0e3 * 1. * betah_cat(V * 1.0e3))
CaTm1h1_m1h0 = smodel.VDepSReac('CaTm1h1_m1h0',
ssys_stoch,
slhs=[CaT_m1h1],
srhs=[CaT_m1h0],
k=lambda V: 1.0e3 * 1. * betah_cat(V * 1.0e3))
CaTm0h1_m0h0 = smodel.VDepSReac('CaTm0h1_m0h0',
ssys_stoch,
slhs=[CaT_m0h1],
srhs=[CaT_m0h0],
k=lambda V: 1.0e3 * 1. * betah_cat(V * 1.0e3))
if cyl160:
OC_CaT = smodel.GHKcurr('OC_CaT',
ssys_stoch,
CaT_m2h1,
Ca_stoch,
virtual_oconc=Ca_oconc,
computeflux=True)
else:
OC_CaT = smodel.GHKcurr('OC_CaT',
ssys_stoch,
CaT_m2h1,
Ca_stoch,
computeflux=True)
OC_CaT.setP(CaT_P)
##### BK channel ####################
BK_C0 = smodel.Spec('BK_C0', mdl_det)
BK_C1 = smodel.Spec('BK_C1', mdl_det)
BK_C2 = smodel.Spec('BK_C2', mdl_det)