//Created by libAntimony v2.4
model cloutier_model_2009__environment(time_)

  // Variable initializations:
  time_ = ;
end

model cloutier_model_2009__NAn(NAn_, Vn_leak_Na, Vn_pump, Vn_stim, time_)

  // Rate Rules:
  NAn_' = (Vn_leak_Na + Vn_stim) - 3 * Vn_pump;

  // Variable initializations:
  NAn_ = 15.533;
  Vn_leak_Na = ;
  Vn_pump = ;
  Vn_stim = ;
  time_ = ;
end

model cloutier_model_2009__GLCn(GLCn, V_en_GLC, Vn_hk, time_)

  // Rate Rules:
  GLCn' = V_en_GLC - Vn_hk;

  // Variable initializations:
  GLCn = 0.2633;
  V_en_GLC = ;
  Vn_hk = ;
  time_ = ;
end

model cloutier_model_2009__G6Pn(G6Pn, Vn_pgi, Vn_hk, time_)

  // Rate Rules:
  G6Pn' = Vn_hk - Vn_pgi;

  // Variable initializations:
  G6Pn = 0.7275;
  Vn_pgi = ;
  Vn_hk = ;
  time_ = ;
end

model cloutier_model_2009__F6Pn(F6Pn, Vn_pgi, Vn_pfk, time_)

  // Rate Rules:
  F6Pn' = Vn_pgi - Vn_pfk;

  // Variable initializations:
  F6Pn = 0.1091;
  Vn_pgi = ;
  Vn_pfk = ;
  time_ = ;
end

model cloutier_model_2009__GAPn(GAPn, Vn_pgk, Vn_pfk, time_)

  // Rate Rules:
  GAPn' = 2 * Vn_pfk - Vn_pgk;

  // Variable initializations:
  GAPn = 0.0418;
  Vn_pgk = ;
  Vn_pfk = ;
  time_ = ;
end

model cloutier_model_2009__PEPn(PEPn, Vn_pgk, Vn_pk, time_)

  // Rate Rules:
  PEPn' = Vn_pgk - Vn_pk;

  // Variable initializations:
  PEPn = 0.0037;
  Vn_pgk = ;
  Vn_pk = ;
  time_ = ;
end

model cloutier_model_2009__PYRn(PYRn, Vn_pk, Vn_ldh, Vn_mito, time_)

  // Rate Rules:
  PYRn' = Vn_pk - (Vn_ldh + Vn_mito);

  // Variable initializations:
  PYRn = 0.0388;
  Vn_pk = ;
  Vn_ldh = ;
  Vn_mito = ;
  time_ = ;
end

model cloutier_model_2009__LACn(LACn, Vn_ldh, Vne_LAC, time_)

  // Rate Rules:
  LACn' = Vn_ldh - Vne_LAC;

  // Variable initializations:
  LACn = 0.3856;
  Vn_ldh = ;
  Vne_LAC = ;
  time_ = ;
end

model cloutier_model_2009__NADHn(NADHn, Vn_pgk, Vn_ldh, Vn_mito, time_)

  // Rate Rules:
  NADHn' = Vn_pgk - (Vn_ldh + Vn_mito);

  // Variable initializations:
  NADHn = 0.0319;
  Vn_pgk = ;
  Vn_ldh = ;
  Vn_mito = ;
  time_ = ;
end

model cloutier_model_2009__ATPn(ATPn, nOP, Vn_hk, Vn_pfk, Vn_pgk, Vn_pk, Vn_ATPase, Vn_pump, Vn_mito, Vn_ck, dAMP_dATPn, time_)

  // Rate Rules:
  ATPn' = ((Vn_pgk + Vn_pk + nOP * Vn_mito + Vn_ck) - (Vn_hk + Vn_pfk + Vn_ATPase + Vn_pump)) * power(1 - dAMP_dATPn, -1);

  // Variable initializations:
  ATPn = 2.2592;
  nOP = ;
  Vn_hk = ;
  Vn_pfk = ;
  Vn_pgk = ;
  Vn_pk = ;
  Vn_ATPase = ;
  Vn_pump = ;
  Vn_mito = ;
  Vn_ck = ;
  dAMP_dATPn = ;
  time_ = ;
end

model cloutier_model_2009__PCrn(PCrn, Vn_ck, time_)

  // Rate Rules:
  PCrn' = - Vn_ck;

  // Variable initializations:
  PCrn = 4.2529;
  Vn_ck = ;
  time_ = ;
end

model cloutier_model_2009__O2n(O2n, NAero, Vn_mito, Vcn_O2, time_)

  // Rate Rules:
  O2n' = Vcn_O2 - NAero * Vn_mito;

  // Variable initializations:
  O2n = 0.0975;
  NAero = ;
  Vn_mito = ;
  Vcn_O2 = ;
  time_ = ;
end

model cloutier_model_2009__GLUn(GLUn, Rng, Vg_gs, Vn_stim_GLU, time_)

  // Rate Rules:
  GLUn' = Vg_gs * (1 / Rng) - Vn_stim_GLU;

  // Variable initializations:
  GLUn = 3.0;
  Rng = ;
  Vg_gs = ;
  Vn_stim_GLU = ;
  time_ = ;
end

model cloutier_model_2009__NAg(NAg, Vg_leak_Na, Vg_pump, Veg_GLU, time_)

  // Rate Rules:
  NAg' = (Vg_leak_Na + 3 * Veg_GLU) - 3 * Vg_pump;

  // Variable initializations:
  NAg = 13.36;
  Vg_leak_Na = ;
  Vg_pump = ;
  Veg_GLU = ;
  time_ = ;
end

model cloutier_model_2009__GLCg(GLCg, Vcg_GLC, Veg_GLC, Vg_hk, time_)

  // Rate Rules:
  GLCg' = (Vcg_GLC + Veg_GLC) - Vg_hk;

  // Variable initializations:
  GLCg = 0.1656;
  Vcg_GLC = ;
  Veg_GLC = ;
  Vg_hk = ;
  time_ = ;
end

model cloutier_model_2009__G6Pg(G6Pg, Vg_pgi, Vg_glys, Vg_glyp, Vg_hk, time_)

  // Rate Rules:
  G6Pg' = (Vg_hk + Vg_glyp) - (Vg_pgi + Vg_glys);

  // Variable initializations:
  G6Pg = 0.7326;
  Vg_pgi = ;
  Vg_glys = ;
  Vg_glyp = ;
  Vg_hk = ;
  time_ = ;
end

model cloutier_model_2009__F6Pg(F6Pg, Vg_pgi, Vg_pfk, time_)

  // Rate Rules:
  F6Pg' = Vg_pgi - Vg_pfk;

  // Variable initializations:
  F6Pg = 0.1116;
  Vg_pgi = ;
  Vg_pfk = ;
  time_ = ;
end

model cloutier_model_2009__GAPg(GAPg, Vg_pgk, Vg_pfk, time_)

  // Rate Rules:
  GAPg' = 2 * Vg_pfk - Vg_pgk;

  // Variable initializations:
  GAPg = 0.0698;
  Vg_pgk = ;
  Vg_pfk = ;
  time_ = ;
end

model cloutier_model_2009__PEPg(PEPg, Vg_pgk, Vg_pk, time_)

  // Rate Rules:
  PEPg' = Vg_pgk - Vg_pk;

  // Variable initializations:
  PEPg = 0.0254;
  Vg_pgk = ;
  Vg_pk = ;
  time_ = ;
end

model cloutier_model_2009__PYRg(PYRg, Vg_ldh, Vg_pk, Vg_mito, time_)

  // Rate Rules:
  PYRg' = Vg_pk - (Vg_ldh + Vg_mito);

  // Variable initializations:
  PYRg = 0.1711;
  Vg_ldh = ;
  Vg_pk = ;
  Vg_mito = ;
  time_ = ;
end

model cloutier_model_2009__LACg(LACg, Vg_ldh, Vge_LAC, Vgc_LAC, time_)

  // Rate Rules:
  LACg' = Vg_ldh - (Vge_LAC + Vgc_LAC);

  // Variable initializations:
  LACg = 0.4651;
  Vg_ldh = ;
  Vge_LAC = ;
  Vgc_LAC = ;
  time_ = ;
end

model cloutier_model_2009__NADHg(NADHg, Vg_pgk, Vg_ldh, Vg_mito, time_)

  // Rate Rules:
  NADHg' = Vg_pgk - (Vg_ldh + Vg_mito);

  // Variable initializations:
  NADHg = 0.0445;
  Vg_pgk = ;
  Vg_ldh = ;
  Vg_mito = ;
  time_ = ;
end

model cloutier_model_2009__ATPg(ATPg, nOP, Vg_hk, Vg_pfk, Vg_pgk, Vg_pk, Vg_ATPase, Vg_pump, Vg_mito, Vg_ck, Vg_gs, dAMP_dATPg, time_)

  // Rate Rules:
  ATPg' = ((Vg_pgk + Vg_pk + nOP * Vg_mito + Vg_ck) - (Vg_hk + Vg_pfk + Vg_ATPase + Vg_pump + Vg_gs)) * power(1 - dAMP_dATPg, -1);

  // Variable initializations:
  ATPg = 2.24;
  nOP = ;
  Vg_hk = ;
  Vg_pfk = ;
  Vg_pgk = ;
  Vg_pk = ;
  Vg_ATPase = ;
  Vg_pump = ;
  Vg_mito = ;
  Vg_ck = ;
  Vg_gs = ;
  dAMP_dATPg = ;
  time_ = ;
end

model cloutier_model_2009__PCrg(PCrg, Vg_ck, time_)

  // Rate Rules:
  PCrg' = - Vg_ck;

  // Variable initializations:
  PCrg = 4.6817;
  Vg_ck = ;
  time_ = ;
end

model cloutier_model_2009__O2g(O2g, NAero, Vg_mito, Vcg_O2, time_)

  // Rate Rules:
  O2g' = Vcg_O2 - NAero * Vg_mito;

  // Variable initializations:
  O2g = 0.1589;
  NAero = ;
  Vg_mito = ;
  Vcg_O2 = ;
  time_ = ;
end

model cloutier_model_2009__GLYg(GLYg, Vg_glys, Vg_glyp, time_)

  // Rate Rules:
  GLYg' = Vg_glys - Vg_glyp;

  // Variable initializations:
  GLYg = 2.5;
  Vg_glys = ;
  Vg_glyp = ;
  time_ = ;
end

model cloutier_model_2009__GLUg(GLUg, Veg_GLU, Vg_gs, time_)

  // Rate Rules:
  GLUg' = Veg_GLU - Vg_gs;

  // Variable initializations:
  GLUg = 0.0;
  Veg_GLU = ;
  Vg_gs = ;
  time_ = ;
end

model cloutier_model_2009__GLCe(GLCe, Reg, Ren, Vce_GLC, Veg_GLC, V_en_GLC, time_)

  // Rate Rules:
  GLCe' = Vce_GLC - (Veg_GLC * (1 / Reg) + V_en_GLC * (1 / Ren));

  // Variable initializations:
  GLCe = 0.3339;
  Reg = ;
  Ren = ;
  Vce_GLC = ;
  Veg_GLC = ;
  V_en_GLC = ;
  time_ = ;
end

model cloutier_model_2009__LACe(LACe, Reg, Ren, Vec_LAC, Vne_LAC, Vge_LAC, time_)

  // Rate Rules:
  LACe' = (Vne_LAC * (1 / Ren) + Vge_LAC * (1 / Reg)) - Vec_LAC;

  // Variable initializations:
  LACe = 0.3986;
  Reg = ;
  Ren = ;
  Vec_LAC = ;
  Vne_LAC = ;
  Vge_LAC = ;
  time_ = ;
end

model cloutier_model_2009__GLUe(GLUe, Reg, Ren, Vn_stim_GLU, Veg_GLU, time_)

  // Rate Rules:
  GLUe' = Vn_stim_GLU * (1 / Ren) - Veg_GLU * (1 / Reg);

  // Variable initializations:
  GLUe = 0.0;
  Reg = ;
  Ren = ;
  Vn_stim_GLU = ;
  Veg_GLU = ;
  time_ = ;
end

model cloutier_model_2009__O2c(O2c, Rcn, Rcg, Vc_O2, Vcn_O2, Vcg_O2, time_)

  // Rate Rules:
  O2c' = Vc_O2 - (Vcn_O2 * (1 / Rcn) + Vcg_O2 * (1 / Rcg));

  // Variable initializations:
  O2c = 7.4201;
  Rcn = ;
  Rcg = ;
  Vc_O2 = ;
  Vcn_O2 = ;
  Vcg_O2 = ;
  time_ = ;
end

model cloutier_model_2009__GLCc(GLCc, Rce, Rcg, Vc_GLC, Vce_GLC, Vcg_GLC, time_)

  // Rate Rules:
  GLCc' = Vc_GLC - (Vce_GLC * (1 / Rce) + Vcg_GLC * (1 / Rcg));

  // Variable initializations:
  GLCc = 4.6401;
  Rce = ;
  Rcg = ;
  Vc_GLC = ;
  Vce_GLC = ;
  Vcg_GLC = ;
  time_ = ;
end

model cloutier_model_2009__LACc(LACc, Rce, Rcg, Vc_LAC, Vec_LAC, Vgc_LAC, time_)

  // Rate Rules:
  LACc' = Vc_LAC + (Vec_LAC * (1 / Rce) + Vgc_LAC * (1 / Rcg));

  // Variable initializations:
  LACc = 0.3251;
  Rce = ;
  Rcg = ;
  Vc_LAC = ;
  Vec_LAC = ;
  Vgc_LAC = ;
  time_ = ;
end

model cloutier_model_2009__CO2c(CO2c, Rcn, Rcg, Vnc_CO2, Vc_CO2, Vgc_CO2, time_)

  // Rate Rules:
  CO2c' = (Vnc_CO2 * (1 / Rcn) + Vgc_CO2 * (1 / Rcg)) - Vc_CO2;

  // Variable initializations:
  CO2c = 2.12;
  Rcn = ;
  Rcg = ;
  Vnc_CO2 = ;
  Vc_CO2 = ;
  Vgc_CO2 = ;
  time_ = ;
end

model cloutier_model_2009__Vv(Vv, Fin_t, Fout_t, time_)

  // Rate Rules:
  Vv' = Fin_t - Fout_t;

  // Variable initializations:
  Vv = 0.0237;
  Fin_t = ;
  Fout_t = ;
  time_ = ;
end

model cloutier_model_2009__dHb(dHb, Fin_t, Fout_t, Vv, O2a, O2c, time_)

  // Rate Rules:
  dHb' = Fin_t * (O2a - O2c) - Fout_t * (dHb / Vv);

  // Variable initializations:
  dHb = 0.0218;
  Fin_t = ;
  Fout_t = ;
  Vv = ;
  O2a = ;
  O2c = ;
  time_ = ;
end

model cloutier_model_2009__Vn_leak_Na(Vn_leak_Na, Sm_n, Vm, Vn, RT, F, NAe, NAn_)

  // Assignment Rules:
  Vn_leak_Na := (Sm_n / Vn) * (gn_NA / F) * ((RT / F) * ln(NAe / NAn_) - Vm);

  // Variable initializations:
  gn_NA = 0.0039;
  Sm_n = ;
  Vm = ;
  Vn = ;
  RT = ;
  F = ;
  NAe = ;
  NAn_ = ;
end

model cloutier_model_2009__Vn_pump(Vn_pump, kpump, Km_pump, Sm_n, Vn, ATPn, NAn_)

  // Assignment Rules:
  Vn_pump := (Sm_n / Vn) * kpump * ATPn * NAn_ * power(1 + ATPn / Km_pump, -1);

  // Variable initializations:
  kpump = ;
  Km_pump = ;
  Sm_n = ;
  Vn = ;
  ATPn = ;
  NAn_ = ;
end

model cloutier_model_2009__Vn_stim(Vn_stim, v_stim)

  // Assignment Rules:
  Vn_stim := v_stim;

  // Variable initializations:
  v_stim = ;
end

model cloutier_model_2009__V_en_GLC(V_en_GLC, GLCe, GLCn)

  // Assignment Rules:
  V_en_GLC := Vm_en_GLC * (GLCe / (GLCe + Km_en_GLC) - GLCn / (GLCn + Km_en_GLC));

  // Variable initializations:
  Km_en_GLC = 5.32;
  Vm_en_GLC = 0.50417;
  GLCe = ;
  GLCn = ;
end

model cloutier_model_2009__Vn_hk(Vn_hk, Km_GLC, GLCn, ATPn, G6Pn, G6P_inh_hk, aG6P_inh_hk)

  // Assignment Rules:
  Vn_hk := Vmax_n_hk * ATPn * (GLCn / (GLCn + Km_GLC)) * (1 - 1 / (1 + exp(- aG6P_inh_hk * (1 * (G6Pn - G6P_inh_hk)))));

  // Variable initializations:
  Vmax_n_hk = 0.0513;
  Km_GLC = ;
  GLCn = ;
  ATPn = ;
  G6Pn = ;
  G6P_inh_hk = ;
  aG6P_inh_hk = ;
end

model cloutier_model_2009__Vn_pgi(Vn_pgi, Km_G6P, Km_F6P_pgi, F6Pn, G6Pn)

  // Assignment Rules:
  Vn_pgi := Vmaxf_n_pgi * (G6Pn / (G6Pn + Km_G6P)) - Vmaxr_n_pgi * (F6Pn / (F6Pn + Km_F6P_pgi));

  // Variable initializations:
  Vmaxf_n_pgi = 0.5;
  Vmaxr_n_pgi = 0.45;
  Km_G6P = ;
  Km_F6P_pgi = ;
  F6Pn = ;
  G6Pn = ;
end

model cloutier_model_2009__Vn_pfk(Vn_pfk, Km_F6P_pfk, Ki_ATP, nH, ATPn, F6Pn)

  // Assignment Rules:
  Vn_pfk := kn_pfk * ATPn * (F6Pn / (F6Pn + Km_F6P_pfk)) * power(1 + power(ATPn / Ki_ATP, nH), -1);

  // Variable initializations:
  kn_pfk = 0.55783;
  Km_F6P_pfk = ;
  Ki_ATP = ;
  nH = ;
  ATPn = ;
  F6Pn = ;
end

model cloutier_model_2009__Vn_pgk(Vn_pgk, ADPn, GAPn, NADn, NADHn)

  // Assignment Rules:
  Vn_pgk := kn_pgk * GAPn * ADPn * (NADn / NADHn);

  // Variable initializations:
  kn_pgk = 0.4287;
  ADPn = ;
  GAPn = ;
  NADn = ;
  NADHn = ;
end

model cloutier_model_2009__Vn_pk(Vn_pk, PEPn, ADPn)

  // Assignment Rules:
  Vn_pk := kn_pk * PEPn * ADPn;

  // Variable initializations:
  kn_pk = 28.6;
  PEPn = ;
  ADPn = ;
end

model cloutier_model_2009__Vn_ldh(Vn_ldh, PYRn, NADHn, LACn, NADn)

  // Assignment Rules:
  Vn_ldh := kfn_ldh * PYRn * NADHn - krn_ldh * LACn * NADn;

  // Variable initializations:
  kfn_ldh = 5.30;
  krn_ldh = 0.1046;
  PYRn = ;
  NADHn = ;
  LACn = ;
  NADn = ;
end

model cloutier_model_2009__Vn_mito(Vn_mito, Km_O2, Km_ADP, Km_PYR, O2n, ADPn, PYRn, ATPn, rATP_mito, aATP_mito)

  // Assignment Rules:
  Vn_mito := Vmax_n_mito * (O2n / (O2n + Km_O2)) * (ADPn / (ADPn + Km_ADP)) * (PYRn / (PYRn + Km_PYR)) * (1 - 1 / (1 + exp(- aATP_mito * (1 * (ATPn / ADPn - 1 * rATP_mito)))));

  // Variable initializations:
  Vmax_n_mito = 0.05557;
  Km_O2 = ;
  Km_ADP = ;
  Km_PYR = ;
  O2n = ;
  ADPn = ;
  PYRn = ;
  ATPn = ;
  rATP_mito = ;
  aATP_mito = ;
end

model cloutier_model_2009__Vne_LAC(Vne_LAC, LACn, LACe)

  // Assignment Rules:
  Vne_LAC := Vmax_ne_LAC * (LACn / (LACn + Km_ne_LAC) - LACe / (LACe + Km_ne_LAC));

  // Variable initializations:
  Vmax_ne_LAC = 0.1978;
  Km_ne_LAC = 0.09314;
  LACn = ;
  LACe = ;
end

model cloutier_model_2009__Vn_ATPase(Vn_ATPase, ATPn)

  // Assignment Rules:
  Vn_ATPase := Vmax_n_ATPase * (ATPn / (ATPn + 0.001));

  // Variable initializations:
  Vmax_n_ATPase = 0.04889;
  ATPn = ;
end

model cloutier_model_2009__Vn_ck(Vn_ck, PCrn, ADPn, CRn, ATPn)

  // Assignment Rules:
  Vn_ck := kfn_ck * PCrn * ADPn - krn_ck * CRn * ATPn;

  // Variable initializations:
  krn_ck = 0.015;
  kfn_ck = 0.0524681;
  PCrn = ;
  ADPn = ;
  CRn = ;
  ATPn = ;
end

model cloutier_model_2009__Vcn_O2(Vcn_O2, Ko2, Vn, HbOP, O2c, O2n)

  // Assignment Rules:
  Vcn_O2 := (PScapn / Vn) * (Ko2 * power(HbOP / O2c - 1.0, -1.0 / nh_O2) - O2n);

  // Variable initializations:
  nh_O2 = 2.7;
  PScapn = 0.2202;
  Ko2 = ;
  Vn = ;
  HbOP = ;
  O2c = ;
  O2n = ;
end

model cloutier_model_2009__Vg_leak_Na(Vg_leak_Na, Sm_g, Vm, Vg, RT, F, NAe, NAg)

  // Assignment Rules:
  Vg_leak_Na := (Sm_g / Vg) * (gg_NA / F) * ((RT / F) * ln(NAe / NAg) - Vm);

  // Variable initializations:
  gg_NA = 0.00325;
  Sm_g = ;
  Vm = ;
  Vg = ;
  RT = ;
  F = ;
  NAe = ;
  NAg = ;
end

model cloutier_model_2009__Vg_pump(Vg_pump, kpump, Km_pump, Sm_g, Vg, ATPg, NAg)

  // Assignment Rules:
  Vg_pump := (Sm_g / Vg) * kpump * ATPg * NAg * power(1 + ATPg / Km_pump, -1);

  // Variable initializations:
  kpump = ;
  Km_pump = ;
  Sm_g = ;
  Vg = ;
  ATPg = ;
  NAg = ;
end

model cloutier_model_2009__Veg_GLC(Veg_GLC, GLCg, GLCe, KO1)

  // Assignment Rules:
  Veg_GLC := KO1 * Vm_eg_GLC * (GLCe / (GLCe + Km_eg_GLC) - GLCg / (GLCg + Km_eg_GLC));

  // Variable initializations:
  Km_eg_GLC = 3.53;
  Vm_eg_GLC = 0.038089;
  GLCg = ;
  GLCe = ;
  KO1 = ;
end

model cloutier_model_2009__Vcg_GLC(Vcg_GLC, GLCg, GLCc)

  // Assignment Rules:
  Vcg_GLC := Vm_cg_GLC * (GLCc / (GLCc + Km_cg_GLC) - GLCg / (GLCg + Km_cg_GLC));

  // Variable initializations:
  Km_cg_GLC = 9.92;
  Vm_cg_GLC = 0.0098394;
  GLCg = ;
  GLCc = ;
end

model cloutier_model_2009__Vg_hk(Vg_hk, Km_GLC, GLCg, ATPg, G6Pg, G6P_inh_hk, aG6P_inh_hk)

  // Assignment Rules:
  Vg_hk := Vmax_g_hk * ATPg * (GLCg / (GLCg + Km_GLC)) * (1 - 1 / (1 + exp(- aG6P_inh_hk * (1 * (G6Pg - G6P_inh_hk)))));

  // Variable initializations:
  Vmax_g_hk = 0.050461;
  Km_GLC = ;
  GLCg = ;
  ATPg = ;
  G6Pg = ;
  G6P_inh_hk = ;
  aG6P_inh_hk = ;
end

model cloutier_model_2009__Vg_pgi(Vg_pgi, Km_G6P, Km_F6P_pgi, F6Pg, G6Pg)

  // Assignment Rules:
  Vg_pgi := Vmaxf_g_pgi * (G6Pg / (G6Pg + Km_G6P)) - Vmaxr_g_pgi * (F6Pg / (F6Pg + Km_F6P_pgi));

  // Variable initializations:
  Vmaxf_g_pgi = 0.5;
  Vmaxr_g_pgi = 0.45;
  Km_G6P = ;
  Km_F6P_pgi = ;
  F6Pg = ;
  G6Pg = ;
end

model cloutier_model_2009__Vg_pfk(Vg_pfk, Km_F6P_pfk, Ki_ATP, nH, ATPg, F6Pg)

  // Assignment Rules:
  Vg_pfk := kg_pfk * ATPg * (F6Pg / (F6Pg + Km_F6P_pfk)) * power(1 + power(ATPg / Ki_ATP, nH), -1);

  // Variable initializations:
  kg_pfk = 0.403;
  Km_F6P_pfk = ;
  Ki_ATP = ;
  nH = ;
  ATPg = ;
  F6Pg = ;
end

model cloutier_model_2009__Vg_pgk(Vg_pgk, ADPg, GAPg, NADg, NADHg)

  // Assignment Rules:
  Vg_pgk := kg_pgk * GAPg * ADPg * (NADg / NADHg);

  // Variable initializations:
  kg_pgk = 0.2514;
  ADPg = ;
  GAPg = ;
  NADg = ;
  NADHg = ;
end

model cloutier_model_2009__Vg_pk(Vg_pk, PEPg, ADPg)

  // Assignment Rules:
  Vg_pk := kg_pk * PEPg * ADPg;

  // Variable initializations:
  kg_pk = 2.73;
  PEPg = ;
  ADPg = ;
end

model cloutier_model_2009__Vg_ldh(Vg_ldh, PYRg, NADHg, LACg, NADg)

  // Assignment Rules:
  Vg_ldh := kfg_ldh * PYRg * NADHg - krg_ldh * LACg * NADg;

  // Variable initializations:
  kfg_ldh = 6.2613;
  krg_ldh = 0.54682;
  PYRg = ;
  NADHg = ;
  LACg = ;
  NADg = ;
end

model cloutier_model_2009__Vg_mito(Vg_mito, Km_O2, Km_ADP, Km_PYR, O2g, ADPg, PYRg, ATPg, rATP_mito, aATP_mito)

  // Assignment Rules:
  Vg_mito := Vmax_g_mito * (O2g / (O2g + Km_O2)) * (ADPg / (ADPg + Km_ADP)) * (PYRg / (PYRg + Km_PYR)) * (1 - 1 / (1 + exp(1 * - aATP_mito * (ATPg / ADPg - 1 * rATP_mito))));

  // Variable initializations:
  Vmax_g_mito = 0.008454;
  Km_O2 = ;
  Km_ADP = ;
  Km_PYR = ;
  O2g = ;
  ADPg = ;
  PYRg = ;
  ATPg = ;
  rATP_mito = ;
  aATP_mito = ;
end

model cloutier_model_2009__Vge_LAC(Vge_LAC, LACg, LACe)

  // Assignment Rules:
  Vge_LAC := Vmax_ge_LAC * (LACg / (LACg + Km_ge_LAC) - LACe / (LACe + Km_ge_LAC));

  // Variable initializations:
  Vmax_ge_LAC = 0.086124;
  Km_ge_LAC = 0.22163;
  LACg = ;
  LACe = ;
end

model cloutier_model_2009__Vgc_LAC(Vgc_LAC, LACg, LACc)

  // Assignment Rules:
  Vgc_LAC := Vmax_gc_LAC * (LACg / (LACg + Km_gc_LAC) - LACc / (LACc + Km_gc_LAC));

  // Variable initializations:
  Vmax_gc_LAC = 0.00021856;
  Km_gc_LAC = 0.12862;
  LACg = ;
  LACc = ;
end

model cloutier_model_2009__Vg_ATPase(Vg_ATPase, ATPg)

  // Assignment Rules:
  Vg_ATPase := Vmax_g_ATPase * (ATPg / (ATPg + 0.001));

  // Variable initializations:
  Vmax_g_ATPase = 0.035657;
  ATPg = ;
end

model cloutier_model_2009__Vg_ck(Vg_ck, PCrg, ADPg, CRg, ATPg)

  // Assignment Rules:
  Vg_ck := kfg_ck * PCrg * ADPg - krg_ck * CRg * ATPg;

  // Variable initializations:
  krg_ck = 0.02073;
  kfg_ck = 0.0243;
  PCrg = ;
  ADPg = ;
  CRg = ;
  ATPg = ;
end

model cloutier_model_2009__Vcg_O2(Vcg_O2, Ko2, nh_O2, Vg, HbOP, O2c, O2g)

  // Assignment Rules:
  Vcg_O2 := (PScapg / Vg) * (Ko2 * power(HbOP / O2c - 1.0, -1.0 / nh_O2) - O2g);

  // Variable initializations:
  PScapg = 0.2457;
  Ko2 = ;
  nh_O2 = ;
  Vg = ;
  HbOP = ;
  O2c = ;
  O2g = ;
end

model cloutier_model_2009__Vc_O2(Vc_O2, Fin_t, Vc, O2a, O2c)

  // Assignment Rules:
  Vc_O2 := 2.0 * (Fin_t / Vc) * (O2a - O2c);

  // Variable initializations:
  Fin_t = ;
  Vc = ;
  O2a = ;
  O2c = ;
end

model cloutier_model_2009__Vc_GLC(Vc_GLC, Fin_t, Vc, GLCa, GLCc)

  // Assignment Rules:
  Vc_GLC := 2.0 * (Fin_t / Vc) * (GLCa - GLCc);

  // Variable initializations:
  Fin_t = ;
  Vc = ;
  GLCa = ;
  GLCc = ;
end

model cloutier_model_2009__Vce_GLC(Vce_GLC, GLCe, GLCc)

  // Assignment Rules:
  Vce_GLC := Vm_ce_GLC * (GLCc / (GLCc + Km_ce_GLC) - GLCe / (GLCe + Km_ce_GLC));

  // Variable initializations:
  Km_ce_GLC = 8.4568;
  Vm_ce_GLC = 0.0489;
  GLCe = ;
  GLCc = ;
end

model cloutier_model_2009__Vc_LAC(Vc_LAC, Fin_t, Vc, LACa, LACc)

  // Assignment Rules:
  Vc_LAC := 2.0 * (Fin_t / Vc) * (LACa - LACc);

  // Variable initializations:
  Fin_t = ;
  Vc = ;
  LACa = ;
  LACc = ;
end

model cloutier_model_2009__Vec_LAC(Vec_LAC, LACe, LACc)

  // Assignment Rules:
  Vec_LAC := Vm_ec_LAC * (LACe / (LACe + Km_ec_LAC) - LACc / (LACc + Km_ec_LAC));

  // Variable initializations:
  Km_ec_LAC = 0.764818;
  Vm_ec_LAC = 0.0325;
  LACe = ;
  LACc = ;
end

model cloutier_model_2009__Vnc_CO2(Vnc_CO2, Vn_mito)

  // Assignment Rules:
  Vnc_CO2 := 3.0 * Vn_mito;

  // Variable initializations:
  Vn_mito = ;
end

model cloutier_model_2009__Vgc_CO2(Vgc_CO2, Vg_mito)

  // Assignment Rules:
  Vgc_CO2 := 3.0 * Vg_mito;

  // Variable initializations:
  Vg_mito = ;
end

model cloutier_model_2009__Vn_stim_GLU(Vn_stim_GLU, R_GLU_NA, Km_GLU, GLUn, KO2, Vn_stim)

  // Assignment Rules:
  Vn_stim_GLU := Vn_stim * R_GLU_NA * KO2 * (GLUn / (GLUn + Km_GLU));

  // Variable initializations:
  R_GLU_NA = ;
  Km_GLU = ;
  GLUn = ;
  KO2 = ;
  Vn_stim = ;
end

model cloutier_model_2009__Vg_gs(Vg_gs, Km_GLU, Km_ATP, GLUg, ATPg)

  // Assignment Rules:
  Vg_gs := Vmax_g_gs * ((GLUg / (GLUg + Km_GLU)) * (ATPg / (ATPg + Km_ATP)));

  // Variable initializations:
  Vmax_g_gs = 0.3;
  Km_GLU = ;
  Km_ATP = ;
  GLUg = ;
  ATPg = ;
end

model cloutier_model_2009__Veg_GLU(Veg_GLU, Km_GLU, GLUe)

  // Assignment Rules:
  Veg_GLU := Vmax_eg_GLU * (GLUe / (GLUe + Km_GLU));

  // Variable initializations:
  Vmax_eg_GLU = 0.0208;
  Km_GLU = ;
  GLUe = ;
end

model cloutier_model_2009__Vc_CO2(Vc_CO2, Fin_t, Vc, CO2a, CO2c)

  // Assignment Rules:
  Vc_CO2 := 2.0 * (Fin_t / Vc) * (CO2c - CO2a);

  // Variable initializations:
  Fin_t = ;
  Vc = ;
  CO2a = ;
  CO2c = ;
end

model cloutier_model_2009__Vg_glys(Vg_glys, G6Pg, GLYg, GLY_inh, aGLY_inh)

  // Assignment Rules:
  Vg_glys := Vmax_glys * (G6Pg / (G6Pg + Km_G6P_glys)) * (1 - 1 / (1 + exp(- aGLY_inh * (1 * (GLYg - GLY_inh)))));

  // Variable initializations:
  Vmax_glys = 0.0001528;
  Km_G6P_glys = 0.5;
  G6Pg = ;
  GLYg = ;
  GLY_inh = ;
  aGLY_inh = ;
end

model cloutier_model_2009__Vg_glyp(Vg_glyp, unitstepSB2, stim, to, to_GLY, tend_GLY, sr_GLY, t1, delta_GLY, time_, GLYg, KO3)

  // Assignment Rules:
  Vg_glyp := Vmax_glyp * (GLYg / (GLYg + Km_GLY)) * deltaVt_GLY;
  deltaVt_GLY := 1 + stim * (delta_GLY * KO3 * (1 / (1 + exp(1 * - sr_GLY * (time_ - (to + to_GLY))))) * (1 - unitstepSB2));

  // Variable initializations:
  Vmax_glyp = 4.922e-5;
  Km_GLY = 1.0;
  unitstepSB2 = ;
  stim = ;
  to = ;
  to_GLY = ;
  tend_GLY = ;
  sr_GLY = ;
  t1 = ;
  delta_GLY = ;
  time_ = ;
  GLYg = ;
  KO3 = ;
end

model cloutier_model_2009__Fin_t(Fin_t, stim, to, tend, sr, t1, deltaf, time_)

  // Assignment Rules:
  Fin_t := CBF0 + (stim * CBF0 * deltaf * (1 / (1 + exp((1 * - sr) * (time_ - ((to + t1) - 3))))) - stim * CBF0 * deltaf * (1 / (1 + exp((1 * - sr) * (time_ - (to + tend + t1 + 3))))));

  // Variable initializations:
  CBF0 = 0.012;
  stim = ;
  to = ;
  tend = ;
  sr = ;
  t1 = ;
  deltaf = ;
  time_ = ;
end

model cloutier_model_2009__Fout_t(Fout_t, Fin_t, CBF0, Vv, Vv0, tv)

  // Assignment Rules:
  Fout_t := CBF0 * ((power(Vv / Vv0, 2) + tv * power(Vv / Vv0, -0.5) * (Fin_t / Vv0)) / (1 + CBF0 * tv * power(Vv / Vv0, -0.5) * (1 / Vv0)));

  // Variable initializations:
  Fin_t = ;
  CBF0 = ;
  Vv = ;
  Vv0 = ;
  tv = ;
end

model cloutier_model_2009__NADn(NADn, NADHn)

  // Assignment Rules:
  NADn := NADH_n_tot - NADHn;

  // Variable initializations:
  NADH_n_tot = 0.22;
  NADHn = ;
end

model cloutier_model_2009__NADg(NADg, NADHg)

  // Assignment Rules:
  NADg := NADH_g_tot - NADHg;

  // Variable initializations:
  NADH_g_tot = 0.22;
  NADHg = ;
end

model cloutier_model_2009__CRn(CRn, PCrn)

  // Assignment Rules:
  CRn := PCrn_tot - PCrn;

  // Variable initializations:
  PCrn_tot = 5.0;
  PCrn = ;
end

model cloutier_model_2009__CRg(CRg, PCrg)

  // Assignment Rules:
  CRg := PCrg_tot - PCrg;

  // Variable initializations:
  PCrg_tot = 5.0;
  PCrg = ;
end

model cloutier_model_2009__ADPn(ADPn, ATPn, ATPtot, qak)

  // Assignment Rules:
  ADPn := (ATPn / 2) * (- qak + root(power(qak, 2) + 4 * qak * (ATPtot / ATPn - 1)));

  // Variable initializations:
  ATPn = ;
  ATPtot = ;
  qak = ;
end

model cloutier_model_2009__ADPg(ADPg, ATPg, ATPtot, qak)

  // Assignment Rules:
  ADPg := (ATPg / 2) * (- qak + root(power(qak, 2) + 4 * qak * (ATPtot / ATPg - 1)));

  // Variable initializations:
  ATPg = ;
  ATPtot = ;
  qak = ;
end

model cloutier_model_2009__u_n(u_n, ATPn, ATPtot, qak)

  // Assignment Rules:
  u_n := power(qak, 2) + 4 * qak * (ATPtot / ATPn - 1);

  // Variable initializations:
  ATPn = ;
  ATPtot = ;
  qak = ;
end

model cloutier_model_2009__u_g(u_g, ATPg, ATPtot, qak)

  // Assignment Rules:
  u_g := power(qak, 2) + 4 * qak * (ATPtot / ATPg - 1);

  // Variable initializations:
  ATPg = ;
  ATPtot = ;
  qak = ;
end

model cloutier_model_2009__dAMP_dATPn(dAMP_dATPn, ATPn, ATPtot, u_n, qak)

  // Assignment Rules:
  dAMP_dATPn := (qak / 2 + qak * (ATPtot / (ATPn * root(u_n)))) - (1 + 0.5 * root(u_n));

  // Variable initializations:
  ATPn = ;
  ATPtot = ;
  u_n = ;
  qak = ;
end

model cloutier_model_2009__dAMP_dATPg(dAMP_dATPg, ATPg, ATPtot, u_g, qak)

  // Assignment Rules:
  dAMP_dATPg := (qak / 2 + qak * (ATPtot / (ATPg * root(u_g)))) - (1 + 0.5 * root(u_g));

  // Variable initializations:
  ATPg = ;
  ATPtot = ;
  u_g = ;
  qak = ;
end

model cloutier_model_2009__AMPn(AMPn, ATPtot, ATPn, ADPn)

  // Assignment Rules:
  AMPn := ATPtot - (ATPn + ADPn);

  // Variable initializations:
  ATPtot = ;
  ATPn = ;
  ADPn = ;
end

model cloutier_model_2009__AMPg(AMPg, ATPtot, ATPg, ADPg)

  // Assignment Rules:
  AMPg := ATPtot - (ATPg + ADPg);

  // Variable initializations:
  ATPtot = ;
  ATPg = ;
  ADPg = ;
end

model cloutier_model_2009__BOLD(Vv0, k1, k2, k3, Vv, dHb, dHb0)

  // Assignment Rules:
  BOLD := Vv0 * ((k1 + k2) * (1 - dHb / dHb0) - (k2 + k3) * (1 - Vv / Vv0));

  // Variable initializations:
  Vv0 = ;
  k1 = ;
  k2 = ;
  k3 = ;
  Vv = ;
  dHb = ;
  dHb0 = ;
end

model cloutier_model_2009__v_stim(v_stim, stim, to, tend, t_n_stim, v1_n, v2_n, time_)

  // Assignment Rules:
  v_stim := stim * (v1_n + v2_n * ((time_ - to) / t_n_stim) * exp(-((time_ - to) * (unitpulseSB / t_n_stim)))) * unitpulseSB;
  unitpulseSB := piecewise( 1 , ( geq(time_, to)) && (time_ <= to + tend ), 0 );

  // Variable initializations:
  stim = ;
  to = ;
  tend = ;
  t_n_stim = ;
  v1_n = ;
  v2_n = ;
  time_ = ;
end

model cloutier_model_2009__unitstepSB(unitstepSB, to, tend, time_)

  // Assignment Rules:
  unitstepSB := piecewise( 1 , geq(time_ - (tend + to), 0) , 0 );

  // Variable initializations:
  to = ;
  tend = ;
  time_ = ;
end

model cloutier_model_2009__model_parameters(Km_PYR, Km_ATP, Ki_ATP, Km_ADP, Km_O2, Km_GLC, Km_GLU, Km_G6P, Km_F6P_pgi, Km_F6P_pfk, Km_pump, nh_O2, Ko2, kpump, ATPtot, NAe, nH, O2a, CO2a, GLCa, LACa, nOP, NAero, Rng, Reg, Ren, Rcn, Rcg, Rce, Sm_n, Vm, RT, F, Vn, G6P_inh_hk, aG6P_inh_hk, rATP_mito, aATP_mito, HbOP, Sm_g, Vg, KO1, Vc, R_GLU_NA, KO2, KO3, GLY_inh, aGLY_inh, CBF0, Vv0, tv, qak, k1, k2, k3, dHb0, stim, to, tend, v1_n, v2_n, t_n_stim, sr, t1, delta_GLY, deltaf, tend_GLY, to_GLY, sr_GLY)

  // Variable initializations:
  Km_PYR = 0.0632;
  Km_ATP = 0.01532;
  Ki_ATP = 0.7595;
  Km_ADP = 0.00107;
  Km_O2 = 0.0029658;
  Km_GLC = 0.105;
  Km_GLU = 0.05;
  Km_G6P = 0.5;
  Km_F6P_pgi = 0.06;
  Km_F6P_pfk = 0.18;
  Km_pump = 0.4243;
  nh_O2 = 2.7;
  Ko2 = 0.089733;
  kpump = 3.17e-7;
  ATPtot = 2.379;
  NAe = 150.0;
  nH = 4.0;
  O2a = 8.34;
  CO2a = 1.2;
  GLCa = 4.8;
  LACa = 0.313;
  nOP = 15.0;
  NAero = 3.0;
  Rng = 1.8;
  Reg = 0.8;
  Ren = 0.4444444444444444;
  Rcn = 0.01222;
  Rcg = 0.022;
  Rce = 0.0275;
  Sm_n = 40500;
  Vm = -70;
  RT = 2577340;
  F = 96500;
  Vn = 0.45;
  G6P_inh_hk = 0.6;
  aG6P_inh_hk = 20.0;
  rATP_mito = 20.0;
  aATP_mito = 5.0;
  HbOP = 8.6;
  Sm_g = 10500;
  Vg = 0.25;
  KO1 = 1.0;
  Vc = 0.0055;
  R_GLU_NA = 0.075;
  KO2 = 1;
  KO3 = 1;
  GLY_inh = 4.2;
  aGLY_inh = 20.0;
  CBF0 = 0.012;
  Vv0 = 0.0236;
  tv = 35.0;
  qak = 0.92;
  k1 = 2.22;
  k2 = 0.46;
  k3 = 0.43;
  dHb0 = 0.064;
  stim = 1;
  to = 200;
  tend = 300;
  v1_n = 0.041;
  v2_n = 2.55;
  t_n_stim = 2;
  sr = 4.59186;
  t1 = 2;
  delta_GLY = 62;
  deltaf = 0.42;
  tend_GLY = 440;
  to_GLY = 83;
  sr_GLY = 4;
end

model cloutier_model_2009__unitstepSB2(unitstepSB2, tend, time_, to, tend_GLY, to_GLY)

  // Assignment Rules:
  unitstepSB2 := piecewise( 1 , geq(time_ - (tend_GLY + to + to_GLY), 0) , 0 );

  // Variable initializations:
  tend = ;
  time_ = ;
  to = ;
  tend_GLY = ;
  to_GLY = ;
end

model *cloutier_model_2009____main()

  // Sub-modules, and any changes to those submodules:
  environment: cloutier_model_2009__environment(time_);
  NAn_: cloutier_model_2009__NAn(NAn_0, Vn_leak_Na0, Vn_pump0, Vn_stim0, time_);
  GLCn: cloutier_model_2009__GLCn(GLCn0, V_en_GLC0, Vn_hk0, time_);
  G6Pn: cloutier_model_2009__G6Pn(G6Pn0, Vn_pgi0, Vn_hk0, time_);
  F6Pn: cloutier_model_2009__F6Pn(F6Pn0, Vn_pgi0, Vn_pfk0, time_);
  GAPn: cloutier_model_2009__GAPn(GAPn0, Vn_pgk0, Vn_pfk0, time_);
  PEPn: cloutier_model_2009__PEPn(PEPn0, Vn_pgk0, Vn_pk0, time_);
  PYRn: cloutier_model_2009__PYRn(PYRn0, Vn_pk0, Vn_ldh0, Vn_mito0, time_);
  LACn: cloutier_model_2009__LACn(LACn0, Vn_ldh0, Vne_LAC0, time_);
  NADHn: cloutier_model_2009__NADHn(NADHn0, Vn_pgk0, Vn_ldh0, Vn_mito0, time_);
  ATPn: cloutier_model_2009__ATPn(ATPn0, nOP, Vn_hk0, Vn_pfk0, Vn_pgk0, Vn_pk0, Vn_ATPase0, Vn_pump0, Vn_mito0, Vn_ck0, dAMP_dATPn0, time_);
  PCrn: cloutier_model_2009__PCrn(PCrn0, Vn_ck0, time_);
  O2n: cloutier_model_2009__O2n(O2n0, NAero, Vn_mito0, Vcn_O20, time_);
  GLUn: cloutier_model_2009__GLUn(GLUn0, Rng, Vg_gs0, Vn_stim_GLU0, time_);
  NAg: cloutier_model_2009__NAg(NAg0, Vg_leak_Na0, Vg_pump0, Veg_GLU0, time_);
  GLCg: cloutier_model_2009__GLCg(GLCg0, Vcg_GLC0, Veg_GLC0, Vg_hk0, time_);
  G6Pg: cloutier_model_2009__G6Pg(G6Pg0, Vg_pgi0, Vg_glys0, Vg_glyp0, Vg_hk0, time_);
  F6Pg: cloutier_model_2009__F6Pg(F6Pg0, Vg_pgi0, Vg_pfk0, time_);
  GAPg: cloutier_model_2009__GAPg(GAPg0, Vg_pgk0, Vg_pfk0, time_);
  PEPg: cloutier_model_2009__PEPg(PEPg0, Vg_pgk0, Vg_pk0, time_);
  PYRg: cloutier_model_2009__PYRg(PYRg0, Vg_ldh0, Vg_pk0, Vg_mito0, time_);
  LACg: cloutier_model_2009__LACg(LACg0, Vg_ldh0, Vge_LAC0, Vgc_LAC0, time_);
  NADHg: cloutier_model_2009__NADHg(NADHg0, Vg_pgk0, Vg_ldh0, Vg_mito0, time_);
  ATPg: cloutier_model_2009__ATPg(ATPg0, nOP, Vg_hk0, Vg_pfk0, Vg_pgk0, Vg_pk0, Vg_ATPase0, Vg_pump0, Vg_mito0, Vg_ck0, Vg_gs0, dAMP_dATPg0, time_);
  PCrg: cloutier_model_2009__PCrg(PCrg0, Vg_ck0, time_);
  O2g: cloutier_model_2009__O2g(O2g0, NAero, Vg_mito0, Vcg_O20, time_);
  GLYg: cloutier_model_2009__GLYg(GLYg0, Vg_glys0, Vg_glyp0, time_);
  GLUg: cloutier_model_2009__GLUg(GLUg0, Veg_GLU0, Vg_gs0, time_);
  GLCe: cloutier_model_2009__GLCe(GLCe0, Reg, Ren, Vce_GLC0, Veg_GLC0, V_en_GLC0, time_);
  LACe: cloutier_model_2009__LACe(LACe0, Reg, Ren, Vec_LAC0, Vne_LAC0, Vge_LAC0, time_);
  GLUe: cloutier_model_2009__GLUe(GLUe0, Reg, Ren, Vn_stim_GLU0, Veg_GLU0, time_);
  O2c: cloutier_model_2009__O2c(O2c0, Rcn, Rcg, Vc_O20, Vcn_O20, Vcg_O20, time_);
  GLCc: cloutier_model_2009__GLCc(GLCc0, Rce, Rcg, Vc_GLC0, Vce_GLC0, Vcg_GLC0, time_);
  LACc: cloutier_model_2009__LACc(LACc0, Rce, Rcg, Vc_LAC0, Vec_LAC0, Vgc_LAC0, time_);
  CO2c: cloutier_model_2009__CO2c(CO2c0, Rcn, Rcg, Vnc_CO20, Vc_CO20, Vgc_CO20, time_);
  Vv: cloutier_model_2009__Vv(Vv0, Fin_t0, Fout_t0, time_);
  dHb: cloutier_model_2009__dHb(dHb1, Fin_t0, Fout_t0, Vv0, O2a, O2c0, time_);
  Vn_leak_Na: cloutier_model_2009__Vn_leak_Na(Vn_leak_Na0, Sm_n, Vm, Vn, RT, F, NAe, NAn_0);
  Vn_pump: cloutier_model_2009__Vn_pump(Vn_pump0, kpump, Km_pump, Sm_n, Vn, ATPn0, NAn_0);
  Vn_stim: cloutier_model_2009__Vn_stim(Vn_stim0, v_stim0);
  V_en_GLC: cloutier_model_2009__V_en_GLC(V_en_GLC0, GLCe0, GLCn0);
  Vn_hk: cloutier_model_2009__Vn_hk(Vn_hk0, Km_GLC, GLCn0, ATPn0, G6Pn0, G6P_inh_hk, aG6P_inh_hk);
  Vn_pgi: cloutier_model_2009__Vn_pgi(Vn_pgi0, Km_G6P, Km_F6P_pgi, F6Pn0, G6Pn0);
  Vn_pfk: cloutier_model_2009__Vn_pfk(Vn_pfk0, Km_F6P_pfk, Ki_ATP, nH, ATPn0, F6Pn0);
  Vn_pgk: cloutier_model_2009__Vn_pgk(Vn_pgk0, ADPn0, GAPn0, NADn0, NADHn0);
  Vn_pk: cloutier_model_2009__Vn_pk(Vn_pk0, PEPn0, ADPn0);
  Vn_ldh: cloutier_model_2009__Vn_ldh(Vn_ldh0, PYRn0, NADHn0, LACn0, NADn0);
  Vn_mito: cloutier_model_2009__Vn_mito(Vn_mito0, Km_O2, Km_ADP, Km_PYR, O2n0, ADPn0, PYRn0, ATPn0, rATP_mito, aATP_mito);
  Vne_LAC: cloutier_model_2009__Vne_LAC(Vne_LAC0, LACn0, LACe0);
  Vn_ATPase: cloutier_model_2009__Vn_ATPase(Vn_ATPase0, ATPn0);
  Vn_ck: cloutier_model_2009__Vn_ck(Vn_ck0, PCrn0, ADPn0, CRn0, ATPn0);
  Vcn_O2: cloutier_model_2009__Vcn_O2(Vcn_O20, Ko2, Vn, HbOP, O2c0, O2n0);
  Vg_leak_Na: cloutier_model_2009__Vg_leak_Na(Vg_leak_Na0, Sm_g, Vm, Vg, RT, F, NAe, NAg0);
  Vg_pump: cloutier_model_2009__Vg_pump(Vg_pump0, kpump, Km_pump, Sm_g, Vg, ATPg0, NAg0);
  Veg_GLC: cloutier_model_2009__Veg_GLC(Veg_GLC0, GLCg0, GLCe0, KO1);
  Vcg_GLC: cloutier_model_2009__Vcg_GLC(Vcg_GLC0, GLCg0, GLCc0);
  Vg_hk: cloutier_model_2009__Vg_hk(Vg_hk0, Km_GLC, GLCg0, ATPg0, G6Pg0, G6P_inh_hk, aG6P_inh_hk);
  Vg_pgi: cloutier_model_2009__Vg_pgi(Vg_pgi0, Km_G6P, Km_F6P_pgi, F6Pg0, G6Pg0);
  Vg_pfk: cloutier_model_2009__Vg_pfk(Vg_pfk0, Km_F6P_pfk, Ki_ATP, nH, ATPg0, F6Pg0);
  Vg_pgk: cloutier_model_2009__Vg_pgk(Vg_pgk0, ADPg0, GAPg0, NADg0, NADHg0);
  Vg_pk: cloutier_model_2009__Vg_pk(Vg_pk0, PEPg0, ADPg0);
  Vg_ldh: cloutier_model_2009__Vg_ldh(Vg_ldh0, PYRg0, NADHg0, LACg0, NADg0);
  Vg_mito: cloutier_model_2009__Vg_mito(Vg_mito0, Km_O2, Km_ADP, Km_PYR, O2g0, ADPg0, PYRg0, ATPg0, rATP_mito, aATP_mito);
  Vge_LAC: cloutier_model_2009__Vge_LAC(Vge_LAC0, LACg0, LACe0);
  Vgc_LAC: cloutier_model_2009__Vgc_LAC(Vgc_LAC0, LACg0, LACc0);
  Vg_ATPase: cloutier_model_2009__Vg_ATPase(Vg_ATPase0, ATPg0);
  Vg_ck: cloutier_model_2009__Vg_ck(Vg_ck0, PCrg0, ADPg0, CRg0, ATPg0);
  Vcg_O2: cloutier_model_2009__Vcg_O2(Vcg_O20, Ko2, nh_O2, Vg, HbOP, O2c0, O2g0);
  Vc_O2: cloutier_model_2009__Vc_O2(Vc_O20, Fin_t0, Vc, O2a, O2c0);
  Vc_GLC: cloutier_model_2009__Vc_GLC(Vc_GLC0, Fin_t0, Vc, GLCa, GLCc0);
  Vce_GLC: cloutier_model_2009__Vce_GLC(Vce_GLC0, GLCe0, GLCc0);
  Vc_LAC: cloutier_model_2009__Vc_LAC(Vc_LAC0, Fin_t0, Vc, LACa, LACc0);
  Vec_LAC: cloutier_model_2009__Vec_LAC(Vec_LAC0, LACe0, LACc0);
  Vnc_CO2: cloutier_model_2009__Vnc_CO2(Vnc_CO20, Vn_mito0);
  Vgc_CO2: cloutier_model_2009__Vgc_CO2(Vgc_CO20, Vg_mito0);
  Vn_stim_GLU: cloutier_model_2009__Vn_stim_GLU(Vn_stim_GLU0, R_GLU_NA, Km_GLU, GLUn0, KO2, Vn_stim0);
  Vg_gs: cloutier_model_2009__Vg_gs(Vg_gs0, Km_GLU, Km_ATP, GLUg0, ATPg0);
  Veg_GLU: cloutier_model_2009__Veg_GLU(Veg_GLU0, Km_GLU, GLUe0);
  Vc_CO2: cloutier_model_2009__Vc_CO2(Vc_CO20, Fin_t0, Vc, CO2a, CO2c0);
  Vg_glys: cloutier_model_2009__Vg_glys(Vg_glys0, G6Pg0, GLYg0, GLY_inh, aGLY_inh);
  Vg_glyp: cloutier_model_2009__Vg_glyp(Vg_glyp0, unitstepSB20, stim, to, to_GLY, tend_GLY, sr_GLY, t1, delta_GLY, time_, GLYg0, KO3);
  Fin_t: cloutier_model_2009__Fin_t(Fin_t0, stim, to, tend, sr, t1, deltaf, time_);
  Fout_t: cloutier_model_2009__Fout_t(Fout_t0, Fin_t0, CBF0, Vv0, Vv00, tv);
  NADn: cloutier_model_2009__NADn(NADn0, NADHn0);
  NADg: cloutier_model_2009__NADg(NADg0, NADHg0);
  CRn: cloutier_model_2009__CRn(CRn0, PCrn0);
  CRg: cloutier_model_2009__CRg(CRg0, PCrg0);
  ADPn: cloutier_model_2009__ADPn(ADPn0, ATPn0, ATPtot, qak);
  ADPg: cloutier_model_2009__ADPg(ADPg0, ATPg0, ATPtot, qak);
  u_n: cloutier_model_2009__u_n(u_n0, ATPn0, ATPtot, qak);
  u_g: cloutier_model_2009__u_g(u_g0, ATPg0, ATPtot, qak);
  dAMP_dATPn: cloutier_model_2009__dAMP_dATPn(dAMP_dATPn0, ATPn0, ATPtot, u_n0, qak);
  dAMP_dATPg: cloutier_model_2009__dAMP_dATPg(dAMP_dATPg0, ATPg0, ATPtot, u_g0, qak);
  AMPn: cloutier_model_2009__AMPn(AMPn0, ATPtot, ATPn0, ADPn0);
  AMPg: cloutier_model_2009__AMPg(AMPg0, ATPtot, ATPg0, ADPg0);
  BOLD: cloutier_model_2009__BOLD(Vv00, k1, k2, k3, Vv0, dHb1, dHb0);
  v_stim: cloutier_model_2009__v_stim(v_stim0, stim, to, tend, t_n_stim, v1_n, v2_n, time_);
  unitstepSB: cloutier_model_2009__unitstepSB(unitstepSB0, to, tend, time_);
  model_parameters: cloutier_model_2009__model_parameters(Km_PYR, Km_ATP, Ki_ATP, Km_ADP, Km_O2, Km_GLC, Km_GLU, Km_G6P, Km_F6P_pgi, Km_F6P_pfk, Km_pump, nh_O2, Ko2, kpump, ATPtot, NAe, nH, O2a, CO2a, GLCa, LACa, nOP, NAero, Rng, Reg, Ren, Rcn, Rcg, Rce, Sm_n, Vm, RT, F, Vn, G6P_inh_hk, aG6P_inh_hk, rATP_mito, aATP_mito, HbOP, Sm_g, Vg, KO1, Vc, R_GLU_NA, KO2, KO3, GLY_inh, aGLY_inh, CBF0, Vv00, tv, qak, k1, k2, k3, dHb0, stim, to, tend, v1_n, v2_n, t_n_stim, sr, t1, delta_GLY, deltaf, tend_GLY, to_GLY, sr_GLY);
  unitstepSB2: cloutier_model_2009__unitstepSB2(unitstepSB20, tend, time_, to, tend_GLY, to_GLY);
end