// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;

// 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 L;
  real Endoplasmic_Reticulum = 1 L;
  real Mitochondria = 1 L;
  real Ca_cyt(time) uM;
  real CaER(time) uM;
  real CaM(time) uM;
  real CaPr(time) uM;
  real Pr(time) uM;
  real v1(time) microkatal;
  real Kch = 4100;
  real K1 = 5;
  real v3(time) microkatal;
  real Kleak = .05;
  real v5(time) microkatal;
  real Kpump = 20;
  real v7(time) microkatal;
  real Kout = 125;
  real K3 = 5;
  real Km = .00625;
  real v9(time) microkatal;
  real Kin = 300;
  real K2 = .8;
  real v11(time) microkatal;
  real Kminus = .01;
  real v12(time) microkatal;
  real Kplus = .1;

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

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