//Created by libAntimony v2.4 model gross_1998__environment(time_) // Variable initializations: time_ = ; end model gross_1998__parameters(f0_1, b0_1, f0_2, b0_2, f0_3, b0_3, f0_4, b0_4, f0_5, b0_5, f0_6, b0_6, delta_z, delta_m, delta_n, alpha_p, alpha_pp, beta_p, beta_pp, C_T, n, m, F, R, T, V, u) // Assignment Rules: b0_6 := (f0_1 * f0_2 * f0_3 * f0_4 * f0_5 * f0_6) / (b0_1 * b0_2 * b0_3 * b0_4 * b0_5); delta_m := (1 - beta_p) - beta_pp; delta_n := (1 - alpha_p) - alpha_pp; u := (F * V) / (R * T); // Variable initializations: f0_1 = 5.9e-1; b0_1 = 40; f0_2 = 2.5e-4; b0_2 = 3.7e3; f0_3 = 1.5e3; b0_3 = 2.4e3; f0_4 = 4.6e3; b0_4 = 5.3e-4; f0_5 = 0.2; b0_5 = 35e-3; f0_6 = 0.2; delta_z = 1.5; alpha_p = 0; alpha_pp = 0; beta_p = 0.11; beta_pp = 0; C_T = 2; n = 1; m = 3; F = 96485.34; R = 8.314; T = 310; V = ; end model gross_1998__ion_concentrations(V, Na_i, Na_o, HCO3_i, HCO3_o) // Variable initializations: V = -60e-3; Na_i = 20; Na_o = 0; HCO3_i = 18; HCO3_o = 0.6; end model gross_1998__rate_constants(Na_i, Na_o, HCO3_i, HCO3_o, f_1, b_1, f_2, b_2, f_3, b_3, f_4, b_4, f_5, b_5, f_6, b_6, f0_1, b0_1, f0_2, b0_2, f0_3, b0_3, f0_4, b0_4, f0_5, b0_5, f0_6, b0_6, delta_z, delta_m, delta_n, alpha_p, alpha_pp, beta_p, beta_pp, C_T, n, m, u) // Assignment Rules: f_1 := f0_1 * power(Na_i, n) * exp((n * alpha_p * u) / 2); b_1 := b0_1 * exp((- n * alpha_p * u) / 2); f_2 := f0_2 * power(HCO3_i, m) * exp((- m * beta_p * u) / 2); b_2 := b0_2 * exp((m * beta_p * u) / 2); f_3 := f0_3 * exp((((delta_z + n * delta_n) - m * delta_m) * u) / 2); b_3 := b0_3 * exp((-((delta_z + n * delta_n) - m * delta_m) * u) / 2); f_4 := f0_4 * exp((- m * beta_pp * u) / 2); b_4 := b0_4 * power(HCO3_o, m) * exp((m * beta_pp * u) / 2); f_5 := f0_5 * exp((n * alpha_pp * u) / 2); b_5 := b0_5 * power(Na_o, n) * exp((- n * alpha_pp * u) / 2); f_6 := f0_6 * exp((- delta_z * u) / 2); b_6 := b0_6 * exp((delta_z * u) / 2); // Variable initializations: Na_i = ; Na_o = ; HCO3_i = ; HCO3_o = ; f0_1 = ; b0_1 = ; f0_2 = ; b0_2 = ; f0_3 = ; b0_3 = ; f0_4 = ; b0_4 = ; f0_5 = ; b0_5 = ; f0_6 = ; b0_6 = ; delta_z = ; delta_m = ; delta_n = ; alpha_p = ; alpha_pp = ; beta_p = ; beta_pp = ; C_T = ; n = ; m = ; u = ; end model gross_1998__kinetic_equations(f_1, b_1, f_2, b_2, f_3, b_3, f_4, b_4, f_5, b_5, f_6, b_6, C_1, C_2, C_3, C_4, C_5, C_6, C_6_temp, C_T, time_) // Assignment Rules: C_6 := C_T - (C_1 + C_2 + C_3 + C_4 + C_5); // Rate Rules: C_1' = (b_1 * C_2 + f_6 * C_6) - (f_1 + b_6) * C_1; C_2' = (f_1 * C_1 + b_2 * C_3) - (b_1 + f_2) * C_2; C_3' = (f_2 * C_2 + b_3 * C_4) - (b_2 + f_3) * C_3; C_4' = (f_3 * C_3 + b_4 * C_5) - (b_3 + f_4) * C_4; C_5' = (f_4 * C_4 + b_5 * C_6) - (b_4 + f_5) * C_5; C_6_temp' = (f_5 * C_5 + b_6 * C_1) - (b_5 + f_6) * C_6; // Variable initializations: f_1 = ; b_1 = ; f_2 = ; b_2 = ; f_3 = ; b_3 = ; f_4 = ; b_4 = ; f_5 = ; b_5 = ; f_6 = ; b_6 = ; C_1 = 0.86691; C_2 = 0.510255; C_3 = 0.004611; C_4 = 0.0041975; C_5 = 0.5678; C_6_temp = 0.046227; C_T = ; time_ = ; end model gross_1998__king_altman_states(f_1, b_1, f_2, b_2, f_3, b_3, f_4, b_4, f_5, b_5, f_6, b_6, C_T, C1, C2, C3, C4, C5, C6) // Assignment Rules: C1_sum := f_2 * f_3 * f_4 * f_5 * f_6 + b_1 * f_3 * f_4 * f_5 * f_6 + b_1 * b_2 * f_4 * f_5 * f_6 + b_1 * b_2 * b_3 * f_5 * f_6 + b_1 * b_2 * b_3 * b_4 * f_6 + b_1 * b_2 * b_3 * b_4 * b_5; C2_sum := f_1 * f_3 * f_4 * f_5 * f_6 + b_2 * f_1 * f_4 * f_5 * f_6 + b_2 * b_3 * f_1 * f_5 * f_6 + b_2 * b_3 * b_4 * f_1 * f_6 + b_2 * b_3 * b_4 * b_5 * f_1 + b_2 * b_3 * b_4 * b_5 * b_6; C3_sum := f_1 * f_2 * f_4 * f_5 * f_6 + b_3 * f_1 * f_2 * f_5 * f_6 + b_3 * b_4 * f_1 * f_2 * f_6 + b_3 * b_4 * b_5 * f_1 * f_2 + b_3 * b_4 * b_5 * b_6 * f_2 + b_1 * b_3 * b_4 * b_5 * b_6; C4_sum := f_1 * f_2 * f_3 * f_5 * f_6 + b_4 * f_1 * f_2 * f_3 * f_6 + b_4 * b_5 * f_1 * f_2 * f_3 + b_4 * b_5 * b_6 * f_2 * f_3 + b_1 * b_4 * b_5 * b_6 * f_3 + b_1 * b_2 * b_4 * b_5 * b_6; C5_sum := f_1 * f_2 * f_3 * f_4 * f_6 + b_5 * f_1 * f_2 * f_3 * f_4 + b_5 * b_6 * f_2 * f_3 * f_4 + b_1 * b_5 * b_6 * f_3 * f_4 + b_1 * b_2 * b_5 * b_6 * f_4 + b_1 * b_2 * b_3 * b_5 * b_6; C6_sum := f_1 * f_2 * f_3 * f_4 * f_5 + b_6 * f_2 * f_3 * f_4 * f_5 + b_1 * b_6 * f_3 * f_4 * f_5 + b_1 * b_2 * b_6 * f_4 * f_5 + b_1 * b_2 * b_3 * b_6 * f_5 + b_1 * b_2 * b_3 * b_4 * b_6; C_sum := C1_sum + C2_sum + C3_sum + C4_sum + C5_sum + C6_sum; C1 := (C_T * C1_sum) / C_sum; C2 := (C_T * C2_sum) / C_sum; C3 := (C_T * C3_sum) / C_sum; C4 := (C_T * C4_sum) / C_sum; C5 := (C_T * C5_sum) / C_sum; C6 := (C_T * C6_sum) / C_sum; // Variable initializations: f_1 = ; b_1 = ; f_2 = ; b_2 = ; f_3 = ; b_3 = ; f_4 = ; b_4 = ; f_5 = ; b_5 = ; f_6 = ; b_6 = ; C_T = ; end model gross_1998__NBC_current(f_3, b_3, f_6, b_6, C_1, C_3, C_4, C_6, C1, C3, C4, C6, delta_z, n, m, F, I_NBC, I_NBC_SS) // Assignment Rules: I_NBC := - F * (delta_z * (b_6 * C_1 - f_6 * C_6) + ((delta_z + n) - m) * (f_3 * C_3 - b_3 * C_4)); I_NBC_SS := - F * (delta_z * (b_6 * C1 - f_6 * C6) + ((delta_z + n) - m) * (f_3 * C3 - b_3 * C4)); // Variable initializations: f_3 = ; b_3 = ; f_6 = ; b_6 = ; C_1 = ; C_3 = ; C_4 = ; C_6 = ; C1 = ; C3 = ; C4 = ; C6 = ; delta_z = ; n = ; m = ; F = ; end model *gross_1998____main() // Sub-modules, and any changes to those submodules: environment: gross_1998__environment(time_); parameters: gross_1998__parameters(f0_1, b0_1, f0_2, b0_2, f0_3, b0_3, f0_4, b0_4, f0_5, b0_5, f0_6, b0_6, delta_z, delta_m, delta_n, alpha_p, alpha_pp, beta_p, beta_pp, C_T, n, m, F, R, T, V, u); ion_concentrations: gross_1998__ion_concentrations(V, Na_i, Na_o, HCO3_i, HCO3_o); rate_constants: gross_1998__rate_constants(Na_i, Na_o, HCO3_i, HCO3_o, f_1, b_1, f_2, b_2, f_3, b_3, f_4, b_4, f_5, b_5, f_6, b_6, f0_1, b0_1, f0_2, b0_2, f0_3, b0_3, f0_4, b0_4, f0_5, b0_5, f0_6, b0_6, delta_z, delta_m, delta_n, alpha_p, alpha_pp, beta_p, beta_pp, C_T, n, m, u); kinetic_equations: gross_1998__kinetic_equations(f_1, b_1, f_2, b_2, f_3, b_3, f_4, b_4, f_5, b_5, f_6, b_6, C_1, C_2, C_3, C_4, C_5, C_6, C_6_temp, C_T, time_); king_altman_states: gross_1998__king_altman_states(f_1, b_1, f_2, b_2, f_3, b_3, f_4, b_4, f_5, b_5, f_6, b_6, C_T, C1, C2, C3, C4, C5, C6); NBC_current: gross_1998__NBC_current(f_3, b_3, f_6, b_6, C_1, C_3, C_4, C_6, C1, C3, C4, C6, delta_z, n, m, F, I_NBC, I_NBC_SS); end