// This model generated automatically from SBML
// WARNING: Renamed variable source to source_ due to JSim namespace restrictions

// unit definitions
import nsrunit;
unit conversion off;

// SBML property definitions
property sbmlRole=string;
property sbmlName=string;
property sbmlCompartment=string;

// SBML reactions
// vinGlc: source 
// vGlcTrans: GlcX => 59Glc 
// vHKPFK: Glc 2ATP => 2trioseP 2ADP 
// vGAPDH: trioseP NAD => BPG NADH 
// vlowpart: BPG 2ADP => ACA 2ATP 
// vADH: ACA NADH => NAD 
// vATPase: ATP => ADP 
// vstorage: Glc 2ATP => 2ADP 
// vglycerol: trioseP NADH => NAD 
// vdifACA: 59ACA => ACAX 
// voutACA: ACAX => sink 

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

  // variable definitions
  real compartment = 1 L;
  real k0 = 2.82529;
  real V1 = 49.1132;
  real Ktrans = .04;
  real V2 = 10.5992;
  real q = 2.64;
  real Ki = 1.24464;
  real k3 = 10.2311;
  real k4 = 126519;
  real k5 = 17.2154;
  real k6 = 3.2;
  real k7 = 1.5115;
  real k8 = 1.87206;
  real k9 = .690539;
  real k10 = .869329;
  real Atot = 3.6;
  real Ntot = .98;
  real yvol = 59;
  real GlcX0 = 24;
  real source_ = 0 M;
  real GlcX(time) M;
  real Glc(time) M;
  real ATP(time) M;
  real trioseP(time) M;
  real ADP(time) M;
  real NAD(time) M;
  real BPG(time) M;
  real NADH(time) M;
  real ACA(time) M;
  real ACAX(time) M;
  real sink = 0 M;
  real vinGlc(time) katal;
  real vGlcTrans(time) katal;
  real vHKPFK(time) katal;
  real vGAPDH(time) katal;
  real vlowpart(time) katal;
  real vADH(time) katal;
  real vATPase(time) katal;
  real vstorage(time) katal;
  real vglycerol(time) katal;
  real vdifACA(time) katal;
  real voutACA(time) katal;

  // equations
  when (time=time.min) GlcX = 6.7194493;
  GlcX:time = (-1*vGlcTrans)/compartment;
  when (time=time.min) Glc = 5;
  Glc:time = (59*vGlcTrans + -1*vHKPFK + -1*vstorage)/compartment;
  when (time=time.min) ATP = 2.1028392;
  ATP:time = (-2*vHKPFK + 2*vlowpart + -1*vATPase + -2*vstorage)/compartment;
  when (time=time.min) trioseP = 7.8474335;
  trioseP:time = (2*vHKPFK + -1*vGAPDH + -1*vglycerol)/compartment;
  when (time=time.min) ADP = 1.49716;
  ADP:time = (2*vHKPFK + -2*vlowpart + vATPase + 2*vstorage)/compartment;
  when (time=time.min) NAD = .65;
  NAD:time = (-1*vGAPDH + vADH + vglycerol)/compartment;
  when (time=time.min) BPG = 2.7551108E-4;
  BPG:time = (vGAPDH + -1*vlowpart)/compartment;
  when (time=time.min) NADH = .33000016;
  NADH:time = (vGAPDH + -1*vADH + -1*vglycerol)/compartment;
  when (time=time.min) ACA = 8.3327608;
  ACA:time = (vlowpart + -1*vADH + -59*vdifACA)/compartment;
  when (time=time.min) ACAX = 1.3121754;
  ACAX:time = (vdifACA + -1*voutACA)/compartment;
  vinGlc = k0/yvol*(GlcX0-GlcX);
  vGlcTrans = V1/yvol*GlcX/(Ktrans+GlcX);
  vHKPFK = V2*Glc*ATP/(1+(ATP/Ki)^q);
  vGAPDH = k3*trioseP*NAD;
  vlowpart = k4*BPG*ADP;
  vADH = k5*ACA*NADH;
  vATPase = k6*ATP;
  vstorage = k7*Glc*ATP;
  vglycerol = k8*trioseP*NADH;
  vdifACA = k9/yvol*(ACA-ACAX);
  voutACA = (k0+k10)/yvol*ACAX;

  // variable properties
  compartment.sbmlRole="compartment";
  k0.sbmlRole="parameter";
  V1.sbmlRole="parameter";
  Ktrans.sbmlRole="parameter";
  V2.sbmlRole="parameter";
  q.sbmlRole="parameter";
  Ki.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";
  Atot.sbmlRole="parameter";
  Ntot.sbmlRole="parameter";
  yvol.sbmlRole="parameter";
  GlcX0.sbmlRole="parameter";
  source_.sbmlRole="species";
  source_.sbmlCompartment="compartment";
  GlcX.sbmlRole="species";
  GlcX.sbmlCompartment="compartment";
  Glc.sbmlRole="species";
  Glc.sbmlCompartment="compartment";
  ATP.sbmlRole="species";
  ATP.sbmlCompartment="compartment";
  trioseP.sbmlRole="species";
  trioseP.sbmlCompartment="compartment";
  ADP.sbmlRole="species";
  ADP.sbmlCompartment="compartment";
  NAD.sbmlRole="species";
  NAD.sbmlCompartment="compartment";
  BPG.sbmlRole="species";
  BPG.sbmlCompartment="compartment";
  NADH.sbmlRole="species";
  NADH.sbmlCompartment="compartment";
  ACA.sbmlRole="species";
  ACA.sbmlCompartment="compartment";
  ACAX.sbmlRole="species";
  ACAX.sbmlCompartment="compartment";
  sink.sbmlRole="species";
  sink.sbmlCompartment="compartment";
  vinGlc.sbmlRole="rate";
  vGlcTrans.sbmlRole="rate";
  vHKPFK.sbmlRole="rate";
  vGAPDH.sbmlRole="rate";
  vlowpart.sbmlRole="rate";
  vADH.sbmlRole="rate";
  vATPase.sbmlRole="rate";
  vstorage.sbmlRole="rate";
  vglycerol.sbmlRole="rate";
  vdifACA.sbmlRole="rate";
  voutACA.sbmlRole="rate";
}