I'm building a synapse model with a neuron and an astrocyte.
I would need to initialise the initial concentrations differently, however, not for the extracellular space.
I defined the ionic species as follow:
Code: Select all
rxd_na = rxd.Species([extracellular, cyt], name = 'na', charge = 1, d = 1)
rxd_k = rxd.Species([extracellular, cyt], name = 'k', charge = 1, d = 1)
rxd_ca = rxd.Species([extracellular, cyt], name = 'ca', charge = 2, d = 1)
na_a = rxd.Species([extracellular, cyt_astro], name='na_a', charge=1, d=1)
k_a = rxd.Species([extracellular, cyt_astro], name='k_a', charge=1, d=1)
ca_a = rxd.Species([extracellular, cyt_astro], name='ca_a', charge=2, d=1)
Code: Select all
h.cao0_ca_ion = 1.8
h.nao0_na_ion = 145
h.ko0_k_ion = 3
h.cai0_ca_ion = 0.006
h.nai0_na_ion = 5
h.ki0_k_ion = 130
h.finitialize(-65*mV)
Code: Select all
rxd_ca[extracellular].initial = 1.8
...
Thanks in advance for any help! :)