// 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
// vHKPFK: 2ATP => 2trioseP 2ADP 
// vGAPDH: trioseP => BPG 
// vlowpart: BPG 2ADP => 2ATP 
// vATPase: ATP => ADP 
// vstorage: 2ATP => 2ADP 
// vglycerol: trioseP => sink 

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

  // variable definitions
  real compartment = 1 volume;
  real k7t(time);
  real V2t(time);
  real q = 2.64;
  real Ki = 1.2446368;
  real k3t = 6.65022;
  real k4 = 1.2651852E5;
  real k6 = 3.2;
  real k8t = .61778;
  real Atot = 3.6;
  real CGlc = 6.38837;
  real CNADH = .33;
  real CNAD = .65;
  private real ATP.amt(time) substance;
  real ATP(time) substance/volume;
  real ATP.init substance/volume;
  private real trioseP.amt(time) substance;
  real trioseP(time) substance/volume;
  real trioseP.init substance/volume;
  private real ADP.amt(time) substance;
  real ADP(time) substance/volume;
  real ADP.init substance/volume;
  private real BPG.amt(time) substance;
  real BPG(time) substance/volume;
  real BPG.init substance/volume;
  private real sink.amt substance;
  real sink substance/volume;
  real sink.init substance/volume;
  real vHKPFK.rate(time) substance/time;
  real vGAPDH.rate(time) substance/time;
  real vlowpart.rate(time) substance/time;
  real vATPase.rate(time) substance/time;
  real vstorage.rate(time) substance/time;
  real vglycerol.rate(time) substance/time;

  // equations
  k7t = CGlc*9.65605/6.38837;
  V2t = CGlc*67.7113/6.38837;
  when (time=time.min) ATP.amt = ATP.init*compartment;
  ATP.amt:time = -2*vHKPFK.rate + 2*vlowpart.rate + -1*vATPase.rate + -2*vstorage.rate;
  ATP = ATP.amt/compartment;
  ATP.init = 2.1028274;
  when (time=time.min) trioseP.amt = trioseP.init*compartment;
  trioseP.amt:time = 2*vHKPFK.rate + -1*vGAPDH.rate + -1*vglycerol.rate;
  trioseP = trioseP.amt/compartment;
  trioseP.init = 5;
  when (time=time.min) ADP.amt = ADP.init*compartment;
  ADP.amt:time = 2*vHKPFK.rate + -2*vlowpart.rate + vATPase.rate + 2*vstorage.rate;
  ADP = ADP.amt/compartment;
  ADP.init = 1.49718;
  when (time=time.min) BPG.amt = BPG.init*compartment;
  BPG.amt:time = vGAPDH.rate + -1*vlowpart.rate;
  BPG = BPG.amt/compartment;
  BPG.init = 2.7550972E-4;
  sink.amt = sink*compartment;
  sink = sink.init;
  sink.init = 0;
  vHKPFK.rate = V2t*ATP/(1+(ATP/Ki)^q);
  vGAPDH.rate = k3t*trioseP;
  vlowpart.rate = k4*BPG*ADP;
  vATPase.rate = k6*ATP;
  vstorage.rate = k7t*ATP;
  vglycerol.rate = k8t*trioseP;

  // variable properties
  compartment.sbmlRole="compartment";
  k7t.sbmlRole="parameter";
  V2t.sbmlRole="parameter";
  q.sbmlRole="parameter";
  Ki.sbmlRole="parameter";
  k3t.sbmlRole="parameter";
  k4.sbmlRole="parameter";
  k6.sbmlRole="parameter";
  k8t.sbmlRole="parameter";
  Atot.sbmlRole="parameter";
  CGlc.sbmlRole="parameter";
  CNADH.sbmlRole="parameter";
  CNAD.sbmlRole="parameter";
  ATP.amt.sbmlRole="speciesAmount";
  ATP.sbmlRole="speciesConcentration";
  ATP.sbmlCompartment="compartment";
  ATP.init.sbmlRole="speciesInitialConcentration";
  trioseP.amt.sbmlRole="speciesAmount";
  trioseP.sbmlRole="speciesConcentration";
  trioseP.sbmlCompartment="compartment";
  trioseP.init.sbmlRole="speciesInitialConcentration";
  ADP.amt.sbmlRole="speciesAmount";
  ADP.sbmlRole="speciesConcentration";
  ADP.sbmlCompartment="compartment";
  ADP.init.sbmlRole="speciesInitialConcentration";
  BPG.amt.sbmlRole="speciesAmount";
  BPG.sbmlRole="speciesConcentration";
  BPG.sbmlCompartment="compartment";
  BPG.init.sbmlRole="speciesInitialConcentration";
  sink.amt.sbmlRole="speciesAmount";
  sink.sbmlRole="speciesConcentration";
  sink.sbmlCompartment="compartment";
  sink.init.sbmlRole="speciesInitialConcentration";
  vHKPFK.rate.sbmlRole="rate";
  vGAPDH.rate.sbmlRole="rate";
  vlowpart.rate.sbmlRole="rate";
  vATPase.rate.sbmlRole="rate";
  vstorage.rate.sbmlRole="rate";
  vglycerol.rate.sbmlRole="rate";
}