- Author:
- WeiweiAi <wai484@aucklanduni.ac.nz>
- Date:
- 2021-07-27 14:55:57+12:00
- Desc:
- 1. Move phiT to the gating_kinetics;
2. Add time constant (time_const) in the model gate_steady;
3. Import cellLib to reuse the modules unit, time, gating_kinetics, phiT (temperature_factor), protocols, etc
4. Move units, time, temperature_factor (renamed to phiT), periodic-IStim-protocol, Voltage_ramp_protocol, gating-variable and gate_steady to cellLib
- Permanent Source URI:
- https://staging.physiomeproject.org/workspace/64f/rawfile/eb2aeec4e8d211b62d370fd78bd1fa1174bb13ab/sed-ml/get_init.py
def init_gate(T, V):
import math
phi = 3**((T - 6.3) / 10)
alpha_m = 0.1 * (V + 25) * phi / (math.exp( (V + 25) / 10) -1 )
beta_m = 4 * math.exp( V / 18 ) * phi
alpha_h = 0.07 * math.exp( V / 20 ) * phi
beta_h = phi / (math.exp( (V + 30) /10 ) + 1)
alpha_n = 0.01 * (V + 10) * phi / (math.exp( (V + 10) / 10) -1 )
beta_n = 0.125 * math.exp( V / 80 ) * phi
m_init = alpha_m / (alpha_m + beta_m)
n_init = alpha_n / (alpha_n + beta_n)
h_init = alpha_h / (alpha_h + beta_h)
return m_init, n_init, h_init
