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

  // Variable initializations:
  time_ = ;
end

model kesmir_1999__B0(B0, pr, mu, rho, delta_B, CT_star, time_)

  // Rate Rules:
  B0' = pr * mu * CT_star - (rho * B0 + delta_B * B0);

  // Variable initializations:
  B0 = 3;
  pr = ;
  mu = ;
  rho = ;
  delta_B = ;
  CT_star = ;
  time_ = ;
end

model kesmir_1999__B1(B1, alpha_B, rho, delta_B, B0, time_)

  // Rate Rules:
  B1' = rho * (1 + alpha_B) * B0 - (rho * B1 + delta_B * B1);

  // Variable initializations:
  B1 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B0 = ;
  time_ = ;
end

model kesmir_1999__B2(B2, alpha_B, rho, delta_B, B1, time_)

  // Rate Rules:
  B2' = rho * (1 + alpha_B) * B1 - (rho * B2 + delta_B * B2);

  // Variable initializations:
  B2 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B1 = ;
  time_ = ;
end

model kesmir_1999__B3(B3, alpha_B, rho, delta_B, B2, time_)

  // Rate Rules:
  B3' = rho * (1 + alpha_B) * B2 - (rho * B3 + delta_B * B3);

  // Variable initializations:
  B3 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B2 = ;
  time_ = ;
end

model kesmir_1999__B4(B4, alpha_B, rho, delta_B, B3, time_)

  // Rate Rules:
  B4' = rho * (1 + alpha_B) * B3 - (rho * B4 + delta_B * B4);

  // Variable initializations:
  B4 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B3 = ;
  time_ = ;
end

model kesmir_1999__B5(B5, alpha_B, rho, delta_B, B4, time_)

  // Rate Rules:
  B5' = rho * (1 + alpha_B) * B4 - (rho * B5 + delta_B * B5);

  // Variable initializations:
  B5 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B4 = ;
  time_ = ;
end

model kesmir_1999__B6(B6, alpha_B, rho, delta_B, B5, time_)

  // Rate Rules:
  B6' = rho * (1 + alpha_B) * B5 - (rho * B6 + delta_B * B6);

  // Variable initializations:
  B6 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B5 = ;
  time_ = ;
end

model kesmir_1999__B7(B7, alpha_B, rho, delta_B, B6, time_)

  // Rate Rules:
  B7' = rho * (1 + alpha_B) * B6 - (rho * B7 + delta_B * B7);

  // Variable initializations:
  B7 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B6 = ;
  time_ = ;
end

model kesmir_1999__B8(B8, alpha_B, rho, delta_B, B7, time_)

  // Rate Rules:
  B8' = rho * (1 + alpha_B) * B7 - (rho * B8 + delta_B * B8);

  // Variable initializations:
  B8 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B7 = ;
  time_ = ;
end

model kesmir_1999__B9(B9, alpha_B, rho, delta_B, B8, time_)

  // Rate Rules:
  B9' = rho * (1 + alpha_B) * B8 - (rho * B9 + delta_B * B9);

  // Variable initializations:
  B9 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B8 = ;
  time_ = ;
end

model kesmir_1999__B10(B10, alpha_B, rho, delta_B, B9, time_)

  // Rate Rules:
  B10' = rho * (1 + alpha_B) * B9 - (rho * B10 + delta_B * B10);

  // Variable initializations:
  B10 = 0;
  alpha_B = ;
  rho = ;
  delta_B = ;
  B9 = ;
  time_ = ;
end

model kesmir_1999__centroblasts_sum(B_sum, B10, B9, B8, B7, B6, B5, B4, B3, B2, B1, B0)

  // Assignment Rules:
  B_sum := B1 + B1 + B3 + B4 + B5 + B6 + B6 + B7 + B8 + B9 + B10;

  // Variable initializations:
  B10 = ;
  B9 = ;
  B8 = ;
  B7 = ;
  B6 = ;
  B5 = ;
  B4 = ;
  B3 = ;
  B2 = ;
  B1 = ;
  B0 = ;
end

model kesmir_1999__C(C, rho, mu, B10, time_)

  // Rate Rules:
  C' = d * rho * B10 * 1 - mu * C;

  // Variable initializations:
  C = 0;
  d = 2;
  rho = ;
  mu = ;
  B10 = ;
  time_ = ;
end

model kesmir_1999__C_star(C_star, mu, CA, time_)

  // Rate Rules:
  C_star' = mu * CA - mu * C_star;

  // Variable initializations:
  C_star = 0;
  mu = ;
  CA = ;
  time_ = ;
end

model kesmir_1999__centrocytes_sum(C_starsum, C, C_star, CT_star, CA)

  // Assignment Rules:
  C_starsum := C + C_star;

  // Variable initializations:
  C = ;
  C_star = ;
  CT_star = ;
  CA = ;
end

model kesmir_1999__M(M, pr, mu, CT_star, time_)

  // Rate Rules:
  M' = (1 - pr) * mu * CT_star;

  // Variable initializations:
  M = 0;
  pr = ;
  mu = ;
  CT_star = ;
  time_ = ;
end

model kesmir_1999__A(A, CA, time_)

  // Assignment Rules:
  log_A := log(A / 1);

  // Rate Rules:
  A' = - z * A - u * CA * 1;

  // Variable initializations:
  A = 500;
  z = 0.02;
  u = 0.15;
  CA = ;
  time_ = ;
end

model kesmir_1999__T(T, alpha_T, CT_star, time_)

  // Rate Rules:
  T' = (sigma * 1 + p * alpha_T * CT_star) - delta_T * T;

  // Variable initializations:
  T = 0;
  p = 2;
  sigma = 5;
  delta_T = 0.8;
  alpha_T = ;
  CT_star = ;
  time_ = ;
end

model kesmir_1999__CA(CA, A, C)

  // Assignment Rules:
  CA := (C * A) / (SA * 1 + A);

  // Variable initializations:
  SA = 500;
  A = ;
  C = ;
end

model kesmir_1999__CT_star(CT_star, T, C_star, time_)

  // Assignment Rules:
  CT_star := (C_star * T) / (ST * 1 + C_star);

  // Variable initializations:
  ST = 50;
  T = ;
  C_star = ;
  time_ = ;
end

model kesmir_1999__alpha_B(alpha_B, B_sum, B0)

  // Assignment Rules:
  alpha_B := KB / (KB + B_sum / 1);

  // Variable initializations:
  KB = 1e4;
  B_sum = ;
  B0 = ;
end

model kesmir_1999__alpha_T(alpha_T, T)

  // Assignment Rules:
  alpha_T := KT / (KT + T / 1);

  // Variable initializations:
  KT = 100;
  T = ;
end

model kesmir_1999__total(B_sum, C_starsum)

  // Assignment Rules:
  total := B_sum + C_starsum;
  log_total := log(total / 1 + 1E-12);

  // Variable initializations:
  B_sum = ;
  C_starsum = ;
end

model kesmir_1999__kinetic_parameters(rho, mu, delta_B, pr)

  // Variable initializations:
  rho = 4;
  mu = 3;
  delta_B = 0.8;
  pr = 0.15;
end

model *kesmir_1999____main()

  // Sub-modules, and any changes to those submodules:
  environment: kesmir_1999__environment(time_);
  B0: kesmir_1999__B0(B00, pr, mu, rho, delta_B, CT_star0, time_);
  B1: kesmir_1999__B1(B11, alpha_B0, rho, delta_B, B00, time_);
  B2: kesmir_1999__B2(B20, alpha_B0, rho, delta_B, B11, time_);
  B3: kesmir_1999__B3(B30, alpha_B0, rho, delta_B, B20, time_);
  B4: kesmir_1999__B4(B40, alpha_B0, rho, delta_B, B30, time_);
  B5: kesmir_1999__B5(B50, alpha_B0, rho, delta_B, B40, time_);
  B6: kesmir_1999__B6(B60, alpha_B0, rho, delta_B, B50, time_);
  B7: kesmir_1999__B7(B70, alpha_B0, rho, delta_B, B60, time_);
  B8: kesmir_1999__B8(B80, alpha_B0, rho, delta_B, B70, time_);
  B9: kesmir_1999__B9(B90, alpha_B0, rho, delta_B, B80, time_);
  B10: kesmir_1999__B10(B100, alpha_B0, rho, delta_B, B90, time_);
  centroblasts_sum: kesmir_1999__centroblasts_sum(B_sum, B100, B90, B80, B70, B60, B50, B40, B30, B20, B11, B00);
  C: kesmir_1999__C(C0, rho, mu, B100, time_);
  C_star: kesmir_1999__C_star(C_star0, mu, CA0, time_);
  centrocytes_sum: kesmir_1999__centrocytes_sum(C_starsum, C0, C_star0, CT_star0, CA0);
  M: kesmir_1999__M(M0, pr, mu, CT_star0, time_);
  A: kesmir_1999__A(A0, CA0, time_);
  T: kesmir_1999__T(T0, alpha_T0, CT_star0, time_);
  CA: kesmir_1999__CA(CA0, A0, C0);
  CT_star: kesmir_1999__CT_star(CT_star0, T0, C_star0, time_);
  alpha_B: kesmir_1999__alpha_B(alpha_B0, B_sum, B00);
  alpha_T: kesmir_1999__alpha_T(alpha_T0, T0);
  total: kesmir_1999__total(B_sum, C_starsum);
  kinetic_parameters: kesmir_1999__kinetic_parameters(rho, mu, delta_B, pr);
end