// 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: DHAP GAP => MGO 
// v2: MGO GSH => SDLGSH 
// v3: SDLGSH => GSH DLac 
// v4: MGO => Acetol 

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

  // variable definitions
  real compartment = 1 volume;
  real k1 = .01636;
  real k2 = 6.6E-4;
  real V1 = 186.45;
  real km1 = 3.56;
  real km2 = 1.64;
  real V2 = 8.09;
  real km3 = .91;
  real kout = 1;
  real V3 = 17.85;
  real NADPH = 1.7;
  real km4 = .65;
  real km5 = .075;
  private real DHAP.amt substance;
  real DHAP substance/volume;
  real DHAP.init substance/volume;
  private real GAP.amt substance;
  real GAP substance/volume;
  real GAP.init substance/volume;
  private real MGO.amt(time) substance;
  real MGO(time) substance/volume;
  real MGO.init substance/volume;
  private real GSH.amt(time) substance;
  real GSH(time) substance/volume;
  real GSH.init substance/volume;
  private real SDLGSH.amt(time) substance;
  real SDLGSH(time) substance/volume;
  real SDLGSH.init substance/volume;
  private real DLac.amt substance;
  real DLac substance/volume;
  real DLac.init substance/volume;
  private real Acetol.amt substance;
  real Acetol substance/volume;
  real Acetol.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;

  // equations
  DHAP.amt = DHAP*compartment;
  DHAP = DHAP.init;
  DHAP.init = 2.5;
  GAP.amt = GAP*compartment;
  GAP = GAP.init;
  GAP.init = .12;
  when (time=time.min) MGO.amt = MGO.init*compartment;
  MGO.amt:time = v1.rate + -1*v2.rate + -1*v4.rate;
  MGO = MGO.amt/compartment;
  MGO.init = 0;
  when (time=time.min) GSH.amt = GSH.init*compartment;
  GSH.amt:time = -1*v2.rate + v3.rate;
  GSH = GSH.amt/compartment;
  GSH.init = 4;
  when (time=time.min) SDLGSH.amt = SDLGSH.init*compartment;
  SDLGSH.amt:time = v2.rate + -1*v3.rate;
  SDLGSH = SDLGSH.amt/compartment;
  SDLGSH.init = 0;
  DLac.amt = DLac*compartment;
  DLac = DLac.init;
  DLac.init = 0;
  Acetol.amt = Acetol*compartment;
  Acetol = Acetol.init;
  Acetol.init = 0;
  v1.rate = k1*GAP+k2*DHAP;
  v2.rate = V1*MGO*GSH/((km1+MGO)*(km2+GSH));
  v3.rate = V2*SDLGSH/(km3+SDLGSH);
  v4.rate = V3*NADPH*MGO/((km4+MGO)*(km5+NADPH));

  // variable properties
  compartment.sbmlRole="compartment";
  k1.sbmlRole="parameter";
  k2.sbmlRole="parameter";
  V1.sbmlRole="parameter";
  km1.sbmlRole="parameter";
  km2.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  km3.sbmlRole="parameter";
  kout.sbmlRole="parameter";
  V3.sbmlRole="parameter";
  NADPH.sbmlRole="parameter";
  km4.sbmlRole="parameter";
  km5.sbmlRole="parameter";
  DHAP.amt.sbmlRole="speciesAmount";
  DHAP.sbmlRole="speciesConcentration";
  DHAP.sbmlCompartment="compartment";
  DHAP.init.sbmlRole="speciesInitialConcentration";
  GAP.amt.sbmlRole="speciesAmount";
  GAP.sbmlRole="speciesConcentration";
  GAP.sbmlCompartment="compartment";
  GAP.init.sbmlRole="speciesInitialConcentration";
  MGO.amt.sbmlRole="speciesAmount";
  MGO.sbmlRole="speciesConcentration";
  MGO.sbmlCompartment="compartment";
  MGO.init.sbmlRole="speciesInitialConcentration";
  GSH.amt.sbmlRole="speciesAmount";
  GSH.sbmlRole="speciesConcentration";
  GSH.sbmlCompartment="compartment";
  GSH.init.sbmlRole="speciesInitialConcentration";
  SDLGSH.amt.sbmlRole="speciesAmount";
  SDLGSH.sbmlRole="speciesConcentration";
  SDLGSH.sbmlCompartment="compartment";
  SDLGSH.init.sbmlRole="speciesInitialConcentration";
  DLac.amt.sbmlRole="speciesAmount";
  DLac.sbmlRole="speciesConcentration";
  DLac.sbmlCompartment="compartment";
  DLac.init.sbmlRole="speciesInitialConcentration";
  Acetol.amt.sbmlRole="speciesAmount";
  Acetol.sbmlRole="speciesConcentration";
  Acetol.sbmlCompartment="compartment";
  Acetol.init.sbmlRole="speciesInitialConcentration";
  v1.rate.sbmlRole="rate";
  v2.rate.sbmlRole="rate";
  v3.rate.sbmlRole="rate";
  v4.rate.sbmlRole="rate";
}