I'm trying to define a MOD file capable of reading the internal concentration of a substance and, with a slim to none delay, "transpose" the same concentration in a ecs.
The ecs is defined with rxd.Extracellular(). The substance is generated by a MOD file in the inside and as an rxd.species on the outside, with a diffusion value to prevent local accumulation.
I'm following a leak current and pump schemes:
Code: Select all
USEION x READ xi WRITE xo CHARGE 1
INITIAL {
imax = 1
}
ASSIGNED {
xo (mM)
xi (mM)
}
BREAKPOINT {
xo = imax *xi
}
When i use vec.record(cell.section[0](0.5)._ref_xo) it stays at 1 for one integration step and than it falls at 0 for the rest of the simulation.