// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = 1E-6 mole;
unit volume = litre;
unit area = metre^2;
unit length = metre;
unit time = second;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// v1: .25CaER <=> Ca_cyt 
// v3: .25CaER <=> Ca_cyt 
// v5: Ca_cyt <=> .25CaER 
// v7: .25CaM <=> Ca_cyt 
// v9: Ca_cyt <=> .25CaM 
// v11: CaPr <=> Pr Ca_cyt 
// v12: Pr Ca_cyt <=> CaPr 

math main {
  realDomain time second;
  time.min=0;
  extern time.max;
  extern time.delta;

  // variable definitions
  real Cytosol = 1 volume;
  real Endoplasmic_Reticulum = 1 volume;
  real Mitochondria = 1 volume;
  private real Ca_cyt.amt(time) substance;
  real Ca_cyt(time) substance/volume;
  real Ca_cyt.init substance/volume;
  private real CaER.amt(time) substance;
  real CaER(time) substance/volume;
  real CaER.init substance/volume;
  private real CaM.amt(time) substance;
  real CaM(time) substance/volume;
  real CaM.init substance/volume;
  private real CaPr.amt(time) substance;
  real CaPr(time) substance/volume;
  real CaPr.init substance/volume;
  private real Pr.amt(time) substance;
  real Pr(time) substance/volume;
  real Pr.init substance/volume;
  real v1.rate(time) substance/time;
  real Kch = 4100;
  real K1 = 5;
  real v3.rate(time) substance/time;
  real Kleak = .05;
  real v5.rate(time) substance/time;
  real Kpump = 20;
  real v7.rate(time) substance/time;
  real Kout = 125;
  real K3 = 5;
  real Km = .00625;
  real v9.rate(time) substance/time;
  real Kin = 300;
  real K2 = .8;
  real v11.rate(time) substance/time;
  real Kminus = .01;
  real v12.rate(time) substance/time;
  real Kplus = .1;

  // equations
  when (time=time.min) Ca_cyt.amt = Ca_cyt.init*Cytosol;
  Ca_cyt.amt:time = v1.rate + v3.rate + -1*v5.rate + v7.rate + -1*v9.rate + v11.rate + -1*v12.rate;
  Ca_cyt = Ca_cyt.amt/Cytosol;
  Ca_cyt.init = .35;
  when (time=time.min) CaER.amt = CaER.init*Endoplasmic_Reticulum;
  CaER.amt:time = -.25*v1.rate + -.25*v3.rate + .25*v5.rate;
  CaER = CaER.amt/Endoplasmic_Reticulum;
  CaER.init = .76;
  when (time=time.min) CaM.amt = CaM.init*Mitochondria;
  CaM.amt:time = -.25*v7.rate + .25*v9.rate;
  CaM = CaM.amt/Mitochondria;
  CaM.init = .29;
  when (time=time.min) CaPr.amt = CaPr.init*Cytosol;
  CaPr.amt:time = -1*v11.rate + v12.rate;
  CaPr = CaPr.amt/Cytosol;
  CaPr.init = 85.45;
  when (time=time.min) Pr.amt = Pr.init*Cytosol;
  Pr.amt:time = v11.rate + -1*v12.rate;
  Pr = Pr.amt/Cytosol;
  Pr.init = 34.55;
  v1.rate = Cytosol*(Kch*Ca_cyt^2*(CaER-Ca_cyt)/(K1^2+Ca_cyt^2));
  v3.rate = Cytosol*Kleak*(CaER-Ca_cyt);
  v5.rate = Endoplasmic_Reticulum*Kpump*Ca_cyt;
  v7.rate = Cytosol*CaM*(Kout*Ca_cyt^2/(K3^2+Ca_cyt^2)+Km);
  v9.rate = Mitochondria*(Kin*Ca_cyt^8/(K2^8+Ca_cyt^8));
  v11.rate = Cytosol*Kminus*CaPr;
  v12.rate = Cytosol*Kplus*Ca_cyt*Pr;

  // variable properties
  Cytosol.sbmlRole="compartment";
  Endoplasmic_Reticulum.sbmlRole="compartment";
  Mitochondria.sbmlRole="compartment";
  Ca_cyt.amt.sbmlRole="speciesAmount";
  Ca_cyt.sbmlRole="speciesConcentration";
  Ca_cyt.sbmlCompartment="Cytosol";
  Ca_cyt.init.sbmlRole="speciesInitialConcentration";
  CaER.amt.sbmlRole="speciesAmount";
  CaER.sbmlRole="speciesConcentration";
  CaER.sbmlCompartment="Endoplasmic_Reticulum";
  CaER.init.sbmlRole="speciesInitialConcentration";
  CaM.amt.sbmlRole="speciesAmount";
  CaM.sbmlRole="speciesConcentration";
  CaM.sbmlCompartment="Mitochondria";
  CaM.init.sbmlRole="speciesInitialConcentration";
  CaPr.amt.sbmlRole="speciesAmount";
  CaPr.sbmlRole="speciesConcentration";
  CaPr.sbmlCompartment="Cytosol";
  CaPr.init.sbmlRole="speciesInitialConcentration";
  Pr.amt.sbmlRole="speciesAmount";
  Pr.sbmlRole="speciesConcentration";
  Pr.sbmlCompartment="Cytosol";
  Pr.init.sbmlRole="speciesInitialConcentration";
  v1.rate.sbmlRole="rate";
  Kch.sbmlRole="parameter";
  K1.sbmlRole="parameter";
  v3.rate.sbmlRole="rate";
  Kleak.sbmlRole="parameter";
  v5.rate.sbmlRole="rate";
  Kpump.sbmlRole="parameter";
  v7.rate.sbmlRole="rate";
  Kout.sbmlRole="parameter";
  K3.sbmlRole="parameter";
  Km.sbmlRole="parameter";
  v9.rate.sbmlRole="rate";
  Kin.sbmlRole="parameter";
  K2.sbmlRole="parameter";
  v11.rate.sbmlRole="rate";
  Kminus.sbmlRole="parameter";
  v12.rate.sbmlRole="rate";
  Kplus.sbmlRole="parameter";
}