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

  // Variable initializations:
  time_ = ;
end

model starbuck_1990__Rs(Rs, kfc, kr, ks, kt, krec_c, krec_i, numcell, fr, Ric, Rii, Lo, Cs, time_)

  // Rate Rules:
  Rs' = (kr * Cs + krec_c * Ric + krec_i * (1.0 - fr) * Rii + ks) - ((kfc / numcell) * Rs * Lo + kt * Rs);

  // Variable initializations:
  Rs = ;
  kfc = ;
  kr = ;
  ks = ;
  kt = ;
  krec_c = ;
  krec_i = ;
  numcell = ;
  fr = ;
  Ric = ;
  Rii = ;
  Lo = ;
  Cs = ;
  time_ = ;
end

model starbuck_1990__Cs(Cs, kfc, kr, kc, ku, kt, krec_c, krec_i, numcell, fr, Ts, Ps, Lo, Rs, Cic, Cii, time_)

  // Rate Rules:
  Cs' = (ku * Ts + krec_c * Cic + krec_i * (1.0 - fr) * Cii) - ((kfc / numcell) * Rs * Lo + kr * Cs + kc * Cs * Ps + kt * Cs);

  // Variable initializations:
  Cs = ;
  kfc = ;
  kr = ;
  kc = ;
  ku = ;
  kt = ;
  krec_c = ;
  krec_i = ;
  numcell = ;
  fr = ;
  Ts = ;
  Ps = ;
  Lo = ;
  Rs = ;
  Cic = ;
  Cii = ;
  time_ = ;
end

model starbuck_1990__Ts(Ts, kc, ku, Cs, Ps, lamda, time_)

  // Rate Rules:
  Ts' = kc * Cs * Ps - (ku * Ts + lamda * Ts);

  // Variable initializations:
  Ts = ;
  kc = ;
  ku = ;
  Cs = ;
  Ps = ;
  lamda = ;
  time_ = ;
end

model starbuck_1990__Ps(Ps, kc, ku, krec_p, Ts, Cs, Pi_, lamda, time_)

  // Rate Rules:
  Ps' = (ku * Ts + krec_p * Pi_) - (kc * Cs * Ps + lamda * Ps);

  // Variable initializations:
  Ps = ;
  kc = ;
  ku = ;
  krec_p = ;
  Ts = ;
  Cs = ;
  Pi_ = ;
  lamda = ;
  time_ = ;
end

model starbuck_1990__Lo(Lo, kr, kfc, krec_c, Rs, Cs, Lic, numcell, time_)

  // Rate Rules:
  Lo' = (numcell * kr * Cs + krec_c * Lic) - kfc * Rs * Lo;

  // Variable initializations:
  Lo = ;
  kr = ;
  kfc = ;
  krec_c = ;
  Rs = ;
  Cs = ;
  Lic = ;
  numcell = ;
  time_ = ;
end

model starbuck_1990__Rii(Rii, kfe, kr_, kdeg, krec_i, numcell, fr, Lii, Cii, time_)

  // Rate Rules:
  Rii' = kr_ * Cii - ((kfe / numcell) * Rii * Lii + krec_i * (1.0 - fr) * Rii + fr * kdeg * Rii);

  // Variable initializations:
  Rii = ;
  kfe = ;
  kr_ = ;
  kdeg = ;
  krec_i = ;
  numcell = ;
  fr = ;
  Lii = ;
  Cii = ;
  time_ = ;
end

model starbuck_1990__Cii(Cii, kfe, kr_, ku_, kdeg, krec_i, numcell, fr, Ti, Lii, Rii, time_)

  // Rate Rules:
  Cii' = ((kfe / numcell) * Rii * Lii + ku_ * Ti) - (kr_ * Cii + fr * kdeg * Cii + krec_i * (1.0 - fr) * Cii);

  // Variable initializations:
  Cii = ;
  kfe = ;
  kr_ = ;
  ku_ = ;
  kdeg = ;
  krec_i = ;
  numcell = ;
  fr = ;
  Ti = ;
  Lii = ;
  Rii = ;
  time_ = ;
end

model starbuck_1990__Ti(Ti, kdeg, ku_, Ts, lamda, time_)

  // Rate Rules:
  Ti' = lamda * Ts - (ku_ * Ti + kdeg * Ti);

  // Variable initializations:
  Ti = ;
  kdeg = ;
  ku_ = ;
  Ts = ;
  lamda = ;
  time_ = ;
end

model starbuck_1990__Pi(Pi_, ku_, krec_p, Ti, Ps, lamda, time_)

  // Rate Rules:
  Pi_' = (ku_ * Ti + lamda * Pi_) - krec_p * Pi_;

  // Variable initializations:
  Pi_ = ;
  ku_ = ;
  krec_p = ;
  Ti = ;
  Ps = ;
  lamda = ;
  time_ = ;
end

model starbuck_1990__Lii(Lii, kr_, kfe, kdeg, Rii, Cii, numcell, time_)

  // Rate Rules:
  Lii' = numcell * kr_ * Cii - (kdeg * Lii + kfe * Rii * Lii);

  // Variable initializations:
  Lii = ;
  kr_ = ;
  kfe = ;
  kdeg = ;
  Rii = ;
  Cii = ;
  numcell = ;
  time_ = ;
end

model starbuck_1990__Ric(Ric, kfv, kr, kt, krec_c, numcell, Rs, Lic, Cic, time_)

  // Rate Rules:
  Ric' = (kr * Cic + kt * Rs) - ((kfv / numcell) * Ric * Lic + krec_c * Ric);

  // Variable initializations:
  Ric = ;
  kfv = ;
  kr = ;
  kt = ;
  krec_c = ;
  numcell = ;
  Rs = ;
  Lic = ;
  Cic = ;
  time_ = ;
end

model starbuck_1990__Cic(Cic, kfv, kr, kt, krec_c, numcell, Cs, Lic, Ric, time_)

  // Rate Rules:
  Cic' = ((kfv / numcell) * Ric * Lic + kt * Cs) - (kr * Cic + krec_c * Cic);

  // Variable initializations:
  Cic = ;
  kfv = ;
  kr = ;
  kt = ;
  krec_c = ;
  numcell = ;
  Cs = ;
  Lic = ;
  Ric = ;
  time_ = ;
end

model starbuck_1990__Lic(Lic, kr, kfv, krec_c, Ric, Cic, numcell, time_)

  // Rate Rules:
  Lic' = numcell * kr * Cic - (krec_c * Lic + kfv * Ric * Lic);

  // Variable initializations:
  Lic = ;
  kr = ;
  kfv = ;
  krec_c = ;
  Ric = ;
  Cic = ;
  numcell = ;
  time_ = ;
end

model starbuck_1990__kinetic_constants(kr, kc, kuc, ks, lamda, kt, kr_, ku_, krec_c, krec_i, krec_p, kdeg, ku, kfc, kfe, kfv, numcell, fr)

  // Assignment Rules:
  ku := kc / Kcp;
  kfc := (kf * numcell) / (Nav * Vplate);
  kfe := kf / (Nav * Vendosome);
  kfv := kf / (Nav * Vvesicle);
  numcell := cell / plate;
  Vplate := L / plate;

  // Variable initializations:
  kr = 3.6E-2;
  kc = 3.8E-7;
  kuc = 1.2E-2;
  ks = 1.0;
  lamda = 1.8;
  kt = 0.05;
  kr_ = 3.6E-1;
  ku_ = 0.1216;
  krec_c = 0.05;
  krec_i = 0.05;
  krec_p = 20.0;
  kdeg = 0.02;
  fr = ;
  kf = 1.2E6;
  K_DC = 3.0E-8;
  K_DT = 3.2E5;
  Kcp = 3.125E-6;
  RT = 1.8E5;
  PT = 6.6E4;
  PT_RT = 0.37;
  cell = ;
  plate = ;
  L = ;
  Vendosome = 5.0E-19;
  Vvesicle = 5.0E-15;
  Nav = 6.0221419E23;
end

model *starbuck_1990____main()

  // Sub-modules, and any changes to those submodules:
  environment: starbuck_1990__environment(time_);
  Rs: starbuck_1990__Rs(Rs0, kfc, kr, ks, kt, krec_c, krec_i, numcell, fr, Ric0, Rii0, Lo0, Cs0, time_);
  Cs: starbuck_1990__Cs(Cs0, kfc, kr, kc, ku, kt, krec_c, krec_i, numcell, fr, Ts0, Ps0, Lo0, Rs0, Cic0, Cii0, time_);
  Ts: starbuck_1990__Ts(Ts0, kc, ku, Cs0, Ps0, lamda, time_);
  Ps: starbuck_1990__Ps(Ps0, kc, ku, krec_p, Ts0, Cs0, Pi_0, lamda, time_);
  Lo: starbuck_1990__Lo(Lo0, kr, kfc, krec_c, Rs0, Cs0, Lic0, numcell, time_);
  Rii: starbuck_1990__Rii(Rii0, kfe, kr_, kdeg, krec_i, numcell, fr, Lii0, Cii0, time_);
  Cii: starbuck_1990__Cii(Cii0, kfe, kr_, ku_, kdeg, krec_i, numcell, fr, Ti0, Lii0, Rii0, time_);
  Ti: starbuck_1990__Ti(Ti0, kdeg, ku_, Ts0, lamda, time_);
  Pi_: starbuck_1990__Pi(Pi_0, ku_, krec_p, Ti0, Ps0, lamda, time_);
  Lii: starbuck_1990__Lii(Lii0, kr_, kfe, kdeg, Rii0, Cii0, numcell, time_);
  Ric: starbuck_1990__Ric(Ric0, kfv, kr, kt, krec_c, numcell, Rs0, Lic0, Cic0, time_);
  Cic: starbuck_1990__Cic(Cic0, kfv, kr, kt, krec_c, numcell, Cs0, Lic0, Ric0, time_);
  Lic: starbuck_1990__Lic(Lic0, kr, kfv, krec_c, Ric0, Cic0, numcell, time_);
  kinetic_constants: starbuck_1990__kinetic_constants(kr, kc, kuc, ks, lamda, kt, kr_, ku_, krec_c, krec_i, krec_p, kdeg, ku, kfc, kfe, kfv, numcell, fr);
end