//Created by libAntimony v2.4
model cooling_2007__all(Pc, Pcg, J10, J11)

  // Assignment Rules:
  L := piecewise( Ls / (1 + exp(-(80) * ((t - ts) - 0.05))) , ( t < ts + 0.15) && (geq(t, ts) ), Ls , geq(t, ts + 0.15) , 0 );
  J1 := kf1 * R * L - kr1 * Rl;
  kr1 := kf1 * Kd1;
  J2 := kf2 * R * Gd - kr2 * Rg;
  kr2 := kf2 * Kd2;
  J3 := kf3 * Rl * Gd - kr3 * Rlg;
  J4 := kf4 * L * Rg - kr4 * Rlg;
  kr4 := kf4 * Kd4;
  J5 := kf5 * Rlg;
  J6 := kf6 * Rlg;
  J7 := kf7 * Gt_;
  J8 := kf8 * P * Ca - kr8 * Pc;
  J9 := kf9 * P * Gt_ - kr9 * Pg;
  J10 := kf10 * Pc * Gt_ - kr10 * Pcg;
  J11 := kf11 * Pg * Ca - kr11 * Pcg;
  kr11 := kf11 * Kd11;
  J12 := kf12 * Pcg;
  J13 := kf13 * Pg;
  J14 := (kf14 * Pc * PIP2) / (Km14 / Cpc + PIP2);
  J15 := (kf15 * Pcg * PIP2) / (Km15 / Cpc + PIP2);
  J16 := kf16 * IP3;
  Cpc := Cc / Cp;
  Cc := 1 / (Vc * 6.022E2);
  Cp := 1 / (Vc * Rpc);

  // Rate Rules:
  Ca' = Cpc * -(1) * (J8 + J11);
  Gd' = (J7 + J13 + J12) - (J2 + J3);
  Gt_' = J6 - (J7 + J9 + J10);
  R' = -(1) * (J1 + J2);
  Rl' = (J1 + J6) - J3;
  Rg' = J2 - J4;
  Rlg' = ((J3 - J5) + J4) - J6;
  Rlgp' = J5;
  IP3' = Cpc * (J14 + J15) - J16;
  Pc' = (J8 + J12) - J10;
  Pcg' = (J10 + J11) - J12;
  P' = J13 - (J9 + J8);
  Pg' = J9 - (J11 + J13);

  // Variable initializations:
  t = ;
  Ca = 0.1;
  Ls = 0.1;
  ts = 10;
  Gd = 10000;
  Gt_ = 0;
  R = 13.9;
  Rl = 0;
  Rg = 5.06;
  Rlg = 0;
  Rlgp = 0;
  IP3 = 0.015;
  PIP2 = 4000;
  Pc = 9.09;
  Pcg = 0;
  P = 90.9;
  Pg = 0;
  kf1 = 0.0003;
  Kd1 = 3e-5;
  kf2 = 2.75e-4;
  Kd2 = 27500;
  kf3 = 1;
  kr3 = 0.001;
  kf4 = 0.3;
  Kd4 = 3e-5;
  kf5 = 0.0004;
  kf6 = 1;
  kf7 = 0.15;
  kf8 = 0.0167;
  kr8 = 0.0167;
  kf9 = 0.0042;
  kr9 = 1;
  kf10 = 0.042;
  kr10 = 1;
  kf11 = 0.0334;
  Kd11 = 0.1;
  kf12 = 6;
  kf13 = 6;
  kf14 = 0.444;
  Km14 = 19.8;
  kf15 = 3.8;
  Km15 = 5;
  kf16 = 1.25;
  Vc = 2550;
  Rpc = 4.61;
end

model *cooling_2007____main()

  // Sub-modules, and any changes to those submodules:
  all: cooling_2007__all(Pc, Pcg, J10, J11);
end