Location: BG_ns_Ca @ 4bcdb0a9ff2f / matlab_parameter_fitting / calc_IGHK.m

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2022-03-30 11:46:41+13:00
Desc:
With channel density checked using patch clamp against kinetic model. Other fixes to be aligned with Clancy. also using Vns offset
Permanent Source URI:
https://staging.physiomeproject.org/workspace/83b/rawfile/4bcdb0a9ff2fc8bbd4d4723ca68b8df2f6cefc8e/matlab_parameter_fitting/calc_IGHK.m
Compare... Download Source 

function I_GHK = calc_IGHK(G_GHK,V,cKi_st,cKo, z)
R = 8.314;
T = 310;
F = 96485;

V_norm = z*F*V/R/T;
GHK_factor = V_norm./(1-exp(-V_norm));
idx_nan = isnan(GHK_factor);
GHK_factor(idx_nan) = 1;
I_GHK = G_GHK*GHK_factor.*(cKi_st - cKo*exp(-V_norm));

end