// 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
// Reaction1: A R <=> C 
// Reaction2: A <=> EmptySet 
// Reaction3: C <=> R 
// Reaction4: R <=> EmptySet 
// Reaction5: A DA <=> DAp 
// Reaction6: DAp <=> A DA 
// Reaction7: DA <=> DA MA 
// Reaction8: DAp <=> DAp MA 
// Reaction9: MA <=> EmptySet 
// Reaction10: MA <=> A MA 
// Reaction11: A DR <=> DRp 
// Reaction12: DRp <=> A DR 
// Reaction13: DR <=> DR MR 
// Reaction14: DRp <=> DRp MR 
// Reaction15: MR <=> EmptySet 
// Reaction16: MR <=> MR R 

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

  // variable definitions
  real deterministicOscillator = 1 L;
  real EmptySet = 0 mol;
  real A(time) mol;
  real C(time) mol;
  real DA(time) mol;
  real DAp(time) mol;
  real DR(time) mol;
  real DRp(time) mol;
  real MA(time) mol;
  real MR(time) mol;
  real R(time) mol;
  real Reaction1(time) mole/hr;
  real gammaC = 2;
  real Reaction2(time) mole/hr;
  real deltaA = 1;
  real Reaction3(time) mole/hr;
  real Reaction4(time) mole/hr;
  real deltaR = .2;
  real Reaction5(time) mole/hr;
  real gammaA = 1;
  real Reaction6(time) mole/hr;
  real thetaA = 50;
  real Reaction7(time) mole/hr;
  real alphaA = 50;
  real Reaction8(time) mole/hr;
  real alphaAp = 500;
  real Reaction9(time) mole/hr;
  real deltaMA = 10;
  real Reaction10(time) mole/hr;
  real betaA = 50;
  real Reaction11(time) mole/hr;
  real gammaR = 1;
  real Reaction12(time) mole/hr;
  real thetaR = 100;
  real Reaction13(time) mole/hr;
  real alphaR = .01;
  real Reaction14(time) mole/hr;
  real alphaRp = 50;
  real Reaction15(time) mole/hr;
  real deltaMR = .5;
  real Reaction16(time) mole/hr;
  real betaR = 5;

  // equations
  when (time=time.min) A = 0;
  A:time = -1*Reaction1 + -1*Reaction2 + -1*Reaction5 + Reaction6 + Reaction10 + -1*Reaction11 + Reaction12;
  when (time=time.min) C = 0;
  C:time = Reaction1 + -1*Reaction3;
  when (time=time.min) DA = 1;
  DA:time = -1*Reaction5 + Reaction6 + -1*Reaction7 + Reaction7;
  when (time=time.min) DAp = 0;
  DAp:time = Reaction5 + -1*Reaction6 + -1*Reaction8 + Reaction8;
  when (time=time.min) DR = 1;
  DR:time = -1*Reaction11 + Reaction12 + -1*Reaction13 + Reaction13;
  when (time=time.min) DRp = 0;
  DRp:time = Reaction11 + -1*Reaction12 + -1*Reaction14 + Reaction14;
  when (time=time.min) MA = 0;
  MA:time = Reaction7 + Reaction8 + -1*Reaction9 + -1*Reaction10 + Reaction10;
  when (time=time.min) MR = 0;
  MR:time = Reaction13 + Reaction14 + -1*Reaction15 + -1*Reaction16 + Reaction16;
  when (time=time.min) R = 0;
  R:time = -1*Reaction1 + Reaction3 + -1*Reaction4 + Reaction16;
  Reaction1 = A*R*gammaC;
  Reaction2 = A*deltaA;
  Reaction3 = C*deltaA;
  Reaction4 = R*deltaR;
  Reaction5 = A*DA*gammaA;
  Reaction6 = DAp*thetaA;
  Reaction7 = DA*alphaA;
  Reaction8 = DAp*alphaAp;
  Reaction9 = MA*deltaMA;
  Reaction10 = MA*betaA;
  Reaction11 = A*DR*gammaR;
  Reaction12 = DRp*thetaR;
  Reaction13 = DR*alphaR;
  Reaction14 = DRp*alphaRp;
  Reaction15 = MR*deltaMR;
  Reaction16 = MR*betaR;

  // variable properties
  deterministicOscillator.sbmlRole="compartment";
  EmptySet.sbmlRole="species";
  EmptySet.sbmlCompartment="deterministicOscillator";
  A.sbmlRole="species";
  A.sbmlCompartment="deterministicOscillator";
  C.sbmlRole="species";
  C.sbmlCompartment="deterministicOscillator";
  DA.sbmlRole="species";
  DA.sbmlCompartment="deterministicOscillator";
  DAp.sbmlRole="species";
  DAp.sbmlCompartment="deterministicOscillator";
  DR.sbmlRole="species";
  DR.sbmlCompartment="deterministicOscillator";
  DRp.sbmlRole="species";
  DRp.sbmlCompartment="deterministicOscillator";
  MA.sbmlRole="species";
  MA.sbmlCompartment="deterministicOscillator";
  MR.sbmlRole="species";
  MR.sbmlCompartment="deterministicOscillator";
  R.sbmlRole="species";
  R.sbmlCompartment="deterministicOscillator";
  Reaction1.sbmlRole="rate";
  gammaC.sbmlRole="parameter";
  Reaction2.sbmlRole="rate";
  deltaA.sbmlRole="parameter";
  Reaction3.sbmlRole="rate";
  Reaction4.sbmlRole="rate";
  deltaR.sbmlRole="parameter";
  Reaction5.sbmlRole="rate";
  gammaA.sbmlRole="parameter";
  Reaction6.sbmlRole="rate";
  thetaA.sbmlRole="parameter";
  Reaction7.sbmlRole="rate";
  alphaA.sbmlRole="parameter";
  Reaction8.sbmlRole="rate";
  alphaAp.sbmlRole="parameter";
  Reaction9.sbmlRole="rate";
  deltaMA.sbmlRole="parameter";
  Reaction10.sbmlRole="rate";
  betaA.sbmlRole="parameter";
  Reaction11.sbmlRole="rate";
  gammaR.sbmlRole="parameter";
  Reaction12.sbmlRole="rate";
  thetaR.sbmlRole="parameter";
  Reaction13.sbmlRole="rate";
  alphaR.sbmlRole="parameter";
  Reaction14.sbmlRole="rate";
  alphaRp.sbmlRole="parameter";
  Reaction15.sbmlRole="rate";
  deltaMR.sbmlRole="parameter";
  Reaction16.sbmlRole="rate";
  betaR.sbmlRole="parameter";
}