// 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
// 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 L;
  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;
  real ATP(time) M;
  real trioseP(time) M;
  real ADP(time) M;
  real BPG(time) M;
  real sink = 0 M;
  real vHKPFK(time) katal;
  real vGAPDH(time) katal;
  real vlowpart(time) katal;
  real vATPase(time) katal;
  real vstorage(time) katal;
  real vglycerol(time) katal;

  // equations
  k7t = CGlc*9.65605/6.38837;
  V2t = CGlc*67.7113/6.38837;
  when (time=time.min) ATP = 2.1028274;
  ATP:time = (-2*vHKPFK + 2*vlowpart + -1*vATPase + -2*vstorage)/compartment;
  when (time=time.min) trioseP = 5;
  trioseP:time = (2*vHKPFK + -1*vGAPDH + -1*vglycerol)/compartment;
  when (time=time.min) ADP = 1.49718;
  ADP:time = (2*vHKPFK + -2*vlowpart + vATPase + 2*vstorage)/compartment;
  when (time=time.min) BPG = 2.7550972E-4;
  BPG:time = (vGAPDH + -1*vlowpart)/compartment;
  vHKPFK = V2t*ATP/(1+(ATP/Ki)^q);
  vGAPDH = k3t*trioseP;
  vlowpart = k4*BPG*ADP;
  vATPase = k6*ATP;
  vstorage = k7t*ATP;
  vglycerol = 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.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  trioseP.sbmlRole="species";
  trioseP.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  BPG.sbmlRole="species";
  BPG.sbmlCompartment="compartment";
  sink.sbmlRole="species";
  sink.sbmlCompartment="compartment";
  vHKPFK.sbmlRole="rate";
  vGAPDH.sbmlRole="rate";
  vlowpart.sbmlRole="rate";
  vATPase.sbmlRole="rate";
  vstorage.sbmlRole="rate";
  vglycerol.sbmlRole="rate";
}