Location: BG_TCC @ bc202e51bde6 / Clancy_matlab_parameter_fitting / calc_IGHK.m

Author:
Shelley Fong <sfon036@UoA.auckland.ac.nz>
Date:
2022-07-06 16:26:28+12:00
Desc:
Updating to LRd 34.4 pL cell parameters
Permanent Source URI:
https://models.cellml.org/workspace/831/rawfile/bc202e51bde6be8ddf88d35d908ef1465dcf067d/Clancy_matlab_parameter_fitting/calc_IGHK.m

function I_GHK = calc_IGHK(G_GHK,V,cNi_st,cNo,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.*(cNi_st - cNo*exp(-V_norm));

end