// Created by libAntimony v2.13.0
model *Overgaard2007_PDmodel_IL21()

  // Compartments and Species:
  compartment COMpartment;

  // Assignment Rules:
  M_c := piecewise(M_night, ((tprime/3600) >= t_night) && ((tprime/3600) < t_day), M_day);
  k := kb + kinc*(T - T_b*(1 + pEtot*BR)) + f2_drug;
  f_prime := delta_high_dose*(1 + exp(-alpha*(time - (tdose1 + t_prime))))^-1;
  E_slow := AMT_dose*pEs2*(gNsTs1 + gNsTs2 + gNsTs3);
  E_fast := pEf2*(gNfTf1 + gNfTf2 + gNfTf3);
  gNsTs1 := piecewise((Ks^Ns/6)*exp(-Ks*X1)*X1^(Ns - 1), X1 > 0, 0);
  gNsTs2 := piecewise((Ks^Ns/6)*exp(-Ks*X2)*X2^(Ns - 1), X2 > 0, 0);
  gNsTs3 := piecewise((Ks^Ns/6)*exp(-Ks*X3)*X3^(Ns - 1), X3 > 0, 0);
  gNfTf1 := piecewise((Kf^Nf/6)*exp(-Kf*X1)*X1^(Nf - 1), X1 > 0, 0);
  gNfTf2 := piecewise((Kf^Nf/6)*exp(-Kf*X2)*X2^(Nf - 1), X2 > 0, 0);
  gNfTf3 := piecewise((Kf^Nf/6)*exp(-Kf*X3)*X3^(Nf - 1), X3 > 0, 0);
  M_night := (1 - f_prime)*M_night_baseline + f_prime*M_day;
  tprime := time*3600*1 - floor(time*3600*1/day_length)*day_length;
  M_day := (kb + kinc*(T_day - T_b))*(T_day - T_a);
  kb := M_b/(T_b - T_a);
  f2_drug := 0;
  T_day := T_b + delta_T/2;
  T_night := T_b - delta_T/2;
  M_night_baseline := (kb + kinc*(T_night - T_b))*(T_night - T_a);
  Ks := Ns/Ts;
  X1 := (time - tdose1)/24;
  X2 := (time - tdose2)/24;
  X3 := (time - tdose3)/24;
  Kf := Nf/Tf;

  // Rate Rules:
  M' = -km*(M - M_c);
  T' = c^-1*(M - k*(T - T_a));
  BR' = f_prime*(E_slow + E_fast)*(1 - BR) - kR*BR;

  // Compartment initializations:
  COMpartment = 1;

  // Variable initializations:
  M = 3.5;
  km = 1.1375;
  T = 38.785;
  c = 3.47;
  T_a = 21;
  BR = 0;
  kR = 5.35;
  AMT_dose = 3;
  pEs2 = 2.43;
  pEf2 = 3.57;
  delta_high_dose = 1;
  alpha = 0.2229166;
  tdose1 = 24;
  t_prime = 45.12;
  T_b = 38;
  delta_T = 1.57;
  kinc = 0.0258;
  tdose2 = 72;
  tdose3 = 120;
  t_night = 6.73;
  t_day = 17.5;
  day_length = 86400;
  pEtot = 0.144;
  pEf1 = 1;
  pEs1 = 0.2;
  pEf3 = 8;
  pEs3 = 50;
  M_b = 3;
  Ns = 4;
  Nf = 4;
  Ts = 2.45;
  Tf = 0.368;

  // Other declarations:
  var M, M_c, T, k, BR, f_prime, E_slow, E_fast, gNsTs1, gNsTs2, gNsTs3, gNfTf1;
  var gNfTf2, gNfTf3, M_night, tprime, M_day, kb, f2_drug, T_day, T_night;
  var M_night_baseline, Ks, X1, X2, X3, Kf;
  const COMpartment, km, c, T_a, kR, AMT_dose, pEs2, pEf2, delta_high_dose;
  const alpha, tdose1, t_prime, T_b, delta_T, kinc, tdose2, tdose3, t_night;
  const t_day, day_length, pEtot, pEf1, pEs1, pEf3, pEs3, M_b, Ns, Nf, Ts;
  const Tf;

  // Display Names:
  M is "Metabolic rate";
  km is "rate constant Metabolism";
  M_c is "circadian rhythm";
  T is "Temperature";
  c is "specific heat constant";
  k is "heat conductance";
  T_a is "ambient temperature";
  BR is "Bound Receptor";
  f_prime is "Priming";
  E_slow is "Slow Effect";
  E_fast is "Fast Effect";
  T_b is "basiline temperature";
  delta_T is "temperature difference";
  kb is "heat conductance baselinevalue";
  Ns is "No. of transit compartment (slow)";
  Nf is "No. of transit compartment (fast)";
  Ts is "mean total delay (slow)";
  Tf is "mena total delay (fast)";
end

Overgaard2007_PDmodel_IL21 is "Overgaard2007_PDmodel_IL21"

Overgaard2007_PDmodel_IL21 model_entity_is "http://identifiers.org/biomodels.db/MODEL0911110000"
Overgaard2007_PDmodel_IL21 model_entity_is "http://identifiers.org/biomodels.db/BIOMD0000000238"
Overgaard2007_PDmodel_IL21 description "http://identifiers.org/pubmed/17009101"
Overgaard2007_PDmodel_IL21 property "http://identifiers.org/mamo/MAMO_0000046"
Overgaard2007_PDmodel_IL21 taxon "http://identifiers.org/taxonomy/9541"
Overgaard2007_PDmodel_IL21 hypernym "http://identifiers.org/go/GO:0060086"
Overgaard2007_PDmodel_IL21 hypernym "http://identifiers.org/omim/605384"
