//Created by libAntimony v2.4

// Warnings from automatic translation:
//    Unable to use the formula "K_mCSQN - (CSQN_max + CSQN_buff + delta_Ca_JSR + Ca_JSR_old)" to set the assignment rule for b1:  Loop detected:  b1's definition (K_mCSQN - (CSQN_max + CSQN_buff + delta_Ca_JSR + Ca_JSR_old)) either includes itself directly (i.e. 's5 := 6 + s5') or by proxy (i.e. 's5 := 8*d3' and 'd3 := 9*s5').
//    Unable to use the formula "K_mCSQN * (CSQN_buff + delta_Ca_JSR + Ca_JSR_old)" to set the assignment rule for c1:  Loop detected:  c1's definition (K_mCSQN * (CSQN_buff + delta_Ca_JSR + Ca_JSR_old)) either includes itself directly (i.e. 's5 := 6 + s5') or by proxy (i.e. 's5 := 8*d3' and 'd3 := 9*s5').
//    Unable to use the formula "TRPN_buff + CMDN_buff + delta_Cai + Cai_old" to set the assignment rule for Ca_total:  Loop detected:  Ca_total's definition (TRPN_buff + CMDN_buff + delta_Cai + Cai_old) either includes itself directly (i.e. 's5 := 6 + s5') or by proxy (i.e. 's5 := 8*d3' and 'd3 := 9*s5').

model riemer_sobie_tung_1998_version01__environment(time_)

  // Variable initializations:
  time_ = ;
end

model riemer_sobie_tung_1998_version01__membrane(V, R, T, F, time_, i_Na, i_Ca_L, i_Ca_T, i_Kr, i_Ks, i_NaCa, i_K1, i_Kp, i_p_Ca, i_Na_b, i_Ca_b, i_NaK, i_ns_Ca, i_s)

  // Rate Rules:
  V' = -(1.0 / Cm) * (i_Na + i_Ca_L + i_Ca_T + i_Kr + i_Ks + i_K1 + i_Kp + i_NaCa + i_p_Ca + i_Na_b + i_Ca_b + i_NaK + i_ns_Ca + i_s + I_st);

  // Variable initializations:
  V = -86.5;
  R = 8314.0;
  T = 310.0;
  F = 96500.0;
  Cm = 1.0;
  I_st = -100.0;
  time_ = ;
  i_Na = ;
  i_Ca_L = ;
  i_Ca_T = ;
  i_Kr = ;
  i_Ks = ;
  i_NaCa = ;
  i_K1 = ;
  i_Kp = ;
  i_p_Ca = ;
  i_Na_b = ;
  i_Ca_b = ;
  i_NaK = ;
  i_ns_Ca = ;
  i_s = ;
end

model riemer_sobie_tung_1998_version01__fast_sodium_current_m_gate(m, V, time_)

  // Assignment Rules:
  alpha_m := (0.32 * (V + 47.13)) / (1.0 - exp(-0.1 * (V + 47.13)));
  beta_m := 0.08 * exp(- V / 11.0);

  // Rate Rules:
  m' = alpha_m * (1.0 - m) - beta_m * m;

  // Variable initializations:
  m = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__fast_sodium_current_h_gate(h, V, time_)

  // Assignment Rules:
  alpha_h := piecewise( 0.135 * exp((80.0 + V) / -6.8) , V < -40.0 , 0.0 );
  beta_h := piecewise( 3.56 * exp(0.079 * V) + 310000.0 * exp(0.35 * V) , V < -40.0 , 1.0 / (0.13 * (1.0 + exp(-(V + 10.66) / 11.1))) );

  // Rate Rules:
  h' = alpha_h * (1.0 - h) - beta_h * h;

  // Variable initializations:
  h = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__fast_sodium_current_j_gate(j, V, time_)

  // Assignment Rules:
  alpha_j := piecewise( (-127140.0 * exp(0.2444 * V) - 0.00003474 * exp(-0.04391 * V)) * ((V + 37.78) / (1.0 + exp(0.311 * (V + 79.23)))) , V < -40.0 , 0.0 );
  beta_j := piecewise( (0.1212 * exp(-0.01052 * V)) / (1.0 + exp(-0.1378 * (V + 40.14))) , V < -40.0 , (0.3 * exp(-0.0000002535 * V)) / (1.0 + exp(-0.1 * (V + 32.0))) );

  // Rate Rules:
  j' = alpha_j * (1.0 - j) - beta_j * j;

  // Variable initializations:
  j = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__fast_sodium_current(i_Na, E_Na, time_, V, R, F, T, Nai, Nao, m, h, j)

  // Sub-modules, and any changes to those submodules:
  fast_sodium_current_m_gate: riemer_sobie_tung_1998_version01__fast_sodium_current_m_gate(m, V, time_);
  fast_sodium_current_h_gate: riemer_sobie_tung_1998_version01__fast_sodium_current_h_gate(h, V, time_);
  fast_sodium_current_j_gate: riemer_sobie_tung_1998_version01__fast_sodium_current_j_gate(j, V, time_);

  // Assignment Rules:
  i_Na := g_Na * power(m, 3.0) * h * j * (V - E_Na);
  E_Na := ((R * T) / F) * ln(Nao / Nai);

  // Variable initializations:
  g_Na = 16.0;
  R = ;
  F = ;
  T = ;
  Nai = ;
  Nao = ;
end

model riemer_sobie_tung_1998_version01__L_type_Ca_channel_d_gate(d, V, time_)

  // Assignment Rules:
  alpha_d := d_infinity / tau_d;
  beta_d := (1.0 - d_infinity) / tau_d;
  d_infinity := 1.0 / (1.0 + exp(-((V + 10.0) / 6.24)));
  tau_d := d_infinity * ((1.0 - exp(-((V + 10.0) / 6.24))) / (0.035 * (V + 10.0)));

  // Rate Rules:
  d' = alpha_d * (1.0 - d) - beta_d * d;

  // Variable initializations:
  d = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__L_type_Ca_channel_f_gate(f, V, time_)

  // Assignment Rules:
  alpha_f := f_infinity / tau_f;
  beta_f := (1.0 - f_infinity) / tau_f;
  f_infinity := 1.0 / (1.0 + exp((V + 35.06) / 8.6)) + 0.6 / (1.0 + exp((50.0 - V) / 20.0));
  tau_f := 1.0 / (0.0197 * exp(-power(0.0337 * (V + 10.0), 2.0)) + 0.02);

  // Rate Rules:
  f' = alpha_f * (1.0 - f) - beta_f * f;

  // Variable initializations:
  f = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__L_type_Ca_channel_f_Ca_gate(f_Ca, V, time_, Cai)

  // Assignment Rules:
  f_Ca := 1.0 / (1.0 + Cai / Km_Ca);

  // Variable initializations:
  Km_Ca = 0.6;
  V = ;
  time_ = ;
  Cai = ;
end

model riemer_sobie_tung_1998_version01__L_type_Ca_channel(i_Ca_L, i_CaCa, i_CaK, i_CaNa, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, time_, V, Cai, R, T, F, Cao, Nao, Ko, Nai, Ki, d, f, f_Ca)

  // Sub-modules, and any changes to those submodules:
  L_type_Ca_channel_d_gate: riemer_sobie_tung_1998_version01__L_type_Ca_channel_d_gate(d, V, time_);
  L_type_Ca_channel_f_gate: riemer_sobie_tung_1998_version01__L_type_Ca_channel_f_gate(f, V, time_);
  L_type_Ca_channel_f_Ca_gate: riemer_sobie_tung_1998_version01__L_type_Ca_channel_f_Ca_gate(f_Ca, V, time_, Cai);

  // Assignment Rules:
  i_Ca_L := i_CaCa + i_CaK + i_CaNa;
  i_CaCa := d * f * f_Ca * I_CaCa;
  i_CaK := d * f * f_Ca * I_CaK;
  i_CaNa := d * f * f_Ca * I_CaNa;
  I_CaCa := P_Ca * power(2.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Cai * Cai * exp((2.0 * V * F) / (R * T)) - gamma_Cao * Cao) / (exp((2.0 * V * F) / (R * T)) - 1.0));
  I_CaK := P_K * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Ki * Ki * exp((1.0 * V * F) / (R * T)) - gamma_Ko * Ko) / (exp((1.0 * V * F) / (R * T)) - 1.0));
  I_CaNa := P_Na * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Nai * Nai * exp((1.0 * V * F) / (R * T)) - gamma_Nao * Nao) / (exp((1.0 * V * F) / (R * T)) - 1.0));

  // Variable initializations:
  gamma_Nai = 0.75;
  gamma_Nao = 0.75;
  gamma_Ki = 0.75;
  gamma_Ko = 0.75;
  P_Ca = 0.00054;
  P_Na = 0.000000675;
  P_K = 0.000000193;
  gamma_Cai = 1.0;
  gamma_Cao = 0.34;
  R = ;
  T = ;
  F = ;
  Cao = ;
  Nao = ;
  Ko = ;
  Nai = ;
  Ki = ;
end

model riemer_sobie_tung_1998_version01__T_type_Ca_channel_b_gate(b, time_, V)

  // Assignment Rules:
  b_infinity := 1.0 / (1.0 + exp(-((V + 14.0) / 10.8)));
  tau_b := (3.7 + 6.1) / (1.0 + exp((25.0 + V) / 4.5));

  // Rate Rules:
  b' = (b_infinity - b) / tau_b;

  // Variable initializations:
  b = ;
  time_ = ;
  V = ;
end

model riemer_sobie_tung_1998_version01__T_type_Ca_channel_g_gate(g, time_, V)

  // Assignment Rules:
  g_infinity := 1.0 / (1.0 + exp((V + 60.0) / 5.6));
  tau_g := piecewise( 12.0 , V > 0.0 , -0.875 * V + 12.0 );

  // Rate Rules:
  g' = (g_infinity - g) / tau_g;

  // Variable initializations:
  g = ;
  time_ = ;
  V = ;
end

model riemer_sobie_tung_1998_version01__T_type_Ca_channel(i_Ca_T, time_, V, Cai, Cao, R, T, F, b, g)

  // Sub-modules, and any changes to those submodules:
  T_type_Ca_channel_b_gate: riemer_sobie_tung_1998_version01__T_type_Ca_channel_b_gate(b, time_, V);
  T_type_Ca_channel_g_gate: riemer_sobie_tung_1998_version01__T_type_Ca_channel_g_gate(g, time_, V);

  // Assignment Rules:
  i_Ca_T := g_Ca_T * power(b, 2.0) * g * (V - E_Ca);
  E_Ca := ((R * T) / (2.0 * F)) * ln(Cao / Cai);

  // Variable initializations:
  g_Ca_T = 0.05;
  Cai = ;
  Cao = ;
  R = ;
  T = ;
  F = ;
end

model riemer_sobie_tung_1998_version01__rapid_time_dependent_potassium_current_Xr_gate(Xr, V, time_)

  // Assignment Rules:
  Xr_infinity := 1.0 / (1.0 + exp(-((V + 21.5) / 7.5)));
  tau_Xr := 1.0 / ((0.00138 * (V + 14.2)) / (1.0 - exp(-0.123 * (V + 14.2))) + (0.00061 * (V + 38.9)) / (exp(0.145 * (V + 38.9)) - 1.0));

  // Rate Rules:
  Xr' = (Xr_infinity - Xr) / tau_Xr;

  // Variable initializations:
  Xr = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__rapid_time_dependent_potassium_current_Rr_gate(Rr, V, time_)

  // Assignment Rules:
  Rr := 1.0 / (1.0 + exp((V + 9.0) / 22.4));

  // Variable initializations:
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__rapid_time_dependent_potassium_current(i_Kr, time_, V, R, T, F, Ko, Ki, Xr, Rr)

  // Sub-modules, and any changes to those submodules:
  rapid_time_dependent_potassium_current_Xr_gate: riemer_sobie_tung_1998_version01__rapid_time_dependent_potassium_current_Xr_gate(Xr, V, time_);
  rapid_time_dependent_potassium_current_Rr_gate: riemer_sobie_tung_1998_version01__rapid_time_dependent_potassium_current_Rr_gate(Rr, V, time_);

  // Assignment Rules:
  i_Kr := g_Kr * Xr * Rr * (V - E_Kr);
  g_Kr := 0.02614 * root(Ko / 5.4);
  E_Kr := ((R * T) / F) * ln(Ko / Ki);

  // Variable initializations:
  R = ;
  T = ;
  F = ;
  Ko = ;
  Ki = ;
end

model riemer_sobie_tung_1998_version01__slow_time_dependent_potassium_current_Xs_gate(Xs, tau_Xs, V, time_)

  // Assignment Rules:
  tau_Xs := 1.0 / ((0.0000719 * (V + 30.0)) / (1.0 - exp(-0.148 * (V + 30.0))) + (0.000131 * (V + 30.0)) / (exp(0.0687 * (V + 30.0)) - 1.0));
  Xs_infinity := 1.0 / (1.0 + exp(-(V + -1.5) / 16.7));

  // Rate Rules:
  Xs' = (Xs_infinity - Xs) / tau_Xs;

  // Variable initializations:
  Xs = ;
  V = ;
  time_ = ;
end

model riemer_sobie_tung_1998_version01__slow_time_dependent_potassium_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, Cai, R, T, F, Xs)

  // Sub-modules, and any changes to those submodules:
  slow_time_dependent_potassium_current_Xs_gate: riemer_sobie_tung_1998_version01__slow_time_dependent_potassium_current_Xs_gate(Xs, tau_Xs, V, time_);

  // Assignment Rules:
  i_Ks := g_Ks * power(Xs, 2.0) * (V - E_Ks);
  g_Ks := (0.057 + 0.19) / (1.0 + exp((-7.2 + P_Ca) / 0.6));
  E_Ks := ((R * T) / F) * ln((Ko + P_NaK * Nao) / (Ki + P_NaK * Nai));
  P_Ca := -log(Cai) + 3.0;

  // Variable initializations:
  P_NaK = 0.01833;
  Ko = ;
  Ki = ;
  Nao = ;
  Nai = ;
  Cai = ;
  R = ;
  T = ;
  F = ;
end

model riemer_sobie_tung_1998_version01__time_independent_potassium_current_K1_gate(K1_infinity, V, time_, E_K1)

  // Assignment Rules:
  K1_infinity := alpha_K1 / (alpha_K1 + beta_K1);
  alpha_K1 := 1.02 / (1.0 + exp(0.2385 * ((V - E_K1) - 59.215)));
  beta_K1 := (0.49124 * exp(((V + 5.476) - E_K1) / 12.45) + exp((V - (E_K1 + 594.31)) / 16.2)) / (1.0 + exp(-0.5143 * (V - (E_K1 + 4.753))));

  // Variable initializations:
  V = ;
  time_ = ;
  E_K1 = ;
end

model riemer_sobie_tung_1998_version01__time_independent_potassium_current(i_K1, E_K1, time_, V, Ko, Ki, R, T, F, K1_infinity)

  // Sub-modules, and any changes to those submodules:
  time_independent_potassium_current_K1_gate: riemer_sobie_tung_1998_version01__time_independent_potassium_current_K1_gate(K1_infinity, V, time_, E_K1);

  // Assignment Rules:
  i_K1 := g_K1 * K1_infinity * (V - E_K1);
  E_K1 := ((R * T) / F) * ln(Ko / Ki);
  g_K1 := 0.75 * root(Ko / 5.4);

  // Variable initializations:
  Ko = ;
  Ki = ;
  R = ;
  T = ;
  F = ;
end

model riemer_sobie_tung_1998_version01__plateau_potassium_current(i_Kp, time_, V, E_K1)

  // Assignment Rules:
  i_Kp := g_Kp * Kp * (V - E_Kp);
  E_Kp := E_K1;
  Kp := 1.0 / (1.0 + exp((7.488 - V) / 5.98));

  // Variable initializations:
  g_Kp = 0.00552;
  time_ = ;
  V = ;
  E_K1 = ;
end

model riemer_sobie_tung_1998_version01__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai)

  // Assignment Rules:
  i_p_Ca := I_pCa * (Cai / (K_mpCa + Cai));

  // Variable initializations:
  K_mpCa = 0.5;
  I_pCa = 1.15;
  time_ = ;
  Cai = ;
end

model riemer_sobie_tung_1998_version01__sodium_background_current(i_Na_b, time_, V, E_Na)

  // Assignment Rules:
  i_Na_b := g_Nab * (V - E_NaN);
  E_NaN := E_Na;

  // Variable initializations:
  g_Nab = 0.00141;
  time_ = ;
  V = ;
  E_Na = ;
end

model riemer_sobie_tung_1998_version01__calcium_background_current(i_Ca_b, time_, V, R, T, F, Cai, Cao)

  // Assignment Rules:
  i_Ca_b := g_Cab * (V - E_CaN);
  E_CaN := ((R * T) / (2.0 * F)) * ln(Cao / Cai);

  // Variable initializations:
  g_Cab = 0.003016;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Cai = ;
  Cao = ;
end

model riemer_sobie_tung_1998_version01__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko)

  // Assignment Rules:
  i_NaK := I_NaK * f_NaK * (1.0 / (1.0 + power(K_mNai / Nai, 1.5))) * (Ko / (Ko + K_mKo));
  f_NaK := 1.0 / ((1.0 + 0.1245 * exp(-0.1 * ((V * F) / (R * T)))) + 0.0365 * sigma * exp(-((V * F) / (R * T))));
  sigma := (1.0 / 7.0) * (exp(Nao / 67.3) - 1.0);

  // Variable initializations:
  I_NaK = ;
  K_mNai = 10.0;
  K_mKo = 1.5;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Nai = ;
  Nao = ;
  Ko = ;
end

model riemer_sobie_tung_1998_version01__non_specific_calcium_activated_current(i_ns_Ca, i_ns_Na, i_ns_K, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, R, T, F, Cao, Nao, Ko, Nai, Ki, time_, Cai, V)

  // Assignment Rules:
  i_ns_Ca := i_ns_Na + i_ns_K;
  i_ns_Na := I_ns_Na * (1.0 / (1.0 + power(K_m_ns_Ca / Cai, 3.0)));
  i_ns_K := I_ns_K * (1.0 / (1.0 + power(K_m_ns_Ca / Cai, 3.0)));
  I_ns_Na := P_ns_Ca * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Nai * Nai * exp((1.0 * V * F) / (R * T)) - gamma_Nao * Nao) / (exp((1.0 * V * F) / (R * T)) - 1.0));
  I_ns_K := P_ns_Ca * power(1.0, 2.0) * ((V * power(F, 2.0)) / (R * T)) * ((gamma_Ki * Ki * exp((1.0 * V * F) / (R * T)) - gamma_Ko * Ko) / (exp((1.0 * V * F) / (R * T)) - 1.0));

  // Variable initializations:
  P_ns_Ca = 1.75E-7;
  gamma_Nai = ;
  gamma_Nao = ;
  gamma_Ki = ;
  gamma_Ko = ;
  R = ;
  T = ;
  F = ;
  Cao = ;
  Nao = ;
  Ko = ;
  Nai = ;
  Ki = ;
  K_m_ns_Ca = 1.2;
  time_ = ;
  Cai = ;
  V = ;
end

model riemer_sobie_tung_1998_version01__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao)

  // Assignment Rules:
  i_NaCa := K_NaCa * (1.0 / (power(K_mNa, 3.0) + power(Nao, 3.0))) * (1.0 / (K_mCa + Cao)) * (1.0 / (1.0 + K_sat * exp((eta - 1.0) * V * (F / (R * T))))) * (exp(eta * V * (F / (R * T))) * power(Nai, 3.0) * Cao - exp((eta - 1.0) * V * (F / (R * T))) * power(Nao, 3.0) * Cai);

  // Variable initializations:
  K_NaCa = 2000.0;
  K_mNa = 87.5;
  K_mCa = 1.38;
  K_sat = 0.1;
  eta = 0.35;
  time_ = ;
  V = ;
  R = ;
  T = ;
  F = ;
  Nai = ;
  Nao = ;
  Cai = ;
  Cao = ;
end

model riemer_sobie_tung_1998_version01__stretch_activated_current(i_s, time_, V)

  // Assignment Rules:
  i_s := g_s * (V - E_s);

  // Variable initializations:
  g_s = ;
  E_s = -20.0;
  time_ = ;
  V = ;
end

model riemer_sobie_tung_1998_version01__calcium_buffers_in_the_myoplasm(Ca_JSR_new, Cai, time_, CSQN_max, K_mCSQN, CSQN_buff)

  // Assignment Rules:
  Ca_JSR_new := (root(power(b1, 2.0) + 4.0 * c1) - b1) / 2.0;
  Cai := (2.0 / 3.0) * root(power(b, 2.0) - 3.0 * c) * cos(arccos((9.0 * b * c - (2.0 * power(b, 3.0) + 27.0 * d)) / (2.0 * power(power(b, 2.0) - 3.0 * c, 3.0 / 2.0))) / 3.0) - b / 3.0;
  TRPN_buff := TRPN_max * (Cai / (Cai + K_mTRPN));
  CMDN_buff := CMDN_max * (Cai / (Cai + K_mCMDN));
  delta_Ca_JSR := Ca_JSR_old - Ca_JSR_new;
  delta_Cai := Cai_old - Cai;
  b := (CMDN_max + TRPN_max + K_mTRPN + K_mCMDN) - Ca_total;
  c := (K_mCMDN * K_mTRPN + TRPN_max * K_mCMDN + CMDN_max * K_mTRPN) - Ca_total * (K_mTRPN + K_mCMDN);
  d := -(K_mTRPN * K_mCMDN) * Ca_total;

  // Variable initializations:
  K_mTRPN = 0.5;
  K_mCMDN = 2.38;
  TRPN_max = 70.0;
  CMDN_max = 50.0;
  Ca_total = ;
  Ca_JSR_old = ;
  Cai_old = ;
  b1 = ;
  c1 = ;
  time_ = ;
  CSQN_max = ;
  K_mCSQN = ;
  CSQN_buff = ;
end

model riemer_sobie_tung_1998_version01__calcium_fluxes_in_the_SR(i_rel, i_up, i_leak, i_tr, CSQN_buff, CSQN_max, K_mCSQN, time_, Cai, Ca_JSR_new, Ca_NSR)

  // Assignment Rules:
  i_rel := G_rel * (Ca_JSR_new - Cai);
  i_up := I_up * (Cai / (Cai + K_mup));
  i_leak := K_leak * Ca_NSR;
  i_tr := (Ca_NSR - Ca_JSR_new) / tau_tr;
  CSQN_buff := CSQN_max * (Ca_JSR_new / (Ca_JSR_new + K_mCSQN));
  G_rel := piecewise( G_rel_max * ((delta_Ca_i2 - delta_Ca_ith) / ((K_mrel + delta_Ca_i2) - delta_Ca_ith)) * (1.0 - exp(-(t / tau_on))) * exp(-(t / tau_off)) , calcium_overload == 0.0 , G_rel_max * (1.0 - exp(-(t / tau_on))) * exp(-(t / tau_off)) );
  G_rel_max := piecewise( piecewise( 0.0 , delta_Ca_i2 < delta_Ca_ith , 60.0 ) , calcium_overload == 0.0 , piecewise( 0.0 , CSQN_buff < CSQN_th , 4.0 ) );
  K_leak := I_up / Ca_NSR_max;

  // Variable initializations:
  CSQN_max = 10.0;
  K_mCSQN = 0.8;
  tau_on = 2.0;
  tau_off = 2.0;
  tau_tr = 180.0;
  t = 0.0;
  K_mrel = 0.8;
  delta_Ca_i2 = ;
  delta_Ca_ith = 0.18;
  CSQN_th = 0.7;
  K_mup = 0.92;
  I_up = 0.005;
  Ca_NSR_max = 0.15;
  calcium_overload = ;
  time_ = ;
  Cai = ;
  Ca_JSR_new = ;
  Ca_NSR = ;
end

model riemer_sobie_tung_1998_version01__ionic_concentrations(Nai, Nao, Cao, Ki, Ko, Ca_NSR, time_, F, i_Na, i_CaNa, i_CaCa, i_Na_b, i_ns_Na, i_NaCa, i_NaK, i_CaK, i_Kr, i_Ks, i_K1, i_Kp, i_ns_K, i_tr, i_leak, i_up)

  // Rate Rules:
  Nai' = -(i_Na + i_CaNa + i_Na_b + i_ns_Na + i_NaCa * 3.0 + i_NaK * 3.0) * (A_cap / (V_myo * F));
  Ki' = -(i_CaK + i_Kr + i_Ks + i_K1 + i_Kp + i_ns_K + -(i_NaK * 2.0)) * (A_cap / (V_myo * F));
  Ko' = (i_CaK + i_Kr + i_Ks + i_K1 + i_Kp + i_ns_K + -(i_NaK * 2.0)) * (A_cap / (V_cleft * F));
  Ca_NSR' = -((i_leak + i_tr) - i_up);
  Ca_foot' = - i_CaCa * (A_cap / (2.0 * V_myo * F)) * R_A_V;

  // Variable initializations:
  Nai = 18.0;
  Nao = 140.0;
  Cao = 1.8;
  Ki = 145.0;
  Ko = 5.4;
  Ca_NSR = 15.0;
  Ca_foot = ;
  A_cap = 0.000153;
  R_A_V = ;
  V_myo = ;
  V_cleft = ;
  time_ = ;
  F = ;
  i_Na = ;
  i_CaNa = ;
  i_CaCa = ;
  i_Na_b = ;
  i_ns_Na = ;
  i_NaCa = ;
  i_NaK = ;
  i_CaK = ;
  i_Kr = ;
  i_Ks = ;
  i_K1 = ;
  i_Kp = ;
  i_ns_K = ;
  i_tr = ;
  i_leak = ;
  i_up = ;
end

model *riemer_sobie_tung_1998_version01____main()

  // Sub-modules, and any changes to those submodules:
  environment: riemer_sobie_tung_1998_version01__environment(time_);
  membrane: riemer_sobie_tung_1998_version01__membrane(V, R, T, F, time_, i_Na, i_Ca_L, i_Ca_T, i_Kr, i_Ks, i_NaCa, i_K1, i_Kp, i_p_Ca, i_Na_b, i_Ca_b, i_NaK, i_ns_Ca, i_s);
  fast_sodium_current: riemer_sobie_tung_1998_version01__fast_sodium_current(i_Na, E_Na, time_, V, R, F, T, Nai, Nao, m, h, j);
  L_type_Ca_channel: riemer_sobie_tung_1998_version01__L_type_Ca_channel(i_Ca_L, i_CaCa, i_CaK, i_CaNa, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, time_, V, Cai, R, T, F, Cao, Nao, Ko, Nai, Ki, d, f, f_Ca);
  T_type_Ca_channel: riemer_sobie_tung_1998_version01__T_type_Ca_channel(i_Ca_T, time_, V, Cai, Cao, R, T, F, b, g);
  rapid_time_dependent_potassium_current: riemer_sobie_tung_1998_version01__rapid_time_dependent_potassium_current(i_Kr, time_, V, R, T, F, Ko, Ki, Xr, Rr);
  slow_time_dependent_potassium_current: riemer_sobie_tung_1998_version01__slow_time_dependent_potassium_current(i_Ks, time_, V, Ko, Ki, Nao, Nai, Cai, R, T, F, Xs);
  time_independent_potassium_current: riemer_sobie_tung_1998_version01__time_independent_potassium_current(i_K1, E_K1, time_, V, Ko, Ki, R, T, F, K1_infinity);
  plateau_potassium_current: riemer_sobie_tung_1998_version01__plateau_potassium_current(i_Kp, time_, V, E_K1);
  sarcolemmal_calcium_pump: riemer_sobie_tung_1998_version01__sarcolemmal_calcium_pump(i_p_Ca, time_, Cai);
  sodium_background_current: riemer_sobie_tung_1998_version01__sodium_background_current(i_Na_b, time_, V, E_Na);
  calcium_background_current: riemer_sobie_tung_1998_version01__calcium_background_current(i_Ca_b, time_, V, R, T, F, Cai, Cao);
  sodium_potassium_pump: riemer_sobie_tung_1998_version01__sodium_potassium_pump(i_NaK, time_, V, R, T, F, Nai, Nao, Ko);
  non_specific_calcium_activated_current: riemer_sobie_tung_1998_version01__non_specific_calcium_activated_current(i_ns_Ca, i_ns_Na, i_ns_K, gamma_Nai, gamma_Nao, gamma_Ki, gamma_Ko, R, T, F, Cao, Nao, Ko, Nai, Ki, time_, Cai, V);
  Na_Ca_exchanger: riemer_sobie_tung_1998_version01__Na_Ca_exchanger(i_NaCa, time_, V, R, T, F, Nai, Nao, Cai, Cao);
  stretch_activated_current: riemer_sobie_tung_1998_version01__stretch_activated_current(i_s, time_, V);
  calcium_buffers_in_the_myoplasm: riemer_sobie_tung_1998_version01__calcium_buffers_in_the_myoplasm(Ca_JSR_new, Cai, time_, CSQN_max, K_mCSQN, CSQN_buff);
  calcium_fluxes_in_the_SR: riemer_sobie_tung_1998_version01__calcium_fluxes_in_the_SR(i_rel, i_up, i_leak, i_tr, CSQN_buff, CSQN_max, K_mCSQN, time_, Cai, Ca_JSR_new, Ca_NSR);
  ionic_concentrations: riemer_sobie_tung_1998_version01__ionic_concentrations(Nai, Nao, Cao, Ki, Ko, Ca_NSR, time_, F, i_Na, i_CaNa, i_CaCa, i_Na_b, i_ns_Na, i_NaCa, i_NaK, i_CaK, i_Kr, i_Ks, i_K1, i_Kp, i_ns_K, i_tr, i_leak, i_up);
end