Location: BG_Kr @ e95514b3ab1f / matlab_parameter_finding / calc_IGHK.m

Author:
Shelley Fong <s.fong@auckland.ac.nz>
Date:
2022-04-01 13:56:25+13:00
Desc:
With new densities found from fitting Pan - K + Kr to Clancy model
Permanent Source URI:
http://models.cellml.org/workspace/82c/rawfile/e95514b3ab1f0bb93bd69dff14346e3e15668e1f/matlab_parameter_finding/calc_IGHK.m

function I_GHK = calc_IGHK(G_GHK,V,cNi_st,cNo)
R = 8.314;
T = 310;
F = 96485;

V_norm = 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.*(cNi_st - cNo*exp(-V_norm));

end