// This model generated automatically from SBML

// unit definitions
import nsrunit;
unit conversion off;
unit item=scalar;
unit substance = 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: s => ACA 
// v2: ACA => p 
// v3: s => PKA 
// v4: PKA => p 
// v5: s => ERK2 
// v6: ERK2 => p 
// v7: s => REGA 
// v8: REGA => p 
// v9: s => cAMP 
// v10: cAMP => p 
// v11: s => excAMP 
// v12: excAMP => p 
// v13: s => CAR1 
// v14: CAR1 => p 

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

  // variable definitions
  real compartment = 1 volume;
  real k1 = 2;
  real k2 = .9;
  real k3 = 2.5;
  real k4 = 1.5;
  real k5 = .6;
  real k6 = .8;
  real k7 = 1;
  real k8 = 1.3;
  real k9 = .3;
  real k10 = .8;
  real k11 = .7;
  real k12 = 4.9;
  real k13 = 23;
  real k14 = 4.5;
  private real s.amt substance;
  real s substance/volume;
  real s.init substance/volume;
  private real ACA.amt(time) substance;
  real ACA(time) substance/volume;
  real ACA.init substance/volume;
  private real CAR1.amt(time) substance;
  real CAR1(time) substance/volume;
  real CAR1.init substance/volume;
  private real p.amt substance;
  real p substance/volume;
  real p.init substance/volume;
  private real PKA.amt(time) substance;
  real PKA(time) substance/volume;
  real PKA.init substance/volume;
  private real cAMP.amt(time) substance;
  real cAMP(time) substance/volume;
  real cAMP.init substance/volume;
  private real ERK2.amt(time) substance;
  real ERK2(time) substance/volume;
  real ERK2.init substance/volume;
  private real REGA.amt(time) substance;
  real REGA(time) substance/volume;
  real REGA.init substance/volume;
  private real excAMP.amt(time) substance;
  real excAMP(time) substance/volume;
  real excAMP.init substance/volume;
  real v1.rate(time) substance/time;
  real v2.rate(time) substance/time;
  real v3.rate(time) substance/time;
  real v4.rate(time) substance/time;
  real v5.rate(time) substance/time;
  real v6.rate(time) substance/time;
  real v7.rate(time) substance/time;
  real v8.rate(time) substance/time;
  real v9.rate(time) substance/time;
  real v10.rate(time) substance/time;
  real v11.rate(time) substance/time;
  real v12.rate(time) substance/time;
  real v13.rate(time) substance/time;
  real v14.rate(time) substance/time;

  // equations
  s.amt = s*compartment;
  s = s.init;
  s.init = 1;
  when (time=time.min) ACA.amt = ACA.init*compartment;
  ACA.amt:time = v1.rate + -1*v2.rate;
  ACA = ACA.amt/compartment;
  ACA.init = 3.39;
  when (time=time.min) CAR1.amt = CAR1.init*compartment;
  CAR1.amt:time = v13.rate + -1*v14.rate;
  CAR1 = CAR1.amt/compartment;
  CAR1.init = 2.45;
  p.amt = p*compartment;
  p = p.init;
  p.init = 1;
  when (time=time.min) PKA.amt = PKA.init*compartment;
  PKA.amt:time = v3.rate + -1*v4.rate;
  PKA = PKA.amt/compartment;
  PKA.init = 2.2;
  when (time=time.min) cAMP.amt = cAMP.init*compartment;
  cAMP.amt:time = v9.rate + -1*v10.rate;
  cAMP = cAMP.amt/compartment;
  cAMP.init = 1.38;
  when (time=time.min) ERK2.amt = ERK2.init*compartment;
  ERK2.amt:time = v5.rate + -1*v6.rate;
  ERK2 = ERK2.amt/compartment;
  ERK2.init = 1.13;
  when (time=time.min) REGA.amt = REGA.init*compartment;
  REGA.amt:time = v7.rate + -1*v8.rate;
  REGA = REGA.amt/compartment;
  REGA.init = 1.24;
  when (time=time.min) excAMP.amt = excAMP.init*compartment;
  excAMP.amt:time = v11.rate + -1*v12.rate;
  excAMP = excAMP.amt/compartment;
  excAMP.init = .48;
  v1.rate = k1*CAR1;
  v2.rate = k2*ACA*PKA;
  v3.rate = k3*cAMP;
  v4.rate = k4*PKA;
  v5.rate = k5*CAR1;
  v6.rate = k6*PKA*ERK2;
  v7.rate = k7;
  v8.rate = k8*ERK2*REGA;
  v9.rate = k9*ACA;
  v10.rate = k10*REGA*cAMP;
  v11.rate = k11*ACA;
  v12.rate = k12*excAMP;
  v13.rate = k13*excAMP;
  v14.rate = k14*CAR1;

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  k3.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k5.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k7.sbmlRole="parameter";
  k8.sbmlRole="parameter";
  k9.sbmlRole="parameter";
  k10.sbmlRole="parameter";
  k11.sbmlRole="parameter";
  k12.sbmlRole="parameter";
  k13.sbmlRole="parameter";
  k14.sbmlRole="parameter";
  s.amt.sbmlRole="speciesAmount";
  s.sbmlRole="speciesConcentration";
  s.sbmlCompartment="compartment";
  s.init.sbmlRole="speciesInitialConcentration";
  ACA.amt.sbmlRole="speciesAmount";
  ACA.sbmlRole="speciesConcentration";
  ACA.sbmlCompartment="compartment";
  ACA.init.sbmlRole="speciesInitialConcentration";
  CAR1.amt.sbmlRole="speciesAmount";
  CAR1.sbmlRole="speciesConcentration";
  CAR1.sbmlCompartment="compartment";
  CAR1.init.sbmlRole="speciesInitialConcentration";
  p.amt.sbmlRole="speciesAmount";
  p.sbmlRole="speciesConcentration";
  p.sbmlCompartment="compartment";
  p.init.sbmlRole="speciesInitialConcentration";
  PKA.amt.sbmlRole="speciesAmount";
  PKA.sbmlRole="speciesConcentration";
  PKA.sbmlCompartment="compartment";
  PKA.init.sbmlRole="speciesInitialConcentration";
  cAMP.amt.sbmlRole="speciesAmount";
  cAMP.sbmlRole="speciesConcentration";
  cAMP.sbmlCompartment="compartment";
  cAMP.init.sbmlRole="speciesInitialConcentration";
  ERK2.amt.sbmlRole="speciesAmount";
  ERK2.sbmlRole="speciesConcentration";
  ERK2.sbmlCompartment="compartment";
  ERK2.init.sbmlRole="speciesInitialConcentration";
  REGA.amt.sbmlRole="speciesAmount";
  REGA.sbmlRole="speciesConcentration";
  REGA.sbmlCompartment="compartment";
  REGA.init.sbmlRole="speciesInitialConcentration";
  excAMP.amt.sbmlRole="speciesAmount";
  excAMP.sbmlRole="speciesConcentration";
  excAMP.sbmlCompartment="compartment";
  excAMP.init.sbmlRole="speciesInitialConcentration";
  v1.rate.sbmlRole="rate";
  v2.rate.sbmlRole="rate";
  v3.rate.sbmlRole="rate";
  v4.rate.sbmlRole="rate";
  v5.rate.sbmlRole="rate";
  v6.rate.sbmlRole="rate";
  v7.rate.sbmlRole="rate";
  v8.rate.sbmlRole="rate";
  v9.rate.sbmlRole="rate";
  v10.rate.sbmlRole="rate";
  v11.rate.sbmlRole="rate";
  v12.rate.sbmlRole="rate";
  v13.rate.sbmlRole="rate";
  v14.rate.sbmlRole="rate";
}