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

  // Variable initializations:
  time_ = ;
end

model raposo_2002__Q_Ca_p(Q_Ca_p, J_bp_Ca, J_pb_Ca, J_pu_Ca, J_ip_Ca, time_)

  // Rate Rules:
  Q_Ca_p' = (J_bp_Ca + J_ip_Ca) - (J_pb_Ca + J_pu_Ca);

  // Variable initializations:
  Q_Ca_p = 1.000;
  J_bp_Ca = ;
  J_pb_Ca = ;
  J_pu_Ca = ;
  J_ip_Ca = ;
  time_ = ;
end

model raposo_2002__Q_P_p(Q_P_p, J_bp_P, J_pb_P, J_pu_P, J_ip_P, J_pc_P, J_cp_P, time_)

  // Rate Rules:
  Q_P_p' = (J_bp_P + J_ip_P + J_cp_P) - (J_pb_P + J_pu_P + J_pc_P);

  // Variable initializations:
  Q_P_p = 1.000;
  J_bp_P = ;
  J_pb_P = ;
  J_pu_P = ;
  J_ip_P = ;
  J_pc_P = ;
  J_cp_P = ;
  time_ = ;
end

model raposo_2002__Q_P_c(Q_P_c, J_pc_P, J_cp_P, time_)

  // Rate Rules:
  Q_P_c' = J_pc_P - J_cp_P;

  // Variable initializations:
  Q_P_c = 3226.0;
  J_pc_P = ;
  J_cp_P = ;
  time_ = ;
end

model raposo_2002__Q_PTH_p(Q_PTH_p, S_PTH, C_PTH_p, k_PTH, time_)

  // Rate Rules:
  Q_PTH_p' = S_PTH - k_PTH * 1.0 * C_PTH_p;

  // Variable initializations:
  Q_PTH_p = 1.000;
  S_PTH = ;
  C_PTH_p = ;
  k_PTH = ;
  time_ = ;
end

model raposo_2002__Q_Ca_b(Q_Ca_b, J_bp_Ca, J_pb_Ca, time_)

  // Rate Rules:
  Q_Ca_b' = J_pb_Ca - J_bp_Ca;

  // Variable initializations:
  Q_Ca_b = 100.0;
  J_bp_Ca = ;
  J_pb_Ca = ;
  time_ = ;
end

model raposo_2002__Q_P_b(Q_P_b, J_bp_P, J_pb_P, time_)

  // Rate Rules:
  Q_P_b' = J_pb_P - J_bp_P;

  // Variable initializations:
  Q_P_b = 1.000;
  J_bp_P = ;
  J_pb_P = ;
  time_ = ;
end

model raposo_2002__Q_E_k(Q_E_k, S_E, k_E, time_)

  // Rate Rules:
  Q_E_k' = S_E - k_E * Q_E_k;

  // Variable initializations:
  Q_E_k = 1.000;
  S_E = ;
  k_E = ;
  time_ = ;
end

model raposo_2002__Q_D_p(Q_D_p, k_D, Q_E_k, time_)

  // Rate Rules:
  Q_D_p' = 1.0 * Q_E_k - k_D * Q_D_p;

  // Variable initializations:
  Q_D_p = 1.000;
  k_D = ;
  Q_E_k = ;
  time_ = ;
end

model raposo_2002__Q_TCa_i(Q_TCa_i, k1_D, k2_D, time_)

  // Rate Rules:
  Q_TCa_i' = (1.0 - Q_TCa_i) * k1_D - Q_TCa_i * k2_D;

  // Variable initializations:
  Q_TCa_i = 1.000;
  k1_D = ;
  k2_D = ;
  time_ = ;
end

model raposo_2002__Q_C_PT(Q_C_PT, k3_D, k4_D, time_)

  // Rate Rules:
  Q_C_PT' = (1.0 - Q_C_PT) * k3_D - Q_C_PT * k4_D;

  // Variable initializations:
  Q_C_PT = 1.000;
  k3_D = ;
  k4_D = ;
  time_ = ;
end

model raposo_2002__Q_TP_k(Q_TP_k, k1_P, k2_P, time_)

  // Rate Rules:
  Q_TP_k' = (1.0 - Q_TP_k) * k1_P - Q_TP_k * k2_P;

  // Variable initializations:
  Q_TP_k = 1.000;
  k1_P = ;
  k2_P = ;
  time_ = ;
end

model raposo_2002__model_fluxes(J_bp_Ca, J_pb_Ca, J_bp_P, J_pb_P, J_pu_Ca, J_pu_P, J_ip_Ca, J_ip_P, J_pc_P, J_cp_P, S_PTH, Q_Ca_b, Q_Ca_p, Q_TCa_i, Q_C_PT, Q_P_c, C_Ca_i, C_Ca_p, C_P_p, k_Ca_b, k1_D, k2_D, k3_D, k4_D, k3_P, k4_P, Y_Ca_i_2plus, Y_Ca_p_1minus)

  // Assignment Rules:
  J_pb_Ca := k_Ca_b * Q_Ca_b;
  J_bp_P := Stoic_Ca_P * J_bp_Ca;
  J_pb_P := Stoic_Ca_P * J_pb_Ca;
  J_ip_Ca := (Y_Ca_i_2plus * (1.0 - Q_TCa_i) * 1.0 * k1_D + k_Ca_i * 1.0 * (C_Ca_i - C_Ca_p)) - Q_TCa_i * k2_D;
  J_ip_P := J_P_ing;
  J_pc_P := k4_P * 1.0 * C_P_p;
  J_cp_P := k3_P * Q_P_c;
  S_PTH := 1.0 * Y_Ca_p_1minus * (k3_D * (1.0 - Q_C_PT) - k4_D * Q_C_PT);

  // Variable initializations:
  J_bp_Ca = 3.3;
  J_pu_Ca = 0.21;
  J_pu_P = 1.27;
  J_P_ing = 1.000;
  Stoic_Ca_P = 0.464;
  k_Ca_i = 1.000;
  Q_Ca_b = ;
  Q_Ca_p = ;
  Q_TCa_i = ;
  Q_C_PT = ;
  Q_P_c = ;
  C_Ca_i = ;
  C_Ca_p = ;
  C_P_p = ;
  k_Ca_b = ;
  k1_D = ;
  k2_D = ;
  k3_D = ;
  k4_D = ;
  k3_P = ;
  k4_P = ;
  Y_Ca_i_2plus = ;
  Y_Ca_p_1minus = ;
end

model raposo_2002__model_parameters(k1_D, k2_D, k3_D, k4_D, k1_P, k2_P, k3_P, k4_P, k_PTH, k_E, k_D, k_Ca_b, C_Ca_p, C_Ca_i, C_P_p, C_PTH_p, C_D_p, Y_Ca_i_2plus, Y_Ca_p_1minus, S_E, Q_TP_k)

  // Assignment Rules:
  k1_D := Y_i_D_p_1plus;
  k2_D := Y_i_D_p_1minus;
  k3_D := Y_PT_D_p_1plus;
  k4_D := Y_PT_D_p_1minus;
  k1_P := Y_k_Jext_1plus;
  Y_Ca_i_2plus := Y_Ca_i_2plus_Max * (C_Ca_i / (C_Ca_i + X_R_i_Ca));
  Y_Ca_p_1minus := Y_Ca_p_1plus_minus_Max * (a * (1.0 - tanh(b_PTH_Ca * (C_Ca_p - X_R_PTH_Ca))) + d);
  Ca_thr := 1.95 + Y_PTH_p_2plus;
  Ca_T := Ca_thr / C_Ca_p;
  P_thr := Y_PTH_p_2minus * Q_TP_k;
  P_T := P_thr / C_P_p;
  Y_i_D_p_1plus := Y_i_D_p_1plus_minus_Max * (a * (1.0 + tanh(b_i_D * (C_D_p - X_R_i_D))) + d);
  Y_PT_D_p_1plus := Y_PT_D_p_1plus_minus_Max * (a * (1.0 + tanh(b_PTH_D * (C_D_p - X_R_PTH_D))) + d);
  Y_k_Jext_1plus := Y_k_Jext_1plus_minus_Max * (a * (1.0 + tanh(b_E_PTH * (1.0 * C_P_p - X_R_E_PTH))) + d);
  Y_PTH_p_2plus := Y_PTH_p_2plus_Max * (C_PTH_p / (C_PTH_p + 1.0 * X_R_P_k_PTH));
  Y_PTH_p_2minus := Y_PTH_p_2minus_Max * (X_R_P_k_PTH / (1.0 * C_PTH_p + X_R_P_k_PTH));
  Y_i_D_p_1minus := Y_i_D_p_1plus_minus_Max * (a * (1.0 - tanh(b_i_D * (C_D_p - X_R_i_D))) + d);
  Y_PT_D_p_1minus := Y_PT_D_p_1plus_minus_Max * (a * (1.0 - tanh(b_PTH_D * (C_D_p - X_R_PTH_D))) + d);

  // Variable initializations:
  k2_P = 1.000;
  k3_P = 1.000;
  k4_P = 51.8;
  k_PTH = 100.0;
  k_E = 0.05;
  k_D = 0.1;
  k_Ca_b = 1.000;
  C_Ca_p = 2.4;
  C_Ca_i = 1.000;
  C_P_p = 1.2;
  C_PTH_p = 3.85;
  C_D_p = 90.0;
  S_E = 1.000;
  Y_Ca_i_2plus_Max = 1.000;
  Y_PTH_p_2plus_Max = 1.000;
  Y_PTH_p_2minus_Max = 1.000;
  Y_Ca_p_1plus_minus_Max = 0.02;
  Y_i_D_p_1plus_minus_Max = 0.02;
  Y_PT_D_p_1plus_minus_Max = 0.01;
  Y_k_Jext_1plus_minus_Max = 0.01;
  X_R_PTH_Ca = 1.0;
  X_R_i_Ca = 1.0;
  X_R_PTH_D = 90.0;
  X_R_E_PTH = 3.85;
  X_R_i_D = 90.0;
  X_R_P_k_PTH = 90.0;
  b_PTH_Ca = 0.05;
  b_PTH_D = 0.03;
  b_E_PTH = 0.55;
  b_i_D = 0.03;
  a = 0.85;
  d = 0.15;
  Q_TP_k = ;
end

model *raposo_2002____main()

  // Sub-modules, and any changes to those submodules:
  environment: raposo_2002__environment(time_);
  Q_Ca_p: raposo_2002__Q_Ca_p(Q_Ca_p0, J_bp_Ca, J_pb_Ca, J_pu_Ca, J_ip_Ca, time_);
  Q_P_p: raposo_2002__Q_P_p(Q_P_p0, J_bp_P, J_pb_P, J_pu_P, J_ip_P, J_pc_P, J_cp_P, time_);
  Q_P_c: raposo_2002__Q_P_c(Q_P_c0, J_pc_P, J_cp_P, time_);
  Q_PTH_p: raposo_2002__Q_PTH_p(Q_PTH_p0, S_PTH, C_PTH_p, k_PTH, time_);
  Q_Ca_b: raposo_2002__Q_Ca_b(Q_Ca_b0, J_bp_Ca, J_pb_Ca, time_);
  Q_P_b: raposo_2002__Q_P_b(Q_P_b0, J_bp_P, J_pb_P, time_);
  Q_E_k: raposo_2002__Q_E_k(Q_E_k0, S_E, k_E, time_);
  Q_D_p: raposo_2002__Q_D_p(Q_D_p0, k_D, Q_E_k0, time_);
  Q_TCa_i: raposo_2002__Q_TCa_i(Q_TCa_i0, k1_D, k2_D, time_);
  Q_C_PT: raposo_2002__Q_C_PT(Q_C_PT0, k3_D, k4_D, time_);
  Q_TP_k: raposo_2002__Q_TP_k(Q_TP_k0, k1_P, k2_P, time_);
  model_fluxes: raposo_2002__model_fluxes(J_bp_Ca, J_pb_Ca, J_bp_P, J_pb_P, J_pu_Ca, J_pu_P, J_ip_Ca, J_ip_P, J_pc_P, J_cp_P, S_PTH, Q_Ca_b0, Q_Ca_p0, Q_TCa_i0, Q_C_PT0, Q_P_c0, C_Ca_i, C_Ca_p, C_P_p, k_Ca_b, k1_D, k2_D, k3_D, k4_D, k3_P, k4_P, Y_Ca_i_2plus, Y_Ca_p_1minus);
  model_parameters: raposo_2002__model_parameters(k1_D, k2_D, k3_D, k4_D, k1_P, k2_P, k3_P, k4_P, k_PTH, k_E, k_D, k_Ca_b, C_Ca_p, C_Ca_i, C_P_p, C_PTH_p, C_D_p, Y_Ca_i_2plus, Y_Ca_p_1minus, S_E, Q_TP_k0);
end