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

  // Variable initializations:
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__Asc(Asc, pG, pDA, T_Asc, time_)

  // Rate Rules:
  Asc' = gamma_Asc * Asc * (C_Asc - Asc) - (pG + pDA) * T_Asc * Asc;

  // Variable initializations:
  Asc = ;
  C_Asc = ;
  gamma_Asc = ;
  pG = ;
  pDA = ;
  T_Asc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__Csc(Csc_, GtnC, pE, pDC, pP, T_Csc, time_)

  // Rate Rules:
  Csc_' = (gamma_Csc * Csc_ * (C_Csc - Csc_) + ((gmax * power(GtnC, 2.0)) / (power(GtnC, 2.0) + power(alpha_Csc, 2.0))) * Csc_) - (pE + pDC + pP) * T_Csc * Csc_;

  // Variable initializations:
  Csc_ = ;
  C_Csc = ;
  gamma_Csc = ;
  alpha_Csc = ;
  gmax = ;
  GtnC = ;
  pE = ;
  pDC = ;
  pP = ;
  T_Csc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__G(G, AC, Fd, alpha_HA, pG, gamma_Fdmax, alpha_Fd, Asc, T_Asc, time_)

  // Rate Rules:
  G' = (pG * T_Asc * Asc + k_gmax * (1.0 - power(AC, 2.0) / (power(AC, 2.0) + power(alpha_HA, 2.0))) * G) - (gamma_Fdmax * (1.0 - power(Fd, 2.0) / (power(Fd, 2.0) + power(alpha_Fd, 2.0))) * G + gamma_Gc * G);

  // Variable initializations:
  G = ;
  k_gmax = ;
  gamma_Gc = ;
  AC = ;
  Fd = ;
  alpha_HA = ;
  pG = ;
  gamma_Fdmax = ;
  alpha_Fd = ;
  Asc = ;
  T_Asc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__DA(DA, AC, Fd, pDA, alpha_HA, gamma_Fdmax, alpha_Fd, Asc, T_Asc, time_)

  // Rate Rules:
  DA' = (pDA * T_Asc * Asc + ((k_dmax * power(AC, 2.0)) / (power(AC, 2.0) + power(alpha_HA, 2.0))) * DA + gamma_Fdmax * (1.0 - power(Fd, 2.0) / (power(Fd, 2.0) + power(alpha_Fd, 2.0))) * DA) - gamma_DA * DA;

  // Variable initializations:
  DA = ;
  k_dmax = ;
  gamma_DA = ;
  AC = ;
  Fd = ;
  pDA = ;
  alpha_HA = ;
  gamma_Fdmax = ;
  alpha_Fd = ;
  Asc = ;
  T_Asc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__DC(DC, pDC, Csc_, T_Asc, time_)

  // Rate Rules:
  DC' = pDC * T_Asc * Csc_ - gamma_DC * DC;

  // Variable initializations:
  DC = ;
  gamma_DC = ;
  pDC = ;
  Csc_ = ;
  T_Asc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__E(E, pE, GtnC, Csc_, T_Csc, time_)

  // Rate Rules:
  E' = (pE * T_Csc * Csc_ - gamma_E * E) + ((k_Emax * power(GtnC, 2.0)) / (power(GtnC, 2.0) + power(alpha_E, 2.0))) * E;

  // Variable initializations:
  E = ;
  k_Emax = ;
  gamma_E = ;
  alpha_E = ;
  pE = ;
  GtnC = ;
  Csc_ = ;
  T_Csc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__P(P, pP, Csc_, T_Csc, time_)

  // Rate Rules:
  P' = pP * T_Csc * Csc_ - gamma_P * P;

  // Variable initializations:
  P = ;
  gamma_P = ;
  pP = ;
  Csc_ = ;
  T_Csc = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__feedback_mechanisms(pG, pDA, pP, pE, pDC, G, DA, DC, E, P)

  // Assignment Rules:
  pG := exp((- f * power(G, 2.0)) / power(N_, 2.0));
  pDA := exp((- f * power(DA, 2.0)) / power(N_, 2.0));
  pP := exp((- f * power(P, 2.0)) / power(N_, 2.0));
  pE := exp((- f * power(E, 2.0)) / power(N_, 2.0));
  pDC := exp((- f * power(DC, 2.0)) / power(N_, 2.0));

  // Variable initializations:
  f = ;
  N_ = ;
  G = ;
  DA = ;
  DC = ;
  E = ;
  P = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__GtnA(GtnA, alpha_Fd, beta_G, AC, G, SA, NC, NE, Fd, time_)

  // Rate Rules:
  GtnA' = -(G * ((KNG1 * NE) / ((NE + alpha_NG1) * (1.0 + SA / kSG) * (1.0 + power(AC, 2.0) / (power(AC, 2.0) + power(kAG, 2.0)))) + (KNG2 * NC) / ((NC + alpha_NG2) * (1.0 + SA / kSG) * (1.0 + power(AC, 2.0) / (power(AC, 2.0) + power(kAG, 2.0)))) + (KFG * Fd) / ((Fd + alpha_Fd) * (1.0 + SA / kSG) * (1.0 + power(AC, 2.0) / (power(AC, 2.0) + power(kAG, 2.0))))) * ((kG + beta_G) * GtnA));

  // Variable initializations:
  GtnA = ;
  KNG1 = 6.28E-17;
  KNG2 = 8.75E-17;
  KFG = ;
  kAG = ;
  kSG = 9.0E-11;
  kG = ;
  alpha_NG1 = 1.0E-10;
  alpha_NG2 = 1.0E-10;
  alpha_Fd = ;
  beta_G = ;
  AC = ;
  G = ;
  SA = ;
  NC = ;
  NE = ;
  Fd = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__GtnC(GtnC, beta_G, GtnA, time_)

  // Rate Rules:
  GtnC' = beta_G * GtnA - kappa_G * GtnC;

  // Variable initializations:
  GtnC = ;
  kappa_G = ;
  beta_G = ;
  GtnA = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__SA(SA, kappa_S, alpha_NS, kNS, kSS, KGS, NC, NE, AA, DA, time_)

  // Rate Rules:
  SA' = -(DA * ((KAS * AA) / ((AA + alpha_AS) * (1.0 + SA / kSS) * (1.0 + NC / kNS)) + (KGS * NE) / ((NE + alpha_NS) * (1.0 + SA / kSS) * (1.0 + NC / kNS))) * (kappa_S * SA));

  // Variable initializations:
  SA = ;
  KAS = 8.04E-15;
  alpha_AS = 0.05;
  kappa_S = ;
  alpha_NS = ;
  kNS = ;
  kSS = ;
  KGS = ;
  NC = ;
  NE = ;
  AA = ;
  DA = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__SC(SC, kappa_S, alpha_NS, kNS, kSS, KGS, NC, NE, DC, GtnC, time_)

  // Rate Rules:
  SC' = -(DC * ((KNS * NE) / ((NE + alpha_NS) * (1.0 + SC / kSS) * (1.0 + NC / kNS)) + (KGS * GtnC) / ((GtnC + alpha_GS) * (1.0 + SC / kSS) * (1.0 + NC / kNS))) * (kappa_S * SC));

  // Variable initializations:
  SC = ;
  KNS = 1.14E-15;
  alpha_GS = 5.20E-12;
  kappa_S = ;
  alpha_NS = ;
  kNS = ;
  kSS = ;
  KGS = ;
  NC = ;
  NE = ;
  DC = ;
  GtnC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__HC(HC, E, SC, NE, NC, GtnC, time_)

  // Rate Rules:
  HC' = -(E * ((KNH * NE) / ((NE + alpha_NH) * (1.0 + SC / kSH)) + (KGH * GtnC) / ((GtnC + alpha_GH) * (1.0 + SC / kSH))) * (kappa_H * HC));

  // Variable initializations:
  HC = ;
  KNH = 7.59E-16;
  KGH = 7.77E-16;
  kSH = 9.0E-10;
  alpha_GH = 3.0E-10;
  alpha_NH = 3.25E-8;
  kappa_H = ;
  E = ;
  SC = ;
  NE = ;
  NC = ;
  GtnC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__AC(AC, P, beta_A, BC, hb, alpha_HA, Fd, alpha_Fd, NC, SC, HC, GtnC, time_)

  // Rate Rules:
  AC' = -(P * ((KNA * NC) / ((NC + alpha_NA) * (1.0 + SC / kSA)) + HC / (HC + alpha_HA) + (KGA * GtnC) / ((GtnC + alpha_GA) * (1.0 + SC / kSA)) + (KHA * HC) / ((HC + alpha_HA) * (1.0 + SC / kSA))) * (hb * AC * BC + ((k_Fdmax * Fd) / (Fd + alpha_Fd)) * AC + beta_A * AC));

  // Variable initializations:
  AC = ;
  KNA = 2.33E-11;
  KHA = 7.96E-10;
  KGA = 4.98E-11;
  kSA = 9.0E-10;
  alpha_NA = 5.0E-6;
  alpha_GA = 1.8E-10;
  k_Fdmax = ;
  P = ;
  beta_A = ;
  BC = ;
  hb = ;
  alpha_HA = ;
  Fd = ;
  alpha_Fd = ;
  NC = ;
  SC = ;
  HC = ;
  GtnC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__AA(AA, beta_A, AC, time_)

  // Rate Rules:
  AA' = beta_A * AC - kappa_A * AA;

  // Variable initializations:
  AA = ;
  kappa_A = ;
  beta_A = ;
  AC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__BC(BC, alpha_NB, AC, hb, NC, time_)

  // Rate Rules:
  BC' = (k_BCmax * NC) / (NC + alpha_NB) - (hb * AC * BC + beta_B * BC);

  // Variable initializations:
  BC = ;
  k_BCmax = ;
  beta_B = ;
  alpha_NB = ;
  AC = ;
  hb = ;
  NC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__BA(BA, alpha_NB, AA, hb, NC, time_)

  // Rate Rules:
  BA' = (k_BAmax * NC) / (NC + alpha_NB) - (hb * AA * BA + kappa_B * BA);

  // Variable initializations:
  BA = ;
  k_BAmax = ;
  kappa_B = ;
  alpha_NB = ;
  AA = ;
  hb = ;
  NC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__NC(NC, kFd, Fd, AC, time_)

  // Rate Rules:
  NC' = ((Nmax1 * Fd) / ((Fd + kFd) * (1.0 + power(AC, 2.0) / (power(AC, 2.0) + power(kAN1, 2.0)))) - kappa_NC * NC) + Bas1;

  // Variable initializations:
  NC = ;
  Nmax1 = ;
  kAN1 = ;
  kappa_NC = ;
  Bas1 = ;
  kFd = ;
  Fd = ;
  AC = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__NE(NE, kFd, AC, Fd, time_)

  // Rate Rules:
  NE' = ((Nmax2 * Fd) / ((Fd + power(kFd, 2.0)) * (1.0 + power(AC, 2.0) / (power(AC, 2.0) + power(kAN2, 2.0)))) - kappa_NE * NE) + Bas2;

  // Variable initializations:
  NE = ;
  Nmax2 = ;
  kAN2 = ;
  kappa_NE = ;
  Bas2 = ;
  kFd = ;
  AC = ;
  Fd = ;
  time_ = ;
end

model joseph_zavros_merchant_kirschner_2003_version01__model_parameters(beta_G, beta_A, kFd, Fd, gamma_Fdmax, alpha_Fd, kappa_S, alpha_NS, alpha_NB, alpha_HA, kNS, kSS, KGS, T_Asc, T_Csc, hb)

  // Variable initializations:
  beta_G = 1.5;
  beta_A = 2.72;
  kFd = ;
  Fd = ;
  gamma_Fdmax = ;
  alpha_Fd = ;
  kappa_S = ;
  alpha_NS = 6.28E-7;
  alpha_NB = ;
  alpha_HA = 2.0E-8;
  kNS = 1.0E-9;
  kSS = 9.0E-11;
  KGS = 2.54E-18;
  T_Asc = ;
  T_Csc = ;
  hb = ;
end

model *joseph_zavros_merchant_kirschner_2003_version01____main()

  // Sub-modules, and any changes to those submodules:
  environment: joseph_zavros_merchant_kirschner_2003_version01__environment(time_);
  Asc: joseph_zavros_merchant_kirschner_2003_version01__Asc(Asc0, pG, pDA, T_Asc, time_);
  Csc_: joseph_zavros_merchant_kirschner_2003_version01__Csc(Csc_0, GtnC0, pE, pDC, pP, T_Csc, time_);
  G: joseph_zavros_merchant_kirschner_2003_version01__G(G0, AC0, Fd, alpha_HA, pG, gamma_Fdmax, alpha_Fd, Asc0, T_Asc, time_);
  DA: joseph_zavros_merchant_kirschner_2003_version01__DA(DA0, AC0, Fd, pDA, alpha_HA, gamma_Fdmax, alpha_Fd, Asc0, T_Asc, time_);
  DC: joseph_zavros_merchant_kirschner_2003_version01__DC(DC0, pDC, Csc_0, T_Asc, time_);
  E: joseph_zavros_merchant_kirschner_2003_version01__E(E0, pE, GtnC0, Csc_0, T_Csc, time_);
  P: joseph_zavros_merchant_kirschner_2003_version01__P(P0, pP, Csc_0, T_Csc, time_);
  feedback_mechanisms: joseph_zavros_merchant_kirschner_2003_version01__feedback_mechanisms(pG, pDA, pP, pE, pDC, G0, DA0, DC0, E0, P0);
  GtnA: joseph_zavros_merchant_kirschner_2003_version01__GtnA(GtnA0, alpha_Fd, beta_G, AC0, G0, SA0, NC0, NE0, Fd, time_);
  GtnC: joseph_zavros_merchant_kirschner_2003_version01__GtnC(GtnC0, beta_G, GtnA0, time_);
  SA: joseph_zavros_merchant_kirschner_2003_version01__SA(SA0, kappa_S, alpha_NS, kNS, kSS, KGS, NC0, NE0, AA0, DA0, time_);
  SC: joseph_zavros_merchant_kirschner_2003_version01__SC(SC0, kappa_S, alpha_NS, kNS, kSS, KGS, NC0, NE0, DC0, GtnC0, time_);
  HC: joseph_zavros_merchant_kirschner_2003_version01__HC(HC0, E0, SC0, NE0, NC0, GtnC0, time_);
  AC: joseph_zavros_merchant_kirschner_2003_version01__AC(AC0, P0, beta_A, BC0, hb, alpha_HA, Fd, alpha_Fd, NC0, SC0, HC0, GtnC0, time_);
  AA: joseph_zavros_merchant_kirschner_2003_version01__AA(AA0, beta_A, AC0, time_);
  BC: joseph_zavros_merchant_kirschner_2003_version01__BC(BC0, alpha_NB, AC0, hb, NC0, time_);
  BA: joseph_zavros_merchant_kirschner_2003_version01__BA(BA0, alpha_NB, AA0, hb, NC0, time_);
  NC: joseph_zavros_merchant_kirschner_2003_version01__NC(NC0, kFd, Fd, AC0, time_);
  NE: joseph_zavros_merchant_kirschner_2003_version01__NE(NE0, kFd, AC0, Fd, time_);
  model_parameters: joseph_zavros_merchant_kirschner_2003_version01__model_parameters(beta_G, beta_A, kFd, Fd, gamma_Fdmax, alpha_Fd, kappa_S, alpha_NS, alpha_NB, alpha_HA, kNS, kSS, KGS, T_Asc, T_Csc, hb);
end